CN109558692B - 预测微粒冲击金属件残余应力和马氏体相变的有限元方法 - Google Patents
预测微粒冲击金属件残余应力和马氏体相变的有限元方法 Download PDFInfo
- Publication number
- CN109558692B CN109558692B CN201811591935.6A CN201811591935A CN109558692B CN 109558692 B CN109558692 B CN 109558692B CN 201811591935 A CN201811591935 A CN 201811591935A CN 109558692 B CN109558692 B CN 109558692B
- Authority
- CN
- China
- Prior art keywords
- stress
- increment
- abaqus
- strain
- strengthening
- 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
- 229910000734 martensite Inorganic materials 0.000 title claims abstract description 57
- 239000002184 metal Substances 0.000 title claims abstract description 54
- 239000002245 particle Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000000463 material Substances 0.000 claims abstract description 26
- 230000008569 process Effects 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims description 32
- 238000005728 strengthening Methods 0.000 claims description 27
- 230000008859 change Effects 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 239000007769 metal material Substances 0.000 claims description 5
- 238000005482 strain hardening Methods 0.000 claims description 4
- 230000005489 elastic deformation Effects 0.000 claims description 3
- 239000012798 spherical particle Substances 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000003116 impacting effect Effects 0.000 claims description 2
- 229910000831 Steel Inorganic materials 0.000 claims 1
- 239000010959 steel Substances 0.000 claims 1
- 230000003014 reinforcing effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000004566 building material Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011236 particulate material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 239000007787 solid Substances 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]
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
本发明公开了一种预测微粒高速冲击下金属构件残余应力和马氏体相变的有限元方法,包括以下步骤:(1)建立微粒冲击金属构件的有限元模型;(2)建立综合考虑了塑性应变强化和马氏体相变强化的率相关弹塑性本构模型;(3)编写ABAQUS‑VUMAT用户动态材料子程序,运用应力补偿更新算法实现提出的率相关弹塑性本构模型;(4)将子程序嵌入到ABAQUS中,对微粒流冲击金属构件的过程进行计算,进一步获得金属构件内部的残余应力以及马氏体相的百分比含量。本发明可以用于任意数量的微粒对任意结构形式的金属构件进行冲击的残余应力和马氏体相变预测。
Description
技术领域
本发明涉及金属构件残余应力和组织相变预测,旨在提供预测微粒高速冲击下金属构件残余应力和马氏体相变的有限元方法。
背景技术
微粒流冲击采用大量微粒反复冲击金属构件,引入有益的残余压应力,残余压应力可促使金属表面裂纹闭合,抑制裂纹扩展,提高疲劳裂纹扩展寿命,改善金属构件的抗疲劳性能。而且表层组织发生塑性变形,发生形变诱发马氏体相变,进而产生马氏体相变强化。因此,研究残余应力产生和马氏体相变对推进微粒冲击技术的应用与发展非常关键。由于实验研究费时费力,而且在实际的工程应用中还会受到一定的制约,因此,有限元法成为主要的研究手段。
构建微粒高速冲击金属构件的数值模型时,准确描述金属材料的本构关系至关重要。近年来研究中材料的本构关系多采用Johnson-Cook模型,该模型考虑了应变硬化和应变率硬化的影响,但是没有考虑到马氏体相变强化的影响,残余应力的预测结果并不准确,更无法对高速冲击下金属材料马氏体相变进行研究。因此,提出一种可以准确预测微粒高速冲击下金属构件内残余应力和马氏体相变的材料率相关弹塑性本构模型尤为必要。
发明内容
本发明的主要目的在于克服现有技术中的不足,提供一种能准确预测微粒高速冲击下金属构件残余应力和马氏体相变的有限元方法。
为解决上述技术问题,本发明的解决方案是:一种预测微粒高速冲击下金属构件残余应力和马氏体相变的有限元方法,包括以下步骤:
(1)建立微粒冲击金属构件的有限元模型;
(2)建立综合考虑了塑性应变强化和马氏体相变强化的率相关弹塑性本构模型;
(3)编写ABAQUS-VUMAT用户动态材料子程序,运用应力补偿更新算法实现提出的率相关弹塑性本构模型;
(4)将子程序嵌入到ABAQUS中,对微粒流冲击金属构件的过程进行计算,进一步获得金属构件内部的残余应力以及马氏体相的百分比含量。
进一步的,所述步骤(1)具体为:
基于ABAQUS建立单个或多个球形微粒、金属构件,分别设置材料属性和划分网格,再进行ASSEMBLY对其组装之后设置分析步和通用接触属性。
进一步的,所述步骤(2)具体包括下述步骤:
步骤(2.1)建立金属材料的本构关系:
Δσ=D:Δεe
其中,Δσ为应力增量张量,D为弹性模量张量,Δεe为弹性应变增量张量;
步骤(2.2)建立塑性模型,具体建立方式为:
考虑塑性应变强化和马氏体相变强化的流动准则为:
其中,f为屈服方程;σe为Mises等效应力;所述K和分别为各向同性强化应力和等效塑性应变;所述σ11,σ22,σ33,σ12,σ23和σ31为六个方向的应力分量。
各向同性强化应力:
kw=hwwn
其中,kε为塑性应变强化应力,kw为马氏体相变强化应力,C为材料参数,为等效塑性应变率,为参考应变率,σ0为初始屈服应力,hε为应变强化系数,a为材料参数,hw为相变强化模量,n为强化指数,w为马氏体的百分比含量,其增量方程由Santacreu模型给出:
其中,wmax为马氏体相的最大含量,D0,D1和m为材料参数,为应力三轴度。
进一步的,所述步骤(3)具体包括下述步骤:
步骤(3.1):通过用户子程序VUMAT的SDV定义第n+1增量步开始时的初始状态变量wn值,同时也是第n步结束时的状态变量值,在第n+1增量步开始时读入;其中,所述n是指第n增量步,wn是指第n增量步结束时的金属构件马氏体的百分比含量;
步骤(3.2):VUMAT由应变增量驱动,首先假设应变增量全为弹性,由步骤(2.1)中公式Δσ=D:Δεe计算试探应力:将试探应力代入到屈服准则中,该公式是步骤(2.2)中的公式在第n+1步的特定计算;
如果fn+1≤0,则材料尚处于弹性变形阶段,将试探应力更新为n+1增量步应力,马氏体相含量更新为n+1增量步状态变量wn+1;
如果fn+1>0,则材料屈服,根据应力补偿显式算法更新增量步结束时的应力,即在试探应力的基础上减去多算的塑性应变增量对应的应力增量部分,对试探应力进行负补偿;
等效塑性应变增量运用牛顿迭代法求解:
其中,所述i表示第i次迭代,表示第i次迭代所得的等效塑性应变增量;所述i-1 表示第i-1次迭代,表示第i-1次迭代所得的等效塑性应变增量;更新等效塑性应变增量,直至满足精度要求,得到马氏体相的百分比含量wn+1,更新第n+1增量步应力,更新内能,塑性功。
进一步的,所述步骤(4)具体为:
将步骤(1)建立的模型主文件和步骤(3)建立的ABAQUS-VUMAT用户子程序联合,使用ABAQUS/EXPLICT方法对微粒高速冲击金属构件进行计算,获得马氏体相变含量;将ABAQUS/EXPLICIT中得到的动态应力状态导入到ABAQUS/STANDARD模块中进行隐式分析,确定静态平衡下的残余应力场,即完成微粒流高速冲击下金属构件内残余应力和马氏体相变的预测。
与现有技术相比,本发明的有益效果是:
(1)建立的材料率相关弹塑性本构模型同时考虑了塑性应变强化和马氏体相变强化的影响。
(2)利用ABAQUS-VUMAT用户子程序来数值实现所建立的材料本构模型,能准确预测微粒高速冲击下金属构件残余应力和马氏体相的百分比含量。
(3)基于通用有限元技术,本发明可以用于任意数量的微粒对任意结构形式的金属构件进行冲击的残余应力和马氏体相变预测。
附图说明
图1为本发明对某一具体实例所采用的有限元模型;
图2为本发明对所提出的率相关弹塑性本构模型的VUMAT数值实现流程图;
图3为对图2中实例进行预测得到的残余应力分布结果图;
图4为对图2中实例进行预测得到的马氏体相含量结果图。
具体实施方式
首先需要说明的是,本发明是计算机技术在残余应力和组织相变预测领域的一种应用。在本发明的实现过程中,会涉及到多个软件功能模块的应用。申请人认为,如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后,在结合现有公知技术的情况下,本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。凡本发明申请文件提及的均属此范畴,申请人不再一一列举。
下面结合附图与具体实施方式对本发明作进一步详细描述:
根据微粒和金属板材的对称性,在ABAQUS/CAE中建立包括四个微粒、一个金属板材构件的高速冲击四分之一对称模型,如图1所示,微粒将依次撞击金属板材的四个拐角区域。金属板材的大小为0.5×0.5×2.0mm,球形微粒的直径为1.0mm。固定金属板材的底面,在其四个侧面分别施加相应的对称约束。微粒和金属面之间设置库伦摩擦,摩擦系数为0.3。金属材料为AISI348,密度为7800kg/m3,弹性模量E=210GPa,泊松比μ=0.3。采用减缩积分三维八节点实体单元C3D8R对金属板材和微粒进行离散,金属板材的单元尺寸为0.05mm,微粒的单元尺寸为0.03mm。微粒材料为二氧化硅,其强度和硬度都高于金属板材,因此微粒的变形可以忽略不计,将微粒模拟成刚体,参考点位于其球心,施加初始速度v=100m/s在参考点上。本实例中所用材料参数如表1所示。
表1实例中所用的材料参数
建立材料AISI348的本构关系如下:
Δσ=D:Δεe
其中,Δσ为应力增量张量,D为弹性模量张量,Δεe为弹性应变增量张量;
建立塑性模型,具体建立方式为:
考虑塑性应变强化和马氏体相变强化的流动准则为:
其中,f为屈服方程;σe为Mises等效应力;所述K和分别为各向同性强化应力和等效塑性应变;所述σ11,σ22,σ33,σ12,σ23和σ31为六个方向的应力分量。
各向同性强化应力:
kw=470*w0.97
其中,kε为塑性应变强化应力,kw为马氏体相变强化应力,为等效塑性应变率,为参考应变率,w为马氏体的百分比含量,其增量方程由Santacreu模型给出:
其中,为应力三轴度。
利用ABAQUS/EXPLICIT计算模拟微粒冲击金属板材的过程,如图2利用用户子程序VUMAT判断材料点是否进入塑性,运用应力补偿更新算法实现提出的率相关弹塑性本构模型,计算得到马氏体相变含量。具体过程如下:通过用户子程序VUMAT的SDV定义第n+1 增量步的开始时的初始状态变量wn值,同时也是第n步结束时的状态变量值,在第n+1增量步开始时读入;其中,所述n是指第n增量步,wn是指第n增量步结束时的金属板材内的马氏体的百分比含量。
VUMAT由应变增量驱动,首先假设应变增量全为弹性,由公式Δσ=D:Δεe计算试探应力:将试探应力代入到屈服准则该公式是步骤(2.2)中的公式在第n+1步的特定计算;
如果fn+1≤0,则材料尚处于弹性变形阶段,将试探应力更新为n+1增量步应力,马氏体相含量更新为n+1增量步状态变量wn+1;如果fn+1>0,则材料屈服,根据应力补偿显式算法更新增量步结束时的应力,即在试探应力的基础上减去多算的塑性应变增量对应的应力增量部分,对试探应力进行负补偿;
等效塑性应变增量运用牛顿迭代法求解:
其中,所述i表示第i次迭代,表示第i次迭代所得的等效塑性应变增量;所述i-1 表示第i-1次迭代,表示第i-1次迭代所得的等效塑性应变增量;更新等效塑性应变增量,直至满足精度要求,得到马氏体相的百分比含量wn+1,更新第n+1增量步应力,更新内能,塑性功。
最后,将ABAQUS/EXPLICIT中得到的动态应力状态导入到ABAQUS/STANDARD模块中进行隐式分析,确定静态平衡下金属板材的残余应力场。
图3为沿金属板材沿厚度方向的残余应力分布,随着深度的增加残余压应力先增大后减小;图4为金属板材表面中心点的马氏体相变含量随时间的变化,可见马氏体相变含量随着冲击次数的增加而增多。所以本发明提出的率相关弹塑性本构模型可以准确预测微粒流高速冲击下金属构件内部的残余应力和马氏体相变。
本发明提出的率相关弹塑性本构模型同时考虑了塑性变形和材料相变的影响,结合该本构模型,利用FORTRAN语言编写子程序并嵌入到ABAQUS中,能准确预测微粒高速冲击下金属构件残余应力和马氏体相变,为深入阐明残余应力产生和马氏体相变的机理提供了技术支撑,推进微粒流冲击在提高关键材料和零部件疲劳性能中的应用。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (4)
1.一种预测微粒冲击金属件残余应力和马氏体相变的有限元方法,其特征在于,包括以下步骤:
(1)建立微粒冲击金属构件的有限元模型;
(2)建立综合考虑了塑性应变强化和马氏体相变强化的率相关弹塑性本构模型;
所述步骤(2)具体包括下述步骤:
步骤(2.1)建立金属材料的本构关系:
Δσ=D:Δae
其中,Δσ为应力增量张量,D为弹性模量张量,Δae为弹性应变增量张量;
步骤(2.2)建立塑性模型,具体建立方式为:
考虑塑性应变强化和马氏体相变强化的流动准则为:
其中,f为屈服方程;σe为Mises等效应力;所述K和分别为各向同性强化应力和等效塑性应变;所述σ11,σ22,σ33,σ12,σ13和σ31为六个方向的应力分量;
各向同性强化应力:
kw=hwwn
其中,kε为塑性应变强化应力,kw为马氏体相变强化应力,C为材料参数,为应变率,为参考应变率,σ0为初始屈服应力,hε为应变强化系数,a为材料参数,hw为相变强化模量,n为强化指数,w为马氏体的百分比含量,其增量方程由Santacreu模型给出:
其中,wmax为马氏体相的最大含量,D0,D1和m为材料参数,为应力三轴度;
(3)编写ABAQUS-VUMAT用户动态材料子程序,运用应力补偿更新算法实现提出的率相关弹塑性本构模型;
(4)将子程序嵌入到ABAQUS中,对微粒流冲击金属构件的过程进行计算,进一步获得金属构件内部的残余应力以及马氏体相的百分比含量。
2.根据权利要求1所述的一种预测微粒冲击金属件残余应力和马氏体相变的有限元方法,其特征在于,所述步骤(1)具体为:
基于ABAQUS建立单个或多个球形微粒、金属构件,分别设置材料属性和划分网格,再进行ASSEMBLY对其组装之后设置分析步和通用接触属性。
3.根据权利要求1所述的一种预测微粒冲击金属件残余应力和马氏体相变的有限元方法,其特征在于,所述步骤(3)具体包括下述步骤:
步骤(3.1):通过用户子程序VUMAT的SDV定义第n+1增量步开始时的初始状态变量wn值,同时也是第n步结束时的状态变量值,在第n+1增量步开始时读入;其中,所述n是指第n增量步,wn是指第n增量步结束时的金属构件马氏体的百分比含量;
步骤(3.2):VUMAT由应变增量驱动,首先假设应变增量全为弹性,由步骤(2.1)中公式Δσ=D:Δεe计算试探应力:将试探应力代入到屈服准则中,该公式是步骤(2.2)中的公式在第n+1步的特定计算;
如果fn+1≤0,则材料尚处于弹性变形阶段,将试探应力更新为n+1增量步应力,马氏体相含量更新为n+1增量步状态变量wn+1;
如果fn+1>0,则材料屈服,根据应力补偿显式算法更新增量步结束时的应力,即在试探应力的基础上减去多算的塑性应变增量对应的应力增量部分,对试探应力进行负补偿;
等效塑性应变增量运用牛顿迭代法求解:
其中,所述i表示第i次迭代,表示第i次迭代所得的等效塑性应变增量;所述i-1表示第i-1次迭代,表示第i-1次迭代所得的等效塑性应变增量;更新等效塑性应变增量,直至满足精度要求,得到马氏体相的百分比含量wn+1,更新第n+1增量步应力,更新内能,塑性功。
4.根据权利要求1所述的一种预测微粒冲击金属件残余应力和马氏体相变的有限元方法,其特征在于,所述步骤(4)具体为:
将步骤(1)建立的模型主文件和步骤(3)建立的ABAQUS-VUMAT用户子程序联合,使用ABAQUS/EXPLICT方法对微粒高速冲击金属构件进行计算,获得马氏体相变含量;将ABAQUS/EXPLICIT中得到的动态应力状态导入到ABAQUS/STANDARD模块中进行隐式分析,确定静态平衡下的残余应力场,即完成微粒流高速冲击下金属构件内残余应力和马氏体相变的预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811591935.6A CN109558692B (zh) | 2018-12-25 | 2018-12-25 | 预测微粒冲击金属件残余应力和马氏体相变的有限元方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811591935.6A CN109558692B (zh) | 2018-12-25 | 2018-12-25 | 预测微粒冲击金属件残余应力和马氏体相变的有限元方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109558692A CN109558692A (zh) | 2019-04-02 |
CN109558692B true CN109558692B (zh) | 2019-07-12 |
Family
ID=65871338
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811591935.6A Active CN109558692B (zh) | 2018-12-25 | 2018-12-25 | 预测微粒冲击金属件残余应力和马氏体相变的有限元方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109558692B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795885B (zh) * | 2019-11-22 | 2023-08-04 | 东北大学 | Trip钢动态变形过程相变诱发塑性的有限元模拟方法 |
CN112131762B (zh) * | 2020-08-07 | 2022-10-11 | 上海大学 | 模拟马氏体相变的网格自适应有限元方法 |
CN112711884B (zh) * | 2020-12-29 | 2024-01-23 | 华南理工大学 | 显隐式算法结合的脆性材料准静态破坏的仿真预测方法 |
CN112989654B (zh) * | 2021-02-25 | 2024-07-12 | 江苏大学 | 预测冲击载荷下激光冲击成形极限的有限元方法 |
CN116611204B (zh) * | 2023-03-23 | 2023-10-20 | 哈尔滨工业大学 | 一种模拟含纳米孔洞的钢组织中马氏体相变的分子动力学建模方法 |
CN117556671B (zh) * | 2023-11-24 | 2024-04-30 | 中国石油大学(华东) | 爆炸条件下临氢金属承压结构动态损伤与断裂预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102122311A (zh) * | 2011-02-21 | 2011-07-13 | 北京航空航天大学 | 一种基于有限元的动力调谐陀螺加速稳定剖面生成方法 |
CN106777769A (zh) * | 2017-01-08 | 2017-05-31 | 浙江大学 | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10303827B2 (en) * | 2016-04-05 | 2019-05-28 | Rolls-Royce Corporation | Predicting cracking in cooled metal or alloy components |
CN106503292B (zh) * | 2016-09-20 | 2018-04-24 | 浙江大学 | 预测低速冲击下复合材料层合板渐进失效的有限元方法 |
CN108932362B (zh) * | 2018-05-07 | 2023-01-06 | 华南理工大学 | 一种预测高压氢系统橡胶密封圈密封特性的有限元方法 |
CN108637525A (zh) * | 2018-05-17 | 2018-10-12 | 中国石油大学(华东) | 一种可免焊后去应力处理的低合金高强度钢用埋弧焊丝 |
-
2018
- 2018-12-25 CN CN201811591935.6A patent/CN109558692B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102122311A (zh) * | 2011-02-21 | 2011-07-13 | 北京航空航天大学 | 一种基于有限元的动力调谐陀螺加速稳定剖面生成方法 |
CN106777769A (zh) * | 2017-01-08 | 2017-05-31 | 浙江大学 | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109558692A (zh) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109558692B (zh) | 预测微粒冲击金属件残余应力和马氏体相变的有限元方法 | |
CN109946006B (zh) | 基于混合硬化模型的微粒流冲击金属材料力学行为预测方法 | |
Wang et al. | Dislocation-based study on the influences of shot peening on fatigue resistance | |
Meo et al. | Finite element analysis of residual stress induced by shot peening process | |
Chaise et al. | Modelling of multiple impacts for the prediction of distortions and residual stresses induced by ultrasonic shot peening (USP) | |
Guagliano | Relating Almen intensity to residual stresses induced by shot peening: a numerical approach | |
Miao et al. | Numerical simulation of the stress peen forming process and experimental validation | |
Shih et al. | A model for hysteretic behavior of rhombic low yield strength steel added damping and stiffness | |
Shivpuri et al. | Elasto-plastic pseudo-dynamic numerical model for the design of shot peening process parameters | |
Kang et al. | Multiple impact modelling for shot peening and peen forming | |
Xiao et al. | Prediction of shot peen forming effects with single and repeated impacts | |
Yang et al. | Effect of initial surface finish on effectiveness of shot peening treatment using enhanced periodic cell model | |
Burchitz | Improvement of springback prediction in sheet metal forming | |
Bhuvaraghan et al. | Optimization of the fatigue strength of materials due to shot peening: A Survey | |
Wu et al. | Effect of shot peening on surface residual stress distribution of SiCp/2024Al | |
Xiao et al. | Estimation of peening effects of random and regular peening patterns | |
Tamarkin et al. | Surface quality assurance and process reliability in the processing with a ball-rod hardener | |
Zhang et al. | Discrete element simulation of surface mechanical attrition treatment with rough-surface sonotrode | |
Zhao et al. | Analytical and experimental study on the surface generation mechanism in shot peening | |
Xiao et al. | Prediction of peen forming stress and curvature with dynamic response of compressively prestressed target | |
Farrahi | Side effects of shot peening on fatigue crack initiation life | |
Balokhonov et al. | Simulation of deformation and fracture of coated material with account for propagation of a Lüders-Chernov band in the steel substrate | |
Zhang et al. | Modeling of multiple shots for analyzing shot peening controlled parameters on formed curvature radius | |
Wang et al. | Mechanical behavior study of asymmetric deformation in double-sided symmetrical sequential shot peening process | |
Yeganeh-Salman et al. | Evaluating the accuracy of mass scaling method in non-linear quasi-static finite element analysis of RC structures |
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 |