UTRACLOAD不收敛问题

sansirosk

一个二维平面应变问题，我想在表面施加一个法向赫兹载荷和一个与之成比例的切向摩擦力，摩擦系数miu，现在是计算不收敛，请大神指教啊，真诚感谢啊！

1.       SUBROUTINE UTRACLOAD(ALPHA,T_USER,KSTEP,KINC,TIME,NOEL,COORDS,NPT
   
2.      1                     DIRCOS,JLTYP,SNAME)
   
3. C
   
4.       INCLUDE 'ABA_PARAM.INC'
   
5. C
   
6.       DIMENSION T_USER(3),TIME(2),COORDS(3),DIRCOS(3,3)
   
7.       CHARACTER*80 SNAME
   
8. C   实际定义载荷(Start)
   
9.       PARAMETER(Pmax=1706,Con_radius=0.0792,ZERO=0.0d0,Miu=0.1d0)
   
10.           x=COORDS(1)
    
11.           P=Pmax*SQRT(abs(1.0-(x/Con_radius)**2))
    
12.           XLEFT=-Con_radius+0.000002     !设置误差范围，积分点坐标和节点坐标不一样
    
13.           XRIGHT=Con_radius-0.000002     !设置误差范围，积分点坐标和节点坐标不一样
    
14.       IF(x.LE.XRIGHT.AND.x.GE.XLEFT)THEN
    
15.                   ALPHA=P*SQRT(1.0+Miu**2)
    
16.                   T_USER(1)=-1.0/SQRT(1.0+Miu**2)
    
17.           T_USER(2)=-Miu/SQRT(1.0+Miu**2)
    
18.           T_USER(3)=0.0
    
19.       ELSE
    
20.           ALPHA=ZERO
    
21.           T_USER(1)=0.0
    
22.           T_USER(2)=-1.0
    
23.           T_USER(3)=0.0
    
24.       END IF
    
25. C   实际定义载荷(End)
    
26.       RETURN
    
27.       END

复制代码