完全拉格朗日-二维-大变形梁-理论开发.docVIP

完全的拉格朗日-二维-大变形梁 假设 1.1.适用范围：大位移、大转动、小应变、线弹性 1.2.模型：铁木辛柯 示意图 基本公式 3.1.沿变形前梁轴线建立局部坐标系，轴线上点的位移为 3.2.轴向应变与轴力 取变形前轴线上一段矢量微元，变形后变为 定义伸长比为 轴向应变为 轴力为 其中分别为弹性模量、截面积 3.3.曲率与弯矩 设截面转角为 曲率为 弯矩为 其中为截面惯性矩 3.4.剪应变与剪力 设轴线转角为，根据轴线微元变形前后的位形，可得与之间的关系 剪应变为 截面剪力为 3.5.关于铁木辛柯模型与伯努利模型的选择 前面是铁木辛柯模型的一般描述。而伯努利模型假设，变形后截面与轴线仍然垂直，即截面的剪应变，因此截面转角与轴线转角相等，即。对于经典小变形问题，可以把转角与位移的关系化简为，然后采用厄米特插值形函数来处理。但大变形问题不能如此化简，所以导不出相应的形函数来引入上述关系式。所以我们采用铁木辛柯理论，将转角与位移看作独立的未知量，进行独立插值。 3.6.关于单元阶数的选择 既然确定了转角与位移独立插值，单元阶数又出现了一个问题。如果选择2节点线单元，那么位移的形函数就线性的，这与实际情况相差太远。所以我们应该选用3节点以上的单元，而选用过高的阶数，又难免数值上的龙格震荡现象。实践表明，3节点单元就是一个不错的选择 弱解形式及其线性化 4.1.根据虚位移原理，内力在虚应变上做的功，等于外力在虚位移上做的功 其中分别为分布力，集中力9，集中力矩 4.2.虚应变展开 记 4.2.1.轴向虚应变 其中 同理 所以 4.2.2.剪切虚应变 不能做独立未知量，应由表示。但关系式涉及三角函数，可考虑将也转换为三角函数形式。由于我们假设小应变，所以 展开得 又由于小应变，所以上式中的分母可替换为 这样就用位移和转角简洁表达出了剪应变 在小应变的假设下，上式并不会带来明显的误差 其虚应变为 记，代入得 4.3.下面进行线性化 原弱解形式记为 其牛顿切线法格式为（参见文档《非线性问题的牛顿切线法.doc》） 其中为线性化后的应变，为位移梯度 4.3.1.正应变部分 左端应变线性化为 右端应变相减为 其中 因此 代入得 4.3.2.剪应变部分 左端应变线性化为 右端应变相减为 4.3.3.代入方程，得到最终的弱解形式 脚本 5.1.nbmull3.ges nbmull3 defi disp u,v,st, var u1 v1 st1 u2 v2 st2 u3 v3 st3 refc rx, coor x, coef un,vn,stn, func epsl,kapa,gama, dord 1,1,1, node 3 mate pe pv pa piz fu fv fm \ 200e9 0.3 0.01 1.e-5 0.0 0.0 0.0 $c6 pg=pe/(2*(1+pv)) shap u= u1=rx*(rx-1.)/2 u3=(1.-rx**2) u2=rx*(1.+rx)/2 v= v1=rx*(rx-1.)/2 v3=(1.-rx**2) v2=rx*(1.+rx)/2 st= st1=rx*(rx-1.)/2 st3=(1.-rx**2) st2=rx*(1.+rx)/2 tran x= x(1)=rx*(rx-1.)/2 x(3)=(1.-rx**2) x(2)=rx*(1.+rx)/2 coef un= un(1)=rx*(rx-1.)/2 un(3)=(1.-rx**2) un(2)=rx*(1.+rx)/2 vn= vn(1)=rx*(rx-1.)/2 vn(3)=(1.-rx**2) vn(2)=rx*(1.+rx)/2 stn= stn(1)=rx*(rx-1.)/2 stn(3)=(1.-rx**2) stn(2)=rx*(1.+rx)/2 gaus = 3 -sqrt(3.0/5.0) 5.0/9.0 0.0 8.0/9.0 sqrt(3.0/5.0) 5.0/9.0 func $cv eux={un/x} $cv evx={vn/x} $cv yita=((1+eux)**2+evx**2)**0.5 $cv cosf=(1+eux)/yita $cv sinf=evx/yita $cv coss=cos(stn) $cv sins=sin(stn) $cv yita2=(1+eux)*coss+evx*sins epsl=+[u/x]*cosf+[v/x]*sinf kapa=+[st/x] gama=+[v/x]*coss-[u/x]*sins-[st]*yita2 stif dist= +[epsl;epsl]*pe*pa +[kapa;kapa]*pe*piz +[gama