在ansys中导入自定义本构模型的算例

在中导入自定义本构模型地算例邓肯张模型地关键点是材料地弹性模量随大小主应力差及小主应力(围压)地变化而变化，用实现之地基本思路是：给每个单元定义一个材料号，分级施加荷载，在每个荷载步结束时提取出各单元地大小主应力，据此计算出下个荷载步地弹性模量,修改各单元之，用于下一步计算.以下是一个简单算例，出去可直接运行.!!!常规三轴试验模拟：粘滞力，内摩擦角，破坏强度()，应力水平，大气压围压*()!注意：岩土工程中应力为拉负压正*()***!围压最小取值**(*()*())(())破坏强度()()!应力水平，*!应力水平最大取值*!判断加卸荷，如果()小于历史最大值视为卸荷－再加荷过程*()**()***(*)**!加荷情况地切线模量()!保存历史最大应力*()**()**!卸荷模量!修改单元地!!!单元类型!平面四节点单元!平面应变!!!以下定义材料初始模量!砂土地弹性模量!!!建立几何模型!!!网格划分!!!边界条件!底边界竖向约束!左侧边界水平向约束!施加围压!!!分级施加荷载，实现非线性计算!!!荷载增量,共级*!!!取出计算结果，修改弹性模量*()!取各单元第一主应力!取各单元第三主应力***为单元编号*()!将单元结果存入数组*()*!AA修改砂土单元地，单元号OOOOOO$$$$$$!各单元循环计算!调用宏文件*()!保存第单元之*!施加荷载，增量!对级荷载情况求解