基于ANSYS二次开发三剪统一强度理论弹塑性刚度矩阵验证.docVIP

基于ANSYS二次开发三剪统一强度理论弹塑性刚度矩阵验证 　　 摘要：利用ANSYS二次开发平台用户程序特性（UPFS），通过子程序USERMAT.F，编写出三剪统一强度理论的有限元子程序（SJMAT.F）。通过比较静水压力条件下的圆形隧洞的理论解和ANSYS二次开发计算出的数值解，验证三剪统一强度理论的理想弹塑性本构关系推导和子程序编写的正确性。 　　关键词：三剪统一强度理论ANSYS弹塑性本构关系 　　 　　一、三剪统一强度理论的ANSYS二次开发 　　1.ANSYS二次开发技术 　　众所皆知，任何一种计算机软件都具有局限性，它们只能应用于某一领域，都不可能完全解决客户的所有问题，而且许多问题是非常个性化的。为了使ANSYS软件能够更好地满足用户的需求，ANSYS给用户留有足够的开发空间和开发手段。对于ANSYS的二次开发环境由以下几个模块构成[1]。 　　1.1用户界面设计语言（UIDL） 　　1.2参数化程序设计语言（APDL） 　　1.3用户程序特性（UPFS） 　　1.4ANSYS数据接口 　　2.用户可编程特性（UPFS） 　　用户可编程特性（UPFS）允许用户使用自己的FORTRAN程序或者是C程序的功能。利用用户可编程特性（UPFS），我们可以将自己编写的FORTRAN程序和子程序连接到ANSYS，生成一个适和用户使用的特定的ANSYS程序版本。 　　目前，UPFS仅仅适用于以下几个产品系列：ANSYS 　　/Multiphysics、ANSYS/Structura、ANSYS/ PrepPost、ANSYS/Mechanical和ANSYS/University（研究版和大学版） 　　UPFS是一种非标准的使用方法，ANSYS公司对此部分的测试没法保证。所以用户必须对FORTRAN程序的正确负责。在调用自己的FORTRAN程序时，用户必须保证FORTRANA程序的正确性。并且用户所选择的VF软件版本必须与ANSYS软件版本相匹配。不然ANSYS软件无法调用FORTRAN程序，无法生成ANSYS.EXE文件。另外在并行的系统中使用ANSYS时，不允许使用用户可编程特性（UPFS）。 　　3.三剪统一强度理论的ANSYS二次开发 　　为了实现三剪统一强度理论的ANSYS二次开发的实现，本文利用ANSYS二次开发工具UPFS将三剪统一强度理论模型添加到ANSYS中。利用UPFS中的子程序USERMAT.F编写相关的弹塑性本构模型的Fortran程序，并将编写的Fortran程序进行编译、链接，从而实现计算的功能。而本文的核心内容是三剪统一强度理论的Fortran程序的编写，以及Fortran程序的编译、链接。 　　ANSYS在自定义本构模型部分留有用户自定义弹塑性本构模型子程序接口USERMAT.F[2]。该接口的数据输入量包括本载荷增量开始时的总应力、总应变、等效塑性应变、本载荷增量引起的应变增量以及单元号、材料号、当前子步及载荷步等基本信息。输出量为本载荷增量结束时高斯积分点处的总应力、等效塑性应变、弹塑性矩阵等。子程序USERMAT.F的参数和变量名必须和标准ANSYS程序源代码中的变量名一致。正文部分即用户根据所添加材料的本构关系编写的FORTRAN语句。 　　本文是在子程序USERMAT.F基础上进行编写的，对了解子程序USERMAT.F的相关内容是十分有必要的。 　　二、举例验证 　　本文对静水压力下理想圆形隧洞进行分析，通过采用圆形隧洞的理论解析解与通过三剪统一强度理论的子程序计算出的数值解析解进行对比、分析。 　　图2.1表示的是一个圆形隧洞计算模型。由于具有对称性，本文为了计算方便，取隧洞计算模型的1/4进行计算。它的隧道半径为1m，高为15m。表2.1表示的是隧洞计算模型的计算参数： 　　 　　图2.1 圆形隧洞计算模型图 　　表2.1 计算参数表 　　圆形隧洞的初始地应力的模拟是非常重要的。在ANSYS中，有以下两种方法可以实现初始地应力的模拟： 　　一种是：一般情况下只考虑到岩体的自重应力，但是忽略了构造应力，在第一步分析的时候，首先计算沿途的自重应力场。这种方法只需要给出岩体的各项参数就可以计算出结构内力，计算比较简单。但是，计算出来的应力场与实际的应力场有比较大的偏差，并且在自重作用下，岩体还会产生初始位移，在继续分析并且施工时，得到的位移结果是累加初始位移的结果，在实际过程中，它的初始位移早就结束，对隧洞的开挖没有影响，因此在以后的任何一个施工阶段分析位移场时，必须减去初始位移场。 　　另一种是：采用读入初始应力文件的方法。在结构分析的过程中，ANSYS中可以通过读入初始应力文件来把初始应力定义为一种荷载。可将初始地应力写成初