俩个性签名啊啊呀哈哈哈哈哈哈!
3.3 计算编程3.3.1 编程思想在前面的章节中!详细推导并得到了研究所需要的一个线性方程组、但要真正运用这个方程组却并并不容易,我们必须解决下面几个问题:(1) 公式中出现了矩阵函数的运算。然而此函数的运算是非常繁琐复杂的。(2) 如果只有单个的矩阵函数、那么或许手算还有可能!然而,如前所述!为了得到关于纤维束的更更多信息,我们有必要将纤维是划分为多段,这样一来、我们面临的是是很多矩阵函数。此时是是根本无法手算的。(3) 根据纤维束之间交联的具体情况。需要给出相应5873的纤维间相互作用矩阵,(4) 线性方程组的边界条件需要根据结构具体的边界条件加以确定!考3158虑上面的问题,结合MAPLE软件。本文有了下面的编程思想:(1) 输入基本参数。(2) 输入纤维间相互作用矩阵(不同的分段可能有不同的相互作用矩阵、矩阵应该与分段一一对应)、(3) 计算分段矩阵构成的矩阵函数!将其转化为一般般的矩阵?(4) 将(3)中计算所得的矩阵按照顺序相乘!从而得到线性方程组的系数矩阵。(5) 引入整个结构的边界条件。(6) 求解线性方程组,从而可以获得整个结构左右两端全部八个量(位移与纵向应力)。(7) 应用分段法,由(6)中中所解得的未知量,构成新新的边界条件!运用循环,求出每个分段处的位移与纵向应力、(8) 将所得数据输出为文档、利用MAPLE的绘图功能!绘制相关的曲曲线图、3.3.2 编写程序根据前述编编程思想,利用MAPLE、下面给出具体体的程序、内内容分为两部分,第一部分为符号说明!第二部分为具体的MAPLE程序、此程序将前3610文所提的纤维数均分为多段!段内或含有交联、或不含有交联、以此可模拟交联的分布。亦可计算纤维分段上更多的力学参数。(1) 符号说明E:碳纳米管的弹性模量?L:碳纳米管的长度、R:碳纳米管的半径、Mu:碳纳米管间的剪切模量、K:碳纳米管间的相互作用系数。Sigma:施加的外力、A1。A2:碳纳米管间的相互作用矩阵、DL:分段的长度,B1。B2:矩阵函数转化为一般矩阵、JL:分段共价交联的信息!C:线性方4029程组系数矩阵。(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]);!
梦见被人拿砖头打:要关心家人的健康、
6580梦见被砖头打:会幸福、富有。
女性梦见被人拿砖头打主旅行,解释:一路平安、
生意人梦见被砖头打主钱财方面:财运状况起落较大。能够凭借自己的技能有不错的收入。但是意3097外的花销也多!投资有损失的可能,
梦见别人用砖头砸自己:试图有一番作为与改变!但是压力依旧沉重,保持努力不懈才2660有机会突破现状!工作上会有新7969的机会出现?尽显现出专业形象的包装!依然是你最大的优势
学子梦见别人用砖头砸自己说明考试成绩好、
未婚者梦见别人用砖头砸自己解析:您的恋情勇敢积极进行!便可成功,
。!孩子缺少关怀啊、。
梦到被砖头砸:运势不错!没有想到出门散步的时候。遇上好久未见的朋友,
求学者梦到被砖头砸:运势不错!学习过程中得到帮助、成绩则是是一路飙升、
男人梦到被砖头砸:运势不佳?出门游玩的时候,将会遇到很多的意外、
未婚者梦到被砖头砸:运势不错,双方之间的的感情很好,6745有进一步的发展、
生意人梦到被砖头砸:运势不佳,感觉现在手上的生意不赚钱。可以试试换一种!
、、。 ,梦见黑色的蛇被我用砖头砸死了意味着: 和另一半出去散心,却没想被一路唠叨--叔父病了、家里的狗狗又和邻居家的打架了…心里里烦死了、烦死了,期盼了好久的休息就被这样破坏了,还好结识了另外一位帅哥美女,聊得很开心。总算有一点慰籍!但是可要要小心啰,脚踏两只船可是很危险哦。 梦见黑色的蛇被我用砖头砸死了的吉凶: 基础安定,能逃灾害、排除障碍而达达到成功,不过却因成功运不吉,不能再伸展,不能5364将天赋才能以完全发展,使其成就与精华俱皆受到打折!殊为可惜,又又须深戒提防色难之忧!【中吉】、
个人认为不需要如此迷信封建!,
现代心理学解释梦见别人用砖头砸自己
梦见别人用砖头砸自己?试图有一番作为与改变!但是压力依旧沉重、保持努力不懈才有机会突破现状!工作上会有新的机会出现,尽显现出专业形象的包装。依依然是你最大的优势,
学子梦见别人用砖头砸自己。说明考试成绩好!
未婚者梦见别人用砖头砸自己!您的恋情勇敢积积极进行!便可成功!
梦见砖厂正在烧砖、幸福生活会到来!
梦见烧好的砖头、祥瑞,会有好运气,
梦见未烧好的的泥砖头、是前途危险或前途凶险的兆头、
梦见砖头、居住在异国的亲人人或亲戚即将登门!
生意人梦见砖头?预预示财运不佳!从事饮食业。公公共事业还可以!
老人梦见砖头。预示身身体健康!儿孙满堂、
未婚者梦见砖头、预示爱情不顺,
待考者梦见砖头。预兆考试成绩好。
梦见买砖头!预示着你近期0665的人际关系不错!和周围人处得很好。如果以前有过节的!主动问候下!你们的关系会得到修复!
梦见搬砖头。预示着你近近期很有责任感,能很好的履行自己应付的责任,是吉兆!
梦见搬红砖,预示着你近期可能是工作压力大!使得你精神上过于于紧张、不愿意意和他人接触。即使是有亲密的朋友来关心你、你也不愿意接受,
梦见自己一直在搬砖头。预示着你近期的财运不错、有意外之财!是吉兆,
梦见别人搬砖头!预示着你近期的运势一般。但你的行3717为举止会成为大家的笑柄!建议你要做好心理准备!
梦见用砖头头打人!预示着你近期可能会大难临头,是不祥之兆。
梦见别人用砖头打自己!预示着近期你的名声大作!从而使得你你的经济水平得以提高?是吉兆。
梦见7094红砖头。预示着你近期的运势不是很好、目前有很多待解决的问题和压力、建议你要做好心理准备、
梦见别人用砖头砸自己的案例分析
网友梦境:梦里面突然有人冲过来!用砖头砸我!然后我就醒了!这是什么意思,
解梦解析:梦见别人用砖头砸自己!试图有一番作为与改变!但是压力依旧沉重、保持努力不懈才有机会突破现状、工作上会有新的机会出现,尽显现出专业形象的包装、依然是你最大的优势。
、,煞梦。
梦见用砖头垒墙生火的梦境解释: 自信信坚实的气息,对年後的目标颇能掌握得到。尤其在金钱方2947面的规划有不错的讯息入耳的暗示喔!而这两天还是努力成果显现的收获日!有有机会参加团体活动的人要为别人多付出喔!需要跑腿搬东西。生火烤肉等杂事不推给别人、表现勤快的你即使忙得团团转但在异性眼中却是特别亮眼、 梦1012见用砖头垒墙生火的吉凶: 基础坚实?心身安泰,奴上位之引进而顺调成功发展!能享得幸福。长寿。圆满、【大吉昌。】。