UTRACLOAD不收敛问题
一个二维平面应变问题，我想在表面施加一个法向赫兹载荷和一个与之成比例的切向摩擦力，摩擦系数miu，现在是计算不收敛，请大神指教啊，真诚感谢啊！ SUBROUTINE UTRACLOAD(ALPHA,T_USER,KSTEP,KINC,TIME,NOEL,COORDS,NPT1 DIRCOS,JLTYP,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION T_USER(3),TIME(2),COORDS(3),DIRCOS(3,3)
CHARACTER*80 SNAME
C 实际定义载荷(Start)
PARAMETER(Pmax=1706,Con_radius=0.0792,ZERO=0.0d0,Miu=0.1d0)
x=COORDS(1)
P=Pmax*SQRT(abs(1.0-(x/Con_radius)**2))
XLEFT=-Con_radius+0.000002 !设置误差范围，积分点坐标和节点坐标不一样
XRIGHT=Con_radius-0.000002 !设置误差范围，积分点坐标和节点坐标不一样
IF(x.LE.XRIGHT.AND.x.GE.XLEFT)THEN
ALPHA=P*SQRT(1.0+Miu**2)
T_USER(1)=-1.0/SQRT(1.0+Miu**2)
T_USER(2)=-Miu/SQRT(1.0+Miu**2)
T_USER(3)=0.0
ELSE
ALPHA=ZERO
T_USER(1)=0.0
T_USER(2)=-1.0
T_USER(3)=0.0
END IF
C 实际定义载荷(End)
RETURN
END
页:
[1]