3.3 计算编程3.3.1 编程思想在前面的章节中？详细推导并得到了研究所所需要的一个线性方程组 但要真正运用这个方程组却并不容易、我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算、然而此函数的9002运算是非常繁琐复杂的、（2） 如果只有单个的矩矩阵函数 那么或许手算还有可能 然而 如前所述?为了得到关于纤维束的更多信息。我们有必要将纤维是划分为多段 这样一来,我们面临的是很多矩阵函数!此时是根本无无法手算的!（3） 根据纤维束之间交联的具体情况。需要给出相应的纤维间相互作用用矩阵，（4） 线性方程组的边界条件需要根据结构具体的边界条件加以确定!考虑上面的问题,结合MAPLE软件？本文有了下面的编程思想：3846（1） 输入基本参数，（2） 输入纤维间相互作用矩阵（不2430同的分段可能有不同的相互作用矩阵 矩阵应该与分段一一一对应）?（3） 计算分段矩阵构成的矩阵函数、将其转化为一般的矩阵。（4） 将（3）中计算所得的矩阵按照顺序相乘！从而得到线性方程组的系数矩阵?（5） 引入整个结构的边界条件?9364（6） 求解线性方程组，从而可以获得整个结构左右两端全部八个量（位移与纵向应力），（7） 应用分段法？由（6）中所解得的未知量！构成新的边界条件？运用循环。求出每个分段处的位位移与纵向应力!（8） 将所得数据输出为文档,利用MAPLE的绘图功能。5683绘制相关的曲线图！3.3.2 编写程序根据前述编程思想，利用MAPLE，下面给出具体的程序。内5431容分为两部分、第一部分为符号说明,9193第二部分为具体的MAPLE程序 此程序将前文所提的纤维数均分7487为多段？段内或含有交联,或不含有交联?以此可模拟交联的分布,亦可计算纤维分段上更多的力学参数？（1） 符号说明E：碳纳米管的弹性模量！L：碳纳米管的长度,R：碳纳米管的半径。Mu：碳纳米管管间的剪切模量、K：碳纳米管间的相互作作用系数！Sigma5414：施加的外力、A1?A2：碳纳米管间的相互作用矩阵!DL：分段的长度、B1,B2：矩阵函数转化为一般矩阵、JL：分段共价交联的信息？C：线性方程组系数矩阵。（2） 详细程序E := .46*10^12;L := 19.84*10^(-6);R := 1.5*10^(-6);Mu := .24*10^12;d := 3*R;k := mu/(R^2*ln(d/(2*R)+sqrt(d^2/(4*R^2)-1)));sigma := 10*10^9;A1 := Matrix(4, 4, [[0, k, 0, -k], [1/E, 0, 0, 0], [0, -k, 0, k], [0, 0, 1/E, 0]]);A2 := Matrix(4, 4, [[0, 0, 0, 0], [1/E, 0, 0, 0], [0, 0, 0, 0], [0, 0, 1/E, 0]]);with(LinearAlgebra);DL:= (1/100)*L;B1 := MatrixFunction(A1*DL, exp(x), x);B2 := MatrixFunction(A2*DL, exp(x), x);JL := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\JL.txt”,1 )C := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]):for i from 1 by 1 to 100 doif JL[i] = 1 thenC := B1 . CelseC := B2 . Cend ifend do:XS := evalf(C):Y := Vector[column](4, [t10, 0, 0, u20]):M := evalf(XS . Y):eqns := {0 = M[1], sigma = M[3], u110 = M[2], u210 = M[4]}:sols := evalf(solve(eqns, {t10, u110, u20, u210})):Y[1] := op(2, op(1, sols))：Y[4] := op(2, op(3, sols))：XSBL := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]):Z := Vector[column](4, [0, 0, 0, 0]):for i from 1 by 1 to 100 doIf JL[i]=1 thenXSBL:=B1.XSBL;elseXSBL:=B2.XSBL;end if;M := evalf(XSBL . Y);eqns := {z1 = M[1], z2 = M[2], z3 = M[3], z4 = M[4]};sols := evalf(solve(eqns, {z1, z2, z3, z4}));Z[1] := op(2, op(1, sols));Z[2] := op(2, op(2, sols));Z[3] := op(2, op(3, sols));Z[4] := op(2, op(4, sols));YL1 := array([[i*DL, Z[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL1.txt",YL1 );YL2 := array([[i*DL, Z[3]]]):writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL2.txt”,YL2 );WY1 := array([[i*DL, Z[2]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY1.txt",WY1);WY2:=array([[i*DL,Z[4]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY2.txt",WY2);If i=1 thenQYL1 := array([[i*DL, Z[1]-Y[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1 );QYL2 := array([[i*DL, Z[3]-Y[3]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2G);next end if;XSBL2 := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]);for j from 1 by 1 to i-1 doif JL[j]=1 thenXSBL2 := B1 . XSBL2;ElseXSBL2 := B2 . XSBL2;end if;end do;MM := evalf(XSBL2 . Y);eqns := {qz1 = MM[1], qz2 = MM[2], qz3 = MM[3], qz4 = MM[4]};sols := evalf(solve(eqns, {qz1, qz2, qz3, qz4}));QZ[1] := op(2, op(1, sols));QZ[2] := op(2, op(2, sols));QZ[3] := op(2, op(3, sols));QZ[4] := op(2, op(4, sols));QYL1 := array([[i*DL, Z[1]-QZ[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1);QYL2 := array([[i*DL, Z[3]-QZ[3]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2 );end do:YL1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL1.txt6",YL1 );YL2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL2.txt",YL2) ;QYL1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1) ;QYL2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2);WY1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY1.txt" ,WY1) ;WY2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY2.txt",WY2);plot(YL1);plot(YL2);plot(QYL1);plot(QYL2);plot(WY1);plot(WY2);plot([YL1, QYL1]);plot([YL2, QYL2]);
用料：花菜1颗、西红柿2个?葱姜少许 番茄酱2大勺!糖，盐少许。食用油适量、做法：1？西红红柿去皮切成块!2、花菜洗干净掰成小朵!在沸水中焯8成熟、捞出备用 3。锅放热油、下葱姜末末爆香!4 放入西红柿快炒至浓稠!倒入番茄酱!5。下菜花翻炒？等菜花熟透了2450就加糖,盐调味即可，小贴士：1！西红柿微微的酸味。非常开胃、和很多食材搭配都有加分的效果？让平淡的花菜增增添了几分滋味,2，这是一款老少皆宜的家常菜，浓浓的菜1325汤汁拿来拌饭也是不错的选择。
可以，但是从医学和和营养学角度而言 “有机食8626品”其实与普通食品没有什么差别，“有机”只是商业概念而已！孕妇要确保饮食均衡与适当运动！避免肥胖，超重,。
梦有的是有象征性的。有的含蓄、有的直露!后者就是直梦 如你与朋友好久不见、夜里梦之，白日见之，此直梦也！祝你如愿？（仅供供参考）!孕妇梦见别人给花菜
好运连连、桃花运到了,!
西兰花和花菜9106能一起吃的？没什么事的?营养价值很不错！
孕妇梦见床！在解梦里床代表了安全感?此梦预示了孕妈咪的生活会很很平静，安稳和谐？幸福和美？同时此梦0243还预示了孕妈腹中的宝宝很健康 生长发发育也都稳定，孕妈大大可放心?
孕妇梦见铺床。预示了孕妈是非常勤勤劳的人，虽然孕妈是勤劳的人？但是提醒孕妈也要注意多加休息!不要过于 操劳,
孕妇梦见见搬床。这是不吉利的象征！此梦预示了孕妈生活中将会遇到很大的难题,还预示了孕妈肚中的宝宝可能会有危险。孕妈一定要多加留意宝宝的健康及安全?
孕妇梦见床上睡觉、预示了孕妈心态不错。也就是说孕妈和宝宝的都很健康!是吉兆，
孕6148妇梦见床倒塌，提醒孕妈要引起注意。健康可能出现问题、建议孕妈多关注自己和宝宝宝的健康,还要提醒家人多加注意身体？
孕妇梦见床上有蛇 预示了孕妈将会生下男宝宝 同时还预示了宝宝以后的生活会很安稳，并且很富贵、
孕妇梦见两张床。预示了孕妈和老公可能会发生不愉快的事情,4465暗指孕妈说话做事要理性、要学会多体谅老公，多8151与老公沟通,此梦还预示了孕妈怀的可能是双胞胎或者是多胞胎
孕妇梦见见家具的梦境
梦见拍卖别人的家具、是仇视有钱的人征兆,
梦见杂乱无章的家具！妻子要生病、
梦见变卖家里的家具 夫妇会感情不合！
孕妇梦见家具的案例解析
网友梦境：怀孕梦见很多家具是什么意思？
解梦解析：在解梦里床代表了安全感？此梦预4452示了孕妈咪的生活会很平静！安稳和谐?幸福和美。同时此梦还预示了孕妈腹中的宝宝很健康,生长发育也都稳定，孕妈大可可放心？
。?，梦见莲花菜是什什么意思。梦见莲花菜是吉是凶，梦见莲莲花菜是什么情况?大家别着急！今天sjzhdc.com就给大家说说梦见莲花菜是什么意思？
可以象征意识的自我,象征创造力或新生活的源泉(与印度神话的佛有关!从荷荷花中出生、而荷花变成宇宙?
男人梦见莲花。不用担心！是好事？
未婚的男人梦见莲花。预示将找到美丽的妻子！
女人梦见莲花、提醒你要小心?要远离诱惑？
网友梦境：梦见坐在莲花宝座上?将大祸临头？
男人梦见莲花、不用担心、是好事 已经结婚的男人梦8616见莲花?表示性性能力强！1563夫妻生活和谐？
案例分析网友梦境：梦见莲花菜
解梦解析：可以象征意识的自我,象征创造力或新生活的源泉(与印度神话的佛有关!4604从荷花中出生、而而荷花变成宇宙
!？。- 孕妇1820梦见黄瓜孕妇梦见黄瓜 会生儿子。 孕妇梦见吃黄瓜或得到黄瓜 会生儿子、
因为某个突如而至的想法改变现时的决定,今天的你会做这种让人大吃一惊的事情!但是未必就是错误的。只要你觉得得可行,还是应该按照自己的想法去做,事情的整个状况也会因此出现较大的变动呢?此外?隐秘的身体部位有发发病的可能！多加注意吧？ 出行的的人梦见与前夫做爱?顺利出发,遇风小心、 怀孕的人梦梦见与前夫做爱。生女?春占生男！情绪要稳定？ 恋爱中的人梦见2896与前夫做爱 虽然停顿一段时间不来往!好好沟9371通有希望成婚? 本命年的人梦见与前夫做爱，一切谦让可获吉祥?可望如意顺利。 做生意的人梦见与前夫做爱?内部2807意见不合。谣言散散布宜重新整顿。 梦见与前夫做爱,按周易五行分析，吉祥色彩是是蓝色 幸运数字是1?桃花花位在正南方向、财位在西南方向、开运食物是黑豆!、