用于计算合金表面划伤缺陷底部应力集中系数的方法技术

本发明专利技术公开了一种用于计算合金表面划伤缺陷底部应力集中系数的方法，首先建立缺口问题的弹性力学模型，得到不同深度缺口对应的应力集中函数；然后使用有限元软件ABAQUS建立合金表面划伤和后继加载仿真计算模型，划伤模型采用断裂应变控制工件被划部分单元的破坏和消除，后继加载在划伤模型基础上计算，自然考虑了划伤处的局部残余应力和局部塑性应变对应力集中的影响；最后根据有限元数值分析得到的局部应力结果，对应力集中系数表达式进行修正，得到应力集中系数修正公式。本发明专利技术提供了评估合金材料表面划伤的应力集中系数的方法和公式，为含表面划伤缺陷的该种材料结构的疲劳寿命分析提供了依据。

【技术实现步骤摘要】
用于计算合金表面划伤缺陷底部应力集中系数的方法
本专利技术涉及金属结构和材料疲劳断裂研究领域，尤其涉及一种用于计算合金表面划伤缺陷底部应力集中系数的方法。
技术介绍
金属(合金)结构的疲劳寿命与其存在的缺陷密切相关，航空结构的设计正从传统的安全寿命方法向损伤容限和缺陷容限法转变。缺陷容限分析依赖于缺陷处的局部应力、应变状态的分析。表面划伤是金属结构最常遇到的缺陷之一，在制造和使用阶段均可能发生，会对其疲劳寿命造成很大影响。表面划伤不仅会产生应力集中，而且在划伤的过程中还造成局部的残余应力和塑性变形，这些都会对应力集中产生影响。目前划伤等缺陷造成的应力集中分析部完善，现有的一些工作是完全基于弹性力学分析的，没有考虑上述局部的残余应力和塑性变形影响，因此实际上不能准确反映真实的情况。
技术实现思路
本专利技术所要解决的技术问题是针对
技术介绍
中所涉及到的缺陷，提供一种用于计算合金表面划伤缺陷底部应力集中系数的方法。本专利技术为解决上述技术问题采用以下技术方案：用于计算合金表面划伤缺陷底部应力集中系数的方法，包含以下步骤：步骤1)，根据缺口附近的弹性应力场函数，进而推出应力分量的方程，通过引入边界条件求得对称拉伸载荷下缺口附近的应力场的通解，得到缺口底部应力集中系数的表达式，得到不同底部半径和深度的缺口对应的应力集中函数；步骤2)，使用ABAQUS软件创建刀具和工件模型，设置刀具为刚体，工件为弹塑性体，通过vumate子程序设置断裂应变控制单元的破坏和消除，刀具与工件之间定义面面接触,划痕的尺寸由底部半径R和深度D表征；非线性动力学计算产生划痕后，采用重启动分析方法导入划伤后的模型，此模型包含了划伤过程的残余应力和塑性应变信息，随后施加垂直于划痕方向的拉伸载荷，划伤模型采用断裂应变控制工件被划部分单元的破坏和消除，后继加载在划伤模型基础上计算，考虑了划伤处的局部残余应力和局部塑性应变对应力集中的影响，得到有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果；步骤3)，对弹性力学缺口底部应力集中系数表达式进行修正，引入修正参数；根据有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果，拟合出修正参数随R、D的变化函数，根据有限元数值分析得到的局部应力结果，对步骤1)中的应力集中系数表达式进行修正，得到修正过后的应力集中系数。作为本专利技术用于计算合金表面划伤缺陷底部应力集中系数的方法进一步的优化方案，步骤1)的详细步骤如下：令缺口附近弹性应力场为：ψ(z)＝bzλ+czμ式中，a、b、c、d是复数形式，λ和μ是实数，λ＞0且λ＞μ；推出应力分量的表达式如下：σθ＝λrλ-1[a1(1+λ)cos(1-λ)θ+b1cos(1+λ)θ+a2(1+λ)sin(1-λ)θ-b2sin(1+λ)θ]+μrμ-1[d1(1+μ)cos(1-μ)θ+c1cos(1+μ)θ+d2(1+μ)sin(1-μ)θ-c2sin(1+μ)θ]σr＝λrλ-1[a1(3-λ)cos(1-λ)θ-b1cos(1+λ)θ+a2(3-λ)sin(1-λ)θ+b2sin(1+λ)θ]+μrμ-1[d1(3-μ)cos(1-μ)θ-c1cos(1+μ)θ+d2(3-μ)sin(1-μ)θ+c2sin(1+μ)θ]τrθ＝λrλ-1[a1(1-λ)sin(1-λ)θ+b1sin(1+λ)θ-a2(1-λ)cos(1-λ)θ+b2cos(1+λ)θ]+μrμ-1[d1(1-μ)sin(1-μ)θ+c1sin(1+μ)θ-d2(1-μ)cos(1-μ)θ+c2cos(1+μ)θ]式中，(r、θ)表示合金表面任一点的极坐标，σθ、σr和τrθ表示该点θ、r以及切向的应力分量，a1、a2分别是a的实部、虚部，b1、b2分别是b的实部、虚部，c1、c2分别是c的实部、虚部，d1、d2分别是d的实部、虚部；引入以下边界条件：式中，(u，v)表示合金表面任一点在辅助坐标系中的坐标，q为从1到2的实数；r0、u0、v0分别为缺口圆弧段深度以及缺口自由边界的坐标；获得缺口在对称拉伸载荷作用下缺口附近应力场的通解，在缺口底部点r＝r0，θ＝0处的表达式为：式中，R表示缺口底部圆弧半径；χb、χc、χd表达式如下：式中，则缺口底部点应力集中系数表达式为：式中，Kθ、Kr、Krθ分别表示任一点θ、r以及切向的的应力集中系数。作为本专利技术用于计算合金表面划伤缺陷底部应力集中系数的方法进一步的优化方案，步骤3)的详细步骤如下：步骤3.1)，对缺口底部应力集中系数表达式进行修正，引入修正参数a1、a2：步骤3.2)，根据有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果，分别拟合出λ、μ、χb、χc、χd和弹性力学缺口的张角2α之间的函数关系式；步骤3.3)，对于有限元模拟的多组不同底部半径R、深度D的划伤，提取其模型底部残余应力以及后继加载划痕底部应力分量结果，求得应力集中系数Kθ、Kr，进而得到多组不同底部半径R、深度D的划伤对应的修正参数a1、a2的值；步骤3.4)，根据有限元模拟的多组不同底部半径R、深度D的划伤对应的修正参数a1、a2的值，采用多项式形式分别对a1、a2进行拟合，得到a1、a2的拟合函数，进而得到修正过后的应力集中系数。本专利技术采用以上技术方案与现有技术相比，具有以下技术效果：1.将V型缺口附近应力场的理论结果应用到合金表面划痕附近应力场的计算，具有针对性；2.结合有限元模拟，考虑了划痕的塑性，更符合实际情况。附图说明图1为缺口弹性理论分析的模型示意图；图2为本专利技术的实现流程图；图3(a)为本专利技术表面划伤模拟中整体网格的有限元模型；图3(b)为本专利技术表面划伤模拟中局部网格的有限元模型；图4为系数λ的拟合曲线；图5为系数μ的拟合曲线；图6为系数χb的拟合曲线；图7为系数χc的拟合曲线；图8为系数χd的拟合曲线；图9为系数a1随划痕深度D、底部半径R的变化曲面图；图10为系数a2随划痕深度D、底部半径R的变化曲面图。具体实施方式下面结合附图对本专利技术的技术方案做进一步的详细说明：本专利技术可以以许多不同的形式实现，而不应当认为限于这里所述的实施例。相反，提供这些实施例以便使本公开透彻且完整，并且将向本领域技术人员充分表达本专利技术的范围。在附图中，为了清楚起见放大了组件。如图1、图2所示，本专利技术公开了一种用于计算合金表面划伤缺陷底部应力集中系数的方法，以TB6钛合金为例，首先推导对称拉伸载荷下V型缺口底部点的应力分量表达式，具体过程为：缺口附近弹性应力场ψ(z)＝bzλ+czμ可推出应力分量的表达式如下σθ＝λrλ-1[a1(1+λ)cos(1-λ)θ+b1cos(1+λ)θ+a2(1+λ)sin(1-λ)θ-b2sin(1+λ)θ]+μrμ-1[d1(1+μ)cos(1-μ)θ+c1cos(1+μ)θ+d2(1+μ)sin(1-μ)θ-c2sin(1+μ)θ]σr＝λrλ-1[a1(3-λ)cos(1-λ)θ-b1cos(1+λ)θ+a2(3-λ)sin(1-λ)θ+b2sin(1+λ)θ]+μrμ-1[d1(3-μ)cos(1-μ)θ-c1cos(1+μ)θ+d2(3-μ)sin(1-μ)θ+c本文档来自技高网...

【技术保护点】
1.用于计算合金表面划伤缺陷底部应力集中系数的方法，其特征在于，包含以下步骤：步骤1)，根据缺口附近的弹性应力场函数，进而推出应力分量的方程，通过引入边界条件求得对称拉伸载荷下缺口附近的应力场的通解，得到缺口底部应力集中系数的表达式，得到不同底部半径和深度的缺口对应的应力集中函数；步骤2)，使用ABAQUS软件创建刀具和工件模型，设置刀具为刚体，工件为弹塑性体，通过vumate子程序设置断裂应变控制单元的破坏和消除，刀具与工件之间定义面面接触,划痕的尺寸由底部半径R和深度D表征；非线性动力学计算产生划痕后，采用重启动分析方法导入划伤后的模型，此模型包含了划伤过程的残余应力和塑性应变信息，随后施加垂直于划痕方向的拉伸载荷，划伤模型采用断裂应变控制工件被划部分单元的破坏和消除，后继加载在划伤模型基础上计算，考虑了划伤处的局部残余应力和局部塑性应变对应力集中的影响，得到有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果；步骤3)，对弹性力学缺口底部应力集中系数表达式进行修正，引入修正参数；根据有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果，拟合出修正参数随R、D的变化函数，根据有限元数值分析得到的局部应力结果，对步骤1)中的应力集中系数表达式进行修正，得到修正过后的应力集中系数。...

【技术特征摘要】
1.用于计算合金表面划伤缺陷底部应力集中系数的方法，其特征在于，包含以下步骤：步骤1)，根据缺口附近的弹性应力场函数，进而推出应力分量的方程，通过引入边界条件求得对称拉伸载荷下缺口附近的应力场的通解，得到缺口底部应力集中系数的表达式，得到不同底部半径和深度的缺口对应的应力集中函数；步骤2)，使用ABAQUS软件创建刀具和工件模型，设置刀具为刚体，工件为弹塑性体，通过vumate子程序设置断裂应变控制单元的破坏和消除，刀具与工件之间定义面面接触,划痕的尺寸由底部半径R和深度D表征；非线性动力学计算产生划痕后，采用重启动分析方法导入划伤后的模型，此模型包含了划伤过程的残余应力和塑性应变信息，随后施加垂直于划痕方向的拉伸载荷，划伤模型采用断裂应变控制工件被划部分单元的破坏和消除，后继加载在划伤模型基础上计算，考虑了划伤处的局部残余应力和局部塑性应变对应力集中的影响，得到有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果；步骤3)，对弹性力学缺口底部应力集中系数表达式进行修正，引入修正参数；根据有限元模拟的多组不同底部半径R、深度D的划伤后继加载底部应力集中系数的数值结果，拟合出修正参数随R、D的变化函数，根据有限元数值分析得到的局部应力结果，对步骤1)中的应力集中系数表达式进行修正，得到修正过后的应力集中系数。2.根据权利要求1所述的用于计算合金表面划伤缺陷底部应力集中系数的方法，其特征在于，步骤1)的详细步骤如下：令缺口附近弹性应力场为：ψ(z)＝bzλ+czμ式中，a、b、c、d是复数形式，λ和μ是实数，λ＞0且λ＞μ；推出应力分量的表达式如下：σθ＝λrλ-1[a1(1+λ)cos(1-λ)θ+b1cos(1+λ)θ+a2(1+λ)sin(1-λ)θ-b2sin(1+λ)θ]+μrμ-1[d1(1+μ)cos(1-μ)θ+c1cos(1+μ)θ+d2(1+μ)sin(1-μ)θ-c2sin(1+μ)θ]σr＝λrλ-1[a1(3-λ)cos(1-λ)θ-b1cos(1+λ)θ+a2(3-λ)sin(1-λ)θ+b2sin(1+λ)θ]+μr...

【专利技术属性】
技术研发人员：倪阳，周储伟，周世友，郝建群，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32