1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
| C
C This is for pure titanium :
C
PROGRAM plasticity
IMPLICIT NONE
REAL sigma
REAL epsilon
REAL alpha
REAL epsilonP
INTEGER N,i
REAL k
REAL epsmax
REAL E
REAL sigmaY
REAL epsilondot
REAL F
REAL Depsilon
REAL Dtime
C
C initialisations :
C
N=500
k=0.4
epsmax=0.13
E=110000
sigmaY=260
epsilondot=10e-2
sigma=0
epsilon=0
alpha=0
epsilonP=0
Depsilon=epsmax/N
Dtime=Depsilon/epsilondot
5 DO i=1,N
epsilon=epsilon+Depsilon
sigma=E*(epsilon-epsilonP)
F=abs(sigma)-sigmaY+k*alpha
IF (F.LT.0) THEN
sigma=sigma
GO TO 5
ELSE
sigma=(1-((F/(E+k))*E)/(abs(sigma)))*sigma
epsilonP=epsilonP+(F/(E+k))*sign(sigma)
alpha=alpha+(F/(E+k))
GO TO 5
END IF
END DO
END |
Partager