玻尔百科复合材料通过组合两种或多种不同物质而制成，其性能优于单一组分，从而在从航空航天到生物力学等领域实现了突破性进展。对于科学家和工程师而言，核心挑战在于如何在不进行昂贵且耗时的物理测试的情况下预测这些复杂材料的行为。本文旨在填补这一知识空白，带领读者全面进入复合材料建模的世界。本文旨在让读者对如何通过数学描述和计算模拟这些材料，从而设计出更好、更轻、更强的结构有一个基本的理解。
接下来的章节将引导您探索这个迷人的主题。首先，在“原理与机理”中，我们将深入探讨用于预测复合材料刚度、强度和失效的基础理论，从简单的平均法则开始，逐步发展到复杂的计算技术。然后，在“应用与跨学科联系”中，我们将探索这些模型如何应用于解决现实世界中的工程问题，并发现“复合材料思维方式”如何为结构力学之外的多个领域提供深刻的见解。
让我们从一个简单的问题开始我们的旅程：到底什么是复合材料？您可能会想到高科技的例子，如一级方程式赛车或现代飞机中的碳纤维。但这个想法本身与文明一样古老。想想用稻草加固的泥砖，或用钢筋增强的现代混凝土。原理是相同的：取两种或多种不同的材料，将它们组合起来，创造出一种性能优于或仅仅是不同于其单一组分的新材料。它不仅仅是盐和胡椒那样的混合物；它是一个团队，其中成员——坚固、刚硬的纤维和将其固定在一起的周围基体——协同工作。
作为科学家和工程师，我们面临的巨大挑战是预测这个团队的行为，而不必构建和测试每一种可能的组合。我们希望建立一套理论，一套规则，能够根据其组成部分的属性告诉我们整体的属性——即有效或均质化属性。
想象一个用于笔记本电脑热管的简单热模型，其设计目的是将热量从处理器中带走。一个初步的、幼稚的模型可能会假设热管由单一、均匀的材料制成。但如果为了节省成本，它实际上是由两段不同金属连接而成的呢？这种仅仅进行平均的简化模型将会是错误的。中点的实际温度将取决于两段材料的热导率以及它们的排列方式。热量必须串联地“流”过它们，而总的热阻是各部分热阻之和。这个简单的例子给了我们一个深刻的教训：在复合材料中，几何构型——即组分的精确排列方式——与组分本身同等重要。
最早且最简单的捕捉这一点的尝试是著名的混合律。让我们想象一种由纤维和基体组成的简单复合材料。它的刚度（即杨氏模量 ）是多少？存在两种极端的、理想化的情景。
首先，想象沿着纤维方向拉伸复合材料。如果我们假设纤维和基体完美粘合且拉伸量相同——即均匀应变条件——那么总刚度就是组分刚度的加权平均。这就是 Voigt 模型，它给出了复合材料刚度的上限。
这里， 和 分别是纤维和基体的体积分数——它们占总体积的百分比。这是“乐观主义者”的估计。
现在，想象垂直于纤维方向拉伸复合材料。如果我们这次假设纤维和基体承受相同的应力——即均匀应力条件——我们会得到一个不同的答案。总柔度（刚度的倒数）是组分柔度的加权平均。这就是 Reuss 模型，它给出了一个下限。
这是“悲观主义者”的估计。复合材料的真实刚度位于这两个界限之间，即所谓的 Voigt 和 Reuss 界限。要使用这些公式，我们需要知道体积分数。在实践中，测量每种组分的重量要容易得多。这意味着我们需要一种方法，将容易测量的重量分数转换为物理上至关重要的体积分数，这个计算依赖于材料的密度，甚至必须考虑像空隙这样的微小制造缺陷。
Voigt 和 Reuss 模型优美而简单，但现实往往更为复杂。虽然 Voigt 模型在预测沿纤维方向的刚度方面表现出色，但对于横跨纤维方向的性能，这两个模型都可能极不准确。为什么呢？因为均匀应变或均匀应力的假设几乎从未成立。
当您垂直于纤维方向加载复合材料时，材料内部的应力场会变成一个旋转的、复杂的图案。应力必须“绕流”过非常刚硬的纤维，从而在基体中产生高应力和低应力区域。忽略了这种错综复杂变化的简单模型无法捕捉到其物理本质。
为了做得更好，我们需要更复杂的模型。一个杰出的例子是 Halpin-Tsai 关系式。Halpin 和 Tsai 没有从一个笼统的假设出发，而是采取了一种更务实、更具物理学家风格的方法。他们的模型是一个巧妙的插值，一个旨在提高准确性的半经验公式。它看起来是这样的：
这里， 是我们想要预测的属性（如横向模量 ）， 和 是基体和纤维的属性，而 是一个取决于纤维与基体属性比值的项。其中的奥妙在于参数 。这是一个可调整的“形状参数”，用于解释增强体的几何形状——无论是圆形纤维、方形纤维，还是短粗的颗粒。
这种形式的天才之处在于，它在数学结构上保证了在已知的物理极限情况下的正确性。对于极少数纤维（）的情况，它正确地再现了无限大基体中单根纤维的精确解。而随着纤维变得越来越密集，分母项 巧妙地以一种唯象的方式捕捉了“多夹杂物相互作用”的效应。增加更多纤维对刚度的提升效果会递减，这是简单混合律所忽略的一种非线性。这是一个绝佳的例子，展示了如何运用一点数学上的巧思——在此例中是使用一种称为 Padé 近似的有理函数——从一个已知的简单案例通向一个复杂的、普遍的现实。
问题的几何形状和约束条件的重要性怎么强调都不为过。思考一下一块非常薄的复合材料片——单层板——和一块非常厚的块体之间的区别。对于薄片，我们通常可以假设它处于平面应力状态，这意味着其厚度方向上没有应力，因为它可以自由膨胀或收缩。然而，对于厚块体，中间的材料受到周围材料的约束，阻止其在厚度方向上变形。这是一种平面应 strain 状态。对于像复合材料这样的各向异性材料，你计算出的有效平面内刚度会根据你所做的假设而根本不同。操作的顺序——应用物理假设与应用几何极限——是不可交换的。这是一个微妙但深刻的提醒：我们的模型的好坏取决于它们所建立的假设。
像 Halpin-Tsai 这样的解析模型功能强大，但它们有其局限性。对于具有真正复杂、随机微观结构的复合材料该怎么办？为此，我们求助于计算机，创建一个“数字实验室”。这个想法被称为多尺度建模。我们不可能逐个原子地模拟整个飞机机翼，所以我们利用了尺度的分离。
我们从材料中分离出一小块重复的、在统计上能代表整体的部分。这被称为代表性体积单元 (RVE)。然后，我们仅对这个微小的 RVE 进行详细的有限元 (FE) 模拟，以计算出其有效性能。但是，你如何对一个从更大材料中切割出来的小盒子施加载荷呢？你不可能直接抓住它！
答案是一个优美而巧妙的数学技巧：周期性边界条件。我们告诉计算机，我们 RVE 盒子的右面“粘合”到它的左面，顶面粘合到底面。一个面上的任何位移都会在对面精确地镜像。通过这样做，我们模拟了一种在所有方向上完美无限重复的材料。我们正在建模的是真正的“体”材料，不受我们微小计算机模型任何人为边界效应的影响。
这种计算方法不仅仅是一个方便的技巧；它建立在物理学中一些最深刻、最美丽的数学之上。次可加遍历定理为均质化方法为何有效提供了严格的基础，即使对于具有随机、非周期性微观结构的材料也是如此。这个定理告诉我们，如果一个随机材料是“遍历的”——意味着一个足够大的样本在统计上与任何其他大样本无法区分——那么，当我们取越来越大的 RVE 时，复杂、随机的微观行为将平均化为一个单一的、可预测的、确定性的宏观行为。本质上，随机性在宏观尺度上被“冲刷”掉了。正是这种魔力让我们能够用一个简单、干净、有效的模型取代那个混乱、复杂的微观世界，从而用于设计现实世界的结构。
到目前IE，我们讨论了刚度和变形。但是当复合材料断裂时会发生什么呢？这是失效准则的领域。与只有一个屈服强度的简单金属不同，复合材料的强度要复杂得多。它的强度取决于载荷的方向，以及不同类型应力——拉伸、压缩和剪切——同时作用的组合。
为了捕捉这一点，我们使用交互式失效准则，其中最著名的是Tsai-Wu 准则。它在应力空间中定义了一个“失效面”。想象一个多维空间，其中每个轴代表一个不同的应力分量（, , 等）。Tsai-Wu 准则描述了这个空间中的一个椭球体。如果材料中任何一点的应力状态位于椭球体内部，则材料是安全的。如果它接触或超出了椭球体的表面，失效就开始了。
其通用方程看起来令人望而生畏，是一个包含许多项的多项式。但在这里，物理学家的朋友——对称性——再次伸出援手。对于一个正交各向异性材料（一种具有三个相互垂直对称平面的材料，如单向单层板），我们可以证明，如果我们翻转剪切应力的符号，材料的响应必须相同。这一物理要求迫使通用多项式中的许多系数为零，从而将准则显著简化为一个可处理的形式。
现在，我们必须做一个关键的区分。当金属在延性意义上“失效”时，它会屈服。它会发生永久的塑性变形，就像弯曲一个回形针一样。它仍然可以承受载荷。而当复合材料铺层失效时，通常是一个脆性事件。基体中形成微裂纹，或纤维断裂。材料受损了；它从根本上失去了刚度以及以同样方式承载载荷的能力。
我们可以在计算机上模拟这个被称为渐进性失效的过程。想象一个由许多不同角度的铺层组成的层合板。我们施加一个小的载荷，然后在我们的模拟中，我们逐个检查每个铺层并询问：“根据你所承受的应力，你是否已根据 Tsai-Wu 准则失效？”如果所有铺层的答案都是“否”，我们就增加载荷。如果有一个铺层回答“是！”，我们不会移除它，而是“惩罚”它。在铺层折扣法中，我们大幅降低其刚度矩阵，使其“变软”且效果降低。现在，整个层合板变弱了。原来由失效铺层承担的载荷必须重新分配给它的邻居。我们重新运行分析。这种重新分配可能会导致相邻的铺层失效。我们重复这个过程——加载、检查、失效、重新分配——观察损伤在层合板中级联传播，直到整个结构再也无法承受载荷。这就是我们如何通过计算来预测复合结构在接近最终失效时复杂的、逐层“解体”的过程。
当然，现实世界还有更多的花样。在层合板的自由边缘——例如，一个切口或面板的侧面——会产生二维理论无法预测的复杂三维应力。这些层间应力可能导致各层剥离，即分层。预测这些情况需要详细的三维有限元模型和一个非常仔细的后处理工作流程，以确保结果在物理上有意义且在数值上是准确的。
这把我们带到了复合材料建模前沿最深刻的挑战之一。当我们试图在我们先进的多尺度（）模拟中实施一个简单的损伤模型——即刚度随应变减小，称为软化——时，会发生一场数学灾难。
数学告诉损伤要集中在一个无限薄的带中——一个体积为零的裂纹。这意味着破坏材料所需的能量为零，这在物理上是荒谬的。在计算机模拟中，这个“裂纹”会盲目地沿着有限元网格的线条发展。改变网格，你就会改变答案。这个问题被称为是不适定的；解决方案具有病态的网格依赖性。
我们如何解决这个问题？我们意识到我们的简单模型缺少了一部分物理原理。断裂并非纯粹的局部事件。一个点的材料状态必须取决于其邻近区域发生的情况。我们必须在我们的模型中引入一个内禀长度尺度来对其进行正则化。
有几种巧妙的方法可以做到这一点。我们可以使用非局部模型，其中一个点的损伤是其周围小体积内应变的平均值。我们可以使用梯度增强模型，其中材料的能量不仅取决于应变，还取决于损伤的空间梯度——产生一个尖锐的损伤前沿是需要能量的。或者我们可以使用更先进的理论，如微形态连续介质，它通过代表微观结构的自由度来丰富宏观模型本身。
所有这些方法都有相同的效果：它们将损伤涂抹在一个微小但有限的宽度上，确保产生裂纹的能量非零，并且随着网格细化，模拟结果收敛到一个单一的、物理上有意义的答案。通过驯服这种不稳定性，我们正在更接近物质的真实本质，它不是一个抽象的、局部的连续体，而是一个复杂的、相互作用的系统，其中一个点的行为总是与其周围环境相联系。这就是复合材料建模的美丽而持续的探索：构建日益逼真的数学图像，以描绘材料错综复杂而又优雅的现实。
在我们之前的讨论中，我们揭示了复合材料的基本原理。我们了解到，通过巧妙地组合不同的物质——例如，将强纤维置于柔顺的基体中——我们可以创造出性能超越其单一组分的新材料。在某种意义上，我们已经成为材料建筑师。但建筑不仅仅是理解砖块和砂浆的特性；它是用它们来建造大教堂、桥梁和房屋。那么，这种新获得的建筑能力将我们引向何方？我们能用它来建造什么？
本章是一次进入广阔且往往令人惊讶的世界的旅程，在那里，复合材料科学得以应用。我们将看到这些原理如何不仅仅是抽象的规则，而是解决我们一些最紧迫工程挑战的实用工具。然后我们将发现，“复合材料思维方式”是一种通用语言，让我们能够理解远超力学范围的现象，从热量和电流的流动到生命本身的结构。
复合材料最直接、或许也是最具标志性的应用，是在追求既极其坚固又惊人轻巧的结构的宏伟探索中。想象一下现代飞机的机翼、一级方程式赛车的底盘，或高性能的自行车车架。在每一种情况下，目标都是最大化性能，而一个关键的衡量标准是刚度重量比。例如，我们如何设计一根梁，使其在尽可能抵抗弯曲的同时，不增加不必要的质量？
对于复合材料来说，这成了一个有趣的谜题，有两个可调的旋钮。首先是零件的宏观形状——它的宽度、高度和横截面几何形状。其次是材料本身的微观“配方”——碳纤维、玻璃纤维和聚合物树脂的体积分数。复合材料设计的真正威力在于，当我们意识到我们可以同时优化这两者时。我们可以问计算机：在给定的预算和可用成分下，创造出最刚、最轻的梁的形状和材料混合物的绝对最佳组合是什么？
解决方案通常涉及一个漂亮的关注点分离。最佳形状由力学定律决定，通常倾向于更高、更薄的横截面以最大化截面二次矩，就像I型梁一样。同时，最佳材料成分通过解决一个不同的问题来找到：在满足成本约束的IE提下，找到能提供最高比刚度（杨氏模量除以密度）的材料混合物。这种双重优化使工程师能够突破可能的界限，创造出用钢或铝等单一金属无法想象的结构。
但第一天的强度是不够的。一个结构必须承受一生的振动、阵风和波动的载荷。这就引出了疲劳的关键问题。任何材料，当经受数百万次的应力循环时，最终都会变弱并失效，即使峰值应力远低于它在单次拉伸中能承受的应力。对于复合材料而言，预测这种疲劳寿命至关重要。例如，直升机旋翼桨叶在每次飞行中都会经受复杂的循环载荷。
工程师们已经开发出复杂的模型来预测复合材料部件的寿命。他们发现，不仅应力循环的幅度很重要，其平均值也很重要。在一个高张力和低张力之间波动的应力循环比一个围绕零点以相同幅度波动的循环更具破坏性。通过创建一个同时考虑幅度和平均值的“损伤驱动应力度量”，然后将其代入一个幂律关系，我们可以估算出失效前的循环次数。这使我们能够设计出可以在任何危险出现之前安全退役的部件，这是现代安全工程的基石。
即使有最好的设计，我们也必须面对失效的可能性。在复合材料中，失效很少是简单、干净的断裂。它是一个复杂的过程，通常由一个微小的、几乎看不见的裂纹引发。断裂力学这门科学致力于理解这些裂纹如何扩展。在一个本身就是由不同材料粘合而成的复合材料中，当裂纹到达一个界面时，它面临一个关键的选择：是穿透到下一种材料中，还是选择阻力最小的路径转向，沿着界面传播？后一个过程，即分层，是层状复合材料的主要失效模式。
预测这种行为是一场高风险的游戏。计算科学家们开发了一整套工具来分析裂纹尖端的能量平衡，使用诸如积分和虚拟裂纹闭合技术（VCCT）等方法。选择正确的工具取决于具体情况——材料是弹性行为，还是有塑性变形？我们处理的是简单的裂纹还是两种不同材料界面处的复杂裂纹？一个合理的决策流程对于可靠的分析至关重要。裂纹本身的选择——是穿透还是偏转——可以通过比较驱动裂纹前进的可用能量与界面的韧性来预测。如果能量释放率超过了界面的断裂韧性，就可能发生分层。这个微观决策具有宏观后果，决定了整个部件的失效路径和最终强度。
到目前为止，我们讨论了分析和优化由固定菜单成分制成的复合材料。但如果我们能够设计成分本身呢？这就是多尺度建模和拓扑优化的革命性前景。在这里，我们不把材料视为给定的，而是视为一个设计空间。一个计算机算法的任务是在给定的设计域内分配有限数量的材料，以实现最大刚度等目标。
其结果往往令人叹为观止，产生出复杂的、类似骨骼的晶格和看起来有机的桁架，这些结构完美地适应了它们的力学环境。要使这种方法奏效，其基本假设是“尺度分离”原则：晶格的微观特征必须远小于大结构中应力变化的尺度。这使我们能够“均质化”微晶格的属性，并在宏观尺度上将其视为一种光滑、连续的材料。
然而，一个好的科学家，就像一个好的工程师一样，必须知道他们工具的局限性。这种尺度分离可能会失效。在尖角或集中点载荷附近，宏观应力场的变化可能非常迅速，以至于可以与微观结构本身的尺寸相媲美。在这些“边界层”中，均质化模型不再有效，需要对微观结构进行更详细、直接的模拟。理解我们的模型在哪里适用、在哪里失效，与模型本身同等重要。
此外，在计算机上运行这些复杂模拟本身就是一门艺术。有时，直接应用数值方法可能会导致奇异的、非物理的结果。复合材料有限元分析中的一个经典例子是“锁定”。当使用简单的单元类型模拟具有非常刚硬纤维的材料时，数值模型可能会变得人为地僵硬，即“锁定”，无法真实地变形。为了克服这个问题，计算科学家们设计了一些巧妙的技巧，例如“选择性减缩积分”。这项技术有意地对材料能量中过硬的部分使用一个不太精确的积分规则，同时对其余部分保持完全的精度。这是一个“少即是多”方法的绝佳例子——计算中某一部分的刻意不精确，却为整体解锁了一个更准确、更具物理意义的结果。
也许从研究复合材料中得到的最深刻的启示是，混合、分层和有效性能这些核心思想并不仅限于结构力学。它们构成了一种通用语言，描述了广泛的物理现象。
考虑热的流动。一个为隔热而设计的、由多层不同材料制成的复合墙，其行为与受载的结构复合材料非常相似。每一层都对热传导呈现出一定的阻力。要找到墙体的总热阻，我们只需将各层的阻力相加，就像我们对电路中的电阻器所做的那样。由此，我们可以推导出一个“复合毕渥数”，这是一个单一的无量纲数，它告诉我们整个墙体是会均匀冷却，还是内部会产生显著的温度梯度。这个单一的数字决定了我们是否可以使用简单的“集总电容”模型，这是热工学中的一个强大简化方法。
这种与电学的联系并不仅仅是一种类比。如果我们用两种不同电介质的交替薄片制造一种层状材料，我们实际上是创造了一叠串联的电容器。当一个电场垂直于这些层施加时，复合材料的有效介电常数可以用与串联电容器完全相同的公式来计算。这一原理是设计用于电子元件的新材料的核心，甚至在嵌段共聚物的自组装纳米结构中也能观察到 [@problem-id:374553]。
“复合”思维方式甚至更深入地延伸到能源和电子材料的世界。考虑一种混合离子电子导体（MIEC），这是现代电池和燃料电池中的关键材料。这些通常是通过混合电子导电相和离子导电相而制成的复合材料。一个关键问题是：导电相的体积分数达到多少时，整个材料才开始导电？这是一个经典的*逾渗理论问题。在临界体积分数（即“逾渗阈值”）以下，导电颗粒是绝缘海洋中的孤岛。在阈值处，第一个连续的、贯穿样品的路径形成，电导率突然飙升。即使在阈值之上，有效电导率也取决于这些路径的几何形状。载流子必须通过一条扭曲、复杂的路径，而这种几何阻抗由一个称为曲折度*的属性来表征。通过对逾渗、曲折度和各相固有属性之间相互作用的建模，我们可以预测和设计我们未来能源材料的性能 [@problem-id:2500623]。
当我们发展这些复杂的工具和概念时，令人谦卑地意识到，在许多方面，我们只是在追赶自然的脚步。生命世界是复合材料设计的终极展示，经过数十亿年的进化而臻于完美。
看看你自己的骨骼就知道了。骨骼是复合工程的杰作：一个由坚韧的胶原蛋白原纤维组成的基体，用坚硬的羟基磷灰石晶体加固。它是一种活的复合材料，不断适应其环境。正如 Wolff 定律所述，骨骼会根据其所受的机械载荷进行自我重塑。在高应力区域，骨骼变得更致密，其内部小梁结构与主应力方向对齐，从而创造出一种自然优化的结构。这种适应通过两个不同的过程发生。塑建（Modeling）改变骨骼的整体尺寸和形状，例如通过在一侧表面增加新骨而在另一侧移除骨骼来矫直弯曲的骨骼。而重塑（Remodeling）则是一个内部维护过程，其中小块的陈旧或受损骨骼被新骨替换，保持整体形态，但确保组织的健康和完整性 [@problem-id:2619272]。
这个原则不仅限于动物。植物界充满了结构复合材料的例子。一根普通草叶的茎之所以能在风中弯曲而不折断，要归功于厚壁组织——长而刚硬的木质化纤维嵌入在较软的薄壁组织基体中。为了理解这种组织的力学，植物学家和生物力学家进行的实验，材料工程师会觉得非常熟悉。他们设计弯曲测试，以仔细地解耦纤维体积分数与纤维方向的影响，旨在建立一个该组织刚度的预测模型。
从喷气式飞机的机翼到蜻蜓的翅膀，从人造压力容器到树木的木质树干，原理都是相同的。通过理解混合的艺术，我们不仅开启了工程学的新前沿，而且对我们周围的世界获得了更深刻、更统一的视角。发现之旅远未结束；建筑的可能性是真正无穷无尽的。