CN113408102A - 一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法 - Google Patents

Abstract

一种空心球‑金属基三相复合材料细观尺度精细化建模仿真方法，属于三相复合材料领域。首先，复现等效体心立方均匀化模型的单一胞元结构，建立单一胞元结构与三相复合材料结构的映射关系；其次，根据目标三相复合材料结构的外形尺寸、三相复合材料的球基分布规律和球基间拓扑结构关系创建目标结构几何模型；最后，将三相复合材料结构几何模型离散为各种单元后进行位移插值，为离散单元赋予材料属性和节点约束条件，获得三相复合材料细观尺度精细化仿真模型。本发明能够得到与真实三相复合材料极为相似的胞元结构，该结构可准确再现空心球复合金属泡沫材料的各组成相结构分布；可研究固态金属基体和空心球界面间的连接失效和基体材料的裂纹拓展等微观变化。

Description

一种空心球-金属基三相复合材料细观尺度精细化建模仿真 方法

技术领域

本发明属于三相复合材料领域，涉及一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法。

背景技术

空心球-金属基三相复合材料是一种先进的多孔材料，由胞元尺寸、形状和壁厚为设定值的金属空心球和固态金属基体组成，具有高强度、高刚度和优异的吸声性能，因此其在汽车、航天、能源工程和医疗器械等领域具有广泛的应用前景。相较于传统泡沫金属材料，空心球-金属基三相复合材料具有均匀分布的胞元结构和牢固的球基结合强度。因此，它作为一种各向同性材料具有更高的强度、刚度和能量吸收能力。然而，空心球-金属基三相复合材料的细观结构和力学性能十分复杂，其精细化建模仿真方法是研究该材料变形模式和损伤机理的基础，也是当前该领域的研究重点。

空心球-金属基三相复合材料的建模分析可分为基于宏观尺度和微观尺度的两类方法，这两类方法分别从不同角度定义和阐述空心球-金属基三相复合材料的材料特性。宏观尺度上的建模方法忽略空心金属球的影响，主要基于空心球-金属基三相复合材料在多种应变率下的应力应变曲线，采用均匀化方法等效具有复杂材料属性的金属基三相复合材料。该方法虽可以在宏观尺度反映材料在目标工况下的力学特性，但无法准确表征材料各组成相的细观结构和变形损伤机理。微观尺度上的建模方法忽略了空心球和基体材料间的物质差异，主要基于基体材料线弹性应力应变曲线，建立空心球-金属基三相复合材料无球壁的仿真模型。该方法可以在一定程度上反映给定工况下金属孔隙变形情况，但因其忽略空心球壁对基体材料的支撑和连接作用，故无法给出各组成相的失效行为，即无法界定基体和空心球界面间的连接失效和基体材料的裂纹拓展路径等微观过程。

外国学者Afsaneh Rabiei等在2019年发表的《Ballistic performance ofcomposite metal foam against large caliber threats Composite Structures》中，基于空心球-金属基三相复合材料在多种应变率下的应力应变曲线，忽略空心金属球和多孔结构的影响，采用均匀化方法等效具有复杂材料属性的空心球-金属基三相复合材料。该方法虽可以在宏观尺度反映材料在目标工况下的力学特性，但其缺点是无法对材料各组成相的细观结构和变形损伤机理准确表征。因此，研究空心球-金属基三相复合材料亟需建立不同孔隙率下能反映各组成相结构分布和损伤演化的细观尺度模型，用于研究基体和空心球界面间的连接失效以及基体材料的裂纹拓展路径，即通过等效体心立方均匀化方法建立空心球-金属基三相复合材料细观尺度模型。

发明内容

本发明提出一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，建立空心球壁、固态金属基体复合材料的等效体心立方均匀化细观尺度模型。该模型可准确再现空心球-金属基三相复合材料的各组成相结构分布，用于研究基体和空心球界面间的连接失效和基体材料的裂纹拓展路径等微观过程，这对提高金属基三相复合材料数值模拟的准确性有深远意义。

为了达到上述目的，本发明采用的技术方案为：

一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，首先，根据空心球-金属基三相复合材料各组成相的几何参数、目标孔隙率复现等效体心立方均匀化模型的单一胞元结构，建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系；其次，根据目标空心球-金属基三相复合材料结构的外形尺寸、空心球-金属基三相复合材料的球基分布规律和球基间拓扑结构关系创建目标结构几何模型；最后，将空心球-金属基三相复合材料结构几何模型离散为各种单元，并在离散的单元节点上进行位移插值，为离散单元赋予材料属性和节点约束条件。具体如下：

第一步，根据空心球-金属基三相复合材料各组成相的几何参数、目标孔隙率复现等效体心立方均匀化模型的单一胞元结构，建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系，具体方法如下：

(1)根据空心球-金属基三相复合材料各组成相的几何参数、目标孔隙率复现等效体心立方均匀化模型的单一胞元结构：

因空心球-金属基三相复合材料的孔隙率介于体心立方堆积孔隙率和简单立方堆积孔隙率之间，且细观尺度上空心球在基体材料中的位置分布与体心立方堆积方式较一致，故以等效均匀化体心立方形式排布空心球-金属基三相复合材料的球基，并通过增加空心球间距以降低体心立方堆积孔隙率，从而得到与真实空心球-金属基三相复合材料极为相似的胞元结构。

空心球-金属基三相复合材料的孔隙率为：

其中，V为目标空心球-金属基三相复合材料的体积，单位mm³；V_m为空心球-金属基三相复合材料的基体材料体积，单位mm³；V_s为空心球-金属基三相复合材料的空心球体积，单位mm³；V_h为空心球-金属基三相复合材料孔隙的体积，单位mm³。

鉴于该空心球-金属基三相材料细观模型为等效体心立方均匀化模型，故从整体模型中提取一个等效体心立方体单胞胞元结构，该立方体单胞胞元结构与空心球-金属基三相材料具有相同的孔隙率：

其中，V_d为立方体胞元的体积，单位mm³；V_dm为立方体胞元的基体体积，单位mm³；V_ds为立方体胞元的空心球体积，单位mm³；V_dh为立方体胞元孔隙的体积，单位mm³。

立方体单胞胞元体积为：

V_d＝a³ (3)

其中，a为立方体胞元的边长，单位mm。

立方体胞元中孔隙的体积为：

其中，R为空心球外部半径，单位mm；t为空心球壁厚，单位mm。

进而，由公式(2)、(3)、(4)可推导出单一立方体单胞胞元边长为：

通过空心球-金属基三相复合材料的目标孔隙率、空心球外径尺寸和空心球壁厚即可计算得到单一立方体胞元的边长。

(2)建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系。

假设目标结构模型为长方体，目标模型在笛卡尔坐标系下放置，为消除空心球-金属基三相复合材料的边界效应，其仿真模型任意方向空心球数量和单一胞元数量应满足：

N_dmin＝N_min-1 (6)

其中，N_dmin为消除边界效应时目标模型任意方向单一胞元数量最小值；N_min为消除边界效应时目标模型任意方向空心球数量最小值。

则该模型几何参数为：

L_x＝m·a(N_dmin≤m≤m_max) (7)

L_y＝n·a(N_dmin≤n≤n_max) (8)

L_z＝p·a(N_dmin≤p≤p_max) (9)

其中，L_x为目标模型在x方向长度，单位mm；L_y为目标模型在y方向长度，单位mm；L_z为目标模型在z方向长度，单位mm；m为目标模型x方向单一立方体胞元的个数；n为目标模型y方向单一立方体胞元个数；p为目标模型z方向单一立方体胞元个数；m_max为目标模型x方向单一立方体胞元的个数最大值；n_max为目标模型y方向单一立方体胞元的个数最大值；p_max为目标模型z方向单一立方体胞元的个数最大值。

即建立单一胞元结构与目标空心球-金属基三相复合材料结构的数量和几何映射关系。

第二步，根据目标空心球-金属基三相复合材料结构的外形尺寸、空心球-金属基三相复合材料的球基分布规律和球基间拓扑结构关系创建目标结构几何模型，具体如下：

(1)在笛卡尔坐标系下，以(0，0，0)、(L_x，0，0)、(0，L_y，0)、(0，0，L_z)为基础点创建目标样件几何实体；

(2)在y＝0的X-O-Z平面创建第一层空心球坐标点，即：

P_i,j,k＝(ia,ja,ka) (i＝0,1,···,m,k＝0,1,···,p,j＝0) (10)

(3)在

的X-O-Z平面创建第二层空心球坐标点，即：

(4)将两层坐标点沿Y轴正方向线性复制阵列，阵列数量为n，阵列间距为a，单位mm；再将第一层坐标点沿Y轴正方向线性复制阵列，阵列数量为2，阵列节距为na，单位mm，最终获得空心球球心位置坐标点。

(5)在空心球球心坐标点处创建实心球体，球体半径为空心球内部半径r：

r＝R-t (12)

通过布尔减法运算去除样件实体中在球心坐标点处创建的半径为r的球形实体，获得多孔实体模型。

(6)由于基体与空心球间的拓扑结构关系，两部分的网格节点需要相互拟合并传递节点力。因此，首先在(4)中生成的球心位置坐标处创建空心球外壁曲面，曲面半径与空心球外部半径尺寸相同；再通过曲面切割实体，利用空心球外壁曲面切割(5)中获得的多孔实体模型。

通过切割可以得到壁厚为(R-r)的空心球和多孔基体，两者间曲面的拓扑状态为共享曲面，以共享曲面为目标面创建的2D网格可以实现网格节点一一对应，以备下一步精准地网格划分。根据以上方法即可得到空心球-金属基三相复合材料结构的几何模型。

第三步，将空心球-金属基三相复合材料结构几何模型离散为各种单元，并在离散的单元节点上进行位移插值，为离散单元赋予材料属性和节点约束条件。具体为：

(1)将空心球-金属基三相复合材料结构几何模型离散为各种单元，并在离散的单元节点上进行位移插值。

将空心球-金属基三相复合材料结构几何模型各部分节点的位置与时间的变化关系用离散单元节点的坐标插值表示为：

其中，i＝1,2,3为空间坐标指标；I＝1,2,···,n表示有限单元节点；其中μ_iI表示单元第I个节点在i方向的位移；N_I(x,y,z)为第I个节点的拉格朗日插值函数。

由于空心球壁、球间距为狭小尺寸，为防止出现局部应力集中，对空心球-金属基三相复合材料进行三角形二维单元面网格处理，使面网格与复杂的几何表面属性一致，且满足：

p≤t (14)

其中，p为二维网格尺寸；t为空心球壁厚。

为保证计算精度，选择十节点四面体单元离散第二步中生成的空心球-金属基三相复合材料几何模型，在结构紧凑区域(其中，在空心球壁及空心球之间的区域为结构紧凑区域)自动生成较小的网格单元，最大程度保持网格表面与第二步中生成的目标结构几何模型的一致性，在离散实体时确保四面体网格与二维网格具有相关性，且满足：

p’≤t (15)

其中，p’为十节点四面体单元尺寸；t为空心球壁厚。

经过离散后，即可得到空心球-金属基三相复合材料离散单元，其中包括以空心球几何结构离散得到的空心球离散单元和以多孔基体几何结构离散得到的固态基体离散单元。

(2)为空心球-金属基三相复合材料离散单元赋予材料属性和节点约束条件。具体如下：

(a)为空心球离散单元和固态基体离散单元设置屈服准则、材料参数和失效准则。

因空心球、基体材料均为各向同性材料，所以遵从各向同性Von Mises屈服准则，因此为空心球离散单元和固态基体离散单元均赋予各向同性Von Mises屈服准则，即无论在任何应力状态下，当变形体内某一点的应力偏张量的第二不变量J'₂达到某一定值时，该点由弹性状态进入塑性状态。主坐标系下的Von Mises屈服准则为：

其中，σ₁、σ₂、σ₃为笛卡尔坐标系下X、Y、Z方向的单元主应力，σ_s为单元屈服应力。其物理意义为，当空心球、固态基体离散单元体积形状改变的弹性模量达到某一常数时，该单元质点发生屈服。

设空心球、基体材料均为各向同性理想弹塑性金属材料，分别为空心球离散单元和固态基体离散单元定义材料的杨氏模量、屈服强度、质量密度、泊松比、切线模量、失效参数和等效应力-应变曲线。

空心球和固态基体的离散单元均以材料失效时的有效塑性应变和最小时间步定义材料失效指标，当离散单元等效塑性应变达到失效值时，离散单元即发生失效。

设置离散单元属性和截面属性，创建十节点四面体实体单元属性并赋予给空心球、固态基体离散单元。

(b)为空心球-金属基三相复合材料离散单元节点赋予节点约束条件，即在空心球-金属基三相复合材料离散单元节点处建立接触约束。

空心球-金属基三相复合材料各组成相间具有复杂的接触特性，因此空心球与固态基体离散单元节点间的接触摩擦系数由多组参数共同决定：

其中，μ_c为接触摩擦系数；FD为动摩擦因数；FS为静摩擦因数；DC为指数衰减系数；v_rel为接触面相对速度。

同时鉴于空心球和基体组成相的失效模式，在仿真两者单元节点间接触关系时需要定义带有失效准则的绑定接触。该绑定接触基于罚函数的接触算法，在接触相离散单元节点间传递压应力和拉应力，对无穿透的从节点施加接触力使从节点和主面段穿透距离为零，当达到接触失效准则后绑定接触失效，空心球和基体单元节点间的接触类型转换为普通接触。接触失效准则为：

其中，σ_n和σ_s分别为接触面正应力和接触面切应力，单位MPa；NFLS与SFLS分别为对应的接触面失效正应力和接触面失效剪切应力。

当接触面单元间达到界面极限分离距离后，接触点完全破坏，空心球与基体间界面完全分离，两者间脆性金属间相完全损伤失效并形成球壁外侧微裂纹。

通过以上三个步骤即可获得空心球-金属基三相复合材料细观尺度精细化仿真模型。

通过显式求解方法计算目标工况下材料的力学响应和损伤模式，建立材料的仿真模型裂纹观测截面，通过等效应力、等效应变、最大剪应力应力云图观测截面裂纹分布，在细观尺度下得到空心球-金属基三相复合材料的界面失效和裂纹形式，即得到空心球-金属基三相复合材料细观尺度仿真模型的力学数值解。

本发明的有益效果是：

(1)采用等效均匀化体心立方排布方式描述空心球-金属基三相复合材料的球基排布，建立空心球壁、固态基体细观尺度模型；通过增加空心球间距以降低体心立方堆积孔隙率，得到与真实空心球-金属基三相复合材料极为相似的胞元结构，该结构模型可以准确再现空心球复合金属泡沫材料的各组成相结构分布。

(2)考虑空心金属球和金属基体间连接方式的影响，仿真空心球-基体连接处接触特性和材料差异，可以研究固态金属基体和空心球界面间的连接失效和基体材料的裂纹拓展等微观变化。

附图说明

图1为本发明等效体心立方均匀化模型单一胞元结构简化图；图1(a)为单一胞元结构侧视图；图1(b)为单一胞元结构正等测图。

图2为本发明单胞结构与目标空心球-金属基三相复合材料结构映射关系图。

图3为本发明单层空心球圆心坐标分布图；图3(a)为y＝0的X-O-Z平面第一层空心球坐标分布图；图3(b)为

的X-O-Z平面第二层空心球坐标分布图。

图4为本发明目标基体结构几何模型正等测图。

图5为本发明目标空心球-金属基三相复合材料结构几何模型正等测图。

图6为本发明空心球-金属基三相复合材料模型离散单元图。

图7为本发明细观尺度下空心球-金属基三相复合材料的界面失效、应力分布裂纹分布图；

图7(a)为模型整体界面失效、应力分布和裂纹分布图；图7(b)、(c)为仿真模型裂纹观测截面不同阶段的界面失效、应力分布和裂纹分布图。

具体实施方式

下面结合附图和实施例详细说明本发明的结构原理和工作原理。

一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，其特征在于：

本发明的基本变量包括空心球-金属基三相复合材料的孔隙率、空心球外部半径、空心球壁厚、目标几何结构尺寸、基准坐标和三相复合材料参数；参照图5，本发明中选取的空心球-金属基三相复合材料的孔隙率为f＝0.39；空心球外部半径R＝1mm；空心球壁t＝0.1mm；目标几何结构尺寸L_x＝7.5mm、L_y＝22.5mm、L_z＝7.5mm；基准坐标为(0，0，0)；空心球-金属基三相复合材料中空心球材料为不锈钢316L，基体材料为AL7075。具体空心球-金属基三相材料等效体心立方均匀化模型相关参数计算如下：

(1)参见图1，等效体心立方均匀化模型单一胞元结构中的孔隙体积V_dh：

根据孔隙率f＝0.39，计算立方体单胞的体积V_d：

立方体单一胞元结构的边长a:

(2)参见图2，建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系。

目标几何模型x方向单一立方体胞元的个数m：

目标几何模型y方向单一立方体胞元的个数n：

目标几何模型z方向单一立方体胞元的个数p：

(3)参见图3，根据空心球-金属基三相复合材料的球基分布规律创建目标结构在笛卡尔坐标系下的三维几何模型。

在笛卡尔坐标系下，以(0，0，0)、(2.5，0，0)、(0，22.5，0)、(0，0，2.5)为基础点创建目标结构几何模型。

在y＝0的X-O-Z平面创建第一层空心球坐标点，即：

P_i,j,k＝(2.5i,2.5j,2.5k)

(i＝0,1,···,m,k＝0,1,···,p,j＝0) (7)

在y＝1.25的X-O-Z平面创建第二层空心球坐标点，即：

将两层坐标点沿Y轴正方向线性复制阵列，阵列数量为n＝9，阵列间距为2.5，单位mm；再将第一层坐标点沿Y轴正方向线性复制阵列，阵列数量为2，阵列节距为22.5，单位mm，最终获得空心球球心位置坐标点。

(4)参见图4，在空心球球心坐标点处创建实心球体，球体半径为空心球内部半径r：

r＝R-t＝1-0.1＝0.9mm (9)

通过布尔减法运算去除样件实体中在球心坐标点处创建的半径为0.9mm的球形实体，获得多孔实体几何模型。

(5)参见图5，根据基体与空心球间的拓扑结构关系，在(3)中生成的球心位置坐标处创建空心球外壁曲面，曲面半径为R＝1mm；再通过曲面切割实体，利用空心球外壁曲面切割(4)中获得的多孔实体模型。

通过切割可以得到壁厚为t＝0.1mm的空心球和多孔基体，即空心球-金属基三相复合材料结构的几何模型。

(6)参见图6，将(5)中生成的空心球-金属基三相复合材料结构几何模型离散为各种单元，选择十节点四面体单元离散空心球-金属基三相复合材料模型，单元尺寸p：

p≤0.1 (10)

(7)空心球-金属基三相复合材料离散单元赋予材料属性和节点约束条件，空心球、基体材料均为各向同性理想弹塑性金属材料，空心球材料为不锈钢316L，基体为AL7075。

表1空心球-金属基三相复合材料模型材料属性

空心球和固态基体的离散单元均以材料失效时的有效塑性应变和最小时间步定义材料失效指标，当离散单元等效塑性应变达到失效值时，离散单元即发生失效。

定义空心球与金属基体离散单元节点间接触关系为带有失效准则的绑定接触。当接触面单元间达到界面极限分离距离后，接触点完全破坏，空心球与基体间界面完全分离，两者间脆性金属间相完全损伤失效并形成球壁外侧微裂纹。

(8)参见图7，获得目标工况下三相材料的力学响应和损伤模式，得到材料的仿真模型裂纹观测截面，在细观尺度下观测空心球-金属基三相复合材料的界面失效和裂纹形式。

以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。

Claims (3)

1.一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，其特征在于，包括以下步骤：

第一步，首先复现等效体心立方均匀化模型的单一胞元结构，建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系，具体方法如下：

(1)根据空心球-金属基三相复合材料各组成相的几何参数、目标孔隙率复现等效体心立方均匀化模型的单一胞元结构：

以等效均匀化体心立方形式排布空心球-金属基三相复合材料的球基，通过增加空心球间距降低体心立方堆积孔隙率，得到与真实空心球-金属基三相复合材料极为相似的胞元结构；从整体模型中提取一个等效体心立方体单胞胞元结构，该立方体单胞胞元结构与空心球-金属基三相材料具有相同的孔隙率，为：

其中，V_d为立方体胞元的体积，单位mm³；V_dm为立方体胞元的基体体积，单位mm³；V_ds为立方体胞元的空心球体积，单位mm³；V_dh为立方体胞元孔隙的体积，单位mm³；

所述单一立方体单胞胞元的边长为：

其中，a为立方体胞元的边长，单位mm；R为空心球外部半径，单位mm；t为空心球壁厚，单位mm；

(2)建立单一胞元结构与目标空心球-金属基三相复合材料结构的映射关系；

假设目标结构模型为长方体，目标模型在笛卡尔坐标系下放置，为消除空心球-金属基三相复合材料的边界效应，其仿真模型任意方向空心球数量和单一胞元数量应满足：

N_dmin＝N_min-1 (6)

其中，N_dmin为消除边界效应时目标模型任意方向单一胞元数量最小值；N_min为消除边界效应时目标模型任意方向空心球数量最小值；

则该模型几何参数为：

L_x＝m·a(N_dmin≤m≤m_max) (7)

L_y＝n·a(N_dmin≤n≤n_max) (8)

L_z＝p·a(N_dmin≤p≤p_max) (9)

其中，L_x为目标模型在x方向长度，单位mm；L_y为目标模型在y方向长度，单位mm；L_z为目标模型在z方向长度，单位mm；m为目标模型x方向单一立方体胞元的个数；n为目标模型y方向单一立方体胞元个数；p为目标模型z方向单一立方体胞元个数；m_max为目标模型x方向单一立方体胞元的个数最大值；n_max为目标模型y方向单一立方体胞元的个数最大值；p_max为目标模型z方向单一立方体胞元的个数最大值；

即建立单一胞元结构与目标空心球-金属基三相复合材料结构的数量和几何映射关系；

第二步，根据目标空心球-金属基三相复合材料结构的外形尺寸、空心球-金属基三相复合材料的球基分布规律和球基间拓扑结构关系创建目标结构几何模型，具体如下：

(1)在笛卡尔坐标系下，以(0，0，0)、(L_x，0，0)、(0，L_y，0)、(0，0，L_z)为基础点创建目标样件几何实体；

(2)在y＝0的X-O-Z平面创建第一层空心球坐标点，即：

P_i,j,k＝(ia,ja,ka)(i＝0,1,···,m,k＝0,1,···,p,j＝0) (10)

(3)在

的X-O-Z平面创建第二层空心球坐标点，即：

(4)将两层坐标点沿Y轴正方向线性复制阵列，阵列数量为n，阵列间距为a，单位mm；再将第一层坐标点沿Y轴正方向线性复制阵列，阵列数量为2，阵列节距为na，单位mm，最终获得空心球球心位置坐标点；

(5)在空心球球心坐标点处创建实心球体，球体半径为空心球内部半径r：

r＝R-t (12)

通过布尔减法运算去除样件实体中在球心坐标点处创建的半径为r的球形实体，获得多孔实体模型；

(6)在(4)中生成的球心位置坐标处创建空心球外壁曲面，曲面半径与空心球外部半径尺寸相同；通过曲面切割实体，利用空心球外壁曲面切割(5)中获得的多孔实体模型；

通过切割得到壁厚为(R-r)的空心球和多孔基体，两者间曲面的拓扑状态为共享曲面，以共享曲面为目标面创建2D网格，得到空心球-金属基三相复合材料结构的几何模型；

第三步，为空心球-金属基三相复合材料的离散单元赋予材料属性和节点约束条件：

(1)将空心球-金属基三相复合材料结构几何模型离散为各种单元，并在离散的单元节点上进行位移插值；

将空心球-金属基三相复合材料结构几何模型各部分节点的位置与时间的变化关系用离散单元节点的坐标插值表示为：

其中，i＝1,2,3为空间坐标指标；I＝1,2,···,n表示有限单元节点；其中μ_iI表示单元第I个节点在i方向的位移；N_I(x,y,z)为第I个节点的拉格朗日插值函数；

对空心球-金属基三相复合材料进行三角形二维单元面网格处理，使面网格与复杂的几何表面属性一致，且满足：

p≤t (14)

其中，p为二维网格尺寸；t为空心球壁厚；

选择十节点四面体单元离散第二步中生成的空心球-金属基三相复合材料几何模型，在结构紧凑区域自动生成较小的网格单元，最大程度保持网格表面与第二步中生成的目标结构几何模型的一致性，在离散实体时确保四面体网格与二维网格具有相关性，且满足：

p′≤t (15)

其中，p’为十节点四面体单元尺寸；t为空心球壁厚；

经过离散后，得到空心球-金属基三相复合材料离散单元，其中包括以空心球几何结构离散得到的空心球离散单元和以多孔基体几何结构离散得到的固态基体离散单元；

(2)为空心球-金属基三相复合材料离散单元赋予材料属性和节点约束条件；具体如下：

(a)为空心球离散单元和固态基体离散单元设置屈服准则、材料参数和失效准则；

设空心球、基体材料均为各向同性理想弹塑性金属材料，分别为空心球离散单元和固态基体离散单元定义材料的杨氏模量、屈服强度、质量密度、泊松比、切线模量、失效参数和等效应力-应变曲线；

空心球和固态基体的离散单元均以材料失效时的有效塑性应变和最小时间步定义材料失效指标，当离散单元等效塑性应变达到失效值时，离散单元即发生失效；

设置离散单元属性和截面属性，创建十节点四面体实体单元属性并赋予给空心球、固态基体离散单元；

(b)为空心球-金属基三相复合材料离散单元节点赋予节点约束条件，即在空心球-金属基三相复合材料离散单元节点处建立接触约束；

空心球与固态基体离散单元节点间的接触摩擦系数由多组参数共同决定：

其中，μ_c为接触摩擦系数；FD为动摩擦因数；FS为静摩擦因数；DC为指数衰减系数；v_rel为接触面相对速度；

同时鉴于空心球和基体组成相的失效模式，在仿真两者单元节点间接触关系时需要定义带有失效准则的绑定接触；当达到接触失效准则后绑定接触失效，空心球和基体单元节点间的接触类型转换为普通接触；

当接触面单元间达到界面极限分离距离后，接触点完全破坏，空心球与基体间界面完全分离，两者间脆性金属间相完全损伤失效并形成球壁外侧微裂纹；

通过以上三步获得空心球-金属基三相复合材料细观尺度精细化仿真模型；计算目标工况下材料的力学响应和损伤模式，建立材料的仿真模型裂纹观测截面，通过观测截面裂纹分布，在细观尺度下得到空心球-金属基三相复合材料的界面失效和裂纹形式。

2.根据权利要求1所述的一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，其特征在于，根据权利要求1所述的所述步骤(2)(a)中，为空心球离散单元和固态基体离散单元均赋予各向同性Von Mises屈服准则，主坐标系下的Von Mises屈服准则为：

其中，σ₁、σ₂、σ₃为笛卡尔坐标系下X、Y、Z方向的单元主应力，σ_s为单元屈服应力；其物理意义为，当空心球、固态基体离散单元体积形状改变的弹性模量达到某一常数时，该单元质点发生屈服。

3.根据权利要求1或2所述的一种空心球-金属基三相复合材料细观尺度精细化建模仿真方法，其特征在于，所述步骤(2)(b)中，接触失效准则为：

其中，σ_n和σ_s分别为接触面正应力和接触面切应力，单位MPa；NFLS与SFLS分别为对应的接触面失效正应力和接触面失效剪切应力。

Images