渡槽流体耦合模态分析命令流,ANSYS

finish
/clear
/filename,ducaoa
!!!!!!!!!!!!渡槽参数!!!!!!!!!!!!!!!
h1=0.6 !上弦杆高
b1=0.3 !上弦杆宽
h2=0.5 !下弦杆高
b2=0.3 !下弦杆宽
h3=0.25 !顶横系梁高
b3=0.25 !顶横系梁宽
h4=0.5 !底横系梁高
b4=0.3 !底横系梁宽
dd=0.3 !竖杆尺寸
d1=0.12 !渡槽侧板厚
d2=0.12 !渡槽底板厚
hh=3.8 !渡槽高
bb=3.04 !渡槽宽
dp=(hh-0.25)*0.75 !渡槽水深

/prep7
et,1,BEAM188,1,,2 !梁单元BEAM188
ET,2,63 !两侧板单元shell181
et,3,fluid30 !流体fluid30流固耦合单元

mp,ex,1,32.5e9 !上弦杆C40混凝土
mp,prxy,1,0.3
mp,dens,1,2600

MP,EX,2,34.5e9 !下弦杆C50混凝土
MP,PRXY,2,0.3
MP,DENS,2,2590

MP,EX,3,30e9 !上横系梁C30混凝土
MP,PRXY,3,0.3
MP,DENS,3,2550

MP,EX,4,30e9 !下横系梁C30混凝土
MP,PRXY,4,0.3
MP,DENS,4,2565

MP,EX,5,30e9 !槽身侧板C30混凝土
MP,PRXY,5,0.3
MP,DENS,5,2200

MP,EX,6,30e9 !槽身底板C30混凝土
MP,PRXY,6,0.3
MP,DENS,6,2540

MP,EX,7,32.5e9 !C40混凝土竖杆
MP,PRXY,7,0.3
MP,DENS,7,2580

MP,EX,8,32.5e9 !C30混凝土竖杆
MP,PRXY,8,0.3
MP,DENS,8,2520

MP,EX,9,32.5e9 !C30墩体混凝土
MP,PRXY,9,0.3
MP,DENS,9,2560

MP,EX,10,32.5e9*1e4 !C30墩体混凝土
MP,PRXY,10,0.3
MP,DENS,10,2560

mp,dens,11,1000 !水的材料属性
mp,sonc,11,1500
mp,mu,11,0.5

!截面定义
sectype,1,beam,rect, !左侧拱
secdata,h1,b1
secoff,user,,b1/2
sectype,2,beam,rect !右侧拱
secdata,h1,b1
secoff,user,,-b1/2
sectype,3,beam,rect !左侧梁
secdata,h2,b2
secoff,user,,b1/2,
sectype,4,beam,rect !右侧梁
secdata,h2,b2
secoff,user,,-b1/2
sectype,5,beam,rect !左侧竖杆
secdata,dd,dd
secoff,user,,b1/2,
sectype,6,beam,rect !右侧竖杆
secdata,dd,dd
secoff,user,,-b1/2,
sectype,7,beam,rect !顶横隔梁
secdata,h3,b3
sectype,8,beam,rect !底横隔梁
secdata,h4,b4
secoff,user,h4/2+d1/2
sectype,9,beam,rect !墩内横隔梁
secdata,dbt,dbt
r,1,d1
!!!!!!!!!!!!!!渡槽模型建立!!!!!!!!!!!!!!!!!
CSYS,0
k, 1 ,-18.500,.000,.0000,
k, 2 ,-16.250,1.691,.0000,
l, 1 ,2
k, 3 ,-13.750,3.312,.0000,
l, 2 ,3
k, 4 ,-11.250,4.664,.0000,
l, 3 ,4
k, 5 ,-8.750,5.745,.0000,
l, 4 ,5
k, 6 ,-6.250,6.555,.0000,
l, 5 ,6
k, 7 ,-3.750,7.096,.0000,
l, 6 ,7
k, 8 ,-1.250,7.366,.0000,
l, 7 ,8
k, 9 ,1.250,7.366,.0000,
l, 8 ,9
k, 10 ,3.750,7.096,.0000,
l, 9 ,10
k, 11 ,6.250,6.555,.0000,
l, 10 ,11
k, 12 ,8.750,5.745,.0000,
l, 11 ,12
k, 13 ,11.250,4.664,.0000,
l, 12 ,

13
k, 14 ,13.750,3.312,.0000,
l, 13 ,14
k, 15 ,16.250,1.691,.0000,
l, 14 ,15
k, 16 ,18.500,.000,.0000,
l, 15 ,16
k, 17 ,-20,.000,.0000,
l, 17 ,1
k, 18 ,-16.250,.000,.0000,
l, 1 ,18
k, 19 ,-13.750,.000,.0000,
l, 18 ,19
k, 20 ,-11.250,.000,.0000,
l, 19 ,20
k, 21 ,-8.750,.000,.0000,
l, 20 ,21
k, 22 ,-6.250,.000,.0000,
l, 21 ,22
k, 23 ,-3.750,.000,.0000,
l, 22 ,23
k, 24 ,-1.250,.000,.0000,
l, 23 ,24
k, 25 ,1.250,.000,.0000,
l, 24 ,25
k, 26 ,3.750,.000,.0000,
l, 25 ,26
k, 27 ,6.250,.000,.0000,
l, 26 ,27
k, 28 ,8.750,.000,.0000,
l, 27 ,28
k, 29 ,11.250,.000,.0000,
l, 28 ,29
k, 30 ,13.750,.000,.0000,
l, 29 ,30
k, 31 ,16.250,.000,.0000,
l, 30 ,31
l, 31 ,16
k, 32 ,20,.000,.000,
l, 16 ,32
l, 2 ,18
l, 3 ,19
l, 4 ,20
l, 5 ,21
l, 6 ,22
l, 7 ,23
l, 8 ,24
l, 9 ,25
l, 10 ,26
l, 11 ,27
l, 12 ,28
l, 13 ,29
l, 14 ,30
l, 15 ,31

k,33,-20,hh,
k,34,-18.50,hh
k,35,-16.25,hh
k,36,-13.75,hh
k,37,20,hh,
k,38,18.50,hh
k,39,16.25,hh
k,40,13.75,hh

l,17,33
l,1,34
l,2,35,
l,3,36
l,16,38
l,15,39
l,14,40
l,14,40
l,32,37
allsel

lgen,2,all,,,,,-bb-2*b1,100,0
lsel,s,,,1,15
lsel,a,,,1+54,15+54
cm,l1,line !上弦杆
lsel,s,,,16,32
lsel,a,,,16+54,32+54
cm,l2,line !下弦杆
lsel,s,,,33,46
lsel,a,,,33+54,46+54
cm,l3a,line !竖杆c40混凝土
lsel,s,,,47,54
lsel,a,,,47+54,54+54
cm,l3b,line !竖杆c30混凝土

lsel,none
*do,i,16,32,1
l,i,i+100
*enddo
l,1,101
wpcsys,1
wpoff,,,-b1-d1/2
lsbw,all
wpoff,,,-bb+d1
lsbw,all
cm,l4,line !底横隔梁
allsel

wpcsys,1
wprota,,-90
wpoff,,,dp !渡槽内水深度
cmsel,s,l3a
lsbw,all
wpoff,,,hh-dp
lsbw,all
cm,l3a,line !竖杆C40
wpoff,,,-hh+dp
cmsel,s,l3b
lsbw,all
cm,l3b,line !竖杆C30

wpcsys,1
ksel,s,,,43,61,18
ksel,a,,,58,76,18
kgen,2,all,,,,hh,,200,,0

wpcsys,1
allsel
a,243,43,58,258
a,43,58,76,61
a,61,76,276,261

wprota,,,90
wpoff,,,-18.50
asbw,all,
wpoff,,,2.25
asbw,all,
*do,i,1,13
wpoff,,,2.5
asbw,all
*enddo
wpoff,,,2.25
asbw,all,
wpcsys,1
wprota,,-90
wpoff,,,dp !渡槽内水深度
asbw,all,
cm,aa,area

asel,s,loc,z,-b1-d1/2
cm,aa1,area,
asel,s,loc,y,0
cm,aa2,area
asel,s,loc,z,-bb+d1/2-b1
cm,aa3,area

wpcsys,1
lsel,none
*do,i,33,40,1
l,i,i+100
*enddo
*do,i,141,150
l,i,i+10
*enddo
wpoff,,,-b1-d1/2
lsbw,all
wpoff,,,-bb+d1
lsbw,all
cm,l5,line !顶横隔梁

cmsel,s,aa2
vdrag,all,,,,,,402
allsel
nummrg,kp
numcmp,all

wpcsys,1
ss=1 !网格比例大小
/eshape,1
!拱结构网格划分
cmsel,s,l1
lesize,all,,,5*ss
lsel,r,loc,z,0
latt,1,,1,,,,1,
lmesh,all
cmsel,s,l1
lsel,r,loc,z,-bb-2*b1
latt,1,,1,,,,2,
lmesh,all

!梁网格划分
cmsel,s,l2
lsel,u,,,16,32,16
lsel,u,,,65,81,16
lesize,all,,,5

*ss
lsel,s,,,16,32,16
lsel,a,,,65,81,16
lesize,all,,,3
cmsel,s,l2
lsel,r,loc,z,0
latt,2,,1,,,,3,
lmesh,all
cmsel,s,l2
lsel,r,loc,z,-bb-2*b1
latt,2,,1,,,,4,
lmesh,all

!顶横隔梁网格划分
cmsel,s,l5
lsel,u,length,,b1+d1/2
lesize,all,,,8*ss
lsel,s,length,,b1+d1/2
lesize,all,,,1*ss
cmsel,s,l5
latt,3,,1,,,,7,
lmesh,all

!底横隔梁网格划分
cmsel,s,l4
lsel,u,length,,b1+d1/2
lesize,all,,,8*ss
lsel,s,length,,b1+d1/2
lesize,all,,,1*ss
cmsel,s,l4
latt,4,,1,,,,8,
lmesh,all

!吊杆网格划分
allsel
cmsel,s,l3a
lesize,all,0.5*ss
lsel,r,loc,z,0
latt,7,,1,,,,5
lmesh,all !C40混凝土竖杆

cmsel,s,l3a
lsel,r,loc,z,-bb-2*b1
latt,7,,1,,,,6
lmesh,all

cmsel,s,l3b
lesize,all,0.5*ss
lsel,r,loc,z,-bb-2*b1
latt,8,,1,,,,6
lmesh,all !C30混凝土竖杆

cmsel,s,l3b
lesize,all,0.5*ss
lsel,r,loc,z,0
latt,8,,1,,,,5
lmesh,all

cmsel,s,aa1
lsla,s,
lesize,all,0.5
aatt,5,1,2,,
amesh,all

cmsel,s,aa2
lsla,s,
lesize,all,0.5
aatt,5,1,2,,
amesh,all

cmsel,s,aa3
lsla,s,
lesize,all,0.5
aatt,5,1,2,,
amesh,all

lsel,s,length,,dp
lesize,all,0.5
mshkey,1
mshape,3d
vatt,11,,3,
vmesh,all

lsel,s,loc,x,-20
lsel,a,loc,x,20
lclear,all
allsel
numcmp,all

!施加约束
allsel
asel,s,loc,y,dp
nsla,s,1
d,all,pres,0
allsel
esel,s,type,,2,3
vsel,s,loc,y,0,dp
asel,s,ext
asel,u,loc,x,20
asel,u,loc,x,-20
asel,u,loc,y,dp
sfa,all,,fsi

nsel,s,loc,x,18.5,20
nsel,a,loc,x,-20,-18.5
nsel,r,loc,y,0
nsel,r,loc,z,0
d,all,ux,,,,,uy,uz
nsel,s,loc,x,18.5,20
nsel,a,loc,x,-20,-18.5
nsel,r,loc,y,0
nsel,r,loc,z,-bb-2*b1
d,all,ux,,,,,uy,uz
asel,s,loc,x,-20
asel,a,loc,x,20
lsel,s,ext
nsll,s,1
nsel,u,loc,z,-b1-d1/2
nsel,u,loc,z,-bb-b1+d1/2
nsel,u,loc,y,0
cm,nnn,node
nsel,s,loc,x,-20
nsel,a,loc,x,20
asel,s,loc,x,-20
asel,a,loc,x,20
nsla,u,
cmsel,u,nnn
dsym,symm,x
allsel
acel,,9.8
/pnum,mat,1
/num,1
/eshape,1
eplot

/solu
antype,modal
upcoord,1,on
prstes,on
modopt,unsym,10 !非对称模态提取方法处理流固耦合问题
eqslv,front
mxpand,10,,,1
solve

/POST1
SET,FIRST
PLNSOL,U,SUM,0,1.0