这是vibinfo的一个程序,运行有一个错误,望指出,谢谢

liforent · 发表于 2006-5-15 14:55

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

这是vibinfo的一个程序,运行有一个错误,望指出,谢谢! ！这是一个金属容器装满水，用桐木做保护包装，由1.2米处摔下的接触碰撞计算，水被分配到容器壁上，为减少计算时间，从木包装接触地面计算开始，到容器反弹时结束。 !金属容器木质包装的撞击 KEYW,LSDYNA,1 /PREP7 /view,,1,2,3 /plopts,info,1 /pbc,all,1 /pnum,mat,1 /num,1 ET,1,SHELL163 KEYOPT,1,1,6 KEYOPT,1,2,0 KEYOPT,1,3,0 KEYOPT,1,4,0, ET,2,SOLID164 ! --- 输入金属容器材料参数 mp,ex,1,1.86e11 mp,nuxy,1,0.29 mp,dens,1,45435 r,1 rmodif,1,1,,,0.006,,,, ! --- 输入木材材料参数 MP,DENS,2,500 MP,EX,2,1.5e8 MP,NUXY,2,0.2 TB,BISO,2,,,, TBDAT,1,2.5e6 TBDAT,2, ! --- 生成并划分容器(1/8)网格 CSWPLA,11,1,297/597,1, CSYS,11, K,,0.597,,, K,,0.597,-45,, K,,0.597,-90,, L,1,2 L,2,3 CSYS,0, k,,0.597,0.45, k,5000,0,0.1, k,6000,0,-0.1, l,5000,6000 l,1,4 FLST,2,3,4,ORDE,3 FITEM,2,1 FITEM,2,-2 FITEM,2,4 FLST,8,2,3 FITEM,8,6000 FITEM,8,5000 AROTAT,P51X,,,,,,P51X,,90,, LESIZE,1,,,6,,,,,1 LESIZE,2,,,10,,,,,1 LESIZE,4,,,10,,,,,1 LESIZE,5,,,10,,,,,1 LESIZE,6,,,6,,,,,1 LESIZE,7,,,10,,,,,1 LESIZE,8,,,10,,,,,1 LESIZE,9,,,10,,,,,1 LESIZE,10,,,10,,,,,1 amesh,1,3,1 ! --- 生成并划分木盒(1/8)网格 block,0,1,-0.7,0.18,-1,0 CSWPLA,11,1,0.5,1, CSYS,11, K,,0.6,,, K,,0.6,-45,, K,,0.6,-90,, L,16,17 L,17,18 CSYS,0 k,,0.6,0.45, l,16,19 FLST,2,3,4,ORDE,2 FITEM,2,23 FITEM,2,-25 FLST,8,2,3 FITEM,8,6000 FITEM,8,5000 AROTAT,P51X,,,,,,P51X,,90,, FLST,3,3,5,ORDE,2 FITEM,3,10 FITEM,3,-12 VSBA,1,P51X !-------布尔操作的警告信息 VDELE,2,,,1 ! ----为网格划分分割体 K,7000,0.5,-0.8,, K,7001,1.5,-0.25,, !-------创建两条分割线 FLST,2,2,3 FITEM,2,17 FITEM,2,7000 FLST,3,1,8 FITEM,3,0.42,-0.4,0 BSPLIN,P51X,P51X FLST,2,2,3 FITEM,2,17 FITEM,2,7001 FLST,3,1,8 FITEM,3,0.6,-0.28,0 BSPLIN,P51X,P51X ! ----分割体 ADRAG,17,,,,,,29 ADRAG,18,,,,,,29 VSBA,3,5 VSBA,2,7 wprot,,-90 VSBW,4 ! ---- 网格划分 TYPE,2 MAT,2 LESIZE,24,,,10,,,,,1 LESIZE,26,,,10,,,,,1 LESIZE,29,,,10,,,,,1 LESIZE,40,,,10,,,,,1 AMESH,10 VSWEEP,1 ACLEAR,10 LESIZE,28,,,6,,,,,1 LESIZE,47,,,6,,,,,1 LESIZE,23,,,6,,,,,1 LESIZE,35,,,4,,,,,1 VSWEEP,3 VSWEEP,2 VSWEEP,5 LDELE,3,,,1 ALLSEL,ALL NUMMRG,ALL,,,,LOW NUMCMP,ALL APLOT ! --- A正撞---复制单元(还有计算其它的撞法.....) wpoff,,,0.45 CSWPLA,11,0,1,1, CSYS,11, ASEL,S,MAT,,1 ARSYM,X,all,,,,0,0 ARSYM,y,all,,,,0,0 ARSYM,z,all,,,,0,0 VSEL,S,MAT,,2 VSYMM,X,all,,,,0,0 VSYMM,y,all,,,,0,0 ALLSEL,ALL NUMMRG,ALL,,,,LOW NUMCMP,ALL CSYS,0, CSDELE,11,,1, WPSTYLE,,,,,,,,0 EPLOT ! ---- 定义接触组合 esel,s,mat,,1 nsle cm,ncase,node cm,ecase,elem esel,s,mat,,2 nsle cm,nwood,node cm,ewood,elem FLST,5,12,5,ORDE,12 FITEM,5,10 FITEM,5,-11 FITEM,5,14 FITEM,5,48 FITEM,5,54 FITEM,5,64 FITEM,5,67 FITEM,5,73 FITEM,5,82 FITEM,5,86 FITEM,5,91 FITEM,5,98 ASEL,S, , ,P51X NSLA,R,1 CM,target,NODE ALLSEL,ALL EPLOT ! --- 定义容器加速度及初始速度 EDVE,VELO,ncase,0,-4.85,0,0,0,0,,,,,, EDVE,VELO,nwood,0,-4.85,0,0,0,0,,,,,, *DIM,time,ARRAY,2,1,1, , , *DIM,accc,ARRAY,2,1,1, , , *SET,TIME(1,1,1),0 *SET,TIME(2,1,1),1 *SET,ACCC(1,1,1),9.8 *SET,ACCC(2,1,1),9.8 EDLOAD,ADD,ACLY,0,ncase,TIME,ACCC,0,,,,, EDLOAD,ADD,ACLY,0,nwood,TIME,ACCC,0,,,,, ! --- 定义木盒的约束 nsel,s,loc,y,-0.7 d,all,,0,,,,uy nsel,s,loc,Y,-0.7 nsel,r,loc,X,0 nsel,r,loc,Z,-1 d,all,ux,0,,,,, nsel,s,loc,Y,-0.7 nsel,r,loc,X,0 nsel,r,loc,Z,1 d,all,ux,0,,,,, nsel,s,loc,Y,-0.7 nsel,r,loc,X,0 nsel,r,loc,Z,0 d,all,ux,0,,,,,uz ! --- 在容器与木盒之间定义的面面接触 edcgen,asts,ncase,target,0,0,0,0,0, , , , ,0,10000000 edcontact,1.0,,,,,,2,,0 ! --- 为后处理建立单元组合half NSEL,S,LOC,Z,-0.0001,-1 ESLN,S CM,half,ELEM nsel,all esel,all ! --- 计算输出设置 /solu time,0.0045 !计算时间 edrst,50 !输出文件数 edhtime,500 EDDUMP,1, EDENERGY,1,1,1,1 ALLSEL,ALL eplot ! --- 计算 FINISH /SOLU /STAT,SOLU SOLVE

liforent · 发表于 2006-5-15 16:21

[em06]

账号		自动登录	找回密码
密码			我要加入