埋地管道动力时程分析没有问题和警告,但应力应变结果很小,位移基本没有怎么回事? 300
埋地管道动力时程分析没有问题和警告,计算结果位移非常小,应力也很小,感觉不太对,我用的是m,pa为单位,有大神能帮忙检查下,看是什么问题不,非常感谢,以下是完整的命令流,附件为地震数据文件。
/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,15600
k,1
k,2,0,0,15
l,1,2
circle,1,0.25
adrag,2,3,4,5,,,1
lsel,r,loc,z,0
lesize,all,,,4
lsel,s,length,,15
lesize,all,,,25
mshape,0,2d
mshkey,1
amesh,all
NGEN,2,10000,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,15
MAT,1
REAL, 5
ESYS, 0
SECNUM,
TSHAP,LINE
EINTF,0.0001, ,LOW,0,0,0,0,0
alls
NSEL,S, , ,10001,15000,1
D,all,all
nsel,s,loc,z,0
NSEL,A,LOC,Z,15
nsel,r,,,1,5000
d,all,ux,,,,,uy,roty,rotz
alls
*SET,NT,1000
*SET,DT,0.02
*DIM,AC,,NT,1
/solu
antype,trans
alphad,
betad,
timint,off
nlgeom,on
time,1e-6
nsubst,4
kbc,1
acel,0,9.8
solve
timint,on
NROPT,FULL,,
CNVTOL,F,,0.05,,0.01
CNVTOL,U,,0.05,,0.01
/input,dizhengbo,txt
*do,i,1,NT
acel,0,AC(i)*9.8,0
time,i*DT
nsubst,4
solve
*enddo