CN116738780B - 考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 - Google Patents
考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 Download PDFInfo
- Publication number
- CN116738780B CN116738780B CN202310609963.0A CN202310609963A CN116738780B CN 116738780 B CN116738780 B CN 116738780B CN 202310609963 A CN202310609963 A CN 202310609963A CN 116738780 B CN116738780 B CN 116738780B
- Authority
- CN
- China
- Prior art keywords
- crack
- deflection
- length
- sample
- stress intensity
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims abstract description 13
- 230000008878 coupling Effects 0.000 claims abstract description 7
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 238000005859 coupling reaction Methods 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 6
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 8
- 229910000861 Mg alloy Inorganic materials 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000009864 tensile test Methods 0.000 description 3
- 238000005336 cracking Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 229910001297 Zn alloy Inorganic materials 0.000 description 1
- PGTXKIZLOWULDJ-UHFFFAOYSA-N [Mg].[Zn] Chemical compound [Mg].[Zn] PGTXKIZLOWULDJ-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910052749 magnesium Inorganic materials 0.000 description 1
- 239000011777 magnesium Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
- G06F21/60—Protecting data
- G06F21/602—Providing cryptographic facilities or services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- Bioethics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,1)建立含偏折裂纹的紧凑拉伸试样二维平面应力模型,在紧凑拉伸试样的两个销孔中心分别建立参考点,建立参考点和销孔内壁间的耦合约束;2)对模型进行网格划分,对裂纹尖端处的网格进行加密;3)通过隐式求解器,计算位移、应变及应力场,得到柔度和等效应力强度因子;4)建立裂纹长度ai/W与μi的关系,feq与ai/W的关系。本发明针对含偏折裂纹的紧凑拉伸试样,给出了不同偏折角度、裂纹扩展长度下的等效应力强度因子计算公式,可有效提升含偏折裂纹结构的裂纹监测及可靠性评价精度,减少因评价不合理造成的安全及经济风险。
Description
技术领域
本发明属于结构力学疲劳裂纹分析技术领域,具体涉及一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法。
背景技术
精准地监测裂纹长度和获取结构材料的疲劳裂纹扩展速率是基于损伤容限理论可靠性评价的前提和关键。然而,各向异性材料,如轧制的镁合金、钛合金及增材制造材料,有显著的材料取向性,在基于紧凑拉伸试样的纯I型加载下也会产生显著的裂纹偏折,从而无法准确监测裂纹扩展长度,也无法精确获取材料的裂纹扩展速率。
目前,根据金属材料疲劳裂纹扩展试验方法,紧凑拉伸试样是常用的试样类型。然而,针对紧凑拉伸试样,当前的基于柔度法的疲劳裂纹扩展长度预测公式及应力强度因子(裂纹扩展驱动力)计算公式均是假设裂纹沿着初始方向扩展。因此,使用直裂纹紧凑拉伸试样公式得到的裂纹扩展长度及疲劳裂纹扩展速率会带来较大的误差,造成无法对各向异性构件进行精确的可靠性评价和寿命管理,引入较大的安全及经济风险。为了解决当前存在的问题,常采用数值分析方法计算偏折裂纹下的柔度及应力强度因子。
然而,目前报道的结果中仅是考虑了几个裂纹偏折角度下的应力强度因子,无法适用于大范围、其他偏折角度下的结果,更没有建立基于柔度法的疲劳偏折裂纹长度预测公式。
综上所述,亟需采用数值分析方法,求解大范围(0°-90°)裂纹偏折角度、不同裂纹扩展长度下的紧凑拉伸试样柔度和应力强度因子,发展基于柔度法的裂纹长度监测公式和应力强度因子计算公式,最终实现对考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率的精准计算。
发明内容
本发明的目的在于克服现有技术的不足,提供一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,针对含偏折裂纹的紧凑拉伸试样,建立不同裂纹偏折角度、裂纹扩展长度下的柔度和应力强度因子的计算公式,实现对考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率的精准计算。
本发明解决其技术问题是通过以下技术方案实现的:
一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其特征在于:所述方法的步骤为:
S1、采用有限元分析软件Abaqus,建立含偏折裂纹的紧凑拉伸(CT)试样二维平面应力模型,输入材料的弹性模量E和泊松比ν,在CT试样的两个销孔中心分别建立参考点,通过运动耦合约束的方式,建立参考点和销孔内壁间的耦合约束,外部载荷P分别施加在上、下两个参考点上;
S2、采用带有缩减积分的8节点四边形单元对上述二维平面应力模型进行网格划分,对裂纹尖端处的网格进行加密并引入奇异单元,采用聚焦环式网格,最小网格尺寸为1μm;
S3、通过Abaqus自带的隐式求解器,计算位移、应变及应力场,得到不同裂纹偏折角度和裂纹扩展长度下的柔度和等效应力强度因子;
S4、对不同的裂纹偏折角度β(0°~90°),分别采用式(1)和式(2)对数据进行多项式拟合,建立裂纹长度与CT试样宽度的比值ai/W与归一化柔度μi的关系,归一化的等效应力强度因子feq(ai/W)与裂纹长度与试样宽度的比值ai/W间的关系:
其中:ai为裂纹长度,数值上等于(a0+Δaicos(β)),但对于偏折角度为90°的数值上等于(a0+Δai);
a0为初始裂纹长度;
Δai为裂纹偏折扩展长度;
β为裂纹偏折角度;
W为试样宽度;
μi为第i个归一化的柔度;
a-i为拟合参数;
feq(αi)=(2+αi)(1-αi)1.5·(a+bαi+cαi 2+dαi 3+eαi 4+fαi 5)(2)
其中:αi=ai/W是投影裂纹长度;
a-f为拟合参数。
而且,所述步骤S1中含偏折裂纹的紧凑拉伸(CT)试样的结构尺寸参考ASTME1820-17标准,具体为:
试样宽度W为50mm;
沿着裂纹方向上的整体试样宽度为1.25W=62.5mm;
垂直于裂纹方向上的整体试样宽度为1.2W=60mm;
裂纹在试样中间,初始裂纹长度a0为0.2W=10mm;
销孔沿着裂纹上下对称分布,销孔中心距离裂纹的长度为0.355W=17.75mm、距离试样开口边界为0.25W=12.5mm,销孔直径为0.25W=12.5mm;
裂纹偏折角度β分别为0°、7.5°、15°、22.5°、30°、37.5°、45°、52.5°、60°、67.5°、75°、82.5°和90°;
裂纹偏折扩展长度Δai在试样尺寸及模拟精度允许的情况下,尽量覆盖较大范围。
而且,所述步骤S3中柔度和等效应力强度因子的具体计算步骤为:
柔度是通过加载线上、缺口上下两个点处的变化位移与反应力的比值而得到,如式(3)所示:
其中:Δdi为变化的位移,ΔPi为变化的反应力;
等效应力强度因子由式(4)计算:
其中:KI和KII分别为I型和II型下应力强度因子,通过围道积分计算,可以直接从Abaqus的dat文件中输出。
而且,所述步骤S4中归一化柔度和归一化等效应力强度因子具体表达式如下:
其中:E为材料的弹性模量;
B为试样厚度(1mm);
C0为初始柔度,数值上等于8.60863;
其中:Keq为等效应力强度因子;
BN为净试样厚度,无侧槽时,BN=B;
Pi为施加的载荷。
本发明的优点和有益效果为:
1、本发明考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,针对含偏折裂纹的紧凑拉伸试样,提出基于柔度法的疲劳裂纹长度预测方法,给出了不同偏折角度、裂纹扩展长度下的等效应力强度因子计算公式,公式简单、易用。
2、本发明考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,疲劳裂纹长度预测公式及等效应力强度因子的求解公式,可有效提升含偏折裂纹结构的裂纹监测及可靠性评价精度,减少因评价不合理造成的安全及经济风险。
附图说明
图1为本发明的紧凑拉伸试样几何示意图;
图2为本发明试验试样在镁合金轧制板上的取样方式示意图;
图3为本发明有限元模型几何示意图;
图4为本发明紧凑拉伸试样有限元网格图;
图5为本发明有限元计算的裂纹长度与柔度之间的关系图;
图6为本发明有限元计算的等效应力强度因子与裂纹长度之间的关系图;
图7为本发明试验和预测的疲劳裂纹扩展长度曲线图;
图8为本发明基于等效应力强度因子建立的疲劳裂纹扩展速率关系图。
具体实施方式
下面通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其创新之处在于:所述方法的步骤为:
S1、采用有限元分析软件Abaqus,建立含偏折裂纹的紧凑拉伸(CT)试样二维平面应力模型,输入材料的弹性模量E和泊松比ν,在CT试样的两个销孔中心分别建立参考点,通过运动耦合约束的方式,建立参考点和销孔内壁间的耦合约束,外部载荷P分别施加在上、下两个参考点上;
含偏折裂纹的紧凑拉伸(CT)试样的结构尺寸参考ASTME1820-17标准,具体为:
试样宽度W为50mm;
沿着裂纹方向上的整体试样宽度为1.25W=62.5mm;
垂直于裂纹方向上的整体试样宽度为1.2W=60mm;
裂纹在试样中间,初始裂纹长度a0为0.2W=10mm;
销孔沿着裂纹上下对称分布,销孔中心距离裂纹的长度为0.355W=17.75mm、距离试样开口边界为0.25W=12.5mm,销孔直径为0.25W=12.5mm;
裂纹偏折角度β分别为0°、7.5°、15°、22.5°、30°、37.5°、45°、52.5°、60°、67.5°、75°、82.5°和90°;
裂纹偏折扩展长度Δai在试样尺寸及模拟精度允许的情况下,尽量覆盖较大范围;
S2、采用带有缩减积分的8节点四边形单元对上述二维平面应力模型进行网格划分,对裂纹尖端处的网格进行加密并引入奇异单元,采用聚焦环式网格,最小网格尺寸为1μm;
S3、通过Abaqus自带的隐式求解器,计算位移、应变及应力场,得到不同裂纹偏折角度和裂纹扩展长度下的柔度和等效应力强度因子;
柔度和等效应力强度因子的具体计算步骤为:
柔度Ci是通过加载线上、缺口上下两个点处的变化位移与反应力的比值而得到,如式(3)所示:
其中:Δdi为变化的位移,ΔPi为变化的反应力。
等效应力强度因子由式(4)计算:
其中:KI和KII分别为I型和II型下应力强度因子,通过围道积分计算,可以直接从Abaqus的dat文件中输出。
S4、对不同的裂纹偏折角度β(0°~90°),分别采用式(1)和式(2)对数据进行多项式拟合,建立裂纹长度与CT试样宽度的比值ai/W与归一化柔度μi的关系,归一化的等效应力强度因子feq(ai/W)与裂纹长度与试样宽度的比值ai/W间的关系:
其中:ai为裂纹长度,数值上等于(a0+Δaicos(β)),对于偏折角度为90°的数值上等于(a0+Δai);
a0为初始裂纹长度;
Δai为裂纹偏折扩展长度;
β为裂纹偏折角度;
W为试样宽度;
μi为第i个归一化的柔度;
a-i为拟合参数,列于表1中;
表1基于柔度的疲劳裂纹长度预测公式的相关系数
feq(αi)=(2+αi)(1-αi)1.5·(a+bαi+cαi 2+dαi 3+eαi 4+fαi 5)(2)其中:αi=ai/W是投影裂纹长度;
a-f为拟合参数,列于表2中。
表2归一化的等效应力强度因子计算公式的相关系数
归一化柔度和归一化等效应力强度因子具体表达式如下:
其中:E为材料的弹性模量;
B为试样厚度(1mm);
C0为初始柔度,数值上等于8.60863;
其中:Keq为等效应力强度因子;
BN为净试样厚度,无侧槽时,BN=B;
Pi为施加的载荷。
本发明以某镁合金制造公司生产的轧制ZK60镁锌合金为例,使用本发明提出的方法,得到其准确的疲劳裂纹扩展长度及疲劳裂纹扩展速率。通过单轴拉伸试验,测得ZK60镁合金的基本力学性能:横向方向(TD)的弹性模量E11为45256MPa,法线方向(ND)的弹性模量E22为44527MPa,泊松比ν12为0.35,剪切模量G12为17000MPa。
下面,将按照本发明的方法,对ZK60镁合金疲劳裂纹扩展长度的准确预测并得到其疲劳裂纹扩展速率。
(1)试样设计及加工
采用紧凑拉伸(CT)试样,其几何示意图如图1所示。试样宽度W为30mm,初始裂纹长度a0为6mm,试样厚度B为3.4mm。CT试样在ZK60镁合金轧制板TD-ND面上切取(如图2),设计CT试样的裂纹方向与TD方向不同的夹角,依次为0°、45°和67.5°,分别命名为试样0°-90°、45°-45°和67.5°-22.5°,其中第一个角度是裂纹方向与TD方向的角度,第二个角度为裂纹方向与ND方向的角度。CT试样加工完毕后,在试验之前对表面进行打磨抛光,以便在高倍摄像机下观察裂纹扩展。
(2)疲劳裂纹扩展实验过程
所有疲劳裂纹扩展试验均在MTS370.02拉扭疲劳试验机上进行,采用Epsilon3541COD引伸计测量裂纹张开位移。首先是疲劳预制裂纹,采用正弦波加载,应力比为0.1,加载频率为10Hz,通过高倍相机实时监测裂纹扩展长度,当疲劳预制裂纹长度达1mm时,正式开始疲劳裂纹扩展试验。在正式的疲劳裂纹扩展实验过程中,采用快循环与慢循环交替加载,快循环加载频率为10Hz,慢循环加载频率为0.05Hz,在慢循环加载过程中相机每间隔5s拍一张试样照片,用于监测裂纹扩展长度,在试样完全断裂之后,停止试验。
(3)试验数据处理
实验过程中会收集两类数据,一类是相机所拍摄的试样表面的照片,一类是MTS试验机采集的文本数据,包括时间、载荷、循环圈数等。首先,采用Matlab软件从文本数据中提取所有慢循环阶段的循环圈数、载荷、位移和时间,再根据时间从相机所拍摄试样照片中找到对应的照片。随后,利用图片处理软件,如ImageJ对照片中的裂纹偏折角及裂纹扩展长度进行测量。通过位移与载荷的比值得到柔度,最终,得到每个时刻对应下的循环圈数、柔度、疲劳裂纹偏折角度及裂纹扩展长度。
(4)有限元建立柔度和等效应力强度因子与裂纹扩展长度间的关系
使用有限元软件Abaqus,建立含偏折裂纹的紧凑拉伸试样的二维平面应力模型(图3所示),其中试样宽度W为50mm,初始裂纹长度a0为0.2W=10mm,销孔直径R为0.25W=12.5mm,裂纹偏折角度β分别为0°、7.5°、15°、22.5°、30°、37.5°、45°、52.5°、60°、67.5°、75°、82.5°和90°,裂纹偏折扩展长度Δa尽量大,使所有试样的a/W在0.9左右。注,对于偏折角度β为0°-82.5°时,a/W=(a0+Δai cos(β))/W,对于偏折角度β为90°时,a/W=(a0+Δai)/W。输入镁合金弹性模量45256MPa,泊松比0.35。采用带有缩减积分的8节点四边形单元对模型进行网格划分,对裂纹尖端处的网格进行加密和引入奇异单元,采用聚焦环式网格,最小网格尺寸为1μm,有限元网格模型如图4所示。计算不同裂纹偏折角度和裂纹扩展长度下的柔度和等效应力强度因子,结果分别如图5和图6所示,拟合的经验公式如式(1)和式(2)。
通过试验,得到试样0°-90°、45°-45°和67.5°-22.5°对应的裂纹偏折角度β分别为0°、36.95°和62.2°。将试验得到的柔度代入式(1),通过插值,就可以得到偏折角度为0°和62.2°对应的裂纹长度,如图7所示。
通过与相机拍摄的对比,发现当前发展的基于柔度的疲劳裂纹扩展长度预测方法是有效的。另一方面,将得到的瞬时的裂纹长度代入式(2)和式(6),通过插值,可以计算试样45°-45°对应的等效应力强度因子Keq,借助Matlab软件拟合疲劳裂纹扩展速率da/dN关系,最终可以得到da/dN-Keq关系曲线,如图8所示。
尽管为说明目的公开了本发明的实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此,本发明的范围不局限于实施例和附图所公开的内容。
Claims (4)
1.一种考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其特征在于:所述方法的步骤为:
S1、采用有限元分析软件Abaqus,建立含偏折裂纹的CT试样二维平面应力模型,输入材料的弹性模量E和泊松比ν,在CT试样的两个销孔中心分别建立参考点,通过运动耦合约束的方式,建立参考点和销孔内壁间的耦合约束,外部载荷P分别施加在上、下两个参考点上;
S2、采用带有缩减积分的8节点四边形单元对上述二维平面应力模型进行网格划分,对裂纹尖端处的网格进行加密并引入奇异单元,采用聚焦环式网格,最小网格尺寸为1μm;
S3、通过Abaqus自带的隐式求解器,计算位移、应变及应力场,得到不同裂纹偏折角度和裂纹扩展长度下的柔度和等效应力强度因子;
S4、对不同的裂纹偏折角度β(0°~90°),分别采用式(1)和式(2)对数据进行多项式拟合,建立裂纹长度与CT试样宽度的比值ai/W与归一化柔度μi的关系,归一化的等效应力强度因子feq(ai/W)与裂纹长度与试样宽度的比值ai/W间的关系:
ai/W=a+b/μi+c/μi 2+d/μi 3+e/μi 4+f/μi 5+g/μi 6+h/μi 7+i/μi 8 (1)
其中:ai为裂纹长度,数值上等于(a0+Δaicos(β)),但对偏折角度为90°的数值上等于(a0+Δai);
a0为初始裂纹长度;
Δai为裂纹偏折扩展长度;
β为裂纹偏折角度;
W为试样宽度;
μi为第i个归一化的柔度;
a、b、c、d、e、f、g、h和i为拟合参数;
feq(αi)=(2+αi)/(1-αi)1.5·(a′+b′αi+c′αi 2+d′αi 3+e′αi 4+f′αi 5) (2)
其中:αi=ai/W是投影裂纹长度;
a′、b′、c′、d′、e′和f′为拟合参数;
S5、对于不同裂纹偏折角度,通过式(1)得到实时的裂纹长度a,通过式(2)计算等效应力强度因子Keq,借助Matlab软件拟合疲劳裂纹扩展长度与循环圈数N的比值da/dN,最终得到da/dN-Keq关系曲线即为疲劳裂纹扩展速率关系。
2.根据权利要求1所述考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其特征在于:所述步骤S1中含偏折裂纹的CT试样的结构尺寸参考ASTME1820-17标准,具体为:
试样宽度W为50mm;
沿着裂纹方向上的整体试样宽度为1.25W=62.5mm;
垂直于裂纹方向上的整体试样宽度为1.2W=60mm;
裂纹在试样中间,初始裂纹长度a0为0.2W=10mm;
销孔沿着裂纹上下对称分布,销孔中心距离裂纹的长度为0.355W=17.75mm、距离试样开口边界为0.25W=12.5mm,销孔直径为0.25W=12.5mm;
裂纹偏折角度β分别为0°、7.5°、15°、22.5°、30°、37.5°、45°、52.5°、60°、67.5°、75°、82.5°和90°。
3.根据权利要求1所述考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其特征在于:所述步骤S3中柔度和等效应力强度因子的具体计算步骤为:
柔度Ci是通过加载线上、缺口上下两个点处的变化位移与反应力的比值而得到,如式(3)所示:
其中:Δdi为变化的位移,ΔPi为变化的反应力;
等效应力强度因子由式(4)计算:
其中:KI和KII分别为I型和II型下应力强度因子,通过围道积分计算,直接从Abaqus的dat文件中输出。
4.根据权利要求3所述考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法,其特征在于:所述步骤S4中归一化柔度和归一化等效应力强度因子具体表达式如下:
其中:E为材料的弹性模量;
B为试样厚度;
C0为初始柔度,数值上等于8.60863;
其中:Keq为等效应力强度因子;
BN为净试样厚度,无侧槽时,BN=B;
Pi为施加的载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310609963.0A CN116738780B (zh) | 2023-05-29 | 2023-05-29 | 考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310609963.0A CN116738780B (zh) | 2023-05-29 | 2023-05-29 | 考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116738780A CN116738780A (zh) | 2023-09-12 |
CN116738780B true CN116738780B (zh) | 2024-01-26 |
Family
ID=87912483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310609963.0A Active CN116738780B (zh) | 2023-05-29 | 2023-05-29 | 考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116738780B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117371271B (zh) * | 2023-09-22 | 2024-06-18 | 天津大学 | 基于销轴式单边缺口拉伸试样的各向异性材料裂纹长度及断裂性能测试方法 |
CN117371272B (zh) * | 2023-09-22 | 2024-04-19 | 天津大学 | 适用于不同各向异性材料及尺寸的夹持式单边缺口拉伸试样裂纹长度及断裂性能计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102564844A (zh) * | 2011-12-30 | 2012-07-11 | 中国船舶重工集团公司第七二五研究所 | 一种采用双引伸计测量紧凑拉伸试样断裂参量的方法 |
CN110672441A (zh) * | 2019-10-21 | 2020-01-10 | 江苏理工学院 | 一种疲劳裂纹扩展过程中裂尖状态的表征方法 |
CN111665132A (zh) * | 2020-06-17 | 2020-09-15 | 天津大学 | 一种开口角隅疲劳裂纹扩展测量方法 |
CN114323940A (zh) * | 2021-12-21 | 2022-04-12 | 中国科学院金属研究所 | 梯型波加载条件下的紧凑拉伸试样疲劳裂纹扩展速率的测试方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104237007B (zh) * | 2013-06-07 | 2016-09-14 | 合肥通用机械研究院 | 一种测试复合型疲劳裂纹扩展速率的装置 |
-
2023
- 2023-05-29 CN CN202310609963.0A patent/CN116738780B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102564844A (zh) * | 2011-12-30 | 2012-07-11 | 中国船舶重工集团公司第七二五研究所 | 一种采用双引伸计测量紧凑拉伸试样断裂参量的方法 |
CN110672441A (zh) * | 2019-10-21 | 2020-01-10 | 江苏理工学院 | 一种疲劳裂纹扩展过程中裂尖状态的表征方法 |
CN111665132A (zh) * | 2020-06-17 | 2020-09-15 | 天津大学 | 一种开口角隅疲劳裂纹扩展测量方法 |
CN114323940A (zh) * | 2021-12-21 | 2022-04-12 | 中国科学院金属研究所 | 梯型波加载条件下的紧凑拉伸试样疲劳裂纹扩展速率的测试方法 |
Non-Patent Citations (2)
Title |
---|
CRO试样的疲劳裂纹扩展行为试验方法研究;赵兴华;蔡力勋;包陈;;《工程力学》(第11期);第25-33页 * |
刘争 ; 孙少南 ; 王昕 ; 陈旭.面内和面外约束效应对3D紧凑拉伸试样的影响.《压力容器先进技术—第九届全国压力容器学术会议论文集》.2017,第1-7页. * |
Also Published As
Publication number | Publication date |
---|---|
CN116738780A (zh) | 2023-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116738780B (zh) | 考虑裂纹偏折的紧凑拉伸试样疲劳裂纹扩展长度及速率计算方法 | |
Peirs et al. | Novel technique for static and dynamic shear testing of Ti6Al4V sheet | |
KR101152323B1 (ko) | 계장화 압입시험법을 이용한 비등방향 잔류응력의 평가방법 | |
CN106644711B (zh) | 一种延性材料单轴本构关系测试方法 | |
CN111665132A (zh) | 一种开口角隅疲劳裂纹扩展测量方法 | |
CN117371271B (zh) | 基于销轴式单边缺口拉伸试样的各向异性材料裂纹长度及断裂性能测试方法 | |
Weeks et al. | Comparison of J-integral from single specimen SE (T) tests on API-5L X100 line pipe steel | |
CN108548720B (zh) | I型裂纹弹塑性理论公式获取延性材料j阻力曲线的方法 | |
Lin et al. | Biaxial tensile behavior of Ti-6Al-4V under proportional loading | |
Esmaeili et al. | Investigation of T-stress and tensile strength effect on crack tip conditions and crack initiation angle in off-axis laminate composite | |
Pantousa et al. | Experimental and numerical investigation of the ductile fracture of structural steel at elevated temperatures | |
Li et al. | Identification of material parameters from punch stretch test | |
Fa-Kun et al. | Numerical investigation of frictional effect on measuring accuracy of different small specimen creep tests | |
Wolf et al. | Crack growth behaviour in biaxial stress fields: Calculation of K-factors for cruciform specimens | |
Han et al. | Relation between ASTM E606 specimen geometry and misalignment in strain-controlled fatigue testing | |
Ranjan et al. | A fracture-mechanics-based methodology for fatigue life prediction of single crystal nickel-based superalloys | |
He et al. | A wide range stress intensity factor solution for an eccentrically cracked middle tension specimen with clamped ends | |
CN113358678B (zh) | α钛变形过程介观应力和织构的半定量预测及可视化方法 | |
Tretyakova et al. | Experimental study of the influence of strain-stress state on the jerky flow in metals and alloys | |
CN117371272B (zh) | 适用于不同各向异性材料及尺寸的夹持式单边缺口拉伸试样裂纹长度及断裂性能计算方法 | |
Mei et al. | Characterization of mechanical property degradation of ion-irradiated materials | |
Qi et al. | Compliance-based testing method for fatigue crack propagation rates of mixed-mode I–II cracks | |
CN108918298B (zh) | 三维裂尖约束效应及三维裂尖等效应力场的表征方法 | |
Jin et al. | The Yield Behavior of Oxygen-Free Copper Under Biaxial Quasi-Static and Dynamic Loadings | |
Arnold et al. | Interpretation of Experimental Tensile Test Results and Its Implication on Damage Modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |