黄永刚教授的子程序的文档。原文档为英文，为方便初学者学习，翻译为中文

给出了正交各向异性材料在材料主轴上的本构关系 ，并导出了定向结晶材料和单晶材料的本构关系

本构模型matlab代码充足 AMPLE-实质性学习环境 AMPLE是MATLAB中材质点方法的准静态隐式实现。 可以从项目webapges获得有关AMPLE的更多信息： AMPLE是具有规则四边形背景网格的弹塑性大变形材料点代码（主要代码是ample.m）。 连续体框架基于更新的拉格朗日结构，其中包括两个不同的本构模型：线性弹性和具有von Mises屈服面的线性弹性完美塑性模型。 建议用户首先查看“ documentation”文件夹中的“ index.html”文件。 html文件已使用M2HTML（）进行了编译，它们应提供代码概述以及它们如何组合在一起。 这些脚本还应该与MATLAB的“帮助”功能兼容。 版本历史： AMPLE-原始的AMPLE版本 AMPLE_1.1-版本1.1，于2020年8月发布-专注于运行时的改进

目的为准确模拟地震作用下钢筋混凝土梁柱构件的弹塑性受力特性，解决ABAQUS软件显式分析模块缺乏用于三维梁单元的混凝土本构模型及钢筋本构模型的问题。方法采用的混凝土单轴本构模型考虑了箍筋对混凝土的约束影响以及混凝土裂面效应的影响，钢筋单轴本构模型合理考虑了Bauschinger效应，利用ABAQUS用户子程序接口VUMAT进行二次开发，编制了用于显式动力分析的梁单元混凝土、钢筋本构模型的计算程序，并对往复荷载下钢筋混凝土柱受力性能进行了数值模拟。结果通过与ABAQUS自带的钢筋本构模型对比，验证了采用笔者

二维平面应力问题中的弹性材料模型 [s，D] = Elastic_Model（材料，e） 输入： -------- 材料：包含材料属性 Material.E（弹性模量）和 Material.v（泊松比）的变量e: (3*1) 应变向量 输出： ----------- s: (3*1) 应力向量D: (3*3) 本构矩阵

本构模型matlab代码FSI求解器 用Matlab和C ++（使用Eigen和Spectra）编写的用于流体-结构相互作用的单片3D解算器。 使用Newmark方案及时求解方程，在此期间，使用弧长方法求解大变形弹性和流体动力学的非线性方程。 目录 快速总结 该代码旨在解决由以下耦合方程组成的3D流固耦合问题： 固体区的大变形弹性：假设材料是各向同性的，并且遵循St Venant-Kirchhoff本构方程； 流体域中的不可压缩的Navier-Stokes：我们假设牛顿不可压缩的流体。 使用任意的拉格朗日-欧拉公式来编写方程式，并且通过流体和固体之间的界面处的速度以及边界条件和初始条件的连续性方程式来封闭系统。 实现了一种基本的方法，该方法可在时间步长中移动ALE网格以跟随界面的移动。 获得新网格作为解决位移边界条件所定义的弹性问题的方法，以确保正确跟踪界面。 因此，通过一对弹性系数来控制失真。 该方法是整体的，在每个时间步都解决了一个单一的大型数值系统。 每个时间增量的更新规则由Newmark方案给出，它具有可调整的beta和gamma参数。 在Newmark方案的每次迭代中，必须解

部分Johnson-Cook+金属本构模

修正剑桥模型子程序，正常固结软粘土的本构模型

针对粘土剑桥模型，开发了显式应力积分算法，用于模拟应力应变关系

应用于J2流动理论—径向返回算法 9 应力更新算法 小应变时的弹—塑性本构关系和框5.6的J2 流动理论， 注意到塑性流动方向是在偏应力的方向，给出为 J2塑性流动理论基于Mises屈服面，它特别适用于金属塑性， 该模型的关键假设是压力对在金属中的塑性流动没有影响；屈服条件和塑性流动方向是基于应力张量的偏量部分。 它也是屈服表面的法向，即 在偏应力空间(π平面)，Mises屈服表面是环状，法向是径向。在塑性流动方向(径向)定义一个单位法向矢量为 * J2的特点是材料进行塑性后体积不变，塑性流动仅与形变有关