s6cdefo3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine

s6cdefc3 (px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vx5, vx6, vy1, vy2, vy3, vy4, vy5, vy6, vz1, vz2, vz3, vz4, vz5, vz6, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz, dhxx, dhxy, dhxz, dhyx, dhyy, dhyz, dhzx, dhzy, dhzz, dh4, dh5, dh6, px1h, px2h, px3h, py1h, py2h, py3h, pz1h, pz2h, pz3h, ji33, b1x, b1y, b2y, b2x, b1122, b1221, b2212, b1121, b1xh, b1yh, b2xh, b2yh, b1122h, b1221h, b2212h, b1121h, ddhv, nu, nel)

Function/Subroutine Documentation

s6cdefc3()

subroutine s6cdefc3	(		px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			d4,
			d5,
			d6,
			wxx,
			wyy,
			wzz,
			dhxx,
			dhxy,
			dhxz,
			dhyx,
			dhyy,
			dhyz,
			dhzx,
			dhzy,
			dhzz,
			dh4,
			dh5,
			dh6,
			px1h,
			px2h,
			px3h,
			py1h,
			py2h,
			py3h,
			pz1h,
			pz2h,
			pz3h,
			ji33,
			b1x,
			b1y,
			b2y,
			b2x,
			b1122,
			b1221,
			b2212,
			b1121,
			b1xh,
			b1yh,
			b2xh,
			b2yh,
			b1122h,
			b1221h,
			b2212h,
			b1121h,
			ddhv,
			nu,
		integer, intent(in)	nel )

Definition at line 28 of file s6cdefo3.F.

C   ---
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
C     REAL
      my_real
     .   vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), 
     .   vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), 
     .   vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), 
     .   px1(*), px2(*), px3(*), px4(*),
     .   py1(*), py2(*), py3(*), py4(*),
     .   pz1(*), pz2(*), pz3(*), pz4(*), 
     .   dxx(*), dxy(*), dxz(*),
     .   dyx(*), dyy(*), dyz(*),
     .   dzx(*), dzy(*), dzz(*), d4(*), d5(*), d6(*),
     .   wxx(*), wyy(*), wzz(*),
C   +++
     .   px1h(*), px2h(*), px3h(*), 
     .   py1h(*), py2h(*), py3h(*), 
     .   pz1h(*), pz2h(*), pz3h(*), 
     .   ji33(*),b1x(mvsiz,2),b1y(mvsiz,2),b2x(mvsiz,2),b2y(mvsiz,2),
     .   b1xh(mvsiz,2),b1yh(mvsiz,2),b2xh(mvsiz,2),b2yh(mvsiz,2),
     .   b1122(*),b1221(*),b2212(*),b1121(*),
     .   b1122h(*),b1221h(*),b2212h(*),b1121h(*),
     .   dhxx(*), dhxy(*), dhxz(*),  
     .   dhyx(*), dhyy(*), dhyz(*), 
     .   dhzx(*), dhzy(*), dhzz(*),dh4(*), dh5(*), dh6(*),
     .   ddhv(*),nu(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I, J
C     REAL
      my_real
     .   vx14(mvsiz), vy14(mvsiz), vz14(mvsiz),
     .   vx25(mvsiz), vy25(mvsiz), vz25(mvsiz),
     .   vx36(mvsiz), vy36(mvsiz), vz36(mvsiz),
     .   vx14n(mvsiz), vy14n(mvsiz), vz14n(mvsiz),
     .   vx25n(mvsiz), vy25n(mvsiz), vz25n(mvsiz),
     .   vx36n(mvsiz), vy36n(mvsiz), vz36n(mvsiz),
     .   vx3614n(mvsiz),vy3614n(mvsiz),vx2514n(mvsiz),vy2514n(mvsiz),
     .   vxhi(mvsiz), vyhi(mvsiz), vzhi(mvsiz), dt1d2, dt1d1
      my_real 
     .  nu1,pxsvx,pysvy,pzsvz,pxavx,pyavy,pzavz,termad
C-----------------------------------------------
      DO i=1,nel
       vx14(i)=vx1(i)+vx4(i)
       vx25(i)=vx2(i)+vx5(i)
       vx36(i)=vx3(i)+vx6(i)
       vxhi(i)=vx4(i)+vx5(i)+vx6(i)-vx1(i)-vx2(i)-vx3(i)
       vy14(i)=vy1(i)+vy4(i)
       vy25(i)=vy2(i)+vy5(i)
       vy36(i)=vy3(i)+vy6(i)
       vyhi(i)=vy4(i)+vy5(i)+vy6(i)-vy1(i)-vy2(i)-vy3(i)
       vz14(i)=vz1(i)+vz4(i)
       vz25(i)=vz2(i)+vz5(i)
       vz36(i)=vz3(i)+vz6(i)
       vzhi(i)=vz4(i)+vz5(i)+vz6(i)-vz1(i)-vz2(i)-vz3(i)
      ENDDO
      DO i=1,nel
       vx14n(i)=-vx1(i)+vx4(i)
       vx25n(i)=-vx2(i)+vx5(i)
       vx36n(i)=-vx3(i)+vx6(i)
       vy14n(i)=-vy1(i)+vy4(i)
       vy25n(i)=-vy2(i)+vy5(i)
       vy36n(i)=-vy3(i)+vy6(i)
       vz14n(i)=-vz1(i)+vz4(i)
       vz25n(i)=-vz2(i)+vz5(i)
       vz36n(i)=-vz3(i)+vz6(i)
C
       vx3614n(i)=vx36n(i)-vx14n(i)
       vy3614n(i)=vy36n(i)-vy14n(i)
       vx2514n(i)=vx25n(i)-vx14n(i)
       vy2514n(i)=vy25n(i)-vy14n(i)
      ENDDO
C                                                                     12
C   ----------constant part--------------------------------------
      DO i=1,nel
        dxx(i)=px1(i)*vx14(i)+px2(i)*vx25(i)+
     .         px3(i)*vx36(i)+px4(i)*vxhi(i)
        dyy(i)=py1(i)*vy14(i)+py2(i)*vy25(i)+
     .         py3(i)*vy36(i)+py4(i)*vyhi(i)
        dzz(i)=pz1(i)*vz14(i)+pz2(i)*vz25(i)+
     .         pz3(i)*vz36(i)+pz4(i)*vzhi(i)
        dxy(i)=py1(i)*vx14(i)+py2(i)*vx25(i)+
     .         py3(i)*vx36(i)+py4(i)*vxhi(i)
        dyx(i)=px1(i)*vy14(i)+px2(i)*vy25(i)+
     .         px3(i)*vy36(i)+px4(i)*vyhi(i)
C   --------shear traitement-------
        dxz(i)=ji33(i)*two*vxhi(i)
     .         -b1122(i)*vx36n(i)+b1221(i)*vx25n(i)
     .         +b2212(i)*(vy25n(i)-vy36n(i))  
        dxz(i)= dxz(i)+ji33(i)*(vxhi(i)-vx14n(i))
     .         -b1x(i,1)*vx3614n(i)+b1x(i,2)*vx2514n(i)  
     .         -b1y(i,1)*vy3614n(i)+b1y(i,2)*vy2514n(i)  
        dyz(i)=ji33(i)*two*vyhi(i)
     .         -b1122(i)*vy25n(i)+b1221(i)*vy36n(i)
     .         +b1121(i)*(vx36n(i)-vx25n(i))  
        dyz(i)= dyz(i)+ ji33(i)*(vyhi(i)-vy14n(i))
     .         +b2x(i,1)*vx3614n(i)-b2x(i,2)*vx2514n(i)
     .         +b2y(i,1)*vy3614n(i)-b2y(i,2)*vy2514n(i)  
C
        dzx(i)=half*(px1(i)*vz14(i)+px2(i)*vz25(i)+
     .         px3(i)*vz36(i))
        dzy(i)=half*(py1(i)*vz14(i)+py2(i)*vz25(i)+
     .         py3(i)*vz36(i))
      ENDDO
C
      DO i=1,nel
        wxx(i)=zero
        wyy(i)=zero
        wzz(i)=zero
      ENDDO
C   ----------non-constant part--------------------------------------
      DO i=1,nel
        nu1 = nu(i)/(one - nu(i))
        pxavx = px1(i)*vx14n(i)+px2(i)*vx25n(i)+px3(i)*vx36n(i)
        pyavy = py1(i)*vy14n(i)+py2(i)*vy25n(i)+py3(i)*vy36n(i)
        pzavz = pz1(i)*vz14n(i)+pz2(i)*vz25n(i)+pz3(i)*vz36n(i)
        pxsvx = px1h(i)*vx14(i)+px2h(i)*vx25(i)+px3h(i)*vx36(i)
        pysvy = py1h(i)*vy14(i)+py2h(i)*vy25(i)+py3h(i)*vy36(i)
        pzsvz = pz1h(i)*vz14(i)+pz2h(i)*vz25(i)+pz3h(i)*vz36(i)
        termad=-nu(i)*pyavy-nu1*pzavz
        ddhv(i)= termad
        dhxx(i)= pxsvx+pxavx+termad
        termad=-nu(i)*pxsvx-nu1*pzsvz
        ddhv(i)= ddhv(i)+termad
        dhyy(i)= pysvy+pyavy+termad
        termad=-nu(i)*(pxsvx+pyavy)-nu1*(pxavx+pysvy)
        ddhv(i)= ddhv(i)+termad
        dhzz(i)= pzsvz+pzavz+termad
        dhxy(i)=py1(i)*vx14n(i)+py2(i)*vx25n(i)+
     .          py3(i)*vx36n(i)+
     .          py1h(i)*vx14(i)+py2h(i)*vx25(i)+
     .          py3h(i)*vx36(i)
        dhyx(i)=px1(i)*vy14n(i)+px2(i)*vy25n(i)+
     .          px3(i)*vy36n(i)+
     .          px1h(i)*vy14(i)+px2h(i)*vy25(i)+
     .          px3h(i)*vy36(i)
C   --------shear traitement-------
        dhxz(i)=-b1122h(i)*vx36n(i)+b1221h(i)*vx25n(i)
     .          +b2212h(i)*(vy25n(i)-vy36n(i))  
        dhxz(i)= dhxz(i)
     .         -b1xh(i,1)*vx3614n(i)+b1xh(i,2)*vx2514n(i)  
     .         -b1yh(i,1)*vy3614n(i)+b1yh(i,2)*vy2514n(i)  
        dhyz(i)=-b1122h(i)*vy25n(i)+b1221h(i)*vy36n(i)
     .          +b1121h(i)*(vx36n(i)-vx25n(i))  
        dhyz(i)= dhyz(i) 
     .         +b2xh(i,1)*vx3614n(i)-b2xh(i,2)*vx2514n(i)
     .         +b2yh(i,1)*vy3614n(i)-b2yh(i,2)*vy2514n(i)  
        dhzx(i)=half*(px1h(i)*vz14(i)+px2h(i)*vz25(i)+
     .          px3h(i)*vz36(i))
        dhzy(i)=half*(py1h(i)*vz14(i)+py2h(i)*vz25(i)+
     .          py3h(i)*vz36(i))
      ENDDO
C
      RETURN