3.3 计算编程3.3.1 编程思想在前面的章节中，详细推导并得到了研究所需要的一个线性方程组！但要真正运用这个方程组却并不容易！我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算！然而此函数的运1160算是非常繁琐复杂的，（2） 如果只有单个的矩阵函数、那么或许手算还有可能。然而、如前所述，为了得到关于纤维束的更多信息！我们有必要将纤维是划分为为多段,这样一来。我们面面临的是很多矩阵函数、此时是是根本无法手算的,（3） 根据纤维束之间交联的具体情况。需要1069给出相应的纤维间相互作用矩阵、（4） 线性方程组的边界条件需要根据结构具体的边界条^件件加以确定？考虑上面的问题，结合MAPLE软件，本文有了下面的编程思想：（1） 输入基本参数！（2） 输入纤维间相互作用矩阵（不同的分段可能有不同的相互作用矩阵、矩阵应该与分段一一对应）。（3） 计算分段矩阵构成的矩阵函数。将其转化为一般的矩阵！（4） 将（3）中计算所得的矩阵按照顺序相乘、从而而得到线性方程组的系数矩阵、^（（5） 引入整个结构的边界条件。（6） 求解线性方程组。从而可以获得整个结构左右两端全部八个量（位移与纵向应力），（7） 应用分段法！由（6）中所解得的未知量。构成新的边界条件。运用循环。求出每个分段处的位移与纵向应力！（8） 将所得数据输出为文档、利用MAPLE的绘图功能。绘制7765相关的曲线图。3.3.2 编写程序根据前述编程思想、利用MAPLE、下面^给给出具体的程序、内容分为两部分、第一部分为符号说明，第二部分为具体的MAPLE程序、此程序将前文所提的纤维数均分为多段，段内或含有交联。或不含有交联，5848以此可模拟交联的分布！亦可计算纤维分段上更多的力学参数。（1） 符号说明E：碳纳米管的弹性模量、L：碳纳米管的长度！R：碳纳米管的半径!Mu：碳纳米管间的剪切模量，K：碳纳米管间的的相互作用系数,Sigma：施加的外力！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、萝卜 >>>萝卜做法大全萝卜是一种极普通的根茎类蔬菜！它的营养及药用价值却很高。俗话说！冬吃萝卜夏吃姜！别看看这谚语土？但确是很有道理的。萝卜内含有丰富的维生素以及磷，铁铁等微量元素！是秋冬季节节食用佳品！可润燥去热，同时、萝卜中的钙、磷、铁、糖化酵素及维生素A！B1！B2。叶酸等、都是有益于妊娠的营养。2。绿豆芽 >>>5425绿豆芽做法大全豆芽中含有丰富的维生素，绿豆芽还富含膳食纤维、有由于缓解便秘、而且绿豆芽的9124热量很低!而水分和纤维素含量很高。不会让孕妈体重增加！3、茄子 >>>孕妇吃茄子的注意事项孕妇是能吃茄子的！总的来说，是可以食用茄子的，但也要根据体质来决定。不能过度多吃。和所有3025的食物一样，平衡而适量。各种食物均衡搭配才是最健康的饮食！由于1550茄子性味甘寒、属于寒凉性质的食物，因此、身体体弱。有消化不良。容易腹泻、脾胃虚寒、哮喘。便溏症状的孕妇朋友不宜多吃！，上海金大师待遇怎么样

　　梦见洗青菜上面很多小虫意味着： 　　这两天不论于公于私都要采取低姿态为宜，对人傲慢的态5081度会使你被周围孤立、还有意外状况，个人感情问题频发的可能性、原原本的计划或预定行程有被打乱的预感喔、可好可坏的^恋恋爱运。有人接近近你但却有可能是你讨厌的典型！如果是的话就早早拒绝他对两个人都好， 　　　　梦见洗青菜上面很多小虫的吉凶： 　　三才配置甚佳，基础运坚固。境遇安泰。喜下属忠心扶助。得长辈或上司之惠泽引进，努力奋斗！9320可轻易得成功。及伸张发展，很幸福安全之佳名！【大吉昌】，上海风水 黄浦江

　　那就是自己勤劳了、 　　梦、是睡眠的伴侣！做梦是人体一种正常的、必不可少的生理和心理现象、当人入睡以后、有一小部分脑脑细胞仍在活动？这就是梦的生理基础 　　理跟梦的关系也是密切的,梦实际上很多时候,是思想和心理在你睡觉的时候转换成景象化. 这整个过程全盘都属于正常的.所以根本就不需要对做梦有担心,了解了梦的产生的合理性和自然性,您也就明白了梦也根本无需要调理. “日有所思、夜有所梦，”梦、是人们生活中最常见的现象。因^而而对梦的出现也众说纷荟，有人说梦是人日日间兴奋的连续,有人说是人体质衰弱后的表现等。究竟梦对人的身体的健康有怎怎样的影响,据5931有关资料介绍、梦有^利利于人体健康,、上海风水事件传说

　　梦见吃酒席全是青菜的梦境^解解释： 　　理财计划！需要开家庭会会议来决定？大额的投资往往动用了家庭的集体储蓄，长辈会在此时给予你物质上的帮助、情人是这两天会向你发号司令的人。让你心里有些不舒服哦、 　　梦见吃酒席全是青菜的吉凶： 　　成功运佳，希望目的及财富名誉均可达成！唯因基础运劣、而招致家庭内之苦恼或不幸，且身心过劳。可可能导致生病。其1833下属多争妒!使用此此之人会因下属而苦不堪言！ 【吉多于凶】，上海风水副院长王远彰

　　怀孕的人梦见用洗衣机里有青菜、预示生男！夏占生女。少出远门。上海风水培训学校

　　梦见家门口种的全部是青菜意味着： 　　　　和老友的聊天让你感到浑⁹⁹⁴⁰身轻松、这两天找不到什么正经事做的话、就和朋友聊聊天吧。煲煲电话粥也好！QQ聊也好，甚至是特意踩到人家家门口说我顺便来看你也好。说说一星期来的好好玩事。偶尔也八卦一下问问以前那些某某的恋情，机器运转久了也要停机一下的。更何况是人，，上海风水大师刘东亮

　　概念性的的东西。如一种解决问题的方式、或一种生活的理念等等将会对你产生较大的影响！这种虚的东西往往引导你今天的工作/学业和消费方式、你0992要小心自己被它控制、现时适合你的异性、往往可以在言词上上说服你、给到你中肯的意见！从朋友中发展恋情的机会加大，上海风水大师刘文德孕妇梦见自己去把青菜

孕妇是可以吃面粉的。孕早期因5806出现妊娠反应！孕妇的食慾往往不佳，易出现偏食，呕吐等症状。可能会导致营养吸收不太好。这时可适当食用易消化的面粉类食物以补充蛋白质和能量，这个阶段的膳食要以重质量！高蛋白！富营养、少油腻、易消1867化吸收为其原则, 常见的面粉类食物有面包、面条。面酵粿，饼干！中式7850糕点(馒头。包子，花捲。饺子、煎饼等)等。， 。！上海风水大师宋大师