Civil Engineering Reference
In-Depth Information
FUNCTION
SigmaK(dxr,r,E,ny,Cdim)
!--------------------------------------------
! FUNDAMENTAL SOLUTION FOR Normal Stresses
! isotropic material (Kelvin solution)
!--------------------------------------------
REAL,INTENT(IN) ::
dxr(:)
! rx/r etc.
REAL,INTENT(IN) ::
r
! r
REAL,INTENT(IN) ::
E
! Young's modulus
REAL,INTENT(IN) ::
ny
! Poisson's ratio
INTEGER,INTENT(IN) ::
Cdim
! Cartesian dimension
REAL
::
SigmaK(Cdim,Cdim) ! Returns array CdimxCdim
INTEGER
::
n,i,j
REAL
::
G,c,c2,c3,c4
! Temps
G= E/(2.0*(1+ny))
SELECT CASE
(Cdim)
CASE
(2) ! Plane strain solution
n= 1
c2= 1.0/(4.0*Pi*(1.0-ny))
c3= 1.0-2.0*ny
c4= 2.0
CASE
(3) ! Three-dimensional solution
n= 2
c2= 1.0/(8.0*Pi*(1.0-ny))
c3= 1.0-2.0*ny
c4= 3.0
CASE DEFAULT
END SELECT
Direction_Pi: &
DO
i=1,Cdim
Direction_Sigma: &
DO
j=1,Cdim
IF
(i == j)
THEN
SigmaK(i,i)= -c2/r**n*(c3*dxr(i)+c4*dxr(i)**3)
ELSE
SigmaK(i,j)= -c2/r**n*(-c3*dxr(i) + c4*dxr(i)**2*dxr(j))
END IF
END DO &
Direction_Sigma
END DO &
Direction_Pi
RETURN
END FUNCTION SigmaK
The discretised form can be written as
E
N
E
N
¦¦
¦¦
ee
ee
³
cu
P
'
T
u
'
U
t
6
P Q
,
H
Q
dV
(13.46)
i
ni
ni
i
0
n
n
en
11
en
11
V
Search WWH ::
Custom Search