要将所有的力转化到节点或者单元上,应该是你加载荷的时候,是加在几何上的 而不是直接加载节点上的
lzkhnu 发表于 2010-10-31 23:00
我认为我是施加到有限元模型上的,不是施加在几何模型上。
这是我的命令流,希望你帮我看看:
fini
/cle
/units,si
h=26.5
h1=30.45
h2=15.48
h3=11
gm=1000
/prep7
et,1,plane42
mp,ex,1,2.55e10
mp,nuxy,1,0.167
mp,dense,1,2400
keyopt,1,3,0
keyopt,1,5,0
keyopt,1,6,0
et,2,plane42
mp,ex,2,1.15e10
mp,nuxy,2,0.16
mp,dense,2,2300
keyopt,2,3,0
keyopt,2,5,0
keyopt,2,6,0
et,3,plane42
mp,ex,3,6e10
mp,nuxy,3,0.167
mp,dense,3,2660
keyopt,3,3,0
keyopt,3,5,0
keyopt,3,6,0
!定义关键点
!坝剖面
K,1,0,0
K,2,5,0
K,3,6,1
K,4,27.6,1
K,5,27.6,9.3
!材料分线
K,10,6,1.5
K,11,27.1,1.5
K,12,27.1,6.2
K,13,17.58198,6.2
K,14,3.38198,25
K,15,0.6,25
K,16,0.6,8
K,17,6,8
K,18,8.51815,18.2
K,19,0.6,18.2
K,20,0,25
K,21,27.6,6.2
!溢流面
!圆弧段
K,6,27.1,9.3
k,7,16.5011,10.16
!中间直线段
K,8,7.87419,21.66615
!顶部曲线段
K,9,0,25.98761
K,30,1.16549,26.5
K,31,5.16549,24.64157
!灌浆廊道
!圆心
K,22,2.5,3
K,23,2.5,5.25
K,24,2,5.25
K,25,5.5,5.25
K,26,5,5.25
K,27,5,3
K,28,3.75,5.25
k,29,22.1811,7.32
circle,28,1.25,,,180
circle,28,1.75,,,180
larc,6,7,29
larc,30,31,14,6.259
larc,9,30,15,1.32
nummrg,all
!地基
K,32,-0.6799,0.97128
K,37,-40.4299,0.97128
K,34,-40.4299,-52.02872
K,35,80.6,-52.02872
K,38,80.6,1
!定义面
a,15,14,18,19
a,19,18,13,12,11,10,17,16
a,9,20,15,14,18,13,12,21,5,6,7,8,31,30
a,20,15,19,16,17,10,11,12,21,4,3,2,1
a,22,27,26,25,36,24,23
asba,4,5
a,22,27,26,25,36,24,23
a,23,33,26,27,22
asba,4,5
a,37,34,35,38,4,3,2,1,32
wprota,,-90
wpoff,6,,8
asbw,2
wpoff,11.51815,,-1.8
asbw,8
wpoff,9.58185,,-4.7
asbw,6
wpoff,-21.1,,-0.5
asbw,10
wpoff,-6,,7
asbw,11
wpcsys,-1
esize,1
mshape,0,2d
mshkey,1
allsel
asel,s,,,5
amap,5,18,19,16,39
amesh,5
allsel
asel,s,,,2
amap,2,39,17,40,13
amesh,2
allsel
asel,s,,,9
amap,9,40,10,11,12
amesh,9
esize,1
mshape,0,2d
mshkey,1
asel,s,,,8
amesh,8
asel,s,,,12
amap,12,42,43,41,4
amesh,12
asel,s,,,6
amap,6,43,1,2,3
amesh,6
asel,s,,,10
amap,10,44,16,15,20
amesh,10
esize,1
mshape,0,2d
mshkey,0
asel,s,,,13
amesh,13
esize,2
mshape,0,2d
mshkey,1
asel,s,,,4
amap,4,37,34,35,38
amesh,4
esize,0.1
mshape,0,2d
mshkey,1
asel,s,,,7
amap,7,24,23,26,25
amesh,7
esize,0.5
mshape,0,2d
mshkey,1
asel,s,,,1
amesh,1
esize,1
mshape,0,2d
mshkey,0
asel,s,,,3
amesh,3
!加载及求解
allsel
FINISH
/SOL
SAVE
!施加约束
nsel,all
nsel,s,loc,x,-40.4299
nsel,r,loc,y,-52.02872,0.97128
d,all,ux
nsel,all
nsel,s,loc,x,80.6
nsel,r,loc,y,-52.02872,0.97128
d,all,ux
nsel,all
nsel,s,loc,x,-40.4299,80.6
nsel,r,loc,y,-52.02872
d,all,all
!施加重力
TIME,1
ALLSEL
acel,,9.8
SOLVE
TIME,2
allsel
sfl,41,pres,gm*9.8*h1
sftran
allsel
sfl,39,pres,gm*9.8*h2
sftran
nsel,all
nsel,s,loc,x,0
sfgrad,pres,,y,0,-gm*9.8
sf,all,pres,gm*9.8*(h1+1)
nsel,all
nsel,s,loc,x,27.6
nsel,r,loc,y,1,9.3
sfgrad,pres,,y,1,-gm*9.8
sf,all,pres,gm*9.8*h2
nsel,all
nsel,s,loc,y,9.3
sfgrad,pres,,x,27.6,1
sf,all,pres,gm*9.8*(h2-9.3)
nsel,all
nsel,s,loc,x,11.7611,16.5011
nsel,r,loc,y,10.16,16.48
sfgrad,pres,,y,16.48,-gm*9.8
sf,all,pres,
nsel,all
nsel,s,loc,x,16.5011,22.1811
nsel,r,loc,y,7.32,10.16
sfgrad,pres,,y,10.16,-gm*9.8
sf,all,pres,gm*9.8*(h2-10.16)
nsel,all
nsel,s,loc,x,22.1811,27.1
nsel,r,loc,y,7.32,9.3
sfgrad,pres,,y,9.3,-gm*9.8
sf,all,pres,gm*9.8*(h2-9.3)
nsel,all
nsel,s,loc,x,0
nsel,r,loc,y,0,11
sfgrad,pres,,y,11,-770*9.8
sf,all,pres,
esel,all
esel,s,,,411,416,1
sfgrad,pres,,x,0,-24777
sfe,all,3,pres,0,308201
esel,all
sfe,411,4,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1524.22*2
esel,all
sfe,405,4,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1524.22*3
esel,all
esel,s,,,381,404,1
sfgrad,pres,,x,0,-1524.22
sfe,all,2,pres,0,0.24*gm*9.8*(h1-h2)+9.8*gm*h2-1524.22*4
SOLVE |
|