4. 材料属性定义：各向同性/各向异性材料、功能材料本构关系、温度/场依赖属性

做多物理场仿真这么多年，我最大的体会是：材料定义决定了仿真的下限。网格剖分得再漂亮，边界条件设得再精确，只要材料参数给错了，结果就是一堆漂亮的垃圾。今天咱们就聊聊材料属性这块，我踩过的坑不少，希望能帮你绕过去。

4.1 各向同性 vs 各向异性：别搞混了方向

先说说最基础的概念。各向同性材料，说白了就是材料在各个方向上性能都一样。比如普通钢材、玻璃、陶瓷，你从哪个方向拉它，弹性模量都差不多。这类材料定义起来最简单，一个标量就搞定。

但功能材料就没这么省心了。压电陶瓷、复合材料、单晶硅……这些都是各向异性的。什么意思？你从不同方向施加电场，它产生的应变完全不同。我记得有一次做PZT-5H压电片的仿真，客户给的参数表里只写了介电常数，我按各向同性算，结果谐振频率差了15%。后来才发现，压电材料的介电常数、弹性常数、压电常数都是张量，方向搞错了全白搭。

在COMSOL里定义各向异性材料时，我习惯用坐标旋转来对齐材料的主轴方向。比如压电陶瓷的极化方向通常是Z轴，但实际贴片时可能转了45度。这时候别傻乎乎地改矩阵数值，用旋转坐标系功能，省事又不容易出错。

4.2 功能材料的本构关系：压电、磁致伸缩、形状记忆

功能材料的本构关系，是仿真的灵魂。咱们挑几个常见的说说。

4.2.1 压电本构：电场和力场的耦合

压电材料的本构方程有两种形式——e型和d型。我个人更常用e型，因为它在有限元里数值稳定性更好。

// e型本构（应力-电荷形式）
T = cE * S - e^T * E
D = e * S + εS * E

// d型本构（应变-电压形式）
S = sE * T + d^T * E
D = d * T + εT * E

这里T是应力，S是应变，E是电场，D是电位移。cE是短路弹性刚度，e是压电应力常数，εS是夹持介电常数。

我曾经犯过一个低级错误：把压电常数e和d搞混了。e的单位是C/m²，d的单位是C/N，差了整整一个数量级。那次仿真出来的位移只有实测的十分之一，排查了两天才找到原因。嗯，现在每次定义材料我都会先检查单位。

4.2.2 磁致伸缩：磁场和力场的耦合

磁致伸缩材料（比如Terfenol-D）的本构和压电类似，只是把电场换成了磁场。核心方程是：

// 磁致伸缩本构
T = cH * S - e^T * H
B = e * S + μS * H

H是磁场强度，B是磁感应强度，μS是恒应变磁导率。你想想看，这和压电本构是不是长得一模一样？只是物理量换了个马甲。所以理解了压电，磁致伸缩也就顺理成章了。

4.2.3 形状记忆合金：热-力耦合的非线性

形状记忆合金（SMA）的本构就复杂多了，涉及马氏体相变。它不是线性的，而是有滞回环。我建议初学者先用简化的Brinson模型，把马氏体体积分数作为内变量，再逐步过渡到更复杂的Auricchio模型。

4.3 温度/场依赖属性：材料不是一成不变的

很多新手会忽略这一点：材料的属性是随温度和场强变化的。比如压电陶瓷的压电常数d33，在高温下会急剧下降，接近居里温度时甚至退极化。再比如磁致伸缩材料的饱和磁致伸缩系数，随磁场强度变化是非线性的。

我做过一个超声换能器的项目，工作温度从25℃升到80℃。一开始我用常温下的材料参数，仿真出来的谐振频率和实测差了2kHz。后来把弹性模量、介电常数都设成温度的函数，结果一下就吻合了。

4.3.1 温度依赖的常见形式

材料属性	温度依赖形式	典型参数
弹性模量E	线性：E(T)=E₀(1+αΔT)	α≈-3e-4/℃（金属）
介电常数ε	Curie-Weiss：ε∝1/(T-Tc)	Tc=130℃（PZT-5H）
压电常数d33	多项式拟合	d33(T)=d33₀(1-βT²)
热导率k	反比：k∝1/T	硅在300K时k≈150 W/mK

在COMSOL里定义温度依赖属性，我通常用解析函数或插值函数。如果实验数据点不多，用线性插值就够了；数据点密集的话，用样条插值更平滑。

4.3.2 场依赖属性：电场/磁场的影响

除了温度，电场和磁场也会改变材料属性。比如压电陶瓷在强电场下会出现电致伸缩效应，介电常数不再是常数。磁致伸缩材料在强磁场下会饱和，磁导率急剧下降。

定义场依赖属性时，我建议用分段函数。举个例子，磁致伸缩材料的磁导率μ随磁场H变化：

// 磁导率的分段定义
if (H < H_sat) {
    μ(H) = μ_initial + (μ_sat - μ_initial) * (H / H_sat)^2;
} else {
    μ(H) = μ_sat;
}

这里H_sat是饱和磁场强度，μ_initial是初始磁导率，μ_sat是饱和磁导率。这种分段定义在数值计算中很稳定，不会出现震荡。

4.4 知识体系总览

下面这张图是我自己梳理的材料属性定义框架，你可以对照着看：

这张图把材料属性定义的三大块串起来了。左边是各向同性/各向异性的选择，中间是功能材料的本构方程，右边是温度和场依赖的处理。三者缺一不可。

最后说一句：材料参数一定要和实验对标。仿真做得再花哨，没有实验验证就是空中楼阁。我每次拿到新材料，都会先做一个简单的单物理场验证——比如只加电场看位移，或者只加热看热膨胀。确认基础参数没问题了，再上多物理场耦合。这样出了问题也好排查，不至于一头雾水。

---

公众号:蓝海资料掘金营，微信deep3321