finish /clear /filname,qiaopan,1 /prep7 et,1,beam189 et,2,mass21 mp,ex,1,3.45e10 mp,nuxy,1,0.3 mp,dens,1,2500 r,1,1,1,1,1,1,1SECTYPE,1,BEAM,MESH,1sec SECOFFSET,ORIG, SECREAD,1sec,SECT, ,MESH SECTYPE,2,BEAM,MESH,2sec SECOFFSET,ORIG, SECREAD,2sec,SECT, ,MESH SECTYPE,3,BEAM,MESH,3sec SECOFFSET,ORIG, SECREAD,3sec,SECT, ,MESH SECTYPE,4,BEAM,MESH,4sec SECOFFSET,ORIG, SECREAD,4sec,SECT, ,MESH SECTYPE,5,BEAM,MESH,5sec SECOFFSET,ORIG, SECREAD,5sec,SECT, ,MESHSECTYPE,6,BEAM,MESH,6sec SECOFFSET,ORIG, SECREAD,6sec,SECT, ,MESH SECTYPE,7,BEAM,MESH,7sec SECOFFSET,ORIG, SECREAD,7sec,SECT, ,MESHSECTYPE,8,BEAM,MESH,8sec SECOFFSET,ORIG, SECREAD,8sec,SECT, ,MESHSECTYPE,9,BEAM,MESH,9secSECOFFSET,ORIG, SECREAD,9sec,SECT, ,MESHSECTYPE,10,BEAM,MESH,10sec SECOFFSET,ORIG, SECREAD,10sec,SECT, ,MESH SECTYPE,11,BEAM,MESH,11sec SECOFFSET,ORIG, SECREAD,11sec,SECT, ,MESH SECTYPE,12,BEAM,MESH,100sec SECOFFSET,ORIG, SECREAD,4sec,SECT, ,MESH SECTYPE,13,BEAM,MESH,13sec SECOFFSET,ORIG, SECREAD,13sec,SECT, ,MESH SECTYPE,14,BEAM,MESH,14sec SECOFFSET,ORIG, SECREAD,14sec,SECT, ,MESH! /pnum,kp,1 /pnum,line,1!创建主梁 k,1 k,2,1.25 k,3,3.75 k,4,5 k,5,9 k,6,13 k,7,17 k,8,21 k,9,25 k,10,29.5 k,11,34 k,12,35*do,i,1,11,1 l,i,i+1 *enddok,100,1 !创建一个定位点allsel gplot!赋予属性*do,i,1,11,1 lsel,s,i latt,1,1,100,i *enddo lsel,s,all lesize,all,4 lmesh,all！镜像处理 lsel,s,1,10,1 lsymm,x,all csys,4 kwpave,12 lsel,s,1,21,1 lplot lsymm,x,all!修改单元坐标系 lsel,s,12,42,1 lreverse,all,0 lsel,s,33,42,1 lreverse,all,0 allsel gplot!增加边跨合拢段和现浇段 kgen,2,23,-2 kgen,2,23,-2.9 kgen,2,23,-4.4 kgen,2,23,-5.9 kgen,2,46,2 kgen,2,46,2.9 kgen,2,46,4.4 kgen,2,46,5.9 l,50,49 l,49,48 l,48,47 l,47,23 l,46,51l,51,52 l,52,53 l,53,54!赋予边跨端部截面属性 lsel,s,45,48,1latt,1,1,100,11 !100 为定位点 lsel,s,44,49,5 latt,1,1,100,13 lsel,s,43,50,7 latt,1,1,100,14lsel,s,45,48,1 lsel,a,44,49,5 lsel,a,43,50,7 lesize,all,4 lmesh,all numcmp,all nummrg,all！ d,155,ux d,155,uy d,155,uz d,1,uy d,1,uz d,78,uy d,78,uz d,176,uy d,176,uz!进行求解 allsel gplot! NSEL,S,LOC,Y,0 *SET,J,155 NMAX=NDINQR(0,13)*DIM,PINODE,NMAX ！定义节点数组，用于存放节点号 *DIM,NODEX,NMAX ！定义节点 X 坐标数组，用于存放各个节点的 X 坐标*SET,PINODE(1),J *SET,NODEX(1),NX(J)*DO,I,2,NMAX*SET,NI,NNEAR(J)*SET,PINODE(I),NINSEL,U,J*SET,J,NI *ENDDO ALLSEL,ALL！*SET,V,8/3.6 *SET,FREQ,49 *SET,PI,3.1415926 *SET,W,2*PI*FREQ *SET,DETT,0.25/V *SET,FMAX,257000 *SET,G,7000*9.8FINISH /SOLU ANTYPE,TRANS TRNOPT,FULL TIMINT,OFF ACEL,10 TIME,1E-5 NSUBST,2 KBC,1 SOLVE TIMINT,ONOUTRES,ALL,ALL*DO,I,2,118,1*IF,I,LT,117,THENOUTRES,ALL,ALLTIME,(I-1)*DETTFDELE,ALL,ALLF,PINODE(I),FY, -GALLSELSOLVE*ELSETIME,400 NSUBST,40 FDELE,ALL,ALLALLSELSOLVE*ENDIF*ENDDO/POST26 /rgb,index,100,100,100,0 /rgb,index,80,80,80,13 /rgb,index,60,60,60,14 /rgb,index,0,0,0,15 /replot!跨中位移 NSOL,2,38,U,Y PLVAR,2!中跨跨中应力 ESOL,3,13,38,SMISC,35 PLVAR,3 /IMAGE,SAVE,ZKKZYL,JPG！中跨 1/4 处应力 ESOL,5,7,17, SMISC,35 PLVAR,5 /IMAGE,SAVE,ZK4YL,JPG！边跨跨中应力 ESOL,6,20,60, SMISC,35 PLVAR,6 /IMAGE,SAVE,BKKZYL,JPG!边跨距边支座 1/4 处应力 ESOL,7,23,66, SMISC,35 PLVAR,7 /IMAGE,SAVE,BK4YL,JPG！中间支点应力随时间 ESOL,8,1,1, SMISC,35 PLVAR,8 /IMAGE,SAVE,ZJZDYL,JPG!输出中跨跨中应力 *CREAT,ANSUITMP*CFOPEN,ZKKZYL,TXT *GET,num_var,VARI,3,NSETS *DIM,JZ1,ARRAY, num_var,2 VGET,JZ1(1,1),1 VGET,JZ1(1,2),3 *VWRITE,JZ1(1,1),JZ1(1,2) (F20.5,5X,F20.5) *CFCLOS *END /INPUT,ANSUITMP！输出中跨 1/4 处应力 *CREAT,ANSUITMP *CFOPEN,ZKKZ4YL,TXT *GET,num_var,VARI,5,NSETS *DIM,JZ2,ARRAY, num_var,2 VGET,JZ2(1,1),1 VGET,JZ2(1,2),5 *VWRITE,JZ2(1,1),JZ2(1,2) (F20.5,5X,F20.5) *CFCLOS *END /INPUT,ANSUITMP!输出边跨跨中 *CREAT,ANSUITMP *CFOPEN,BKKZYL,TXT *GET,num_var,VARI,6,NSETS *DIM,JZ3,ARRAY, num_var,2 VGET,JZ3(1,1),1 VGET,JZ3(1,2),6 *VWRITE,JZ3(1,1),JZ3(1,2) (F20.5,5X,F20.5) *CFCLOS *END /INPUT,ANSUITMP!输出边跨距边支座 1/4 应力 *CREAT,ANSUITMP *CFOPEN,BKKZ4YL,TXT *GET,num_var,VARI,7,NSETS *DIM,JZ4,ARRAY, num_var,2 VGET,JZ4(1,1),1 VGET,JZ4(1,2),7*VWRITE,JZ4(1,1),JZ4(1,2) (F20.5,5X,F20.5) *CFCLOS *END /INPUT,ANSUITMP！输出中间支点应力*CREAT,ANSUITMP *CFOPEN,ZJZDYL,TXT *GET,num_var,VARI,8,NSETS *DIM,JZ5,ARRAY, num_var,2 VGET,JZ5(1,1),1 VGET,JZ5(1,2),8 *VWRITE,JZ5(1,1),JZ5(1,2) (F20.5,5X,F20.5) *CFCLOS *END /INPUT,ANSUITMP