CN113609731B - 一种预判复杂应力状态下镁合金变形临界开裂的方法 - Google Patents
一种预判复杂应力状态下镁合金变形临界开裂的方法 Download PDFInfo
- Publication number
- CN113609731B CN113609731B CN202110877719.3A CN202110877719A CN113609731B CN 113609731 B CN113609731 B CN 113609731B CN 202110877719 A CN202110877719 A CN 202110877719A CN 113609731 B CN113609731 B CN 113609731B
- Authority
- CN
- China
- Prior art keywords
- magnesium alloy
- strain
- deformation
- critical cracking
- cracking
- 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
Images
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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- 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)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (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)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明公开了一种预判复杂应力状态下镁合金变形临界开裂的方法,涉及镁合金塑性成形领域。具体方法包括:1)实验获取镁合金的临界开裂应变εf;2)建立临界开裂应变εf与应力三轴度σ*、相对应变速率相对温度T*之间的数学模型;3)比较临界开裂应变εf和模拟应变ε;4)模拟确定镁合金的临界开裂损伤Cf;5)建立临界开裂损伤Cf与变形温度T、应变速率之间的数学模型;6)确定复杂应力状态下镁合金变形的临界开裂判据。本发明制定一种能够准确预判复杂应力状态下镁合金变形临界开裂的判据,可为冲压、旋压等处于复杂应力状态下的深加工工艺优化提供指导,能有效提高镁合金深加工产品成形质量和生产效率,降低生产成本,节约资源。
Description
技术领域
本发明属于镁合金塑性成形领域,具体涉及一种预判复杂应力状态下镁合金变形临界开裂的方法。
背景技术
航空航天、轨道交通、新能源汽车等领域装备构件的轻量化制造需求越发紧急和迫切。镁合金作为最轻的工程金属材料,在卫星适配器锥桶件、厢体盖板、车身覆盖件等异型壁板构件减重上应用前景广阔。此类构件通常具有大尺寸、变曲率和异形深腔的形状特征且服役在承载/次承载状态,采用镁合金板材深加工(如冲压、旋压等)成形是实现其快速低成本一体化制备的重要途径。但镁合金为密排六方晶体结构,复杂应力状态下的塑性变形协调能力差,稳定变形区间窄,导致塑性变形受工艺条件和变形状态强烈影响,控制不当极易产生开裂缺陷,成为限制其发展的一个重要因素。
通过模拟手段无法观察到材料的宏观开裂现象,同时,仅通过实验手段获取镁合金临界开裂状态的方法,不但会消耗大量原材料使得实验成本高,而且费时费力。此外,与轴向拉伸、压缩不同,材料在冲压、旋压等深加工成形时处于复杂应力状态,传统基于单轴加载数据制定开裂判据的方法,不适用于复杂应力状态下板材深加工开裂预判。本发明结合了物理实验和有限元模拟,考虑三轴应力状态、应变速率和变形温度的综合影响,建立了三者与临界开裂应变之间的数学关系,进而制定一种能够准确预判复杂应力状态下镁合金变形开裂的判据,据此实现工艺参数匹配关系预调整,对规避开裂缺陷发生,提高板材深加工成形质量具有重要的意义。
发明内容
针对镁合金板材在复杂应力状态下塑性成形时极易产生裂纹缺陷的问题,本发明提供了一种预判复杂应力状态下镁合金变形临界开裂的方法。该方法结合了轴向拉伸实验和有限元模拟,通过建立适合描述473K~723K温度、0.01s-1~20s-1应变速率、0.33~1.67应力三轴度范围内的临界开裂应变数学模型和用于计算模拟过程中临界损伤的数学模型,进而建立包含累积损伤因子的临界开裂判据,据此对复杂应力状态下热塑性加工过程中镁合金的开裂趋势做出准确预判,以此优化工艺过程。
本发明中临界开裂应变εf数学模型的建立方法为。
步骤1:确定实验状态下镁合金的临界开裂应变εf:通过对具有不同缺口半径的镁合金矩形截面试样进行轴向拉伸,缺口半径R≥0.5mm,变形温度T为 473K~723K、拉伸速度v为1.5~30000mm/min,获得应力σ和应变ε数据,在应力-应变曲线上选取临近应力卸载的加工硬化率骤降点,骤降点处的应变,即为临界开裂应变εf。
步骤2:对步骤1中的镁合金矩形截面试样进行相同实验条件下的拉伸变形有限元模拟。
(1)建立临界开裂应变εf对应力三轴度σ*的依赖关系:跟踪变形区的应变ε分布状态,确定模拟中应变首先达到临界开裂应变εf时的位置,即起裂处,提取该处的应力状态,通过计算该处的应力三轴度σ*,其中σ1、σ2、σ3为起裂处的三个主应力,对εf-σ*数据进行非线性拟合得到/>式中,χ1是与镁合金材料牌号相关的常数,取值与镁质量分数呈负相关,范围是-13.8554~-11.7532。
(2)建立临界开裂应变εf对相对应变速率的依赖关系:通过/>计算得到相对应变速率/>式中,L为拉伸试样标距,进而对/>数据进行非线性拟合得到式中,χ2是与温度有关的常数,T<573K 时,χ2的取值范围是3.20~3.22,573K≤T≤723K时,χ2取值范围是2.05~2.28,χ3是与板材表面粗糙度有关的常数,取值与粗糙度呈负相关,范围是-0.12~0.05。
(3)建立临界开裂应变εf对相对温度T*的依赖关系:通过计算得到相对温度T*,式中,Tm为镁合金熔点,进而对εf-T*数据进行线性拟合得到εf=(χ4+χ5T*),χ4是与温度有关的常数,温度T<573K时,χ4的取值范围是1.480~1.533,温度573K≤T≤723K时,χ4的取值范围是-0.501~-0.473;χ5是与镁合金材料牌号相关的常数,取值与镁质量分数呈正相关,范围是 -1.90~2.61。
步骤三:通过上述数学模型,计算镁合金的临界开裂应变εf,比较临界开裂应变εf和有限元模拟变形过程中提取的应变ε,当ε≥εf时,说明发生开裂,当ε<εf时,进一步建立临界开裂判据。
本发明中临界开裂判据的建立方法为。
步骤四:获取有限元模拟中镁合金开裂处的临界损伤值Cf:采用断裂损伤准则,式中,C为损伤值,n为镁合金的加工硬化指数,对相同实验条件下试样的拉伸变形进行有限元模拟,跟踪试样在变形过程中的应变状态和损伤分布,确定模拟中应变首先达到临界开裂应变εf的点的位置,获得该点损伤值,即为临界开裂损伤Cf。
步骤五:建立临界开裂损伤Cf与变形温度T和应变速率之间的数学模型:温度和应变速率对镁合金变形临界开裂损伤Cf的综合影响,通过引入 Zener-Hollomon因子,来表示,式中,Q为镁合金变形激活能,应变速率/>由/>计算得到,对Cf和数据进行线性拟合得到/>式中,M和N是与材料相关的常数。
通过上述过程,发明一种预判复杂应力状态下镁合金变形临界开裂的方法,将上述临界开裂判据编译到有限元模拟软件中,能够准确预判镁合金在复杂应力状态下的实际变形过程中是否发生开裂。
本发明的优点在于:结合物理实验和有限元模拟,建立了镁合金累积损伤与应力状态、变形温度、应变速率之间的关系,来定量预判在复杂应力状态下镁合金变形过程中的临界开裂状态,具有参数易于获得、方法简单易行等优点,可为镁合金的冲压、旋压等处于复杂应力状态下的深加工工艺优化提供指导,据此实现工艺参数匹配关系预调整,对规避开裂缺陷发生、提高镁合金深加工产品成形质量和生产效率、降低生产成本和节约资源具有重要意义。
附图说明
图1是预判复杂应力状态下镁合金变形临界开裂方法的实现流程。
图2是拉伸实验的试样示意图。
具体实施方式
本发明提供一种预判复杂应力状态下镁合金变形临界开裂的方法。
以厚度宽度是7mm*300mm的铸轧态AZ31镁合金板材为例,其化学成分组织如表1所示。
表1实施例铸态AZ31镁合金的化学成分(wt%)
采用万能拉伸机、Gleeble 3800热模拟试验机对具有不同缺口的镁合金矩形截面试样进行轴向拉伸,变形温度T=573K~673K、拉伸速度v =1.5~1500mm/min;拉伸试样为光滑试样和缺口半径R为1mm、4mm的缺口试样,试样采用相同的标距段L=25mm。
热拉伸实验中试样的升温速率为10K/s,达到实验温度后保温3min,同时在试样侧面焊上热电偶用于实时温度监测和采集,拉伸直至试样断裂。
本实例中临界开裂应变εf数学模型的建立过程如下。
表2试样在起裂时对应的临界开裂应变与应力三轴度值
(2)对相同实验条件下,相同的铸轧态AZ31镁合金试样进行拉伸变形的有限元模拟,跟踪变形区的应变分布状态,确定模拟中应变首先达到临界开裂应变εf时的位置,即起裂处,提取该处的应力状态,通过计算该处的应力三轴度σ*,如表2,对εf-σ*数据进行非线性拟合,拟合相关系数为0.94503;通过计算得到相对应变速率/>对/>数据进行非线性拟合,拟合相关系数为 0.97382;通过/>计算得到相对温度T*,如表3,对εf-T*数据进行线性拟合,拟合相关系数为0.96473。
表3不同温度时刻的断裂应变值
(3)通过上述数学模型,计算镁合金的临界开裂应变εf,经比较,本实例中铸轧态AZ31镁合金的临界开裂应变εf大于有限元模拟变形过程中提取的应变ε,即ε<εf。
进一步建立临界开裂判据,建立方法为。
(4)分别对变形温度T=573K,应变速率缺口半径R=4mm 的缺口试样和变形温度T=673K,应变速率/>的光滑试样进行损伤模拟分析,据轴向拉伸实验结果得到本实例中铸轧态AZ31镁合金的加工硬化指数n =0.3,模拟采用断裂损伤准则,跟踪试样在变形过程中的应变状态和损伤分布。在试样标距段中心沿垂直于拉伸方向的宽向等距离选取跟踪点,光滑试样选取6个点,缺口试样选取4 个点,利用有限元软件的点跟踪模块记录拉伸过程中各跟踪点损伤值累积,确定模拟中应变首先达到临界开裂应变εf的点位置,获得该点损伤值,即为临界开裂损伤Cf。
(5)由应力应变曲线计算得到材料的变形激活能Q=236.543KJ/mol,通过计算不同温度和应变速率下的Zener-Hollomon因子,该参数表征了温度和应变速率对镁合金变形临界开裂损伤Cf的综合影响,对Cf和/>的数据进行线性拟合得到/>拟合相关系数为0.96378。
计算得到D=-0.37<0,说明在变形温度T=573K~673K,拉伸速度v =1.5~1500mm/min下,本实例的铸轧态AZ31镁合金在变形过程中不发生开裂。
Claims (2)
1.一种预判复杂应力状态下镁合金变形临界开裂的方法,其特征在于,该方法结合了物理实验和有限元模拟,包括以下步骤:
步骤1:对具有不同缺口半径的镁合金矩形截面试样进行轴向拉伸,缺口半径R≥0.5mm,变形温度T为473~723K、拉伸速度v为1.5~30000mm/min,获取应力-应变曲线上临近应力卸载的加工硬化率骤降点,其中σ为拉伸应力,ε为拉伸应变,确定骤降点应变为临界开裂应变εf;
步骤2:对步骤1中的镁合金矩形截面试样进行相同实验条件下的拉伸变形有限元模拟,跟踪变形区的应变ε分布状态,确定模拟中应变首先达到临界开裂应变εf时的位置,即起裂处,提取该处应力状态,通过计算该处的应力三轴度σ*,其中σ1、σ2、σ3为起裂处的三个主应力,建立临界开裂应变εf与应力三轴度σ*、相对应变速率/>相对温度T*之间的数学模型,其中相对应变速率/>通过/>计算得到,相对温度T*通过/>计算得到,式中,L为拉伸试样标距,Tm为镁合金熔点;
步骤3:比较步骤2中模型计算得到的镁合金临界开裂应变εf和有限元模拟变形过程中提取的应变ε,当ε≥εf时,说明发生开裂,当ε<εf时,进一步建立临界开裂判据;
步骤4:采用断裂损伤准则,式中,C为损伤值,n为镁合金的加工硬化指数,对相同实验条件下试样的拉伸变形进行有限元模拟,跟踪试样在变形过程中的应变ε分布状态和损伤C分布状态,确定模拟中应变首先达到临界开裂应变εf的点位置,获得该点损伤值,即为临界开裂损伤Cf;
式中,χ1是与镁合金材料牌号相关的常数,取值与镁质量分数呈负相关,范围是-13.8554~-11.7532;χ2是与温度有关的常数,温度T<573K时,χ2的取值范围是3.20~3.22,温度573K≤T≤723K时,χ2的取值范围是2.05~2.28;χ3是与板材表面粗糙度有关的常数,取值与粗糙度呈负相关,范围是-0.12~0.05;χ4是与温度有关的常数,温度T<573K时,χ4的取值范围是1.480~1.533,温度573K≤T≤723K时,χ4的取值范围是-0.501~-0.473;χ5是与镁合金材料牌号相关的常数,取值与镁质量分数呈正相关,范围是-1.90~2.61。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110877719.3A CN113609731B (zh) | 2021-08-01 | 2021-08-01 | 一种预判复杂应力状态下镁合金变形临界开裂的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110877719.3A CN113609731B (zh) | 2021-08-01 | 2021-08-01 | 一种预判复杂应力状态下镁合金变形临界开裂的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113609731A CN113609731A (zh) | 2021-11-05 |
CN113609731B true CN113609731B (zh) | 2023-06-27 |
Family
ID=78338948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110877719.3A Active CN113609731B (zh) | 2021-08-01 | 2021-08-01 | 一种预判复杂应力状态下镁合金变形临界开裂的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113609731B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504185A (zh) * | 2014-12-10 | 2015-04-08 | 西北工业大学 | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 |
RU2634800C1 (ru) * | 2016-11-22 | 2017-11-03 | Акционерное общество "Чепецкий механический завод" | Способ определения порога напряжений коррозионного растрескивания стали или сплава при постоянной деформации |
CN111241721A (zh) * | 2019-11-12 | 2020-06-05 | 太原科技大学 | 一种镁合金薄板轧制边部开裂准则及深度预判方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3558628B2 (ja) * | 2002-06-05 | 2004-08-25 | 住友電工スチールワイヤー株式会社 | マグネシウム合金板およびその製造方法 |
-
2021
- 2021-08-01 CN CN202110877719.3A patent/CN113609731B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504185A (zh) * | 2014-12-10 | 2015-04-08 | 西北工业大学 | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 |
RU2634800C1 (ru) * | 2016-11-22 | 2017-11-03 | Акционерное общество "Чепецкий механический завод" | Способ определения порога напряжений коррозионного растрескивания стали или сплава при постоянной деформации |
CN111241721A (zh) * | 2019-11-12 | 2020-06-05 | 太原科技大学 | 一种镁合金薄板轧制边部开裂准则及深度预判方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113609731A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Constitutive analysis of homogenized 7005 aluminum alloy at evaluated temperature for extrusion process | |
Srivastava et al. | Computer modeling and prediction of thermal fatigue cracking in die-casting tooling | |
Hockauf et al. | Mechanical properties and corrosion behaviour of ultrafine-grained AA6082 produced by equal-channel angular pressing | |
Gao et al. | Investigation on mechanical behavior and plastic damage of AA7075 aluminum alloy by thermal small punch test: Experimental trials, numerical analysis | |
Wang et al. | Hot deformation of austenite and prediction of microstructure evolution of cross-wedge rolling | |
RU2618016C2 (ru) | Титановый сплав для крепежа | |
CN107937842B (zh) | 一种Al-Zn-Mg-Cu铝合金的热处理方法 | |
Vijayanandh et al. | Comparative approaches for fatigue life estimation of aluminium alloy for aerospace applications | |
Yang et al. | Enhancing hoop strength of titanium alloy tube by cross spinning | |
CN105677949A (zh) | Tc6钛合金锻件微观组织参数的预测方法 | |
Wang et al. | Deformation efficiency, homogeneity, and electrical resistivity of pure copper processed by constrained groove pressing | |
Wang et al. | Tensile properties and a modified s-Johnson-Cook model for constitutive relationship of AA7075 sheets at cryogenic temperatures | |
CN113609731B (zh) | 一种预判复杂应力状态下镁合金变形临界开裂的方法 | |
Chen et al. | A study on static properties of high-performance steel after corrosion damage | |
Wang et al. | Constitutive modeling for high temperature flow behavior of a high-strength manganese brass | |
Joo et al. | Finite element and experimental analysis of closure and contact bonding of pores during hot rolling of steel | |
Zhang et al. | Effect of different temperatures on the metallographic structure and tensile property of 2024‐T4 alloy in integral heating single point incremental forming | |
WANG et al. | Deformation behavior of TC1 titanium alloy sheet under double-sided pressure | |
Scarabello et al. | FE modelling of large ingot hot forging | |
CN113486458A (zh) | 基于deform模拟预测各向异性对钢折弯开裂影响的方法 | |
Gerin et al. | Effect of cold forming on the high cycle fatigue behaviour of a 27MnCr5 steel | |
CN111366606A (zh) | 一种基于膨胀曲线的热处理变形预测方法 | |
Zhuang et al. | Experimental investigation of the effect of the material damage induced in sheet metal forming process on the service performance of 22MnB5 steel | |
Remshev et al. | Investigation of the influence of modification of the working tool for the manufacture of niobium case products | |
Graf et al. | Aluminium-Matrix-Composites (AMC) for Hot Forged Components |
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 |