《《3D打印金属植入物有限元分析方法》》

ICS11.040.40

中国医疗器械行业协会团体标准

T/CAMDI045—2020

3D打印金属植入物有限元分析方法

FiniteElementAnalysisof3DPrintedMetalImplants

2020-06-18发布2020-07-01实施

中国医疗器械行业协会发布

T/CAMDI045--2020

3D打印金属植入物有限元分析方法

1范围

本标准规定了使用有限元分析（FEA）技术对3D打印金属植入物进行数值模拟分析的要求、方

法、报告内容等。

本标准适用于对3D打印金属植入物在产品设计阶段和验证阶段围绕产品力学性能、运动功能和

生物功能所开展的相关终产品的结构性能进行的有限元分析。

依据本标准在临床或实验等效条件下对3D打印金属植入物的性能评价，可作为产品设计阶段和

验证阶段的参考依据。如获得产品真实可靠性的最终评价依据，应添加相应的体外实验、动物实验

等资料。

2规范性引用文件

下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本

文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。

定制式医疗器械监督管理规定（试行）

GB/T31054机械产品计算机辅助工程有限元数值计算术语

GB/T33582机械产品结构有限元力学分析通用规则

GB/T35022增材制造主要特性和测试方法零件和粉末原材料

GB/T36984外科植入物用多孔金属材料X射线CT检测方法

T/CAMDI025定制式医疗器械力学等效模型

ASTMF2996非组合式金属髋关节股骨柄有限元分析标准方法

3术语和定义

GB/T33582机械产品结构有限元力学分析通用规则中的术语与定义同样适用于本标准。

3.1有限单元法finiteelementmethod

将连续的求解域离散为有限个单元，并在给定边界条件下，利用有限单元的近似解逼近真实物

理系统的数值分析方法。

3.2有限元分析finiteelementanalysis

基于有限单元法的结构性能分析。

3.3有限元模型finiteelementmodel

以有限数量单元的组合逼近产品结构的模型。

3.4有限元建模finiteelementmodeling

构建有限元模型的过程，包括几何模型构建和（或）处理、材料属性定义、网格划分、边界条

件施加等步骤。

1

T/CAMDI045--2020

3.5单元element

具有几何、物理属性的最小求解域。

3.6节点node

单元之间的铰接点。每个单元仅在节点处和相邻单元及外部发生联系。

3.7边界条件boundaryconditions

在给定工况下，求解域边界上的几何、物理条件。

3.8材料属性materialproperty

材料的物理性能参数，如弹性模量、泊松比、密度等。

3.9各向同性isotropy

是指物体的物理、化学性质不因方向而有所变化的特性，即在不同方向所测得的性能数值是相同

的。

3.10各向异性anisotropy

物体的全部或部分物理、化学等性质随方向的不同而有所变化的特性。

3.11敏感性分析sensitivityanalysis

在存在不确定性的情况下，通过改变某一变量测试该变量对模型或系统结果的稳健性影响。

3.12应力集中stressconcentration

结构局部区域，刚度或形状急剧变化引起的应力数值明显增高的现象。

4要求

4.1基本要求

对于建立的有限元模型，应考虑风险最高的产品几何和实际使用情况。建议采用SI国际单位作为

有限元仿真的标准单位。对有限元模型的几何、材料边界条件等的简化应不影响植入物性能的评价。

4.2模型数据

被分析产品的几何外形和细节特征的设计过程充分考虑临床医师和工程师的意见。

4.2.1宏观结构

有限元模型的宏观结构应符合被分析结构的几何外形，其几何细节可以从图纸、几何模型、或任

何其它与所定义模型几何形状一致的资料获取。有限元模型输入几何结构的简化应不影响植入物性能

的评价。

4.2.2多孔微结构

3D打印金属植入物根据其功能不同，可能存在不同种类的多孔微结构。多孔微结构通常有减重或

骨长入的作用。在对3D打印金属植入物进行有限元分析的过程中应对这一特征进行充分考虑。在计算

能力许可的条件下，几何模型应能充分反映3D打印形成的孔隙结构。

2

T/CAMDI045--2020

4.3网格模型

4.3.1网格单元、尺寸的确定应进行网格敏感性分析

4.3.2非均一化的网格划分，应对植入物性能评价关注区域的细小特征、孔结构、接触面等做网格

细化处理。

4.4材料属性

4.4.1各向同性材料

为了计算有限元模型的应力和应变，至少需要输入至有限元模型中的所需材料属性为弹性模量

（E）和泊松比（μ）。

4.4.2各向异性材料

在无法证明材料为各向同性的情况下，在进行有限元分析前应充分评估材料在各个方向的力学

属性，3D打印金属植入物应考虑按照各向异性材料进行属性赋值，并对各向异性材料属性对结果分

析的影响进行评估。

4.5边界条件

边界条件的施加应充分考虑被分析结构在人体内的约束环境，使有限元分析模型能充分反映植

入物的约束状况。必要时边界条件的输入条件应同临床医师进行充分沟通。

加载条件应与被分析结构在体情况相符。加载条件分为载荷控制和位移控制，应根据实际情况

合理选择加载方式。

4.6模型计算

将上述几何模型、材料参数、边界条件等设置完成，经离散化处理，从而建立完整的有限元分

析模型。根据模型可能的失效模式，应合理选择强度理论和分析类型，并对选择的合理性做出说明。

4.7模型验证

应采用有效手段对有限元模型的分析方法、输入参数和计算结果的有效性进行验证。

4.8软件

对3D打印金属植入物进行有限元分析的软件，应包含完整的静力学分析与动力学分析功能或模

块。使用行业中常用并公认的有限元软件进行分析时，应提供软件名称和版本号。当使用其他软件

时需提供软件分析可靠性的说明。

注：可通过使用不同软件对相同算例进行运算结果的对比验证软件的可靠性。可靠性说明应单

独以验证报告的形式提供。

5方法

5.1模型数据

核查导入有限元分析软件的几何模型数据，几何模型数据应满足5.2的要求，尺寸精度满足植入

物性能评价要求，并说明合理性。

仅包含植入物部件的有限元分析，其几何模型从产品的CAD几何文件中获取。

包含植入物以外多个其他结构的有限元分析，如模拟体外试验装置或模拟人体预期生理环境，

需要对其他部件几何模型数据的来源及合理性进行说明，并明确植入物及部件间的相对位置关系。

3

T/CAMDI045--2020

在构建植入物有限元模型时，可以对几何模型的细小特征，如：小倒圆、倒边特征、非重要区域的

薄、锐特征和难以网格划分的多孔结构等进行简化，并对简化部位及简化理由做出合理性说明。

5.2网格模型

网格模型的建立应满足5.3的要求。

网格划分推荐优先选用六面体单元，其次选四面体单元，优先选用高阶，网格单元大小应根据网

格敏感性分析确定，满足差异阈值为5%以内的敏感性要求。

在网格尺寸上建议由大到小多次划分，对比计算结果。结果差异小于5%时，认定网格收敛。在需

要网格细化的部位，应选取适合的细化准则，并标明单元细化结果。可以通过自动划分网格和手动划

分网格、或是二者结合的方法对模型进行有限元网格划分。最重要的考虑因素在于，所使用单元的类

型、尺寸和形状必须能代表预期的行为，且没有重大的数值局限性。对于单元划分的质量，大多数有

限元分析软件具有检查单元形状的内置程序工具。若无此工具，则需要进行额外的检查。

存在局部网格细化的情况，指出其具体位置和局部网格尺寸。

5.3材料属性

材料属性应满足5.4的要求。各向同性材料获取材料参数的方法，可参照经验或文献获取，或按照

相应标准测试方法进行试验。根据分析目的做线弹性分析时须至少获取如下参数：材料的弹性模量（E）、

泊松比（μ）等线弹性段参数；对于材料非线性分析时还应该获取应力-应变曲线或者屈服拉伸强度

（YTS）、极限拉伸强度（UTS）、断裂伸长率（A），断面收缩率（Z）等非线性段参数。各向异性材

料在不同方向按照上述方式分别获取。

其中，应变可用引伸计、应变片、光学全场测量等方式获取。注意因应变获取方式不同而带来的，

工程应力-应变和真实应力-应变之间的转化。

5.3.1各向同性的或实体部分金属，材料属性可以通过材料检验报告获取。判断材料是否为各向同性

材料属性的方法：采用在X\Y\Z以及45度摆放打印的试棒进行拉伸测试，当结果差异小于15%时，认为

材料为各向同性材料，反之则各向异性。

5.3.2当材料属性不明确时，应采用相同工艺下打印的试样进行相应的力学测试获得材料属性，并

进行留样。提供的留样，应能充分反映材料的拉伸、压缩等力学性能。

注1：可参照GB/T35022中推荐的试验方法对增材制造样件进行机械性能检测。

注2：在对各向异性材料进行赋值时，建议应首先考虑将材料的各向异性，并按照实测材料参数进行

赋予；或可使用各向异性材料中性能最弱的方向的材料参数作为整个材料的属性赋值。材料属性的选

取和赋值过程需在报告中完整论述。

5.3.3对于多孔结构，建议按照其自身几何特征进行有限元网格划分和简化。如果多孔结构可以简化

为具有代表性结构的单元体进行阵列填充方式形成的多孔材料，可将多孔部分按照其宏观几何外轮廓

建立模型，其材料属性通过对代表性结构单元体进行力学实验获得，并按照实体材料进行属性赋值。

对于梯度孔径或随机排布的多孔结构，尽量不进行简化处理。

如果多孔结构在整个植入物中占比极小，且对整体力学性能无影响（例如骨整合面多孔层）。可

将这种多孔结构简化为实体，材料属性按照实体金属的百分比赋予。同时分析人员需提供该简化过程

的合理性说明。

4

T/CAMDI045--2020

注：对于随机生成的多孔结构，应先根据GB/T36984中的方法先确定孔隙分布的均匀性。

5.3.4对于有限元模型中涉及的其他材料应尽量反映其真实材料属性。

5.4边界条件

边界条件的设置应满足5.5的要求。

5.4.1在医工交互的条件下，根据临床实际情况，对有限元模型应施加符合在体情况的位移、约束、

力或其他必要的边界条件。明确有限元模型载荷、位移的施加方式、作用位置、大小和方向，模拟

植入物在体的功能状况。

5.4.2边界条件的施加应充分考虑模型中是否存在与其他部件或自身发生表面接触的区域。对该类

型区域施加接触边界条件。

5.4.3对于设计上存在对称性结构的产品，如果其边界条件也是对称的，则可以使用部分对称结构，

对其施加对称边界条件，进行简化计算。否则不得做简化处理。

5.4.4边界条件的设置应注意力学等效模型的建立，可以参考T/CAMDI025-2019定制式医疗器械力

学等效模型中的相关要求进行，明确各部分的装配关系，并确认其合理性。

注：当3D打印金属植入物的类型和受力状态与现有标准或方法类似时，可参考使用类似产品的测试

标准。有统计学意义的数据也可以作为加载载荷的条件。

5.5模型计算

模型的计算应满足5.6的要求，并定义及关注如下要素。

5.5.1定义分析步、分析类型和输出参数。

5.5.2应对有限元分析的计算过程进行检查，确保几何模型、网格模型、材料参数和边界条件等在

分析模型中被正确合理的定义，并能正确代表所分析的情况。确保有限元分析结果无单元应力应变

不协调的情况。

5.6结果输出

计算结果通过图表等形式呈现，并对结果进行文字描述、解释和分析总结。通过云图呈现结果，

须包括图例、最大值应力、最大位移值及最大值发生的位置。要有必要的局部放大图以及关键截面

图，多参数结果比较时采用统一图例。明确输出参数与植入物性能评价指标的关系。

注：可参照附录A中实例的样式进行结果输出。

5.7模型验证

5.7.1实验方法。应采用有效手段对有限元模型的有效性进行验证，有限元与物理测试应尽可能体

现实际测试的工况，可以进行适当简化。将实验结果同有限元结果进行对比，通过两者结果的一致

程度，对模型的有效性进行验证。对于力学实验建议采用引伸计、应变片或非接触式光学测量的方

式进行数据测量。

5.7.2文献对比。在实验很难复现有限元模拟的情况下，建议通过采用文献对比的方式，通过查阅

相关研究文献中对于类似产品的研究结果，作为支持性数据，证明有限元模型的有效性。

5

T/CAMDI045--2020

6报告

应在报告中全面的记录3D打印金属植入物的有限元分析过程。报告格式应符合可接受的工程报告

稿的格式要求。报告应包含但不限于以下信息：

框架内容

a)有限元分析目标或需求，包含分析对象的完整描述，包括分析目的、产品介绍、分析方法等；

b)分析软件，使用的商用软件应提供软件的名称及版本号。非商用软件提供软件分析可行性说明。

使用软件子程序及特殊模块等，需提供验证分析方法的详细信息；

c)参考标准，包含力学实验和有限元仿真工况中涉及的全部参考标准；

d)几何结构，包含模型的几何数据，包含外观尺寸、多孔结构等，对几何模型的简化情况及其合理

性描述；

e)网格模型，包含网格的类型、大小及网格数量，计算网格收敛性的方法和过程；

f)材料参数，提供材料参数的数据来源，当材料性能参数来源于力学实验时，应提供相应的力学实

验报告作为附件材料；

g)边界条件，对于模型中使用的边界条件应提供数据来源的文献或标准；

h)输出结果，应给出典型的图表结果，如应力/应变云图、数据表等，图表应简明易懂，不应有无

关信息；

i)模型验证，提供参照的文献或进行的力学实验过程，模型验证过程中的对比数据分析结果，模型

有效性的分析过程与结果；

j)结论与建议

7附录

附录A建议的报告样式

附录B全实体结构的典型产品的操作实例

附录C含有多孔结构的典型产品的操作实例

6

T/CAMDI045--2020

附录A

（规范性）

3D打印人工颞下颌关节有限元分析报告

A.1有限元分析目标

颞下颌关节（TMJ）是颅面复合体的重要组成部分，对于言语，咀嚼，吞咽和表达情绪等活动必

不可少。涉及TMJ的疾病是比较普遍的，其疾病常可累及一个或多个附近部位，主要表现为关节疼

痛(开口受限及咀嚼能力受影响）、严重者可见不同程度的面部畸形，对患者的正常生活、学习等带

来不利影响。16％‐59％的患病率报告的症状和33‐86％的临床体征表示,修复TMJ的颞下颌关节置换

（TJR）手术已被广泛接受为治疗TMJ病情严重的可靠方法。当保守和非手术疗法无效时，其中包括

晚期退行性关节疾病，肿瘤，发育异常，强直，外伤和骨折等疾病，采用颞下颌关节置换（TJR）手

术是一种有效且可靠的关节重建方法。

随着3D打印技术和医学影像技术的迅速发展，通过三维重建和3D打印技术制造的定制化关节

假体更符合患者的解剖结构。使用仿真力学分析软件分析设计完成后的关节假体在受力情况下的应

力应变等值以及解决应力集中等问题，可以为颞下颌关节修复重建的设计提供更可靠的技术保证。

本次有限元分析采用的是静力学分析，主要目的是对根据特定病例定制的个性化颞下颌人工关

节进行分析验证，并对比两种不同髁突设计构型。

A.2分析软件

采用ANSYS17.0软件进行仿真计算

A.3参考标准

GB∕T13810-2017外科植入物用钛及钛合金加工材

A.4几何结构

模型名称几何尺寸（单位：mm）是否简化及说明

上颌骨148.44x141.69x105.33是

下颌骨123.47x98.87x84.00是

关节窝17.15x23.74x17.98否

仿生关节柄15.59x27.38x55.59否

人工关节柄16.70x27.03x55.78否

螺钉7.39x3.85x5.23是。为了减少计算量，螺钉未设计螺纹

关节盘25.98x16.34x11.85是

A.5网格模型

部件名称网格类型网格数节点

上颌骨Quadratictetrahedral79191130941

下颌骨Quadratictetrahedral4083263417

关节窝假体Lineartetrahedral271487957

仿生关节柄Lineartetrahedral3823410932

7

T/CAMDI045--2020

人工关节柄Lineartetrahedral337649590

钛钉Quadratictetrahedral8511624

关节盘Quadratictetrahedral25664634

A.6材料参数

材料类型杨氏模量（GPa）松泊比

上颌骨/下颌骨Bone-C150.3

关节窝假体Polyethylene0.80.4

仿生关节柄Ti-6AI-4V1100.3

人工关节柄Ti-6AI-4V1100.3

钛钉Ti-6AI-4V1100.3

关节盘DiskAnterior44.1MPa0.4

参考文献：

HuiGao,XianLi,ChunjuanWang,PingJi,ChaoWang,Mechanobiologicallyoptimizationofa3D

titanium-meshimplantformandibularlargedefect:Asimulatedstudy,MaterialsScience&

EngineeringC.104(2019)109934.

A.7边界条件

本次分析给定的载荷是咬合力，具体值见下表：

最大咬合力（N）最大咬合力（N）

颞前肌222.3

颞后肌158.6

左侧咬肌浅层196.3

咬肌中层37.9

咬肌深层44.9

翼内肌170.8

翼外肌上头65.5

翼外肌下头5.1

固定部位：

8

T/CAMDI045--2020

接触关系：

接触物体接触关系

上颌骨关节盘Frictionless

上颌骨钛钉Bonded

上颌骨下颌骨No-Separation

上颌骨人工关节窝No-Separation

下颌骨关节柄Bonded

下颌骨钛钉Bonded

下颌骨关节盘Frictionless

关节柄钛钉No-Separation

关节窝钛钉No-Separation

关节柄关节窝Frictional

参考文献：

[1]Y.W.Lee,H.J.You,J.A.Jung,D.W.Kim,Mandibularreconstructionusingcustomized

three-dimensionaltitaniumimplant,Arch.Craniofac.Surg.19(2018)152–156.

[2]HuiGao,XianLi,ChunjuanWang,PingJi,ChaoWang,Mechanobiologicallyoptimizationofa3D

titanium-meshimplantformandibularlargedefect:Asimulatedstudy,MaterialsScience&EngineeringC.

104(2019)109934.

[3]DavidAcklanda,DaleRobinsona,PeterVeeSinLeea,GeorgeDimitroulis,Designandclinicaloutcome

ofanovel3D-printedprostheticjointreplacementforthehumantemporomandibularjoint,Clinical

Biomechanics56(2018)52–60.

A.8输出结果

A.9网格敏感性验证

9

T/CAMDI045--2020

有限元网格尺寸设置在1mm左右最为合理，可以保证足够的计算精度，又不至于消耗过大的计

算时间。收敛性如下图所示：

A.10结论与建议

从有限元分析结果可以得出，人工髁突的接触应力大于仿生髁突设计颞下颌关节，因此在设计定

制化关节假体时，应优先选择仿生髁突头。

10

T/CAMDI045--2020

附录B

（资料性）

股骨柄有限元分析及实验验证

B.1分析目标

使用本标准的方法，对定制式3D打印金属股骨柄植入产品进行静态和动态力学性能分析，验证

所设计的结构和3D打印制造的产品是否可靠，为临床应用提供设计和制造依据。

B.2分析软件

软件名称：ANSYS版本号：18.2

B.3参考标准

YY/Y0809.4-2018,ASTMF2996-13,GB/T13810-2017

B.4几何结构

图1为特征尺寸，股骨柄假体长度（CT）=150.53mm。

图2根据ASTMF2996-13对模型进行处理，红色为固定约束位置，绿色区域为过渡区域，避

免由于约束区域刚度过大产生的判断误差。

图3为有限元分析整体模型。

图1图2图3

B.5网格模型

（1）网格类型：四面体网格

（2）网格大小：1mm

（3）网格数量：223276

B.6材料参数

股骨柄球头/压块

材料TC4钴铬钼

抗拉强度1100MPa1350MPa

屈服强度1000MPa1050MPa

11

T/CAMDI045--2020

弹性模量110000MPa210000MPa

泊松比0.30.3

TC4S-N曲线（图4）：

图4TC4疲劳曲线图

B.7边界条件

静态

载荷：2300N（图5标注B）

约束：按ASTMF2996-13方式对CT=90mm[2]以下区域进行全约束（图5标注A）；约束压块X、Z

方向的位移为0mm，Y方向free（图5标注C）

疲劳

载荷：恒定载荷（Ratio），载荷比𝑅2300.1

2300

图5边界条件

B.8输出结果

B.8.1静态分析

12

T/CAMDI045--2020

静态分析

位移应力

结果：分析获得最大位移0.150mm。结果：分析获得颈部最大应力为136.31MPa。

应力随位置变化规律

结果：从位置①-位置②应力变化曲线，最大应结果：从位置①到位置②的应力变化曲线，最大应

力为92.24MPa出现在颈部。力为135.58MPa出现在颈部。

13

T/CAMDI045--2020

B.8.2疲劳分析

疲劳应力

结果：颈部最大交变应力为76.67MPa。

安全系数疲劳寿命

结果：模拟分析获得该模型尺寸植入件的最小安

结果：模拟分析获得的疲劳寿命能达到1000万次。

全系数为8.6。

B.9模型验证

B.9.1股骨柄实验验证

采用SLM工艺打印出与模拟分析模型一致的钛合金股骨柄植入件，打印后的植入件随基板进行

整体热处理，打印工艺参数及热处理规范如表1所示，完成热处理的植入件经打磨去除表面台阶和粗

糙纹路，获得光滑表面。

表1钛合金股骨柄打印工艺参数与热处理规范

名称材料激光功率扫描速度扫描间距层厚

股骨柄TC4ELI280W1200mm/s0.14mm30

热处理规范：处理温度800℃；保温2h；空冷（气保护炉进行）

股骨柄植入件按YY/T0809.4-2018方法要求进行疲劳测试，加载载荷2300N，循环疲劳次数500

14

T/CAMDI045--2020

万次，试验装置如图6所示；股骨柄静态载荷位移数据如图7所示，2300N时位移变化为0.162mm，

有限元模拟结果为0.150mm，吻合度为92.59%；疲劳试验数据采集如图8所示，500次循环过程中系统

没有出现失稳显现，试验完成后取出植入件进行检测，无任何缺陷出现。

图6股骨柄植入件疲劳试验装置图7股骨柄静态载荷位移数据

图8疲劳试验数据采集

B.9.2力学性能验证

使用与植入件同炉打印的试样，并经过相同的热处理后进行力学性能检测和孔隙率测定，获得

的力学性能和孔隙率如表2及图9、图10所示。

表2随炉样件力学性能

抗拉强度（MPa）屈服强度（MPa）延伸率（%）

TC4ELI1100±501000±5015±2

图9TC4样件的拉伸曲线图10TC4样件的孔隙度

15

T/CAMDI045--2020

随炉钛合金试样经截面法金相检测孔隙率为0.006%，致密度达到99.96%，最大孔洞尺寸小于20.14

μm；抗拉强度1100±50MPa，屈服强度1000±50MPa，延伸率15±2%，满足GB/T13810-2017的

要求。

B.10结论与建议

有限元分析可获得优化的几何模型和应力分布状况，为临床应用植入件提供了设计和制造的依据，

以目前的SLM工艺水平，采用较佳工艺参数生产的定制式/个性化植入产品，其力学性能能够满足临

床植入的要求，鉴于植入件安全性能考虑，需要进行同炉试样的相关性能验证，至少包括力学性能和

孔隙率检测。

参考文献：

[1]YY/T0809.4-2018外科植入物部分和全髋关节假体第四部分：带柄股骨部件疲劳性能试验和性能要

求.[S]

[2]ASTMF2996-13StandardPracticeforFiniteElementAnalysis(FEA)ofNon-ModularMetallic

OrthopaedicHipFemoralStems.[S]

[2]GB/T13810-2107外科植入物用钛及钛合金加工材.[S]

[3]ChethanK.N.MohammadZuber，Staticstructuralanalysisofdifferentstemdesignsusedintotalhip

arthroplastyusingfiniteelementmethod.

[4]傅彦棉，彭伟，盛浩，等.基于ABAQUS及疲劳试验的烧结股骨柄设计研究[J],2019,16（3）：71-73.

16

T/CAMDI045--2020

附录C

（资料性）

3D打印金属椎间融合器有限元分析及实验验证

C.1软件信息

软件名称和版本：MagicsV23（打印前处理），ABAQUS2017（有限元前处理和计算），

SOLIDWORKS2014（几何模型建立），Mimics21（CT扫描逆向重建），Hypermesh2017（有限元前

处理）。

C.23D打印实体标准件对标

C.2.13D打印实体件产品测试

实体拉伸标准件拉伸测试，按照《GB∕T228.1-2010金属材料拉伸试验第1部分室温试验方法》

进行设计，并采用引伸计采集应变。获取试验的应力应变曲线。下面以1号样品为例作为展示（图1），

实验照片和试验曲线如图2所示。

图1样品设计

60000

50000

40000

30000

载荷(N)

20000

10000

0

-1012345

拉伸位移(mm)

图2实验照片及样品实验曲线

表1实验测得的材料属性数据

弹性模量(自动)屈服强度(偏移量0.002mm/mm)抗拉强度

(MPa)(MPa)(MPa)

1111,327.609441,029.591111,087.72729

17

T/CAMDI045--2020

图3试验后样品

C.2.23D打印实心件产品有限元分析

拉伸测试的有限元模型如图4（a）所示，材料属性设置如图4（b）所示，分为弹性段和塑性段。

弹性阶段输入材料的杨氏模量和泊松比，塑性阶段输入参数第一列为应力，第二列为塑性应变。

（a）（b）

图4拉伸测试有限元模型：（a）有限元模型；（b）材料属性

有限元模型的边界条件设置如图5所示。移动夹持端，使得标距段产生4mm的位移。

图5有限元模型的边界条件设置

计算后输出结果中的应力云图和标距段的力‐位移曲线，如图6和图7所示：

18

T/CAMDI045--2020

图6结果应力云图

图7标距段力‐位移曲线

实验后样品形貌同有限元模型分析后样品形貌进行对比，对比如图8所示：

图8实验样品同有限元模型的形貌对比

实验所得的力‐位移曲线同有限元分析的结果进行对比，如图9所示。可以看出，2条曲线在弹

性区、塑性区均非常重合，说明在该区域对标成功。不考虑破坏区，因此未进行对标。

19

T/CAMDI045--2020

图9实验所得的力‐位移曲线同有限元分析的对比结果

C.33D打印多孔件产品对标

C.3.13D打印多孔件产品测试

多孔随炉样块样品设计及打印的样品如图10所示，控制尺寸为10*10*10。

图10多孔随炉样块样品设计及打印的样品

将多孔样块进行研所实验，所得的实验曲线和结果数据如图11所示。

20

T/CAMDI045--2020

1000

0

-1000

-2000

-3000

载荷(N)载荷

-4000

-5000

-6000

-1012345

压缩位移(mm)

最大压缩载荷弹性模量(自动杨氏)抗压强度

(N)(MPa)(MPa)

15,641.547361,807.1579757.56681

25,822.810061,545.5310958.22810

平均值5,732.178711,676.3445357.89745

图11多孔样块实验数据：斜率和峰值载荷基本一致，斜率和峰值载荷基本一致，说明样品重复

性好，实验误差小。

C.3.23D打印多孔件产品有限元分析

建立多孔样件的有限元模型，采用梁单元B31,Truss单元没有旋转自由度只能承受轴向载荷，

BEAM单元可以承受弯矩，更符合实际。几何模型如图12所示。丝径横截面积为圆形，半径0.12mm。

图12多孔样件的几何模型

边界条件，如图13所示。底部完全固定，在顶端施加垂直载荷，直至失效。

21

T/CAMDI045--2020

图13边界条件设定

对不同单元尺寸进行收敛性分析分析单元尺寸表2所示，结果如图14所示。

表2梁单元不同单元大小对应的模型数量对比表格：

梁单元长度单元数量节点数量最大载荷

0.8mm7772332112252.9

0.4mm1432498739521.1

0.2mm29498250478117.3

0.1mm60494560437639.8

图14梁单元不同长度的载荷位移结果对比。

22

T/CAMDI045--2020

可以看出，梁单元长度从0.2mm降至0.1mm时，最大载荷变化约5%，符合收敛性要求，因此

后续模型采用0.2mm的梁单元进行划分。

实验结果曲线同有限元分析结果进行对比，对比结果如图15所示：

图15实验结果曲线同有限元分析结果

可以看出，半径0.12mm，划分0.2mm的梁单元，其计算结果与实验曲线的最大值相近。斜率比

试验的略大，但比较接近。此处斜率的差异可以解释为测试样品与有限元模型的细节差异导致，结

果如表3所示。

表3实验结果

最大压缩载荷杨氏模量

(N)（MPa）

R=0.125217.3825754.36

TEST5641.5519286.49

C.43D打印