ABAQUS有一个非常大的单元库，含有多种类型单元，可以满足各种复杂的力学分析。相比于写一个完整的有限元求解程序，在一些现有程序的基础上编写用户单元可以大大降低开发成本、缩短开发时间，并且可以充分利用ABAQUS提供的强大的前后处理能力。此外，ABAQUS内置的求解器的效率非常高，求解非线性问题具有很好的收敛性，这也使得用户单元子程序具有非常广的应用前景，是一个很好的助力科研和工程的工具。......
2023-11-03
C++调用外部库的方法ABAQUS有限元子程序开发及应用
【摘要】：相较于Fortran语言，C＋＋语言的一大优势是含有丰富的外部库，可以满足各种计算需求，能大大提高开发子程序的效率。第1章已经简单介绍了使用外部库的好处，并通过一个简单的例子说明了使用外部库可以提高代码的可读性。本节以数值计算库armadillo为例来介绍如何在C＋＋语言中调用外部库。如果没有进行这一步骤，则在运行中会弹出报错消息框并退出，而且弹出的是BLAS库的消息框。不过，在armadillo库中使用一维数组对矩阵进行初始化时是列优先的。
相较于Fortran语言，C＋＋语言的一大优势是含有丰富的外部库，可以满足各种计算需求，能大大提高开发子程序的效率。第1章已经简单介绍了使用外部库的好处，并通过一个简单的例子说明了使用外部库可以提高代码的可读性。本节以数值计算库armadillo为例来介绍如何在C＋＋语言中调用外部库。
通常情况下，只需要正常安装armadillo库就可以使用了。但如果想要使用一些好用的函数（如矩阵乘法、解方程的函数等），则需要引入BLAS库和LAPACK库。要想引入这两个库，就需要在链接器中提高它们的lib，具体操作步骤如下：
第1步，修改环境文件abaqus_v6.env中的link_sl变量，在其末尾添加两个文件名：blas_win64_MT.lib和lapack_win64_MT.lib。
第2步，将文件blas_win64_MT.lib和lapack_win64_MT.lib复制到Visual Studio（VS）的cpp库中（选择amd64），具体目录为C：/ProgramFiles（x86）/Microsoft Visual Studio 11.0/VC/lib/amd64。
第3步，将对应的两个.dll文件复制到system32目录中。如果没有进行这一步骤，则在运行中会弹出报错消息框并退出，而且弹出的是BLAS库的消息框。(www.chuimin.cn)
通过以上三步，就可以配置出非常方便使用且功能强大的C＋＋语言数值计算环境。在配置armadillo库后，就可以更加直观和方便地定义矩阵和向量，具体可以通过下面的代码来定义：
使用上面的代码定义时，需要在子程序接口中把DDSDDE变量（上面的程序中为DDSDDEori）声明为double*或void*类型。在第2条代码中，false表示不复制数据，直接用原始位置（内存地址）；true表示遵守strict-aliasing规则，即不添加空间，保证总是使用这一段空间。
注意：这里不需要对矩阵进行转置。这是因为，对于二维数组，虽然在Fortran中是列优先的，而在C＋＋中是行优先的。不过，在armadillo库中使用一维数组对矩阵进行初始化时是列优先的。
有关基于ABAQUS的有限元子程序开发及应用的文章
- 详细阅读
-
基本数据类型-基于ABAQUS的有限元子程序详细阅读
Fortran语言提供了5种内在的数据类型，也可以根据需要来自定义数据类型。此外，在声明整型变量时，还可以指定变量使用的字节数。在Fortran语言中，用两个连续的实数来存储复数的这两部分。......
2023-11-03
-
蠕变材料模型概述-基于ABAQUS的有限元子详细阅读
蠕变是固体材料在保持应力不变的条件下，应变随时间延长而增加的现象。由于蠕变，材料在某瞬时的应力状态一般不仅与该瞬时的变形有关，而且与该瞬时以前的变形过程有关[24]。虽然Abaqus/Standard提供了丰富的蠕变本构模型，但Abaqus/Explicit还没有关于蠕变的本构模型。然而，对于一些特定的问题，采用Abaqus/Standard很难求解收敛，此时不得不借助Abaqus/Explicit进行求解，如果其中涉及的材料具有较为明显的蠕变特性，则需要借助VUMAT在Abaqus/Explicit中实现蠕变本构模型[26]。......
2023-11-03
-
基于线程的并行计算：ABAQUS有限元子程序详细阅读
在基于线程的并行计算中，对于公共空间的变量或者共享的资源，需要进行额外的保护。它们的编号为1～10，开发者可以通过它们的编号进行使用。初始化互斥锁的最佳位置在子程序uexternaldb中。示例如下：经过上面代码的计算，1号互斥锁（＃1）就可以在其他要被调用的用户子程序中使用了，以保护对共享公用块中变量的访问，防止其被多个互相竞争的线程同时访问而造成破坏。......
2023-11-03
-
ABAQUS有限元子程序UELMAT接口详细阅读
用户子程序UELMAT的接口和用户子程序UEL的接口相似，具体如下：与用户子程序UEL不同，用户子程序UELMAT可以直接调用ABAQUS内置的材料模型，这样我们可以在编程时重点考虑在单元层面的实现上，而无须考虑复杂的材料实现。为了在模型中使用用户子程序UELMAT，需要在模型的inp文件中添加接口，以指明哪些单元需要使用用户子程序UELMAT计算其单元刚度矩阵和右手边残差向量。......
2023-11-03
-
子程序编程及调用方法详解详细阅读
子程序调用可通过辅助机能M98代码进行，在M98指令中，子程序号由地址P规定。如“N15 M98 P30100；”为调用子程序O0100三次，而“N15 M98 P3100；”则表示调用子程序O3100一次。......
2023-06-25
-
基于ABAQUS的有限元子程序开发与应用：M详细阅读
不同于基于线程的并行计算，MPI实现的并行是进程级的。基于MPI的用户子程序的并行计算，总是通过一些特定的函数来实现，先简单介绍一下这些函数：●MPI_Init：告知MPI系统进行所有必要的初始化设置。它写在启动MPI并行计算的最前面。以上只是表达了作为一个MPI并行计算的基本结构，并没有真正涉及进程之间的通信，为了更好地进行并行，必然需要在进程间通信。下面给出一个具体的基于MPI的并行计算的例子：......
2023-11-03
-
Fortran语言程序开发及应用在ABAQU详细阅读
Fortran语言主程序的语法示例如下：接下来，以一个简单的Fortran语言程序为例，介绍Fortran语言程序的基本结构。完整的程序代码如下：在上面的Fortran语言程序中，代码后添加了相应的注释，这有助于我们很好地理解程序的逻辑和流程。除了字符串常量外，Fortran语言不区分字母的大写和小写。对此，很多习惯使用其他程序语言的编写者在转而使用Fortran语言编程时会有不适应，导致耗费很多时间用于调试一些低级错误。......
2023-11-03
相关推荐