波形钢腹板箱梁命令流

!简支SPC40m-----------工况1正载

finish$/clear$/prep7 !*************************************************** ***************************进入前处理

/title,the analysis of SPC

/nerr,0,90000000 !修改限制的警告数量

!---------------------------------------------------------------------参数化

*SET,T1,0.012

et,1,solid45

et,2,shell63

et,3,link10 !

mp,alpx,1,1e-5

mp,ex,2,2.1e11

mp,dens,2,7800

mp,nuxy,2,0.3

mp,alpx,2,1.2e-5

mp,ex,3,1.95e11

mp,dens,3,7800

mp,prxy,3,0.3

mp,alpx,3,1.2e-5

r,1,T1

!---------------------------------------------------------------------预应力钢束

*SET,areagjx,1.4e-4 !单根面积

*dim,gjarea,array,4 !定义钢束

*SET,gjarea(1),12*areagjx

*SET,gjarea(2),12*areagjx

*SET,F,1860*140*0.6 !单根张拉力

*SET,E,1.95e11

*SET,L,1

*SET,detL,F*L/(E*areagjx) !每根钢绞线单元的伸长量

*SET,intrn,detL/(detL+L) !每根钢绞线的拉应变

r,4,gjarea(1),istrn !定义体外束M15-19实常数

!---------------------------------------------------------------------几何模型

local,11,0,-3.0065 !定义局部坐标系11

k,1,3.0065

k,2,2.6185

k,3,2.0185

k,4,1.5485

k,5,1.2425

k,6

k,7,,-0.36

k,8,1.2425,-0.36

k,9,1.5485,-0.471

k,10,2.0185,-0.642

k,11,2.6185,-0.36

k,12,3.0065,-0.36

k,13,1.55,-1.93

k,14,1.08,-1.759

k,15,0.84,-2

k,16,,-2

k,17,,-2.28

k,18,0.84,-2.28

k,19,1.08,-2.28

k,20,1.55,-2.28

k,21,-2.0185

k,22,-3.3685

k,23,-3.3685,-0.28

k,24,-2.0185,-0.642

*do,j,0,4

v,1+i-100+j,2+i-100+j,11+i-100-j,12+i-100-j,1+i+j,2+i+j,11+i-j,12+i-j

*enddo

*enddo

*do,i,100,5800,100

*do,j,0,2

v,13+i-100+j,14+i-100+j,19+i-100-j,20+i-100-j,13+i+j,14+i+j,19+i-j,20+i-j v,21+i-100,22+i-100,23+i-100,24+i-100,21+i,22+i,23+i,24+i

*enddo

*enddo

!---------------------------------------钢腹板建模

vsel,all

wpoffs,,,0.35

*do,i,1,58

vsel,s,loc,z,(i-1)*0.63,i*0.63

vsbw,all

wpoffs,,,0.63

*enddo

!-------------------复制1/4直波

vsbw,all

!-------------------切割底板结合部--直段内侧

wpstyl,defa$wprota,,,90$wpoffs,,,-1.790245

vsel,s,loc,x,-1.45,-1.92$vsel,r,loc,y,-1.5,-2.2

*do,i,0,29

vsel,u,loc,z,-0.28+i*1.26,0.63+i*1.26

*enddo

vsbw,all

!-------------------切割底板结合部--斜段1

kwplan,,2456,2294,2455

*do,i,0,28

vsel,s,loc,x,-1.45,-1.92$vsel,r,loc,y,-1.5,-2.2$vsel,r,loc,z,0.35+i*1.26,0.7+i*1.26 vsbw,all$wpave,-1.790245,-2.2,0.63+(i+1)*1.26

*enddo

!-------------------切割底板结合部--斜段2

kwplan,,2296,2295,2458

*do,i,0,28

vsel,s,loc,x,-1.45,-1.92$vsel,r,loc,y,-1.5,-2.2$vsel,r,loc,z,0.98+i*1.26,1.26+i*1.26 vsbw,all$wpave,-1.592755,-2.2,1.26+(i+1)*1.26

*enddo

!-------------------------------------------------顶板结合部切割

!-------------------切割顶板结合部--直段外侧

wpstyl,defa$wprota,,,90$wpoffs,,,-1.124255

vsel,s,loc,x,-0.98,-1.45$vsel,r,loc,y,-1.5,0

*do,i,0,28

vsel,u,loc,z,0.35+i*1.26,1.26+i*1.26

*enddo

vsbw,all

!-------------------切割顶板结合部--直段内侧

wpstyl,defa$wprota,,,90$wpoffs,,,-1.321745

vsel,s,loc,x,-0.98,-1.45$vsel,r,loc,y,-1.5,0

*do,i,0,29

vsel,u,loc,z,-0.28+i*1.26,0.63+i*1.26

*enddo

vsbw,all

!-------------------切割顶板结合部--斜段1

kwplan,,2765,2763,2766

*do,i,0,28

vsel,s,loc,x,-0.98,-1.45$vsel,r,loc,y,-1.5,0$vsel,r,loc,z,0.35+i*1.26,0.7+i*1.26 vsbw,all$wpave,-1.321745,0,0.63+(i+1)*1.26

*enddo

!-------------------切割顶板结合部--斜段2

kwplan,,2597,2767,2598

*do,i,0,28

vsel,s,loc,x,-0.98,-1.4$vsel,r,loc,y,-1.5,0$vsel,r,loc,z,0.98+i*1.26,1.26+i*1.26 vsbw,all$wpave,-1.124255,-2.2,1.26+(i+1)*1.26

*enddo

!-------------------建立支点处几何体

a,10,2762,2291,13$a,2762,9,14,2291

a,9,8,15,14$a,8,7,16,15

asel,s,loc,z,0

vext,all,,,,,-1.555

vsel,s,loc,z,0,-1.555

vgen,2,all,,,,,38.445

!-------------------建立横隔板几何体