ansys确定边坡安全系数

/FILNAME,Slope Stability,0
/PREP7
!进行D-P模型参数转化
*AFUN,DEG
fai=5
c=1.2E4
f=1.7
a=2*sin(fai)/(sqrt(3)*(3-sin(fai)))
k=6*c*cos(fai)/(sqrt(3)*(3-sin(fai)))

sfai=sqrt(a**(2)*3.1413*9/(2*3**(1/2)+a**(2)*3.1415))
ss=2**(1/2)*3**(1/4)*sfai/sqrt(3.1415*(9-sfai**(2)))
*if,ss-a,le,1e-4,then
sfai=sfai
*else
sfai=-sfai
*endif
dfai=asin(sfai)
dc=k*sqrt(2*3.1415*(9-sfai**(2)))/6/3**(1/4)/cos(dfai)

dfai=atan(tan(dfai)/f) !转换后的内摩擦角
dc=dc/f !转换后的粘聚力

ET,1,PLANE42 !定义1号单元
KEYOPT,1,3,2 !设置单元关键字为平面应变

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2E7
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,1827
TB,DP,1,,,
TBMODIF,1,1,dc
TBMODIF,1,2,dfai
TBMODIF,1,3,0
!定义关键点
K,1,0,0,0
K,2,15,0,0
K,3,15,8,0
K,4,10,8,0
K,5,5,4,0
K,6,0,4,0
A,1,2,3,4,5,6

ESIZE,0.25 !设置总体单元尺寸
AMESH,ALL !划分所有的面

/SOL
ANTYPE,0
NLGEOM,1 !考虑大变形
NROPT,FULL, ,OFF
!定义位移约束
DL,6, ,UX,
DL,2, ,UX,
DL,1, ,ALL,
ACEL,0,9.8,0, !施加重力加速度
OUTRES,ALL,ALL, !输出所有子步
TIME,1
NSUBST,100 !设置求解子步数为100

SOLVE
/POST1
*GET,U110,NODE,110,U,SUM
*GET,U100,NODE,100,U,SUM
*GET,U92,NODE,92,U,SUM