/filname,ch08
/title,static analysis of compressor structuren with hole
/prep7
et,1,solid45 mp,ex,1,1.15e5 mp,prxy,1,0.30782 mp,dens,1,4.48e-9
k,1,226,,208.8 k,2,226,,258.7 k,3,157,,258.7 k,4,237.5,,220.3 k,5,229.2,,220.3 k,6,237.5,,208.8 k,7,126,,276.7 k,8,138,,276.7 k,9,102.5,,263 k,10,102.5,,248.7 k,11,237.5,,273.8 k,12,237.5,,264.1 k,13,135,,248.7 k,14,243.85,,273.8 k,15,243.85,,254.8 k,16,229.2,,254.8 k,17,162.5,,264.1
/pnum,kp,1 /pnum,line,1 /number,2 /view,1,,,-1 /auto,1
lstr,1,2 lstr,2,3 lstr,1,6 lstr,6,4 lstr,4,5 lstr,5,16 lstr,16,15 lstr,15,14 lstr,14,11 lstr,11,12 lstr,12,17
lstr,8,7 lstr,7,9 lstr,9,10 lstr,10,13 lplot
ltan,11,8 !创建切线 ltan,2,13 al,all
/pnum,line,0 /pnum,kp,0
k,18,,,, k,19,,,300
vrotat,all,,,,,,18,19,60 /view,1,1,1,1 /auto,1 save
wprot,30
cyl4,200,0,10,,,,300
!注意圆柱的深度方向是沿工作平面坐标系 !z向的； vsbv,1,2 save
/view,1,,,-1 /auto,1
wpcsys,-1,0 lstr,18,19 flst,2,2,8 fitem,2,0,0,0 fitem,2,220,0,0 circle,p51x,,,,60 flst,2,2,8 fitem,2,0,0,0 fitem,2,180,0,0 circle,p51x,,,,60
lsla,u !将与面关联的线从选择集中去除 /pnum,line,1 adrag,53,,,,,,52 adrag,54,,,,,,52 /pnum,line,0 allsel
aslv,u
vsba,all,all allsel wprot,,90 wprot,,,21 vsbw,all wprot,,,9 vsbw,all wprot,,,9 vsbw,all ldele,52,,,1 !删除轴线 csys,1
!激活坐标系转为总体柱坐标系 vsel,s,loc,x,180,220 cm,holevol,volu allsel vplot save
!定义周期对称分析选项
!定义高低角度组件
!即周期对称结构的对称边界
!组件命名时，前面的cyclic即为组件根名 asel,s,loc,y,0
cm,cyclic_m01l,area asel,s,loc,y, 60
cm,cyclic_m01h,area allsel
cyclic,6,60,1,’cyclic’
!此步为指定周期对称分析选项
!连接线和面以便于映射网格划分 esize,3
cmsel,s,holevol vsel,r,loc,y,21,30 aslv,s
wpcsys,-1,0 wprot,30 wpoff,200 cswpla,11,1
asel,u,loc,z,264.1 asel,u,loc,z,258.7
asel,u,loc,x,10 csys,1
asel,u,loc,y,30 accat,all lsla,s
lsel,r,loc,z,264.1 lccat,all lsla,s
lsel,r,loc,z,258.7 lccat,all
!生成网格，采用扫略网格划分 type,1
mshape,0,3d vsel,s,loc,y,0,21 vsweep,all
vsel,s,loc,y,21,30 vsel,u,loc,x,180,220 vsweep,all cmsel,s,holevol vsel,r,loc,y,21,30 vsweep,all
vsel,s,loc,y,0,30 asel,s,acca adele,all lsel,s,lcca ldele,all
csys,4
vsymm,y,all,,,,0,0 allsel
nummrg,all numcmp,all csys,1 nrotat,all save /solu
antype,static nsel,s,loc,x,237.5
nsel,r,loc,z,220.3,208.8 d,all,,,,,,ux,uy allsel
nsel,s,loc,z,208.8
d,all,,,,,,uy,uz allsel eplot save
omega,,,1191.11
nsel,s,loc,x,243.5
*get,no_nodes,node,,count f,all,fx,628232/no_nodes allsel save solve finish
/post1
plnsol,u,x,2,1 plnsol,u,y,2,1 plnsol,u,z,2,1 plnsol,u,sum,2,1 plnsol,s,x,0,1 plnsol,s,y,0,1 plnsol,s,z,0,1 plnsol,s,eqv,0,1 finish