Civil Engineering Reference
In-Depth Information
DO iel=1,nels; at=at+pt5*ell(iel)*(loads(iel)+loads(iel+1)); END DO
IF(j==ntime)press(1:)=loads(1:)
IF(j/npri*npri==j)WRITE(11,'(3e12.4)')time,(a0-at)/a0,loads(nres)
END DO timesteps
WRITE(11,'(/a,e10.4,a)')"
Depth
Pressure (time=",ntime*dtim,")"
WRITE(11,'(3e12.4)')0.0,press(1)
WRITE(11,'(2e12.4)')(SUM(ell(1:i)),press(i+1),i=1,nels)
STOP
END PROGRAM p81
Scalar integers:
fixed freedoms
number of fixed nodes
i
simple counter
iel
simple counter
iwp
SELECTED REAL KIND(15)
j
simple counter
nels
number of elements
number of degrees of freedom in the mesh
neq
number of nodes per element
nod
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 time steps required
nstep
time step number at which spatial distribution is to be printed
ntime
Scalar reals:
at
holds area beneath isochrone by Trapezoid Rule at time
t
holds area beneath isochrone by Trapezoid Rule at time
t
=
0
a0
calculation time step
dtim
set to 0
.
5
pt 5
10
20
set to 1
×
penalty
time integration weighting parameter
theta
holds elapsed time
t
time
set to 0
.
0
zero
Dynamic integer arrays:
etype
element property types
diagonal term locations
kdiag
nodes with fixed values
node
element node numbers
num
Dynamic real arrays:
bp
global “mass” matrix
ell
element lengths
kc
element conductivity matrix
kv
global conductivity matrix
loads
excess pore pressure values
