1、/prep7ET,1,SOLID45ET,2,COMBIN39,1 ET,3,COMBIN39,2MP,EX,1,2.10E11MP,NUXY,1,0.3MP,DENS,1,7800MP,DAMP,2,0.3R,1,-2,-15486.62,-0.06904,-15486.62,0,0 !定义土弹簧刚度曲线（水平方向）RMORE,0.06904,15486.62,2,15486.62,R,2,-2,-10000,-0.06904,-10000,0,0 !定义土弹簧刚度曲线（竖直向下） RMORE,0.06904,10000,2,10000,R,3 *SET,Ri,0.381 *SET,T,0.
2、0238*SET,P,10e6*SET,chang,1*T*SET,h,50cyl4,0,0,ri,0,ri+t,90,-hcyl4,0,0,ri,90,ri+t,180,-hcyl4,0,0,ri,180,ri+t,270,-hcyl4,0,0,ri,270,ri+t,360,-hwpoffs,z-changVSBW,ALLwpcsys,-1,0wpoffs,z-2VSEL,S,5,11,2VSBW,ALLwpcsys,-1,0vsel,allCSYS,0VSYMM,z,all, , , ,0,0vsel,allvglue,allnumcmp,allALLSEL,ALL LSEL,S,32,
3、35,1LSEL,a,41,42,1LSEL,a,48,49,1LSEL,a,54,61,1LSEL,a,75,86,1LSEL,a,129,132,1Lesize,all,1LSEL,S,12LSEL,a,16LSEL,a,9LSEL,a,10LSEL,a,21LSEL,a,22LSEL,a,51LSEL,a,52LSEL,a,67LSEL,a,68LSEL,a,43LSEL,a,44LSEL,a,45LSEL,a,46LSEL,a,63LSEL,a,64Lesize,all,24LSEL,S,99LSEL,a,100LSEL,a,5LSEL,a,7LSEL,a,103LSEL,a,104L
4、SEL,a,115LSEL,a,116Lesize,all,18ALLSEL,ALL vsel,all TYPE,1 MAT,1Real,3VSWEEP,all nsel,s,loc,x,0,0.4048csys,1nsel,r,loc,x,0.4048csys,0ngen,2,100000,all,3nsel,s,loc,x,0,0.4048csys,1nsel,r,loc,x,0.4048csys,0*get,nnub,node,count*get,nmin,node,num,min*dim,A,nnubn0=nminA(1)=n0*do,j,2,nnubn0=ndnext(n0)A(j)
5、=n0*enddotype,2mat,2real,1*do,i,1,nnub,1e,A(i),A(i)+100000*enddonsel,s,loc,x,0,-0.4048csys,1nsel,r,loc,x,0.4048csys,0ngen,2,200000,all,-3nsel,s,loc,x,0,-0.4048csys,1nsel,r,loc,x,0.4048csys,0*get,nnub,node,count*get,nmin,node,num,min*dim,B,nnubn0=nminB(1)=n0*do,j,2,nnubn0=ndnext(n0)B(j)=n0*enddotype,
6、2mat,2real,1*do,i,1,nnub,1e,B(i),B(i)+200000*enddonsel,s,loc,y,0,-0.4048csys,1nsel,r,loc,x,0.4048csys,0ngen,2,300000,all,-3nsel,s,loc,y,0,-0.4048csys,1nsel,r,loc,x,0.4048csys,0*get,nnub,node,count*get,nmin,node,num,min*dim,C,nnubn0=nminC(1)=n0*do,j,2,nnubn0=ndnext(n0)C(j)=n0*enddotype,3mat,2real,2*d
7、o,i,1,nnub,1e,C(i),C(i)+300000*enddoALLSEL,ALLCSYS,1Asel,s,loc,x,riNsla,s,1SF,ALL,PRES,pALLSEL,ALLcsys,1nsel,s,loc,x,2,4d,all,allcsys,0nsel,s,loc,z,-hnsel,a,loc,z,hdsym,symm,zallsfini/nerr,-1,10000000/solucsys,0ANTYPE,TRANS TRNOPT,FULL NROPT,full,on timint, off time, 1e-8 acel,9.8 cnvtol,f,1e5,0.1,2
8、 nsubst, 2 KBC, 1 Lswrite, 1 Solve timint, on nsubst,1neqit,30 nlgeom,on *dim,acez,TABLE,500,1,1,time,accel*create,ansuitmp*tread,acez,el_centro,dat,1*end/input,ansuitmpOUTRES,ALL,ALL,t1=0.02t2=10dt=0.02*do,tm,t1,t2,dttime,tmalpha,0 beta,0 acel,acez(tm),9.8,solve*enddoFini/POST26 NUMVAR,2e7 NSOL,2,4
9、233,U,x, XVAR,1 PLVAR,2, , , , , , , , , , DERIV,3,2,1,1PLVAR,3, , , , , , , , , ,DERIV,4,3,1,1PLVAR,5, , , , , , , , , ,Ansol,6,4233,s,y,sy_6XVAR,1 PLVAR,6,Ansol,7,4233,epel,y,epely_7XVAR,1 PLVAR,7,A B CChang 4233 6920 69382 3721 4701 51514 3720 4702 51506 3719 4700 51498 3718 4699 514810 3717 4698 5147
