ansys命令流error求解? 20
浏览:1208 回答:2
壳模型周围均布三向弹簧约束,求解总是出现 *** WARNING *** CP = 5.148 TIME= 11:15:35 Input time should be greater than zero. *** ERROR *** CP = 5.647 TIME= 11:15:35 The displacements input as real constants must be in ascending order for element 1937. *** ERROR *** CP = 5.647 TIME= 11:15:35 The displacements input as real constants must be in ascending order for element 961. *** ERROR *** CP = 5.647 TIME= 11:15:35 The last input displacement for element 1937 must be positive. 求解答!以下是命令流: /prep7 et,1,shell63 r,1,0.016,0.016,0.016,0.016,,,,, et,2,combin39 keyopt,2,3,1 et,3,combin39 keyopt,3,3,2 et,4,combin39 keyopt,4,3,3 et,5,combin39 keyopt,5,3,3 r,2,-5*e-3,-5*e4,-1.5*e-3,-5*e4,0,0!定义横向弹簧本构 rmore,1.5*e-3,5*e4,5*e-3,5*e4 r,3,-5*e-3,-5*e4,-1.5*e-3,-5*e4,0,0!定义轴向弹簧本构 rmore,1.5*e-3,5*e4,5*e-3,5*e4 r,4,-4*e-3,-6.9*e4,-2.5*e-3,-6.9*e4,0,0!定义垂直弹簧本构 rmore,1.5*e-3,1.667*e4,5*e-3,1.667*e4 r,5,0,0,0.05,600000,0.3,1.4*e6!定义等效弹簧本构 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,dens,1,7800 k,1 k,2,0,0,30 l,1,2 circle,1,0.5 adrag,2,3,4,5,,,1 lsel,r,loc,z,0 lesize,all,,,4 lsel,s,length,,30 lesize,all,,,60 mshape,0,2d mshkey,1 amesh,all NGEN,2,1000,all, , , , , ,1, TYPE, 2 MAT, 1 REAL, 2 ESYS, 0 SECNUM, TSHAP,LINE EINTF,0.0001, ,LOW,0,0,0,0,0 TYPE,3 MAT, 1 REAL,3 ESYS, 0 SECNUM, TSHAP,LINE EINTF,0.0001, ,LOW,0,0,0,0,0 TYPE,4 MAT, 1 REAL,4 ESYS,0 SECNUM, TSHAP,LINE EINTF,0.0001, ,LOW,0,0,0,0,0 alls !两端弹簧 NSEL,S,LOC,Z,0 NSEL,A,LOC,Z,30 TYPE, 5 MAT,1 REAL, 5 ESYS, 0 SECNUM, TSHAP,LINE EINTF,0.0001, ,LOW,0,0,0,0,0 alls NSEL,S, , ,1001,1976,1 D,all,all nsel,s,loc,z,0 NSEL,A,LOC,Z,30 nsel,r,,,1,1000 d,all,ux,,,,,uy,roty,rotz alls /solu antype,trans alphad, betad, timint,off nlgeom,on time,1*e-6 nsubst,4 kbc,1 acel,0,-9.8 solve
是你对 combin39 单元定义的实常数有问题。按照帮助文件,该实常数是 位移-力 曲线,其中的位移应该是升序排列的,而你的实常数定义,比如:
r,2,-5*e-3,-5*e4,-1.5*e-3,-5*e4,0,0 !定义横向弹簧本构
rmore,1.5*e-3,5*e4,5*e-3,5*e4
其中的实数的格式不符合要求,比如:
-5*e-3, -1.5*e-3, ,0
前两个实数中间多了一个 * 号,应该改为:
-5e-3, -1.5e-3
类似的实数都需要修改,包括:
time,1*e-6
需要改为:
time,1e-6
按照原来的命令流,生成的模型中实常数 2 如下:
显然是错误的。
修改后的命令流如下:
! 以下是修改后的命令流:
fini
/clear
/filn,3D_springs,1
/prep7
et,1,shell63
r,1,0.016,0.016,0.016,0.016,,,,,
et,2,combin39
keyopt,2,3,1
et,3,combin39
keyopt,3,3,2
et,4,combin39
keyopt,4,3,3
et,5,combin39
keyopt,5,3,3
r,2,-5e-3,-5e4,-1.5e-3,-5e4,0,0 !定义横向弹簧本构
rmore,1.5e-3,5e4,5e-3,5e4
r,3,-5e-3,-5e4,-1.5e-3,-5e4,0,0 !定义轴向弹簧本构
rmore,1.5e-3,5e4,5e-3,5e4
r,4,-4e-3,-6.9e4,-2.5e-3,-6.9e4,0,0 !定义垂直弹簧本构
rmore,1.5e-3,1.667e4,5e-3,1.667e4
r,5,0,0,0.05,600000,0.3,1.4e6! 定义等效弹簧本构
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7800
k,1
k,2,0,0,30
l,1,2
circle,1,0.5
adrag,2,3,4,5,,,1
lsel,r,loc,z,0
lesize,all,,,4
lsel,s,length,,30
lesize,all,,,60
mshape,0,2d
mshkey,1
amesh,all
NGEN,2,1000,all, , , , , ,1,
TYPE, 2
MAT, 1
REAL, 2
ESYS, 0
SECNUM,
TSHAP,LINE
EINTF,0.0001, ,LOW,0,0,0,0,0
TYPE,3
MAT, 1
REAL,3
ESYS, 0
SECNUM,
TSHAP,LINE
EINTF,0.0001, ,LOW,0,0,0,0,0
TYPE,4
MAT, 1
REAL,4
ESYS,0
SECNUM,
TSHAP,LINE
EINTF,0.0001, ,LOW,0,0,0,0,0
alls !两端弹簧
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,30
TYPE, 5
MAT,1
REAL, 5
ESYS, 0
SECNUM,
TSHAP,LINE
EINTF,0.0001, ,LOW,0,0,0,0,0
alls
NSEL,S, , ,1001,1976,1
D,all,all
nsel,s,loc,z,0
NSEL,A,LOC,Z,30
nsel,r,,,1,1000
d,all,ux,,,,,uy,roty,rotz
alls
/solu
antype,trans
alphad,
betad,
timint,off
nlgeom,on
time,1e-6 ! 原为 1*e-6, 被当作 -6
nsubst,4,10,2 ! 原为 nsubst,4
kbc,1
acel,0,-9.8
solve