3.3 计算编程3.3.1 编程思想在前面的章节中、详细推导并得到了研究所需要的一个线性方程组！但要真正运用这这个方程组却并不容易？我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算，然而此函数的运算是非常繁琐琐复杂的！（2） 如果2473只有单个的矩阵函数。那么或1640许手算还有可能，然而！如前所述，为了得到关于纤维束的更多信息！我们有必要将纤维是划分为多段。这样一来?我们面临的是很多多矩阵函数。此时是根本无法手算的。（3） 根据纤维束之间交联的具体情况 需要给出相应的纤维间相相互作用矩阵？（4） 线性方程组的边界条件需要根据结构具具体的边界条件加以确定！考虑上面面的问题！结合MAPLE软件？本文有了下面的编程思想：（1） 输入基本参数!（2） 输入纤维间相互作用矩阵（不同的分段可能有不同的相互作用矩阵,矩阵应该与分段一一对应）、（3） 计算分段矩阵构成的矩阵函数。将其转转化为一般的矩阵？（4） 将（3）中计算所得的矩阵按照顺序相乘、从而得到线性方程组的系数矩阵!（5） 引引入整个结构的边界条件!（6） 求解线性方程组?从而可以获得整个结构左右两端全部八个量（位移与纵向应力）?（7） 应应用分段法 由（6）中所解得的未知量！构成新的边界条件，运用循环？求出每个分段处6241的位移与纵向应力、（8） 将所得数据输出为文档!利用MAPLE的绘图功能！绘制相相关的曲线图!3.3.2 编写程序根据前述编程思想，利用MAPLE,下面给出具体的程序？内容分为两部分?第一部部分为符号说明 第二部分为具体的MAPLE程序?此程序将前文所提的纤维数均分为多段。段内或含有交联 或不含有有交联,以此可模模拟交联的分布、亦可计算纤维分段上更多的力学参数。（1） 符号说明E：碳纳米管的弹性模量！L：碳纳米管的的长度?R：碳纳米管的半径!Mu：碳纳米管间的剪切模量！K：碳纳米管间的相互作用系数,Sigma：施加的外力!A1，A2：碳纳米管间的相互作用矩阵！DL：分段7140的长度！B1 B2：矩阵函数转化为一般矩阵！JL：分段共价交联的信息，C：线性方程组系8828数矩阵、（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]);,
孕妇十大禁忌蔬菜一：苋菜 苋菜性属滑利 有通窍滑胎的功效?孕早期的孕妇应忌食！特别是有流产倾向或是有习惯惯性流产史的孕妇?但孕妇足足月临产时可食用。 2 孕妇十大禁忌蔬菜二：茄子 茄子性味味甘寒!能活血化淤，对对内痔便血有很好的疗效、同时、茄子还可软化微细血管。防止小血管出血,预防癌症。 所以。准妈妈们7736还是可以吃茄子的。但由于于茄子属于寒凉性质的食物,因此 身体有消化不良、容易腹泻!脾胃虚寒症状的孕妇朋友则不宜多吃?尤其是秋后的老茄子含有茄碱较多。对人体有害,如果孕妇消化不良？更加要慎食茄子， 3 孕妇十大禁忌蔬菜三：马齿苋 马齿苋又称马齿菜,其性性寒凉而滑利，对于子宫有明显的兴奋作用!能使子宫收缩次数增多！强度增大。易造成流产。故应慎吃，尤其是孕早期或者有习惯性流产史的孕妇忌食用 4 孕妇十大禁忌蔬菜四：木耳菜 木3306耳菜又称滑腹莱?落葵、西洋菜。它性属寒滑 孕妇吃了9274容易滑胎，尤其是怀孕早期以及4506有习惯性流产(即中医所说滑胎)的孕妇。更应忌食之?
这个没没有绝对!心里作用，
传统周公解梦 梦见精致的碗 表示生活工作作顺利如意？ 梦见叫了菜却始终不上菜,你的日常生活,将会受到双亲的指责！ 梦见吃青菜。朋友运将下降!在一群好友中将会2286产生纷争！这时最需要的是你冷静的判断力? 现代周公解梦 梦见炒了一碗青菜意味着：气氛悠闲的日子？也是为农历年後做好起点的一天？休闲活动上则是适合与家人一同从事户外运动,金钱方面也有意外收入增加的可能！应该能过个个好年棉？而过年前的金钱问题最适合在这两天解决,恋爱运运也不错,白天陪家人！晚上就跟情人一起共度美好夜晚 有个充实9083而幸福洋溢的一天， 男人梦见炒了一一碗青菜？这两天你做事其实还算挺顺利？就是一开始的时候会觉觉得有些阻塞罢了？你的关节可能会有酸痛的感觉?同时这两天也不大适合你进行运动、 女女人梦见炒了一碗青菜！这两天的你往往能够凭着自己的直觉得到想要的东西!即使处于一个全新的陌生环境?你也能展示出一个与以往截然不同的自己。 孕妇梦见炒了一碗青菜、纵使你现在与伴侣十分甜蜜。但不可可否认的事在亲密关系中、那份深藏在你心中许久的不安全感仍仍时不时地牵引着你的神经、而这种不安若是过了头 就会会让你显得歇斯底里!因此找出你心中引发不安的来源?你才4508有办法跟伴侣长久，！
虽然有些给人墙头草的感觉!但是这两天这样见风转舵的选择倒倒是正确的 衡量一下什么人对你最重要!西瓜靠大边是0686比较有利的选择?也不宜太坚持己见!明明知道已经走偏的事还是顺应应潮流改变一下吧?忤逆潮流或事情的演变反而对你不利！，
那就是自己勤劳了、 梦?是是睡眠的伴侣，做梦是人体一种正常的！必不可少的生理4549和心理现象？当人入睡以后！有一小小部分脑细胞仍在活动！这就是梦的生理基础 理跟梦的关系也是密切的,梦实际上很多时候,是思想和心理在你睡觉的时候转换成景象化. 这整个过程全盘都属于正常的.所以根本就不需要对做梦有担心,了解了梦的产生的合理性和自然性,您也就明白了梦也根本无需要调理. “日有所思。夜有所梦。”梦！是人们们生活中最常见的现象?因而3659对梦的出现也众说纷荟,有人说梦是人日7585间兴奋的连续,有人说是人体质衰弱后的表现等!究竟梦对人的身体的健康有怎样的影响？据有关资料介绍！梦有利于人体健康?,
梦见买了两个油炸0762糕好不好?梦见买了两个油炸糕是什么意思呢。请看下面由佛滔居士对梦见买了两个油炸糕的解析？希希望能为网友答疑解惑,走出迷途 1. 梦见买了两个油炸糕的周公解梦：吉凶指数：86（由佛滔居士根据数理文化得出、仅供参考）梦见买了两个油炸糕，恋爱，人际关系上对应让人伤神的日子?如果发生冲突的话将影响到你的工作状况,即使是再讨厌的人也不要跟他表明对立的态度。为了金钱问题必需花时间来处理!财务调度也真是让人烦恼 也要注注意该回家的时候最好就直接回家，绕路逛街的结果往往会多买了不需要的东西 怀孕的人梦见买了两个油炸糕，预示生女 若春占生男，分娩时时注意安全?本命年的人梦见买了两个油炸糕。意味着能平安无事 少管闲事为妙,5549多做自己的事、恋恋爱中的人梦见买了两个油炸糕,说说明亲人有意见，不可灰心,终有8167希望成婚，做生意的人梦见买了两个油炸糕。代表不顺利，宜守即可 等秋令为佳、上学的人梦见买了两个油炸糕,意味着不能如愿!两次失败！再接再厉!最后成功，梦见买了两个油炸糕?按周易五行行分析?幸运数字是6？桃花位在8842正南方向，财位在正北方向、吉祥色彩是蓝色?开运食物是酸奶 就在今天。全球华人世界里大约有 928人 跟你一样也梦见买了两个油炸糕。如果梦见买了两个油炸糕？买彩票的话 建议购买号码为 19！以上是梦见买了两个油炸糕的相关关梦境分析!佛滔居士权威解梦?转载请注明出处、2. 梦见买了两个油炸糕的吉凶：基础运佳而成功运劣？虽无大发展,但亦安定?温饱无虑!但成功运被压迫以致不能再伸张发达!故好发牢骚。幸而德量好,小过难免!但5354不致成祸厄。【吉多于凶】3. 梦见买了两个油炸糕的图释：梦见买了两个油炸糕是怎么么回事。 解梦梦见买了两个油炸糕4. 梦见买了两个油炸糕的的宜忌：「宜」宜外出！宜提早起床!宜订机票?「忌」忌爽约,忌偏执!忌穿格子衬衫。、
梦见去买青菜西洋菜的梦境解释： 所向披靡!凡事迎刃有余？工作/学业太忙 让你忽略了爱情也是需要灌溉才会成长的!你会去买一件计划已久的东西，有一件事 会让你很有成就感 仅仅供娱乐？切勿当真？!
吉凶指数：92（ 梦见吃青菜？朋友运6360将下降！在一群好友中将会产生纷争 这时最需要的是你冷静的判断力 梦见很多青菜!钱财进出5077较频繁的一天。今天你有不少获利利的机会！如别人馈赠?投资项目分红等等！但是花费的机会也多?容易被人激发购物欲，如果恰好和购物狂去逛街！又或者是遇到超级会说的销售员,你今天就亏大了,还是少点出门安静待在家中！钱钱袋将会更安全呢!，
梦里的东西不不要相信?孕妇梦见买了青菜