断裂分析实例(J积分和应力强度因子)Word文档下载推荐.docx

断裂分析实例(J积分和应力强度因子)Word文档下载推荐.docx

《断裂分析实例(J积分和应力强度因子)Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《断裂分析实例(J积分和应力强度因子)Word文档下载推荐.docx（9页珍藏版）》请在冰豆网上搜索。

对圆盘的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,,,,

!

*

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

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

sbctran

lswrite,1

NLGEOM,0

lsel,,,,2

sfl,2,pres,-1.2+009,-1.2+009,

sftran

NSUBST,15,1000,5

ARCLEN,1,0,0

lswrite,2

sfl,2,pres,0.0,0.0,

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