2. 各向异性来源与材料分类

各向异性不是材料的一种“附加属性”,而是其内部组织结构在宏观力学响应上的直接投影。要建立可靠的本构模型,第一步不是急着写刚度矩阵,而是回答两个问题:各向异性从何而来,以及它属于哪一类对称性。前者决定模型需要哪些内变量或方向参数,后者决定独立弹性常数有几个、刚度矩阵长什么样。

一、各向异性的四类主要来源

晶体织构。金属与陶瓷的单晶本身即各向异性,例如面心立方铜单晶的弹性模量沿〈111〉方向约为〈100〉方向的2.9倍。多晶体若晶粒取向随机,宏观上近似各向同性;但经过轧制、拉拔、挤压后,晶粒会形成择优取向,即织构。织构程度通常用取向分布函数(ODF)或极图描述,工程上常用织构系数或Lankford系数r值间接表征。轧制铝板的r值在0°、45°、90°方向分别为0.6、0.9、0.7,这种面内差异正是塑性各向异性的来源。

纤维增强。连续纤维复合材料中,纤维承担主要轴向载荷,基体传递剪力和横向应力。单向板的纵向模量E₁由混合律近似给出:E₁ ≈ V_f·E_f + (1−V_f)·E_m;而横向模量E₂受基体控制,通常只有E₁的1/10~1/20。以T300碳纤维/环氧为例,V_f=0.6时E₁≈135 GPa,E₂≈9 GPa,G₁₂≈5 GPa,ν₁₂≈0.30。这种“一强两弱”的极端差异,使单向板天然具有横观各向同性。

层合结构。将不同铺层角的单层板按顺序叠合,各向异性被“设计”出来。层合板的面内刚度由经典层合理论给出:

A_ij = Σ Q̄_ij^(k) · (z_k − z_{k−1})
D_ij = (1/3) Σ Q̄_ij^(k) · (z_k³ − z_{k−1}³)

其中Q̄_ij是第k层在整体坐标系下的折算刚度。铺层顺序改变,A、D矩阵随之改变,因此层合板可以呈现正交各向异性、一般各向异性甚至耦合拉伸—弯曲。这正是“结构各向异性”与“材料各向异性”的叠加。

生物组织。骨、肌腱、软骨、血管壁等由胶原纤维、弹性蛋白和矿物相按特定方向排列而成。密质骨沿骨干轴向的弹性模量约17~20 GPa,横向约10~12 GPa,呈现明显的正交各向异性;肌腱几乎只沿纤维方向承载,可近似为横观各向同性;而心脏瓣膜等软组织在大变形下还会出现纤维重排,各向异性随应变演化。

各向异性来源 → 对称性分类 → 独立常数 晶体织构 轧制/拉拔择优取向 ODF、极图、r值 面内塑性差异 纤维增强 E₁≫E₂,G₁₂低 V_f、纤维方向角 单向板横观各向同性 层合结构 铺层角、顺序、厚度 A、D矩阵耦合 可设计各向异性 生物组织 胶原纤维定向排列 骨、肌腱、血管壁 大变形下重排 对称性分类(按材料对称面数量) 各向同性 → 横观各向同性 → 正交各向异性 → 单斜 → 三斜 对称类型 对称面 独立常数 典型材料 各向同性 无穷多 2(E, ν) 随机取向多晶、铸态金属 横观各向同性 1个对称轴(面内无穷多) 5 单向板、肌腱、拉拔棒材 正交各向异性 3个互相垂直对称面 9 正交铺层板、密质骨、木材 单斜:1个对称面,13个独立常数;三斜:无对称面,21个独立常数

二、对称性分类与独立常数

广义胡克定律写作σ = C : ε,C是四阶刚度张量。三维情况下最多有21个独立常数,但材料内部对称性会大幅削减这个数目。判断对称性的操作方法是:找出材料中所有“力学等价”的方向,即沿这些方向做相同实验得到相同响应。

  • 各向同性:任意方向响应相同,只有2个独立常数,如E和ν。随机取向的多晶体、铸态金属接近此类。
  • 横观各向同性:存在一个特殊轴,垂直于该轴的平面内所有方向等价。独立常数为5个:E₁、E₂、ν₁₂、ν₂₃、G₁₂。单向纤维板、拉拔棒材、肌腱属于此类。
  • 正交各向异性:存在三个互相垂直的对称面,独立常数为9个:E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。正交铺层板、密质骨、木材是典型代表。
  • 单斜:仅有一个对称面,独立常数13个。当铺层角不是0°或90°、或存在单一剪切耦合时出现。
  • 三斜:无任何对称面,21个独立常数,一般各向异性。短纤维随机但带取向、复杂编织复合材料可能接近此类。
判断对称性时,先看材料是否存在一个“制造方向”或“生理方向”。纤维方向、轧制方向、骨干轴向通常就是特殊轴。若只有一个特殊轴,优先按横观各向同性建模;若存在三个正交的特殊方向,则按正交各向异性处理。

三、从来源到模型的对应关系

实际建模时,来源与分类之间并非一一对应,而是“来源决定方向,分类决定常数个数”。例如:

材料/结构 主要来源 推荐对称类型 需标定的常数
单向碳纤维/环氧板 纤维增强 横观各向同性 E₁、E₂、ν₁₂、ν₂₃、G₁₂
[0/90]s正交层合板 层合结构 正交各向异性 E_x、E_y、ν_xy、G_xy(面内)
轧制铝板 晶体织构 正交各向异性 0°、45°、90°拉伸曲线
密质骨 生物组织 正交各向异性 轴向、径向、环向模量
±45°缠绕管 层合+纤维 单斜或一般各向异性 拉伸—剪切耦合项

以单向板为例,若把1轴取为纤维方向、2轴取为面内垂直纤维方向、3轴取为厚度方向,则2-3平面是各向同性面,刚度矩阵中C₂₂=C₃₃,C₁₂=C₁₃,C₅₅=C₆₆,独立常数降为5个。若把铺层角改为30°,折算刚度Q̄中会出现Q̄₁₆、Q̄₂₆,拉伸与剪切耦合,对称性从正交降为单斜,独立常数增加到13个。这就是“方向定义改变对称性”的直接体现。

坐标系约定必须全文一致。建议:1轴沿纤维或轧制方向,2轴面内垂直,3轴为厚度法向。若中途旋转坐标系,必须重新计算Q̄并说明旋转角的正负定义,否则刚度矩阵和实验标定结果会对不上。

四、工程判断流程

面对一个新对象,可按以下顺序快速定位:

  1. 观察是否存在明显的定向组织:纤维、晶粒流线、胶原束、铺层方向。
  2. 确定特殊方向的数量:一个 → 横观各向同性;三个正交 → 正交各向异性;一个对称面 → 单斜;无 → 三斜。
  3. 用独立常数个数反推需要做哪些实验。横观各向同性至少需要5个独立测试,正交各向异性至少9个。
  4. 检查是否存在拉—剪耦合。若存在,说明对称性低于正交,需保留耦合项。
  5. 在模型中显式记录材料主轴与全局坐标系的夹角,作为本构矩阵旋转的依据。