本帖最后由 proejhy 于 2017-10-25 15:59 编辑
!请问该例子中如何考虑自重变形下的响应,重力直接加入acel,0,9.8,0不能起作用,也就是预变形不能进行
FINISH$/CLEAR$/PREP7
LSP=20$B=0.3$H=0.6$V=360/3.6$GRA=9.81
M1=6000$EL=0.2$ET,1,BEAM189
ET,2,MASS21,,,4
ET,3,CONTA175
ET,4,TARGE169
KEYOPT,3,12,2
MP,EX,1,2.1E11$MP,PRXY,1,0.3$MP,DENS,1,7800
R,1$R,2,M1$R,3,,,0.1$SECTYPE,1,BEAM,RECT$SECDATA,B,H
K,1$K,2,LSP$K,3,LSP/2,1$L,1,2$LATT,1,1,1,,,3,1$ESIZE,EL$LMESH,ALL
*GET,NMAX,NODE,,NUM,MAX
NMAX=NMAX+1
N,NMAX
TYPE,2$REAL,2$E,NMAX
TYPE,3$REAL,3$E,NMAX$TYPE,4$ESEL,S,TYPE,,1$ESURF$ALLSEL,ALL
D,1,UX,,,,,UY$D,2,UY$D,ALL,UZ$FINISH
ACEL,,9.8,
/SOLU$ANTYPE,TRANS$OUTRES,ALL,LAST$AUTOTS,ON$NLGEOM,ON
TIMINT,OFF$TIME,0.05$NSUBST,10$KBC,0$D,NMAX,UX
F,NMAX,FY,-M1*GRA$SOLVE
TIME,0.1$SOLVE
TIMINT,ON$KBC,0$NSUBST,5
TI0=0.0005
*DO,TI,TI0,LSP/V,TI0
TIME,TI+0.1
D,NMAX,UX,TI*V
SOLVE$*ENDDO
/POST26$N1=NODE(LSP/2,0,0)$NSOL,2,N1,U,Y$NSOL,3,NMAX,U,Y$PLVAR,2,3
NSOL,4,N1,A,Y$NSOL,5,NMAX,A,Y$PLVAR,4,5
DERIV,6,2,1$DERIV,6,6,1$DERIV,7,3,1$DERIV,7,7,1$PLVAR,6,7
ESOL,8,102,,CONT,STAT$ESOL,9,102,,CONT,PRES$PLVAR,8$PLVAR,9 |