关于同时具有接触非线性和材料非线性的子弹冲击问题

panwang

/title,copper bullet impacting a rigid wall <br>/prep7 !进入前处理 <br>et,1,visco106 !选择单元，visco106 <br>et, 2,contac26 !选择二维点面接触单元 <br>keyopt,1,3,1 !轴对称问题 <br>!下面定义所需的材料属性 <br>mp,ex,1,117e9 !材料弹性模量 <br>mp,dens,1,8930 !材料密度 <br>mp,nuxy,1,0.35 !泊松比 <br>r,2,2e15 ! 定义刚性壁常数 <br>!材料属性定义完毕 <br>!下面建立有限元网格 <br>csys,1 ！建立柱坐标 <br>k,1,0,0 <br>k,2,0.002,0 <br>k,3,0.004,0 <br>k,4,0.004,-45 <br>k,5,0.004,-90 <br>k,6,0.0025,-90 <br>k,7,0.0025,-45 <br>l,3,4 <br>l,4,5 <br>csys,0 <br>a,2,3,4,7 <br>a,7,4,5,6 <br>a,1,6,7,2 <br>k,8,0,0.010 <br>k,9,0.0025,0.010 <br>k,10,0.004,0.010 <br>a,1,2,9,8 <br>a,2,3,10,9 <br>type, 1 <br>real,1 <br>esize,,3 <br>amesh,all <br>save <br>！生成接触单元 <br>n,1001,-2*0.004,-0.004 <br>n,1002,2*0.004,-0.004 <br>csys,1 <br>type,2 <br>real,2 <br>*do,j,1,7 <br>nodc=node(0.004,-90* (j-1)/6,0) <br>e,nodc,1001,1002 <br>*enddo <br>csys,0 <br>!下面定义约束 <br>d,all,,0,,,,uz <br>nsel,s,loc,x,0 <br>d,all,,0,,,,ux <br>allsel <br>!约束定义完毕，设置求解选项 <br>tb,biso,1,1,,0 !设置Bilinear 应力应变曲线 <br>tbmodif,2,1,4e8 !屈服强度 <br>tbmodif,3,1,1e8 !材料的剪切模量 <br>tbplot,biso,1 !显示材料应力应变曲线 <br>fini !退出前处理 <br>/solu !进入求解状态 <br>solcontrol,1 !求解控制设为on <br>antype,4 !瞬态分析 <br>cnvtol,f,,,,-1 <br>trnopt,full !Full <br>lumpm,o <br>nlgeom,1 !考虑大变形因素 <br>ic,all,uy,,-227 !定义初始速度 <br>time,1.1e-5 !时间长度 <br>autots,1 !打开自动时间步长设定 <br>deltim,1e-7 !时间步长增量 <br>kbc,1 <br>outres,all,4 !输出所有子步结果 <br>solve !求解 <br>fini !退出求解 <br>/post1 !进入通用后处理 <br>set,last !导入最后一步计算结果 <br>!下面定义参数提取节点变形量 <br>top_node=29 <br>*get,deform,node,top_node ,u,y <br>!下面显示应变 <br>/dscale,1,1 <br>pldisp,2 <br>plnsol,epto,eqv <br>fini <br>!退出通用后处理，进入时间-历程后处理，显示速度时间变化曲线 <br>/post26 <br>nsol,2,top_node,u,y,displace <br>deriv,3,2,1,,velocity,,,1 <br>/axlab,y,velocity <br>plvar,3 <br>fini !分析结束

[此贴子已经被VibInfo于2006-1-10 20:47:37编辑过]

lxy0910

xiexiele[em07]谢谢