断裂分析实例(J积分和应力强度因子).docx
《断裂分析实例(J积分和应力强度因子).docx》由会员分享，可在线阅读，更多相关《断裂分析实例(J积分和应力强度因子).docx（9页珍藏版）》请在冰豆网上搜索。
分享一个断裂分析实例（J积分和应力强度因子）（最终版）
有些地方没有加注释。
加起来比较麻烦。
后面不再修改了。
建议参考help关于fracture部分以及vm143部分。
有更新[attach]92985[/attach]以上附件没有删除。
一、问题描述
如图所示为一等厚度空心圆盘（参看图1.1），厚度4mm，内径r=10+Δ，Δ=22mm，
外径R=500mm，材料属性数据：
双线性（参看图1.3）E=2.1×1011Pa，ET=6.0×109Pa，屈服极限σs=500MPa，μ=0.3，密度ρ=8500kg/m3，采用Mises屈服准则，裂纹初始长度为a0=2.0mm。
裂纹如图1.2，裂纹位于0度90度180度和270度的位置，图中粗短线表示裂纹。
要求用有限元解
（1）载荷为均布拉力q=150MPa时的KI。
（参看图1.1）
（2）载荷为均布拉力q=1200MPa时的J积分值。
（参看图1.1）
（3）当载荷为转速n=600r/min的KI，并计算沿裂纹尖端不同路径的积分值，与
KI比较。
[attach]93106[/attach][attach]93107[/attach][attach]93108[/attach]
二、求解
求解中统一采用国际单位制，长度m，压力、应力与弹性模量Pa，密度Kg/m3，转速rad/s。
对圆盘的1/4进行ANSYS建模，网格划分如图2.1。
单元类型为6节点三角形单元plane2。
裂纹附近单元边长为0.0002m。
载荷施加如图2.2，扇形两条半径（裂纹处除外）上施加对称位移边界条件，弧上加均布拉力。
裂纹处无位移约束。
[attach]93109[/attach][attach]93110[/attach]
计算的应力强度因子和J积分结果如表2.1。
前面给出的J积分由于坐标系错误，做
法不对，现在改正1200Mpa下的J积分结果。
但是，结果差别不大。
因为，即使在局部坐标
系下，J积分中用到的XG,YG,ZG的坐标还是在全球坐标系下的。
坐标系的改变只对J积分的第一项有一点影响（对Y坐标积分这部分，Y坐标变了）。
新加两个附件：
改正后求解1200MPa下J积分的命令流，及两张路径定义的图片命
令流里面涉及到路径定义的命令是GUI方式进行的。
所以要分开看。
[attach]93704[/attach]。
表2.1应力强度因子和J积分结果
载荷项目不同路径下的结果
150MPaKI3.1359e72.986e7
3.3068e72.7695e7
（此处J积分没改）J-Integral3193.416293182.079173245.22042
3167.23352
1200MPaKI5.4207e95.3398e9
5.382e95.285e9
J-Integral4701408.684777189.414674756.71
4709657.93
600r/min
KI
55011
54901
5467254463
（此处J积分没改）J-Integral1.56498021.57193831.53405748
1.57113766
三、一些需要讨论的问题（后面补的帖解决了大部分问题，这些留在这里供参考）：
1、求解方法选择
1200MPa下采用默认的牛顿-拉普生法老是遇到收敛问题，经常不收敛，或者单元高
度扭曲。
而且需要采用的载荷子步也要很多。
（有人做的时候很少遇到这种情况，这与网格划分质量及网格多少有关，网格越少越容易收敛）
采用弧长法能解决这个问题。
内径参数r=10+Δ，Δ=22mm这种情况直接就求出来了。
Δ参数不同时可能也不能直接就求出来，也不收敛。
但是可以在450，550MPa的地方加两个载荷步就能收敛。
如Δ=6mm时，直接用一个载荷步不收敛，我分成150，450，550，1200MPa
几个载荷步（依次存为载荷步1，2，3，4）y，再求解150到1200MPa就成功了（lssolve,1,4,1）。
不收敛的原因在于500MPa是屈服极限。
后来验证过，网格画稀一点，不用弧长法直接用默认设置就可以收敛。
2、J积分解法的疑问
我是按照help上的介绍做的，但是对其求解的做法有些怀疑。
个人认为求解
δuy/δy偏导数的话，应该是把路径沿y方向移动才是对的。
Help里面δux/δx的求法在数学上正确的。
另外，ANSYS里面给了一个路径项求导的操作：
generalpostproc->path
operation->differentiate。
（differentiate是不是求导，请指教）那这个东西用来求
偏导数不行吗？
为什么help里面要那么来求偏导数。
3、这个实例的建模
这个实例的建模，我是建的1/4模型。
1/8模型也可以建出来，但是我对于1/8模型
还能不能用对称边界条件有怀疑。
1/4模型用对称边界条件是绝对正确的。
另外，对称边界条件得到的约束条件在载荷步里面查看到，约束是发生在环向的（柱坐标系下），环向位移约束为零。
加约束的时候直接加环向位移为零，也是可以的。
[attach]92985[/attach]
命令流部分
四、命令流（log文件另附）[attach]93111[/attach]
1、建模和求解部分（这里的建模网格划分比较密，可能不是很实用,这里的网格划
分不好，在裂纹尖端第一行单元没有奇异性，最好还是用kscon来做）：
/prep7
/COM,Structural
ET,1,PLANE2
KEYOPT,1,3,3
R,1,0.004,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1+011
MPDATA,PRXY,1,,0.3
MPDATA,DENS,1,,8500
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,5+008,6+009,,,,
R,1,0.004,
!
*
wpstyle,0.001,0.1,-1,1,0.003,0,2,,5
k,1,0,0,0
circle,1,0.032,,,90
/PNUM,KP,1
/PNUM,LINE,1
/PNUM,AREA,1
/NUMBER,0
/REPLOT
k,4,0.034,0,0
k,5,0,0.034,0
k,6,0.042,0,0
k,7,0,0.046,0
k,8,0.060,0,0
k,9,0,0.060,0
k,10,0.080,0,0
k,11,0,0.080,0
k,12,0.12,0,0
k,13,0,0.12,0
k,14,0.18,0,0
k,15,0,0.18,0
k,16,0.32,0,0
k,17,0,0.32
circle,1,0.5,,,90
l,2,4
l,4,6
l,6,8
l,8,10
l,10,12
l,12,14
l,14,16
l,16,18
l,3,5
l,5,7
l,7,9
l,9,11
l,11,13
l,13,15
l,15,17
l,17,19
SAVE
al,all
/COLOR,NUM,DGRA,1
/COLOR,NUM,BMAG,2
/COLOR,NUM,RED,3
/COLOR,NUM,CBLU,4
/COLOR,NUM,MRED,5
/COLOR,NUM,GREE,6
/COLOR,NUM,ORAN,7
/COLOR,NUM,MAGE,8
/COLOR,NUM,YGRE,9
/COLOR,NUM,BLUE,10
/COLOR,NUM,GCYA,11
/REPLOT
aplot
lesize,1,0.0004
lesize,2,0.04
lesize,3,0.0002
lesize,4,0.0002
lesize,5,0.0004
lesize,6,0.001
lesize,7,0.002
lesize,8,0.004
lesize,9,0.008
lesize,10,0.016
lesize,11,0.0002
lesize,12,0.0002
lesize,13,0.0004
lesize,14,0.001
lesize,15,0.002
lesize,16,0.004
lesize,17,0.008
lesize,18,0.016
MSHAPE,1,2D
MSHKEY,0
CM,_Y,AREA
CMSEL,S,_Y
AMESH,_Y
CMDELE,_Y
save
/SOLU
NSUBST,5,50,2
AUTOTS,1
sfl,2,pres,-1.5+008,-1.5+008
allsel,all
SFTRAN
lsel,,,,4,10,1
lsel,a,,,12,18,1
lplot
dl,all,,symm
allsel,all
sbctran
lswrite,1
!
*
NLGEOM,0
lsel,,,,2
sfl,2,pres,-1.2+009,-1.2+009,
allsel,all
sftran
NSUBST,15,1000,5
ARCLEN,1,0,0
lswrite,2
!
*
allsel,all
save
!
*
NLGEOM,0
sfl,2,pres,0.0,0.0,
allsel,all
sftran
OMEGA,0,0,62.831852,0
NSUBST,1,1,1
ARCLEN,0,0,0
lswrite,3
/GST,1
lssolve,1,2,1
2、J积分部分（对于半边裂纹，如果你已经定义了路径的话，直接把这部分命令流
输入进去就可以了）
/post1
!
local,11,0,0.034,0,0!
这里不应该建立局部坐标系。
只有计算应力强度因子才需要。
这
里只需要保证全局坐标系的X方向与裂纹平行就是了。
csys,0
！
这里应该有一个定义path，这里没有写出。
etable,volu,volu,
etable,sene,sene,
sexp,wden,sene,volu,1,-1,
pdef,wden,etab,w
升级会员 