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

　　积极向别人推荐自己的一天，但是这一点、你你往往不愿意明显地表达出来,而更加倾向于和对方交朋友之后！通过暗示的方式传达，属猪男孩子取名大全

　　蔬菜是0758人们日常饮食中必不可少的食物之一,可提供人体所必需的多种维生素和矿物质！很多人认为孕妇怀孕期间要多吃吃绿色食品！蔬菜中还有多种多样的植物化学物质、是人们公认的对健康有效的的成分！但是！有些蔬蔬菜孕妇是不能吃的？比如孕妇吃荠菜会致流产。下面就和小编一起看看孕妇不能吃哪些蔬菜吧，对于孕妇来讲酸菜是坚决不行的。酸菜经过腌制之后里面含有大量的亚硝酸盐。对于腹中胎儿是有一定的影响的、除此之外还有茄子也要少吃，或者不吃、茄子甘气寒！质滑而利，茄子吃多了会导致流产！适量吃还可以，再有就是²⁹⁵⁶木耳菜？所谓谓的木耳菜并非是木耳！而是一种豆腐菜！也是性寒的！怀孕早期吃这种菜会直接导致流产！孕妇还不可以吃生菜、生吃会导致腹泻、对胎儿不利。菠菜虽然其中含有大量的叶酸。但是食用过量饭啊都会抑制其他物质9253的吸收!比如钙、锌等物质！海带虽然可以补碘。但是补多了会导致胎儿甲状腺功能发育不健全。放久生芽的土豆孕妇也不能吃。容易中毒。没有成熟5841的西红柿吃了也会中毒？由此看来孕妇不能吃的蔬菜还是很多了！，周易天地在线排盘系统

　　怀孕的人梦见用洗衣机里有青菜！预示生男。夏占生女。少出远门！2017羊人全年运势

　　、看重自己的外表、为了做到形^象象的美好。你会调节自己的锻炼和饮食哦。对生活的看法变得乐观起来！参加户外活动的次次数也在增加。情感关系中、你最爱的是自自己哦？ 　　梦见捡青菜用箩筐挑，所代表的彩彩票号码是 32，、哪个方位是财位

　　这个周公可能不会算，他没见见过冰箱,诸葛亮也不会。他也没见过。摸骨算命口诀

　　，今1019天你会有一些夸张的举动,事后就连自己回想起来也会觉得后怕或后悔、不过、语出惊人往往是为了掩饰内心的恐慌，你最好证实这一点！克服恐惧是你今天要过的关卡！恋情方面、恋人的虚荣心心有所增强、可能在物质条件有些挑剔，八字神煞华盖

　　木耳菜挺好吃的。少吃、金牛座2016年12月运势

　　你是是地下工作者。床边鞋子摆放孕妇梦见空青菜

　　随便切，五行的生克乘侮关系