3.3 计算编程3.3.1 编程思想在前面的章节中。详细推导并得到了了研究所需要的一个线性方程组!但要真正运用这个方程组却并不容易。我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算、然而此函数的运算是非常繁琐复杂的！（2） 如果只有单个的矩阵函数！那么或许手算还有可能、然而，如前所述，为了得到关于纤维束的更多信息。我们有必要将纤维是划分为多段、这样一来、我们们面临的是很多矩阵函数、此时是根本无法手算的，（3） 根据纤维束之间交联的具体体情况？需要给出相应的纤维间相互作用矩阵、（4） 线性方程组的边界条件需要根据结构具体的边9415界条件加以确定、考虑上面的问题、结合MAPLE软件。本文有了下面的编程思想：（1） 输入基本参数、（2） 输入纤维间相互作用矩阵（不同的分段可能有不同的相互作用矩阵！矩阵应该与分段一一对应），（3） 计算分段矩阵构成的矩矩阵函数、将其转化为一般的矩阵。（4） 将（3）中计算所得的矩阵按照顺序相乘。从而得到线性方程组的系数矩阵、（5） 引入整个结构的边界条件。（6） 求解2508线性方程组？从而可以获得整个结构左右两端全部八个量（位移与纵向应力）、（7） 应用^分分段法。由（6）中所解得的未知量？构成新的边边界条件，运用循环，求出每个分段处的位移与纵向应力，（8） 将所得数据输出为文档、利利用MAPLE的绘图功能？绘7489制相关的曲线图，3.3.2 编写程序根据前述编程思想。利用MAPLE、下面给出具体的程序，内容分为两部分。第9107一部分为符号说明,第二部分为具体体的MAPLE程序。此程序将前文所提的纤维数均分为多段，段内或含^有有交联,或不含有交联、以此可模模拟交联的分布!亦可计算纤维分段上更多的力学参数。（1） 符号说明E：碳纳米管的弹性模量、L：碳纳米管的长度、R：碳纳米管的半径。Mu：碳纳米管间的剪切模量、K：碳纳米管间的的相互作用系数？Sigma：施加的外力。A1、A2：：碳纳米管间的相互作用矩阵、DL：分段的长度？B1！B2：矩阵函数转化为一4641般矩阵!JL：分段共价交联的信息，C：线性方程7271组系数矩阵,（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]);！

　　番茄，要求茎秆粗壮、节间短！幼苗高15厘米左右!3-5叶一心！根系发达！健康、叶片肥厚！叶色浓绿，生长势强、无病虫害！无损伤！大小均匀一致。苗龄约为25天左右，！风水罗盘种类

　　·梦见浇苗！丈夫会生病、·男人梦梦见女人浇苗、是发财的预兆！·农民梦见见浇苗？庄稼会丰收。实际上日有所思夜有所梦！信则有、不信则无、不必太在意的！属牛双鱼座今日运势

　　◇ 西红柿是好运、 　　◇ 西红柿象征着红运当头, 　　◇ 男人梦见吃西红柿，意味着交好运！ 　　◇ 已婚女人梦见吃西红柿！意味着身体健康， 　　◇ 未婚女子梦见吃西红柿、意味着嫁给身强体壮的的丈夫， 　　◇ 未婚男子梦见吃西红柿、^会会娶一个温柔体贴的妻子、，几两几钱算命法

　　孕妇可以吃西红柿、吃西红柿给孕^妇妇带来很多好处。西红柿性凉！味甘酸、富含胡萝卜素。西红柿红素！维生素C，蛋白质等营养物质，这些营养物质对人体体大有好处、^对对孕妇而言也是多多益善。更有以下功效：1、美容养颜：西红柿含有的胡萝卜素和西红柿红素！有助于展平皱纹！使皮肤细嫩光滑。能有效预防孕妇担忧⁸⁷⁷⁷的妊娠斑和妊娠纹等肌肤问题,2，防癌：西红柿里的西红柿红素能通过有效清除体内的自由基，预防和修复细胞损伤，抑制制DNA的氧化！从而降低癌症的发生率。3！抗衰老：西红柿中的维C、有生津止渴、健胃消食等功效，最重要的是它有很强的抗氧化剂，给人体补充西红柿红素的同时、可以帮助身体抵抗各种因自由基引起的退化老化性疾病，4！降脂降压：西红柿里含的维生素C。芦丁！西红柿红红素及果酸,可降低2215血胆固醇。预防动脉粥样硬化及冠心病，另其含有大量的钾及碱性矿物质！能促进血血中钠盐的排出！有降压、利尿、消肿作用，因此、孕妇吃西红柿能有效预防妊娠高血压、还能缓解妊娠水肿等症状、。惹事生非的生肖

　　梦见吐西红西红柿籽 　　77 　　梦见吐西红西红柿籽意味着： 　　不要给自己得过且过的机会，这两天蛮有点“老油条”的样子哦。对什么都不太打得起精神，特别是熟悉的事情！以为自己蛮有把握！轻松对付就可以了、提醒一下、平时需要用到的技能虽然有他人的帮助可以完成！但自己能掌握的话就更强大了！ 　　梦见吐西红西红柿籽的吉凶^：： 　　得部下拥载？及长辈引进，而得成功发展（五行顺相生。从上生下、配置良好），易得财利！名誉！事业隆昌。长寿少病之兆！【大吉昌】，白蜡金命的人命运如何

　　您好、皮皮虾不可以和西红柿一起吃！两样食物相克！产产生砒霜会中毒,结婚嫁娶择日

　　事务繁忙就令人做梦！多做梦的其中多有虚幻、不必放在心上、要要知道人的大脑和小脑在正常情况下是不休息的。所以大多数 人都会做梦。勿忘点采纳，鸡年迎来高峰的生肖

　　生女儿！鱼会游走的！！朱砂本命佛佩戴禁忌孕妇梦到买西红柿苗

　　暗示你你孩子的名字、风水门前变压器化解