ansys土木工程应用实例命令流.doc

第三章回弹/filname,huitan/units,si/prep7et,1,plane42keyopt,1,3,2mp,ex,1,2.160687961e11mp,nuxy,1,0.3TB,MISO,1,1,25, TBTEMP,0TBPT,0.000814,175.88e6,TBPT,0.001531,324.98e6,TBPT,0.002819,591.26e6,TBPT,0.003846,801.457e6,TBPT,0.004239,873.22e6,TBPT,0.004594,924.399e6,TBPT,0.004955,956.88e6,TBPT,0.00531,976.99e6,TBPT,0.005699,990.2138e6,TBPT,0.006261,1000e6,TBPT,0.006957,1008.34e6,TBPT,0.007804,1015.5e6,TBPT,0.008779,1021.185e6,TBPT,0.009827,1026.47e6,TBPT,0.010929,1031.289e6,TBPT,0.012066,1035.444e6,TBPT,0.0132,1039.693e6,TBPT,0.014354,1044.886e6,TBPT,0.015494,1049.796e6,TBPT,0.016631,1053.951e6,TBPT,0.017785,1058.106e6,TBPT,0.019039,1062.827e6,TBPT,0.020276,1067.643e6,csys,1k,1,0.00512k,2,0.0064kgen,36,1,2,1,10*do,j,1,69,2 a,j,j+1,j+3,j+2*enddoa,71,72,2,1savelesize,all,4amesh,allnummrg,allnumcmp,allsavecp,1,uy,166,170,171,172,173finish/solunsel,s,loc,y,270d,all,allnplotallselnplotnsel,s,loc,y,90d,all,uxnplotallseloutres,all,alltime,1NSUBST,100nlgeom,on d,166,uy,-0.00128solve*get,def1,node,166,U,Ynsel,s,loc,y,90ddele,all,uyallself,166,fy,0solve*get,def2,node,166,U,Yfinish/post1 !进入后处理器set,list !显示所有计算的结果set,1,last,1 !读取第一个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状plnsol,s,eqv,1,1 !显示节点的等效应力set,2,last,1 !读取第二个加载步最后子步的结果pldisp,1 !同时显示变形前与变形后结构的形状。plnsol,s,eqv,1,1 !显示节点的等效应力自适应划分/PREP7SMRT,OFFANTYPE,STATIC ET,1,PLANE42,3 ! 定义平面应力单元有单元厚度 MP,EX,1,210E9 ! 定义材料属性MP,NUXY,1,.3R,1,0.1 ! 给定厚度LOCAL,11,1,0.5 ! 定义区域坐标系统为圆柱坐标系，Y轴是X的0.5倍。K,1,2,90 K,2,2,0 L,1,2 LOCAL,12,1,0.8461585 ! 定义区域坐标系统为圆柱坐标系，Y轴是X的0.846倍。K,3,3.25,90K,4,3.25,0.0L,3,4CSYS,0 !回到卡式坐标系下。L,2,4L,1,3AL,2,4,1,3 ! 创建面积DL,4,1,SYMM ! 创建对称面约束 DL,3,1,SYMMSFL,2,PRES,-10E6 ! 定义线上的分布载荷MSHK,2 ! 混合网格划分MSHA,0,2D ! 四边形网格SAVE ! 备份数据库FINISHADAPT,7,7,1 ! 求解次数为nsoln=4次，能量误差百分比stargt=7,单元尺寸改变最大比例因子facmx=1/POST1EPLOT ! 绘制单元PRERR ! 打印能量误差百分比NSEL,S,LOC,Y,0.0NSEL,R,LOC,X,2.0*GET,MNODE,NODE,NUM,MAX*GET,SY_D,NODE,MNODE,S,Y ! GET DESIRED STRESS SY VALUENSEL,ALLESEL,ALLFINISH第四章钢筋混凝土板!在四个角点增加钢性垫片，在中心点处直接施加-2mm，压碎关闭计算完毕/clear,nostart/config,nres,5000/prep7!定义单元类型et,1,solid65 !混凝土单元et,2,solid45keyopt,1,6,3keyopt,1,5,1*afun,degR,1,2,0.01,90,0,2,0.01,0,0 !整体模型，实常数!定义材料属性mp,ex,1,24000!混凝土材料属性mp,prxy,1,0.2tb,conc,1,1,9tbdata,0.35,1,3.1125,-1mp,ex,2,2e5!纵向受拉钢筋，钢支座材料属性mp,prxy,2,0.25tb,bkin,2,1,2,1tbdata,360,20000blc4,1000,1000,100blc4,100,100,-50wpoff,1000,0,0blc4,-100,100,-50wpoff,0,1000,0blc4,-100,-100,-50wpoff,-1000,0,0blc4,100,-100,-50/view,1,1,1,1/ANG,1,240,ZS,1 /replot /REP,FASTlsel,s,loc,z,50lesize,all,2lsel,alllsel,s,loc,x,500lesize,all,20lsel,alllsel,s,loc,y,500lesize,all,20vsel,s,loc,z,50vmesh,alltype,2mat,2lsel,s,loc,z,-25lesize,all,1lsel,alllsel,s,loc,z,-50lesize,all,2vsel,s,loc,z,-25vmesh,allnummrg,allnumcmp,allfini/soluasel,s,loc,z,-50da,all,allnsel,allnsel,s,loc,z,100nsel,r,loc,y,500nsel,r,loc,x,500d,all,uz,-2nlgeom,onNSUBST,100 OUTRES,ALL,ALL NEQIT,40pred,onCNVTOL,F, ,0.05,2allselSOLVE fini /post1PLNSOL,U,z,0,1Esel,s,type,1PLNSOL,S,1,0,1!画第一主应力ETABLE,rebar_1,NMISC,43 !获取rebar1平均等效应力PLETAB,REBAR_1,AVGETABLE,rebar_2,NMISC,47 !获取rebar2平均等效应力PLETAB,REBAR_2,AVGetable,seqv_i,nmisc,5!获取solid单元i节点的平均等效应力PLETAB,seqv_i,AVGETABLE,11,NMISC,53!获取第一个积分点的状态ETABLE,22,NMISC,60!获取第二个积分点的状态ETABLE,33,NMISC,67!获取第三个积分点的状态ETABLE,44,NMISC,74!获取第四个积分点的状态ETABLE,55,NMISC,81!获取第五个积分点的状态ETABLE,66,NMISC,88!获取第六个积分点的状态ETABLE,77,NMISC,95!获取第七个积分点的状态ETABLE,88,NMISC,102!获取第八个积分点的状态PLETAB,11,AVGPLETAB,22,AVGPLETAB,33,NOAVPLETAB,44,AVGPLETAB,55,NOAV/DEVICE,VECTOR,1PLCRACK,0,0!显示所有开裂位置PLCRACK,0,1!显示第一开裂位置钢筋混凝土梁!采用均匀网格划分，同时取消两种垫片,DMX=6.415,最大第一主应力263.629，28个子部计算完毕/clear,nostart/config,nres,5000/prep7!定义单元类型et,1,solid65 !混凝土单元et,2,pipe20!钢筋单元et,3,plane42!用于拉伸成体单元之用!定义实常数r,1,18,8.99!纵向受拉钢筋实常数r,2,8,3.99!纵向受压钢筋，箍筋实常数r,3!混凝土!定义材料属性mp,ex,1,24000!混凝土材料属性mp,prxy,1,0.2tb,conc,1,1,9tbdata,0.4,1,3,-1mp,ex,2,2e5!纵向受拉钢筋mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,350mp,ex,3,2e5!横向箍筋，受压钢筋材料属性mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,200/pnum,node,1/pnum,elem,1!产生所有的节点*!n,1n,9,150fill,1,9ngen,11,9,1,9,1,30ngen,19,1000,1,99,1,-75/view,1,1,1,1!箍筋，受压钢筋*!type,2real,2mat,3!水平箍筋*do,ii,11,16,1 e,ii,ii+1*enddo*do,ii,83,88,1 e,ii,ii+1*enddo!竖直箍筋*do,ii,11,74,9 e,ii,ii+9*enddo*do,ii,17,80,9 e,ii,ii+9*enddo!产生整个模型的箍筋egen,19,1000,1,28,1!此时最大的单元编号为532!纵向受压钢筋*do,ii,83,17083,1000 e,ii,ii+1000*enddo*do,ii,89,17089,1000 e,ii,ii+1000*enddo!纵向受拉钢筋*!type,2real,1mat,2*do,ii,11,17011,1000 e,ii,ii+1000*enddo*do,ii,17,17017,1000 e,ii,ii+1000*enddo!此时最大单元编号为604/view,1,1,1,1/pnum,elem,0/pnum,node,0/eshape,1eplot!*生成面单元，以便拉伸成体单元!k,1k,2,150,k,3,150,300k,4,300a,1,2,3,4lsel,s,loc,y,0lsel,a,loc,y,300lesize,all,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,150lesize,all,10type,3amesh,alltype,1!拉伸成混凝土单元real,3mat,1extopt,esize,18extopt,aclear,1vext,all,-1350/pnum,mat,1/pnum,node,0/pnum,elem,0/view,1,1,1,1eplotallselnummrg,allnumcmp,alleplot/solunsel,s,loc,y,-0!右边界约束nsel,r,loc,z,-1275d,all,uyd,all,uxasel,s,loc,z,0!对称面约束da,all,symmnsel,all!施加外部载荷nsel,s,loc,y,300nsel,r,loc,z,-600d,all,uy,-5allselgplotnlgeom,onNSUBST,200 OUTRES,ALL,ALL NEQIT,50pred,onCNVTOL,F, ,0.05,2,0.5allselSOLVE FINISH /post1allsel/DEVICE,VECTOR,1PLCRACK,0,1!显示第一开裂位置PLCRACK,0,2!显示第二开裂位置!/DEVICE,VECTOR,0!PLNSOL,U,Y,0,1!ESEL,S,TYPE,2!PLNSOL,S,1,0,1预应力钢筋混凝土/PREP7 ET,1,LINK8 ET,2,SOLID65ET,3,PLANE42R,1,6126 MP,EX,1,2e5MP,PRXY,1,0.3 MP,ALPX,1,2e-5 MP,EX,2,3.354e MP,PRXY,2,0.2 BLC4, , ,780,210BLC4,275,210,230,1690 BLC4,-570,1900,1920,200 wpoff,390,191 wprot,0,90ASBW,1 wprot,0,0,90ASBW,4ASBW,5/PNUM,AREA,1/PNUM,LINE,1gplotaglue,allnummrg,allnumcmp,allSAVElesize,8,20lesize,21,8lesize,22,8lesize,5,4lesize,9,2lesize,7,2lsel,s,21,22,1lsel,a,5LCCAT,allallselamesh,6LSEL,s,LCCAlDELE,allallsellesize,4,20 lesize,6,20lesize,18,2lesize,19,2lsel,s,18,19,1LCCAT,allallselamesh,4LSEL,s,LCCAlDELE,allallsellesize,2,6lesize,3,8lesize,14,2lesize,12,2lsel,s,2,18,16LCCAT,allallselamesh,2LSEL,s,LCCAlDELE,allallsellesize,20,6lesize,17,8lesize,13,2lsel,s,19,20,1LCCAT,allallselamesh,5LSEL,s,LCCAlDELE,allallsellesize,10,6lesize,15,6lesize,11,6lesize,1,8lesize,16,8allselamesh,1amesh,3SAVETYPE, 2 EXTOPT,ESIZE,40,0, EXTOPT,ACLEAR,1 EXTOPT,ATTR,0,0,0 MAT,2 REAL,1 VEXT,all, , ,0,0,4000, /VIEW, 1 ,1,1,1 /ANG, 1 EPLOT TYPE,1 MAT,1REAL,1 lmesh,29FINISH /SOLDL,23, ,ALL,DL,38, ,ALL, DL,1, ,UY, DL,16, ,UY, DK,1, , , ,0,UX,UY SAVEBFL,29,TEMP,-375 CNVTOL,F, ,0.05,2, ,NSUBST,10,0,0 OUTRES,ERASEOUTRES,ALL,ALL TIME,5 SOLVE4.6混凝土开裂模拟/TITLE,Crack/PREP7ET,1,SOLID65KEYOPT,1,5,2ET,2,LINK8ET,3,PIPE16R,1R,2,0.15R,3,1,0.5MP,EX,1,2E6MP,PRXY,1,0TB,CONC,1TBDATA,-1MP,EX,2,30E6 MP,PRXY,2,0.3!BLC4, , ,5,6,1.5/VIEW, 1 ,1,1,1 /ANG, 1 /REP,FAST!/PNUM,LINE,1LSEL,LOC,Z,0.1,1.4LSEL,R,LOC,Y,0TYPE, 2 MAT, 2REAL, 2LESIZE,ALL,1LMESH,ALL ALLSELLSEL,LOC,Z,1.5LSEL,R,LOC,Y,1,4TYPE, 3 MAT, 1REAL, 3LESIZE,ALL,4LMESH,ALLALLSELTYPE, 1 MAT, 1REAL, 1LSEL,LOC,X,1,4LESIZE,ALL,1ALLSELMSHAPE,0,3DMSHKEY,1VMESH,ALLFINISH!/PREP7 CE,1,0,4,UZ,1,7,UZ,-1,7,ROTX,-3,CE,2,0,6,UZ,1,7,UZ,-1,7,ROTX,-1.5, CE,3,0,8,UZ,1,7,UZ,-1,7,ROTX,1.5, CE,4,0,5,UZ,1,7,UZ,-1,7,ROTX,3, CE,5,0,1,UZ,1,11,UZ,-1,11,ROTX,-3, CE,6,0,12,UZ,1,11,UZ,-1,11,ROTX,-1.5, CE,7,0,10,UZ,1,11,UZ,-1,11,ROTX,1.5,CE,8,0,9,UZ,1,11,UZ,-1,11,ROTX,3, NSEL,LOC,Z,0D,ALL,ALLALLSELNSEL,LOC,Z,1.5D,ALL,ROTYD,ALL,ROTZALLSELF,7,MX,300F,11,MX,300/SOLNSUBST,5,0,0AUTOTS,1SOLVE!/POST1SET,LAST *GET,SCON,NODE,17,S,ZETAB,ST,LS,1ESORT,ST*GET,STL,SORT,MAX*status,parm /DEVICE,VECTOR,1PLCRACKSAVEFINISH4.7混凝土浇筑模拟/TITLE,CASTCONCRETE/PREP7ET,1,SOLID70MP,KXX,1,300.89 MP,C,1,0.94500 MP,DENS,1,2400 MP,KXX,2,300.89 MP,C,2,0.94500 MP,DENS,2,2600 BLOCK,50,2,40, BLOCK,-25,75,-10,0,-20,60, /VIEW, 1 ,1,1,1 /ANG, 1 LSEL,LOC,Z,0LSEL,R,LOC,X,1,49LESIZE,ALL,10ALLSELLSEL,LOC,Z,40LSEL,R,LOC,X,1,49LESIZE,ALL,10ALLSELLSEL,LOC,X,0LSEL,R,LOC,Z,1,39LESIZE,ALL,10ALLSELLSEL,LOC,X,50LSEL,R,LOC,Z,1,39LESIZE,ALL,10ALLSELLSEL,LOC,Z,0LSEL,R,LOC,Y,1,1.5LESIZE,ALL,2ALLSELLSEL,LOC,Z,40LSEL,R,LOC,Y,1,1.5LESIZE,ALL,2ALLSELLSEL,LOC,Z,-20LSEL,R,LOC,X,1,49LESIZE,ALL,20ALLSELLSEL,LOC,Z,60LSEL,R,LOC,X,1,49LESIZE,ALL,20ALLSELLSEL,LOC,X,-25LSEL,R,LOC,Z,1,39LESIZE,ALL,20ALLSELLSEL,LOC,X,75LSEL,R,LOC,Z,1,39LESIZE,ALL,20ALLSELLSEL,LOC,Z,-20LSEL,R,LOC,Y,-9,-1.3LESIZE,ALL,10ALLSELLSEL,LOC,Z,60LSEL,R,LOC,Y,-9,-1.3LESIZE,ALL,10ALLSELVSEL,LOC,Y,1,1.5MAT,1VMESH,1ALLSELVSEL,LOC,Y,-9,-1.5MAT,2VMESH,2ALLSELNUMMRG,ALLNUMCMP,NODENROPT,FULLESEL,S,MAT,2NSLE,S IC,ALL,TEMP,16.5ASEL,LOC,X,-25ASEL,A,LOC,X,75ASEL,A,LOC,Y,-10ASEL,A,LOC,Z,-20ASEL,A,LOC,Z,60NSLA,SSF,ALL,HFLUX,0ALLSELNSEL,LOC,X,49,76NSEL,R,LOC,Y,0NSEL,A,LOC,X,-26,1NSEL,R,LOC,Y,0NSEL,A,LOC,Z,39,61NSEL,R,LOC,Y,0NSEL,A,LOC,Z,-21,1NSEL,R,LOC,Y,0SF,ALL,CONV,2016,20ESEL,MAT,1EKILL,ALLNSLE,SIC,ALL,TEMP,40NSEL,LOC,Y,0,1ESLN,S,1EALIVE,ALLALLSELNSEL,LOC,X,1,49NSEL,R,LOC,Z,1,39NSEL,R,LOC,Y,0CM,DD,NODE ALLSELNSEL,LOC,Y,0,1ESLN,S,1NSLE,SNSEL,U,DD SF,ALL,CONV,2016,20HE00=47880/1.862BFE,ALL,HGEN,HE00ALLSEL/SOLUANTYPE,TRANSTIME,1AUTOTS,0DELTIM,1KBC,0SOLVE/POST1PLNSOL,TG,SUM/PREP7NSEL,LOC,Y,0,1ESLN,S,1HE00=47880*(2/2.862- 1/1.862)BFE,ALL,HGEN,HE00ALLSEL/SOLUANTYPE,TRANSTIME,2AUTOTS,0DELTIM,1KBC,0SOLVE/POST1PLNSOL,TG,SUM/PREP7ESEL,MAT,1NSLE,SSFEDELE,ALL,ALL,CONVALLSELNSEL,LOC,X,1,49NSEL,R,LOC,Z,1,39NSEL,R,LOC,Y,0,1CM,KK,NODE NSEL,LOC,Y,0,2ESLN,S,1NSLE,SNSEL,U,KK SF,ALL,CONV,2016,20ALLSELNSEL,LOC,Y,0,1ESLN,S,1HE00=47880*(3/3.862-2/2.862)BFE,ALL,HGEN,HE00ALLSELNSEL,LOC,Y,1,2ESLN,S,1HE00=47880/1.862BFE,ALL,HGEN,HE00ALLSEL/SOLUANTYPE,TRANSTIME,3AUTOTS,0DELTIM,1KBC,0SOLVE/POST1PLNSOL,TG,SUM/PREP7NSEL,LOC,Y,0,1ESLN,S,1HE00=47880*(4/4.862-3/3.862)BFE,ALL,HGEN,HE00ALLSELNSEL,LOC,Y,1,2ESLN,S,1HE00=47880*(2/2.862-1/1.862)BFE,ALL,HGEN,HE00ALLSEL/SOLUANTYPE,TRANSTIME,4AUTOTS,0DELTIM,1KBC,0SOLVE/POST1PLNSOL,TG,SUMSAVEFINISH!/PREP7ET,1,SOLID45MP,EX,1,0.325E10MP,DENS,1,2400MP,PRXY,1,0.2MP,ALPX,1,0.9E-5MP,REFT,1,15MP,EX,2,2.08E10MP,DENS,2,2600MP,PRXY,2,0.2MP,ALPX,2,0.9E-5ESEL,MAT,2MPCHG,2,ALLALLSEL/SOLUASEL,LOC,X,-25ASEL,A,LOC,X,75ASEL,A,LOC,Y,-10ASEL,A,LOC,Z,-20ASEL,A,LOC,Z,60NSLA,SD,ALL,ALLALLSELESEL,MAT,1EKILL,ALLNSEL,LOC,Y,0,1ESLN,S,1EALIVE,ALLMPCHG,1,ALLALLSELACEL,0,9.8,0LDREAD,TEMP,1,CASTCONCRETE,rth,ANTYPE,TRANSNLGEOM,ONCNVTOL,F,0.2,2,200TIME,1AUTOTS,0DELTIM,1KBC,0SOLVE/SOLULDREAD,TEMP,2,CASTCONCRETE,rth,TIME,2SOLVE/SOLUNSEL,LOC,Y,1,2ESLN,S,1EALIVE,ALLMPCHG,1,ALLALLSELLDREAD,TEMP,3,CASTCONCRETE,rth,TIME,3SOLVE/SOLULDREAD,TEMP,4,CASTCONCRETE,rth,TIME,4SOLVEFINISHSAVE!/POST1PLNSOL,U,YPLNSOL,S,YPLESOL,S,Y5.2行架桥/TITLE,TRUSS/PREP7ET,1,BEAM4ET,2,LINK8ET,3,SHELL63R,1,3,0.25,2.25,3,1,-0.0001,1R,2,3,0.25,2.25,3,1,1R,3,2,1/6,2/3,2,1,2/3R,4,1R,5,1,-1.5E-4R,6,0.3MP,EX,1,35E10MP,DENS,1,3000MP,PRXY,1,1/6MP,EX,2,21E10MP,DENS,2,3500MP,PRXY,2,0.3MP,EX,3,3.5E10MP,DENS,3,2500MP,PRXY,3,1/6MP,EX,4,3.5E10MP,DENS,4,2500MP,PRXY,4,1/6!N,1,0,-5,0NGEN,4,1,1,12N,5,0,-5,12NGEN,3,1,5,12NSYM,Y,7,ALLTYPE,1REAL,1MAT,1E,1,2E,2,3E,3,4E,8,9E,9,10E,10,11TYPE,1REAL,2MAT,1E,5,6E,6,7E,12,13E,13,14E,5,12E,6,13E,7,14TYPE,1REAL,3MAT,2E,7,4E,14,11TYPE,2REAL,4MAT,3E,5,2E,2,7E,12,9E,9,14TYPE,2REAL,5MAT,3E,5,1E,6,2E,7,3E,12,8E,13,9E,14,10TYPE,3REAL,6MAT,4E,1,2,9,8E,2,3,10,9E,3,4,11,10NSYM,X,14,ALLESYM,14,ALLNUMMRG,ALLNUMCMP,ALLSAVE/SOLUD,17,UXD,17,UYD,17,UZD,22,UXD,22,UYD,22,UZD,4,UYD,4,UZD,11,UYD,11,UZACEL,10F,1,FZ,-10000F,8,FZ,-10000SAVESOLVE!/POST1PLVECT,UPRNSOL,U,ZPRESOL,ELEMFINISHSAVE5.3悬索桥/TITLE,CABLE/PREP7ET,1,LINK10ET,2,BEAM3R,1,1,1E-5R,2,0.02,1E-5R,3,0.5,1/24,1MP,EX,1,25E10MP,PRXY,1,0.1MP,DENS,1,1000MP,EX,2,30E10MP,PRXY,2,0.1MP,DENS,2,1000A=250*DO,I,1,10 X=10*(I-1) Y=A*COSH(X/A)-A N,I,X,Y*ENDDON,11,100,20*DO,I,1,10 X=10*I Y=0 N,11+I,X,Y*ENDDOtype,1real,1mat,1*DO,I,1,10 E,I,I+1*ENDDOtype,1real,2mat,1*DO,I,2,10 E,I,I+10*ENDDOtype,2real,3mat,2E,1,12*DO,I,1,9 E,I+11,I+12*ENDDONSYM,X,21,ALLESYM,21,ALLNUMMRG,ALLNUMCMP,ALLSAVE!/SOLUD,41,UXD,41,UYD,21,UYD,31,UXD,31,UYD,11,UXD,11,UYNSEL,U,NODE,41NSEL,U,NODE,21D,ALL,UZALLSELACEL,10,NROPT,FULLESEL,REAL,2ESEL,A,REAL,3EKILL,ALLALLSELNLGEOM,ONTIME,1NSUBST,10,20,5OUTRES,ERASEOUTRES,ALL,ALLSSTIF,ONSOLVE/POST1ESEL,LIVEPLVECT,UALLSEL/PREP7UPGEOM,1,LAST,LAST,Cable,RST/SOLUSOLVE/POST1ESEL,LIVEPLVECT,UALLSELSAVE/SOLUF,1,FY,-50000F,2,FY,-25000F,22,FY,-25000TIME,2SOLVE/POST1ESEL,LIVEPLVECT,UALLSEL/SOLUDOFSEL,S,FYFCUM,ADD,1DOFSEL,ALL*DO,I,1,9XX0=I*10XX=(I+1)*10NSEL,LOC,X,-XX0NSEL,R,LOC,Y,0.1,30NSEL,A,LOC,X,XX0NSEL,R,LOC,Y,0.1,30F,ALL,FY,-25000ALLSELNSEL,LOC,X,XXNSEL,R,LOC,Y,0.1,30NSEL,A,LOC,X,-XXNSEL,R,LOC,Y,0.1,30F,ALL,FY,-25000ALLSELTIME,I+2SOLVE*ENDDO/POST1ESEL,LIVEPLVECT,UALLSELSAVE!/SOLUNSEL,LOC,Y,0.1,30NSEL,A,NODE,1FDELE,ALL,FYALLSELESEL,U,LIVEEALIVE,ALLALLSELTIME,12SOLVE/POST1ESEL,LIVEPLVECT,UALLSELFINISHSAVE5.4.1!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。/config,nres,2000/PREP7 !建模ET,1,BEAM3 !梁单元MP,EX,1,207E9 !弹性模量MP,DENS,1,2000 !密度MP,NUXY,1,0.3R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高NELEM=100LBEAM=32 !梁长度deltL=LBEAM/NELEMV=120f=1000w=10 !移动荷载速度NNODE=NELEM+1deltt=deltL/V*3.6*do,i,1,NNODE,1N,i,(i-1)*deltL !节点点矩为0.3*enddo*do,i,1,NELEME,i,i+1*enddofinish/solud,1,alld,NNODE,allantype,trans *do,ii,1,101+8,1outres,all,all !输出所有加载子步的所有结果time,ii*deltt !每个荷载步的最终求解时间nsubst,5 !每个荷载步包含五个荷载子步*if,ii,lt,8,then fdele,all,all !删除所有的集中力荷载 nsel,s,loc,x,(ii-1)*deltl !选择集中力作用节点 f,all,FY,-f*cos(w*ii*deltt) !施加集中荷载 allsel !选择所有 solve*elseif,ii,lt,101 fdele,all,all nsel,s,loc,x,(ii-1)*deltl nsel,a,loc,x,(ii-7-1)*deltl f,all,fy,-f*cos(w*ii*deltt) allsel solve*else fdele,all,all nsel,s,loc,x,(ii-7-1)*deltl f,all,fy,-f*cos(w*ii*deltt) allsel solve*endif*enddofinish/post26NSOL,2,nelem/2+1,U,Y,UY_2 !定义第二个变量为UY_2,其值为101号节点的Y方向的位移xvar,1 !定义时间变量为坐标横轴plvar,2 !输出位移图deriv,3,2,1,1 !将位移对时间求导，得到速度，存为变量3plvar,3 !输出速度图fini5.4.2!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。/config,nres,2000/PREP7 !建模ET,1,BEAM3 !梁单元MP,EX,1,207E9 !弹性模量MP,DENS,1,2000 !密度MP,NUXY,1,0.3R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高NELEM=100LBEAM=32 !梁长度deltL=LBEAM/NELEMf=1000V=120 !移动荷载速度NNODE=NELEM+1deltt=deltL/V*3.6*do,i,1,NNODE,1N,i,(i-1)*deltL !节点点矩为0.3*enddo*do,i,1,NELEME,i,i+1*enddofinish/solud,1,alld,NNODE,allantype,trans *do,ii,1,nnode+8,1outres,all,all !输出所有加载子步的所有结果time,ii*deltt !每个荷载步的最终求解时间nsubst,5 !每个荷载步包含五个荷载子步*if,ii,lt,8,then fdele,all,all !删除所有的集中力荷载 nsel,s,