| 查看: 2116 | 回复: 1 | ||
[求助]
跪求大神看看这个命令流哪里有问题，好人一生平安！！！！ 已有1人参与
|
|
命令流如下： finish /clear /PREP7 !* ET,1,SOLID45 ET,2,SOLID65 et,3,SOLID65 !* MP,EX,1,4800 MP,EY,1,4800 MP,EZ,1,51500 MP,PRXY,1,0.02 MP,PRYZ,1,0.02 MP,PRXZ,1,0.27 MP,GXY,1,2353 MP,GYZ,1,2353 MP,GXZ,1,20275 mp,ex,2,20000 mp,nuxy,2,0.25 tb,concr,2 tbdata,,,,2,20 mp,ex,3,1.5e4 mp,nuxy,3,0.2 tb,dp,3 tbdata,1,4,30,10 r,1$r,2$r,3 ! CYLIND,16,0,-300,1500,0,90, CYLIND,50,16,0,1500,0,90, CYLIND,500,50,0,1500,0,90, ! allsel,all vsel,s,,,1 aslv,s lsla,s lplot lesize,1,,,4 lesize,2,,,4 lesize,3,,,6 lesize,4,,,6 lesize,5,,,4 lesize,6,,,4 lesize,7,25 lesize,8,25 lesize,9,25 type,1 mat,1 real,1 mshape,0,3d mshkey,0 vsweep,all !* ! allsel,all vsel,s,,,2 aslv,s lsla,s lplot lesize,10,,,4 lesize,11,,,6 lesize,12,,,4 lesize,13,,,6 lesize,14,,,6 lesize,15,,,4 lesize,16,,,6 lesize,17,,,4 lesize,18,25 lesize,19,25 lesize,20,25 lesize,21,25 type,2 mat,2 real,2 mshape,0,3d mshkey,0 vsweep,all !* ! allsel,all vsel,s,,,3 aslv,s lsla,s lplot lesize,23,,,6 lesize,25,,,6 lesize,26,,,6 lesize,28,,,6 lesize,22,50 lesize,24,50 lesize,27,50 lesize,29,50 lesize,30,50 lesize,31,50 lesize,32,50 lesize,33,50 type,3 mat,3 real,3 mshape,0,3d mshkey,0 vsweep,all !* allsel,all ! save ! ! mp,mu,4,0.4 ! 修改 mat,4 r,4,4 real,4 et,4,170 et,5,174 KEYOPT,5,9,0 KEYOPT,5,10,2 nropt,unsym ! 不对称矩阵完全牛顿—拉普森法 !*generate the target sunface allsel,all asel,s,,,9 cm,target1,area type,4 nsla,s,1 esln,s,0 esurf !*generate the contact sunface asel,s,,,3 cm,contact1,area type,5 nsla,s,1 esln,s,0 esurf ! allsel,all eplot ! !*CONTACT PAIR CREATION - END mp,mu,5,0.7 ! 修改 mat,5 r,5,4 real,5 et,6,170 et,7,174 KEYOPT,7,9,0 KEYOPT,7,10,2 nropt,unsym !*generate the target sunface allsel,all asel,s,,,15 cm,target2,area type,6 nsla,s,1 esln,s,0 esurf !*generate the contact sunface asel,s,,,8 cm,contact2,area type,7 nsla,s,1 esln,s,0 esurf ! allsel,all eplot !*CONTACT PAIR CREATION - END save ! !* /solu da,4,symm !对称约束 da,5,symm da,10,symm da,11,symm da,16,symm da,17,symm da,13,all da,14,all da,7,all ! antype,0 nlgeom,1 !大变形 autots,on !自动时间步 neqit,50 !子步最大迭代次数 pred,on !非线性分析预测器（加速收敛） ! solcontrol,on,on ! nropt,full,,on !定义H法选项（270） nsubst,20,10000,2 !子步数 time,100 !时间 allsel,all solve !* nsubst,100,10000,10 outres,all,all autots,1 time,200 asel,s,loc,z,-300 sfa,all,,pres,-10 !修改荷载值 allsel,all solve 错误提示为： *** ERROR *** CP = 51.231 TIME= 12:26:52 The value of UZ at node 67 is 1.342765716E+09. It is greater than the current limit of 1000000. This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model is properly constrained. *** ERROR *** CP = 51.387 TIME= 12:30:06 *** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION *** If one or more parts of the model are held together only by contact verify that the contact surfaces are closed. You can check contact status in the SOLUTION module for the converged solutions using CNCHECK. *** ERROR *** CP = 51.480 TIME= 12:30:22 *** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION *** Rigid body motion can also occur when net section yielding has occurred resulting in large displacements for small increments of load or when buckling has occurred. You can plot the time history curve for node 67 in the UZ direction to check for stiffness (slope of the curve) approaching zero. |
回复此楼
2楼2016-01-10 23:00:02
50