资源描述：
《LS-DYNA材料的二次开发.ppt》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、LS-DYNA材料的二次开发ANSYS/LS-DYNA专题培训内容●二次开发环境●主程序及入口条件●开发材料的本构、子程序及求解输入文件描述●编译、运行新的求解器●开发Kelvin_voigt粘弹材料●用新材料模式做大变形分析LS-DYNA二次开发基于FORTRAN环境在PC和UNIX平台下都需要进行连接编译，生成新的求解器Pc平台需安装digitalvisualfortran5.0或Microsoftpowerstation4.0提供的资源包括：Ls-dyna.f用户自定义本构子程序Ls-dyna.lib静
2、态连接库Ls-dyna.dspdigitalFORTRANworkspace文件或MAKEFILE用于（包括主程序）Readme.txt说明文件二次开发环境在UNIX平台提供如下资源：Makefile执行编译批处理文件（？）object文件（内含主程序）dyn21.f用户定义本构子程序支持的平台：940.2版DECNECIBMHPSGISUN950d版COMPAQCRAYSGIIBMLINUXSUN用户平台需安装FORTRAN77二次开发如何实现？用户自定义的本构代替ls-dyna.f或dyn21.f中的相关
3、本构描述LS-DYNA共提供10种*user_defined_material_model,由这些输入数据为自定义本构提供参数，完成分析在程序中使用的自定义subroutine要和Jobname.K中指定的相同在一次分析中，用户最多可同时使用10种自定义材料本构主程序及入口条件c******************************************************************c
4、LS-DYNAmainprogramentry
5、c************************
6、******************************************programlsdyna3dcalldyna3dstopendc******************************************************************入口条件参数传递ccm(1)=young'smodulusccm(2)=poisson'sratioCcm(n)=用户在点K中给定的新本构参数ceps(1)=x应变增量ceps(2)=y应变增量ceps(3)=z应变增量ceps(4
7、)=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(3)=z应力csig(4)=xy应力csig(5)=yz应力csig(6)=zx应力cchisv(1)=历史变量1chisv(2)=历史变量2…chisv(n)=历史变量n每个积分步、主程序提供如下这些已知量：
8、6个应变增量可能涉及的历史变量hisv（n）单元类型的字符串当前时间当前时间步长用户在点K文件中给定如下参数：弹性模量波松比其它参数cm(n)Ls-dyna.f或dyn21.f应完成的工作：求出6个应力增量求出其它可能涉及的历史变量hisv（n）参数说明由主程序提供的所有参数基于单元坐标系，计算得到的应力显然如此，之后由主程序将其转换到整体坐标系所有的历史变量在初始调用子程序时将置零能量计算完全由主程序完成子程序举例c*******************************************
9、***********************c
10、user-definedsubroutineexample
11、c******************************************************************subroutineumat41(cm,eps,sig,hisv,dt1,capa,etype,time)cisotropicelasticmaterial(sampleusersubroutine)cvariablesccm(1)=young'smodulusccm(
12、2)=poisson'sratiocceps(1)=localxstrainincrementceps(2)=localystrainincrementceps(3)=localzstrainincrementceps(4)=localxystrainincrementceps(5)=localyzstrainincrementceps(6)=localzxstrainincrementcsig(1
