CN105574282A - 3d裂纹起裂位置的预测方法 - Google Patents
3d裂纹起裂位置的预测方法 Download PDFInfo
- Publication number
- CN105574282A CN105574282A CN201510988163.XA CN201510988163A CN105574282A CN 105574282 A CN105574282 A CN 105574282A CN 201510988163 A CN201510988163 A CN 201510988163A CN 105574282 A CN105574282 A CN 105574282A
- Authority
- CN
- China
- Prior art keywords
- crackle
- crack
- ref
- finite element
- location
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种3D裂纹起裂位置的预测方法。因为在预测裂纹的起裂位置时已纳入拘束效应(即考虑了裂纹的尺寸与实际工况的影响),所以,本发明所涉及的3D裂纹起裂位置的预测方法可以准确的评定3D裂纹起裂位置。本发明提供的方法简单、成本低,便于工程实践中应用,可在实际测量中准确的预测3D裂纹起裂位置,为下一步结构完整性评定准确地提供依据。
Description
技术领域
本发明涉及一种评定结构完整性的方法,具体涉及一种三维(3D)裂纹起裂位置的预测方法。
背景技术
在实际结构的制造和使用过程中,总会不可避免的产生3D裂纹。3D裂纹一旦产生,将严重影响结构的使用寿命。为了避免事故的发生,在发现3D裂纹后,需要对裂纹进行结构完整性评定与跟踪监测,这就需要对3D裂纹的起裂位置(最危险位置)进行预测。3D裂纹的形状一般用半椭圆描述。
一般而言,3D裂纹的J积分最大处断裂驱动力最大,所以在工程实践中,一般选用J积分最大处作为最危险位置进行重点监测。然而,3D裂纹的起裂不仅和驱动力有关,还和3D裂纹的断裂阻力有关。大量实验和研究表明,裂纹的断裂阻力受试样/结构几何、裂纹尺寸和加载方式等因素的影响,这种影响通常被称为“拘束效应”,拘束的增加导致材料断裂阻力的降低。截至目前,仍未有一种简单的方法在预测3D裂纹的起裂位置时纳入拘束效应。
发明内容
本发明是为了解决上述课题而进行的,目的在于提供一种纳入拘束效应的3D裂纹起裂位置的预测方法,该方法综合考虑裂纹驱动力与断裂阻力。
本发明提供了一种3D裂纹起裂位置的预测方法,具有这样的特征,包括如下步骤:
(1)通过ABAQUS等有限元软件针对含裂纹的结构建立有限元模型,对有限元模型的3D裂纹的不同位置处的J积分进行模拟计算,其中J积分为弹塑性断裂力学中一个与路径无关的积分;同时,针对J积分引入Jref参数,Jref为材料或结构的断裂韧性;
(2)对有限元模型的3D裂纹的不同位置处的进行计算,Ap=APEEQ/Aref,
其中,Ap为描述拘束效应的拘束参数,APEEQ为等效塑性应变(PEEQ)等值线所围绕区域的面积,Aref为Jref下PEEQ等值线所围绕区域的面积;
(3)对有限元模型的3D裂纹前沿不同位置处的进行计算,计算得到的的最大值处为3D裂纹的起裂位置。
在本发明提供的3D裂纹起裂位置的预测方法中,还可以具有这样的特征:其中,同一材料的Jref、Aref均设置为常数1。
在本发明提供的3D裂纹起裂位置的预测方法中,还可以具有这样的特征:其中,将3D裂纹近似为半椭圆形,并根据裂纹结构对称性,选取二分之一或四分之一的3D裂纹结构建立有限元模型。
发明的作用和效果
本发明所涉及的3D裂纹起裂位置的预测方法,因为在预测裂纹的起裂位置时已纳入拘束效应(即考虑了裂纹的尺寸与实际工况的影响),所以,本发明所涉及的3D裂纹起裂位置的预测方法可以准确的评定3D裂纹起裂位置。本发明提供的方法简单、成本低,便于工程实践中应用,可在实际测量中准确的预测3D裂纹起裂位置,为下一步结构完整性评定准确地提供依据。
附图说明
图1(a)是本发明的实施例3D裂纹起裂位置的预测方法中板结构的拉伸图,图1(b)是本发明的实施例3D裂纹起裂位置的预测方法中板结构的内部裂纹;
图2(a)是本发明的实施例3D裂纹起裂位置的预测方法中建立的有限元模型,图2(b)是本发明的实施例3D裂纹起裂位置的预测方法中有限元模型的裂纹面划分网格,图2(c)是本发明的实施例3D裂纹起裂位置的预测方法中有限元模型图的裂尖部分c的局部网格放大图;
图3(a)是本发明的比较例3D裂纹起裂位置的预测方法中沿裂纹前沿的计算得到的J积分分布图,图3(b)是本发明的实施例3D裂纹起裂位置的预测方法中沿裂纹前沿的计算得到的J积分分布图;
图4是本发明的实施例3D裂纹起裂位置的预测方法中沿裂纹前沿测出的拘束参数的分布;以及
图5是本发明的实施例3D裂纹起裂位置的预测方法中测得的纳入拘束效应后裂纹的起裂位置。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,以下实施例结合附图对本发明3D裂纹起裂位置的预测方法作具体阐述。
图1(a)是本实施例3D裂纹起裂位置的预测方法中板结构的拉伸图,图1(b)是本实施例3D裂纹起裂位置的预测方法中板结构的内部裂纹。
如图1所示,本实施例以板结构进行模拟。3D裂纹的材料选用D6AC钢。a、c分别为半椭圆裂纹的短、长半轴,φ为过裂纹周线上任一点的径向线与长轴间的夹角;t为板结构的厚度,W为板结构的宽度;S为拉伸载荷。
图2(a)是本实施例3D裂纹起裂位置的预测方法中建立的有限元模型,图2(b)是本实施例3D裂纹起裂位置的预测方法中有限元模型的裂纹面划分网格,图2(c)是本实施例3D裂纹起裂位置的预测方法中有限元模型图的裂尖部分c的局部网格放大图。
如图2所示,采用ABAQUS/Standard软件针对含3D裂纹的板结构进行建模,根据3D裂纹的对称性,将含裂纹的板结构均分成四份,选取其中一份板结构建立有限元模型,裂纹前端采用较细的网格划分。
图3(a)是比较例中3D裂纹起裂位置的预测方法中沿裂纹前沿的计算得到的J积分分布图,图3(b)是本实施例3D裂纹起裂位置的预测方法中沿裂纹前沿的计算得到的J积分分布图。
图3(a)中,比较例中裂纹深度a/t=0.65,形状比a/c=0.61,拉伸载荷S=1006.3MPa,现有方法不纳入拘束效应预测的裂纹起裂位置为2φ/π=0.30,与实际测量的裂纹起裂位置2φ/π=0.34相差0.04弧度即2.29度,相差较小。
可在本实施例(图3(b))中,裂纹深度a/t=0.19,形状比a/c=0.76,拉伸载荷S=1300.8MPa,现有方法不纳入拘束效应预测的裂纹起裂位置为2φ/π=1,与实际测量位置2φ/π=0.68相差0.32弧度即18.33度,相差较大,如果在实际工程中很可能会埋下巨大的安全隐患。
图4是本实施例3D裂纹起裂位置的预测方法中沿裂纹前沿测出的拘束参数的分布;图5是本实施例3D裂纹起裂位置的预测方法中测得的纳入拘束效应后裂纹的起裂位置。
如图4和图5所示,为减小预测值与实际测量的误差,本实施例采用以下步骤:
(1)通过ABAQUS软件针对含裂纹的结构建立有限元模型,对有限元模型的3D裂纹的不同位置处的J积分进行模拟计算,其中J积分为弹塑性断裂力学中一个与路径无关的积分;同时,针对J积分引入Jref参数,Jref为材料或结构的断裂韧性;
(2)对有限元模型的3D裂纹的不同位置处的进行计算,Ap=APEEQ/Aref,
其中,Ap为描述拘束效应的拘束参数,APEEQ为等效塑性应变(PEEQ)等值线所围绕区域的面积,Aref为Jref下PEEQ等值线所围绕区域的面积;
(3)对有限元模型的3D裂纹前沿不同位置处的进行计算,计算得到的的最大值处为3D裂纹的起裂位置。
通过本实施例的方法计算预测的裂纹起裂位置为2φ/π=0.6,与实际测量位置2φ/π=0.68相差0.08弧度即4.58度,预测结果与实际测量结果相差较小,比之现有的方法大幅度的提高了精准度。本实施例中,引入Jref与Aref参数的目的是消除量纲。因为同一材料而言,Jref、Aref为一个定量。为简化计算,可将同一材料的Jref与Aref参数均作为常数1,且不影响预测结果。
实施例的作用与效果
本发明所涉及的3D裂纹起裂位置的预测方法,因为在预测裂纹的起裂位置时已纳入拘束效应(即考虑了裂纹的尺寸与实际工况的影响),所以,本发明所涉及的3D裂纹起裂位置的预测方法可以准确的评定3D裂纹起裂位置。本发明提供的方法简单、成本低,便于工程实践中应用,可在实际测量中准确的预测3D裂纹起裂位置,为下一步结构完整性评定准确地提供依据。
以上实施例仅为本发明构思下的基本说明,不对本发明进行限制。而依据本发明的技术方案所作的任何等效变换,均属于本发明的保护范围。
Claims (3)
1.一种3D裂纹起裂位置的预测方法,用于预测结构中不规则3D裂纹的起裂位置,其特征在于,包括如下步骤:
(1)通过ABAQUS有限元软件针对含裂纹的结构建立有限元模型,对所述有限元模型的3D裂纹的不同位置处的J积分进行模拟计算,其中J积分为弹塑性断裂力学中一个与路径无关的积分;同时,针对J积分引入Jref参数,Jref为材料或所述结构的断裂韧性;
(2)对有限元模型的3D裂纹的不同位置处的进行计算,Ap=APEEQ/Aref,
其中,Ap为描述拘束效应的拘束参数,APEEQ为等效塑性应变(PEEQ)等值线所围绕区域的面积,Aref为Jref下PEEQ等值线所围绕区域的面积;
(3)对有限元模型的3D裂纹前沿不同位置处的进行计算,计算得到的所述的最大值处为所述3D裂纹的起裂位置。
2.根据权利要求1所述的3D裂纹起裂位置的预测方法,其特征在于:
其中,同一材料的Jref、Aref均设置为常数1。
3.根据权利要求1所述的3D裂纹起裂位置的预测方法,其特征在于:
其中,将所述3D裂纹近似为半椭圆形,并根据裂纹的结构对称性,选取二分之一或四分之一的3D裂纹结构建立所述有限元模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510988163.XA CN105574282A (zh) | 2015-12-24 | 2015-12-24 | 3d裂纹起裂位置的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510988163.XA CN105574282A (zh) | 2015-12-24 | 2015-12-24 | 3d裂纹起裂位置的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105574282A true CN105574282A (zh) | 2016-05-11 |
Family
ID=55884412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510988163.XA Pending CN105574282A (zh) | 2015-12-24 | 2015-12-24 | 3d裂纹起裂位置的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105574282A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108918298A (zh) * | 2018-08-10 | 2018-11-30 | 西南交通大学 | 三维裂尖约束效应、等效应力场的表征方法及评估金属结构延性断裂韧性的方法 |
CN111539144A (zh) * | 2020-04-20 | 2020-08-14 | 中车青岛四方机车车辆股份有限公司 | 一种含裂纹结构件的断裂韧性计算方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102136018A (zh) * | 2011-03-17 | 2011-07-27 | 哈尔滨工业大学 | 焊缝含中心裂纹受拉伸载荷对接接头实现等承载的设计方法及k因子的应用 |
CN102176212A (zh) * | 2010-12-30 | 2011-09-07 | 哈尔滨工业大学 | 含i型中心裂纹对接接头的应力强度因子的推导方法及应用 |
CN103116682A (zh) * | 2013-03-13 | 2013-05-22 | 空军工程大学 | 一种基于有限元软件的孔边三维裂纹建模方法 |
CN103471940A (zh) * | 2013-09-27 | 2013-12-25 | 华东理工大学 | 面内及面外统一拘束与材料延性断裂韧性关系的确定方法 |
CN104504185A (zh) * | 2014-12-10 | 2015-04-08 | 西北工业大学 | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 |
-
2015
- 2015-12-24 CN CN201510988163.XA patent/CN105574282A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176212A (zh) * | 2010-12-30 | 2011-09-07 | 哈尔滨工业大学 | 含i型中心裂纹对接接头的应力强度因子的推导方法及应用 |
CN102136018A (zh) * | 2011-03-17 | 2011-07-27 | 哈尔滨工业大学 | 焊缝含中心裂纹受拉伸载荷对接接头实现等承载的设计方法及k因子的应用 |
CN103116682A (zh) * | 2013-03-13 | 2013-05-22 | 空军工程大学 | 一种基于有限元软件的孔边三维裂纹建模方法 |
CN103471940A (zh) * | 2013-09-27 | 2013-12-25 | 华东理工大学 | 面内及面外统一拘束与材料延性断裂韧性关系的确定方法 |
CN104504185A (zh) * | 2014-12-10 | 2015-04-08 | 西北工业大学 | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 |
Non-Patent Citations (5)
Title |
---|
YANG J, WANG G Z, XUAN F Z, ET AL.: "Unified correlation of in‐plane and out‐of‐plane constraints with fracture toughness", 《FATIGUE & FRACTURE OF ENGINEERING MATERIALS & STRUCTURES》 * |
关凯书,徐一飞等: "带环形裂纹的微试样小冲杆试验法测试材料断裂韧性探讨", 《塑性工程学报》 * |
张国尚,荆洪阳等: "焊接接头脆性断裂评定的局部法研究进展", 《压力容器》 * |
杨杰,王国珍,轩福贞,涂善东: "基于裂尖等效塑性应变的面内与面外统一拘束参数的研究", 《核技术》 * |
杨杰: "面内_面外统一拘束参数及其与材料及焊接接头断裂韧性的关联", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108918298A (zh) * | 2018-08-10 | 2018-11-30 | 西南交通大学 | 三维裂尖约束效应、等效应力场的表征方法及评估金属结构延性断裂韧性的方法 |
CN108918298B (zh) * | 2018-08-10 | 2020-03-06 | 西南交通大学 | 三维裂尖约束效应及三维裂尖等效应力场的表征方法 |
CN111539144A (zh) * | 2020-04-20 | 2020-08-14 | 中车青岛四方机车车辆股份有限公司 | 一种含裂纹结构件的断裂韧性计算方法及*** |
CN111539144B (zh) * | 2020-04-20 | 2023-09-22 | 中车青岛四方机车车辆股份有限公司 | 一种含裂纹结构件的断裂韧性计算方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Livieri | Use of J-integral to predict static failures in sharp V-notches and rounded U-notches | |
Liu et al. | Buckling failure mode analysis of buried X80 steel gas pipeline under reverse fault displacement | |
Lee et al. | Fracture prediction of thin plates under hemi-spherical punch with calibration and experimental verification | |
TWI391657B (zh) | A rupture determination method, a rupture determination apparatus, a program, and a computer-readable recording medium | |
Abbasi et al. | New attempt to wrinkling behavior analysis of tailor welded blanks during the deep drawing process | |
Kramberger et al. | Assessing the cracking behavior of auxetic cellular structures by using both a numerical and an experimental approach | |
CN105115821A (zh) | 一种基于有限元的材料断裂韧性确定方法 | |
Nourpanah et al. | Ductile crack growth and constraint in pipelines subject to combined loadings | |
Wei et al. | Modeling of mixed-mode crack growth in ductile thin sheets under combined in-plane and out-of-plane loading | |
Kibey et al. | Tensile strain capacity equations for strain-based design of welded pipelines | |
UA115980C2 (uk) | Система і спосіб прогнозу прогину з проклацуванням панелей з формованих сталевих листів | |
Liao et al. | Applicability of damage models for failure analysis of threaded bolts | |
CN105574282A (zh) | 3d裂纹起裂位置的预测方法 | |
Hu et al. | Failure analysis based on J-integral values: A case study of hydrogen blistering defect | |
Park et al. | Evaluation of stress and crack behavior using the extended finite element method in the composite layer of a type III hydrogen storage vessel | |
Hradil et al. | Ductility limits of high strength steels | |
Islam et al. | Experimental verification of automated design of reinforced concrete deep beams | |
Costa et al. | Numerical simulation of quasi-static compression behavior of the toe cap component for safety footwear | |
Peng et al. | Three-dimensional fracture analysis of large opening box girder with crack damage under bending and torsion loads | |
Zhang et al. | Analytical method for evaluating the impact response of stiffeners in a ship side shell subjected to bulbous bow collision | |
Guo et al. | Research on deformation control technology of back plate based on vibration aging | |
Nicak et al. | Atlas plus: Design of large scale fracture mechanics tests on a ferritic pipe | |
Slamene et al. | Damage analysis of functionally graded materials: A finite element investigation utilizing the Gurson–Tvergaard–Needleman (GTN) model for notched plates | |
Hertelé et al. | Sensitivity of plastic response of defective pipeline girth welds to the stress-strain behavior of base and weld metal | |
Meleshko | Hybrid model for determining the stiffness of a spatial rod with consideration geometric nonlinearity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160511 |