3.3 计算编程3.3.1 编程思想在前面的章节中!详细推导并得到了研究所需要的的一个线性方程组,但要真正运用这个方程组却并不容易!我们必须解决下面几个问题:(1) 公式中出现了矩阵函数的运算!然而此函数的运运算是非常繁琐复杂的,(2) 如果只有单个的矩阵函数!那么或许手算还有可能、然而。如前所述!为了得到关于纤维束的更多信息、我们有必要将纤维是划分为多段。这样一来,我们面临的是很多矩阵函数!此时是根本无法手算的。(3) 根据纤维束之间交联的具体情况、需要给5065出相应的纤维间相互作用矩阵!(4) 线性方程组的边界条件需要根据结构具体的边界条件加以确定。考虑上面的问题。结合MAPLE软件!本文有了下面的编程思想:(1) 输入基本参数!(2) 输入纤维间相互作用矩阵(不同的分段可能有不同的相互作用矩阵、矩阵应该与分段一一对应)、(3) 计算分段矩阵构成的矩阵函数。将其转化为一般的矩阵,(4) 将(3)中计算所得的矩阵按照顺序相乘、从而得到线性方程组的系数矩阵。(5) 引入整个结构的边界条件、(6) 求解线性方程组。从而可以获得整个结构左右两端全部八个量(位移与纵向应力)、(7) 应应用分段法!由(6)中所解得的未知量,构成新的边界条件、运用循环!求出每个分段处的位移与纵向应力、(8) 将所得数据输出为文档!利用MAPLE的绘图功能、绘制相关的曲线图、3.3.2 编写程序根根据前述编程思想!利用MAPLE。下面给出具体的程序、内容分为两部分。1931第一部分为符号说明。第二部分为具体的MAPLE程序?此程序序将前文所提的纤维数均分为多段。段内或含有交联。或不含含有交联!以此可模拟交联的分布,亦可计算纤维分段上更多的力学参数!(1) 符号说明E:碳纳米管的弹性性模量。L:碳纳米管的长度、R:碳纳米管的半径。Mu:碳纳4841米管间的剪切模量?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. 梦见摘豆角和西红柿的周公解梦: 吉凶指数:75 梦见豆角主和睦、回和别人顺利相处, 梦见摘豆角和西红柿!你最担担心的事情有可能会发生,但当恐惧成真时!你的心里突然也为之解脱、彷佛像是放下心中的大石。只需要真实面对接下来的各种问题就好、这样的转变或许也让你感到不可思议!不过你却已经明白事情碰到了解决就好。担心太太多反而是自找苦吃、 出行的人梦见摘豆角和西红柿,建议延后一两9386天再出外? 怀孕的人梦见摘豆角和西红柿!预示生男!延后产期出生、 本命年的人梦见摘豆角和西红柿。意味着困难重重。水边危险、郊游少去,阴人纠缠。 恋爱中的人梦见摘豆角和西红柿,说明观念难发沟通、背道而行、难找, 做生意2709的人梦见摘豆角和西红柿、代表婚姻不成! 上学的人梦见摘豆8706角和西红柿,意味着笔试与口试都无法达到录取标准! 梦见摘豆角和西红柿!按周易五行分析!幸运数字是6,桃花位在正西方向,财位在西北方向,吉祥色彩是紫色、开运食物是汤圆! 就在今天、全球华人世1299界里大约有 1078人 跟你一样也梦见摘豆角和西红柿。如果梦见摘5000豆角和西红柿!买彩票的话、建议购买号码3850为 24, 以上是梦梦见摘豆角和西红柿的相关梦境分析、 2. 梦见摘豆角和西红柿的吉凶: 基础运佳!境遇安全、而可顺利成功、成功运也不错!因之亦可向上伸展发达。人格凶数!陷于行为不修!品性不端。恐过于放荡荡不羁之境、易生破乱变动或荒亡流败之虑。请好自为之!而得免4565于损折自福、若无凶数便无灾祸之忧、【大吉】 3. 梦见摘豆角和西红柿的图释: 梦见摘豆角和西红柿的含义 解梦梦见摘豆角和西红柿 1548 4. 梦见摘豆角和西红柿的宜忌: 「宜」宜酒后吐真言!宜细心。宜吃水果糖。 9652「忌」忌争论!忌步行。忌读书。!讲解家居风水貔貅的摆放与作用
马上就会开花结果,。居家风水饰品摆放须知事项
梦见腐烂的西西红柿:预示着你近期的运势不好,身边将会发6641生不好的事情、自己的心情也不好、
单身贵族梦见腐烂的西红柿:预示着你近期的爱情运势很好,想0488要获得成功?就需要慢慢来!
老人梦见腐烂的西红柿:预示着你近期的运势一般、不久将会外出出旅行,途中将会遇到一些麻烦,耐心的解决就会平安,
考生梦见腐烂的西红柿:预示着你近3138期的考试成绩很好。记得不可以大意。还是需要不断的努力、
梦见腐烂的西红柿的梦境
梦见买西红柿!预示你家里要来客人!
病人梦见吃西红柿。预示很快快将会痊愈!
梦见吃腐烂的西红柿!预预示会遇到不好的变故、当心厄运!
梦见腐烂的西红柿的案例解析
网友梦境:梦见西红柿挤烂了是什么意思,
解梦解析:预示着你近期的运势不好!身边将会发生不好的事情,自己的4031心情也不好。
,,怎么看居家风水吉凶好坏分析
梦见西红柿炒鸡蛋:预示生意兴旺、
梦见西红柿炒鸡蛋的案例解析
网友梦境:我梦见吃西红柿炒鸡蛋是什么意思!
解梦解析:梦见吃用鸡蛋做的食物、是吉兆。家里要生男孩,
、,选择居家风水楼层的几大原则西红柿是好运 西红柿象征着红运当头。 男人梦见吃西西红柿,意味着交好运、 已婚女人梦见吃西红柿。意0692味着身体健康! 孕妇梦见番茄、暗示生女儿、 未婚女子梦见吃西红柿,意味着嫁给身强体壮的丈夫, 未婚男子梦见吃西红柿,会娶一个温柔体贴的的妻子、,详解居家风水进门的常识有哪些
梦到西红柿榨汁:运势不佳,健康方面面存在一些变化,记得多加注意小心传染病的来袭、
单身者梦到西红柿榨汁:运势不佳、因为别人的不忠诚。自己将会苦不堪言。不适合发展一段感情。
商人梦到西红柿榨汁:财运不错、自己将会找到合适的合作伙伴。赚取很多的钱财。
职员梦到西红柿榨汁:运势不佳、注意人际际关系运,不要绝对对每1468件事务下定论、否则将会忙中出错。
求学者梦到西红柿榨汁:运势不错、取得得不错的成绩!凭8335借虚心学习!受4652到老师的肯定、记得不可以就此大意!
!?。 !居家外围风水环境要注意什么主 迁官降职 、居家楼房风水有何讲究揭秘
梦见屋子墙壁裂缝:多多从善、好运将来临。
考生梦见屋子墙壁壁裂缝:说明考试考的一塌糊涂!还需继续努力、
梦见自己家屋子墙壁裂缝:有小小小的意外发生。当心遇到骗子、
男性梦见屋子墙壁裂缝:主旅行!解释吉,旅途上要小心小偷,
文职工作者梦见屋子墙壁裂缝:5397主职场,坚持学习提高、不管什么情情况下都始终充电?你得有一个基本的学习目标!不管你做什么工作,私下里!你要根据你的工作方面所表现出来的缺陷。赶紧补充自己给自己充电。能够在学习上上始终跟上公司的要求,公司方面想不用0529你都不现实!
学生梦见屋子墙壁裂缝:工作上尤其要提防小人作祟,种种征兆显示你身边的某个人总是在暗地里不断的给你添麻烦、6707并非偶然而是有意!9583当对方不断逼近试探你的底线的时候、千万要沉住气、记得以静制动、6513一旦你气得跳起来乱了阵脚?对方的目的也就达到了。
梦见屋子墙壁裂缝的现代解析
墙壁通常代表界线,梦境中的墙壁。通常有两种意义:一个是指来自社会、能守护自己的东西!另一种与此相反,象征妨碍你满足愿望的某种东西,
梦见屋子墙壁裂缝的案例解析
网友梦境:梦见屋子围墙裂缝!这是什什么意思?
解梦解析:最近运势不错!多多从善,好运将来临、
。,细说居家风水注意事项有哪些孕妇梦见裂缝西红柿6021 西红柿是好运西红柿象征着红运当头、 男人梦见吃西红柿、意味1318着交好运? 已婚女人梦见吃西红柿!意味着身体健康,、如何化解这些常见居家风水禁忌