cdkderi3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cdkderic3 (jft, jlt, x2, y2, x3, y3, area2, alpe, aldt, px2, py2, px3, py3, px, py, pxy, pyy, vol, volg, nu, thk2)
subroutine	cdkderi3 (jft, jlt, px2, py2, px3, py3, px, py, pxy, pyy, bz1, bz2, bz3, brx1, brx2, brx3, bry1, bry2, bry3, ksi, eta)

Function/Subroutine Documentation

cdkderi3()

subroutine cdkderi3	(	integer	jft,
		integer	jlt,
			px2,
			py2,
			px3,
			py3,
			px,
			py,
			pxy,
			pyy,
			bz1,
			bz2,
			bz3,
			brx1,
			brx2,
			brx3,
			bry1,
			bry2,
			bry3,
			ksi,
			eta )

Definition at line 104 of file cdkderi3.F.

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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT
      my_real
     .     px2(*), py2(*), px3(*), py3(*), 
     .     px(mvsiz,3), py(mvsiz,3), pxy(mvsiz,3), pyy(mvsiz,3),
     .     bz1(mvsiz,2),bz2(mvsiz,2),bz3(mvsiz,2), brx1(mvsiz,3),
     .     brx2(mvsiz,3),brx3(mvsiz,3),bry1(mvsiz,3),bry2(mvsiz,3),bry3(mvsiz,3),
     .     ksi, eta
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real :: xr1,xr2,xr3,yr1,yr2,yr3,xs1,xs2,xs3,ys1,ys2,ys3,x2y,y2x,x6,y6,xy6,rxy1(3),rxy2(3)
C-----------------------------------------------
      x2y =one-two*ksi-eta
      y2x =one-two*eta-ksi
      x6  = six*ksi-one
      y6  = six*eta-one
      xy6 = x6+y6-three      
      DO i=jft,jlt
       xr1= px(i,3)*x2y
       xr2= px(i,2)*eta
       xr3= px(i,1)*eta
       xs1= -px(i,2)*y2x
       xs2= -px(i,3)*ksi
       xs3= -px(i,1)*ksi
       bz1(i,1)=px2(i)*(-xr1+xr3)-px3(i)*(xs3+xs2)
       bz1(i,2)=-px2(i)*(xr2+xr3)+px3(i)*(-xs1+xs3)
C       BZ1(1,I)=-BZ1(2,I)-BZ1(3,I)
C
       yr1= py(i,3)*x2y
       yr2= py(i,2)*eta
       yr3= py(i,1)*eta
       ys1= -py(i,2)*y2x
       ys2= -py(i,3)*ksi
       ys3= -py(i,1)*ksi
       bz2(i,1)=py2(i)*(-yr1+yr3)-py3(i)*(ys3+ys2)
       bz2(i,2)=-py2(i)*(yr2+yr3)+py3(i)*(-ys1+ys3)
C       BZ2(1,I)=-BZ2(2,I)-BZ2(3,I)
C
       bz3(i,1)=px2(i)*(-yr1+yr3)-px3(i)*(ys3+ys2)+
     .           py2(i)*(-xr1+xr3)-py3(i)*(xs3+xs2)
       bz3(i,2)=-px2(i)*(yr2+yr3)+px3(i)*(-ys1+ys3)
     .           -py2(i)*(xr2+xr3)+py3(i)*(-xs1+xs3)
C       BZ3(1,I)=-BZ3(2,I)-BZ3(3,I)
      ENDDO
C
      DO i=jft,jlt
       xr1= pxy(i,3)*x2y
       xr2= -pxy(i,2)*eta
       xr3= pxy(i,1)*eta
       xs1= pxy(i,2)*y2x
       xs2= -pxy(i,3)*ksi
       xs3= pxy(i,1)*ksi
       brx1(i,1)=px2(i)*(xr1+xr2)+px3(i)*(xs1+xs2)
       brx1(i,2)=px2(i)*(xr1+xr3)+px3(i)*(xs3+xs2)
       brx1(i,3)=px2(i)*(xr2+xr3)+px3(i)*(xs1+xs3)
C
       yr1= pyy(i,3)*x2y
       yr2= one-pyy(i,2)*eta
       yr3= pyy(i,1)*eta-one
       ys1= pyy(i,2)*y2x
       ys2= one-pyy(i,3)*ksi
       ys3= pyy(i,1)*ksi-one
       brx2(i,1)=py2(i)*(yr1+yr2)+py3(i)*(ys1+ys2)
       brx2(i,2)=py2(i)*(yr1+yr3)+py3(i)*(ys3+ys2)
       brx2(i,3)=py2(i)*(yr2+yr3)+py3(i)*(ys1+ys3)
C
       rxy1(1)=px2(i)*(yr1+yr2)+px3(i)*(ys1+ys2)
       rxy2(1)=py2(i)*(xr1+xr2)+py3(i)*(xs1+xs2)
       rxy1(2)=px2(i)*(yr1+yr3)+px3(i)*(ys3+ys2)
       rxy2(2)=py2(i)*(xr1+xr3)+py3(i)*(xs3+xs2)
       rxy1(3)=px2(i)*(yr2+yr3)+px3(i)*(ys1+ys3)
       rxy2(3)=py2(i)*(xr2+xr3)+py3(i)*(xs1+xs3)
       brx3(i,1)=rxy1(1)+rxy2(1)
       brx3(i,2)=rxy1(2)+rxy2(2)
       brx3(i,3)=rxy1(3)+rxy2(3)
C
       bry1(i,1)=rxy1(1)+(px2(i)+px3(i))*xy6
       bry1(i,2)=rxy1(2)+px2(i)*x6
       bry1(i,3)=rxy1(3)+px3(i)*y6
C
       bry2(i,1)=-rxy2(1)
       bry2(i,2)=-rxy2(2)
       bry2(i,3)=-rxy2(3)
C
       bry3(i,1)=(py2(i)+py3(i))*xy6+brx2(i,1)-brx1(i,1)
       bry3(i,2)=py2(i)*x6+brx2(i,2)-brx1(i,2)
       bry3(i,3)=py3(i)*y6+brx2(i,3)-brx1(i,3)
      ENDDO
C
C---------------------------------------------------------
      RETURN
C

cdkderic3()

subroutine cdkderic3	(	integer	jft,
		integer	jlt,
			x2,
			y2,
			x3,
			y3,
			area2,
			alpe,
			aldt,
			px2,
			py2,
			px3,
			py3,
			px,
			py,
			pxy,
			pyy,
			vol,
			volg,
			nu,
			thk2 )

Definition at line 28 of file cdkderi3.F.

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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT
      my_real
     .     x2(*), y2(*), x3(*), y3(*),area2(*), 
     .     px2(*), py2(*), px3(*), py3(*), 
     .     px(mvsiz,3), py(mvsiz,3), pxy(mvsiz,3), pyy(mvsiz,3), 
     .     aldt(*),alpe(*),vol(*),volg(*),nu(*),thk2(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
      my_real :: areai(mvsiz),al1, al2, al3, almax,almin,fac
      my_real :: x32(mvsiz), y32(mvsiz),al4, al5, al6
C-----------------------------------------------
      DO i=jft,jlt
       areai(i)=one/area2(i)
       px2(i)=y3(i)*areai(i)
       py2(i)=-x3(i)*areai(i)
       px3(i)=-y2(i)*areai(i)
       py3(i)=x2(i)*areai(i)
       x32(i) =x3(i)-x2(i)
       y32(i) =y3(i)-y2(i)
      ENDDO
C
      DO i=jft,jlt
       al1 = x32(i)*x32(i) + y32(i)*y32(i)       
       al2 = x3(i) * x3(i) + y3(i) * y3(i)
       al3 = x2(i) * x2(i) + y2(i) * y2(i)
       almax = max(al1,al2,al3)
       almin = min(al1,al2,al3)
       fac =1.0+0.6*(1+nu(i))*thk2(i)/almin
       almax = almax*fac
       aldt(i)= area2(i) / sqrt(almax)
       alpe(i)=one
       al4 = three/max(al1,em20)
       al5 = three/max(al2,em20)
       al6 = three/max(al3,em20)
       pxy(i,1)=x32(i)*y32(i)*al4
       pxy(i,2)=x3(i)*y3(i)*al5
       pxy(i,3)=x2(i)*y2(i)*al6
       pyy(i,1)=y32(i)*y32(i)*al4
       pyy(i,2)=y3(i)*y3(i)*al5
       pyy(i,3)=y2(i)*y2(i)*al6
       al4 = two*al4
       al5 = two*al5
       al6 = two*al6       
       px(i,1)=x32(i)*al4
       px(i,2)=-x3(i)*al5
       px(i,3)=x2(i)*al6
       py(i,1)=y32(i)*al4
       py(i,2)=-y3(i)*al5
       py(i,3)=y2(i)*al6
       vol(i) = third*volg(i)
      ENDDO
C---------------------------------------------------------
      RETURN
C