马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
<P>!采用shell99单元(线性层合单元),模拟[0/45/-45/90/90/-45/45/0]层合板受单向拉伸,拉伸应力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> |