pmat3.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    pmat3   ../engine/source/elements/beam/pmat3.F
!||--- called by ------------------------------------------------------
!||    pke3    ../engine/source/elements/beam/pke3.F
!||====================================================================
      SUBROUTINE pmat3(JFT    ,JLT    ,PM     ,GEO     ,YM      ,
     2                 G      ,A1     ,B1     ,B2      ,B3      ,
     3                 AL     ,ALI    ,SH1    ,SH2     ,MTN     ,
     4                 MAT    ,PID    )
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-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT  ,MTN  
      INTEGER MAT(*), PID(*)
C     REAL
      my_real
     .   geo(npropg,*), pm(npropm,*), 
     .    ym(*),g(*),al(*),ali(*),sh1(*),sh2(*),
     .    a1(*),b1(*),b2(*),b3(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX,IPID,J
C     REAL
      my_real
     .   shf(mvsiz),sh,yma2,ali2,sh10,sh0,sh20
C------------- ---------------------------------
C
      DO i=jft,jlt
       g(i)   =pm(22,mat(i))
       ym(i)  =pm(20,mat(i))
       a1(i)  =geo(1,pid(i))
       b1(i)  =geo(2,pid(i))
       b2(i)  =geo(18,pid(i))
       b3(i)  =geo(4,pid(i))
       shf(i) =geo(37,pid(i))
      ENDDO
C
      DO i=jft,jlt
       sh=five_over_6*g(i)*a1(i)
       ali2 =ali(i)*ali(i) 
       yma2=twelve*ym(i)*ali2
       sh10=yma2*b1(i)
       sh20=yma2*b2(i)
       sh0=(one-shf(i))*sh
       sh1(i)=sh0*sh10/(sh+sh10) + shf(i)*sh10
       sh2(i)=sh0*sh20/(sh+sh20) + shf(i)*sh20
      ENDDO
C
      RETURN
      END