好像是应力集中，算不过去了，有什么好办法？

hunanyangdong

我做的一个边坡，左上角浅色区域为即将开挖掉的，在开挖形成的坡面上用子程序加分布荷载，可是算到半途算不过去了，其实力并不大，比自重都小很多。附inp文件
这是加力部分的inp
** STEP: Step-3
**
*Step, name=Step-3, nlgeom=YES, inc=1000, unsymm=YES
*Static
0.001, 1., 1e-05, 1.
*DSLOAD
Surf-1, TRVECNU
**
** PREDEFINED FIELDS
**
** Name: Predefined Field-1   Type: Temperature


这是力的子程序
      SUBROUTINE UTRACLOAD(ALPHA,T_USER,KSTEP,KINC,TIME,NOEL,NPT,
     1 COORDS,DIRCOS,JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION T_USER(3), TIME(2), COORDS(3), DIRCOS(3,3)
      CHARACTER*80 SNAME
      y=COORDS(2)
      ALPHA=14000.0*(40-y)*TIME(1)
      T_USER(1)=0.996194698
      T_USER(2)=-0.087155743
      T_USER(3)=0.0
      RETURN
      END