LS-DYNA材料的二次开发

《LS-DYNA材料的二次开发》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、LS-DYNA材料的二次开发ANSYS/LS-DYNA专题培训内容●二次开发环境●主程序及入口条件●开发材料的本构、子程序及求解输入文件描述●编译、运行新的求解器●开发Kelvin_voigt粘弹材料●用新材料模式做大变形分析LS-DYNA二次开发基于FORTRAN环境在PC和UNIX平台下都需要进行连接编译，生成新的求解器Pc平台需安装digitalvisualfortran5.0或Microsoftpowerstation4.0提供的资源包括：Ls-dyna.f用户自定义本构子程序Ls-dyna.lib静态连接库Ls-dyna.dspdigitalFORTRANworkspace

2、文件或MAKEFILE用于（包括主程序）Readme.txt说明文件二次开发环境在UNIX平台提供如下资源：Makefile执行编译批处理文件（？）object文件（内含主程序）dyn21.f用户定义本构子程序支持的平台：940.2版DECNECIBMHPSGISUN950d版COMPAQCRAYSGIIBMLINUXSUN用户平台需安装FORTRAN77二次开发如何实现？用户自定义的本构代替ls-dyna.f或dyn21.f中的相关本构描述LS-DYNA共提供10种*user_defined_material_model,由这些输入数据为自定义本构提供参数，完成分析在程序中使用的自

3、定义subroutine要和Jobname.K中指定的相同在一次分析中，用户最多可同时使用10种自定义材料本构主程序及入口条件c******************************************************************c

4、LS-DYNAmainprogramentry

5、c******************************************************************programlsdyna3dcalldyna3dstopendc***************************************

6、***************************入口条件参数传递ccm(1)=young'smodulusccm(2)=poisson'sratioCcm(n)=用户在点K中给定的新本构参数ceps(1)=x应变增量ceps(2)=y应变增量ceps(3)=z应变增量ceps(4)=xy应变增量ceps(5)=yz应变增量ceps(6)=zx应变增量c单元类型etype:ceq.“brick”实体单元ceq.“shell”壳单元ceq.“beam”梁单元cctime=当前时间cdt1=当前时间步长ccapa=纵向剪切缩减因子csig(1)=x应力csig(2)=y应力csig(

7、3)=z应力csig(4)=xy应力csig(5)=yz应力csig(6)=zx应力cchisv(1)=历史变量1chisv(2)=历史变量2…chisv(n)=历史变量n每个积分步、主程序提供如下这些已知量：6个应变增量可能涉及的历史变量hisv（n）单元类型的字符串当前时间当前时间步长用户在点K文件中给定如下参数：弹性模量波松比其它参数cm(n)Ls-dyna.f或dyn21.f应完成的工作：求出6个应力增量求出其它可能涉及的历史变量hisv（n）参数说明由主程序提供的所有参数基于单元坐标系，计算得到的应力显然如此，之后由主程序将其转换到整体坐标系所有的历史变量在初始调用子程

8、序时将置零能量计算完全由主程序完成子程序举例c******************************************************************c

9、user-definedsubroutineexample

10、c******************************************************************subroutineumat41(cm,eps,sig,hisv,dt1,capa,etype,time)cisotropicelasticmaterial(sampleusersubroutine)cvaria

11、blesccm(1)=young'smodulusccm(2)=poisson'sratiocceps(1)=localxstrainincrementceps(2)=localystrainincrementceps(3)=localzstrainincrementceps(4)=localxystrainincrementceps(5)=localyzstrainincrementceps(6)=localzxstrainincrementcsig(1