yuqiangliu 发表于 2005-10-7 10:04

这个例子那里表明层间的不同角度?

<P>!采用shell99单元(线性层合单元),模拟层合板受单向拉伸,拉伸应力100MPa<BR>!层合板尺寸:长90mm,宽25mm,厚1mm</P>
<P>/prep7<BR>ET,1,SHELL99,,,,,2,4   <BR>KEYOPT,1,8,1<BR>R,1,8,1   <BR>RMODIF,1,13,1,0,0.125,1,45,0.125   <BR>RMODIF,1,19,1,-45,0.125,1,90,0.125 <BR>MP,EX,1,135000 <BR>MP,EY,1,8000   <BR>MP,EZ,1,8000   <BR>MP,PRXY,1,0.34 <BR>MP,PRYZ,1,0.34 <BR>MP,PRXZ,1,0.34 <BR>MP,GXY,1,4500<BR>MP,GYZ,1,3700<BR>MP,GXZ,1,4500</P>
<P>!建立层合板模型<BR>rectng,-45,45,-12.5,12.5</P>
<P>!定义网格划分<BR>LSEL,S,LENGTH,,90    !选择所有长度为90mm的线段(长度方向上的线段)<BR>LESIZE,all, , ,10, , , , ,1   !长度方向等长划分10个单元<BR>lsel,s,length,,25   !选择所有长度为25mm的线段(宽度方向上的线段)<BR>lesize,all, , ,20,-10, , , ,1!宽度方向划分20个单元,中间疏,两边密</P>
<P>!映射划分层合板<BR>MSHAPE,0,3D <BR>MSHKEY,1<BR>AMESH,all</P>
<P>!定义约束<BR>/solu<BR>nsel,s,loc,x,0!选择所有x=0的节点<BR>D,ALL, , , , , ,UX !给选中的节点施加x方向的约束(无x方向的位移)<BR>nsel,s,loc,y,0!选择所有y=0的节点<BR>D,ALL, , , , , ,UY !给选中的节点施加y方向的约束(无y方向的位移)<BR>allsel,all</P>
<P><BR>!施加拉伸载荷<BR>lsel,s,loc,x,45<BR>lsel,a,loc,x,-45<BR>SFL,ALL,PRES,-100,!施加100MPa的拉应力(pressure)</P>
<P>!求解设置<BR>allsel,all   !求解前选择所有元素<BR>/solu<BR>solve</P>
<P>!计算结果分析(因为是shell单元,故只能计算面内应力,不能考虑层间应力,在边界处没有边缘效应)<BR>/post1<BR>layer,1   !取第一层(从下往上),即0度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,2   !取第二层(从下往上),即45度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,3   !取第三层(从下往上),即-45度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,4   !取第四层(从下往上),即90度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,5   !取第五层(从下往上),即90度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,6   !取第六层(从下往上),即-45度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,7   !取第七层(从下往上),即45度层结果<BR>PRNSOL,S,COMP!list 结果<BR>layer,8   !取第八层(从下往上),即0度层结果<BR>PRNSOL,S,COMP!list 结果<BR>!注:以上所得结果与理论解完全符合</P>
<P><BR>我看了这个例子 看不懂rmodif 的意思,那里表面了0/45/-45/90/90/-45/45/0<BR>谢谢</P>

zp4655958 发表于 2011-10-27 19:13

命令流太乱了
页: [1]
查看完整版本: 这个例子那里表明层间的不同角度?