forked from sanjayg0/feappv
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathumatl1.f
50 lines (38 loc) · 1.75 KB
/
umatl1.f
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
49
50
!$Id:$
subroutine umatl1(eps,theta,td,d,ud,hn,h1,nh,ii,istrt, sig,dd,isw)
! * * F E A P * * A Finite Element Analysis Program
!.... Copyright (c) 1984-2024: Regents of the University of California
! All rights reserved
!-----[--.----+----.----+----.-----------------------------------------]
! Purpose: User Constitutive Model 1
! Input:
! eps(*) - Current strains at point (small deformation)
! - Deformation gradient at point (finite deformation)
! theta - Trace of strain at point
! - Determinant of deforamtion gradient
! td - Temperature change
! d(*) - Program material parameters (ndd)
! ud(*) - User material parameters (nud)
! hn(nh) - History terms at point: t_n
! h1(nh) - History terms at point: t_n+1
! nh - Number of history terms
! ii - Current point number
! istrt - Start state: 0 = elastic; 1 = last solution
! isw - Solution option from element
! Output:
! sig(*) - Stresses at point.
! N.B. 1-d models use only sig(1)
! dd(6,*) - Current material tangent moduli
! N.B. 1-d models use only dd(1,1) and dd(2,1)
!-----[--.----+----.----+----.-----------------------------------------]
implicit none
integer :: nh,istrt,isw, ii
real (kind=8) :: td
real (kind=8) :: eps(*),theta(*),d(*),ud(*),hn(nh),h1(nh)
real (kind=8) :: sig(*),dd(6,*)
! Set initial values to history values in both hn(*) and h1(*)
if(isw.eq.14) then
! Compute and output stress (sig) and (moduli)
else
endif
end subroutine umatl1