Scribd
Upload
53 views

Di Ensio (3,3), GAUSS-R (2), GAUSS-S (2), B (3,8), EB (3,8), BT - B (4,8, Di Ension

This document provides code to compute the stiffness for a 4-node isoparametric surface element in FORTRAN. It defines element geometry, material properties, and loops over 4 Gauss points to calculate shape functions, the Jacobian, and stiffness matrix for the element. The code uses Gaussian quadrature to numerically integrate stiffness contributions from each Gauss point.

Uploaded by

aerosanth
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
53 views

Di Ensio (3,3), GAUSS-R (2), GAUSS-S (2), B (3,8), EB (3,8), BT - B (4,8, Di Ension

This document provides code to compute the stiffness for a 4-node isoparametric surface element in FORTRAN. It defines element geometry, material properties, and loops over 4 Gauss points to calculate shape functions, the Jacobian, and stiffness matrix for the element. The code uses Gaussian quadrature to numerically integrate stiffness contributions from each Gauss point.

Uploaded by

aerosanth
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

Parametric Elements 365

...........................................................................
FORTRAN to compute stiffness for 4-node, isoparametric,CO surface element
...........................................................................
DI~ENSIO ~
(3,3),GAUSS-R(2),GAUSS-S(2),B(3,8),EB(3,8),BT-B(4,8,
DI~ENSION PA~TS~F(4,8,8),S~FF(8,8)
INTEGER GPOINT
REAL Jll,Jl2,J2l,J22
1000 F O R ~ A T ( '',8(F9.5,1X))
* Define global coordinates, the material matrix, and thickness
DATA XA,XB,XC,XD,YA,YB,YC,YD/10,12,12,10,0,0,2,2/
DATA E(l,l),E(2,2),(3,3) /1,1,0.5/
TWICK=l
* Begin loop to evaluate integpand at each of 4 Gauss points
DO 8 0 1-1'2
DO 70 J=1,2
GPOINT=GPOINT+l
* Define the r and S coordinate values of the Gauss integration points
DATA GAUSS-R(l), GAUSS-R(2) /-0.5773502692,0.5~~3502692/
DATA GAUSS-S(l), GAUSS"S(2) /-0.5773502692,0,5773502692/
pe functions for current Gauss Point

C=-0.25
DRl=C*(l-S)
DR2=-lXC*(1-S)
DR3=-1*CX(1+S)
DR4=C*(l+S)

DS2=C*(l+R)
DS3=-l*CX(1+R)
DS4=-1*C*(l-R)
* Calc Jacobian components and determinant for current G.P.
Jll=DRl*XA+DR2*XB+DR3*XC+DR4*XD
Jl2=DRl*YA+DR2*YB+DR3*YC+DR4*YD
J2l=DSl*XA+DS2*XB+DS3*XC+DS4*XD
J22=DSl*YA+DS2*Y~~DS3*YC+DS4*YD
DTJ=(Jll*J22)-(Jl2*J21)

You might also like