Civil Engineering Reference
In-Depth Information
IF(fixed_freedoms/=0)newlo(node)=storbp*value
CALL spabac(bp,newlo,kdiag); loads=newlo
IF(nod==4.AND.j==ntime)THEN; READ(10,*)nci
CALL contour(loads,g_coord,g_num,nci,13); END IF
IF(j/npri*npri==j)WRITE(11,'(2e12.4)')time,loads(nres)
END DO timesteps
STOP
END PROGRAM p82
New scalar integers:
nci
number of contour intervals
number of dimensions
ndim
number of integrating points
nip
number of nodes in the mesh
nn
number of elements in
x(r)
-direction
nxe
number of elements in
y(z)
-direction
nye
New scalar reals:
det
determinant of the Jacobian matrix
.
one
set to 1
0
Scalar characters:
dir
element and node numbering direction
element type
element
type of 2D analysis (
'plane'
or
'axisymmetric'
)
type 2d
New dynamic integer arrays:
g num
global element node numbers matrix
New dynamic real arrays:
coord
element nodal coordinates
shape function derivatives with respect to local coordinates
der
shape function derivatives with respect to global coordinates
deriv
shape functions
fun
integrating point coordinates
gc
nodal coordinates for all elements
g coord
Jacobian matrix
jac
permeability matrix
kay
points
integrating point local coordinates
weights
weighting coefficients
x coords
x(r)
-coordinates of mesh layout
y coords
y(z)
-coordinates of mesh layout
This program is for the analysis of 2D (
ndim=2
) first-order transient problems under
plane or axisymmetric conditions, and is closely based on Program 7.2 in Chapter 7. In
order to simplify the data however, the examples presented here use 4-node rectangular ele-
ments only (
element='quadrilateral'
and
nod=4
). The program includes graphics
subroutines
mesh
and
contour
which generate PostScript files containing, respectively,
images of the finite element mesh (held in
fe95.msh
), and a contour map of the excess
