基于APDL语言实现邓肯－张模型

fini
/clear
/title,nummerical simulation of foundation settlement
/prep7
*dim,Str_max,array,1020  
!偏应力(P1-P3)的历史最大值，1020为单元数，P1、P3为第一、第三主应力
!土力学中，压力为正，拉力为负，P1为竖向(压）应力，P2、P3为侧压应力
*dim,S_max,Array,1020   
!应力水平历史最大值,1020为单元数
*dim,Settlement,array,12
!记录某点沉降值的12×1矩阵，12代表12级加载。（包括下面几个数组，没多大意义）
*dim,M_et,array,12
!最后一个单元，每次加载后的弹性模量
*dim,M_u,array,12
!最后一个单元，每次加载后的泊松比
!此处定义的数组是用于后面的备用计算