dmasanif.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasanif (x, d, elbuf_tab, geo, iparg, ixt, ixp, ixr, mas, pm, el2fa, nbf)

Function/Subroutine Documentation

dmasanif()

subroutine dmasanif	(		x,
			d,
		type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			geo,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixt,*)	ixt,
		integer, dimension(nixp,*)	ixp,
		integer, dimension(nixr,*)	ixr,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf )

Definition at line 29 of file dmasanif.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod
      use element_mod , only : nixt,nixp,nixr
C-----------------------------------------------                        
C   I m p l i c i t   T y p e s                                         
C-----------------------------------------------                        
#include      "implicit_f.inc"                                          
C-----------------------------------------------                        
C   C o m m o n   B l o c k s                                           
C-----------------------------------------------                        
#include      "com01_c.inc"                                             
#include      "com04_c.inc"                                             
#include      "param_c.inc"                                             
C-----------------------------------------------                        
C   D u m m y   A r g u m e n t s                                       
C-----------------------------------------------                        
C     REAL                                                    
      my_real                                                           
     .   mas(*) ,pm(npropm,*),geo(npropg,*),x(3,*),            
     .   d(3,*)                                                         
      INTEGER IPARG(NPARG,*),                                           
     .        IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),EL2FA(*),NBF          
C
      TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
C-----------------------------------------------                        
C   L o c a l   V a r i a b l e s                                       
C-----------------------------------------------                        
C     REAL                                                    
      my_real                                                           
     .   off, p, vonm2, vonm, s1, s2, s12, s3, VALUE,a0,al0,       
     .   rho0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3                       
      INTEGER I, NG, NEL, NFT, ITY, LFT, NPT, ISS, ISC,            
     .        IADD, N, J, LLT, MLW,        
     .        ISTRAIN,NN, K1, K2,JTURB,MT,    
     .        N1,N2,N3,N4,NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,              
     .        OFFSET,NEL_OLD,ITY_OLD,NFT_FA,N_FA,NUVAR                                                     
      REAL R4                                                           
C
      TYPE(G_BUFEL_)  ,POINTER :: GBUF
C-----------------------------------------------
      nn1 = 1                                                           
      nn2 = 1                                                           
      nn3 = 1                                                           
      nn4 = nn3                                                         
      nn5 = nn4                                                         
      nn6 = nn5                                                         
      nn7 = nn6 + numelt                                                
      nn8 = nn7 + numelp                                                
C-----------------------------------------------                        
C                                                                       
      nel_old = 0                                                       
      ity_old = 0                                                       
      DO 490 ng=1,ngroup                                                
        mlw   =iparg(1,ng)                                              
        nel   =iparg(2,ng)                                              
        ity   =iparg(5,ng)                                              
        nft   =iparg(3,ng)                                              
        lft=1                                                           
        llt=nel                                                         
        nft_fa = nft                                                    
C
        gbuf => elbuf_tab(ng)%GBUF
C-----------------------------------------------                        
C       TRUSS                                                           
C-----------------------------------------------                        
        IF(ity==4)THEN                                                
         DO i=lft,llt                                                   
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            rho0 = pm(1,ixt(1,n))                                       
            a0   = geo(1,ixt(4,n))                                      
            n1 = ixt(2,n)                                               
            n2 = ixt(3,n)                                               
            xx1 = x(1,n2)-d(1,n2)-x(1,n1)+d(1,n1)                       
            yy1 = x(2,n2)-d(2,n2)-x(2,n1)+d(2,n1)                       
            zz1 = x(3,n2)-d(3,n2)-x(3,n1)+d(3,n1)                       
            al0  = half*sqrt(xx1*xx1 + yy1*yy1 + zz1*zz1)             
            mas(el2fa(nn6+n_fa)) = rho0*al0*a0                          
         ENDDO                                                          
C-----------------------------------------------                        
C       POUTRES                                                         
C-----------------------------------------------                        
        ELSEIF(ity==5)THEN                                            
         DO i=lft,llt                                                   
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            rho0 = pm(1,ixp(1,n))                                       
            a0   = geo(1,ixp(5,n))                                      
            n1 = ixp(2,n)                                               
            n2 = ixp(3,n)                                               
            xx1 = x(1,n2)-d(1,n2)-x(1,n1)+d(1,n1)                       
            yy1 = x(2,n2)-d(2,n2)-x(2,n1)+d(2,n1)                       
            zz1 = x(3,n2)-d(3,n2)-x(3,n1)+d(3,n1)                       
            al0  = half*sqrt(xx1*xx1 + yy1*yy1 + zz1*zz1)             
            mas(el2fa(nn7+n_fa)) = rho0*al0*a0                          
         ENDDO                                                          
C-----------------------------------------------                        
C       RESSORTS                                                        
C-----------------------------------------------                        
        ELSEIF(ity==6)THEN                                            
          IF(mlw==3)THEN                                              
           DO i=lft,llt                                                 
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
              mas(el2fa(nn8+n_fa))   = half*geo(1,ixr(1,n))           
              mas(el2fa(nn8+n_fa)+1) = half*geo(1,ixr(1,n))           
           ENDDO                                                        
          ELSEIF(mlw==5) THEN                                         
cc           NB1 =IAD - 1                                                 
cc           NB2 =NB1  + NEL                                              
           nuvar =  nint(geo(25,ixr(1,1+nft)))                          
cc           NB3 =NB2  + 3*NEL                                            
cc           NB4 =NB3  + NEL                                              
cc           NB5 =NB4  + 3*NEL                                            
cc           NB6 =NB5                                                     
cc           NB7 =NB6                                                     
cc           NB8 =NB7                                                     
cc           NB9 =NB8  + 3*NEL                                            
cc           NB10=NB9  + 3*NEL                                            
cc           NB11=NB10                                                    
cc           NB12=NB11                                                    
cc           NB13=NB12                                                    
cc           NB14=NB13                                                    
cc           NB15 = NB14 + 3*NEL                                          
cc           NB16 = NB15 + NUVAR*NEL                                      
           DO i=lft,llt                                                 
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            mas(el2fa(nn8+n_fa)) = gbuf%MASS(i)
cc            MAS(EL2FA(NN8+N_FA)) = BUFEL(NB16+I)                        
           ENDDO                                                        
          ELSE                                                          
           DO i=lft,llt                                                 
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            mas(el2fa(nn8+n_fa)) = geo(1,ixr(1,n))                      
           ENDDO                                                        
          ENDIF                                                         
        ELSE                                                            
        ENDIF                                                           
C-----------------------------------------------                        
C       END OF LOOP ON OFFSETS
C-----------------------------------------------                        
 490  CONTINUE                                                          
C-----------------------------------------------                        
C                                                                       
      RETURN