弹性力学中的微观结构——使用结构张量

《Mechanics Research Communications》：Microstructure in elasticity - using fabric tensors

　　

摘要

本文建立了以构织张量形式表述的连续介质弹性理论与基于微观结构影响域导出的弹性理论之间的直接理论联系。文章考察了两种本构表述：Voyiadjis–Kattan表述，其中四阶弹性张量仅用标量构织和二阶偏量构织张量表示；以及一种微观结构驱动表述，其中弹性源自影响域内材料纤维的加权取向矩。通过将微观结构表述展开为显式的二阶和四阶取向矩，并引入最小的二次闭包，弹性张量被简化为仅依赖于单个有效二阶张量。这一简化使得在横向各向同性情形下，两种理论之间可以严格地逐分量等价。推导出了闭式关系，将微观结构取向矩直接映射到控制弹性响应的标量和偏量构织张量。分析进一步得出了相容性条件，量化了两种弹性描述之间何时可以精确识别，并阐明了高阶微观结构统计在控制各向异性刚度中的作用。所得框架为微观结构力学与构织张量连续介质建模之间提供了统一桥梁，为各向异性弹性、微观结构以及非均质材料中的结构-性能关系提供了新的见解，对生物组织尤为相关。

引言

理解微观结构如何决定宏观弹性行为仍然是固体力学、材料科学和生物力学中的核心挑战。经典连续介质弹性理论通过四阶弹性张量来表示刚度，其对称性和不变性性质已得到充分确立；然而，该张量与底层材料结构（如纤维取向、孔隙分布或其他形式的内部微观结构组织）之间的联系，往往需要超越传统各向同性或唯象各向异性模型的额外本构结构。因此，建立微观结构统计与连续介质弹性之间的严格桥梁，对于非均质材料和生物材料的预测性建模至关重要。

已有两个互补的理论方向来应对这一挑战。第一个是构织张量方法，其中各向异性由从取向分布或微裂纹几何导出的标量和张量型内部变量来表示。损伤力学中的奠基性工作表明，二阶和四阶张量可以在热力学一致的框架内编码各向异性刚度和方向性材料响应[1,2]。后续工作表明，弹性本身可直接用构织张量表示，从而允许从随变形或微观结构重组而演化的低阶结构描述符构建四阶刚度张量[3]。该表述提供了一种紧凑的连续介质描述，同时保留了与微观结构取向的有意义联系。

第二个方向从显式的微观结构平均出发推导弹性。在这种表述中，刚度源自结构元素（如纤维或微组成相）在影响域内分布的加权取向矩。基于取向张量和生命加权函数的平均过程产生了依赖于二阶和四阶取向统计的弹性张量，从而将微观结构信息直接嵌入到连续介质响应中[4]。取向矩方法长期以来被用于复合材料、颗粒介质和生物组织中，以连接纤维分布与各向异性刚度[[5], [6], [7]]。这些方法提供了强有力的物理解释，但往往导致高阶张量依赖关系，使本构建模复杂化。

尽管构织张量方法和微观结构平均方法的目标都是将结构与刚度联系起来，但它们基本上是平行发展的。构织表述给出一个仅依赖于单个有效二阶张量的四阶弹性张量，而基于微观结构的弹性通常需要二阶和四阶取向矩。调和这些描述因此是一个基本的理论问题：在什么条件下，微观结构导出的弹性张量可以被构织张量弹性形式精确表示？回答这个问题将提供一个统一框架，连接损伤力学、复合材料理论和生物组织力学，并阐明连续介质构织变量的物理意义。

本文通过直接比较两种代表性表述来解决这一问题。第一种是在Voyiadjis和Kattan损伤力学框架内[3]开发的用标量和偏量构织张量表示的弹性表述。第二种是通过在影响域上对取向相关量进行加权平均得到的基于微观结构的弹性张量，如关于微观结构的近期研究[4]中所表述。通过将微观结构表述展开为显式的二阶和四阶取向矩，并引入将四阶矩与二阶矩联系起来的最小二次闭包，弹性张量被简化为仅依赖于单个有效二阶张量。

为获得解析见解，比较在横向各向同性下进行，这是一种在纤维增强材料、地质材料和软组织中广泛遇到的对称类型[[5], [6], [7], [8]]。推导出的闭式表达式将微观结构取向矩直接映射到控制弹性的标量和偏量构织张量。分析进一步揭示了两种表述精确等价所需的相容性约束，从而确定构织张量连续介质模型何时完全捕获微观结构引起的各向异性，以及何时高阶统计是不可避免的。

除理论统一之外，这一联系具有重要含义。构织张量弹性提供了通过连续介质力学中演化的结构描述符表示各向异性刚度的自然途径[1,2,9]，而基于微观结构的平均提供了对复合材料、多孔固体和生物组织[4,5,6]具有物理依据的描述。建立这些框架之间的等价性能够实现跨学科的结果转移，并支持连接微观力学与连续介质响应的多尺度本构建模策略[10]。

因此，本文的目标是在弹性理论中推导构织张量与微观结构取向矩之间的直接解析关系。所得表述提供了各向异性刚度的统一结构解释，并阐明了高阶微观结构信息在连续介质本构定律中的作用。

大量研究致力于建立微观结构组织与宏观弹性响应之间的严格联系，特别是在颗粒材料、复合材料和非均质固体中。在此背景下，Misra及其同事开发的颗粒微观力学框架[[11], [12], [13]]为本研究提供了最直接的先例之一。他们在颗粒集合高阶弹性常数的研究中，连续介质刚度张量明确地从颗粒对相互作用刚度以及描述接触取向和几何的构织参数导出，从而证明宏观弹性可以用统计定义的结构描述符表示。后续扩展到微结构和伪颗粒连续介质模型的进一步研究表明，当保留微观结构运动学信息时，额外的内部变量和长度尺度效应自然出现，突出了将微观结构简化为低阶张量的优势和局限。这些发展有力支持了本文的概念前提——即只要高阶结构统计函数依赖于二阶取向度量，连续介质弹性就可能允许构织张量表示。

互补见解来自颗粒系统中微观结构演化的研究，特别是与Rosato及其合作者[14]相关的研究。他们对被敲击或振动颗粒堆积的模拟和理论分析强调了取向分布、接触网络和密度弛豫如何在外部激励下演化。尽管这些工作没有直接推导弹性张量，但它们通过展示可测量的微观结构统计量——即潜在的构织描述符——如何从颗粒尺度动力学中产生，提供了必要的物理基础。这些结果对于本构建模很重要，因为它们阐明了取向矩何时足以表征宏观行为，以及何时必须保留更复杂的结构信息。

综合来看，先前的微观力学弹性理论和颗粒微观结构演化研究指向一种统一解释，即构织张量作为统计平均微观结构信息的连续介质载体。本文通过推导构织张量弹性表述与微观结构平均弹性张量之间的显式解析桥梁来推进这一视角，并确定两种描述在何种相容性条件下变为等价。这种综合将该理论定位为更广泛的努力的一部分，旨在跨非均质和各向异性材料系统连接颗粒尺度结构、内部变量和连续介质力学响应。

本文考察的基于微观结构的弹性表述基于Lu、Yang和Ren的近期工作[4]，其取向矩框架作为本分析的起点。

段落节选

连续介质弹性与张量表述

设连续介质体的小应变线弹性响应由本构关系描述：

σ?? = C???? ε??

其中 σ?? 和 ε?? 分别表示柯西应力张量和无穷小应变张量，C???? 是具有经典次对称性和主对称性的四阶弹性张量：

C???? = C???? = C???? = C????

对于各向同性弹性，C???? 用Lamé参数 λ 和 μ 表示如下：

C????^iso = λ δ?? δ?? + μ(δ?? δ?? + δ?? δ??)

其中弹性张量在……

构织张量与微观结构之间的解析对应

第2节的数学推导建立了微观结构取向矩 (m?, m?) 与控制连续介质弹性张量的构织张量表述 (G???, G?2?) 之间的闭式映射。

关键结果是，在横向各向同性和四阶取向矩的最小二次闭包下，从微观结构平均导出的四阶弹性张量变得可以表示为与构织张量弹性相同的结构形式。

式(14)的推导

假设关于 x? 轴的横向各向同性。二阶微观结构张量为：

M?? = diag(m?, m?, m?)

对应的构织张量为：

A?? = diag(a?, a?, a?)

基于构织的弹性张量由下式给出：

C????^F = λ A?? A?? + μ(A?? A?? + A?? A??)

对于横向各向同性构织张量，相关的刚度分量为：

C??^F = λ a?2
C??^F = (λ + 2μ) a?2
C??^F = μ a? a?

微观结构弹性张量给出：

C??^M = α(3/2 m?) + β(15/2 m?2 ? 3/4 m?)
C??^M = α(3 m?) + β(15/2 m?2 ? 3/2 m?)
C??^M = β(15/2 m? m?)

其中 α = E/(1+ν) 且 β = E(1+ν)/(1?2ν)

数值算例

考虑两个代表性数值算例以说明所提出的微观结构-构织映射的适用性和局限性。第一个算例代表强烈不相容的微观结构，第二个算例接近相容条件选取。

结论

本工作发展了基于构织张量的弹性与微观结构平均弹性之间的统一理论联系。

主要发现总结如下：

1. 闭式等价：在横向各向同性和四阶取向矩的二次闭包下，微观结构导出的弹性张量精确简化为构织张量弹性的结构形式。

2. 直接微观结构-构织映射：获得了显式的解析表达式……

CRediT作者贡献声明

George Z. Voyiadjis：写作 – 审阅与编辑，写作 – 原始草稿，可视化，软件，资源，项目管理，调查，形式分析，概念化。
Peter I. Kattan：写作 – 审阅与编辑，写作 – 原始草稿，可视化，验证，软件，方法论，调查，形式分析，数据整理。

George Z. Voyiadjis | Peter I. Kattan