fric子程序求助
浏览:80233 回答:1
新手求助,有两个问题,第一个是OMEGA=78.54D0*MAX(1.0D0-HTIME/4.5D0,0.0D0)中,78.54DO是什么意思?第二个问题TAU(1)括号中的1是摩擦方向的数值,那1代表什么方向呢?谢谢各位了,帮帮忙。
C
IF (LM .EQ. 2) RETURN
C
C SPEED OF ROTATION, RADIUS AND CIRCUMFERENTIAL SLIP INCREMENT
C
HTIME=TIME(2)-0.5*DTIME
OMEGA=78.54D0*MAX(1.0D0-HTIME/4.5D0,0.0D0)
R=COORDS(1)
DSLIP(1)=OMEGA*DTIME*R
C
C COEFFICIENT OF FRICTION
C
AMU=PROPS(1)
C
C PLASTIC DISSIPATION
C
SPD=AMU*PRESS*DSLIP(1)
LM=0
TAU(1)=0.
DDTDDG(1,1)=0.
DDTDDP(1)=0.
DDTDDT(1,1)=0.
C
RETURN
END
C
IF (LM .EQ. 2) RETURN
C
C SPEED OF ROTATION, RADIUS AND CIRCUMFERENTIAL SLIP INCREMENT
C
HTIME=TIME(2)-0.5*DTIME
OMEGA=78.54D0*MAX(1.0D0-HTIME/4.5D0,0.0D0)
R=COORDS(1)
DSLIP(1)=OMEGA*DTIME*R
C
C COEFFICIENT OF FRICTION
C
AMU=PROPS(1)
C
C PLASTIC DISSIPATION
C
SPD=AMU*PRESS*DSLIP(1)
LM=0
TAU(1)=0.
DDTDDG(1,1)=0.
DDTDDP(1)=0.
DDTDDT(1,1)=0.
C
RETURN
END
do 是终止好像 1 应该是x方向吧 不是很确定