Civil Engineering Reference
In-Depth Information
sigma=MATMUL(dee,MATMUL(bee,eld))
! IF(j/npri*npri==j)WRITE(11,'(I5,6E12.4)')iel,gc,sigma
END DO gauss_pts_2
END DO elements_4
END DO time_steps
CALL dismsh(loads,nf,0.05_iwp,g_coord,g_num,13)
CALL vecmsh(loads,nf,0.05_iwp,0.1_iwp,g_coord,g_num,14)
STOP
END PROGRAM p93
New scalar integers:
j
simple counter
loaded nodes
number of loaded nodes
ndof
number of displacement degrees of freedom per element
number of load function points
nlfp
maximum number of load steps
nls
output printed every
npri
time steps
npri
number of material properties
nprops
number of different property types
np types
node number at which time history is to be printed
nres
number of stress/strain terms
nst
nstep
number of time steps required
New scalar reals:
dtim
calculation time step
theta
time integration weighting parameter
holds time elapsed
t
time
tot load
accumulated load
New dynamic integer arrays:
etype
element property types
diagonal term locations
kdiag
New dynamic real arrays:
al
load steps at resolution of calculation time step
rhs “load” increment vector
ans
strain-displacement matrix
bee
coupling matrix
c
stress-strain matrix
dee
element nodal displacements
eld
integrating point coordinates
gc
element stiffness matrix
km
element conductivity matrix
kc
global stiffness matrix
kv
input load/time function
lf
phi0
used in element-by-element “gather” algorithm
phi1
used in element-by-element “scatter” algorithm
sigma
stress terms
