我写了个命令流 可以算的 应力也有的 没具体看 参数随便编的 你自己改一下好了
/filename,Q1
/title,JD
/prep7
et,1,solid65,,,,,,,1
mp,ex,1,24320
mp,prxy,1,0.2
mp,dens,1,20e-6
fc=25.6
ft=2.6
qc=7.8
qt=0.28
tb,concr,1
tbdata,,0.5,0.95,ft,-1
tb,mkin,1,
tbtemp,,strain
tbdata,,0.0002,0.0008,0.0014,0.002,0.0033
tbtemp,0.0
tbdata,,0.19*fc,0.64*fc,0.91*fc,fc,0.9*fc
blc4,,,500,300,500
elemsiz=100
vsel,s,,,1
vatt,1,1,1
mshkey,1
esize,elemsiz
vmesh,all
nsel,s,loc,y,0
d,all,all
allsel,all
*SET,NT,50
*DIM,TIME,,NT
*DIM,AC_X,,NT
*CREATE,ansuitmp
*VREAD,TIME,'TIME','txt',' ',1, , , , , ,
(F6.2)
*END
/INPUT,ansuitmp
*CREATE,ansuitmp
*VREAD,AC_X,'AC_X','txt',' ',1, , , , , ,
(F6.2)
*END
/INPUT,ansuitmp
FINISH
/SOL
ANTYPE,TRANS
nlgeom,off
nropt,full,,on
ncnv,1
neqit,50
lnsrch,on
pred,on
autots,on
nsubst,1,
*DO,K,1,NT
ACEL,AC_X(K),AC_X(K),0
TIME,TIME(K)
SOLVE
*ENDDO
FINISH |