psumg3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	psumg3 (jft, jlt, ym, g, area, b1, b2, b3, al, ali, sh1, sh2, r11, r12, r13, r21, r22, r23, r31, r32, r33, ke11, ke12, ke22)

Function/Subroutine Documentation

psumg3()

subroutine psumg3	(	integer	jft,
		integer	jlt,
			ym,
			g,
			area,
			b1,
			b2,
			b3,
			al,
			ali,
			sh1,
			sh2,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
			ke11,
			ke12,
			ke22 )

Definition at line 28 of file psumg3.F.

C-----------------------------------------------
#include      "implicit_f.inc"
#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 ym(*),g(*),al(*),ali(*),sh1(*),sh2(*),area(*),b1(*),b2(*),b3(*),
     .        r11(*),r12(*),r13(*),r21(*),r22(*),r23(*),r31(*),r32(*),r33(*)
       my_real ke11(6,6,*),ke22(6,6,*),ke12(6,6,*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I,J,EP,MI,MJ
      my_real q(3,3,mvsiz),k11(3,mvsiz),m11(3,mvsiz),mf32(mvsiz),m12(3,mvsiz),mf23(mvsiz),al2,eli,q1,q2,q3,mf11(3,3,mvsiz)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
C--- B1 : Iyy ;B2 : Izz ;B3 : Ixx ;
       DO i=jft,jlt
         al2= al(i)*al(i)
         eli = ym(i)*ali(i)
         k11(1,i)=eli*area(i)
         k11(2,i)=sh2(i)*ali(i)
         k11(3,i)=sh1(i)*ali(i)
         m11(1,i)=g(i)*b3(i)*ali(i)
         m11(2,i)=fourth*al(i)*sh1(i)+eli*b1(i)
         m11(3,i)=fourth*al(i)*sh2(i)+eli*b2(i)
         mf23(i)=half*sh2(i)
         mf32(i)=half*sh1(i)
         m12(1,i)=-m11(1,i)
         m12(2,i)=al(i)*mf32(i)-m11(2,i)
         m12(3,i)=al(i)*mf23(i)-m11(3,i)
         mf32(i)=-mf32(i)
       ENDDO
C---------------------------------------
C   TRANS LOCAL-->GLOBAL 
C---------------------------------------
       DO i=jft,jlt
         q(1,1,i)=r11(i)
         q(1,2,i)=r21(i)
         q(1,3,i)=r31(i)
         q(2,1,i)=r12(i)
         q(2,2,i)=r22(i)
         q(2,3,i)=r32(i)
         q(3,1,i)=r13(i)
         q(3,2,i)=r23(i)
         q(3,3,i)=r33(i)
       ENDDO 
C---------------------------------------
C   ASSEMBLAGE
C---------------------------------------
C---------KII ----Keij=QkiQkj*Kkk---M12 est diag- 
       DO i=1,3 
         mi=i+3
         DO j=i,3
           mj=j+3
           DO ep=jft,jlt
             q1 =q(1,i,ep)*q(1,j,ep)
             q2 =q(2,i,ep)*q(2,j,ep)
             q3 =q(3,i,ep)*q(3,j,ep)
             ke11(i,j,ep)=q1*k11(1,ep)+q2*k11(2,ep)+q3*k11(3,ep)
             ke11(mi,mj,ep)=q1*m11(1,ep)+q2*m11(2,ep)+q3*m11(3,ep)
             ke12(mi,mj,ep)=q1*m12(1,ep)+q2*m12(2,ep)+q3*m12(3,ep)
             ke22(i,j,ep)=ke11(i,j,ep)
             ke22(mi,mj,ep)=ke11(mi,mj,ep)
           ENDDO
         ENDDO
       ENDDO
C---------K23>0,K32>0-------------------------
       DO i=1,3 
         DO j=1,3
           DO ep=jft,jlt
             mf11(i,j,ep)=q(2,i,ep)*mf23(ep)*q(3,j,ep)+q(3,i,ep)*mf32(ep)*q(2,j,ep)
           ENDDO
         ENDDO
       ENDDO
 
       DO i=1,3 
         DO j=1,3
           mj=j+3
           DO ep=jft,jlt
             ke11(i,mj,ep)=mf11(i,j,ep)
             ke22(i,mj,ep)=-mf11(i,j,ep)
           ENDDO
         ENDDO
       ENDDO
 
       DO i=1,6 
         DO j=i,6
           DO ep=jft,jlt
             ke11(j,i,ep)=ke11(i,j,ep)
             ke22(j,i,ep)=ke22(i,j,ep)
           ENDDO
         ENDDO
       ENDDO
       !---------KIJ --------
       DO i=1,3 
         mi=i+3
         DO j=1,3
           mj=j+3
           DO ep=jft,jlt
             ke12(i,j,ep)=-ke11(i,j,ep)
             ke12(i,mj,ep)=mf11(i,j,ep)
             ke12(mi,j,ep)=-mf11(j,i,ep)
           ENDDO
         ENDDO
       ENDDO
       DO ep=jft,jlt
          ke12(5,4,ep)=ke12(4,5,ep)
          ke12(6,4,ep)=ke12(4,6,ep)
          ke12(6,5,ep)=ke12(5,6,ep)
       ENDDO
 
      RETURN