surfmas.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	surfmas (ms, ibufn, ity, amasu, x, id, addmas, admid, titr)

Function/Subroutine Documentation

surfmas()

subroutine surfmas	(		ms,
		integer, dimension(*)	ibufn,
		integer	ity,
			amasu,
			x,
		integer	id,
			addmas,
		integer	admid,
		character(len=nchartitle)	titr )

Definition at line 32 of file surfmas.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE message_mod
        USE names_and_titles_mod , ONLY : nchartitle
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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IBUFN(*),ITY,ID,ADMID
      my_real ms(*),amasu,x(3,*),addmas
      CHARACTER(LEN=NCHARTITLE)::TITR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER K
      my_real
     .   x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4,
     .   x21,y21,z21,x31,y31,z31,x42,y42,z42,
     .   x32,y32,z32,e3x,e3y,e3z,smass,sum,area,
     .   a2,b2,c2,aa,bb,cc,ang1,ang2,ang3
C=======================================================================
      IF(ity /= 7)THEN
C----------------------------------------------
C     MASSES ELEMENTS /4
C----------------------------------------------
        x1=x(1,ibufn(1))
        y1=x(2,ibufn(1))
        z1=x(3,ibufn(1))
        x2=x(1,ibufn(2))
        y2=x(2,ibufn(2))
        z2=x(3,ibufn(2))
        x3=x(1,ibufn(3))
        y3=x(2,ibufn(3))
        z3=x(3,ibufn(3))
        x4=x(1,ibufn(4))
        y4=x(2,ibufn(4))
        z4=x(3,ibufn(4))
C
        x21=x2-x1
        y21=y2-y1
        z21=z2-z1
        x31=x3-x1
        y31=y3-y1
        z31=z3-z1
        x42=x4-x2
        y42=y4-y2
        z42=z4-z2
C
        e3x=y31*z42-z31*y42
        e3y=z31*x42-x31*z42
        e3z=x31*y42-y31*x42
        sum=sqrt(e3x*e3x+e3y*e3y+e3z*e3z)
        area=half*sum
C
C        EMS = RHO * THK * AREA
C
        smass = amasu * area  ! THK is already included in the input AMASU ( = MASS/AREA = RHO * THK)
C nodal mass
        ms(ibufn(1))=ms(ibufn(1)) + smass * fourth
        ms(ibufn(2))=ms(ibufn(2)) + smass * fourth
        ms(ibufn(3))=ms(ibufn(3)) + smass * fourth
        ms(ibufn(4))=ms(ibufn(4)) + smass * fourth
C
        addmas = addmas + smass
C
      ELSE IF(ity == 7)THEN
C----------------------------------------------
C     MASSES ELEMENTS * ANGLE
C----------------------------------------------
        x1=x(1,ibufn(1))
        y1=x(2,ibufn(1))
        z1=x(3,ibufn(1))
        x2=x(1,ibufn(2))
        y2=x(2,ibufn(2))
        z2=x(3,ibufn(2))
        x3=x(1,ibufn(3))
        y3=x(2,ibufn(3))
        z3=x(3,ibufn(3))
C
        x21=x2-x1
        y21=y2-y1
        z21=z2-z1
        x31=x3-x1
        y31=y3-y1
        z31=z3-z1
        x32=x3-x2
        y32=y3-y2
        z32=z3-z2
C
        a2 = x21**2 + y21**2 + z21**2
        b2 = x32**2 + y32**2 + z32**2
        c2 = x31**2 + y31**2 + z31**2
        aa = sqrt(a2)
        bb = sqrt(b2)
        cc = sqrt(c2)
C
        ang1 = acos((a2 + c2 - b2)/(two * aa * cc)) / pi
        ang2 = acos((a2 + b2 - c2)/(two * aa * bb)) / pi
        ang3 = acos((b2 + c2 - a2)/(two * bb * cc)) / pi
C
        IF ( ( (a2 + c2 - b2)/(2. * aa * cc) <= -one ) .OR.
     .       ( (a2 + c2 - b2)/(2. * aa * cc) >=  one ) .OR.
     .       ( (a2 + b2 - c2)/(2. * aa * bb) <= -one ) .OR.
     .       ( (a2 + b2 - c2)/(2. * aa * bb) >=  one ) .OR.
     .       ( (b2 + c2 - a2)/(2. * bb * cc) <= -one ) .OR.
     .       ( (b2 + c2 - a2)/(2. * bb * cc) >=  one ) ) THEN      
          CALL ancmsg(msgid=880,
     .                msgtype=msgerror,
     .                anmode=aninfo,
     .                i1=admid,
     .                c1=titr,i2=id)
        ENDIF
C
        e3x=y21*z31-z21*y31
        e3y=z21*x31-x21*z31
        e3z=x21*y31-y21*x31
        sum=sqrt(e3x*e3x+e3y*e3y+e3z*e3z)
        area=half*sum
C
        smass = amasu * area
C
        ms(ibufn(1))=ms(ibufn(1)) + smass * ang1
        ms(ibufn(2))=ms(ibufn(2)) + smass * ang2
        ms(ibufn(3))=ms(ibufn(3)) + smass * ang3
C
C        ADDMAS = ADDMAS + SMASS * (ANG1 + ANG2 + ANG3)
C    ANG1 + ANG2 + ANG3 = 1.
        addmas = addmas + smass
      END IF
C=======================================================================
      RETURN