ADINA自定义材料开发原理与实际操作.pdfVIP

ADINA 自定义材料开发原理与实际操作 武汉大学 lynnxiong(十方) wanbiao-007 huaizhi1983 目录 1 前言 1 2 ADINA 调用子程序过程 1 2.1 荷载增量步与应变小增量步 1 2.2 非线性有限元的基本解法 2 2.3 ADINA 与子程序的实际调用过程 5 3 二次开发实际操作 8 3.1 自定义材料开发综述 8 3.2 开发环境 8 3.3 开发步骤 9 3.4 注意事项 9 参考文献 9 附录 10 ADINA 中二次开发 1 前言 关于ADINA 自定义材料的二次开发资料较少，其自带的英文理论手册，以及 长江水利委员会在二十世纪八十年代翻译的ADINA 用户手册中，虽然有一定的介 绍和解释，但大多都是一个荷载增量步的计算过程，很容易跟应力计算的小增量 步（subincrement ）混淆，给实际程序编制和调试带来困难。因此，在有限元计算 理论的基础上，结合实际程序调试过程，以及参考了ABAQUS 二次开发的一些资 料，我们在此尝试着给出整个计算调用过程的思路。需要说明的是，以下仅仅是 相关开发人员迄今为止所理解的一个计算过程，可能有错误，仅供大家参考。此 外，限于我们目前学习和研究的范围，在没有另外说明的情况下，下文的计算假 设均是小变形和小应变，且为了简化起见，计算中不包含如热应变等非弹性应变。 对于其它的情况，我们没有研究过，不便多做说明。 2 ADINA 调用子程序过程 2.1 荷载增量步与应变小增量步 首先对计算中的两种增量步做一个简介，一种是荷载增量步，就是在非线性 计算中经常看到的“增量法”中的增量，本质是以一段段的折线去近似曲线，这 种增量步的时间含义在ADINA 中就是step ，在资料中相关的表示是t，Δ t，在应 变方面的相关的表示是STRAIN ；另一种是小增量步，由一个step 的时间增量除以 某个整数以后所得的小增量，它的用途将在下文中叙述，在资料中相关的表示是 τ ，Δ τ 之类，在应变方面的表示是DEPS （子程序中），每个step 的小增量个数 在二次开发中是由用户输入的，在自定义材料中number of strain subdivision 处指 定 （下文记为INTER），ADINA 缺省为10。 1 ADINA 中二次开发 图 2-1 小增量步示意图 图 2-1 为ADINA 结构理论手册中的示意图，Δ eij 即是程序中用到的DEPS。按 照上述说明，两种增量步之间可以理解为以下关系： t S T R A I N = D E P =S INTER I N T E R 2.2 非线性有限元的基本解法 然后再来回顾下有限元法的求解过程。把要求解的区域介质分成有限个单元， 节点位移是基本的未知量。通过位移函数 （形函数），用节点位移唯一地表示单 元内任一点的应变；再利用广义虎克定律，用节点位移唯一地表示单元内任一点 的应力。根据能量原理，找到与单元内部应力状态等效的节点力，从而建立等效 节点力与节点位移的关系。然后将每一单元所承受的荷载，按静力等效原则移置 到节点上。在每一节点建立平衡方程，得到一个线性方程组：解出这个方程组， 求出节点位移，然后可求得每个单元的应力。 （顺便提一下积分点的概念，在求解刚度矩阵和节点荷载时，需计算单元区 域内的积分。但这种积分一般是很复杂的，通常难以用显式表示，往往都用数值 积分方法计算积分值，即在单元内选出某些点，称为积分点，求出被积函数在这 些点的值，然后根据这些数值求出积分值。数值积分有两类方法，一类方法积分 点是等间距的，