当一个含有用户子程序的模型被提交给ABAQUS分析时，正确的编译和链接命令应该被自动执行。对于不同的运行平台，ABAQUS正确的编译和链接命令默认存储在环境文件中，这个文件位于abaqus_dir/site目录下，这里的abaqus_dir是ABAQUS的安装目录。如果想将子程序提供给他人使用，但不希望他人看到子程序的源代码，在这种情况下，可将子程序的源码编译成obj文件提供给他人。这个目标文件可以运行子程序，但不会显示源代码。......
2025-09-30
ABAQUS有限元子程序开发：混合硬化塑性材料本构方程
【摘要】：本节给出混合硬化塑性本构的控制方程，用户子程序要根据这些控制方程来编写。混合硬化塑性本构的控制方程涉及两部分，分别是材料在屈服之前的弹性部分和屈服后的塑性部分。式在Jaumann率下的率形式为将式在共轴旋转框架下进行时间积分，得到增量形式的弹性本构方程：用户子程序中编写代码时依据式进行。如果弹性预测应力超过了屈服应力σY，材料就会发生塑性流动。
本节给出混合硬化塑性本构的控制方程，用户子程序要根据这些控制方程来编写。混合硬化塑性本构的控制方程涉及两部分，分别是材料在屈服之前的弹性部分和屈服后的塑性部分。
1.弹性部分的控制方程
全量形式的弹性本构方程为
式中，上标el表示弹性。
式（5.46）在Jaumann率（共轴旋转框架下）下的率形式为
将式（5.47）在共轴旋转框架下进行时间积分，得到增量形式的弹性本构方程：
用户子程序中编写代码时依据式（5.48）进行。
2.塑性部分的控制方程
von Mises屈服面函数为
式中，Sij——积分点上偏应力张量的分量；
αij——背应力张量的分量；
σY——材料的单轴等效屈服应力，是一个输入参数。
当材料某个积分点上的应力满足式（5.49）时，该点进入屈服状态。von Mises屈服面是偏应力空间中的一个圆柱，其半径为(https://www.chuimin.cn)
对于运动硬化的塑性模型，其屈服面的半径R是一个常数。von Mises屈服面的法线可以写为
定义等效塑性应变率如下：
式中，上标pl表示塑性。
定义塑性流动法则如下：
背应力张量通过Prager-Ziegler线性动态硬化法则进行更新：
式（5.50）～式（5.54）决定了材料的塑性演化行为。但是，我们无法根据这几个式子直接进行编程，需要先进行时间积分，将其转换成可编程的增量形式。积分程序如下：首先进行弹性预测，也就是基于纯的弹性行为来计算等效应力，计算公式如下：
式中，上标o表示上一步的计算结果；上标pr表示预测。
经过一些推导，可以得到封闭形式的等效塑性应变增量的表达式，如下：
从而可以得到背应力、应力张量和塑性应变的更新表达式：
此外，还可以同时从下面的表达式中得到材料的一致雅可比矩阵：
上面的这种算法通常被称为弹性预测-径向返回算法[23]，因为在有效塑性加载条件下，对试验应力的校正会使得应力状态沿着由从屈服面中心向量所定义的方向返回到屈服面上。关于这一算法更加详细的公式推导请参考ABAQUS用户子程序手册1.2.22[1]。
相关文章
- 详细阅读
-
Fortran语言程序开发及应用在ABAQU详细阅读
Fortran语言主程序的语法示例如下：接下来，以一个简单的Fortran语言程序为例，介绍Fortran语言程序的基本结构。完整的程序代码如下：在上面的Fortran语言程序中，代码后添加了相应的注释，这有助于我们很好地理解程序的逻辑和流程。除了字符串常量外，Fortran语言不区分字母的大写和小写。对此，很多习惯使用其他程序语言的编写者在转而使用Fortran语言编程时会有不适应，导致耗费很多时间用于调试一些低级错误。......
2025-09-30
-
颗粒流本构关系的探讨详细阅读
采用分子混沌假设，并认为颗粒之间以二体碰撞为主。上述本构关系式由Gidaspow 针对气固两相流率先提出，在气固两相流模拟中获得了较广泛的应用。因此，上述本构关系严格地说，仅适用于速度梯度较小、颗粒近乎弹性的颗粒流。......
2025-09-29
-
使用ABAQUS有限元程序进行模拟分析详细阅读
ABAQUS 是一套功能强大的工程模拟的有限元软件［6～8］，其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。ABAQUS 被广泛地认为是功能最强的有限元软件，可以分析复杂的固体力学结构力学系统，特别是能够驾驭庞大复杂的问题和模拟高度非线性问题。由于ABAQUS 优秀的分析能力和模拟复杂系统的可靠性，其被各国的工业和研究中所广泛采用。ABAQUS 产品在大量的高科技产品研究中都发挥着巨大的作用。ABAQUS 主要分析功能如表2－1 所示。......
2025-09-29
-
二次开发黏聚裂纹子程序详细阅读
有限元法是目前发展比较完善的数值分析方法，目前常见的通用软件ABAQUS、LS－DYNA、ANSYS 提供了很多应用程序接口和用户子程序接口，方便用户根据专业问题定义自己所需的模型，并通过二次开发建立所需要的本构模型。黏聚裂纹子程序是通过LS－DYNA 二次开发实现的。根据第3 章介绍的黏聚裂纹模型理论，有限元程序在每一时间步提供的应变增量为Δεa，在每次计算前需要提取前一步存储的历史变量。表4－1UMAT 中的子程序......
2025-09-29
-
基于Chapman-Enskog法的颗粒相本详细阅读
在稠密气体分子运动的动理学理论 的分析方法中，主要有Chapman-Enskog迭代法和Grad的13矩方法。首先从Chapman-Enskog法来考察颗粒相脉动量的本构关系，由于问题的复杂性，仅探讨了颗粒近乎弹性 的情形，并与现有的近乎弹性的快速颗粒流理论成果作比较。......
2025-09-29
-
水体分层影响下的垂向混合模拟：控制方程解析详细阅读
为了有效模拟水体分层对垂向混合强度的影响，EFDC模型采用阶的Mellor-Yamada紊流模型［6，7］来求解垂向紊动扩散项，具体方程如下：式中：q 2/2为紊动强度；l为紊动混合长尺度；B 1、E 1、E 2和E 3为经验常数；Qq、Ql为源、汇项。垂向扩散系数Aq一般取和动量方程中的垂向扩散系数A v相等，φv和φb为稳定性函数，反映垂向密度分层对垂向混合的促进和抑制作用。......
2025-09-29
-
材料本构关系与流动应力实验的理论基础详细阅读
为了较真实地模拟材料的切削过程，应把材料的流动应力视为应变、应变率和温度的函数关系，表示为：。表2-1 工件化学成分 描述材料流动应力的模型通常有多种形式。Johnson-Cook本构关系形式简单、待定系数少、适应性强，是应用最为广泛的一种模型[103]。使用较多的流动应力实验方法主要有Hopkinson高速冲击法[104]和正交槽铣法[105]。根据上述情况和国内实验资源，拟采用Hopkinson高速冲击法进行测试，以获取材料的JC本构模型。......
2025-09-29
相关推荐