用户子程序UELMAT的接口和用户子程序UEL的接口相似,具体如下:与用户子程序UEL不同,用户子程序UELMAT可以直接调用ABAQUS内置的材料模型,这样我们可以在编程时重点考虑在单元层面的实现上,而无须考虑复杂的材料实现。为了在模型中使用用户子程序UELMAT,需要在模型的inp文件中添加接口,以指明哪些单元需要使用用户子程序UELMAT计算其单元刚度矩阵和右手边残差向量。......
2023-11-03
用户子程序USDFLD是用于Abaqus/Standard中的场变量定义子程序,其具有以下特性:
(1)允许将材料点处的场变量定义为时间函数或输出变量标识符表中列出的任何可用的材料点处的量的函数,但用户定义的输出变量UVARM和UVARMn除外。
(2)可用于引入解相关的材料特性,因为此类特性可以很容易地定义为场变量的函数。
(3)将在材料定义中包含了用户子定义场变量的单元的所有材料点上被调用。
(4)可以调用实用子程序GETVRM访问材料点的数据。
(5)可以使用和更新解依赖的状态变量(SDV)。
(6)可以与用户子程序UFIELD一起使用,以指定预定义的场变量。
用户子程序USDFLD的Fortran程序接口如下:
下面分别解释子程序的输入输出参数及其意义,这些参数主要分为以下三类。(www.chuimin.cn)
1)必须定义和更新的变量
field(nfield):场变量数组,包含了当前积分点上的所有场变量,它们传进来时是有值的,其值是当前增量步开始时节点上值的插值。更新后的场变量的值会传入以下子程序:CREEP,HETVAL,UEXPAN,UHYPEL,UMAT,UMATHT和UTRS。
2)可以定义和更新的变量
statev(nstatv):状态变量数组,包含了所有解依赖的状态变量(SDV),它们传进来时是有值的,其值是当前增量步开始时的状态变量的值。状态变量数组的长度(nstatv)通过关键字*DEPVAR进行设置。
pnewdt:时间比,建议的时间增量与当前使用的时间增量(dtime)的比值,如果该比值为1,则表明建议的时间增量与当前使用的时间增量一致。注意,这个变量的更新可以和ABAQUS输入文件中设置的自动更新时间增量一起使用,并不冲突。
3)只能使用的变量
除了上面提到的三个变量外,其他变量可被使用,但不能被定义和更新,否则会出现错误。
有关基于ABAQUS的有限元子程序开发及应用的文章
用户子程序UELMAT的接口和用户子程序UEL的接口相似,具体如下:与用户子程序UEL不同,用户子程序UELMAT可以直接调用ABAQUS内置的材料模型,这样我们可以在编程时重点考虑在单元层面的实现上,而无须考虑复杂的材料实现。为了在模型中使用用户子程序UELMAT,需要在模型的inp文件中添加接口,以指明哪些单元需要使用用户子程序UELMAT计算其单元刚度矩阵和右手边残差向量。......
2023-11-03
当一个含有用户子程序的模型被提交给ABAQUS分析时,正确的编译和链接命令应该被自动执行。对于不同的运行平台,ABAQUS正确的编译和链接命令默认存储在环境文件中,这个文件位于abaqus_dir/site目录下,这里的abaqus_dir是ABAQUS的安装目录。如果想将子程序提供给他人使用,但不希望他人看到子程序的源代码,在这种情况下,可将子程序的源码编译成obj文件提供给他人。这个目标文件可以运行子程序,但不会显示源代码。......
2023-11-03
ABAQUS有一个非常大的单元库,含有多种类型单元,可以满足各种复杂的力学分析。相比于写一个完整的有限元求解程序,在一些现有程序的基础上编写用户单元可以大大降低开发成本、缩短开发时间,并且可以充分利用ABAQUS提供的强大的前后处理能力。此外,ABAQUS内置的求解器的效率非常高,求解非线性问题具有很好的收敛性,这也使得用户单元子程序具有非常广的应用前景,是一个很好的助力科研和工程的工具。......
2023-11-03
而由Dharmasena等[90]提出的非典型应变率硬化的材料模型则需要通过编写用户子程序VUHARD来实现。采用用户子程序VUHARD实现上述非典型应变率硬化的材料模型代码如下:为了在模型中使用用户材料塑性硬化子程序VUHARD,需要在inp文件中定义如下语句:上述语句中,第6行指定了需要使用用户子程序VUHARD,第7行的5个数对应于程序中的props变量,其含义分别为屈服应力σY、参考应变率、指数m、切线模量Et和杨氏模量E。......
2023-11-03
本节给出混合硬化塑性本构的控制方程,用户子程序要根据这些控制方程来编写。混合硬化塑性本构的控制方程涉及两部分,分别是材料在屈服之前的弹性部分和屈服后的塑性部分。式在Jaumann率下的率形式为将式在共轴旋转框架下进行时间积分,得到增量形式的弹性本构方程:用户子程序中编写代码时依据式进行。如果弹性预测应力超过了屈服应力σY,材料就会发生塑性流动。......
2023-11-03
Fortran语言提供了5种内在的数据类型,也可以根据需要来自定义数据类型。此外,在声明整型变量时,还可以指定变量使用的字节数。在Fortran语言中,用两个连续的实数来存储复数的这两部分。......
2023-11-03
蠕变是固体材料在保持应力不变的条件下,应变随时间延长而增加的现象。由于蠕变,材料在某瞬时的应力状态一般不仅与该瞬时的变形有关,而且与该瞬时以前的变形过程有关[24]。虽然Abaqus/Standard提供了丰富的蠕变本构模型,但Abaqus/Explicit还没有关于蠕变的本构模型。然而,对于一些特定的问题,采用Abaqus/Standard很难求解收敛,此时不得不借助Abaqus/Explicit进行求解,如果其中涉及的材料具有较为明显的蠕变特性,则需要借助VUMAT在Abaqus/Explicit中实现蠕变本构模型[26]。......
2023-11-03
相较于Fortran语言,C++语言的一大优势是含有丰富的外部库,可以满足各种计算需求,能大大提高开发子程序的效率。第1章已经简单介绍了使用外部库的好处,并通过一个简单的例子说明了使用外部库可以提高代码的可读性。本节以数值计算库armadillo为例来介绍如何在C++语言中调用外部库。如果没有进行这一步骤,则在运行中会弹出报错消息框并退出,而且弹出的是BLAS库的消息框。不过,在armadillo库中使用一维数组对矩阵进行初始化时是列优先的。......
2023-11-03
相关推荐