Civil Engineering Reference
In-Depth Information
c1
c2
intermediate reals
c3
c4
period
period of oscillation
element volume
volume
New scalar logicals:
consistent
.TRUE.
if element mass is “consistent”
New dynamic real arrays:
d1x0_pp
distributed old velocity vector
d2x0_pp
distributed old acceleration vector
d1x1_pp
distributed new velocity vector
d2x1_pp
distributed new acceleration vector
store_mm_pp
distributed element mass storage
temp_pp
distributed intermediate vector
distributed intermediate vector
vu_pp
distributed old displacement vector
x0_pp
distributed new displacement vector
x1_pp
The now familiar consistency of loops
elements_1
and
elements_2
with embed-
ded
gauss_pts_1
appears again. The parallel (3D) version has the lumped mass
matrix for 20-node elements hand-coded in place of the
elmat
routine used in the
serial (2D) case which used 8-node elements. The time stepping loop involves compa-
rable loop
elements_3
while the pcg section involves comparable loop
elements_4
.
Data are listed as Figure 12.40 with results as Figure 12.41 and performance statistics as
Figure 12.42.
nels nxe nze nip
2560 8 8 27
aa bb cc rho
0.125 0.125 0.125 1.0
e v
1.0 0.3
alpha1 beta1
0.0008 0.5
nstep npri theta
32 1 1.0
omega tol limit
0.01 0.0001 3000
Figure 12.40
Data for Program 12.9 example