1 基本概念
1.1 定义与作用
本构方程是描述材料内部响应与外部作用之间关系的数学表达式，通常把应力、应变、温度、时间以及加载历史等变量联系起来。它的核心任务，是用可计算、可分析的形式刻画材料在受载条件下如何变形、传力和演化。
在工程与材料科学中，本构方程具有基础性地位。它既用于解释材料的弹性、塑性、黏弹性等行为，也用于支撑结构分析、工艺模拟和失效预测。许多数值方法之所以能够反映真实材料的力学特征，关键就在于本构模型是否合理。
1.2 发展历史
本构理论的发展，起初主要来源于经典弹性理论。胡克定律的提出，为应力与应变之间的线性关系提供了早期框架。随着材料种类增多和工程问题复杂化，单一线性关系已无法满足需要，于是塑性理论、黏弹性理论和损伤力学等分支逐步形成。
进入计算机辅助工程时代后，本构模型的发展进一步加快。研究者开始将实验数据、理论推导与数值计算结合起来，构建更适合复杂工况的模型。现代本构研究不仅关注宏观力学行为，也重视微观结构、加载路径和多场耦合作用。
1.3 适用范围与局限性
本构方程适用于多种材料体系，如金属、聚合物、陶瓷、复合材料以及岩土等。不同材料的组织结构和变形机制不同，因此其本构形式也会有明显差异。对于小变形、低速加载和温和环境条件，较简单的模型往往就能取得较好效果。
但本构模型通常具有一定局限。首先，模型往往基于特定假设，离开适用条件后准确性会下降。其次，许多材料表现出路径依赖、尺度效应或强非线性特征，单一模型难以全面描述。此外，参数获取、模型选择和数值稳定性也会影响其实际应用效果。
2 力学基础量
2.1 应力
应力是材料内部单位面积上所承受的内力，是描述受力状态的基本量。根据力的作用方式，可分为法向应力和剪应力；从研究角度看，又常用主应力、等效应力等形式来简化问题。
应力不仅反映材料受载程度，也与材料是否进入屈服、是否发生破坏密切相关。在本构方程中，应力通常作为响应量之一，与应变、温度等量共同决定材料状态。
2.2 应变
应变用于描述材料在受力后发生的相对变形程度。对于小变形问题，常采用线应变、剪应变等近似描述；当变形较大时，则需引入有限应变概念，以更准确地反映几何变化。
应变并不只是位移的简单变化，还包含变形方向、变形幅度以及局部形变分布等信息。不同的应变定义适用于不同的理论框架，是建立本构关系时不可缺少的基础变量。
2.3 变形与位移
位移描述物体各点相对于初始位置的移动情况，而变形则更强调材料内部形状和尺寸的改变。位移场是求取应变场的基础，二者在连续介质力学中紧密相关。
在工程分析中，位移常用于边界条件和结果评估，例如结构挠度、沉降或翘曲。变形则更直接地关联材料的局部响应，反映其是否出现拉伸、压缩、剪切或扭转等行为。
2.4 温度与时间变量
温度和时间是本构关系中常见的重要变量。温度会影响材料的刚度、强度、膨胀和相变行为；时间则常与蠕变、松弛、老化等现象相关。
在很多材料中，响应并非只取决于当前载荷，还与加载历程有关。因此，本构方程常需要引入时间作为状态变量，用以描述材料的历史依赖性。对于受热环境中的材料，温度变化还可能引发额外的热应力和结构演化。
3 本构关系的基本形式
3.1 线性本构方程
线性本构方程假设材料响应与外部作用之间近似成正比关系，常用于小变形、低应变和近似均匀材料。由于形式简洁、便于求解，它在工程上应用十分广泛。
这类模型的优点是可解释性强、参数较少、计算稳定性较好；不足之处在于难以描述大变形、强非线性和复杂路径效应。
3.1.1 胡克定律
胡克定律是线性弹性理论的经典表达，说明应力与应变在一定范围内呈正比关系。对于拉伸、压缩和剪切等基本变形模式，它给出了最早、最直观的定量描述。
在理想情况下，材料卸载后能够恢复原状，胡克定律可用来刻画这种可逆行为。它是许多更复杂理论的基础，也是工程力学中的重要起点。
3.1.2 广义胡克定律
广义胡克定律将简单的一维比例关系推广到三维应力状态，通常通过弹性常数矩阵把应力与应变联系起来。对于各向同性材料，表达式可进一步简化；对于各向异性材料，则需要更多独立参数。
这种形式适合描述实际工程材料在小变形条件下的弹性响应，也是有限元分析中常见的基础模型之一。
3.2 非线性本构方程
非线性本构方程用于描述应力与应变之间不再呈简单比例的情况。材料在大变形、高温、高速加载或复杂环境下，常表现出明显的非线性特征。
这类模型能够更真实地反映屈服、硬化、软化和损伤演化等过程，但建模与求解也更复杂，通常需要更多实验数据支撑。
3.3 各向同性与各向异性形式
各向同性本构模型假设材料在各个方向上的性质一致，因此其参数与方向无关。这种假设适用于组织较均匀、结构近似随机的材料。
各向异性本构模型则考虑材料性能随方向变化的情况，常见于晶体材料、纤维增强复合材料和轧制板材等。此类模型更贴近真实材料，但参数识别和理论推导都更复杂。
3.4 微分形式与积分形式
微分形式本构关系强调变量变化率之间的局部联系，常用于描述瞬时响应和增量更新过程。它在数值计算中较为常见，便于逐步积分。
积分形式则从整体历史出发，将当前状态与过去加载历程联系起来，更适合描述具有记忆效应的材料，如黏弹性介质。两种形式在不同理论和计算框架下各有优势。
4 典型材料本构模型
4.1 弹性材料模型
弹性材料模型描述材料在卸载后能够恢复原始形状的行为。其重点在于应力与应变之间的可逆关系，通常不考虑永久变形。
4.1.1 理想弹性模型
理想弹性模型假定材料在加载和卸载过程中完全可逆，不产生耗散和残余应变。它常用于近似描述小范围内的弹性响应。
这种模型简单明了，适合作为理论分析和初步计算的基础，但对真实材料的适应范围有限。
4.1.2 超弹性模型
超弹性模型主要用于大变形弹性材料，如橡胶类材料。它以应变能函数为核心，通过能量导数获得应力响应。
该模型可以较好描述大拉伸和大剪切条件下的力学行为，因而在柔性结构、生物组织和弹性体分析中应用广泛。
4.2 塑性材料模型
塑性材料模型用于描述材料在超过某一阈值后产生不可恢复变形的过程。其关键特征是存在屈服现象和路径依赖性。
4.2.1 屈服准则
屈服准则用于判断材料何时从弹性状态进入塑性状态。常见准则会将应力状态与某一临界条件进行比较，从而确定材料是否发生屈服。
它是塑性理论的起点，也是塑性计算中最重要的判据之一。
4.2.2 加工硬化模型
加工硬化模型描述材料在塑性变形过程中强度逐渐提高的现象。随着塑性应变增加，继续变形所需的应力可能上升。
这类模型能更真实地反映金属在成形和服役中的强化过程，对预测后续承载能力很有帮助。
4.2.3 理想塑性模型
理想塑性模型假设材料一旦屈服，便以近似恒定应力继续变形，不再体现硬化效应。它是分析塑性流动的简化形式。
虽然该模型较粗略，但在某些近似分析和教学场景中仍具有典型意义。
4.3 黏弹性材料模型
黏弹性材料兼具弹性恢复和黏性耗散特征，响应往往与时间有关。聚合物、沥青和部分生物材料常呈现这种行为。
4.3.1 Maxwell模型
Maxwell模型由弹簧与黏壶串联构成，适合描述应力松弛现象。它能体现材料在持续载荷下逐渐流动的特征。
该模型在理论上较为经典，但对某些长期恢复行为的描述能力有限。
4.3.2 Kelvin-Voigt模型
Kelvin-Voigt模型由弹簧与黏壶并联组成，适合刻画蠕变过程。它能反映材料在恒定载荷下逐步变形的趋势。
这一模型常用于近似描述迟滞响应，但在瞬时加载条件下并不总是足够准确。
4.3.3 标准线性固体模型
标准线性固体模型综合了弹性与黏性元件，能够同时描述蠕变和松弛行为。它在实际应用中比单一元件模型更灵活。
由于兼顾了较多现象，这一模型常被视为黏弹性分析中的基础框架之一。
4.4 损伤与断裂模型
损伤模型用于描述材料内部微裂纹、孔隙或界面退化导致的承载能力下降。断裂模型则进一步关注裂纹扩展、失稳和最终破坏过程。
这类模型强调材料性能随载荷与时间不断劣化，常用于评估寿命、疲劳和失效风险。它们在高可靠性结构设计中十分重要。
5 热-力耦合本构
5.1 热膨胀效应
材料受热后通常会发生体积或长度变化，这种现象称为热膨胀。若受约束，热膨胀还会转化为附加应力。
在本构关系中，热膨胀常作为温度引起的附加应变项加入，用于解释温度变化带来的变形偏移。
5.2 热弹性理论
热弹性理论研究温度变化与弹性变形之间的耦合关系。它适用于温度场变化不剧烈、材料仍保持弹性响应的情形。
该理论在精密仪器、热环境结构和热致变形分析中应用较多，可帮助预测温差引起的变形与应力集中。
5.3 热塑性理论
热塑性理论用于描述高温条件下材料的塑性流动行为。随着温度升高，材料屈服强度、硬化能力和流动特性可能明显改变。
在热加工和成形模拟中，热塑性模型通常是不可缺少的组成部分，因为它直接关系到成形载荷、组织变化和最终尺寸精度。
5.4 传热与力学响应耦合
传热与力学耦合意味着温度场和应力应变场彼此影响。材料变形会改变热传导条件，而温度变化又会反过来影响材料刚度和强度。
在一些工况下，这种耦合效应十分显著，例如快速加热、局部冷却或周期性热载荷条件下，不能把热学与力学完全分开处理。
6 连续介质力学中的表述
6.1 张量形式
本构关系在连续介质力学中通常以张量形式表示，以便统一描述多方向、多分量的物理量。张量表示法能够自然处理三维应力和应变状态。
这种表述方式具有良好的数学严谨性，也便于与坐标变换、主方向分析和数值计算结合。
6.2 主方向与本征值问题
主方向是指张量在某些特定方向上的分量可简化为对角形式的方向，对应的量称为主值。通过求解本征值问题，可以找出应力或应变的主方向与主大小。
在工程分析中，主应力和主应变常用于判断材料是否可能发生开裂、滑移或屈服，因为它们更直接反映局部极值状态。
6.3 坐标变换与不变量
坐标变换用于描述同一物理状态在不同参考系中的表示方式。一个合理的本构模型应尽量满足坐标无关性，也就是其物理结论不应依赖坐标选取。
不变量是张量在坐标变换下保持不变的量，常用于构造客观、稳定的本构公式。借助不变量，可以使模型更符合物理规律。
7 参数识别与实验测定
7.1 拉伸与压缩试验
拉伸与压缩试验是获取本构参数的基础手段之一。通过测量载荷与变形关系，可以提取弹性模量、屈服强度、硬化参数等信息。
这类试验直观、操作成熟，常作为材料性能评定和模型标定的首选方法。
7.2 剪切与扭转试验
剪切与扭转试验用于研究材料在切向载荷作用下的力学响应。它们对评估剪切模量、剪切屈服特性和抗扭性能具有重要意义。
对于某些薄壁构件、轴类零件和层状材料，这类试验尤其有价值。
7.3 松弛与蠕变试验
松弛试验通常保持应变恒定，观察应力随时间的变化；蠕变试验则保持应力恒定，观察应变随时间的发展。二者是识别黏弹性参数的重要手段。
通过这些试验，可以较好地确定材料的时间依赖性，并验证模型对长期行为的预测能力。
7.4 反演识别方法
反演识别方法是借助实验结果、优化算法和数值模拟，反向求取本构参数的技术。它特别适用于直接测量困难或参数之间耦合较强的情况。
这类方法能够提高模型拟合精度，但通常对初值、噪声和计算成本较为敏感。
8 数值实现
8.1 有限元中的本构积分
在有限元计算中，本构积分是指在每个增量步内更新材料状态变量，并计算新的应力响应。对于非线性材料，这一步直接决定模拟结果的可靠性。
本构积分需要同时兼顾精度、效率和稳定性，因此常采用增量式更新策略。
8.2 返回映射算法
返回映射算法是塑性与相关非线性本构计算中的常用方法。它用于将试算应力“映射”回满足屈服条件的可行域，从而得到修正后的应力状态。
该算法在工程计算中应用广泛，因为它结构清晰、数值性能较好，适合大规模有限元程序。
8.3 稳定性与收敛性
数值实现不仅要给出结果，还要保证计算过程稳定收敛。若本构更新过于刚性、步长过大或参数设置不当，程序可能出现发散或振荡。
因此，在材料模拟中，算法设计、步长控制和一致切线刚度的处理都十分关键。
8.4 计算程序中的材料子程序
材料子程序是将本构模型嵌入通用计算平台的常用方式。用户可通过编写子程序定义材料响应、状态变量更新和切线矩阵输出。
这种方式提高了模型扩展性，便于将实验拟合结果或自定义理论直接用于复杂结构分析。
9 应用领域
9.1 结构设计与安全评估
本构方程在结构设计中用于预测构件承载能力、变形水平和失效风险。通过合理的材料模型，可以更准确地评估安全裕度和使用寿命。
在桥梁、机械设备、压力容器和航空构件等领域，本构分析是设计验证的重要组成部分。
9.2 材料成形与加工模拟
在轧制、锻造、挤压、冲压和增材制造等过程中，本构模型可用于预测材料流动、应力分布和最终成形质量。它有助于优化工艺参数，减少缺陷。
成形模拟中常常涉及大变形、摩擦、温升和组织演化，因此对本构模型的准确性要求较高。
9.3 复合材料分析
复合材料由不同组分构成，通常具有明显的方向性和界面效应。本构方程用于描述其层合行为、各向异性特征以及损伤扩展过程。
在复合材料设计中，合理的本构模型有助于分析纤维取向、铺层顺序和局部失效对整体性能的影响。
9.4 微纳尺度材料研究
在微纳尺度下，材料行为可能不再遵循传统宏观假设，表面效应、尺寸效应和离散结构影响会更加突出。相应地，本构模型也需要进行修正或重建。
这类研究常见于微电子器件、薄膜、纳米结构和生物微结构材料中，对高精度建模提出了更高要求。
10 相关问题
10.1 本构方程的唯一性
本构方程并不总是唯一的。对于同一材料与同一现象，往往可以构建多种不同形式的模型，只要它们在一定范围内都能拟合实验结果。
因此，模型选择通常取决于研究目标、适用范围、参数可得性以及数值实现难度，而不只是形式是否“唯一正确”。
10.2 本构关系的客观性
客观性要求本构关系不应依赖观察者所选取的参考系或坐标表达方式。也就是说，材料的物理响应应是坐标无关的。
这一要求是本构理论的重要原则之一。若模型不满足客观性，其预测结果可能随坐标变换而改变，从而失去物理意义。
10.3 尺度效应与多尺度建模
尺度效应是指材料在不同尺寸层级上表现出不同力学特性的现象。随着特征尺寸减小，微结构、缺陷分布和界面影响往往变得更加显著。
多尺度建模试图把微观机制与宏观响应联系起来，通过不同层级的模型协同，提升对复杂材料行为的解释能力。
10.4 经验模型与机理模型的比较
经验模型主要依赖实验数据拟合，优点是简洁、实用、易于标定；机理模型则尽量从物理机制出发，强调结构、演化和因果关系。
两者各有优势。经验模型适合快速应用，机理模型更有助于解释现象和外推预测。实际工程中，常需要结合二者，以获得平衡的准确性与可用性。