et,1,plane42 keyopt,1,3,2 mp,ex,1,220e9 mp,nuxy,1,0.3 TB,MISO,1,1,25, TBTEMP,0
TBPT,,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.2e6, TBPT,,0.012066,1035.444e6, TBPT,,0.0132,1039.693e6, TBPT,,0.0143,1044.886e6, TBPT,,0.0194,1049.796e6, TBPT,,0.016631,1053.951e6, TBPT,,0.017785,1058.106e6, TBPT,,0.019039,1062.827e6, TBPT,,0.020276,1067.3e6, csys,1
k,1,0.00512 k,2,0.00
kgen,36,1,2,1,,10 *do,j,1,69,2 a,j,j+1,j+3,j+2 *enddo a,71,72,2,1 save
lesize,all,,,4 amesh,all
nummrg,all numcmp,all save
cp,1,uy,166,170,171,172,173 finish /solu
nsel,s,loc,y,270 d,all,all nplot allsel nplot
nsel,s,loc,y,90 d,all,ux nplot allsel
outres,all,all time,1
NSUBST,100 nlgeom,on
d,166,uy,-0.00128 solve
*get,def1,node,166,U,Y nsel,s,loc,y,90 ddele,all,uy allsel f,166,fy,0 solve
*get,def2,node,166,U,Y finish
/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 !显示节点的等效应力 自适应划分 /PREP7 SMRT,OFF
ANTYPE,STATIC
ET,1,PLANE42,,,3 ! 定义平面应力单元有单元厚度 MP,EX,1,210E9 ! 定义材料属性 MP,NUXY,1,.3
R,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,90 K,4,3.25,0.0 L,3,4
CSYS,0 !回到卡式坐标系下。 L,2,4 L,1,3
AL,2,4,1,3 ! 创建面积
DL,4,1,SYMM ! 创建对称面约束 DL,3,1,SYMM
SFL,2,PRES,-10E6 ! 定义线上的分布载荷 MSHK,2 ! 混合网格划分 MSHA,0,2D ! 四边形网格 SAVE ! 备份数据库 FINISH
ADAPT,7,7,,,1 ! 求解次数为nsoln=4次，能量误差百分比stargt=7,单元尺寸改变最大比例因子facmx=1 /POST1
EPLOT ! 绘制单元
PRERR ! 打印能量误差百分比 NSEL,S,LOC,Y,0.0 NSEL,R,LOC,X,2.0
*GET,MNODE,NODE,,NUM,MAX
*GET,SY_D,NODE,MNODE,S,Y ! GET DESIRED STRESS SY VALUE NSEL,ALL ESEL,ALL FINISH 第四章
钢筋混凝土板
!在四个角点增加钢性垫片，在中心点处直接施加-2mm，压碎关闭计算完毕 /clear,nostart /config,nres,5000 /prep7 !定义单元类型
et,1,solid65 !混凝土单元 et,2,solid45 keyopt,1,6,3 keyopt,1,5,1 *afun,deg
R,1,2,0.01,90,0,2,0.01,0,0 !整体模型，实常数 !定义材料属性 mp,ex,1,24000 mp,prxy,1,0.2 tb,conc,1,1,9
tbdata,,0.35,1,3.1125,-1 mp,ex,2,2e5 mp,prxy,2,0.25 tb,bkin,2,1,2,1 tbdata,,360,20000 blc4,,,1000,1000,100 blc4,,,100,100,-50 wpoff,1000,0,0 blc4,,,-100,100,-50 wpoff,0,1000,0 blc4,,,-100,-100,-50 wpoff,-1000,0,0 blc4,,,100,-100,-50 /view,1,1,1,1 /ANG,1,240,ZS,1 /replot /REP,FAST lsel,s,loc,z,50 lesize,all,,,2 lsel,all
lsel,s,loc,x,500 lesize,all,,,20 lsel,all
lsel,s,loc,y,500 lesize,all,,,20 vsel,s,loc,z,50 vmesh,all type,2 mat,2
lsel,s,loc,z,-25 lesize,all,,,1 lsel,all
lsel,s,loc,z,-50 lesize,all,,,2 vsel,s,loc,z,-25 vmesh,all nummrg,all numcmp,all fini
!纵向受拉钢筋，钢支座材料属性
!混凝土材料属性
/solu
asel,s,loc,z,-50 da,all,all nsel,all
nsel,s,loc,z,100 nsel,r,loc,y,500 nsel,r,loc,x,500 d,all,uz,-2 nlgeom,on NSUBST,100 OUTRES,ALL,ALL NEQIT,40 pred,on
CNVTOL,F, ,0.05,2 allsel SOLVE fini /post1
PLNSOL,U,z,0,1 Esel,s,type,,1 PLNSOL,S,1,0,1 PLETAB,REBAR_1,AVG
ETABLE,rebar_2,NMISC,47 !获取rebar2平均等效应力 PLETAB,REBAR_2,AVG etable,seqv_i,nmisc,5 PLETAB,seqv_i,AVG ETABLE,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 PLETAB,11,AVG PLETAB,22,AVG PLETAB,33,NOAV PLETAB,44,AVG PLETAB,55,NOAV /DEVICE,VECTOR,1 PLCRACK,0,0 PLCRACK,0,1
!显示所有开裂位置 !显示第一开裂位置 !获取第一个积分点的状态 !获取第二个积分点的状态 !获取第三个积分点的状态 !获取第四个积分点的状态 !获取第五个积分点的状态 !获取第六个积分点的状态 !获取第七个积分点的状态
!获取solid单元i节点的平均等效应力
!画第一主应力
ETABLE,rebar_1,NMISC,43 !获取rebar1平均等效应力
ETABLE,88,NMISC,102 !获取第八个积分点的状态
钢筋混凝土梁
!采用均匀网格划分，同时取消两种垫片,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.2 tb,conc,1,1,9 tbdata,,0.4,1,3,-1 mp,ex,2,2e5 mp,prxy,2,0.3 tb,bkin,2,1,2,1 tbdata,,350 mp,ex,3,2e5 mp,prxy,3,0.25 tb,bkin,3,1,2,1 tbdata,,200 /pnum,node,1 /pnum,elem,1
!!产生所有的节点**********!!! n,1 n,9,150 fill,1,9
ngen,11,9,1,9,1,,30 ngen,19,1000,1,99,1,,,-75 /view,1,1,1,1 type,2 real,2 mat,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,,170,1000 e,ii,ii+1000 *enddo
!!!!!纵向受拉钢筋************!!!!!!!!!! type,2 real,1 mat,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,1 eplot
!!!**********生成面单元，以便拉伸成体单元!!!!!! k,1 k,2,150, k,3,150,300 k,4,,300 a,1,2,3,4 lsel,s,loc,y,0 lsel,a,loc,y,300
lesize,all,,,8 lsel,all lsel,s,loc,x,0 lsel,a,loc,x,150 lesize,all,,,10 type,3 amesh,all type,1 real,3 mat,1
extopt,esize,18 extopt,aclear,1 vext,all,,,,,-1350 /pnum,mat,1 /pnum,node,0 /pnum,elem,0 /view,1,1,1,1 eplot allsel nummrg,all numcmp,all eplot /solu
nsel,s,loc,y,-0 nsel,r,loc,z,-1275 d,all,uy d,all,ux
asel,s,loc,z,0 da,all,symm nsel,all
!施加外部载荷
nsel,s,loc,y,300 nsel,r,loc,z,-600 d,all,uy,-5 allsel gplot nlgeom,on NSUBST,200 OUTRES,ALL,ALL NEQIT,50 pred,on
CNVTOL,F, ,0.05,2,0.5 allsel SOLVE FINISH
!对称面约束 !右边界约束
!拉伸成混凝土单元
/post1 allsel
/DEVICE,VECTOR,1 PLCRACK,0,1 PLCRACK,0,2 !PLNSOL,U,Y,0,1 !ESEL,S,TYPE,,2 !PLNSOL,S,1,0,1
!显示第一开裂位置 !显示第二开裂位置
!/DEVICE,VECTOR,0
预应力钢筋混凝土
/PREP7 ET,1,LINK8 ET,2,SOLID65 ET,3,PLANE42 R,1,6126 MP,EX,1,2e5 MP,PRXY,1,0.3 MP,ALPX,1,2e-5 MP,EX,2,3.3e MP,PRXY,2,0.2 BLC4, , ,780,210
BLC4,275,210,230,1690 BLC4,-570,1900,1920,200 wpoff,390,191 wprot,0,90 ASBW,1 wprot,0,0,90 ASBW,4 ASBW,5 /PNUM,AREA,1 /PNUM,LINE,1 gplot aglue,all nummrg,all numcmp,all SAVE
lesize,8,,,20 lesize,21,,,8 lesize,22,,,8 lesize,5,,,4 lesize,9,,,2 lesize,7,,,2 lsel,s,,,21,22,1
lsel,a,,,5 LCCAT,all allsel amesh,6 LSEL,s,LCCA lDELE,all allsel
lesize,4,,,20 lesize,6,,,20 lesize,18,,,2 lesize,19,,,2 lsel,s,,,18,19,1 LCCAT,all allsel amesh,4 LSEL,s,LCCA lDELE,all allsel lesize,2,,,6 lesize,3,,,8 lesize,14,,,2 lesize,12,,,2 lsel,s,,,2,18,16 LCCAT,all allsel amesh,2 LSEL,s,LCCA lDELE,all allsel lesize,20,,,6 lesize,17,,,8 lesize,13,,,2 lsel,s,,,19,20,1 LCCAT,all allsel amesh,5 LSEL,s,LCCA lDELE,all allsel lesize,10,,,6 lesize,15,,,6 lesize,11,,,6 lesize,1,,,8 lesize,16,,,8
allsel amesh,1 amesh,3 SAVE
TYPE, 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,1 REAL,1 lmesh,29 FINISH /SOL DL,23, ,ALL, DL,38, ,ALL, DL,1, ,UY, DL,16, ,UY, DK,1, , , ,0,UX,UY SAVE
BFL,29,TEMP,-375 CNVTOL,F, ,0.05,2, , NSUBST,10,0,0 OUTRES,ERASE OUTRES,ALL,ALL TIME,5 SOLVE
4.6混凝土开裂模拟 /TITLE,Crack /PREP7 ET,1,SOLID65 KEYOPT,1,5,2 ET,2,LINK8 ET,3,PIPE16 R,1 R,2,0.15 R,3,1,0.5 MP,EX,1,2E6
MP,PRXY,1,0 TB,CONC,1 TBDATA,,,,,-1 MP,EX,2,30E6 MP,PRXY,2,0.3
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> BLC4, , ,5,6,1.5 /VIEW, 1 ,1,1,1 /ANG, 1 /REP,FAST
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /PNUM,LINE,1 LSEL,,LOC,Z,0.1,1.4 LSEL,R,LOC,Y,0 TYPE, 2 MAT, 2 REAL, 2 LESIZE,ALL,,,1 LMESH,ALL ALLSEL
LSEL,,LOC,Z,1.5 LSEL,R,LOC,Y,1,4 TYPE, 3 MAT, 1 REAL, 3 LESIZE,ALL,,,4 LMESH,ALL ALLSEL TYPE, 1 MAT, 1 REAL, 1 LSEL,,LOC,X,1,4 LESIZE,ALL,,,1 ALLSEL MSHAPE,0,3D MSHKEY,1 VMESH,ALL FINISH
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /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,0 D,ALL,ALL ALLSEL
NSEL,,LOC,Z,1.5 D,ALL,ROTY D,ALL,ROTZ ALLSEL F,7,MX,300 F,11,MX,300 /SOL
NSUBST,5,0,0 AUTOTS,1 SOLVE
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /POST1 SET,LAST
*GET,SCON,NODE,17,S,Z ETAB,ST,LS,1 ESORT,ST
*GET,STL,SORT,,MAX *status,parm /DEVICE,VECTOR,1 PLCRACK SAVE FINISH
4.7混凝土浇筑模拟 /TITLE,CASTCONCRETE /PREP7 ET,1,SOLID70 MP,KXX,1,300. MP,C,1,0.94500 MP,DENS,1,2400 MP,KXX,2,300. 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,0
LSEL,R,LOC,X,1,49 LESIZE,ALL,,,10 ALLSEL
LSEL,,LOC,Z,40 LSEL,R,LOC,X,1,49 LESIZE,ALL,,,10 ALLSEL LSEL,,LOC,X,0 LSEL,R,LOC,Z,1,39 LESIZE,ALL,,,10 ALLSEL
LSEL,,LOC,X,50 LSEL,R,LOC,Z,1,39 LESIZE,ALL,,,10 ALLSEL LSEL,,LOC,Z,0 LSEL,R,LOC,Y,1,1.5 LESIZE,ALL,,,2 ALLSEL
LSEL,,LOC,Z,40 LSEL,R,LOC,Y,1,1.5 LESIZE,ALL,,,2 ALLSEL
LSEL,,LOC,Z,-20 LSEL,R,LOC,X,1,49 LESIZE,ALL,,,20 ALLSEL
LSEL,,LOC,Z,60 LSEL,R,LOC,X,1,49 LESIZE,ALL,,,20 ALLSEL
LSEL,,LOC,X,-25 LSEL,R,LOC,Z,1,39 LESIZE,ALL,,,20 ALLSEL
LSEL,,LOC,X,75 LSEL,R,LOC,Z,1,39 LESIZE,ALL,,,20 ALLSEL
LSEL,,LOC,Z,-20 LSEL,R,LOC,Y,-9,-1.3 LESIZE,ALL,,,10 ALLSEL
LSEL,,LOC,Z,60
LSEL,R,LOC,Y,-9,-1.3 LESIZE,ALL,,,10 ALLSEL
VSEL,,LOC,Y,1,1.5 MAT,1 VMESH,1 ALLSEL
VSEL,,LOC,Y,-9,-1.5 MAT,2 VMESH,2 ALLSEL NUMMRG,ALL NUMCMP,NODE NROPT,FULL ESEL,S,MAT,,2 NSLE,S IC,ALL,TEMP,16.5 ASEL,,LOC,X,-25 ASEL,A,LOC,X,75 ASEL,A,LOC,Y,-10 ASEL,A,LOC,Z,-20 ASEL,A,LOC,Z,60 NSLA,S
SF,ALL,HFLUX,0 ALLSEL
NSEL,,LOC,X,49,76 NSEL,R,LOC,Y,0 NSEL,A,LOC,X,-26,1 NSEL,R,LOC,Y,0 NSEL,A,LOC,Z,39,61 NSEL,R,LOC,Y,0 NSEL,A,LOC,Z,-21,1 NSEL,R,LOC,Y,0 SF,ALL,CONV,2016,20 ESEL,,MAT,,1 EKILL,ALL NSLE,S
IC,ALL,TEMP,40 NSEL,,LOC,Y,0,1 ESLN,S,1 EALIVE,ALL ALLSEL
NSEL,,LOC,X,1,49 NSEL,R,LOC,Z,1,39
NSEL,R,LOC,Y,0 CM,DD,NODE ALLSEL
NSEL,,LOC,Y,0,1 ESLN,S,1 NSLE,S
NSEL,U,,,DD SF,ALL,CONV,2016,20 HE00=47880/1.862 BFE,ALL,HGEN,,HE00 ALLSEL /SOLU ANTYPE,TRANS TIME,1 AUTOTS,0 DELTIM,1 KBC,0 SOLVE /POST1 PLNSOL,TG,SUM /PREP7
NSEL,,LOC,Y,0,1 ESLN,S,1
HE00=47880*(2/2.862- 1/1.862) BFE,ALL,HGEN,,HE00 ALLSEL /SOLU ANTYPE,TRANS TIME,2 AUTOTS,0 DELTIM,1 KBC,0 SOLVE /POST1 PLNSOL,TG,SUM /PREP7 ESEL,,MAT,,1 NSLE,S
SFEDELE,ALL,ALL,CONV ALLSEL
NSEL,,LOC,X,1,49 NSEL,R,LOC,Z,1,39 NSEL,R,LOC,Y,0,1 CM,KK,NODE
NSEL,,LOC,Y,0,2 ESLN,S,1 NSLE,S
NSEL,U,,,KK SF,ALL,CONV,2016,20 ALLSEL
NSEL,,LOC,Y,0,1 ESLN,S,1
HE00=47880*(3/3.862-2/2.862) BFE,ALL,HGEN,,HE00 ALLSEL
NSEL,,LOC,Y,1,2 ESLN,S,1
HE00=47880/1.862 BFE,ALL,HGEN,,HE00 ALLSEL /SOLU ANTYPE,TRANS TIME,3 AUTOTS,0 DELTIM,1 KBC,0 SOLVE /POST1 PLNSOL,TG,SUM /PREP7
NSEL,,LOC,Y,0,1 ESLN,S,1
HE00=47880*(4/4.862-3/3.862) BFE,ALL,HGEN,,HE00 ALLSEL
NSEL,,LOC,Y,1,2 ESLN,S,1
HE00=47880*(2/2.862-1/1.862) BFE,ALL,HGEN,,HE00 ALLSEL /SOLU ANTYPE,TRANS TIME,4 AUTOTS,0 DELTIM,1 KBC,0 SOLVE /POST1
PLNSOL,TG,SUM SAVE FINISH
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /PREP7 ET,1,SOLID45 MP,EX,1,0.325E10 MP,DENS,1,2400 MP,PRXY,1,0.2 MP,ALPX,1,0.9E-5 MP,REFT,1,15 MP,EX,2,2.08E10 MP,DENS,2,2600 MP,PRXY,2,0.2 MP,ALPX,2,0.9E-5 ESEL,,MAT,,2 MPCHG,2,ALL ALLSEL /SOLU
ASEL,,LOC,X,-25 ASEL,A,LOC,X,75 ASEL,A,LOC,Y,-10 ASEL,A,LOC,Z,-20 ASEL,A,LOC,Z,60 NSLA,S D,ALL,ALL ALLSEL ESEL,,MAT,,1 EKILL,ALL NSEL,,LOC,Y,0,1 ESLN,S,1 EALIVE,ALL MPCHG,1,ALL ALLSEL ACEL,0,9.8,0
LDREAD,TEMP,,,1,,'CASTCONCRETE','rth','' ANTYPE,TRANS NLGEOM,ON
CNVTOL,F,,0.2,2,200 TIME,1 AUTOTS,0 DELTIM,1 KBC,0 SOLVE
/SOLU
LDREAD,TEMP,,,2,,'CASTCONCRETE','rth','' TIME,2 SOLVE /SOLU
NSEL,,LOC,Y,1,2 ESLN,S,1 EALIVE,ALL MPCHG,1,ALL ALLSEL
LDREAD,TEMP,,,3,,'CASTCONCRETE','rth','' TIME,3 SOLVE /SOLU
LDREAD,TEMP,,,4,,'CASTCONCRETE','rth','' TIME,4 SOLVE FINISH SAVE
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /POST1 PLNSOL,U,Y PLNSOL,S,Y PLESOL,S,Y 5.2行架桥 /TITLE,TRUSS /PREP7 ET,1,BEAM4 ET,2,LINK8 ET,3,SHELL63
R,1,3,0.25,2.25,3,1,,-0.0001,1 R,2,3,0.25,2.25,3,1,,,1 R,3,2,1/6,2/3,2,1,,,2/3 R,4,1
R,5,1,-1.5E-4 R,6,0.3 MP,EX,1,35E10 MP,DENS,1,3000 MP,PRXY,1,1/6 MP,EX,2,21E10 MP,DENS,2,3500 MP,PRXY,2,0.3 MP,EX,3,3.5E10 MP,DENS,3,2500
MP,PRXY,3,1/6 MP,EX,4,3.5E10 MP,DENS,4,2500 MP,PRXY,4,1/6
!>>>>>>>>>>>>>>>>>>>>> N,1,0,-5,0 NGEN,4,1,1,,,12 N,5,0,-5,12 NGEN,3,1,5,,,12 NSYM,Y,7,ALL TYPE,1 REAL,1 MAT,1 E,1,2 E,2,3 E,3,4 E,8,9 E,9,10 E,10,11 TYPE,1 REAL,2 MAT,1 E,5,6 E,6,7 E,12,13 E,13,14 E,5,12 E,6,13 E,7,14 TYPE,1 REAL,3 MAT,2 E,7,4 E,14,11 TYPE,2 REAL,4 MAT,3 E,5,2 E,2,7 E,12,9 E,9,14 TYPE,2 REAL,5 MAT,3
E,5,1 E,6,2 E,7,3 E,12,8 E,13,9 E,14,10 TYPE,3 REAL,6 MAT,4 E,1,2,9,8 E,2,3,10,9 E,3,4,11,10 NSYM,X,14,ALL ESYM,,14,ALL NUMMRG,ALL NUMCMP,ALL SAVE /SOLU D,17,UX D,17,UY D,17,UZ D,22,UX D,22,UY D,22,UZ D,4,UY D,4,UZ D,11,UY D,11,UZ ACEL,,,10 F,1,FZ,-10000 F,8,FZ,-10000 SAVE SOLVE
!>>>>>>>>>>>>>>>>>>>> /POST1 PLVECT,U PRNSOL,U,Z PRESOL,ELEM FINISH SAVE 5.3悬索桥 /TITLE,CABLE /PREP7 ET,1,LINK10
ET,2,BEAM3 R,1,1,1E-5 R,2,0.02,1E-5 R,3,0.5,1/24,1 MP,EX,1,25E10 MP,PRXY,1,0.1 MP,DENS,1,1000 MP,EX,2,30E10 MP,PRXY,2,0.1 MP,DENS,2,1000 A=250 *DO,I,1,10 X=10*(I-1) Y=A*COSH(X/A)-A N,I,X,Y *ENDDO N,11,100,20 *DO,I,1,10 X=10*I Y=0 N,11+I,X,Y *ENDDO type,1 real,1 mat,1 *DO,I,1,10 E,I,I+1 *ENDDO type,1 real,2 mat,1 *DO,I,2,10 E,I,I+10 *ENDDO type,2 real,3 mat,2 E,1,12 *DO,I,1,9 E,I+11,I+12 *ENDDO NSYM,X,21,ALL ESYM,,21,ALL NUMMRG,ALL
NUMCMP,ALL SAVE
!>>>>>>>>>>>>>>>>>>>> /SOLU D,41,UX D,41,UY D,21,UY D,31,UX D,31,UY D,11,UX D,11,UY
NSEL,U,NODE,,41 NSEL,U,NODE,,21 D,ALL,UZ ALLSEL ACEL,,10,, NROPT,FULL ESEL,,REAL,,2 ESEL,A,REAL,,3 EKILL,ALL ALLSEL NLGEOM,ON TIME,1
NSUBST,10,20,5 OUTRES,ERASE OUTRES,ALL,ALL SSTIF,ON SOLVE /POST1 ESEL,,LIVE PLVECT,U ALLSEL /PREP7
UPGEOM,1,LAST,LAST,'Cable','RST' /SOLU SOLVE /POST1 ESEL,,LIVE PLVECT,U ALLSEL SAVE /SOLU
F,1,FY,-50000 F,2,FY,-25000
F,22,FY,-25000 TIME,2 SOLVE /POST1 ESEL,,LIVE PLVECT,U ALLSEL /SOLU DOFSEL,S,,FY FCUM,ADD,1 DOFSEL,ALL *DO,I,1,9 XX0=I*10 XX=(I+1)*10 NSEL,,LOC,X,-XX0 NSEL,R,LOC,Y,0.1,30 NSEL,A,LOC,X,XX0 NSEL,R,LOC,Y,0.1,30 F,ALL,FY,-25000 ALLSEL
NSEL,,LOC,X,XX NSEL,R,LOC,Y,0.1,30 NSEL,A,LOC,X,-XX NSEL,R,LOC,Y,0.1,30 F,ALL,FY,-25000 ALLSEL TIME,I+2 SOLVE *ENDDO /POST1 ESEL,,LIVE PLVECT,U ALLSEL SAVE
!>>>>>>>>>>>>>>>>>> /SOLU
NSEL,,LOC,Y,0.1,30 NSEL,A,NODE,,1 FDELE,ALL,FY ALLSEL ESEL,U,LIVE EALIVE,ALL ALLSEL TIME,12
SOLVE /POST1 ESEL,,LIVE PLVECT,U ALLSEL FINISH SAVE 5.4.1
!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。
!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。
/config,nres,2000
/PREP7 !建模 ET,1,BEAM3 !梁单元 MP,EX,1,207E9 !弹性模量 MP,DENS,1,2000 !密度 MP,NUXY,1,0.3
R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高 NELEM=100
LBEAM=32 !梁长度 deltL=LBEAM/NELEM V=120 f=1000
w=10 !移动荷载速度 NNODE=NELEM+1 deltt=deltL/V*3.6 *do,i,1,NNODE,1
N,i,(i-1)*deltL !节点点矩为0.3 *enddo *do,i,1,NELEM E,i,i+1 *enddo finish /solu d,1,all d,NNODE,all antype,trans *do,ii,1,101+8,1 outres,all,all
!输出所有加载子步的所有结果 !每个荷载步的最终求解时间
time,ii*deltt *if,ii,lt,8,then
fdele,all,all
!删除所有的集中力荷载
nsubst,5 !每个荷载步包含五个荷载子步
nsel,s,loc,x,(ii-1)*deltl f,all,FY,-f*cos(w*ii*deltt) 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 *enddo finish /post26
!选择集中力作用节点
!施加集中荷载
allsel !选择所有
NSOL,2,nelem/2+1,U,Y,UY_2 !定义第二个变量为UY_2,其值为101号节点的Y方向的位移 xvar,1 !定义时间变量为坐标横轴 plvar,2 !输出位移图
deriv,3,2,1,,,,,1 !将位移对时间求导，得到速度，存为变量3 plvar,3 !输出速度图 fini 5.4.2
!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。
!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。 /config,nres,2000
/PREP7 !建模 ET,1,BEAM3 !梁单元 MP,EX,1,207E9 !弹性模量 MP,DENS,1,2000 !密度 MP,NUXY,1,0.3
R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高 NELEM=100
LBEAM=32 !梁长度 deltL=LBEAM/NELEM f=1000
V=120 !移动荷载速度 NNODE=NELEM+1
deltt=deltL/V*3.6 *do,i,1,NNODE,1
N,i,(i-1)*deltL !节点点矩为0.3 *enddo *do,i,1,NELEM E,i,i+1 *enddo finish /solu d,1,all d,NNODE,all antype,trans *do,ii,1,nnode+8,1 outres,all,all
!输出所有加载子步的所有结果 !每个荷载步的最终求解时间
time,ii*deltt *if,ii,lt,8,then
fdele,all,all f,all,FY,-f solve
*elseif,ii,lt,nnode fdele,all,all
nsel,s,loc,x,(ii-1)*deltl nsel,a,loc,x,(ii-7-1)*deltl f,all,fy,-f allsel solve *else
fdele,all,all
nsel,s,loc,x,(ii-7-1)*deltl f,all,fy,-f allsel solve *endif *enddo finish /post26
NSOL,2,nelem/2+1,U,Y,UY_2 !定义第二个变量为UY_2,其值为101号节点的Y方向的位移 xvar,1 !定义时间变量为坐标横轴 plvar,2 !输出位移图
deriv,3,2,1,,,,,1 !将位移对时间求导，得到速度，存为变量3 plvar,3 !输出速度图
!删除所有的集中力荷载
!选择集中力作用节点
nsel,s,loc,x,(ii-1)*deltl
nsubst,5 !每个荷载步包含五个荷载子步
!施加集中荷载
allsel !选择所有
fini 5.5 /prep7
/VIEW, 1 ,1,1,1 et,1,shell63 et,2,link8
r,1,0.25 !顶层与底层的厚度 r,2,0.6 !腹板的厚度 r,3,1 !V形支撑厚度
r,4,0.02,-0.005 !预应力钢筋的面积初始与初始应变 mp,ex,1,3.25e10 mp,prxy,1,0.2 mp,dens,1,2700 mp,ex,2,207e9 mp,prxy,2,0.3 mp,dens,2,7800 !桥面关键点 k,1 k,2,4 k,3,10 k,4,16 k,5,20 k,6,24 k,7,30 k,8,36 k,9,40 SAVE
!箱形梁底面关键点 k,10,4,-1.8 k,11,10,-1.8 k,12,16,-1.8 k,13,24,-1.8 k,14,30,-1.8 k,15,36,-1.8 SAVE SAVE
kgen,2,1,15,1,,,-16 kgen,2,1,15,1,,,-85 kgen,2,1,15,1,,,-1 kgen,2,1,15,1,,,-170 SAVE
k,101,,,-36 k,102,4,,-36 k,103,10,,-36
k,104,16,,-36 k,105,20,,-36 k,106,24,,-36 k,107,30,,-36 k,108,36,,-36 k,109,40,,-36
kgen,2,102,104,1,,-3,,8 kgen,2,106,108,1,,-3,,7 SAVE
kgen,2,101,115,1,,,-12,15 kgen,2,116,130,1,,,-12,15 kgen,2,131,145,1,,,-50,15 kgen,3,146,160,1,,,-12,15 SAVE
!顶层面建模并网格划分 *do,i,1,8,1 a,i,i+15,i+16,i+1 *enddo *do,i,16,23,1 a,i,i+85,i+86,i+1 *enddo
*do,i,101,108,1 a,i,i+15,i+16,i+1 *enddo
*do,i,116,123,1 a,i,i+15,i+16,i+1 *enddo *do,i,31,38,1 a,i,i+100,i+101,i+1 *enddo *do,i,31,38,1 a,i,i+115,i+116,i+1 *enddo
*do,i,146,153,1 a,i,i+15,i+16,i+1 *enddo
*do,i,161,168,1 a,i,i+15,i+16,i+1 *enddo *do,i,46,53,1 a,i,i+130,i+131,i+1 *enddo *do,i,46,53,1 a,i,i+15,i+16,i+1
*enddo save
lesize,all,4 aatt,1,1,1 amesh,all
!lsel,s,loc,x,-0.1,0.1 !lsel,a,loc,x,41.9,42.1 !lesize,all,,,10 !allsel !lplot
!lsel,s,loc,z,-0.1,0.1 !lsel,a,loc,z,-165,-166 !lesize,all,,,2 !lplot
asel,u,loc,y,-0.1,0.1 aplot !底层
*do,i,10,11,1 a,i,i+15,i+16,i+1 *enddo *do,i,13,14,1 a,i,i+15,i+16,i+1 *enddo *do,i,25,26,1 a,i,i+85,i+86,i+1 *enddo *do,i,28,29,1 a,i,i+85,i+86,i+1 *enddo
*do,i,110,111,1 a,i,i+15,i+16,i+1 *enddo
*do,i,113,114,1 a,i,i+15,i+16,i+1 *enddo
*do,i,125,126,1 a,i,i+15,i+16,i+1 *enddo
*do,i,128,129,1 a,i,i+15,i+16,i+1 *enddo save
*do,i,40,41,1 a,i,i+100,i+101,i+1
排除顶面面元素!
*enddo *do,i,43,44,1 a,i,i+100,i+101,i+1 *enddo *do,i,40,41,1 a,i,i+115,i+116,i+1 *enddo *do,i,43,44,1 a,i,i+115,i+116,i+1 *enddo
*do,i,155,156,1 a,i,i+15,i+16,i+1 *enddo
*do,i,158,159,1 a,i,i+15,i+16,i+1 *enddo
*do,i,170,171,1 a,i,i+15,i+16,i+1 *enddo
*do,i,173,174,1 a,i,i+15,i+16,i+1 *enddo *do,i,55,56,1 a,i,i+130,i+131,i+1 *enddo *do,i,58,59,1 a,i,i+130,i+131,i+1 *enddo *do,i,55,56,1 a,i,i+15,i+16,i+1 *enddo *do,i,58,59,1 a,i,i+15,i+16,i+1 *enddo save
lesize,all,4 aatt,1,1,1 amesh,all
asel,u,loc,z,-170.1,0.1 /pnum,kp,1 kplot aplot save !腹板
*do,i,2,4,1 a,i,i+8,i+23,i+15 *enddo *do,i,6,8,1 a,i,i+7,i+22,i+15 *enddo *do,i,17,19,1 a,i,i+8,i+93,i+85 *enddo *do,i,21,23,1 a,i,i+7,i+92,i+85 *enddo
*do,i,102,104,1 a,i,i+8,i+23,i+15 *enddo
*do,i,106,108,1 a,i,i+7,i+22,i+15 *enddo
*do,i,117,119,1 a,i,i+8,i+23,i+15 *enddo
*do,i,121,123,1 a,i,i+7,i+22,i+15 *enddo *do,i,32,34,1 a,i,i+8,i+108,i+100 *enddo *do,i,36,38,1 a,i,i+7,i+107,i+100 *enddo *do,i,32,34,1 a,i,i+8,i+123,i+115 *enddo *do,i,36,38,1 a,i,i+7,i+122,i+115 *enddo
*do,i,147,149,1 a,i,i+8,i+23,i+15 *enddo
*do,i,151,153,1 a,i,i+7,i+22,i+15 *enddo
*do,i,162,1,1 a,i,i+8,i+23,i+15
*enddo
*do,i,166,168,1 a,i,i+7,i+22,i+15 *enddo *do,i,47,49,1 a,i,i+8,i+138,i+130 *enddo *do,i,51,53,1 a,i,i+7,i+137,i+130 *enddo *do,i,47,49,1 a,i,i+8,i+23,i+15 *enddo *do,i,51,53,1 a,i,i+7,i+22,i+15 *enddo save aplot
!!先对y方向杆设置划分数目为1和2份，然后在对Z!方向的杆设置单元长度为4 lsel,s,loc,z,-0.1,0.1 lsel,a,loc,z,-17,-15 lsel,a,loc,z,-86,-84 lsel,a,loc,z,-155,-153 lsel,a,loc,z,-171,-169 lesize,all,,,2 lsel,all
lsel,s,loc,z,-37,-35 lsel,a,loc,z,-61,-59 lsel,a,loc,z,-111,-109 lsel,a,loc,z,-135,-133 lesize,all,,,2 !esize,all,4 aatt,1,2,1 amesh,all save
!!!!!!!!!!!!!!!!!!!!!!!!!!!!预应力杆设置 lsel,s,loc,x,3.9,4.1 lsel,a,loc,x,9,11 lsel,a,loc,x,15,17 lsel,a,loc,x,23,25 lsel,a,loc,x,29,31 lsel,a,loc,x,35,37 lsel,u,loc,y,-1.7,-0.1 lplot
lesize,all,4 latt,2,4,2 lmesh,all save !V型支座
asel,u,loc,y,-3.1,0.1 aplot
k,1001,4,-15,-48 k,1002,16,-15,-48 k,1003,24,-15,-48 k,1004,36,-15,-48 k,2001,4,-15,-122 k,2002,16,-15,-122 k,2003,24,-15,-122 k,2004,36,-15,-122 a,110,112,1002,1001 a,112,113,1003,1002 a,113,115,1004,1003 a,140,142,1002,1001 a,142,143,1003,1002 a,143,145,1004,1003 a,155,157,2002,2001 a,157,158,2003,2002 a,158,160,2004,2003 a,185,187,2002,2001 a,187,188,2003,2002 a,188,190,2004,2003 save
lesize,all,4 aatt,1,3,1 amesh,all save allsel nummrg,all numcmp,all /eshape,1 eplot !施加约束
nsel,s,loc,z,-0.1,0.1 nsel,a,loc,z,-171,-169 nsel,r,loc,y,-1.9,-1.7 nplot d,all,uy allsel
nsel,s,loc,y,-16,-14 nplot d,all,all allsel gplot /solu acel,,9.8 pstres,on allsel solve finish /post1 plnsol,u,y /eshape,0
/DSCALE,1,user,10 ETABLE, ,LS,1 PLETAB,LS1,NOAV fini !模态求解 /solu antype,2 upcoord,1,on pstres,on modopt,subsp,10 mxpand,10,,,1 allsel solve FINISH !后处理 /POST1 SET,LIST SET,FIRST PLDISP,0 !SET,NEXT !PLDISP,0 !SET,NEXT !PLDISP,0 !finish 5.7
/TITLE,cable-stayed bridge /prep7 et,1,beam4 et,2,link10
r,1,25.6,20,2000,16,1.6,,,21.8
r,2,16,39.7,77.7,4.7,3.4,,,61.6 r,3,40,213.3,83.3,5,8,,,333.3 r,4,,450,200,6,9,,,8
r,5,1,0.08333333,0.08333333,1,1,,,1/3 r,6,0.12,0 mp,ex,1,3.5e10 mp,prxy,1,0.17 mp,dens,1,2500 mp,ex,2,10e15 mp,prxy,2,0 mp,dens,2,0 mp,ex,3,1.9e10 mp,prxy,3,0.25 mp,dens,3,1200 *do,i,1,59 x=-174*2+(i-1)*6 y1=-14 y2=14
n,3*(i-1)+1,x ! 1,4,7,... n,3*(i-1)+2,x,y1 ! 2,5,8,... n,3*i,x,y2 ! 3,6,9,... *end do type,1 real,1 mat,1 *do,i,1,58 j=3*(i-1)+1 e,j,j+3 *end do type,1 real,5 mat,2 *do,i,1,59 j=3*(i-1)+1 j1=3*(i-1)+2 j2=3*i e,j,j1 e,j,j2 *end do i=59*3
n,i+1,-174,-10,-30 n,i+2,-174,10,-30 n,i+3,-174,-15 n,i+4,-174,15
*do,j,1,5 k=i+4+j
n,k,-174,0,60+(j-1)*18 *end do type,1 real,4 mat,1 e,i+1,i+3 e,i+2,i+4 type,1 real,3 mat,1 e,i+3,i+5 e,i+4,i+5 type,1 real,2 mat,1 *do,j,1,4 k=i+4+j e,k,k+1 *end do
!>>>>>>>>>>>>>>>>>>>>>>>>>> type,2 real,6 mat,3 e,i+5, e,i+5,90 type,2 real,6 mat,3 *do,j,1,8 e,i+6,+3*j e,i+6,-3*j e,i+6,90+3*j e,i+6,90-3*j *end do *do,k,1,3 *do,j,1,7
e,i+6+k,113+(k-1)*21+3*j e,i+6+k,65-(k-1)*21-3*j e,i+6+k,114+(k-1)*21+3*j e,i+6+k,66-(k-1)*21-3*j *end do *end do
i=i+9 nsym,x,i,all esym,,i,all nummrg,all type,1 real,5 mat,2 j=i i=i-9 n,1000,-174 e,1000,i+3 e,1000,i+4 n,2000,174 e,2000,i+3+j e,2000,i+4+j cpintf,uy cpintf,uz cpintf,rotx cpintf,rotz nsel,,loc,z,-30 d,all,all allsel
nsel,,loc,x,-348 nsel,r,loc,y,0 d,all,uy d,all,uz allsel
nsel,,loc,x,348 nsel,r,loc,y,0 d,all,uy allsel NUMCMP,all /solu solve /prep7
r,6,0.012,0.012 /solu lumpm,on pstres,on solve finish
upcoord,1.0,on save /solu
antype,modal modopt,lanb,20 mxpand,10,0,100 lumpm,on pstres,on solve finish
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> /resume
!将生成的地震数据文件放到工作目录下 *dim,AAZ1,array,2,50 *dim,AAY1,array,2,50 *CREATE,ansuitmp
*vread,AAZ1(1,1),ACELZ,txt,,50 (e9.3 , e11.3) *END
/INPUT,ansuitmp *CREATE,ansuitmp
*vread,AAY1(1,1),ACELY,txt,,50 (e9.3 , e11.3) *END
/INPUT,ansuitmp N=117
*dim,fcar,array,N *do,i,1,N tt=(i-1)*0.1
fcar(i)=1000*cos(10*tt) *enddo
*dim,fwind,array,N,N *do,i,1,N tt=(i-1)*0.1 p=50*sin(1.5*tt) fwind(1,i)=3*3*p *do,j,1,N-2
fwind(j+1,i)=6*3*p *enddo
fwind(N,i)=fwind(1,i) *enddo FINISH
/config,nres,2000 /solu antype,trans M=20 *do,i,1,M
time,i*0.1 auto,off deltim,0.05 kbc,0 outres,basic cnvtol,f,,0.01 *if,i,le,50,then
acel,,-AAY1(2,i),10-AAZ1(2,i) *ELSE
ACEL,0,0,10 *endif
*if,i,eq,1,then nsel,s,loc,x,-174*2 nsel,r,loc,y,0 nsel,r,loc,Z,0 f,all,fz,-fcar(j) allsel *else *do,j,2,M
nsel,s,loc,x,-174*2+(j-1)*6 nsel,r,loc,y,0 nsel,r,loc,Z,0 f,all,fz,-fcar(j) allsel
nsel,s,loc,x,-174*2+(j-2)*6 nsel,r,loc,y,0 nsel,r,loc,Z,0 fdele,all,fz allsel *enddo *endif *do,j,1,N
nsel,s,loc,y,-14
nsel,s,loc,x,-174*2+(j-1)*6 f,all,fy,fwind(j,i) allsel *enddo lswrite,i *enddo lssolve,1,M finish 5.7acely
1.000E-01 2.597E-02 2.000E-01 1.065E-01
3.000E-01 1.134E-01 4.000E-01 2.142E-02 5.000E-01 -2.171E-01 6.000E-01 -9.088E-02 7.000E-01 2.191E-01 8.000E-01 3.830E-01 9.000E-01 1.2E-01 1.000E+00 1.537E-01 1.100E+00 -1.436E-01 1.200E+00 1.553E-01 1.300E+00 7.4E-02 1.400E+00 2.819E-01 1.500E+00 8.650E-02 1.600E+00 2.937E-02 1.700E+00 1.975E-01 1.800E+00 4.950E-02 1.900E+00 4.609E-02 2.000E+00 -3.230E-01 2.100E+00 5.031E-02 2.200E+00 -2.851E-01 2.300E+00 -1.761E-01 2.400E+00 -1.056E-01 2.500E+00 1.240E-01 2.600E+00 2.218E-01 2.700E+00 2.256E-01 2.800E+00 -2.678E-01 2.900E+00 1.386E-01 3.000E+00 -5.203E-01 3.100E+00 4.560E-02 3.200E+00 -7.798E-01 3.300E+00 -3.473E-02 3.400E+00 1.902E-01 3.500E+00 2.131E-01 3.600E+00 -4.8E-01 3.700E+00 -3.686E-01 3.800E+00 1.171E+00 3.900E+00 2.259E-01 4.000E+00 6.5E-02 4.100E+00 -5.975E-01 4.200E+00 1.457E-01 4.300E+00 -3.363E-01 4.400E+00 -4.106E-02 4.500E+00 2.952E-01 4.600E+00 -3.781E-02
4.700E+00 1.349E-01 4.800E+00 -3.978E-01 4.900E+00 1.107E-01 5.000E+00 5.493E-01 5.7acelz
1.000E-01 1.298E-02 2.000E-01 5.323E-02 3.000E-01 5.672E-02 4.000E-01 1.071E-02 5.000E-01 -1.086E-01 6.000E-01 -4.4E-02 7.000E-01 1.095E-01 8.000E-01 1.915E-01 9.000E-01 9.461E-02 1.000E+00 7.684E-02 1.100E+00 -7.181E-02 1.200E+00 7.765E-02 1.300E+00 3.822E-02 1.400E+00 1.409E-01 1.500E+00 4.325E-02 1.600E+00 1.469E-02 1.700E+00 9.875E-02 1.800E+00 2.475E-02 1.900E+00 2.304E-02 2.000E+00 -1.615E-01 2.100E+00 2.515E-02 2.200E+00 -1.426E-01 2.300E+00 -8.804E-02 2.400E+00 -5.282E-02 2.500E+00 6.199E-02 2.600E+00 1.109E-01 2.700E+00 1.128E-01 2.800E+00 -1.339E-01 2.900E+00 6.930E-02 3.000E+00 -2.601E-01 3.100E+00 2.280E-02 3.200E+00 -3.9E-01 3.300E+00 -1.736E-02 3.400E+00 9.510E-02 3.500E+00 1.065E-01 3.600E+00 -2.427E-01 3.700E+00 -1.843E-01 3.800E+00 5.856E-01 3.900E+00 1.130E-01
4.000E+00 3.294E-02 4.100E+00 -2.988E-01 4.200E+00 7.287E-02 4.300E+00 -1.681E-01 4.400E+00 -2.053E-02 4.500E+00 1.476E-01 4.600E+00 -1.1E-02 4.700E+00 6.743E-02 4.800E+00 -1.9E-01 4.900E+00 5.534E-02 5.000E+00 2.747E-01 6框架模态分析 /prep7
et,1,beam188 !梁的单元类型
et,2,shell63 !筒体,楼板及外墙单元类型 mp,ex,1,3.25e10 !C40混凝土 mp,prxy,1,0.2 !泊松比 mp,dens,1,2700 !密度 mp,ex,2,3.0e10 !C30混凝土 mp,prxy,2,0.2 !泊松比 mp,dens,2,2700 !密度 sectype,1,beam,rect !框架柱截面 secdata,1.1,1.1
sectype,2,beam,rect !外环梁截面 secdata,0.4,0.6
sectype,3,beam,rect !内框架梁截面 secdata,0.5,0.8
sectype,4,beam,rect !次梁 secdata,0.3,0.5
r,2,0.2 !楼板厚度 r,3,0.3 !筒体厚度
k,5000,22,8,72 !设置beam188单元节点的方向 k,1 k,12,44 kfill,1,12 kgen,5,1,12,,,4 kgen,19,1,60,,,,3 /view,1,1,1,1 /ANGLE,1,270,XM,0 /REPLOT !框架柱共288根 *do,i,1,1021,60 l,i,i+60 *enddo
lgen,3,1,18,1,8,,,2 lgen,2,1,,1,,16,,48 lgen,2,1,18,1,,8,,24 lgen,2,37,,1,,8,,24 lgen,2,1,144,1,28,,,7 latt,1,,1,,5000,,1 lesize,all,1.5 lmesh,all lsel,u,,,all
!线单元的最大编号为288 !外环梁KL2共144根 l,61,65 l,65, l,,92 l,92,68 l,68,72 l,72,120 l,120,109 l,109,61
lgen,18,2,296,1,,,3,60 latt,1,,1,,5000,,2 lesize,all,2 lmesh,all lsel,u,,,all
!线单元的最大编号为432 !内框架梁KL1共108根 l,63,111 l,85, l,,113 l,70,118 l,92,96 l,92,116
lgen,18,433,438,1,,,3,60 latt,1,,1,,5000,,3 lesize,all,2 lmesh,all lsel,u,,,all
!线单元的最大编号为0 !次梁共162根 l,62,110 l,,112 l,69,117 l,71,119 l,97,108
l,73,77 l,80,84 l,90,102 l,91,103
lgen,18,1,9,1,,,3,60 latt,1,,1,,5000,,4 lesize,all,2 lmesh,all allsel
!线单元的最大编号为702 !楼层板 *do,i,61,,1 a,i,i+1,i+13,i+12 *enddo
agen,4,1,4,1,,4,,12 *do,i,68,71,1 a,i,i+1,i+13,i+12 *enddo
agen,4,17,21,1,,4,,12 *do,i,,91,1 a,i,i+1,i+13,i+12 *enddo
agen,18,1,35,1,,,3,60 a,1133,1136,1124,1121 !墙体 !632-635 *do,i,61,,1 a,i,i+1,i+61,i+60 *enddo !636-638 *do,i,,91,1 a,i,i+1,i+61,i+60 *enddo !639-2 *do,i,68,71,1 a,i,i+1,i+61,i+60 *enddo !3-6 *do,i,109,112,1 a,i,i+1,i+61,i+60 *enddo !7-650 *do,i,116,119,1 a,i,i+1,i+61,i+60
*enddo !651-6 *do,i,61,97,12 a,i,i+12,i+72,i+60 *enddo !655-656 *do,i,65,77,12 a,i,i+12,i+72,i+60 *enddo !657-658 *do,i,68,80,12 a,i,i+12,i+72,i+60 *enddo !659-662 *do,i,72,108,12 a,i,i+12,i+72,i+60 *enddo
agen,17,632,662,1,,,3,60 aatt,2,2,2 aesize,all,2 amesh,all asel,u,mat,,2 aplot !筒体
!*do,i,53,1124,60 !a,i,i+3,i+63,i+60 !*enddo
!*do,i,41,1112,60 !a,i,i+3,i+63,i+60 !*enddo
!*do,i,41,1112,60 !a,i,i+12,i+72,i+60 !*enddo
!*do,i,44,1115,60 !a,i,i+12,i+72,i+60 !*enddo
a,53,56,1136,1133 a,41,44,1124,1121 aatt,1,3,2 aesize,all,2 amesh,all asel,u,mat,,1 aplot
a,41,53,1133,1121
a,44,56,1136,1124 aatt,1,3,2 aesize,all,1.5 amesh,all nummrg,all numcmp,all allsel finish /solu
nsel,s,loc,z,0 d,all,all allsel antype,modal modopt,subsp,5 mxpand,5,,,yes solve finish
6网壳非线性分析
!先施加一个初始扰动，然后再用弧长法计算至完成 /prep7 et,1,beam188 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,gxy,1,8e10 !环杆与径杆的界面 sectype,1,beam,ctube,,0 secdata,0.09,0.0925,36 !斜杆截面
sectype,2,beam,ctube,,0 secdata,0.08,0.083,36 !建立模型 csys,2 fi=5 r=50 theta=15 !一半模型节点 *do,i,1,6 *do,j,1,24,1
n,24*(i-1)+j,r,(j-1)*theta,(i-1)*fi+60 *enddo *enddo
n,2000,100,0,90 type,1 mat,1
secnum,1 !环杆 *do,i,1,6 *do,j,1,23,1
e,24*(i-1)+j,24*(i-1)+j+1,2000 *enddo
e,i*24,i*24-23,2000 *enddo !径杆 *do,i,1,5 *do,j,1,24,1
e,24*(i-1)+j,24*(i-1)+j+24,2000 *enddo *enddo
n,1000,50,0,90 *do,i,121,144,1 e,i,1000,2000 *enddo !/pnum,node,1 eplot type,1 real,1 mat,1 secnum,2 !斜杆 *do,i,1,5 *do,j,1,23,1
e,24*(i-1)+j,24*(i-1)+j+25,2000 *enddo
e,i*24,i*24+1,2000 *enddo nummrg,all numcmp,all finish /solu csys,2
nsel,s,loc,z,60 nplot d,all,ux d,all,uy d,all,uz allsel nlgeom,on
f,145,fy,-1000 !施加初始扰动
allsel /OUT,SCRATCH, solve /OUT
nsel,u,loc,z,60 f,all,fz,-100000 outres,all,all arclen,on nsubst,10 allsel /OUT,SCRATCH, solve /OUT finish /post1 pldisp,1 /post26 nsol,2,145,u,z nsol,3,121,u,z nsol,4,97,u,z nsol,5,73,u,z
prod,6,1,,,load,,,12000 prod,7,2,,,,,,-1 prod,8,3,,,,,,-1 prod,9,4,,,,,,-1 prod,10,5,,,,,,-1 !prvar,2,3,4 /axlab,x,deflection /axlab,y,load /grid,1 xvar,7 plvar,6
6网壳特征值屈曲分析
!先施加一个初始扰动，然后再用弧长法计算至完成 /prep7 et,1,beam188 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,gxy,1,8e10 !环杆与径杆的界面 sectype,1,beam,ctube,,0 secdata,0.09,0.0925,36 !斜杆截面
sectype,2,beam,ctube,,0
secdata,0.08,0.083,36 !建立模型 csys,2 fi=5 r=50 theta=15 !一半模型节点 *do,i,1,6 *do,j,1,24,1
n,24*(i-1)+j,r,(j-1)*theta,(i-1)*fi+60 *enddo *enddo
n,2000,100,0,90 type,1 mat,1 secnum,1 !环杆 *do,i,1,6 *do,j,1,23,1
e,24*(i-1)+j,24*(i-1)+j+1,2000 *enddo
e,i*24,i*24-23,2000 *enddo !径杆 *do,i,1,5 *do,j,1,24,1
e,24*(i-1)+j,24*(i-1)+j+24,2000 *enddo *enddo
n,1000,50,0,90 *do,i,121,144,1 e,i,1000,2000 *enddo !/pnum,node,1 eplot type,1 real,1 mat,1 secnum,2 !斜杆 *do,i,1,5 *do,j,1,23,1
e,24*(i-1)+j,24*(i-1)+j+25,2000
*enddo
e,i*24,i*24+1,2000 *enddo nummrg,all numcmp,all finish /solu csys,2
nsel,s,loc,z,60 nplot d,all,ux d,all,uy d,all,uz allsel eplot
nsel,u,loc,z,60 nplot f,all,fz,-1 allsel solve finish /solu
antype,buckle bucopt,subsp,3 mxpand,3 solve finish /post1 set,list SET,FIRST PLDISP,0 SET,NEXT PLDISP,0 SET,NEXT PLDISP,0 finish 6正四角锥网架
!最大的索直径为10～15 cm 间 !支座竖向反力为1300 kN
!构竖向挠度为150 mm ,控制在l/ 400 = 170 mm 以内 !文件名为net D1=0.03
A1=3.1415926*D1*D1/4 D2=0.02
A2=3.1415926*D2*D2/4 ISTR=0.005 ro=7800 /filname,net /prep7 et,1,link8 et,2,link10 r,1,A1 r,2,A2,ISTR
SSTIF,ON mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,dens,1,ro !上层 n,1 n,2,2 n,3,4 n,4,6 n,5,8 n,6,10
ngen,21,6,1,6,1,,,-2 save *do,i,1,5,1 e,i,i+1 *enddo
egen,21,6,1,5,1 save
/view,1,1,1,1 eplot
*do,i,1,115,6 e,i,i+6 *enddo
egen,6,1,106,125,1 save !下层
n,1001,1,-2,-1 n,1002,3,-2,-1 n,1003,5,-2,-1 n,1004,7,-2,-1 n,1005,9,-2,-1
ngen,20,6,1001,1005,1,,,-2 *do,i,1001,1004,1 e,i,i+1 *enddo
egen,20,6,226,229,1 save
*do,i,1001,1109,6 e,i,i+6 *enddo
egen,5,1,306,324,1 save !底层部分 n,2001,3,-6,-10 n,2002,7,-6,-10
ngen,3,30,2001,2002,1,,,-10 save e,2001,2002 egen,3,30,401 e,2001,2031 e,2031,2061 e,2002,2032 e,2032,2062 eplot
!斜杆(上，下层间) *do,i,1001,1005,1 e,i,i-1000 e,i,i-999 e,i,i-994 e,i,i-993 *enddo
egen,20,6,408,427,1 !斜杆（下，底层间） e,2001,1019 e,2001,1021 e,2001,1037 e,2001,1039 e,2002,1021 e,2002,1023 e,2002,1039 e,2002,1041
egen,3,30,808,815,1 !索
n,3001,3,-6,-10 n,3002,7,-6,-10
ngen,3,30,3001,3002,1,,,-10 type,2 real,2 e,1002,3001
e,3001,3031 e,3031,3061 e,3061,1116 e,1004,3002 e,3002,3032 e,3032,3062 e,3062,1118
!索与底层弦梁节点Y向位移的约束 cp,1,uy,2001,3001 cp,2,uy,2002,3002 cp,3,uy,2031,3031 cp,4,uy,2032,3032 cp,5,uy,2061,3061 cp,6,uy,2062,3062 save finish /solu
nsel,s,loc,z,-1 nsel,a,loc,z,-39 nplot d,all,all allsel eplot
!前面的都没有问题 acel,,9.8
f,2031,fy,-500000 f,2032,fy,-500000 time,2 nsub,5 solve
!插入一步，即保存到net.lgw文件中 finish /post1 pldisp,1
etable,evol,volu ssum
*get,vtot,ssum,,item,evol wt=ro*vtot
*GET,uy,NODE,2032,U,Y *status,parm !前面的正确
! LGWRITE,net,lgw,G:\\yezi\\ANSYSC~1\\check\\6-1\\,COMMENT finish
/opt !进入优化器
OPANL,'net','lgw',' ' !指定分析文件 opvar,D1,DV,0.02,0.04 !定义设计变量 opvar,D2,DV,0.01,0.02 opvar,ISTR,DV,0.002,0.006
opvar,uy,SV,-0.17,-0.1 !定义状态变量 opsave,netvar,opt !存储数据 opvar,wt,obj,,,0.1, !定义目标函数 optype,subp !定义一阶方法 opsubp,45 !最大45次跌代 opexe !oplist,all !7桥墩
/TITLE,Buckling /config,nres,2500 /PREP7
ET,1,PLANE42,,,3,,2 ET,2,COMBIN14,,,2 R,1,4 R,2,1E6 MP,EX,1,25E6 MP,PRXY,1,0.17 MP,DENS,1,2500 K,1,-2 K,2,-2,100 K,3,2 K,4,2,100 A,1,2,4,3
LSEL,S,LOC,X,-1,1 LESIZE,ALL,,,4 ALLSEL
LSEL,S,LOC,Y,1,99 LESIZE,ALL,,,50 ALLSEL AMESH,ALL N,1000,-5,100 TYPE,2 REAL,2 E,1000,2
save，bucklingmodel /SOLU
NSEL,S,LOC,Y,100 NSEL,R,LOC,X,-1.5,1.5 F,ALL,FY,-25 ALLSEL
开始优化
列出所有设计序列 NSEL,S,LOC,Y,100 NSEL,R,LOC,X,-2 F,ALL,FY,-12.5 ALLSEL
NSEL,S,LOC,Y,100 NSEL,R,LOC,X,2 F,ALL,FY,-12.5 ALLSEL NSEL,S,LOC,Y D,ALL,ALL ALLSEL D,1000,ALL PSTRES,ON SOLVE finish
save,bucklingrst1 /solu
ANTYPE,BUCKLE BUCOPT,SUBSP,1 MXPAND,1 SOLVE /POST1 SET,,,1,,1 finish
save,bucklingrst2 /PREP7 TB,MISO,1 TBPT,,0.0002,5000 TBPT,,0.0004,9250 TBPT,,0.0006,13000 TBPT,,0.0008,16250 TBPT,,0.001,19000 TBPT,,0.0012,21250 TBPT,,0.0014,23000 TBPT,,0.0016,24250 TBPT,,0.0018,25000 TBPT,,0.0020,25250 TBPT,,0.0022,24312.5 TBPT,,0.0024,23375 TBPT,,0.0026,22437.5 TBPT,,0.0028,21500 TBPT,,0.003,20562.5 /SOLU fac=106
NSEL,S,LOC,Y,100 NSEL,R,LOC,X,-1.5,1.5 F,ALL,FY,-25*fac ALLSEL
NSEL,S,LOC,Y,100 NSEL,R,LOC,X,-2 F,ALL,FY,-12.5*fac ALLSEL
NSEL,S,LOC,Y,100 NSEL,R,LOC,X,2 F,ALL,FY,-12.5*fac ALLSEL /solu UPCOORD,0.1 SAVE,bucklingrst3 ANTYPE,STATIC NLGEOM,ON SSTIF,ON NEQIT,10 NSUBST,80 OUTRES,ALL,ALL ARCLEN,ON,2 SOLVE FINISH SAVE /POST1
SET,,,1,,,,875 PLDISP,1 PLNSOL,S,1,0,1 PLNSOL,U,X,0,1 PLNSOL,NL,SEPL,0,1 PLNSOL,EPEL,EQV,0,1 SAVE,bucklingrst4 FINISH 7隧道
!隧道模拟开挖命令流!希望能够得到加分! !土体的弹性摸量增加 /com fini /cle *set,x1,-12 *set,y1,-12 *set,w1,28.9 *set,h1,30.15
*set,x2,-25 *set,y2,-12 *set,w2,13 *set,h2,30.15 *set,x3,16.9 *set,y3,-12 *set,w3,13 *set,h3,30.15 *set,x4,-25 *set,y4,-30 *set,w4,.9 *set,h4,18 *set,th,0.4 /prep7
et,1,mesh200,2 et,2,mesh200,6 et,3,shell63 et,4,solid45 r,1,th
!壳的厚度
!壳的材料,C30混凝土
mp,ex,1,3.0e10 mp,prxy,1,0.2 mp,dens,1,2500 mp,ex,2,4.5e8 mp,prxy,2,0.32 !tb,dp,2
!tbdata,1,20,30, mp,dens,2,2700 mp,ex,3,4.51e8 mp,prxy,3,0.32 !tb,dp,3
!tbdata,1,20,30, mp,dens,3,2700 k,,0,0 k,,0,3.85 k,,0.88,5.5 k,,2.45,6.15 k,,4.02,5.5 k,,4.9,3.85 k,,4.9,0
larc,1,2,6,8.13 !定义两点之间的圆弧线，larc,p1,p2,pc,rad larc,2,3,6,3.21 larc,3,4,6,2.22 larc,4,5,2,2.22
!挖去岩石的材料 !保留岩石的材料 !用于3-D的2节点线 !用于3-D的4节点四边形 !厚度
*set,length_z,50 !洞的进深
larc,5,6,2,3.21 larc,6,7,2,8.13 larc,7,1,4,6 a,1,2,3,4,5,6,7
!产生面1
!blc4,x1,y1,w,h1 !block4,xcorner,ycorner,width,height,depth !blc4,x2,y2,w,h2 !block4,xcorner,ycorner,width,height,depth 产生面3 blc4,x1,y1,w1,h1 blc4,x2,y2,w2,h2 blc4,x3,y3,w3,h3 blc4,x4,y4,w4,h4 /pnum,area,1
aovl,1,2,3,4,5 !布尔操作重叠，得到面3 nummrg,all,,,,low numcmp,all l,1,8 l,7,9 l,6,10 l,2,11
lsel,s,line,,21,22,1 !用线分割面 lsel,a,line,,7 asbl,5,all
lsel,s,line,,21,24,3 lsel,a,line,,1 asbl,7,all
lsel,s,line,,22,23,1 lsel,a,line,,6 asbl,8,all nummrg,all,,,,low numcmp,all
lsel,s,line,,2,5,1 LCCAT,all lesize,all,,,3 lsel,s,line,,9,11,2 lsel,a,line,,6 lsel,a,line,,1 lesize,all,,,8 lsel,s,line,,8,10,2 lsel,a,line,,7 lesize,all,,,12 lsel,s,line,,21,24,1 lesize,all,,,10,2 type,2
asel,s,area,,5,8,1 amesh,all
!从四个角点上连接出四条线
!产生面2 !产生面3 !产生面4
!产生面5
!mopt,split,on 形单元
!mopt,split,err !mopt,split,warn !mopt,qmesh,main asel,s,area,,1 amesh,1
lsel,s,line,,12,13,1 lesize,all,,,8 lsel,s,line,,15,18,1 lesize,all,,,6,2 asel,s,area,,2,3,1 amesh,all lsel,s,line,,14 lesize,all,,,24 lsel,s,line,,19,20,1 lesize,all,,,6,2 lsel,s,line,,15,17,2 lsel,a,line,,8 LCCAT,all asel,s,area,,4 amesh,all LSEL,s,LCCA LDELE,all
nummrg,all,,,,low numcmp,all allsel
!指定网格划分选项，四边形单元如果违背了单元形状则将其分解为三角
!选用四边形网格划分器，当这种操作失败时，不调用备用的网格划分器
!以下开始拉伸成实体单元 !首先拉伸成壳单元 k,1000,,,-length_z l,1,1000 /view,1,1,1,1 /replot
EXTOPT,ESIZE,10,0, LSEL,S,LINE,,1,7,1 ADRAG,all,,,,,,25 gplot type,3 real,1 mat,1
ASEL,S,loc,z,-25 APLOT
lsel,s,loc,z,-25 lesize,all,,,10
MSHAPE,0,2D MSHKEY,1 amesh,all !拉伸岩石的实体 ASEL,invert aplot
EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1 TYPE,4 MAT,2
asel,r,area,,2,8,1 VDRAG,all,,,,,,25 allsel
!挖去部分岩石的实体 MAT,3
VDRAG,1,,,,,,25 EPLOT
nummrg,all,,,,low numcmp,all
!约束两侧面的X方向的约束 asel,s,loc,x,x2 asel,a,loc,x,x2+w4 da,all,ux,0 alls
!约束地面的Y方向的约束 asel,s,loc,y,y4 da,all,uy,0 alls
asel,s,loc,z,-length_z asel,a,loc,z,0 da,all,uz,0 allsel acel,,10 fini /solu
antype,static deltim,0.1,0.05,0.2 autots,on pred,on lnsrch,on nlgeom,on nropt,full
!使用自动时间步 !打开线性搜索 !打开大位移效果 !设定牛顿-拉普森选项
!选择梁单元，杀死
!打开时间步长预测器
cnvtol,f,,0.02,2,0.5 esel,s,type,,3
ekill,all esel,all esel,s,live nsle,s
nsel,invert d,all,all,0 nsel,all esel,all /PBC,ALL,,1 gplot solve
*do,ii,1,10,1 esel,s,mat,,3 nsle,s
nsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5) esln,r,1 ekill,all
esel,s,type,,3 nsle,s
nsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5) esln,r,1 ealive,all nsle,s ddele,all,all esel,all esel,s,live 元
nsle,s nsel,invert d,all,all,0 nsel,all esel,all solve *enddo /post1
/DEVICE,VECTOR,1 ESEL,S,type,,3 SET,1,LAST,1, PLNSOL,U,Y,0,1 SET,2,LAST,1, PLNSOL,U,Y,0,1 SET,6,LAST,1, PLNSOL,U,Y,0,1
!反向选择，将死单元上的节点约束所有自由度
!选择活单元，此时应该包含两部份，一是梁单元，二是未挖去的岩石单
!激活挖去的岩石对应的壳单元,并将其节点上的约束删除
!选择挖去的岩石，杀死
!选择活的单元 !选择活单元上的节点
!反向选择，即选择了死单元上的节点
!将死单元上的节点约束所有位移，使其不参与矩阵运算
SET,11,LAST,1, PLNSOL,U,Y,0,1 !*dim ,VUY_1 !*dim ,VUY_2 !*dim ,VUY_6 !*dim ,VUY_11 !SET,1,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_1,U,Y !PLETAB,EUY_1,AVG !PRETAB,EUY_1
!*vget,VUY_1,node,387,U,Y !SET,2,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_2,U,Y !PLETAB,EUY_2,AVG !PRETAB,EUY_2
!*vget,VUY_2,node,387,U,Y !SET,6,LAST,1, !PLNSOL,U,Y,0,1 !PLNSOL,U,SUM,0,1 !ETABLE,EUY_6,U,Y !PLETAB,EUY_6,AVG !PRETAB,EUY_6
!*vget,VUY_6,node,323,U,Y !SET,3,LAST,1, !PLNSOL,U,Y,0,1 !PLESOL,SMISC,4 !PLESOL,SMISC,5 !PLESOL,SMISC,6
!MX !MY !MXY
!ETABLE,p_i,SMISC, 13 !ETABLE,p_j,SMISC, 14 !ETABLE,p_k,SMISC, 15 !ETABLE,p_l,SMISC, 16 !ETABLE,p_i_3,SMISC, 17 !ETABLE,p_j_3,SMISC, 18 !ETABLE,p_k_4,SMISC, 20 !ETABLE,p_l_4,SMISC, 19 !ETABLE,p_k_5,SMISC, 21 !ETABLE,p_l_5,SMISC, 22 !fini !/post1
!etable,if,smisc,1 !etable,jf,smisc,7 !etable,im,smisc,6 !etable,jm,smisc,12
大家可以在这个的基础上进一步发掘 7.桩基节点单元 !定义相关参量
*DIM,MK,ARRAY,1,4,, , , !定义水平抗力系数比例常数m的矩阵 *SET,MK(1,1,1) , 6000 !第四层土的m值 *SET,MK(1,2,1) , 5000 !*SET,MK(1,3,1) , 4000 !*SET,MK(1,4,1) , 2000 !*DIM,FQ,ARRAY,1,4,, , , !*SET,FQ(1,1,1) , 40.53 !*SET,FQ(1,2,1) , -33.78 !*SET,FQ(1,3,1) , -16.14 *SET,FQ(1,4,1) , -4.8
*SET,QX, 100 !*SET,QZ, 100 *SET,b0, 1 *SET,h0, 0.5 !进入前处理器 /prep7
et,1,solid45 !et,2,combin14 !et,3,surf1 !keyopt,1,1,0
keyopt,2,3,0 !keyopt,3,4,1 !r,1
mp,ex,1,3.25e10 !mp,nuxy,1,0.3 mp,dens,1,2700 !建立整个桩模型的节点 n,1 n,6,1 fill,1,6
ngen,6,6,1,6,1,,,0.2 ngen,21,36,1,36,1,,0.5
/view,1,1,1,1 !/ang,1 /REP,FAST GPLOT
!建立桩的几何实体模型
第三层土的m值 第二层土的m值 第一层土的m值
定义桩周单位面积摩擦力 中性点以下的单位面积正摩擦力 中性点以上第三层的单位面积负摩擦力水平荷载 桩身单元 弹簧单元 表面效应单元 作为纵向三维弹簧 不包含中间节点 桩的材料属性 以正等侧显示 k,1 k,2,1 k,3,1,,1 k,4,,,1 k,5,,10, k,6,1,10, k,7,1,10,1 k,8,,10,1 /PNUM,kp,1 kplot
v,1,2,3,4,5,6,7,8 lsel,s,loc,y,0 lsel,a,loc,y,10 lesize,all,,,5 allsel
lsel,s,loc,y,5 lesize,all,,,20 vmesh,1
nummrg,all !压缩重复元素 numcmp,all !合并重复元素 !建立第一组弹簧模型,共7*16=102个弹簧 nsel,s,loc,x,1 nsel,r,loc,y,-0.1,8.1 nplot
ngen,2,1000,all,,,2 type,2 *do,i,1,17,1
*if,i,le,4,then !判断属于哪一层土 kk=1
*elseif,i,le,8,then kk=2
*elseif,i,le,12,then kk=3
*elseif,i,le,17,then kk=4 *endif *do,j,6,36,6
*get,vy,node,j+(i-1)*36,loc,y !获取节点的Y坐标值
r,kk+1,mk(1,kk,1)*abs(vy-8)*b0*h0/6 !根据Y值计算水平抗力系数，即弹簧的刚度 real,kk+1
e,j+(i-1)*36,1000+j+(i-1)*36 *enddo *enddo allsel
!最大单元号602
!第二组弹簧单元,共7*16=102个弹簧 allsel
nsel,s,loc,x,0 nsel,r,loc,y,-0.1,8.1 nsel,r,loc,z,-0.1,1.1 nplot
ngen,2,2000,all,,,-2 type,2 *do,i,1,17,1 *if,i,le,4,then kk=1
*elseif,i,le,8,then kk=2
*elseif,i,le,12,then kk=3
*elseif,i,le,17,then kk=4 *endif *do,j,1,31,6
*get,vy,node,j+(i-1)*36,loc,y
r,kk+10,mk(1,kk,1)*abs(vy-8)*b0*h0/6 real,kk+10
e,j+(i-1)*36,2000+j+(i-1)*36 *enddo *enddo allsel
!最大单元号704
!第三组弹簧单元,共7*16=102个弹簧 local,11,0,,,,,90 csys,0 allsel
nsel,s,loc,z,1 nsel,r,loc,y,-0.1,8.1 nsel,r,loc,x,-0.1,1.1 nplot
ngen,2,3000,all,,,,,2 type,2 *do,i,1,17,1 *if,i,le,4,then kk=1
*elseif,i,le,8,then kk=2
*elseif,i,le,12,then kk=3
*elseif,i,le,17,then kk=4 *endif *do,j,31,36,1
*get,vy,node,j+(i-1)*36,loc,y
r,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6 real,kk+20 esys,0
e,j+(i-1)*36,3000+j+(i-1)*36 *enddo *enddo allsel
!最大单元号码806
!第四组弹簧单元,共7*16=102个弹簧 allsel
nsel,s,loc,z,0 nsel,r,loc,y,-0.1,8.1 nsel,r,loc,x,-0.1,1.1 nplot
ngen,2,4000,all,,,,,-2 type,2 *do,i,1,17,1 *if,i,le,4,then kk=1
*elseif,i,le,8,then kk=2
*elseif,i,le,12,then kk=3
*elseif,i,le,17,then kk=4 *endif *do,j,1,6,1
*get,vy,node,j+(i-1)*36,loc,y
r,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6 real,kk+20 esys,0
e,j+(i-1)*36,4000+j+(i-1)*36 *enddo *enddo allsel
!最大单元号码为908 !底层弹簧
allsel csys,0
ngen,2,5000,1,36,1,,-1 type,2
r,1000,mk(1,1,1)*8*b0*h0/6*20 real,1000 *do,i,1,36,1 e,i,5000+i *enddo allsel gplot
!第一组表面效应单元 allsel r,1001 real,1001 local,13,0,,,,90 type,3 *do,i,1,16,1 *do,j,6,30,6
nsel,s,node,,j+(i-1)*36 nsel,a,node,,j+6+(i-1)*36 nsel,a,node,,j+42+(i-1)*36 nsel,a,node,,j+36+(i-1)*36 esys,13 esurf *enddo *enddo
!第二组表面效应单元 allsel
local,13,0,,,,90 type,3 *do,i,1,16,1 *do,j,1,29,6
nsel,s,node,,j+(i-1)*36 nsel,a,node,,j+6+(i-1)*36 nsel,a,node,,j+42+(i-1)*36 nsel,a,node,,j+36+(i-1)*36 esys,13 esurf *enddo *enddo
!第三组表面效应单元 allsel
local,14,0,,,,90
type,3 *do,i,1,16,1 *do,j,1,5,1
nsel,s,node,,j+(i-1)*36 nsel,a,node,,j+1+(i-1)*36 nsel,a,node,,j+37+(i-1)*36 nsel,a,node,,j+36+(i-1)*36 esys,14 esurf *enddo *enddo csys,0
!第四组表面效应单元 allsel
local,15,0,,,,90 type,3 *do,i,1,16,1 *do,j,31,35,1
nsel,s,node,,j+(i-1)*36 nsel,a,node,,j+1+(i-1)*36 nsel,a,node,,j+37+(i-1)*36 nsel,a,node,,j+36+(i-1)*36 esys,15 esurf *enddo *enddo save finish !进入求解器 /solu csys,0 allsel
nsel,s,loc,x,3 nsel,a,loc,x,-2 nsel,a,loc,z,3 nsel,a,loc,z,-2
d,all,all !约束边界条件 nsel,a,loc,y,-1 d,all,all nsel,a,loc,y,0 d,all,ux d,all,uz
ACEL,0,9.81,0, !施加重力加速度 !桩顶表面平均分布载荷
allsel
nsel,s,loc,y,10 nsel,r,loc,x,-0.1,1.1 nsel,r,loc,z,-0.1,1.1 nplot esln,s
sfe,all,6,pres,,5000 !桩顶表面分布荷载 nsel,r,loc,x,1
f,all,fx,QX !水平载荷 nsel,s,loc,y,10 nsel,r,loc,z,1
f,all,fz,QZ !allsel gplot
!第四层土的表面正摩擦力 csys,0 allsel
nsel,s,loc,y,0,1.5 esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,1,1) /PSF,PRES,TANx,2,0,1 gplot
!第三层土的表面负摩擦力 csys,0 allsel
nsel,s,loc,y,2.5,3.5 esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,2,1) /PSF,PRES,TANx,2,0,1 gplot
!第二层土的表面负摩擦力 csys,0 allsel
nsel,s,loc,y,4.5,5.5 esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,3,1) /PSF,PRES,TANx,2,0,1 gplot
!第一层土的表面负摩擦力 csys,0 allsel
水平载荷 nsel,s,loc,y,6.5,7.5 esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,4,1) /PSF,PRES,TANx,2,0,1 gplot allsel gplot time,1 autot,on nlgeom,on nropt,full csys,0 solve finish /POST1
!PLNSOL,U,Y,1,1 !PLNSOL,S,Z,1,1 PATH,path1,2,30,20, PPATH,P51X,1 PATH,STAT /PBC,PATH,1 !/REPLOT /PBC,PATH,0 AVPRIN,0,0, PDEF,SY,S,Y,AVG /PBC,PATH, ,0 PRANGE,1,1,10,YG, PLPATH,SY 7.桩接触分析 /clear /prep7 et,1,plane42 et,2,solid45
mp,ex,1,2.5e10 !mp,nuxy,1,0.2 !mp,dens,1,2500 !mp,ex,2,2.5e8 !mp,nuxy,2,0.4 !mp,dens,2,2000 !tb,dp,2
tbdata,1,19,32,30 !RECTNG,0,1,0,8, !RECTNG,1,5,0,6, !桩的弹性模量 桩的泊松比 桩的密度 土的弹性模量 土的泊松比 土的密度 粘聚力c为19,摩擦角为32度,膨胀角为30 面1 面2
RECTNG,0,5,0,-16, !面3 /pnum,area,1 /pnum,line,1 asel,s,,,2,3,1 aglue,all numcmp,all allsel aplot
lsel,s,,,1,3,2 !划分面1 lesize,all,,,2 lsel,s,,,2,4,2 lesize,all,,,16 amesh,1
lsel,s,,,6,8,2 !lesize,all,0.5 lsel,s,,,5 lesize,all,0.5 lsel,s,,,7 lesize,all,0.5 amesh,2
lsel,s,,,11 !lesize,all,0.5 lsel,s,,,12 lesize,all,0.5 lsel,s,,,10 lesize,all,0.5 lsel,s,,,9 lesize,all,0.5 amesh,3
EXTOPT,ESIZE,2,0, !EXTOPT,ACLEAR,1 type,2 mat,1
VEXT,1,,,,,-1
EXTOPT,ESIZE,2,0, !EXTOPT,ACLEAR,1 mat,2
VEXT,2,3,1,,,-1 allsel /view,1,1,1,1 eplot
k,1001,5,6,-5 k,1002,5,0,-5 k,1003,5,-16,-5
划分面2 划分面3 拉伸成体1 拉伸成体2，3 l,18,1001 l,1001,1002 l,1002,17 l,1002,1003 l,1003,21
al,22,36,37,38 !面4和面5 al,30,38,39,40
EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1
VEXT,19,20,1,-5 !vplot vsel,s,,,4,5,1 vplot
lsel,s,,,36,52,1 lesize,all,0.5 type,2 mat,2
vmesh,all !allsel
vsel,s,,,2,5,1 !VPLOT ESLV,S NSLV,S,1 EPLOT nummrg,all numcmp,all numcmp,all allsel /pnum,area,1 /view,1,1,1,1 /ang,1 /pnum,mat,1,2 eplot
!接触单元设置 allsel et,3,170 et,4,173 keyopt,4,9,0 keyopt,4,12,2 r,10 real,10 mp,mu,2,0.2 vsel,s,,,1 vplot
拉伸成体4和5 划分体4和5 合并体2，3，4，5上的重复单元及节点 asel,s,,,4,6,1 aplot
cm,target,area type,3 nsla,s,1 esurf,all esel,s,type,,3 eplot
/psyms,esys,1 allsel
vsel,s,,,2,5,1 vplot asel,s,,,13 asel,a,,,17 asel,a,,,22 aplot
cm,contact,area type,4 nsla,s,1 esurf,all esel,s,type,,4 eplot
/psyms,esys,1 allsel gplot finish /solu
asel,s,loc,x,0 aplot DA,all,symm asel,s,loc,x,5 DA,all,all aplot
asel,s,loc,z,0 da,all,symm asel,s,loc,z,-5 da,all,all asel,s,loc,y,-16 da,all,all allsel ACEL,0,9.8,0, nsel,s,loc,y,8 esln,s eplot
sfe,all,3,pres,,50 antype,static nlgeom,on time,1 nsub,100 outres,all,all allsel solv finish 8.反应谱分析
!注意：在进行反应普分析时，若想要得到单元的应力状态，必须将solid65单元的非线性属性去掉 /prep7
resume,dam_model,db keyopt,2,5,2 /solu
!施加边界条件 csys,0 dsys,0
nsel,s,loc,x,(2+0.9*0.75)*h nsel,a,loc,x,-1.5*h nplot d,all,ux allsel
nsel,s,loc,z,h
nsel,a,loc,z,-(1+1.5)*h nplot d,all,uz allsel
nsel,s,loc,y,-2*h nplot d,all,uy allsel gplot !模态分析 antype,modal modopt,subsp,10 solve
save,dam_dynamic_rst11,db finish !谱分析 /solu
antype,spectr spopt,sprs,10,yes
!输出solid65单元的节点应力
svtyp,2 sed,1,1,
freq,0.3444,0.3502,0.36,0.3739,0.4175,0.4613,0.5114,0.5121,0.5811,0.5837 sv,1.226,1.2079,1.1945,1.13,1.0312,0.9427,0.859,0.858,0.7657,0.7627 solve
save,dam_dynamic_rst2,db finish !模态扩展 /solu antype,modal expass,on
mxpand,10,,,yes,0.005 solve
save,dam_dynamic_rst3,db finish !合并模态 /solu
antype,spectr srss,0.15,disp solve
save,dam_dynamic_rst4,db finish !查看结果 /post1 /input,,mcom esel,s,type,,2 eplot SET,FIRST PLNSOL,U,SUM,1,1 set,next
PLNSOL,U,SUM,1,1 set,next
PLNSOL,U,SUM,1,1 SET,FIRST PLNSOL,S,1,0,1 PLNSOL,EPTO,1,0,1 SET,next PLNSOL,S,1,0,1 PLNSOL,EPTO,1,0,1 finish /post26 csys,0
nsel,s,loc,y,h nsel,r,loc,z,-0.75*h
nsel,r,loc,x,0 nplot
NSOL,2,1466,U,X,nux NSOL,3,1466,U,Y,nuy NSOL,4,1466,U,z,nuz XVAR,1 PLVAR,2,3,4
ESOL,2,1,19,S,1,s1_19 !变量2代表着：上游坝底侧面，1号单元上19号节点的第一主应力 ESOL,6,33,20,S,1,s1_20 应力
ESOL,7,761,9,S,1,s1_9 应力
ESOL,8,746,7,S,1,s1_7 应力
/axlab,x,time /axlab,y,stress1[pa] XVAR,1 PLVAR,5,6,7,8 8.静态分析 h=180
!大坝高度
!完成大坝浇筑的总的天数
!每天完成的大坝高度 !上游水位高度 !下游水位高度
!坝体不同混凝土材料的分界线，为100m !gm为水的容重 !计算坐标旋转角度 !辅助划分网格之用 !混凝土单元 !岩石单元
!输出solid65单元的节点应力
!100m以下混凝土的材料特性
fini_days=30 h1=100 h2=80
h3=100 gm=1000
!变量5代表着：下游坝顶侧面，746号单元上7号节点的第一主
!变量4代表着：上游坝顶侧面，761号单元上9号节点的第一主!变量3代表着：下游坝底侧面，33号单元上20号节点的第一主
ht=h/fini_days
fia1=90-atan(0.75)*180/3.1415926 !计算下游斜面夹角 fia2=90-fia1 /prep7 et,1,plane42 et,2,solid65 et,3,solid45 keyopt,2,5,2 !1号材料
mp,ex,1,2.85e10 MP,PRXY,1,0.167 TB,CONC,1,1,9,
TBDATA,,0.3,1,1.96e6,22e6 fs=22e6pa, mp,dens,1,2400 !2号材料
mp,ex,2,2.6e10 MP,PRXY,2,0.167 TB,CONC,2,1,9,
!100m以上混凝土的材料特性
!张开剪切传递系数为0.3，抗拉强度fc=1.96e6pa,抗压强度
TBDATA,,0.3,1,1.2e6,17.5e6 !张开剪切传递系数为0.3，抗拉强度fc=1.2e6pa,抗压强度fs=17.5e6pa, mp,dens,2,2400 !3号材料 mp,ex,3,2.9e10 mp,PRXY,3,0.3 mp,dens,3,2600 !创建几何模型 !创建关键点 k,1
k,2,0.9*0.75*h+0.1*h k,3,0.1*h,0.9*h k,4,0,0.9*h !创建线 l,1,2 l,2,3 l,3,4 l,4,1 !创建面 al,1,2,3,4
RECTNG,0,0.1*h,0.9*h,h RECTNG,-1.5*h,0,-2*h,0
!画坝顶矩形 !坝基左面矩形 !坝基正下方矩形
!总矩形
RECTNG,0,0.9*0.75*h+0.1*h,-2*h,0
RECTNG,0.9*0.75*h+0.1*h,(2+0.9*0.75)*h,-2*h,0 !坝基右面矩形 RECTNG,-1.5*h,(2+0.9*0.75)*h,-2*h,h AOVLAP,all /pnum,line,1 /pnum,area,1 nummrg,all numcmp,all aplot
!划分3号面的网格,坝体的顶部 lsel,s,,,3,5,2 lesize,all,,,5 lsel,s,,,12,13,1 lesize,all,,,2 amesh,3
!划分1号面的网格，坝体的下部 lsel,s,,,2,4,2 lesize,all,,,18 lsel,s,,,1 lesize,all,,,5 amesh,1 eplot
!合并重复元素
!将面单元进行粘贴布尔操作 !打开面，线的号码开关
!划分2号面的网格，坝体的正下方基岩 lsel,s,,,11 lesize,all,,,5 lsel,s,,,9,10,1 lesize,all,,,8,4 amesh,2
!划分4号面的网格 lsel,s,,,14 lesize,all,,,5,4 lsel,s,,,7 lesize,all,,,8,4 lsel,s,,,6
lesize,all,,,5,0.25 amesh,4 eplot
!划分5号面的网格 lsel,s,,,15,16,1 lesize,all,,,8,4 lsel,s,,,8
lesize,all,,,8,0.25 amesh,5 eplot
!连接4号线和13号线，并划分7号面的网格 lsel,s,,,4,13,9 LCCAT,all lsel,s,,,19 lesize,all,,,5,4 lsel,s,,,20,21,1 lesize,all,,,20 amesh,7 eplot
!连接2号线和12号线，并划分6号面的网格 lsel,s,,,2,12,10 lccat,all lsel,s,,,17 lesize,all,,,8,4 lsel,s,,,18,22,4 lesize,all,,,20 amesh,6 eplot
!删除前面连接的线元素 allsel lsel,r,lcca ldele,all
!拉伸成坝体单元，采用solid65单元和1号材料，随后改变h3以上坝体的材料属性 EXTOPT,ESIZE,8,0, type,2 mat,1
VEXT,1,3,2,,,-0.75*h /view,1,1,1,1 vplot eplot
!坝体底部单元，拉伸成岩土单元，采用solid45号单元和3号材料 EXTOPT,ESIZE,8,0, type,3 mat,3
VEXT,2,,,,,-0.75*h VEXT,4,5,1,,,-0.75*h eplot
!拉伸生成大坝两侧的岩石单元 allsel
EXTOPT,ESIZE,5,4, EXTOPT,ACLEAR,1 type,3 mat,3
VEXT,1,7,1,,,h /pnum,mat,1 eplot
!利用对称方法完成整个有限元模型 local,11,0,,,-0.75*h csys,11 DSYS,11
/PSYMB,CS,1
!定义局部坐标系 !激活局部坐标系
!显示不同的局部坐标系符号 !显示局部坐标系
!删除面单元
!将1号和3号面拉伸成一半的模型
!拉伸的份数
NSYM,Z,20000,all ENSYM,30000, ,20000,all nummrg,all numcmp,all eplot
!改变材料属性 esel,s,mat,,1 eplot nsle,s nplot
nsel,r,loc,y,h3+1,h+1 nplot esln,s
!选择h3以上的单元，以便转换材料属性
nsel,r,loc,z,-0.75*h+0.1,0.75*h-0.1
!选择单元上的节点
!选择材料编号为1的单元
!合并重复节点和单元
eplot MPCHG,2,all, /REPLOT allsel eplot
save,dam_model,db finish /solu
!施加边界条件 csys,0 dsys,0
nsel,s,loc,x,(2+0.9*0.75)*h nsel,a,loc,x,-1.5*h nplot d,all,ux allsel
nsel,s,loc,z,h
nsel,a,loc,z,-(1+1.5)*h nplot d,all,uz allsel
nsel,s,loc,y,-2*h nplot d,all,uy allsel gplot !施加重力 ACEL,0,9.8,0,
!选择上游坝面单元，施加水平静水压力 ESEL,S,TYPE,,2 eplot
nsel,s,loc,x,0
nsel,r,loc,z,-1.5*h+0.1,-0.1 nsel,r,loc,y,0.1,h1-5
esln,s nplot eplot
/PSF,PRES,NORM,2,0,1 sfgrad,pres,0,y,0,-gm sfe,all,2,pres,,gm*h1 /replot
!下游坝面单元，施加水平静水压力 allsel
!设定将显示压力的方向
!给定载荷梯度为-gm
!若采用(h1-0.1)将导致选择过多的单元,
!使得施加的面分布荷载出现负数的情况
!保存模型文件
!改变材料属性为2号材料
ESEL,S,TYPE,,2 nsle,s eplot
local,12,0,0.9*0.75*h+0.1*h,,,fia2 csys,12 /PSYMB,CS,1 DSYS,12
nsel,s,loc,y,0.1,h2/sin(0.75)-25 !若采用(h2/sin(0.75)-1)将导致选择过多的单元,
!使得施加的面分布荷载出现负数的情况
nsel,r,loc,z,-1.5*h+0.1,-0.1 nsel,u,loc,x,-1000,-2 esln,s nplot eplot
/PSF,PRES,NORM,2,0,1 sfgrad,pres,0,y,0,-gm sfe,all,4,pres,,gm*h2 /replot
!施加上游水位h1产生的渗透压力 !同时施加下游水位h2产生的浮托力 !选择第一部分坝底单元 csys,0 dsys,0 /psymb,cs,1 allsel
nsel,s,loc,y,0
nsel,r,loc,z,-135+0.1,0.1 esln,s
esel,r,type,,2 eplot
SFCUM,PRES,ADD
!设置载荷是叠加的
!下游水位产生的浮托力
sfe,all,5,pres,,gm*h2 sfgrad,pres,0,x,0,-p0 sfe,all,5,pres,,gm*h1 /replot
!选择第二部分坝底单元 allsel
nsel,s,loc,y,0
nsel,r,loc,z,-270+0.1,-135-0.1 esln,s
esel,r,type,,2 eplot p0=0
!设定将显示压力的方向
!给定载荷梯度为-gm
p0=gm*h1/(0.9*0.75*h+0.1*h)
sfgrad,pres,0,x,0,-p0 sfe,all,3,pres,,gm*h2 sfgrad,pres,0,x,0,-p0 sfe,all,3,pres,,gm*h1 SFCUM,PRES,repl
!检查坝底的载荷是否正确 nsel,s,loc,y,0 esln,s
esel,r,type,,2 /replot
save,dam_static_solu,db allsel
outres,all,all autots,on nsubst,20 solve
save,dam_static_rst,db /post1
esel,s,type,,2 set,1,last PLNSOL,U,X PLNSOL,U,Y PLNSOL,EPTO,1,0,1 PLNSOL,S,1,0,1 PLNSOL,S,3,0,1 /DEVICE,VECTOR,1 PLCRACK,0,0
!令斜率为0
!下游水位产生的浮托力
p0=gm*h1/(0.9*0.75*h+0.1*h)
!设置载荷是替代的
!获取第一载荷步，最后子步的结果
!画第一主应变 !画第一主应力 !画第三主应力
!显示开裂位置
!显示开裂的情况，包括裂缝的张开，闭合，压碎
!察看在各个积分点上的开裂情况 ETABLE,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,33,AVG PLETAB,44,AVG finish
因篇幅问题不能全部显示，请点此查看更多更全内容
Copyright © 2019- haog.cn 版权所有 赣ICP备2024042798号-2
违法及侵权请联系：TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务