三维裂纹J积分求解

一、建part
草图尺寸100,建草图为矩形(-10,50)、(50,-10)
拉伸厚度10,ok,得到六面体板
二、材料属性,具体参数自己设置吧,一步步都应该清楚。

3、装配(creat instance)
装配要用第二个(独立的),装配之后要切出一条缝,为后面定义seam,用下图右边的那个
这个是选择分割草图的一条边这里应该选择红色面的右边, 这个是分割需要的草图,黄色的线就是分割线,然后,分割tool-partition
选择cell-extrude/sweep edges
选择红色的圈
这个选择
选中的这个然后选择z轴,看好方向如果方向正确就选ok,如果不正确就选flip然后再选
按住shfit选择涂黄的两条线(就是裂纹线)
(1)先定义seam,special-creat-assign选择红色的区域,就定义了seam
(2)定义crake,special-creat-crake,用围线积分继续,这里选择
的时候容易选择的只是一个面
所以先隐藏那个小圆面,选择的那个就是隐藏的命令,隐藏小圆面后就可以选择crake front了。

Crake Front 亮红色的就是选择的裂纹线
学裂纹尖端有奇异性,所以这儿要按下图设置,不懂的看理论吧。

此时crake定义好了
5、荷载步
按静力学分析,继续默认设置就行,谈后在历史输出步中要编辑
Domain这个要选择之前定义的crake-1下面的围线数我一般用的是8,这个我感觉8左右就行,但是不能就2,3个,也不必要太多(路径无关性)。

6、load
这个我简单加载两端均布拉力,做好加一个约束边界条件,这里我加了一个右边的UZ
7、划分网格,
感觉断裂最麻烦的就是划分网格了,划分网格之前必须对实体分割能比较规则的,而且在裂纹尖端需要时放射状,其实尖端是四边形也能算,但是所有的例子都是放射状的。

因为划分网格比较麻烦我就补一步步的说了那个之前做的例子来示意一下
最好是按照下图进行分割,
矩形里面最好再分一个圆这样可以有效保证网格效果比较好。

分割之后就是布种子,布种子时可以约束选择第三个。

并且在厚度方向的线也要布种子,最好都布,而且数目相同,要不然会划分出来的网格不是一层一层的,求J积分的时候就会提示错误不能求J积分,但是不影响力学分析。

然后就是mesh contral黄色区域设置扫略,绿色区域设置结构话,再选择单元类型,这里用二阶单元,得到的网格就是下图这样。

当然了如果会用hypermesh的可以用hypermesh,本人还不会hypermesh所以只能这么笨的方法划分网格了。

然后就可以提交计算了。

辛苦了一中午为了保护自己的权益,写个自己的名字吧!
2015年5月13日星期三于大连理工大学
张鹏。

合集下载

含裂纹缺陷结构J积分计算方法研究

含裂纹缺陷结构J积分计算方法研究

含裂纹缺陷结构J积分计算方法研究裂纹缺陷是工程结构中常见的问题,在使用中,可能因为各种原因产生缺陷导致结构裂纹,从而影响其使用性能和寿命,甚至危及人员安全。

因此,需要对裂纹缺陷进行研究,并采用计算方法来评估其安全性,以达到保证工程质量和安全的目的。

其中,J积分计算方法是一种常用的结构裂纹评估方法。

J积分计算方法是基于线弹性理论的一种结构裂纹评估方法。

该方法将结构中的裂纹终端处碎片化为无限多的微小裂纹,然后根据线弹性理论建立裂纹尖端应力场解析模型,从而得到J积分的闭合积分式。

J积分是表示单位长度上的断裂能耗散率的物理量,是衡量结构裂纹扩展能力的重要参数。

当结构中有裂纹出现时,只要J值达到特定的值,就说明裂纹达到了稳定扩展状态,从而需要采取相应的应对措施。

J积分的计算方法主要分为两种:数值积分法和半解析法。

数值积分法是将J积分的积分式离散化成有限个节点,然后采用数值方法将积分转化为求和,最终得到计算结果。

数值积分法优点是计算精度高,适用范围广,对于复杂的结构修复方案可以进行精确的评估,同时计算也可以通过计算机程序实现自动化。

但是数值积分法也有缺点,例如计算复杂性高、计算量大、计算时间长等,因此对于大型的结构,一般采用半解析法来进行计算。

半解析法是对J积分的积分式进行一定程度的简化,将其拆分成不同的组份,然后将其进行勾股定理、互变性等数学变换,最终将不定积分转化为常数的形式,并进行简化。

这种方法不仅可以降低计算难度,还可以在一定程度上提高计算速度。

半解析法适用于简单结构的定量评估,例如板、梁的断裂评估,但对于复杂结构的计算,依然需要采用数值积分法进行计算。

总的来说,J积分计算方法是一种可靠、准确的结构裂纹评估方法。

该方法不仅可以定量评估结构裂纹的安全性和扩展性,而且还能够为工程设计提供有效的指导和支持。

同时,随着计算机技术的不断提高和发展,J积分计算方法的计算精度和速度也将不断提高,为工程结构的安全性和可靠性提供更为有力的保障。

05--积分变换方法与三维裂纹问题

05--积分变换方法与三维裂纹问题

u 1
2 E
(1
)
d 2 dy2
2
ei
x
i
d
7 10
v 1
2 E
(1
)
d 3 dy3
(2
) 2
d dx
ei
x
1
2
d
7 11
若求得 ( , y) ,可得 u(x, y) , v(x, y) 。 6
二、半无限弹性平面的位移解
现讨论受分布压力 q(x)的半无限弹性平面问题
边界条件为:
(1) y 0
F
(
)
2
q0a
1/
2
J1(
)
可得位移:
v0
(x)
2q0
(1
E
2
)
a2 x2
在均布压力作用下,裂纹会扩大张开成椭圆形状。 利用这种方法可解许多种裂纹尖端的应力位移场。
15
• 三维裂纹问题的求解
16
受均匀拉伸的椭圆盘状裂纹,Green-Sneddon解
z 0,
x2 a2
y2 b2
1
边界条件:
x2 y2 z2 ij ~ 0
7 22
q( ) cos( x) d 0 ( x a)
0
7 23

引入代换: , q( )
aa
F() , cos()
2
J
1 2
(
)
• 式中 J 是贝塞尔函数 12
三、裂纹问题的对偶积分方程
利用上述代换,边界条件(7-22)、(7-23)写为:
0 F ()J1 ()d g() 2
z 0, (x, y) D zz p0 , xz yz 0

三维裂纹J积分研究

三维裂纹J积分研究
1 三维 裂纹 J积 分 围线积 分 定义 式
最早 由 Ri 提 出的 . 分 围线 积 分 定 义式 是 c e ,积
基 于平 面问题 的 ( 1 : 图 )
过程 中占有重 要地 位 。
J积 分一 般 通 过 有 限元 分 析 求 得 , 了适 合 工 为
-j d ( , 『 塞+ y ( 一- W T ] 1 塞) )
算在 制订 弹 塑性评 定规 范 和具体 缺 陷 的弹塑 性评 定
裂 纹 的 其 它 形 式 以 及 长 屈 服 平 台 的 三 维 - 分 ,积
解 。但 是 , 些 都 是 用 虚 裂 纹 扩 展 法计 算 三 维裂 ] 这 纹 -积分 , , 而笔 者 则借 助 于 三 维裂 纹 J积 分 的 围线 积分 定 义式 编制 程序并 计 算三 维裂 纹 - 分 。 ,积
三维 裂 纹 积 分 研 究
吴龙 平 明斐卿 , ,李 国成 。 ,黄
(. 国石 油 天 然 气 管 道 工 程 有 限 公 司 ,河 北 廊 坊 1中
河北 廊坊

276 ) 5 0 1
0 5 0 ;2 河 北 石 油 职业 技 术 学 院 , 60 0 .
0 5 0 }3 中 国 石 油 大 学 ,山 东 东 营 600 .
维普资讯
第 3 5卷
第 3期






V0 5 N0 3 L 3 .
Ma 2 6 v 00
20 0 6年 5月
PETR(Байду номын сангаас EM I )CH CAI QUI ENT E PM
文 章 编 号 :10 —4 6 2 0 ) 30 1 -4 0 07 6 ( o6 O —0 40

基于MARC的某火炮身管内表面裂纹三维J积分有限元模型

基于MARC的某火炮身管内表面裂纹三维J积分有限元模型

基 于 MARC 的某 火 炮 身 管 内表 面裂 纹 三 维 J积 分 有 限元模 型
周敏 华 。 。李 晓谦 。韦
(.中 南 大 学 机 电工 程 学 院 ,湖 南 长 沙 1

4 08 ;2 1 0 3 .中 国人 民 解放 军 7 7 0部 队 ,广 西 鹿 寨 5 5 0 ) 55 4 6 6
g a , t r tc 1sz ft r c s c lul t d u e h on to fma i u c m b rpr s ur . Th s r l he e iia ieo hec a k wa a c a e nd rt e c dii n o x m m ha e e s e i me h d ha e t r e i e i g a plc ton va u sf r s c iy e l a i nd pr d c i n o he g r e. t o sgr a e ng ne rn p ia i l e o e urt va u ton a e i to ft un ba r 1
通过一条与裂纹前缘重合的曲线设定其周围节点到这条曲线的距离容差保证在裂纹前缘的一个位置只能有一个节点将裂纹前缘附近的节点粘附在这条曲线上然后利用这条曲线选取附着在其上边的节点定义为裂纹前缘离散节点如图所示
火 炮 发 射 与控 制 学 报 第 1期
J OU RNAL OF GUN LAUNCH & CONTROL ・ 9 ・ 3
元 模 型 的建 立 方 法 。通 过 J积 分 值 随 裂 纹 长 度 及 深 度 的 变 化 规 律 。根 据 某
火炮 身管 裂 纹 扩 展 的 临界 J 分 值 ,求 出 了其 在 最 大 膛 压 作 用 下 临 界 裂 纹 尺 寸 ,这 对 火 炮 身 管 的 安 全 性 评 估 积

任意三维裂纹扩展分析

任意三维裂纹扩展分析

任意三维疲劳裂纹扩展分析1.前言在工程实际中,真实的结构总是存在众多缺陷或裂纹,对于一个含裂纹或缺陷的构件,多在其服役荷载远低于容许强度的情况下就发生了破坏。

实际工程结构在经受长时间多因素综合作用下,产生变形、裂纹等缺陷,从而导致整个结构的失效。

结构的失效主要由疲劳引起,其最终失效形式即为断裂,有大约80%以上的工程结构的断裂与疲劳有关,由疲劳引起的巨大经济损失及灾难性的后果不胜枚举。

我们通常不能仅仅因为某个构件出现了裂纹就简单的认为该构件不安全或不可靠,尤其是对于大型设备的重要构件,因为这将使企业耗费高昂的成本。

对于出现的裂纹,以往多采用以下几种处理办法:一是对出现裂纹的构件进行更换,这对于含裂纹但仍能工作的构件是一个巨大的浪费。

二是强行停止使用进行维修,这样会带来巨大的经济损失;三是冒险继续使用,但这样会带来巨大风险,甚至会造成人员伤亡。

所以,人们更想知道,出现的裂纹是否会在既定载荷(包括疲劳载荷在内的任意载荷)下扩展成不安全或失效的临界尺寸,因此,出现了疲劳裂纹扩展分析。

疲劳裂纹扩展分析是采用断裂力学的理论和方法对含裂纹等缺陷构件的失效过程进行分析,以评估产品的安全性和可靠性,可以进行损伤容限评估和剩余寿命预测等,已经在化工机械、飞行器、核工业等各个工程领域得到了广泛应用,并得到了世界各国政府及学术机构的重视。

2.疲劳裂纹扩展分析软件在工程实践中,疲劳裂纹扩展分析已成为评估产品性能、改良产品设计和提高服役寿命的一个重要工具。

目前,疲劳裂纹扩展分析主要有解析法和数值法这样两种方法,下面分别介绍这两种方法。

1)解析法解析法主要依据相应的规范和经验公式,将复杂的三维问题简化为二维问题,并对复杂的裂纹形状和荷载状态进行简化,然后用经验的方法对裂纹安全性进行评估。

但对于大量结构复杂的工程实际问题却无能为力,况且其简化后的分析准确度及是否真实逼近服役情况也值得探讨。

目前,工程上有几款基于解析法而开发的裂纹扩展分析软件,它们主要应用于航空标准结构的裂纹扩展分析,包括DARWIN、NASGRO、AFGROW等。

断裂力学参量[整理版]

断裂力学参量[整理版]

ANSYS求解断裂力学参量的理论方法工程上,线弹性断裂力学中常用应力强度因子K、J积分、G能量释放率这三个参量来描述裂纹场。

ANSYS软件能较好地计算裂纹周围区域的应力分布,并能计算裂纹的应力强度因子K、J积分以及能量释放率G等,其特点是简单、经济、精度高。

下面主要介绍在ANSYS中如何求解应力强度因子K和J积分。

(1)求解应力强度因子ANSYS软件中提供了所谓的“位移外推”法(displacement extrapolation) 来计算应力强度因子[5]。

在线弹性范围内,对于三维裂纹,裂纹尖端的局部位移场与应力强度因子的关系为[6]:)2)22IIIIIIKu kGKv kGKwG⎧=+⎪⎪⎪⎪=+⎨⎪⎪⎪=⎪⎩式中: u、v、w—如图2.5所示裂纹尖端局部直角坐标系下裂纹前端位移;r—如图2.5所示裂纹尖端局部柱坐标系下坐标;G—材料剪切模量;K I、K II、K III—应力强度因子;v—为泊松比;34()3()1vk vv-⎧⎪=⎨-⎪+⎩平面应变或轴对称平面应力当利用裂纹尖端节点的位移进行计算时,应力强度因子和裂纹面节点的位移差存在下列关系:IIIIIIKKK⎧=⎪⎪⎪⎪=⎨⎪⎪⎪=⎪⎩三维裂纹的局部坐标在使用有限元法进行应力强度因子计算时,由于常规单元在裂纹尖端存在奇异性,为使计算准确,必须在裂纹尖端使用细小的单元;如果使用奇异元,即使用二次三角(或五面体)单元,并将靠近裂纹尖端的中间节点置于1/4处,则位于沿裂纹尖端的单元边上的应力和应变与1/消除了奇异性,也就是说,可以用相对比较稀疏的单元得到精度较高的结果。

(2)求解J积分J积分定义为一个围绕裂尖的线积分(二维) 或一个围绕裂纹前沿的面积分。

它用计算裂纹尖端的奇异应力和应变,与积分路径无关。

为了避开裂纹尖点的奇异性,取得较好的精度,积分路径一般取得离裂纹尖点较远。

J积分形式如图2.6所示,其表达式如下:()yxx yuuJ Wdy t t dsx yΓΓ∂∂=-+∂∂⎰⎰式中:W—应变能密度(单位体积应变能);Г—围绕裂纹尖点任意路径;xt—X 方向的作用向量,x x xy yt nσσ=+;yt—Y方向的作用向量,y y xy xt nσσ=+;n—积分路径的外法向向量;s —积分路径距离;围绕裂纹尖端的任意一条J 积分路径在ANSYS 中,为了计算位移向量的偏导数x u x ∂∂与y u y ∂∂,将积分路径向x 正负方向分别移动Δx/2,并求出路径Γ+Δx/2上各点的位移u x1和u y 1以及路径Γ-Δx/2上各点的u x 1和u y 1,则:2121()()x x x y y y u x u u xu y u u y∂∂=-∆⎧⎪⎨∂∂=-∆⎪⎩ ANSYS 具有强大的后处理功能,利用此功能,在求解后可以通过ANSYS 通用后处理器中的单元列表功能,很方便地把各变量映射到自定义的路径中去。

三维J积分在ANSYS中的计算

采取有效的检验方法,其检验周期可在目前维修周期基础上延长。

由于设备的检验方法仅对设备失效可能性产生影响,而对失效后果没有影响,因此,除了降低设备的失效概率之外,还要采取其它措施减少失效后果,以降低风险。

参考文献:[1] 马秋宁,苗立新,蔡有军.炼油化工装置风险评价模式的探讨[J ].当代化工,2001,30(3):1572159.[2] 杨铁成,陈学东,魏安安,等.基于半定量风险分析的加氢装置安全评估[J ].压力容器,2002,19(12):43245.[3] 陈学东,杨铁成,艾志斌,等.基于风险的检测(RBI )在实践中若干问题讨论[J ].压力容器,2005,22(7):36244.[4] American Petroleum Institute ,Risk 2Based Inspection ,BaseResource Document API 581[S].[5] 赵 涛1基于风险的乙烯裂解炉装置维修规划研究[D ]1南京:南京工业大学,2005.[6] Soudabeh A Noori ,John W H Price.Case Study of t he use ofAPI581on H K and HP Material Furnace Tubes[J ].Journal of Pressure Vessel Techonlogy ,2005,127(1):49253.(杜编) 收稿日期:2005210209作者简介:宋启明(19742),男(汉族),辽宁抚顺人,在职硕士研究生,从事压力容器安全评定研究工作。

文章编号:100027466(2006)022*******三维J 积分在ANSYS 中的计算宋启明,李国成(中国石油大学机电学院,山东东营 257061)摘要:利用有限元软件ANSYS 中的二次开发工具A PDL 语言编程,实现对多种典型含裂纹构件的三维弹塑性断裂计算和分析的自动化。

固体火箭发动机药柱裂纹的J积分分析

固体火箭发动机药柱裂纹的J积分分析蒙上阳;胡光宇;刘兵;彭威【摘要】为了探讨某固体火箭发动机药柱纵向裂纹在点火增压时的稳定性,采用三维粘弹性有限元法,在三维J积分圆柱围道曲面内裂纹尖端,构建奇异三维裂纹元,提高模拟精度,分别计算了随裂纹扩展所对应裂纹深度的J积分,并根据J积分随裂纹度变化规律,探讨裂纹的稳定性.研究表明,发动机点火发射时药柱前翼槽出现的纵向裂纹最为危险.【期刊名称】《固体火箭技术》【年(卷),期】2010(033)006【总页数】4页(P646-649)【关键词】固体火箭发动机;有限元法;奇异裂纹元;三维J积分【作者】蒙上阳;胡光宇;刘兵;彭威【作者单位】中国人民解放军63961部队,北京,100012;中国人民解放军63961部队,北京,100012;中国人民解放军63961部队,北京,100012;军械工程学院,石家庄,050003【正文语种】中文【中图分类】V4350 引言在固体火箭发动机点火时,药柱主要承受燃气内压与轴向过载的作用,含有裂纹等缺陷的发动机在点火时易出现事故。

因为发动机点火增压时,燃气可能进入裂纹腔内,导致原有裂纹扩展;若裂纹失稳扩展,会导致发动机发生穿火或轰爆等灾难性事故。

因此,获取控制推进剂药柱裂纹开裂的物理参量,判断裂纹的稳定性,一直是固体推进技术领域最为关心的研究课题之一。

在点火增压过程中,发动机药柱始终处于高温、高压的严酷环境,且这一过程难以重现,试图通过点火试验测试发动机药柱裂纹的稳定性十分困难。

随着数值仿真技术的发展,采用数值方法计算推进剂药柱裂纹开裂的物理参量取得了长足进步。

数值方法不受发动机的结构形状、承受载荷种类及药柱裂纹分布等方面限制,可大大减少物理样机试验,节约研发经费。

固体火箭发动机药柱裂纹数值仿真技术发展重要方向是提高计算精度,不仅通过构建平面奇异单元提高计算平面裂纹的应力强度因子的精度[1-2],而且为进一步模拟发动机的真实响应,通过构建发动机的三维有限元模型获取推进剂药柱三维裂纹的断裂参量[3-6]。

浅谈J积分法在测定材料断裂韧性中的应用

浅谈J积分法在测定材料断裂韧性中的应用J积分法是一种用于测定材料断裂韧性的重要方法。

它是基于应变能耗散的原理,通过对材料断裂过程中应变能的积分来评估其断裂韧性。

本文将从J积分法的基本原理、测定方法以及应用实例三个方面进行论述,以期对J积分法在测定材料断裂韧性中的应用进行深入探讨。

首先,我们来介绍J积分法的基本原理。

J积分法是建立在线弹性力学理论基础之上的,假设材料的断裂过程是一个无限小的裂纹扩展过程。

在无限小裂纹扩展中,裂纹尖端附近的应变场可以近似看作是一个二维弹性问题。

根据线弹性力学理论,材料在拉伸过程中的应变能耗散可以表示为应变能密度的分布累积,即应变能耗散密度的积分。

J积分即是该应变能耗散密度积分的结果。

其次,我们来讨论J积分法的测定方法。

J积分的计算通过对裂纹尖端附近应变场的数值模拟来实现。

一种常用的计算方法是有限元方法,通过构建适当的数值模型来计算应变能耗散密度的积分。

这需要将裂纹尖端的应力和应变信息输入到数值模型中,然后进行数值计算得到J积分的结果。

J积分法的测定还可以通过实验手段实现,常用的方法是通过测量断裂试样的载荷-位移曲线,并根据曲线的降低率计算出断裂韧性参数。

具体而言,可以通过剪切试验或拉伸试验等实验方法,先测定断裂试样的载荷-位移曲线,然后通过求解该曲线的斜率来获得J积分的结果。

最后,我们来探讨J积分法在测定材料断裂韧性中的应用。

J 积分法具有很大的应用潜力,尤其在评估材料的断裂性能方面具有重要意义。

通过测定材料的J积分值,可以得到较为准确的断裂韧性参数,对于评估材料的断裂性能、设计结构的可靠性以及材料的损伤机理研究等方面具有重要意义。

J积分法的应用非常广泛,例如在航空航天领域,评估材料的断裂韧性是确保飞机结构安全可靠的关键因素之一。

通过J积分法可以对航空材料的断裂性能进行准确评估,并指导材料的选用和结构设计。

此外,在汽车工业、核能工业、石油化工等领域,J积分法也被广泛应用于材料断裂韧性的测定,以保证关键设备的安全运行。

Ⅰ型和Ⅲ型复合裂纹J积分计算


J I n t e g r a l C a l c u l a t i o n o f I 。 Ⅲ Mi x e d Mo d e Cr a c k
C A I Y o n g - me i ~, WA N G We i , X I E Y u  ̄ u n
在实 际工 程结 构 中 , 裂 纹 多 数 处 于 复 合 型 变形 状态 , 例如 , 大型 旋转构 件 的裂纹 往 往是张 开 型和扭
裂纹 复合 和裂 尖 单 元 的 奇异 行 为 , J积 分 不 同于 其 它 方法 的 特 点 在 于 积 分 具 有 与 所 取 回路 无 关 的 特 性 ] , 从 而可 以避免 去分 析裂 纹尖 端 附近塑 性 区 的
第4 2 卷 第 2 期 2 0 1 3年 3月
文 章 编 号 :1 0 0 0 — 7 4 6 6 ( 2 0 1 3 ) 0 2 — 0 0 3 5 — 0 4






Vo 1 . 42 NO .2
Ma r .2 O1 3
Hale Waihona Puke P ETRO— CHEM I CAL EQUI P M ENT
t e g r a l e v a l u a t i o n i n I +l l I mi x e d c r a c k i s b a s e d o n t h e c o n s e r v a t i o n i n t e g r a l me t h o d wh i c h s h o ws
关 键 词 :复合裂纹 ;J积分 ;应力强度 因子 ;数值计算
中图分 类号 :O 3 4 6 . 1 ;TB 1 2 5
文献标 志 码 : A
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档