本章深入探讨3D打印材料的物理化学特性，建立材料行为与打印参数之间的定量关系。通过流变学、热力学和力学模型，您将学会预测和优化打印过程中的材料响应，从而实现高质量、高可靠性的打印结果。
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⁻¹，远高于传统成型工艺。
聚合物粘度强烈依赖于温度，可用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%，这解释了为什么温度控制对打印质量至关重要。
聚合物离开喷嘴后会发生弹性恢复（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之间，具体取决于：
实际打印中，通过调整挤出倍率（extrusion multiplier）补偿胀大效应，典型值0.9-1.1。
在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$是丝材截面积。
系统总压降由多个组分构成：
热端入口收缩损失（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$是收缩系数
熔化区压降（温度梯度区）： \(\Delta P_2 = \int_0^{L_m} \frac{8\eta(T(z))Q}{\pi R^4(z)}dz\) 考虑粘度的温度依赖性和可能的径向温度分布
喷嘴压降（主要阻力）： \(\Delta P_3 = \frac{8\eta Q L_n}{\pi R_n^4} \cdot \left(\frac{3n+1}{n}\right)^n\) 对于幂律流体的修正Hagen-Poiseuille方程
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。
玻璃化转变是聚合物从玻璃态到高弹态的二级相变，涉及链段运动的解冻。$T_g$不是固定点而是一个转变区间（通常5-10°C），其值决定了材料的使用温度上限和加工温度下限。
常见3D打印材料的热性能参数：
| 材料 | $T_g$ (°C) | $T_m$ (°C) | 打印温度 (°C) | 热床温度 (°C) | HDT (°C) |
|---|---|---|---|---|---|
| PLA | 55-65 | 150-180 | 190-220 | 20-60 | 50-55 |
| ABS | 105 | 无定形 | 220-250 | 80-110 | 96-98 |
| PETG | 80 | 245 | 220-250 | 60-80 | 70-73 |
| PC | 147 | 无定形 | 250-310 | 100-120 | 138-143 |
| PA6 | 47 | 220 | 240-270 | 60-80 | 55-80 |
| PA12 | 70 | 178 | 240-270 | 60-80 | 55-80 |
| TPU | -35到-5 | 150-220 | 210-240 | 20-60 | 60-80 |
$T_g$的影响因素包括：
对于共聚物和共混物，$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方程，部分相容）
半结晶聚合物的结晶行为对打印件的力学性能、尺寸稳定性和透明度有决定性影响。等温结晶动力学遵循Avrami方程：
\[X_c(t) = 1 - \exp(-Kt^n)\]其中：
Avrami指数$n$的物理意义： | $n$值 | 成核类型 | 生长维度 | 典型情况 | |——-|———|———|———| | 4 | 均相成核 | 3D球晶 | 过冷熔体 | | 3 | 异相成核 | 3D球晶 | 有核剂时 | | 2+1 | 瞬时成核 | 2D片晶+1D加厚 | 薄膜 | | 2 | 瞬时成核 | 2D圆盘 | 受限空间 | | 1 | 瞬时成核 | 1D纤维 | 拉伸结晶 |
结晶速率的温度依赖性呈钟形曲线：
\[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{（异相成核）}\]打印过程的复杂热历史（加热-冷却-再加热循环）显著影响最终微观结构。非等温结晶采用Ozawa扩展：
\[X_c(T,t) = 1 - \exp\left[-\left(\frac{K(T)}{\phi}\right)^m\right]\]其中$\phi = dT/dt$是冷却速率。
实际打印中的热历史可分为四个阶段：
最终结晶度的积分模型：
\[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\]典型结果：
打印过程的体积变化源于热收缩和结晶收缩：
\[\varepsilon_v = \varepsilon_{thermal} + \varepsilon_{crystallization} + \varepsilon_{relaxation}\]各项贡献：
热收缩（占60-70%）： \(\varepsilon_{thermal} = \int_{T_{print}}^{T_{room}} \alpha(T) dT\) 其中$\alpha(T)$在$T_g$附近有突变
结晶收缩（占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$
应力松弛（占5-10%）： \(\varepsilon_{relaxation} = \frac{\sigma_0}{E} \left(1 - \exp(-t/\tau)\right)\)
总收缩率预测：
补偿策略：
层间结合通过聚合物链的相互扩散实现，扩散深度$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$是气体常数。
新层沉积时的界面温度$T_i$决定了熔合程度：
\[T_i = \frac{T_{new}\sqrt{\alpha_{new}} + T_{old}\sqrt{\alpha_{old}}}{\sqrt{\alpha_{new}} + \sqrt{\alpha_{old}}}\]其中$\alpha$是热扩散率，下标表示新/旧层。
定义熔合度$\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%以上。
层间拉伸强度$\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: 熔合完成
━━━━━━━━━━━━
============= (熔合界面)
──────────────
打印过程中的温度梯度导致差异收缩，产生残余应力。对于约束收缩，应力为：
\[\sigma_{res} = E \cdot \alpha \cdot \Delta T \cdot f_{constraint}\]其中$E$是弹性模量，$\alpha$是热膨胀系数，$\Delta T$是温度差，$f_{constraint}$是约束因子（0-1）。
对于长度$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}\]考虑粘弹性松弛的应力演化：
\[\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)\]短纤维增强复合材料的有效模量：
\[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)}}\]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$是纤维长径比。
纤维-基体界面强度通过Kelly-Tyson模型评估：
\[\tau_i = \frac{\sigma_f d_f}{2l_c}\]其中$l_c$是临界纤维长度，$d_f$是纤维直径。对于有效增强，需要$l > 15l_c$。
喷嘴磨损率（Archard方程）：
\[\frac{dV_{wear}}{dt} = K_{wear} \frac{F_n \cdot v_{filament}}{H_{nozzle}}\]建议：
定义目标函数向量：
\[\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$是参数向量。
典型目标包括：
物理约束： \(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\)
使用NSGA-II算法寻找Pareto最优解集：
构建二阶响应面模型：
\[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打印材料行为的定量模型：
关键公式汇总：
习题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设计