最近做课题,用压电做为激发器,给混凝土梁激励,然后看梁的响应!命令流编好后计算时候不收敛,不知哪里出错了,有可能是网格划分问题,有可能是激励有问题。命令流列出来大家帮忙找错!
/PREP7
ET,1,SOLID185
et,2,solid5
L=0.100
W=0.100
DEPTH=0.400
l1=0.015
d1=0.001
BLC4,0,0,L,W,DEPTH
BLOCK,0.0425,0.0575,0.0425,0.0575,0.05,0.051
vptn,1,2
VGLUE,2,3
MP,EX,1,30E9
MP,NUXY,1,0.2
MP,DENS,1,2.5E3
mp,damp,1,0.05
mp,dens,2,7600
mp,perx,2,7.124
mp,pery,2,7.124
mp,perz,2,5.841
tb,PIEZ,2
TBDATA,16,10.5
TBDATA,14,10.5
TBDATA,3,-4.1
TbDATA,6,-4.1
TbDATA,9,14.1
TB,ANEL,2
TBDATA,1,12.8E10,6.8E10,6.6E10
TBDATA,7,12.8E10,6.6E10
TBDATA,12,11E10
TBDATA,16,2.1E10
TBDATA,19,2.1E10
TBDATA,21,2.1E10
vsel,s,,,3
vatt,1,,1
vsel,s,,,2
vatt,2,,2
allsel
wpoffs,0,0,0.05
vsbw,all
wpcsys
wpoffs,0,0,0.051
vsbw,all
wpcsys
wprota,0,90,0
wpoffs,0,0,-0.0425
vsbw,all
wpoffs,0,0,-0.015
vsbw,all
wpcsys
wprota,0,0,90
wpoffs,0,0,0.0425
vsbw,all
wpoffs,0,0,0.015
vsbw,all
wpcsys
allsel
vglue,all
ESIZE,0.010
vmesh,all
DA,114,ALL,0
DA,16,ALL,0
DA,119,ALL,0
DA,124,ALL,0
DA,49,ALL,0
DA,128,ALL,0
DA,70,ALL,0
DA,36,ALL,0
DA,80,ALL,0
WPSTYLE,,,,,,,,0
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'sin2'
*SET,_FNCCSYS,0
*DIM,%_FNCNAME%,TABLE,6,10,1,,,,%_FNCCSYS%
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 100, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 10, 0, 0, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, 99, 0, 1, -3, 0, 0
FINISH
/SOLU
M=10
!N1=NODE(0.050,0.0500,0.10)
ANTYPE,4
TRNOPT,FULL
TINTP,'',0.25,0.5
KBC,0
*DO,I,1,M,1
NSUBST,10
OUTRES,NSOL,ALL
AUTOTS,off
D,42,VOLT,%sin2%
allsel
TIME,0.01*i
SOLVE
*ENDDO
FINI
/post26
NUMVAR,30
N=400
n2=node(0.050,0.0500,0.10)
n3=node(0.050,0.0500,0.2)
n4=node(0.050,0.050,0.250)
n5=node(0.050,0.050,0.30)
n6=node(0.050,0.050,0.350)
nsol,2,N2,U,z,100
nsol,3,n3,u,z,200
nsol,4,n4,u,z,250
nsol,5,n5,u,z,300
nsol,6,n6,u,z,350
plvar,2,3,4,5,6
|