LSDYNA材料的二次开发.pptVIP

* ANSYS/LS-DYNA Training 2001 LS-DYNA材料的 二次开发 ANSYS/LS-DYNA专题培训 内 容 ● 二次开发环境 ● 主程序及入口条件 ● 开发材料的本构、子程序及求解输入文件描述 ● 编译、运行新的求解器 ● 开发Kelvin_voigt粘弹材料 ● 用新材料模式做大变形分析 LS-DYNA二次开发基于FORTRAN环境 在PC和UNIX平台下都需要进行连接编译，生成新的求解器 Pc平台需安装 digital visual fortran 5.0或Microsoft power station 4.0 提供的资源包括： Ls-dyna.f 用户自定义本构子程序 Ls-dyna.lib 静态连接库 Ls-dyna.dsp digital FORTRAN workspace 文件或MAKEFILE 用于（包括主程序） Readme.txt 说明文件 二次开发环境 在UNIX平台提供如下资源： Makefile 执行编译批处理文件（？） object 文件（内含主程序） dyn21.f 用户定义本构子程序 支持的平台： 940.2 版 DEC NEC IBM HP SGI SUN 950d 版 COMPAQ CRAY SGI IBM LINUX SUN 用户平台需安装FORTRAN77 二次开发如何实现？ 用户自定义的本构代替ls-dyna.f或dyn21.f中的相关本构描述 LS-DYNA共提供10种*user_defined_material_model,由这些输入数据为自定义本构提供参数，完成分析 在程序中使用的自定义subroutine要和Jobname.K中指定的相同 在一次分析中，用户最多可同时使用10种自定义材料本构 主程序及入口条件 c****************************************************************** c| LS-DYNA main program entry | c****************************************************************** program lsdyna3d call dyna3d stop end c****************************************************************** 入口条件 参数传递 c cm(1)=youngs modulus c cm(2)=poissons ratio C cm(n)=用户在点K中给定的新本构参数 c eps(1)= x应变增量 c eps(2)= y应变增量 c eps(3)= z应变增量 c eps(4)= xy应变增量 c eps(5)= yz应变增量 c eps(6)= zx应变增量 c 单元类型 etype: c eq.“brick” 实体单元 c eq.“shell” 壳单元 c eq.“beam” 梁单元 c c time=当前时间 c dt1=当前时间步长 c capa=纵向剪切缩减因子 c sig(1)= x应力 c sig(2)=y应力 c sig(3)=z应力 c sig(4)=xy应力 c sig(5)=yz应力 c sig(6)=zx应力 c c hisv(1)=历史变量1 c hisv(2)=历史变量2 … c hisv(n)=历史变量n 每个积分步、主程序提供如下这些已知量： 6个应变增量 可能涉及的历史变量hisv（n） 单元类型的字符串 当前时间 当前时间步长 用户在点K文件中给定如下参数： 弹性模量 波松比 其它参数cm(n) Ls-dyna.f或dyn21.f应完成的工作： 求出6个应力增量 求出其它可能涉及的历史变量hisv（n） 参数说明 ?由主程序提供的所有参数基于单元坐标系，计算得到的应力显然如此，之后由主程序将其转换到整体坐标系 ?所有的历史变量在初始调用子程序 时将置零 ?能量计算完全由主程序完成 子程序举例 c****************************************************************** c| user-defin