LS-DYNA材料的二次开发.docx

ANSYS/LS-DYNA专题培训

LS-DYNA材料的

二次开发

ANSYS/LS-DYNATraining2001

内容

●二次开发环境

●主程序及入口条件

●开发材料的本构、子程序及求解输入文件描述

●编译、运行新的求解器

●开发Kelvin_voigt粘弹材料

●用新材料模式做大变形分析

ANSYS/LS-DYNATraining2001

二次开发环境

LS-DYNA二次开发基于FORTRAN环境

在PC和UNIX平台下都需要进行连接编译，

生成新的求解器

Pc平台需安装digitalvisualfortran5.0或

Microsoftpowerstation4.0

提供的资源包括：

Ls-dyna.f用户自定义本构子程序

Ls-dyna.lib静态连接库

Ls-dyna.dspdigitalFORTRANworkspace文件或MAKEFILE用于（包括主程序）

Readme.txt说明文件

ANSYS/LS-DYNATraining2001

用户平台需安装支持的平台：

FORTRAN77

940.2版

DECNECIBM

HPSGISUN

在UNIX平台提供如下资源：

950d版

Makefile执行编译批处理文件（？）

COMPAQCRAYSGI

IBMLINUXSUN

object文件（内含主程序）

dyn21.f用户定义本构子程序

ANSYS/LS-DYNATraining2001

二次开发如何实现？

用户自定义的本构代替ls-dyna.f或dyn21.f中的相关本构描述

LS-DYNA共提供10种

*user_defined_material_model,由这些输入数据为自定义本构提供参数，完成分析

在程序中使用的自定义subroutine要和Jobname.K

中指定的相同

在一次分析中，用户最多可同时使用10种自定义材料本构

ANSYS/LS-DYNATraining2001

主程序及入口条件

c******************************************************************

c|LS-DYNAmainprogramentry|

c******************************************************************

programlsdyna3d

calldyna3d

stop

end

c******************************************************************

ANSYS/LS-DYNATraining2001

入口条件

csig(1)=x应力参数传递

csig(2)=y应力

ccm(1)=youngsmodulus

csig(3)=z应力

ccm(2)=poissonsratio

csig(4)=xy应力Ccm(n)=用户在点K中给定的新本构参数

ceps(1)=x应变增量csig(5)=yz应力

ceps(2)=y应变增量

csig(6)=zx应力

ceps(3)=z应变增量

c

ceps(4)=xy应变增量

ceps(5)=yz应变增量

chisv(1)=历史变量1

chisv(2)=历史变量2

ceps(6)=zx应变增量

…

c单元类型etype:

chisv(n)=历史变量n

ceq.“brick”实体单元

ceq.“shell”壳单元

ceq.“beam”梁单元

c

ctime=当前时间

cdt1=当前时间步长

ccapa=纵向剪切缩减因子

用户在点K文件中

给定如下参数：

弹性模量

Ls-dyna.f或dyn21.f应

完成的工作：

波松比

求出6个应力增量

其它参数cm(n)

求出其它可能涉及的历史变量hisv（n）

每个积分步、主程序提供

如下这些已知量：

6个应变增量

可能涉及的历史变量hisv（n）

单元类型的字符串

当前时间

当前时间步长

ANSYS/LS-DYNATraining2001

参数说明

?由主程序提供的所有参数基于单元坐标

系，计算得到的应力显然如此，之后由

主程序将其转换到整体坐标系

?所有的历史变量在初始调用子程序

时将置零

?能量计算完全由主程序完成

ANSYS/LS-DYNATraining2001

子程序举例

c******************************************************************

c|user-definedsubro