第5章：材料科学与打印参数

本章深入探讨3D打印材料的物理化学特性，建立材料行为与打印参数之间的定量关系。通过流变学、热力学和力学模型，您将学会预测和优化打印过程中的材料响应，从而实现高质量、高可靠性的打印结果。

5.1 聚合物流变学基础

5.1.1 熔体流动的本构方程

FDM打印的核心是控制熔融聚合物通过喷嘴的流动。聚合物熔体表现出非牛顿流体特性，其剪切应力$\tau$与剪切速率$\dot{\gamma}$的关系遵循幂律模型：

$$\tau = K \dot{\gamma}^n$$ 其中$K$是稠度系数（Pa·s^n），$n$是流动指数（无量纲）。对于大多数热塑性材料，$n < 1$表现为剪切变稀（pseudoplastic）行为。典型值：PLA的$n \approx 0.3-0.4$，ABS的$n \approx 0.35-0.45$。这种剪切变稀特性对打印有利，因为高剪切速率下粘度降低，减小了挤出阻力。

在圆形喷嘴中，通过积分Navier-Stokes方程，可得体积流率$Q$与压降$\Delta P$的关系（Rabinowitsch方程）： $$Q = \frac{\pi R^3}{2K^{1/n}} \left(\frac{n}{3n+1}\right) \left(\frac{\Delta P}{2L}\right)^{1/n}$$ 其中$R$是喷嘴半径，$L$是喷嘴长度。值得注意的是，流率与压降呈非线性关系，指数为$1/n > 1$，意味着压力的小幅增加可显著提高流率。

壁面剪切速率可估算为： $$\dot{\gamma}_w = \frac{4Q}{\pi R^3} \cdot \frac{3n+1}{4n}$$ 典型打印条件下（$Q = 10$ mm³/s，$R = 0.2$ mm），剪切速率可达$10^3-10^4$ s⁻¹，远高于传统成型工艺。

5.1.2 温度依赖性：WLF方程

聚合物粘度强烈依赖于温度，可用Williams-Landel-Ferry (WLF)方程描述： $$\log a_T = \log\frac{\eta(T)}{\eta(T_r)} = -\frac{C_1(T-T_r)}{C_2+(T-T_r)}$$ 其中$a_T$是时温叠加因子，$T_r$是参考温度（通常取$T_g + 50K$），$C_1$和$C_2$是材料常数。对于大多数无定形聚合物，当$T_r = T_g + 50K$时，普适常数为$C_1 \approx 17.44$，$C_2 \approx 51.6K$。

WLF方程在$T_g$到$T_g + 100K$范围内有效，超出此范围应使用Arrhenius方程： $$\eta(T) = \eta_0 \exp\left(\frac{E_a}{RT}\right)$$ 其中$E_a$是流动活化能。实际应用中，温度每升高10°C，粘度约降低30-50%，这解释了为什么温度控制对打印质量至关重要。

5.1.3 挤出胀大现象

聚合物离开喷嘴后会发生弹性恢复（die swell），这是由于流动过程中储存的弹性能释放。挤出物直径$D_e$与喷嘴直径$D_n$的比值称为胀大比： $$B = \frac{D_e}{D_n} = 1 + \frac{1}{2}\left(\frac{\tau_w}{\sigma_y}\right)^2 + \alpha \cdot De$$ 其中$\tau_w$是壁面剪切应力，$\sigma_y$是材料屈服应力，$De = \lambda \dot{\gamma}$是Deborah数（$\lambda$为松弛时间）。典型的胀大比$B$在1.1-1.3之间，具体取决于：

• 剪切速率：高速挤出增大胀大
• 温度：高温降低弹性，减小胀大
• 分子量：高分子量增大胀大
• L/D比：短喷嘴胀大更明显

实际打印中，通过调整挤出倍率（extrusion multiplier）补偿胀大效应，典型值0.9-1.1。

5.1.4 压力驱动流动模型

在Bowden挤出机中，丝材通过PTFE管输送至热端，驱动压力由挤出齿轮产生： $$P_{drive} = \frac{F_{gear}}{A_{filament}} = \frac{\mu F_n}{A_{filament}}$$ 其中$F_{gear}$是齿轮推力，$\mu$是齿轮-丝材摩擦系数（钢齿轮对PLA约0.3-0.4），$F_n$是弹簧预紧力（典型20-50N），$A_{filament} = \pi d^2/4$是丝材截面积。

系统总压降由多个组分构成：

1. 热端入口收缩损失（Bernoulli效应）： $$\Delta P_1 = \frac{1}{2}\rho v^2\left(1-\beta^4\right) + K_c \frac{1}{2}\rho v^2$$ 其中$\beta = d_{nozzle}/d_{filament}$，$K_c \approx 0.5$是收缩系数

2. 熔化区压降（温度梯度区）： $$\Delta P_2 = \int_0^{L_m} \frac{8\eta(T(z))Q}{\pi R^4(z)}dz$$ 考虑粘度的温度依赖性和可能的径向温度分布

3. 喷嘴压降（主要阻力）： $$\Delta P_3 = \frac{8\eta Q L_n}{\pi R_n^4} \cdot \left(\frac{3n+1}{n}\right)^n$$ 对于幂律流体的修正Hagen-Poiseuille方程

4. Bowden管摩擦（远程挤出）： $$\Delta P_4 = \frac{4f L_{tube}}{d_{tube}} \cdot \frac{1}{2}\rho v^2$$ 其中$f$是Darcy摩擦因子

设计准则：$\sum \Delta P_i < 0.8 P_{drive}$，留20%安全裕度防止打滑。典型系统总压降2-5 MPa。

5.2 玻璃化转变与结晶动力学

5.2.1 玻璃化转变温度$T_g$

玻璃化转变是聚合物从玻璃态到高弹态的二级相变，涉及链段运动的解冻。$T_g$不是固定点而是一个转变区间（通常5-10°C），其值决定了材料的使用温度上限和加工温度下限。

常见3D打印材料的热性能参数：

| 材料 | $T_g$ (°C) | $T_m$ (°C) | 打印温度 (°C) | 热床温度 (°C) | HDT (°C) |

$T_g$的影响因素包括：

• 分子结构：刚性基团提高$T_g$（苯环>脂肪链）
• 分子量：Fox-Flory方程：$T_g = T_{g,\infty} - K/M_n$
• 增塑剂：每1%增塑剂降低$T_g$约3-5°C
• 交联度：轻度交联提高$T_g$约10-20°C
• 测试条件：升温速率每增加10°C/min，$T_g$升高2-3°C

对于共聚物和共混物，$T_g$预测： $$\frac{1}{T_g} = \frac{w_1}{T_{g1}} + \frac{w_2}{T_{g2}}$$ （Fox方程，相容体系）$$T_g = w_1 T_{g1} + w_2 T_{g2} + w_1 w_2 \Delta T_{12}$$ （Gordon-Taylor方程，部分相容）

5.2.2 结晶动力学：Avrami方程

半结晶聚合物的结晶行为对打印件的力学性能、尺寸稳定性和透明度有决定性影响。等温结晶动力学遵循Avrami方程：

$$X_c(t) = 1 - \exp(-Kt^n)$$ 其中：

• $X_c$：相对结晶度（0-1）
• $K$：结晶速率常数，$K = K_g \exp(-\Delta G^*/kT)$
• $n$：Avrami指数，反映成核和生长机制

Avrami指数$n$的物理意义： | $n$值 | 成核类型 | 生长维度 | 典型情况 |

结晶速率的温度依赖性呈钟形曲线： $$G(T) = G_0 \exp\left(-\frac{E_d}{R(T-T_g)}\right) \exp\left(-\frac{K_g T_m^2}{T\Delta T \Delta h_f}\right)$$ 第一项是扩散项（WLF型），第二项是成核项（Hoffman-Lauritzen理论）。最大结晶速率温度： $$T_{max} \approx 0.85 T_m \text{（均相成核）或 } 0.80 T_m \text{（异相成核）}$$

5.2.3 热历史效应

打印过程的复杂热历史（加热-冷却-再加热循环）显著影响最终微观结构。非等温结晶采用Ozawa扩展： $$X_c(T,t) = 1 - \exp\left[-\left(\frac{K(T)}{\phi}\right)^m\right]$$ 其中$\phi = dT/dt$是冷却速率。

实际打印中的热历史可分为四个阶段：

1. 初始沉积（$t < 0.1s$）：$T = T_{nozzle}$，熔体状态
2. 快速冷却（$0.1s < t < 1s$）：$\dot{T} \approx -100°C/s$，形成过冷熔体
3. 层间再加热（新层沉积时）：$\Delta T_{reheat} = 20-50°C$，退火效应
4. 缓慢冷却（$t > 60s$）：$\dot{T} \approx -0.5°C/s$，二次结晶

最终结晶度的积分模型： $$X_c^{final} = \int_0^{t_{total}} \frac{dX_c}{dt}dt = \int_{T_{nozzle}}^{T_{room}} \frac{dX_c}{dT} \cdot \frac{dT}{dt} dt$$ 典型结果：

• 薄壁快冷：$X_c \approx 5-15\%$，韧性好，半透明
• 厚壁缓冷：$X_c \approx 30-45\%$，刚性高，不透明
• 退火处理：$X_c$可增加10-20\%

5.2.4 体积收缩模型

打印过程的体积变化源于热收缩和结晶收缩： $$\varepsilon_v = \varepsilon_{thermal} + \varepsilon_{crystallization} + \varepsilon_{relaxation}$$ 各项贡献：

1. 热收缩（占60-70%）： $$\varepsilon_{thermal} = \int_{T_{print}}^{T_{room}} \alpha(T) dT$$ 其中$\alpha(T)$在$T_g$附近有突变

2. 结晶收缩（占20-30%）： $$\varepsilon_{crystallization} = \beta \cdot \Delta X_c = \beta(X_c^{final} - X_c^{initial})$$ 典型值：PLA的$\beta \approx 0.02-0.025$，PA6的$\beta \approx 0.035-0.04$

3. 应力松弛（占5-10%）： $$\varepsilon_{relaxation} = \frac{\sigma_0}{E} \left(1 - \exp(-t/\tau)\right)$$ 总收缩率预测：

• PLA：1.5-3.5%（低结晶度）
• ABS：0.7-0.9%（无定形）
• PA6：1.5-2.5%（吸湿后降低）
• PC：0.6-0.8%（低收缩）

补偿策略：

• 几何放大：$L_{CAD} = L_{target}/(1-\varepsilon_{predicted})$
• 工艺控制：降低$\Delta T$，使用加热腔
• 材料改性：添加成核剂加速结晶，减少后收缩

5.3 层间结合强度模型

5.3.1 熔合理论：爬行动力学

层间结合通过聚合物链的相互扩散实现，扩散深度$d$遵循： $$d(t,T) = \sqrt{D(T) \cdot t}$$ 其中扩散系数$D(T)$为： $$D(T) = D_0 \exp\left(-\frac{E_a}{RT}\right)$$ $E_a$是活化能（典型值50-100 kJ/mol），$R$是气体常数。

5.3.2 接触温度模型

新层沉积时的界面温度$T_i$决定了熔合程度： $$T_i = \frac{T_{new}\sqrt{\alpha_{new}} + T_{old}\sqrt{\alpha_{old}}}{\sqrt{\alpha_{new}} + \sqrt{\alpha_{old}}}$$ 其中$\alpha$是热扩散率，下标表示新/旧层。

5.3.3 熔合度定量评估

定义熔合度$\psi$为： $$\psi = \frac{1}{t_c} \int_0^{t_c} \exp\left(-\frac{E_a}{R}\left(\frac{1}{T(t)} - \frac{1}{T_{ref}}\right)\right)dt$$ 其中$t_c$是接触时间，$T_{ref}$是参考温度（通常取打印温度）。

当$\psi > 0.3$时，层间强度达到材料本体强度的80%以上。

5.3.4 强度预测模型

层间拉伸强度$\sigma_{inter}$与熔合度的关系： $$\sigma_{inter} = \sigma_{bulk} \cdot (1 - e^{-a\psi^b})$$ 其中$\sigma_{bulk}$是材料本体强度，$a \approx 3$，$b \approx 1.5$是拟合参数。

        层间结合示意图

    t=0:  新层沉积
    ━━━━━━━━━━━━ 250°C
    ────────────── 150°C (旧层表面)

    t=0.1s: 界面扩散
    ━━━━━━━━━━━━ 
    ≈≈≈≈≈≈≈≈≈≈≈≈ 200°C (混合区)
    ──────────────

    t=1s: 熔合完成
    ━━━━━━━━━━━━ 
    ============= (熔合界面)
    ──────────────

5.4 翘曲应力分析与预测

5.4.1 残余应力形成机理

打印过程中的温度梯度导致差异收缩，产生残余应力。对于约束收缩，应力为： $$\sigma_{res} = E \cdot \alpha \cdot \Delta T \cdot f_{constraint}$$ 其中$E$是弹性模量，$\alpha$是热膨胀系数，$\Delta T$是温度差，$f_{constraint}$是约束因子（0-1）。

5.4.2 翘曲变形预测：Timoshenko双金属片模型

对于长度$L$、厚度$h$的打印件，翘曲曲率$\kappa$为： $$\kappa = \frac{6\alpha\Delta T}{h} \cdot \frac{(1+m)^2}{3(1+m)^2 + (1+mn)(m^2+\frac{1}{mn})}$$ 其中$m=h_2/h_1$是层厚比，$n=E_2/E_1$是模量比。

最大翘曲高度： $$\delta_{max} = \frac{\kappa L^2}{8}$$

5.4.3 应力演化方程

考虑粘弹性松弛的应力演化： $$\frac{d\sigma}{dt} = E\left(\alpha\frac{dT}{dt} - \frac{\sigma}{\tau_r}\right)$$ 其中$\tau_r$是松弛时间，温度依赖性遵循Arrhenius关系： $$\tau_r(T) = \tau_0 \exp\left(\frac{\Delta H}{RT}\right)$$

5.4.4 翘曲控制策略

1. 温度控制：维持$\Delta T < 30°C$
2. 几何补偿：预变形量$\delta_{pre} = -\delta_{predicted}$
3. 材料改性：添加填料降低$\alpha$（碳纤维可降低80%）
4. 工艺优化：采用对称打印路径减小$f_{constraint}$

5.5 复合材料与纤维增强

5.5.1 混合定律与有效性能

短纤维增强复合材料的有效模量： $$E_{composite} = \eta_0 \eta_l V_f E_f + (1-V_f)E_m$$ 其中$V_f$是纤维体积分数，$\eta_0$是纤维取向因子（随机取向$\eta_0=0.375$，单向$\eta_0=1$），$\eta_l$是长度效率因子： $$\eta_l = 1 - \frac{\tanh(\beta l/2)}{\beta l/2}$$

$$\beta = \sqrt{\frac{2G_m}{E_f r^2 \ln(R/r)}}$$

5.5.2 纤维取向动力学

Jeffery方程描述剪切流中纤维取向演化： $$\frac{D\mathbf{p}}{Dt} = \mathbf{\Omega} \cdot \mathbf{p} + \lambda(\mathbf{D} \cdot \mathbf{p} - \mathbf{p}(\mathbf{p} \cdot \mathbf{D} \cdot \mathbf{p}))$$ 其中$\mathbf{p}$是纤维方向向量，$\mathbf{D}$和$\mathbf{\Omega}$分别是变形率和涡度张量，$\lambda = (r_e^2-1)/(r_e^2+1)$，$r_e$是纤维长径比。

5.5.3 界面剪切强度

纤维-基体界面强度通过Kelly-Tyson模型评估： $$\tau_i = \frac{\sigma_f d_f}{2l_c}$$ 其中$l_c$是临界纤维长度，$d_f$是纤维直径。对于有效增强，需要$l > 15l_c$。

5.5.4 打印参数对复合材料的影响

喷嘴磨损率（Archard方程）： $$\frac{dV_{wear}}{dt} = K_{wear} \frac{F_n \cdot v_{filament}}{H_{nozzle}}$$ 建议：

• 碳纤维含量<20%：使用硬化钢喷嘴
• 碳纤维含量>20%：使用红宝石或碳化钨喷嘴
• 玻璃纤维：降低打印速度30-50%

5.6 参数空间优化方法

5.6.1 多目标优化框架

定义目标函数向量： $$\mathbf{F}(\mathbf{x}) = [f_1(\mathbf{x}), f_2(\mathbf{x}), ..., f_n(\mathbf{x})]^T$$ 其中$\mathbf{x} = [T_{nozzle}, T_{bed}, v_{print}, h_{layer}, ...]^T$是参数向量。

典型目标包括：

• $f_1$：打印时间最小化
• $f_2$：表面粗糙度最小化
• $f_3$：机械强度最大化
• $f_4$：尺寸精度最大化

5.6.2 约束条件

物理约束： $$g_1: T_{nozzle} \in [T_g + 100, T_d - 20]$$ $$g_2: v_{print} \cdot h_{layer} \cdot w_{line} = Q_{extrude}$$ $$g_3: \frac{v_{print}}{a_{max}} < t_{segment}$$ 质量约束： $$g_4: Ra < Ra_{max}$$ $$g_5: \sigma_{inter}/\sigma_{bulk} > 0.8$$

5.6.3 Pareto前沿求解

使用NSGA-II算法寻找Pareto最优解集：

1. 初始化种群$P_0$，大小$N$
2. 对每代$t$： - 生成子代$Q_t$通过交叉变异 - 合并$R_t = P_t \cup Q_t$ - 非支配排序和拥挤度计算 - 选择$P_{t+1}$
3. 收敛判据：$\Delta_{hypervolume} < \epsilon$

5.6.4 响应面方法

构建二阶响应面模型： $$y = \beta_0 + \sum_{i=1}^k \beta_i x_i + \sum_{i=1}^k \beta_{ii} x_i^2 + \sum_{i<j} \beta_{ij} x_i x_j$$ 使用Box-Behnken或中心复合设计进行实验规划，最小实验次数： $$N_{exp} = 2^k + 2k + n_c$$ 其中$k$是因子数，$n_c$是中心点重复次数。

本章小结

本章建立了3D打印材料行为的定量模型：

1. 流变学基础：幂律流体模型预测挤出行为，WLF方程描述温度依赖性
2. 相变动力学：Avrami方程预测结晶度，热历史决定最终性能
3. 界面熔合：扩散动力学控制层间结合，熔合度$\psi>0.3$确保强度
4. 应力变形：Timoshenko模型预测翘曲，残余应力$\sigma_{res}=E\alpha\Delta T$
5. 复合材料：混合定律预测有效模量，纤维取向影响增强效果
6. 参数优化：多目标优化寻找Pareto前沿，响应面方法减少实验

关键公式汇总：

• 挤出流率：$Q \propto (\Delta P)^{1/n}$
• 层间强度：$\sigma_{inter} = \sigma_{bulk}(1-e^{-a\psi^b})$
• 翘曲高度：$\delta_{max} = \kappa L^2/8$
• 复合模量：$E_c = \eta_0\eta_l V_f E_f + (1-V_f)E_m$

练习题

基础题

习题5.1 某PLA材料在220°C时粘度为500 Pa·s，已知其WLF常数$C_1=17.44$，$C_2=51.6K$，参考温度$T_r=110°C$。计算200°C时的粘度。

Hint: 使用WLF方程计算时温叠加因子

习题5.2 喷嘴直径0.4mm，长度2mm，要达到12mm³/s的挤出流率，幂律指数n=0.3，稠度系数K=1000 Pa·s^n。计算所需压降。

Hint: 使用圆管幂律流动公式

习题5.3 打印件从250°C冷却到25°C，材料热膨胀系数$\alpha=70 \times 10^{-6}$/K，弹性模量E=2.5 GPa。计算完全约束下的残余应力。

Hint: 使用热应力公式

挑战题

习题5.4 设计一个实验确定PLA的Avrami参数。已知等温结晶实验在130°C下进行，用DSC测得不同时间的结晶度数据。如何从数据提取n和K值？

Hint: 对Avrami方程取双对数

习题5.5 碳纤维增强PLA，纤维体积分数15%，长度2mm，直径7μm，纤维模量230 GPa，基体模量3.5 GPa。考虑打印时纤维沿挤出方向部分取向（取向因子0.6），计算复合材料纵向模量。

Hint: 使用修正的混合定律，考虑长度效率因子

习题5.6 双材料打印，底层ABS（$T_g=105°C$，打印温度240°C），顶层PC（$T_g=147°C$，打印温度280°C）。评估界面兼容性并提出改进方案。

Hint: 考虑界面温度、热膨胀系数差异和化学相容性

习题5.7 设计一个响应面实验优化PETG打印参数。目标：最大化层间强度，最小化翘曲。因子：喷嘴温度(230-250°C)、床温(70-90°C)、打印速度(30-60mm/s)。

Hint: 使用Box-Behnken设计

常见陷阱与错误

材料存储错误

• 陷阱：忽视材料吸湿性，直接使用
• 后果：打印时水汽爆裂，表面起泡，强度下降50%
• 解决：Nylon需80°C烘干4-6小时，PETG需65°C烘干2-4小时

温度设置误区

• 陷阱：盲目提高温度追求流动性
• 后果：材料降解（分子量下降），发黄，释放有害气体
• 诊断：拉丝严重，层间渗透过度，细节丢失
• 规则：保持在$T_d - 20°C$以下，优先调整速度而非温度

冷却策略失误

• 陷阱：PLA过度冷却，ABS冷却不足
• 影响：PLA过冷导致层间脆弱，ABS欠冷导致变形
• 平衡：桥接需100%冷却，大面积实体需0-30%冷却

复合材料处理

• 陷阱：使用标准喷嘴打印研磨性材料
• 磨损速率：碳纤维PLA可在100m内磨损0.4mm黄铜喷嘴至0.6mm
• 监控：每50小时校准挤出倍率，偏差>5%需更换喷嘴

参数耦合忽视

• 误区：独立调整参数
• 实际：层高改变影响冷却时间，进而影响层间结合
• 方法：维持层时间常数=$h_{layer}/v_{print}$

最佳实践检查清单

新材料导入流程

• [ ] DSC测试确定$T_g$、$T_m$、$T_d$
• [ ] 流变测试获得粘度-温度曲线
• [ ] 打印温度塔（5°C梯度）
• [ ] 回抽测试确定最优回抽距离
• [ ] 桥接测试验证冷却需求
• [ ] 拉伸试样测试层间强度
• [ ] 记录最优参数窗口

打印前材料准备

• [ ] 检查材料干燥度（Nylon<0.2%水分）
• [ ] 确认喷嘴材质匹配
• [ ] 预热热床至$T_g - 10°C$
• [ ] 首层温度提高5-10°C
• [ ] 准备适当脱模剂

质量监控点

• [ ] 首层线宽=喷嘴直径×1.2
• [ ] 层间无明显分离线
• [ ] 拐角无过度渗出
• [ ] 支撑易剥离（间隙0.2mm）
• [ ] 尺寸偏差<±0.2mm
• [ ] 表面粗糙度Ra<层高×0.3

故障诊断顺序

1. [ ] 检查挤出倍率（±5%内）
2. [ ] 验证温度实际值（热电偶）
3. [ ] 测量丝材直径（±0.05mm）
4. [ ] 检查喷嘴磨损
5. [ ] 评估环境温度/湿度
6. [ ] 审查切片参数一致性

文档记录标准

• [ ] 材料批次号和供应商
• [ ] 环境温度湿度
• [ ] 所有打印参数
• [ ] 失败模式和解决方案
• [ ] 成功率统计
• [ ] 成本时间记录