请教，AWE极限载荷求法？

asj78

压力容器，内压。能说下大致的步骤吗，详细更好，谢谢。

pfsr

施加一个极大荷载，用很小的荷载步，做非线性分析。最后一个收敛的荷载值就是极限荷载。

asj78

谢Pfsr版主了，有无具体的例子或文章，呵呵，看懂了，但手麻！

pfsr

版内有类似的例子

koalabear

我也正在学习极限分析，但是始终没学会，在网上流传的接管补强的例题我跑了一上午都没有算完。

koalabear

! ***************环境设置************************
/clear
/filn, E54  
/title, FEA of connecting zone of nozzle to cylinder
/units,si                        !采用国际单位制
! ********* 参数设定*********
Rci=1000                         ! 筒体内半径
tc=30                            ! 筒体厚度
Rco=Rci+tc                       ! 筒体外半径
Lc=4000                          ! 筒体长度
Rno=530                          ! 接管外半径
tn=15                            ! 接管厚度
Rni=Rno-tn                       ! 接管内半径
Li=193                           ! 接管内伸长度
Ln=500                           ! 接管外伸长度
rr1=30                           ! 焊缝外侧过渡圆角半径
rr2=15                           ! 焊缝内侧过渡圆角半径
pi=6                             ! 内压
pc=pi*Rci**2/(Rco**2-Rci**2)     ! 筒体端部轴向平衡面载荷
!****************前处理***************************
/prep7
et,1,95                          ! 定义单元类型   
et,2,45                          ! 定义单元类型         
mp,ex,1,2e5                      ! 定义材料的弹性模量
mp,nuxy,1,0.3                    ! 定义材料的泊松比
TB,BISO,1,1,2,  
TBTEMP,0
TBDATA,,235,,,,,
!****************建立模型***************************
cylind,Rco,Rci,0,-Lc/2,90,270,   ! 生成筒体
wpoff,0,0,-Lc/2                  ! 将工作面沿-Z向移动Lc/2
wprot,0,90,                      ! 将工作面沿yz旋转90度
cylind,Rno,Rni,-Ln-Rci-tc,-Rci+Li,90,180,  ! 生成接管
vovlap,all                       ! 体overlap布尔运算
vsel,s,,,7                       ! 选择筒体
*afun,deg                        ! 设定角度函数中单位为角度
ang1=2*nint(asin(Rno/Rci))       ! 计算接管区切割角度
wprot,0,0,-90+ang1               ! 旋转坐标系
vsbw,all                         ! 切割筒体
afillt,21,12,rr1                 ! 筒体与接管外表面圆角  
afillt,23,35,rr2                 ! 筒体内表面与接管外表面圆角
afillt,14,25,rr2                 ! 生成下辅助过渡圆角
afillt,13,19,rr1                 ! 生成上辅助过渡圆角
alls
askin,91,64                      ! 根据接管外过渡圆角在接管内外表面上的交线蒙面
vsba,4,13                        ! 切割外伸接管
askin,83,72                      ! 根据接管内过渡圆角在接管内外表面上的交线蒙面
vsba,5,18                        ! 切割内伸接管
askin,67,75                      ! 根据内外圆角边界蒙皮生成切割面24
vsba,2,24                        ! 切割筒体
wprot,0,,90-ang1                 ! 旋转坐标系
wprot,0,90,0                     ! 旋转坐标系
wpoff,0,0,-2*Rni                 ! 移动坐标系至接管区轴向切割位置
vsbw,all                         ! 切割筒体
VDELE,       3, , ,1             ! 删除筒体上开孔失去的部分
ADELE,      25, , ,1             ! 删除辅助面
ADELE,      35, , ,1             ! 删除辅助面
al,33,65,5                       ! 生成外圆角区域所需的端面
al,49,73,40                      ! 生成外圆角区域所需的端面
al,66,7,37                       ! 生成内圆角区域所需的端面
al,41,74,50                      ! 生成内圆角区域所需的端面
asel,s,,,31                      ! 选择外圆角区域各面
asel,a,,,38,39
asel,a,,,4,6,2
va,all                           ! 生成外圆角区域体
asel,s,,,21                      ! 选择内圆角区域各面
asel,a,,,41,42
asel,a,,,5,9,4
va,all                           ! 生成内圆角区域体
vsel,s,,,1,5,2                   ! 选择过渡区的体素
vsel,a,,,8,9
vsel,a,,,6
vadd,all                         ! 合并体素
alls
aadd,4,43,29,26,5,51             ! 合并过渡区端面面素
aadd,50,9,52,28,45,6             ! 合并过渡区端面面素
LESIZE,67, , ,80, , , , ,1       ! 指定圆周方向的剖分数
LESIZE,5, , ,5, , , , ,1         ! 指定过渡圆弧的剖分数
vsweep,10                        ! 扫略剖分该区域
alls
type,2
MSHAPE,0,3D                      ! 设定三维映射剖分
MSHKEY,1
VMESH,7                          ! 映射剖分接管上段
VMESH,4                          ! 映射剖分接管上段
accat,56,65                      ! 筒体开孔相邻区区域规则化连面
LESIZE,108, , ,40,0.5, , , ,1    ! 指定筒体接管相邻区的剖分数
VMESH,13                         ! 扫略剖分该区域
vsweep,11                        ! 扫略剖分筒体其它部分
vsweep,12
vsweep,2
alls
nummrg,all                        ! 合并所有相同项
numcmp,all                        ! 压缩
fini                             
!****************求解***************************
/solu
asel,s,loc,z,0                    ! 选择筒体端面
SFA,all,1,PRES,-pc                ! 施加端面平衡面载荷
asel,s,loc,x,0                    ! 选择对称面
asel,a,loc,z,-Lc/2                ! 选择对称面
DA,all,SYMM                       ! 施加对称约束
asel,s,loc,y,Ln+Rci+Tc            ! 选择接管端面
DA,all,UY,                        ! 约束轴向位移
asel,s,,,38                       ! 选择内表面面素
asel,a,,,34
asel,a,,,9
asel,a,,,41
asel,a,,,12
asel,a,,,10
asel,a,,,7
asel,a,,,22
asel,a,,,25
asel,a,,,15
asel,a,,,27
asel,a,,,23
cm,ainner,area
outres,all,all
pred,on
KBC,0
autots,on
NLGEOM, on
nsubst,1000
time,1                                 ! 载荷步1
SFA,all,1,PRES,pi                      ! 施加内压
alls
solve
fini
!****************后处理***************************
/post1
PLNSOL,S,INT,0,1                  ! 显示应力云图
fini