MAKE_STRESS_TENSOR.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_stress_tensor (nn, lambda, mu, duxdx, duydx, duzdx, duxdy, duydy, duzdy, duxdz, duydz, duzdz, sxx, syy, szz, syz, szx, sxy)
	Computes the stress tensor.

Function/Subroutine Documentation

make_stress_tensor()

subroutine make_stress_tensor	(	integer*4	nn,
		real*8, dimension(nn,nn,nn)	lambda,
		real*8, dimension(nn,nn,nn)	mu,
		real*8, dimension(nn,nn,nn)	duxdx,
		real*8, dimension(nn,nn,nn)	duydx,
			duzdx,
		real*8, dimension(nn,nn,nn)	duxdy,
		real*8, dimension(nn,nn,nn)	duydy,
			duzdy,
		real*8, dimension(nn,nn,nn)	duxdz,
		real*8, dimension(nn,nn,nn)	duydz,
			duzdz,
		real*8, dimension(nn,nn,nn)	sxx,
		real*8, dimension(nn,nn,nn)	syy,
		real*8, dimension(nn,nn,nn)	szz,
		real*8, dimension(nn,nn,nn)	syz,
		real*8, dimension(nn,nn,nn)	szx,
		real*8, dimension(nn,nn,nn)	sxy
	)

Computes the stress tensor.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nn	number of 1-D Legendre nodes
[in]	lambda	nodal values of Lame coefficient lambda
[in]	mu	nodal values of Lame coefficient mu
[in]	duxdx	nodal values for spatial derivatives of the displacement
[in]	duydx	nodal values for spatial derivatives of the displacement
[in]	duzdx	nodal values for spatial derivatives of the displacement
[in]	duxdy	nodal values for spatial derivatives of the displacement
[in]	duydy	nodal values for spatial derivatives of the displacement
[in]	duzdy	nodal values for spatial derivatives of the displacement
[in]	duxdz	nodal values for spatial derivatives of the displacement
[in]	duydz	nodal values for spatial derivatives of the displacement
[in]	duzdz	nodal values for spatial derivatives of the displacement
[out]	sxx	nodal values for the stress tensor
[out]	syy	nodal values for the stress tensor
[out]	szz	nodal values for the stress tensor
[out]	syz	nodal values for the stress tensor
[out]	szx	nodal values for the stress tensor
[out]	sxy	nodal values for the stress tensor

Definition at line 43 of file MAKE_STRESS_TENSOR.f90.

      
      
      implicit none
      
      integer*4 :: nn
      integer*4 :: p,q,r
      
      real*8, dimension(nn,nn,nn) :: lambda,mu
      real*8, dimension(nn,nn,nn) :: sxx,syy,szz,syz,szx,sxy
      real*8, dimension(nn,nn,nn) :: duxdx,duxdy,duxdz,duydx,duydy,duydz,duzdx,duzdy,duzdz
      
 
!     STRESS CALCULATION
      
      do r = 1,nn
         do q = 1,nn
            do p = 1,nn
               
               sxx(p,q,r) = lambda(p,q,r) * (duxdx(p,q,r) +duydy(p,q,r) +duzdz(p,q,r)) &
                                        +2.0d0*mu(p,q,r) * duxdx(p,q,r)
               syy(p,q,r) = lambda(p,q,r) * (duxdx(p,q,r) +duydy(p,q,r) +duzdz(p,q,r)) &
                                        +2.0d0*mu(p,q,r) * duydy(p,q,r)
               szz(p,q,r) = lambda(p,q,r) * (duxdx(p,q,r) +duydy(p,q,r) +duzdz(p,q,r)) &
                                        +2.0d0*mu(p,q,r) * duzdz(p,q,r)
               
               syz(p,q,r) = mu(p,q,r) * (duydz(p,q,r) + duzdy(p,q,r))
               szx(p,q,r) = mu(p,q,r) * (duzdx(p,q,r) + duxdz(p,q,r))
               sxy(p,q,r) = mu(p,q,r) * (duxdy(p,q,r) + duydx(p,q,r))
               
            enddo
         enddo
      enddo
      
      return
      

Referenced by compute_energy_error().

Here is the caller graph for this function: