DYNA6 Manual With Examples
DYNA6 Manual With Examples
DYNA6
Dynamic Analysis of Foundations for the Effects of
Harmonic, Transient and Impact Loadings
Contents
0.0 SYSTEM REQUIREMENTS AND INSTALLATION .............................. 5
0.1 System Requirements ...................................................................... 5
0.2 Installation and activating the program ............................................. 5
0.2.1 Standalone version .................................................................... 5
0.2.2 Network version ......................................................................... 6
1.0 PURPOSE OF THE PROGRAM .......................................................... 8
2.0 TYPES OF FOUNDATIONS AND BACKGROUND THEORY ............. 8
2.1 Footing on Piles (Pile)....................................................................... 8
Single Pile Analysis ........................................................................... 10
Low frequency range......................................................................... 10
Pile tip condition ................................................................................ 14
Weak zone ........................................................................................ 15
Free length ........................................................................................ 16
Battered piles .................................................................................... 16
Pile groups ........................................................................................ 16
Static interaction factors .................................................................... 17
Frequency variation of interaction factors ......................................... 19
Groups Containing Battered Piles ..................................................... 22
Option for group stiffness output ....................................................... 22
Flexible Caps .................................................................................... 23
2.2 Footing on Homogeneous Half-Space (half-space) ........................ 24
2.3 Footing on Homogeneous Layer (stratum) ..................................... 25
2.4 Footing on Layered Medium (rigid body) ........................................ 25
2.5 Footing on Layer Overlaying Half-Space (composite medium) ....... 28
2.6 Flexible Rectangular Mat on Half-Space (mat) ............................... 31
2.7 Soil nonlinearity .............................................................................. 32
2.8 Soil material Damping..................................................................... 32
3.0 TYPES OF DYNAMIC LOADING ....................................................... 34
3.1 Loading Type 1: Transient .............................................................. 34
3.2 Loading Type 2: Random ............................................................... 36
50 GB of hard-drive space,
Note: Before the activation DYNA may produce a different product numbers
at each run, however, since they all refer to the same product, you dont need to
re-send the product number.
that the user read this part in order to understand the capabilities and limitations
of the underlying theory and the program.
Figure 2 .1a: Pile Footing (Rigid Cap)
taken as constant and the damping is taken as linear. The dynamic soil reactions
are shown schematically in Figure 2.1.2 with S1 and S2 representing stiffness and
damping, respectively, and the dimensionless frequency a 0 = R/vs where =
frequency is 0.3 for both stiffness and damping. The accuracy of the approach in
the low frequency range depends on the pile to soil, respectively). Figures 2.1.3
and Figure 2.1.4 show the vertical and horizontal pile stiffness compared to other
approaches.
b) Endbearing
Weak zone
The piles may have a weakened zone around them (see Figure 2.1.6),
which makes it possible to approximately account for the effects of imperfect bond
between the pile and the soil, pile sleeves and pile slippage. It also takes into
account the reduction of soil stiffness and increase of soil material damping due to
high strain. The weakened zone effect is calculated using the theory due to Novak
and Sheta (8). The original theory assumes the weak zone to be massless to avoid
wave reflections at the artificial boundary between the weak zone and outside soil
(28). The program allows the use of a weak zone mass participation factor (M.P.F
= 0 to 1) that represents the fraction of the weak zone mass to be added to the pile
mass at each layer. The M.P.F. should increase with the decrease of weak zone
thickness and the increase of weak zone shear modulus ratio (Gm/G) the
recommended maximum value for the M.P.F. is 0.75, with 0.25 to 0.5 being typical.
The weak zone reduces damping much more than stiffness.
Free length
The pile head may protrude from the ground or complete pile-soil separation
at the first layer may be assumed. The free length is accounted for by considering
the topmost layer as void (G1= 0,R1= 0, see Figure 2.1.1, case A). A free length
reduces stiffness, damping and group effects and may be used to account for
gapping, which has similar effects (Figs. 2.1.7a,b).
Battered piles
The effect of batter is handled approximately by calculating the pile
stiffnesses referenced to its local axes using the projected pile length on a vertical
plane. These stiffnesses are then transformed pile length on a vertical plane.
These stiffnesses are then transformed to the global axes system. Such
approximation is valid for slightly battered piles.
Pile groups
No Interaction Option
If the no interaction option is specified, the piles are assumed to act
independently of one another. Vertical and horizontal forces are distributed equally
on individual piles. The group stiffness is calculated by summing pile stiffnesses
using the appropriate transformation to the C.G. of the system. The No Interaction
option is of interest only for comparison with the Interaction option and for
academic purposes.
Interaction Option
If pile-soil-pile interaction is considered, as it should be for closely spaced
piles, the group stiffness and damping are calculated using the superposition
method described in El Naggar and Novak (24). The interaction between each two
piles is used to formulate the group complex flexibility matrix from which the group
stiffnesses and damping are deduced. The group effect is considered separately
for the vertical translation, the horizontal translation in the XZ plane and the
horizontal translation in the YZ plane (Z being the vertical axis). The complex
interaction factor is given by equation 2.1.1.
ij
Figure 2.1.7b: Development of Gap Separating Pile from the Soil with
Number of Cycles for Different Load Amplitudes. Steel Tube Pile with
610mm O.D., Stiff Clay (Swane & Poulos, 1984)
pile values. Negative stiffness may occur depending on the frequency and the
fundamental pile spacing to diameter ratio. Figure 2.1.8 shows an example of
dynamic group stiffness and damping in the vertical and horizontal directions. (The
differences observed are not of great practical significance. Most practical cases
feature very small a0', 0.4 or less.) For soils with high vertical nonhomogeneity (e.g.
Gibson soil), the above mentioned approximations yield less accurate results for
the group stiffness in the horizontal direction.
Figure 2.1.9: Pile Group with Flexible Cap or No Cap. Stiffness and
Flexibility Matrices are Returned Referring to all Pile Heads.
Flexible Caps
When the pile cap is thin relative to its plan dimensions, it may bend
significantly under vertical loads or moments in the vertical planes; however, the
cap horizontal and torsional response can still be calculated employing the rigid
cap option, CAPRIGID.
The flexible cap analysis for vertical loads and rocking moments can be
conducted in two ways:
(a) The user can calculate the complete stiffness matrix of the piles without
a cap using DYNA6 and introduce it into his own finite element
subroutine for the analysis of the superstructure (or flexible cap)
(b) The user can employ the DYNA6 option CAPFLEXIBLE that analyzes
the pile group together with a finite element model of the flexible cap.
The loads are harmonic vertical forces or moments in the vertical planes
Cap contact with the soil and its embedment are not considered.
the footing in which the soil properties may differ from those of the outer part of the
layer. This feature makes it possible to take an approximate account of the
increased confining pressure under the footing.
This option is treated in the program as a special case of a thick pile as
shown in Fig. 2.4.1 and outlined in (11) and should not be applied to shallow
foundation with E/r less than about 5 (E = embedment depth). Deformations of the
footing due to shear, bending and torsion are considered. Also included is the
inertia effect of the body, which can result in strong variation of stiffness with
frequency and even in negative stiffness constants. If this option is used just to
obtain the stiffness and damping of the body without the effect of its mass, the
program can be run with a very small, nominal, value of mass (i.e. with unit weight
approaching zero). The body may protrude from the ground. The free length is
accounted for by considering the natural frequency of the layer applies to the
RIGID-BODY option. The tip condition refers to the soil below the lowest layer. The
stiffness and damping constants returned refer to the top of the footing if the height
of the center of gravity is input as zero. For all options except RIGID-BODY, the
height of the center of gravity is measured from the base of the footing to the centre
of gravity (C.G.) of the whole machine-footing system and is positive if the C.G.
lies above the base. For the RIGID-BODY option, the height of the center of gravity
is measured from the top surface of the footing to the C.G. of the superstructure
(see example 11), as with piles.
a) Surface Footing
b) Embedded Footing
d) Non-Uniform Profile
Vs Gs / s
Where Gs is shear modulus and s is mass density, unit weight R, and
material damping. In DYNA6, material damping is defined as D = tan = 2 where
is the loss angle and is the damping ratio. Hysteretic, frequency independent
material damping is assumed, as it is more realistic than viscous damping. This
results in strong growths of equivalent viscous damping with frequency
approaching zero as is schematically depicted in Fig. 2.8.1. Thus, the role of
material damping can be dominant in the low frequency range, particularly for soil
strata and endbearing piles, vibrating with frequencies lower than the fundamental
frequency of the soil layer (Fig. 2.8.1a), a situation where no geometric damping
exists.
b) Halfspace
The program returns the complete time history of the response in six
degrees of freedom as well as its maxima.
ui g *u
3.
1
In which g is the peak factor and Su is the standard deviation (root-meansquare value) of the response. The peak factor is calculated from the spectra. The
period T may typically be 20 to 60 seconds for earthquakes and 600 to 3600
seconds for machine or wind loading.
components the response has to be calculated for each component separately and
the results added.
In DYNA6, the harmonic excitation can be of two types: the constant
amplitude excitation and the frequency dependent or quadratic excitation. The
constant amplitude excitation is defined as:
P(t ) P0 * sin(t )
3.
P(t ) me * e * 2 * sin(t )
3.
P0 me * e * 2
3.
4
Conversely, Eq. 3.4 can be used to evaluate the excitation product m ee if
the force amplitude, P0, is given for a certain operating frequency. The difference
between the two types of harmonic excitation is further clarified in section 7.4.1.
When the blows are centric, both of these types of hammer foundations can
be represented by a two mass model with two degrees of freedom as shown in
Fig. 3.3a. The vertical displacements v1(t) and v2(t) are return by DYNA6. For the
foundation with the inertial block (Fig. 3.2a), they describe the motion of the anvil
(v1) and the block (v2) respectively. For the directly sprung hammer (Fig. 3.2b),
v1(t) is the motion of the anvil and v2(t) the motion of the trough.
For a hammer with the inertial block, the anvil pad is described by its area,
thickness, Youngs modulus and material damping ratio, which is presumed to be
hysteretic; constants k1 and c1 are evaluated in the program. For the directly
sprung hammer (Fig. 3.2b), total stiffness constant of the springs, k1, and the
damping constant of the dampers, c1, have to be provided. Constants k2, c2
describe the soil or pile properties and are calculated for both types of hammers
within DYNA6 for any of the foundation types shown in Figs. 2.1 to 2.5.
Two types of shock loading are available. For a shock whose duration t is
much shorter than the fundamental period of the foundation, the response is
treated as vibration due to initial velocity issued to the anvil by the blow of the head
(Fig. 3.4a). The initial velocity C of the anvil is to be established by (23):
C (1 kr ) *
m0
*C0
m0 m1
3.
b)
b) 3 or 6 DOF
a) Short Shock
b) Half Sine
4.0 UNITS
Any units can be used but they must be consistent. The output will come
out in the same units in which the input is provided. For the English system, the
units are lb, slug, ft, s with mass in slug = weight/gravity (lb/ft/s2) and gravity = 32.2
ft/s2. For the SI system, the units are N (for force), kg (for mass), m, s or kN, Mg =
1000 kg (tonne), m, s; gravity = 9.81 m/s2.
For random loading, the power spectrum is read in as (force) 2 at circular
frequencies in rad/sec.
When the excitation is due to ground motion xg(t) , the force is defined as
mass(m) x ground acceleration (xg) for harmonic or transient loading: for random
loading, the load power spectrum is m2S
xg()
acceleration.
Imperial System
Time
s (second)
s (seconds)
Length, distance
ft
Gravity constant
9.81 m/s2
32.2 ft/s2
Mass
kg
slug = lb/32.2
kg.m2
slug.ft2
Unit weight
N/m3
lb/ft3
m/s
ft/s
Frequency
rad/s
rad/s
Forces
lb
Moments
N.m
lb.ft
Displacements (output)
ft
English (Imperial): ft, lb, slug or ft, kip (= 1,000 lb), kslg (= 1,000 slug) with
these options the compatible gravity acceleration is entered automatically as either
9.81 m/s2 for SI or 32.2 ft/s2 for the English system. With the user-selected units,
it has to be entered separately.
5.0 NOTATION
For loading types 1 to 3, the input forces and moments and the response
components follow the sequence of labels shown in Fig. 5.1. For example, labels
1 and 2 indicate horizontal directions along the X and Y-axis. The force and
moments are defined with regard to the C.G. of the system.
The forces and translations are positive if they follow the positive direction
of the corresponding axis. Moments and rotations are positive if they act in the
clockwise direction when watched from the origin in the positive direction of the
pertinent axis. Thus, e.g., two positive horizontal forces p 1, p2 acting above the
C.G. in the directions +X, +Y result in moments M4 < 0 and M5> 0. Stiffness
constants, being by definition external forces, follow the same sign convention;
consequently, the cross-stiffness constants of symmetrical foundations differ in
signs.
For the two mass foundations (Fig. 3.3a), the vertical displacements are
positive when they act downward.
6.0 REFERENCES
1. Veletsos, A.S. and Wei, Y.T. (1971) - Lateral and Rocking Vibration of
Footings, J. Soil Mech. And Found. Div., ASCE, SM9, September, pp. 1227- 1248.
2. Veletsos, A.S. and Verbic, B. (1973) - Vibration of Viscoelastic
Foundation, J. Earthquake Engrg. And Struct. Dyn., Vol. 2, pp. 87-102.
3. Veletsos, A.S. and Nair, V. V. D.(1974) - Torsional of Vibration of
Viscoelastic Foundaation, J. Geotech. Div., ASCE, Vol. 100, No. GT3, March, pp.
225-246.
4. Novak, M. And Beredugo, Y.O. (1972) - Vertical Vibration of Embedded
Footings, J. Soil Mechanics and Foundation Division, ASCE , SM12, December,
pp.1291-1310.
5.Beredugo,Y.O. and Novak, M. (1972) - coupled Horizontal and Rocking
Vibration of Embedded Footings, Canadian Geotechnical Journal, Vol. 9, No. 4,
pp. 477-97.
6. Novak, M. And Sachs, K. (1973) - Torsional and Coupled Vibrations of
Embedded Footings, Inter. J. Earthquake Engrg. And Struct. Dyn., Vol.2 No.11,
p.33.
7. Novak, M. (1974) - Effect of Soil on Structural Response to Wind and
Earthquake, Inter. J. Earthquake Engineering and Struct. Dyn., Vol. 3, No.1, pp.
79-96.
8. Novak, M. And Sheata, M. (1980) - Approximate Approach to Contact
Problems of Piles, Proc. Geotech. Engrg. Div. ASCE National Convention
Dynamic Response of Pile Foundations: Analytical Aspects, Florida, October, pp.
53-79.
9. Kausel, E. And Ushijima, R. (1979) - Vertical and Torsional Stiffness of
cylindrical
Footing,
Civil
Eng.
Dept.
Report
R79-6,
MIT,
Cambridge,
Massachusetts.
10 Novak, M and Aboul-Ella, F. (1978a) - Impedance Functions of Piles in
Layered Media, Journal of the Engineering Mechanics Division, ASCE, Vol. 104,
No EM6, Proc. Paper 13847, June, pp. 643-661.
International
Journal
of
Soil
Dynamics
and
Earthquake
24. El Naggar, M.H. and Novak, M., 1996. Nonlinear analysis for dynamic
lateral pile response. J. of Soil Dynamics and Earthquake Engineering, Vol. 15,
No. 4, pp. 233-244.
25 Randolph, M.F. and Poulos, H.G. (1982). Estimating the Flexibility of
Offshore Pile Groups, Numerical Methods in Offshore Piling Proceedings of the
2nd International Conference, University of Texas, Austin, TX.
26. El Sharnouby, H. and Novak M. (1986). Flexibility coefficients and
Interaction Factors for Pile Group Analysis, Canadian Geotechnical Journal, Vol.
23, No. 4, pp. 441-450.
27. Kaynia, A.M. and Kausel, E. (1982). Dynamic Behaviour of Pile
Groups, Conference on Numerical Methods in Offshore Piling. Univ. of Texas,
Austin, TX, pp. 509-532.
28. Novak, M. And Han, Y. (1990). Impedance of Soil Layer With Boundary
Zone, Journal of Geotechnical Engineering, Vol. 116, No. 6, June, pp. 1008-1014.
29. Whittaker, W.L. and Christiano, P. Dynamic Response of Flexible
Plates Bearing on An Elastic Half-Space, RP-125-9-79, Dept. Of Civil Eng.,
Carnegie-Mellon University.
30. El Sharnouby, B. and Novak, M. (1985). Static and Low Frequency
Response of Pile Groups, Canadian Geotechnical Journal, vol. 22, No. 1, pp. 7994.
31. Swane, I.C. and Poulos, H.G. (1984). Shakedown Analysis of Laterally
Loaded Pile Tested in Stiff Clay, Proc. 4th Australian-New Zealand Conf. On
Gemech., Perth, Vol. I, pp. 165-169.
32 Novak, M. (1985) Experiments With Shallow and Deep Foundations,
Proc. Of ASCE Symposium on Vibration Problems in Geotechnical Engineering,
Detroit, Mich., pp. 1-26.
33. Crouse, C.B., Price, T. And Mitchell, R., (1992). Evalutaion of Methods
to Estimate Pile Foundations Stiffnesses for Bridges, Proc. 8th U.S.-Japan Bridge
Engineering Workshop, Chicago, Illinois, pp.14.
34. Han, Y and Vaziri, H. (1992). Dynamic response of Pile Groups under
Lateral Loading, Journal of Soil Dynamics and Earthquake Engineering 11, pp.
87-99.
35 El-Marsafawi, Han, Y. And Novak, M., (1992). Dynamic Experiments on
Two Pile Groups, Journal of Geotechnical Engineering, Vol. 118, No. 4, pp. 576592.
correct
1.0E09
correct
1.0E 9
incorrect
Problem Description:
The following pages outline the various keywords used to describe the
problem to the DYNA6 system.
Problem Title:
The problem may be given a title by clicking on the Title on the menu bar.
The title can be a maximum of 80 characters long. The title is used on the all the
printed output of DYNA6.
The limitations on the title are:
(1) The first character of the title must be alphabetic (i.e. no digits)
(2) The title must not contain equals (=), commas (.) or dollar sign ($).
Settings:
Next to the Title menu is the Settings menu. Under this menu the user can
set the following options:
1) Output Flags
2) Analysis Constants
3) Units
4) Frequency Units
1) Output Flags:
Data Echo:
DYNA6 will generate an echo indicating the interpretation of the input data.
Clicking on the Echo Input Data box may enable this echo.
Dimensionless Results:
The response amplitudes calculated for Harmonic Loading may be nondimensioned by clicking on the Dimensionless Results box.
Note: This command will produce dimensionless amplitudes that are equal to
dynamic amplification factors and are defined by the following formulae:
u 0( )
u 0( min)
Dimensionless rotations:
Quadratic excitation:
F 0( )
A
F 0( min )
m
u 0( )
A
* u0
mee
u 0 ( )
Dimensionless translations:
Dimensionless rotations:
I
Mi
mee * ( )
Pi
*F0
ustat
4*G * r0
(1 )
where r 0
area
(For the Flexible Cap option, G and are the value at the pile tip.)
Note: If the dimensionless amplitude cannot be evaluated because Ust is zero, true
amplitude cannot be evaluated because Ust is zero, true amplitudes are
returned.
This keyword will instruct the DYNA6 system to output the constants stored
in the Foundation Stiffness and Damping Matrices. These constants are referred
to the C.G. of the footing specified in the input and incorporate the mass of soil
and piles. Because of the inertia effect of these masses, the stiffness constants
can be negative, particularly for higher frequencies, soil Poisons ratio ->0.5 and
massive piles. For the option RIGID-BODY, the mass of the embedded footing is
also incorporated in the stiffness and damping constants.
NOTE: For the flexible mat (MAT option) and the flexible cap
(CAPFLEXIBLE option), the keyword MATRIX is not applicable.
2) Analysis Constants:
Damping Safety Factor:
If the applied frequency is close to the resonant frequency of the foundation
or, in the case of foundations supporting a superstructure, close to the
superstructure resonant frequency, the response strongly depends on the
magnitude of the damping constants. To reduce the possibility of damping
overestimation and thus response underestimation, an optional user specified
safety factor may be applied to the damping constants by entering a value in the
Damping Safety Factor box. Where Damping Safety Factor (S) is any
positive number greater than 1.0 defined by:
c(calculated )
c(used )
Example:
Damping Safety Factor = 2.0
This will divide the damping constant by 2.0. A value of S between 1.0 and
2.0 is recommended, with 2.0 suitable especially for pile foundations and deeply
embedded foundations.
Notes:
1. The damping safety factor is not applicable to the flexible mat
(MAT) and the CAPFLEXIBLE options.
2. The damping safety factor is not incorporated in the pile group
stiffness matrices at the pile heads printed using the keywords FLEXIBLE,
DYNFLEX, LATERAL OR DYNLAT (see Section 7.2.1). It is incorporated in
the rigid pile cap matrix printed using keywords MATRIX, TABULATE and
DYNSTIF.
Gravitational Constant:
The gravitational constant is automatically selected when a set of units is
chosen, for the SI units the value is 9.81 m/s2 and for the Imperial system the value
is 32.2 ft/sec2. However, a different value may be entered by the user in the Gravity
Acceleration box. The dimension of the gravitational constant must be consistent
with the dimensions of the other input data and determines the dimensions of the
output. For example with gravity = 32.2 ft/sec2, the displacement amplitudes of the
response are returned in ft; with gravity = 9.81 m/s2 the displacement amplitudes
are returned in m. Amplitudes of rotations are always in radians.
3) Units:
By default, DYNA6 does not include any units in the echo of the input data
or the results. To indicate that a given set of units is to be used in the output file
the user clicks on either the SI Units or Imperial Units button. If it becomes
necessary to use other alternate units, click on the Other button. Once this option
is selected, a set of four boxes becomes available so that the user may choose:
1) Force Units
2) Length Units
3) Mass Units
4) Frequency Units (see next section)
The user is given the option to choose either the standard SI or Imperial unit
for that particular parameter, or to input a new one.
Maximum lengths:
Length: 2 characters (example m or ft)
Force: 3 characters (example N or lb)
Mass: 4 characters (example kg or slug)
Note: The inclusion of units specification for output purposes has no effect on the
numerical values of the results which are governed by the gravity constant.
The user must ensure that the input units are consistent and that they
correspond to the chosen value for the gravity constant.
4) Frequency Units:
The default units for frequency are rad/s. The program also supports the
use of input and output frequencies in Hertz (cycles per second) and rpm
(revolutions per minute.
Notes:
1. The LOAD option RANDOM (SEE SECTION 7.4.3) IS ONLY ALLOWED UNITS
OF rad/s.
2. The plots of the stiffness and damping constants that are produced by the
GRAPH utility would always use RAD/S for the frequency axis, irrespective of
the units used in input and output. The program handles such transformation of
units internally.
Figure 7.1: Notations for Pile Foundation Option with Rigid Cap
(CAPRIGID)
1) PILE CAP
If the Flexible Cap option is selected, clicking on the Pile Cap button will
bring up the Mat Foundation information window. Clicking on the Data button
opens up the Flexible Mat Data window. There are a number of parameters to be
entered by the user:
1) Mat Properties
a) Length LX
b) Length LY
c) Thickness
d) Mass Density
e) Poissons Ratio
f) Youngs Modulus
2) Response at Nodes (maximum of 5)
3) Output Options
a) FLEXIBLE
b) DYNFLEX
4) Number Of Nodes
a) In the X direction NX
b) In the Y direction NY
5) Node Spacing in X
a) Equal
If the Rigid Cap option is selected, then clicking the Pile Cap button opens
the Footing Base Dimensions and Data window. There are three parameters one
can input in this window. For a circular foundation, the radius is entered and for a
rectangular foundation, the lengths in the x and y directions (Lx,Ly) are entered.
The coordinates of the base centre (Xc,Yc,Zc) are also entered. The mass moments
of inertia (Ixx,Iyy,Izz,Ixy,Ixz,Iyz) can be entered by the user if they have been previously
calculated or the user can click on the Calculate box and use the 3DVIEW utility to
calculate them (see ch.9) Once all the data has been entered, click on the OK box.
Notes:
a) If the foundation is not rectangular (or oddly shaped) then a
representative length, Lx or Ly, should be selected and then the other
length should be calculated so that the area of the rectangular footing
would be equal to the area of the actual footing.
b) The coordinates of the base centre should be calculated after the
determination of the coordinates of the C.G. of the system, either by the
user or through the 3DVIEW utility.
c) The last three items in mass (Ixy,Ixz,Iyz ) are products of inertia. They are
important only for grossly asymmetrical foundations; in normal cases,
they can be taken as zeroes.
d) For rigid-body option the masses and mass moments of inertia are those
of the foundation itself are accounted for through the unit weight of the
rigid body.
2) SIDE LAYERS
This option is only available if the Embedded option is chosen and is
unavailable for the Flexible Cap pile foundation and mat foundation option. Clicking
on the Side Layers button brings up the Data for Cap/Footing Side Layers window.
By default, the weak zone around a footing is not considered, but can be
incorporated into the calculations by selecting Weak Zone in this window. The user
can enter the following side layer parameters in this window:
1) Side Layer Thickness
2) Side Layer Shear Wave Velocity
3) Side Layer Unit Weight
4) Poissons Ratio
5) Damping
Once all the data has been entered click OK to return to the Pile Data
window.
3) PILES
Clicking this button brings up the Pile Head Arrangement for Rigid
Pile Cap window. In this window, the user can input the following parameters by
clicking on the appropriate box:
1) Pile Head Condition
a) Fixed
b) Pinned
Note: Pinned head piles are not allowed to carry moments about X
or Y-axes.
2) Pile Tip Condition
a) Floating Tip
b) Endbearing
3) Pile Properties:
a) Pile Length
b) Pile Unit Weight
c) Static Load
4) Pile Material Properties:
a) Youngs Modulus
b) Poissons Ratio
c) Damping
d) Rigidity Coefficient
The user can manually input all the pile coordinates by clicking on Individual
on the menu bar and typing in the coordinates (or pasting an outside file) in the
Coordinates for Individual Piles window. By clicking on the menu bar Mesh, the
Mesh Generation window appears and the user can simply type in the coordinates
of the two opposite corners in the appropriate boxes and the program will formulate
the proper pile mesh. The mesh can be titled and saved for later use. The user has
the following options when developing the mesh:
1) Shape of the Mesh
a) Rectangular
b) Circular
2) Total Number of Piles
a) in X,Y Directions (rectangular mesh)
b) In R, Directions (circular mesh)
3) Pile Spacing (in X,Y or R, directions)
a) Equal (default)
b) Unequal (the user can manually set the spacing by
clicking on the Edit button)
4) Generation Direction
-
Once all the information about the mesh has been entered, click on the
Apply box to return to the Pile Head Arrangement window. If necessary, more
individual piles can be added to the newly created mesh by using the Individual
window, these piles will be added to the end of the list of piles.
Note : The maximum number of piles allowed by the program is 86 for the PC
version operating within the conventional memory (640 KB), and 1000 for
the extended/virtual memory version. The actual number that may be
analyzed depends on the available free memory. Refer to Section 8 for more
details.
If the Flexible cap option has been chosen, clicking on the Piles button will
open the Pile Head Arrangement for Flexible Cap window. The user does not need
to input the locations of the piles but only the material properties for the piles as
follows:
5) Pile Head Condition
a) Fixed
b) Pinned
Note: Pinned head piles are not allowed to carry moments about X
or Y-axes.
6) Pile Tip Condition
a) Floating Tip
b) Endbearing
7) Pile Properties:
a) Pile Length
b) Pile Unit Weight
c) Static Load
8) Pile Material Properties:
a) Youngs Modulus
b) Poissons Ratio
c) Damping
d) Rigidity Coefficient
If any of the piles in the mesh are battered, click on the menu bar Batter to
bring up the Data for Battered Piles window. In this window the user can input the
following information:
1) Pile Number (from mesh)
2) Angle
3) Angle
Click OK to return to the Pile Head Arrangement window.
Clicking Figure on the menu bar brings up a graphical representation of the
mesh. Holding the mouse pointer over one the piles will highlight it and bring up its
coordinates; this diagram can also be printed out. Clicking on the menu bar Print
Coord. will allow the user to print out the pile list along with their respective
coordinates.
4) PILE SOIL
Clicking on the Pile Soil button opens the Data for Pile-Soil System window.
In this window the user inputs information about the pile and soil elements. The
user has the option of choosing:
1) Soil Medium Type
a) Layered
b) Parabolic
2) Pile-Soil-Pile Interaction
a) No Interaction
b) Interaction
3) Weak Zone Interaction (default is no interaction)
For a layered soil (maximum of 30 layers), the user needs to calculate and
input the following parameters:
1) Layer Depth
2) Pile X-radius
3) Pile Y-Radius
4) Pile Area
5) Pile Y Inertia
6) Pile X Inertia
7) Pile Z Inertia
If the Parabolic option is chosen, then the parameters need only be entered
for one layer. The user then clicks on the Soil Elements tab at the bottom of the
window and enters the following information into the spreadsheet:
1) Soil Shear Wave Velocity
2) Soil Unit Weight
3) Poissons Ratio
4) Damping (material damping of soil is defined here as D = tan=
2 where = loss angle and = damping ratio)
5) SOIL BELOW
Clicking on the Soil Below button allows the user to enter the soil properties
below the pile. The following parameters are then typed in the appropriate boxes:
1) Soil Shear Wave Velocity
2) Soil Unit Weight
3) Poissons Ratio
4) Damping (material damping of soil is defined here as D = tan=
2 where = loss angle and = damping ratio)
6) OUTPUT
Clicking on the Output button opens up the Output Options for Pile
Foundation window. The user can select the following output options:
1) Print Piles Forces Due to:
a) Unit Displacements (selecting this box opens up the Data
for Applied Unit Displacements window)
b) Actual Loads (selecting this box opens up the Data for
Applied Actual Loads window)
2) Piles Vertical Stiffness
a) Print in Output File
b) Print in DYNFLEX File
3) Piles Lateral Stiffness
a) Print in Output File
Notes: 1) The keyword FLEXIBLE prompts the printout of the vertical complex
stiffness matrix, n X n, of the group of piles, referred to the individual pile
heads, for the analysis of flexible caps or superstructures. The keyword
DYNFLEX places the same stiffness matrix into a separate file called
DYNFLEX in the current directory. This file is not a part of the output file.
The vertical stiffness matrix is relevant for flexible cap vibration modes
such as the one depicted in Fig.2.6.
2) The keyword LATERAL prompts the printout of the 2n X 2n group
complex stiffness matrix, referred to individual pile heads, pertinent to the
response in the lateral XZ plane (horizontal, rocking and coupling terms),
followed by the corresponding matrix for the YZ plane.
3) The keyword DYNLAT places the same matrices in a separate file called
DYNLAT in the current directory. If the piles are numbered 1,2,3,...n and
the horizontal translation and the rotation at the pile head of pile number i
are denoted by Ui and Fi, respectively, the arrangement of the degrees of
freedom in the lateral stiffness matrices is as follows:
{U1x,F1y,U2x,F2y,...,Uix,Fiy,...,Unx,Fny}
for the XZ plane, with each item having a real and an imaginary part, and
{U1y,F1x,U2y,F2x,...,Uiy,Fix,...,Uny,Fnx}
for the YZ plane
4) The keywords FLEXIBLE, DYNFLEX, LATERAL and DYNLAT should be
used with caution because for a large pile group and a large number of
frequencies, the size of the files and printouts may become excessive.
Clicking on the Resultants button allows the user to input the coordinates of
up to 10 points other than the C.G. The program then calculates the translations
of the these points and places them in the output file.
b) Embedded
c) Poissons ratio
d) Material Damping
3) Resultants (as described in Section 7.2.1, this allows the user to
track the translation of up to 10 points, other than the C.G.)
Notes:
b) Embedded
be subject to bending deformations. The data is entered into the data windows in
the same way as described in Section 7.2.1. Clicking on the Structure button opens
up the Structure and Foundation Data window, this window is quite similar to the
Footing Base Dimensions and Data window. The user can input the following
information into the data window:
1) Structure Mass and Inertia
a) Total Mass
b) Moments of Inertia (Ixx,Iyy,Izz,Ixy,Ixz,Iyz)
2) Rigid Body Top Surface
a) Length Lx
b) Length Ly
3) Top Surface Centre (Xc,Yc,Zc)
Next, the user clicks on the Rigid Body button to bring up the Rigid Body
Foundation Data window. The following data is entered in this window:
1) Rigid Body Properties
a) Total Depth (depth of footing plus soil column underneath)
b) Unit Weight
c) Static Load (approximate)
2) Rigid Body Material
a) Youngs Modulus
b) Poissons Ratio
c) Damping Coefficient
d) Rigidity Coefficient (See footnote 3 in section 7.2.1)
3) End Condition
a) Floating (default)
b) Endbearing
ii) Embedded
Note On Limitations:
The impedance functions are exact for the radio of layer thickness to
halfwidth of the square footing (H/a) equal to 0.5, 1,2,3 and 4 for uniform layers
(Fig. 7.6b(i)) and equal to 2,3,4,5 and 10 for non-uniform layers (Fig. 7.6b(ii)). If
the ratio (H/a) doesnt coincide with one of the above values the program chooses
the closest (H/a) ratio available (interpolation is not implemented because of the
strong non-monotonic variations at high frequencies).
In the composite-medium option, accurate values of stiffness and damping
are used at frequencies equal to 0.10, 0.25,0.50..., 4.75 and 5.0 times (Vs /a)
where Vs is the shear wave velocity at footing base level and a is halfwidth of the
square base(or the equivalent square base). For a frequency less than 0.10 Vs/a,
the program uses the minimum value (0.10 Vs/a) and for frequencies in the range
(0.10-5.0) Vs/a, a linear interpolation is implemented. If the frequency is greater
than 5 (Vs/a) the program uses the maximum value of 5 (Vs/a).
The response phase shift, F, relative to the excitation force may be printed
using the optional keyword
Notes on Soils:
1) Poisson's ratio of the half-space is assumed 0.33 in this option.
2) Material damping of soil is assumed 0.03 and 0.05 for the layer and the
half-space, respectively.
3) Two values for Poisson's ratio of the layer are available 0.33 and 0.45. If
a different value is entered the program sets it to the closest one. (Interpolation is
not implemented in the program because of non-monotonic variations)
4) Vs' is the layer shear wave velocity at footing base level. (Fig. 7.6b)
5) Three values for the shear wave velocity ratio are available 0.8, 0.6, and
0.3. If a different value is entered the program sets it to the closest one.
6) The ratio of unit weight of the half-space to that of the layer (p/p') is
assumed 1.13.
Figure 7.7: Notations for Flexible Mat on Halfspace (Option MAT) and
Pile Foundation with Flexible Cap (CAPFLEXIBLE Option)
a) Elevation
Plan
Notes:
1) The nodes should be labeled from left to fight and from bottom to
top as shown in Fig. 7.7. If the response is required at more than five
nodes, more than one run must be used.
2) The response has real and imaginary parts. Returned are
the absolute value and the phase (optional) from which both the real
and imaginary parts may be obtained.
3) The total number of nodes (NX.NY) allowed by the program
is 144 for the PC version operating within the conventional memory
(640 KB), and 1650 for the extended/virtual memory version. The
actual number that may be analyzed depends on the available free
memory. Refer to Section 8 for more details.
4) The keywords DYNFLEX and FLEXIBLE should be used
with caution as they may produce an excessive amount of data.
a) In Phase (default)
b) Out of Phase (click on the Phase Shift button to edit the
phase shift)
Notes:
forcei (mee)i
pi
op
where i x, y, z
momenti (meer mi )
Mi
op
frequency op. Usually, the magnitudes of pi and Mi are given for the
operating frequency op and have to be divided by op to get the
input. Finally, r = the respective arm of the excitation force with
regard to the centre of gravity of the system and mi is the unbalanced
moment independent of pi. (Often, there are no unbalanced forces;
pi but there is unbalanced moments mi.)
c) All the phase sifts are entered in degrees. For the quadratic
case, the vertical forces are usually 90 degrees out of phase from
the horizontal forces.
d) This option to specify input forces phase shifts is available
only for harmonic loads and rigid footings or rigid caps.
e) Due to numerical round off errors, some phase shifts that
should be zero show random fluctuations in sign and value; these
phase shifts should ignored. The phase shifts are returned in the
range 0 to 2or 0 to -2 as shown in the figure.
f) The number of frequencies should not exceed 1024.The
resultant translation at any point on the machine-foundation system
is due to a combination of both translations and rotations, at the
system's centre of gravity.
g) If the option Dimensionless is used, the resultant
translations are not returned in dimensionless form but are given in
actual values. The option Resultants does not apply to flexible mats
and caps (Mat, Capflexible). All output displacement amplitudes are
"single amplitudes," i.e. displacement amplitudes from the mean
position, not peak-to-peak values.
Each of the above files has an associated text file with the same name but
with an extension .DTX that is also used by the GRAPH utility. Use the keyword
ALL with care; the length of the output file may be increased by up to 1300 lines.
Note: a) If the number of points is a power of 2 (e.g. 256 or 512 or 1024), the
trailing zero load values that are needed to alleviate the problems
associated with the periodicity implied in the Fourier Transform approach,
should be included in the data provided by the user. Otherwise add an extra
zero point (e.g. number of points = 257 or 513 etc.) and the program would
automatically add zeros until the number of points is equal to the next power
of 2.
b) If there is a significant static (DC) component in the transient input
force, as for example in the short moments of a power generator, it should
be separated from the signal for the sake of FFT accuracy and treated
separately as a static or very low frequency harmonic effect.
NOTE:
(a) Keyword PAD is used for two mass foundations with anvil pad
(Fig. 3.2a); Keyword SPRING is used for directly sprung hammers
with anvil springs and dampers (Fig. 3.2b).
(b) Initial velocity of anvil is calculated from Eq. 3.5. This is the
natural frequency of the foundation
3DVIEW
A Program For Machine - Foundation Block Calculation for Use in DYNA6
3DVIEW INTRODUCTION
In DYNA foundations are assumed to be rigid bodies. A single, rigid body
replaces the machine, its foundations, and other attached components with the
reference point being the centre of gravity of the whole system. The excitation
forces, which may act at any point, must be transferred to the system centre of
gravity.
UNITS:
Any system of units can be used but it must be consistent and the same as
in the program, i.e. force, mass, and length in SI units N, kg, m or in Imperial units
lb, slug (weight/g, lb/ft/s2, with 9 = 32.2 ft/s2), ft.
Rectangular Blocks: Each rectangular block is defined by its three side lengths,
three coordinates of its centre of mass, and its density. The
sides of a rectangular block should be parallel to the three
coordinate axes.
Circular Cylinders: Each circular cylinder is defined by its diameter, length, three
coordinates of its centre of mass, density, and the direction of
its axis. The axis of the cylinder should be parallel to one of
the coordinate axes X, Y, Z.
Lumped Masses:
Accounting for Holes: To account for holes for both lumped masses and complex
bodies, the mass is entered with a negative sign. For blocks
and cylinders the density is entered with a negative sign.
Excitation Forces:
The sign convention for the forces shown in Fig. 9.1 is the same as that
used in the main program, i.e. the right hand rule controls the moments. Automatic
change in some signs accompanies the transformation from 3DVIEW axes to
DYNA axes accomplished by 3DVIEW.
Rectangular Blocks, 7 numbers for each block: Input data for rectangular
blocks. The format is:
Lxi Lyi Lzi Xi Yi Zi p; for i = 1,NB (No. of rectangular blocks)
where
Lxi Lyi, Lzi = lengths of rectangular block sides parallel to the X,Y and Z axes,
respectively;
xi, yi, zi = X, Y, and Z coordinates of the block centre of mass.
= mass density of block material.
- Each block data should be entered on a separate line.
- If NB = 0 skip this data block.
Cylindrical Blocks, 6 real numbers and one integer for each cylinder: Input
data for circular cylinders. The format is di li pi xi yi zi NA
where
di, li, pi = diameter, length, and mass density of the circular cylinder,
respectively.
xi, yi, zi = X, Y, and Z coordinates of the cylinder centre of mass.
NA = Axis code = 1 cylinder axis is parallel to X - axis.
= 2 cylinder axis is parallel to Y - axis.
= 3 cylinder axis is parallel to Z - axis.
-Each cylinder data should be entered on a separate line.
- If NY = 0 skip this data block.
Lumped Masses, 10 real numbers: Input data for lumped masses. The
format is:
mi xi Yi zi Ixxi Iyyi Izzi Ixyi Ixzi
where
mi = mass
xi, yi, zi = X, Y and Z-coordinates of the body centre of mass.
Ixxi , Iyyi , Izzi = mass moments of inertia about X, Y and Z-axes, respectively.
Ixyi, Ixzi, Iyzi= mass products of inertia about XY, XZ, and YZ-axes
respectively.
Z-directions,
respectively.
Mxi , Myi , Mzi = moments acting at point i about the X, Y and Z-axes,
respectively.
xi, yi, zi = X, Y, and Z-coordinates of point I.
3 DVIEW OUTPUT
The output of 3DVIEW includes:
The total mass and mass coordinates of the system mass centre (C.G.) in
the 3DVIEW coordinates.
The mass moments of inertia and product of inertia relative to DYNA6 axes.
Resultant excitation forces (including moments relative to DYNA6 axes).
The total mass, mass moments of inertia, and excitation force are input
directly into DYNA6 without any changes in magnitude or signs.
Example
Using the 3DVIEW subroutine, calculate the total mass, position of the
centre of gravity and mass moments of inertia for the whole foundation system
shown in Fig. 8.2. The system comprises four lumped masses, one body that can
be subdivided into rectangular blocks and one cylinder. The specific mass is 2500
kg.
The auxiliary 3DVIEW axes are placed at the base level as shown in Fig.
8.2. The body is subdivided in three rectangular blocks. Block 3 represents a hole
(depression) and is input as a rectangular block with negative mass.
Figure 8.2
DYNA6:
This is the main program, which calculates the stiffness and dumping
constants of the foundation as well as the system response to various load types.
Plot 2D:
This is the post processor. It allows the user to graphically depict the results
of DYNA6 on the screen and to print them.
FILE: By clicking on this option, a menu appears, listing the following file
options.
Import Dyna5 File: to import a data file generated by the older version
(DYNA5)
characters)
Settings: This option allows the user to input the settings including the units
of the problem, the output format and the damping safety factor
Foundation: This option allows the user to select the foundation type. There
are 6 foundation types available to choose from as described in Sec 2.1 listed here.
Pile: This option is for piled foundation. It allows the user to specify the footing
conditions (embedment, shape to input the piles, pile cap properties. It also
allows the user to input the properties of the soil along the pile (Pile-Soil), below
the pile tip (soil below) and adjacent to the pile cap for embedded footing (side
layers). Finally it allows the user to input the output options including the
coordinates for points where the resultant translations are required to be
calculated.
Half-Space: This option is for rigid foundation on a homogeneous soil (halfSpace). It allows the user to input the footing properties (Footing) or calculate
them (using 3Dview by clicking on Calculate). It also allows the user to input
the properties of the half-space (Below) and the adjacent embed footing (Side
Layers) Finally, It allows the user t input the coordinates of the points where the
resultant translations are required (Resultants).
Stratum Foundation: This option is for footing sitting on a layer under lain by
bedrock. It allows the user to input the properties (Footing) or calculate them
using 3Dview (by clicking on calculate) It allows the user to input the
properties of the soil layer (Stratum). Finally, it allows the user to input the
coordinate of the points where the resultant translations are required
(resultants).
allows the user to input the coordinates of the points where the resultant
translations are required (Resultants).
Rigid: This option is for rigid foundations (i.e. piers or caissons). It allows the
user to input the properties of rigid super structure (Structure). It also allows the
user to input the variation of the pier cross-section and soil along the pier and
below its tip (Elements), this includes the geometrical properties of the pier
(Rigid body elements). The properties of the soil along the shaft (Elements for
side Soil Layers) and the properties of soil below the tip of the pier (Elements
for soil Layers Below). And the properties of the soil layer below the foundation
to the bedrock (Below). Finally, it allows the user to input the coordinates of the
points where the resultant translations are required (Resultants)
Mat: This option is for flexible raft foundations. It allows the user to input the
mat data (Data) and the masses of the machines and their Pedestals (Masses)
and the soil properties below the foundation (Soil). It also allows the user to
view the mesh representing the foundation and the lumped masses locations
(Figure). After selecting one of the foundation options, the icons for Loads
and Run will be activated.
Loads: Allows the user to choose one from the following four potions
Transient: this option allows the user to input the characteristics of the
transient load (Transient) where the numbers of data points and time
interval as well as the force amplitudes are required. This option allows the
user to define the format of the FFT of the output load and response. Details
of transient loading are given in Section 3.1.
Random: This option allows the user to input the characteristics of a random
load (Random) including the number of data points observation time and
the minimum and maximum frequency required. Details of random loading
are given in Sec. 3.2.
Impact (Hammer): This option allows the user to input the properties of the
supporting system of the hammer block which could be a pad (Pad) for
Run: After the input process is complete, this option allows the user to run
Dyna6 (Execute Dyna6). After the run is completed successfully the user is
prompted to input the name of the output file. When the user inputs the name
of the output file, the control is returned to the main screen and the Post
Processor option is active.
Post Processor: This option activates the utility Plot 3D where the user can
depict the output graphically or in a spreadsheet. The output depends on the
problem and the format required. This output may include
Vertical Stiffness: This option allows the user to view and print the variation
of the vertical stiffness with frequency.
Vertical Damping: This option allows the user to view and print the variation
of the vertical damping with frequency.
Horizontal stiffness: This option allows the user to view and print the
variation of the vertical damping with frequency
Horizontal Damping: This option allows the user to view and print the
variation of the horizontal damping with frequency
Rocking Stiffness: This option allows the user to view and print the variation
of the Rocking stiffness with frequency.
Rocking Damping: This option allows the user to view and print the variation
of the rocking damping with frequency.
Coupled Hill Rocking Damping: This option allows the user to view and print
the variation of the Coupled Hill Rocking Damping with frequency.
Torsional Stiffness: This option allows the user to view and print the
variation of the torsional stiffness with frequency.
Torsional Damping: This option allows the user to view and print the
variation of the torsional damping with frequency.
Fourier Transform for Translations: It allows the user to view and print the
absolute value of the complex FFT of translation at C.G. (available for
transient loading option only).
Fourier transform for Rotation: It allows the user to view and print the
absolute value of the complex FFT of rotation at C.G. (available for transient
loading option only)
Fourier transform for Forces: It allows the user to view and print the absolute
value of the complex FFT of forces at C.G. (available for transient loading
option only)
Fourier Transform for Moments: It allows the user to view and print the
absolute value of the complex FFT of the moment, for positive frequencies
only (available for transient loading option only).
output
files
ASCII file, i.e., it may be print or edited using any editor. This can be done from the
file option and then using the 0
Windows file manager to find the file and open it
using NotePad and WordPad. .
3
9.4 Graphic Depiction of DYNA
V 6 Results
This option is invoked through the post processor menu or through
0
.
4
PostProcessor utility allows the used to view, edit and/or print the data and
graphs for the output
It has the following options:
File: This option allows the user to open a new file (New) open an existing chart
(Open Chart), save an existing chart (Save chart), import data for plotting from
Excel or Formula one files, or print preview and print data and /or chart.
Data: This option allows the user to change the format of data.
Chart: This option allows the user to edit the chart (title, footnote, legend, exit,
ect.)
View: This option allows the user to change the format for the Tool Bar.
Window: This option allows the user to change the format for the screen
(Cascade, File Horizontal, title, vertical and Arrange icons)
1/2
cannot work. Is
this correct? Negative diagonal stiffnesses may indeed occur under some
conditions; they are correct and cause no problems. They are most likely to be
returned for higher frequencies, rigid footings vibrating in the vertical direction
particularly with soil Poisson's ratio close to 0.5 (see Ref.2 and Fig. 2.5.1), heavy
single piles in weak soil or pile groups (Fig. 1.1.8). Negative stiffnesses pose no
problems in response calculations because the total soil resistance is a resultant
of its real part, depending on stiffness, and imaginary part, depending on damping.
At low frequencies, the diagonal stiffness constants are always positive.
Under large
(1)
comparisons were reported in Refs. 5, 6 and 32, for pile foundations in Refs. 32 to
35 and a few others. Some comparisons of DYNA with other theories are also
described in Chapter 2. More experimental validation is, of course, desirable but
in this regard the theories incorporated in DYNA are no less validated than most
other theories used in practice.
Pile (Single,
Battered)
Harmonic
(Quadratic)
Half-Space
Transient
Half-Space Harmonic
(Embedded/ (Quadrat
Weak)
Ic)
Pile
(Group)
Harmonic
(Quadratic)
Pile
(Group)
Harmonic
Half-Space
Random
Schematic
Stratum
Random
10
Pile
(Parabolic
Soil)
Harmonic
(Quadratic)
11
Rigid-Body
Transient
Response
12
CompositeMedium
Harmonic
(Quadratic)
13
Pile
(Group,
Battered,
Embedded)
Harmonic
(Quadratic)
Response
14
Matf
(Mat
Foundation)
Harmonic
(Quadratic)
15
Pile
(Battered,
General)
Harmonic
(Nonquadratic)
16
Pile
Capflexible
Harmonic
(Quadratic)
8.1
Example 1:
(a)
Evaluate stiffness and damping constants of a single fixed head, vertical pile with
constant cross-section, embedded in a layered soil medium comprising ten layers
adjacent to the pile; consider the soil properties and the weakened zone with
properties given in the Figure below.
(b)
Calculate the response of the pile to harmonic excitation which has constant
amplitudes equal to 1, 1, 1, 1, -1, 1 acting in all six directions and whose frequency
is = 100 rad/s. (The signs of forces and moments have to be consistent with the
sign convention shown in Figure 5.1.)
Units to use: lb, slug, ft.; mp = weak zone mass participation.
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. You will now be brought to a screen with several options, only one of which, you
can currently select. Click on Edit Settings under the Project Control header to proceed.
3. On the project settings menu, select the Imperial Units, check off Echo Input
Data and Stiffness/Damping Matrices and adjust the Damping Safety Factor to 2.
Click OK.
4. Next, select Edit Foundation from the Project properties menu and a selection
of foundations will appear in a new window. Select Pile and under Footing Shape click
Rectangle and OK
5. Now a window titled Pile Foundation will open up. To begin, select the Pile
Cap button to set the values for the pile cap.
6. Input the data for the footing and base dimensions and click OK.
7. Next select the Piles button and input the Pile Head Condition, Pile Material
Properties, Pile Tip Condition, and Pile Properties. Then select the option Individual
from the top of the menu and input the X and Y coordinates of the pile. Click OK twice
to return to the Pile Foundation menu.
8. Select the Pile-Soil button and proceed to input the data for Layer Depth and
Pile X-Radius. Then click Calculate geometrical data and the rest of the table should fill
out. Check off Weak Zone under Soil Medium Type and click on the tab at the bottom
of the Data for Pile-Soil System window, labeled Soil Elements to proceed.
9. Input all of the soil data and click OK to return to the Pile Foundation menu.
10. Select Soil Below on the Pile Foundation menu and input the requested
values. Click OK.
11. Select Output on the Pile Foundation menu and proceed to check off Unit
displacements and input the frequency labels and pile numbers then check off Actual
Displacements and do the same. Click OK three times to return to the main menu and
select Edit Load.
12. On the Load Type Choice menu, select Harmonic and click OK.
13. Select Non-Quadratic from the drop down menu and input all of the other
values accordingly. Click OK when finished to return to the main menu.
Evaluate stiffness and damping constants of a single pile battered in the xz - plane
whose angle to the vertical is = 10o. The pile is embedded in a layered soil medium
comprising ten layers adjacent to the pile;
(b)
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. You will now be brought to a screen with several options, only one of which, you
can currently select. Click on Edit Settings under the Project Control header to proceed.
3. On the project settings menu, select Imperial Units, check off Echo Input
Data and Stiffness/Damping Matrices. Click OK.
4. Next, select Edit Foundation from the Project properties menu and a selection
of foundations will appear in a new window. Select Pile and under Footing Shape click
Rectangle and OK
5. Now a window titled Pile Foundation will open up. To begin, select the Pile
Cap button to set the values for the pile cap.
6. Input the values for the Footing Base and Dimensions of the base and click
OK.
7. Next, select Piles and proceed to input the values for Pile Head/Tip Condition
and Pile/Material Properties. Then, click Individual and input the coordinates for the
single pile. Click OK.
8. On the same menu select Batter and input the batter angle Alpha. Click OK
twice to return to the Pile Foundation menu.
9. Now click the Pile-Soil button and check off Weak Zone and input the Layer
Depth and Pile X-Radius information and click the Calculate button to fill in the rest of
the table. Then click the Soil Elements tab at the bottom of the window.
10. Input all of the required data (Soil Layer Shear Wave Velocity, Soil Unit Wt.,
Poissons Ratio, etc.) Click OK to return to the previous menu.
11. On the Pile Foundation menu, select Soil Below and input the values in the
empty spaces. Click OK twice to return to the project properties menu.
12. Now select the Edit Load button to proceed to the Loading Type Choice menu.
Select the Harmonic load type and click OK to proceed to the next menu.
13. Input the values for the Frequency and Force and Moment Amplitudes. Click
OK.
14. Finally, select Run Project to process the inputs and receive the output values.
Click the Enter bar on the keyboard two times to view the Output File or select Output
File from the Project Properties menu.
The transient response of a 22 m high rigid silo (Figure 8.1) is calculated for
horizontal excitation due to the San Fernando Earthquake, 1971, component S90W with
peak acceleration of 0.11 g (Figure 8.2) acting in the x-direction.
The silo is supported by a 10.5 m diameter foundation resting on a half-space.
Use SI units: kN, Mg = 1000 kg, m.
Note: Better plots than those returned by DYNA can be obtained by generating the output
file DYNAPLT and using it together with the user's plotting subroutine. (See
example 12, p. 8.109.)
C
0
C
s
z
1
0
Half-space
VS = 152
m/sec
19.2
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click
Okay to continue.
2. You will now be brought to a screen with several options, only one of which, you
can currently select. Click on Edit Settings under the Project Control header to proceed.
3. On the Project Settings window check off Other. Under Frequency Units
select Hertz. Under Force Units, check off Other and input kN in the blank space.
Under Mass Units, check off Other and input Mg in the blank space. Finally, under
Length Units check off Meter and click OK.
4. Now select Edit Load on the Project Properties menu and the Foundation
Type Choice window will appear. On this window, select Half Space and under the
heading Footing Shape check off Circle. Click Ok.
6. Input the Footing Mass & Inertia, the Dimensions of Base, and Coord. Of Base
Center values. Click OK to return to the previous menu.
8. Now select the Edit Load option to proceed. The Loading Type Choice
window will appear, select Transient and click OK.
9. Now on the Transient Load Data menu, input the values under Data Points
Properties and check off ALL under Output of FFT for Load and Output of FFT for
Response. Under Amplitude of Forces click the Time History button.
10. The Time History window should be open. Input the adjacent values next to the
appropriate time in the empty space under the Value header. Click OK twice to return to
the Project Properties menu.
11. Finally, select the RUN PROJECT button to run the analysis. Click the Enter
key a few times until the analysis window closes. The Output File window should open up.
If it does not or the Output File was accidentally closed, click the Output File button under
the Project control label.
(a)
Evaluate the stiffness and damping constants of a symmetrical two mass hammer
foundation with an anvil pad. The foundation analyzed is described in detail in Ref.
23 and is shown in Figure 5. It is assumed that separation between the soil and
footing sides may occur due to heavy vibration of the hammer. Consequently, an
effective embedment depth of 5.2 ft (smaller than the actual embedment depth of
8.2 ft) is used in the analysis. The initial estimate of the first natural frequency of
the two mass system is 58 rad/s.
(b)
Calculate the response of this foundation to the impact of the hammer head assuming
initial velocity of the anvil equal to 1.52 ft/s.
(c)
Calculate the maximum stresses on the anvil pad and the soil.
s
PAD
(E = 216 x 104 lb/ft2, D = 0.1)
8
C
VS
f
(
5.2 ft
= 408 ft/s
(
1
VS =
500 ft/s
.10
3.12
.92
4.00)
.10
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should now be open. Click the Edit settings
option to proceed.
4. Select the Edit Foundation option to proceed. The Foundation Type Choice
window should open up. Check off Embedded under Embedment. Click OK.
6. Input the values in the empty spaces indicated under Footing Mass & Inertia
and Dimensions of Base. Click OK.
7. Select Side Layers on the Half-Space Foundation menu and the Side Layer
window should open up. Input the side layer properties in the allotted space and click
OK to proceed.
8. Select Below on the Half-Space Foundation menu and the Properties of Soil
Below window should open up. Input the values for the soil properties in the appropriate
spaces. Click OK twice to return to the Project Properties window.
9. On the Project Properties window, select Edit Load and the Loading Type
Choice window should appear. Select Hammer and click OK.
10. The Type of Hammer Load window should be open. Under Supporting
System, select Pad and under Hammer Load Type, select Short Duration. Click
OK.
11. The Pad System Data window should now be open. Input the values for
Hammer Foundation Properties, Short Duration Pulse data and Pad Supporting
System. Click OK.
12. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 5: Half-Space With Weak Zone and Harmonic Load (resultant translations
and stiffness tabulation are required)
(b)
(c)
The foundation is embedded in halfspace with Vs = 220 ft/s, = 100 lb/ft3, = .25,
D = 10% and is backfilled with weaker soil (side layer: depth = 2 ft, Vs = 180 ft/s, = 75
lb/ft3, = .25, D = .1); weak zone ratios: Gin/Gout = .25, in = .25, Din = 20%, t/Ro = 0.2,
M.P.F. = 0.25.
P
oint 1
Side
Layer
0 ft
VS
Weak Zone:
Gin/Gout = 0.25
.G
180ft/s
t/Ro = 0.2
= 75
D = 0.2
lb/ft3
= 0.25
D = 0.1
= 0.25
ft
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should now be open. Select Edit Settings
under the Project Control heading to continue.
3. In the Project Settings window, select Imperial Units under the Units
heading. Also, under the Frequency Units heading, select Hertz. Finally, under Output
Flags check off Echo Input Data and Output Draft Plots and click OK to proceed.
4. Select Edit Load on the Project Properties window and the Foundation Type
Choice window should open up. Now, select Half-Space, Embedded under
Embedment, and Rectangle under Footing Shape. Click OK to proceed.
6. The Footing Base Dimensions and Data window should open up. Fill out the
missing data beneath the headings Footing Mass & Inertia, Dimensions of Base, Coord.
Of Base Center and Coord. Of CG. Click OK to return to the Half-Space Foundation
menu.
8. Back on the Half-Space Foundation menu click Below to proceed. Fill in the
blank data and click OK.
10. Back on the Project Properties menu select Edit Load to open up the Loading
Type Choice window. Select Harmonic and click OK to proceed.
12. The Form Load Harmonic Phase Shift window should open up. Input the
Phase Shift data in the allotted space and click OK.
13. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Given:
The Machine
Total mass = 9600 kg
Exciting forces are due to rotor unbalances (quadratic) and act in vertical as well as
horizontal dire
2,
meer = 10 kgm2.
The height of the horizontal excitation = 3.66 which is also the height of the machine
centroid.
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window should open up. Under Units, check off SI
Units. Under Output Flags check off Echo Input Data and Output Draft Plots and click
OK to continue.
4. On the Project Properties window select Edit Load below the Project Control
heading. The Foundation Type Choice window will open up. Click Pile, check off
Rectangle blow the Footing shape header and click OK to continue.
5. The Pile Foundation menu should now be open. To continue, click the Pile Cap
button.
6. . The Footing Base Dimensions and Data window should open up. Fill out the
missing data beneath the headings Footing Mass & Inertia, Dimensions of Base, Coord.
Of Base Center and Coord. Of CG. Click OK to return to the Pile Foundation menu.
7. In the Pile Foundation menu click Piles to proceed. The Pile Head
Arrangements for Rigid Piles window will now be open. Check End-Bearing Tip below
the Pile Tip Condition header. Also, input the properties for the steel.
8. Next select Mesh in the tabs to open the Mesh Generation window. Check off
Rectangular Mesh. Input all of the remaining missing data. Click Apply to move on.
Finally, click OK to return to the Pile Foundation menu.
9. Click Pile-Soil to proceed. The Data for Pile-Soil System window will open up.
Fill out the blank windows with the correct data. Click Soil Elements at the bottom of the
window to continue.
10. In the window, fill in the missing data and click OK to proceed.
11.
Finally, on the Pile Foundation menu select Soil Below to open the
Properties of Soil Below the Foundation window. Fill out the blank data and click OK to
continue.
12. The Project Properties window should open up. Select Edit Load beneath the
Project Control heading. The Loading Type Choice window should open up. Select
Harmonic and click OK to proceed.
13. The Harmonic Load Data window should open up. Check off Output Phase
Angle. Also fill in the Frequency data, the Amplitude of Forces/Moments data, and check
off out of phase below the Applied Loads Phase heading. Finally, click Phase Shift to
proceed.
14. The Form Load Harmonic Phase Shift window should open up. Fill in the
Phase Shift data and click OK to finish.
15. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
(b)
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will now be open. Under Units, check off SI
Units. Also, check off Echo Input Data, Stiffness/Damping Matrices, and Output Draft
Plots under the Output Flags heading. Click OK. Click Edit Foundation under the
Project Settings header.
4. The Foundation Type Choice window will open. Click Pile and Rectangle
below the heading Footing Shape. Click OK.
5. The Pile Foundation window will open. First, click Pile Cap to proceed.
6. . The Footing Base Dimensions and Data window should open up. Fill out the
missing data beneath the headings Footing Mass & Inertia, Dimensions of Base, Coord.
Of Base Center and Coord. Of CG. Click OK to return to the Pile Foundation menu.
7. Back on the Pile Foundation window, select Piles to proceed. The Pile Head
Arrangements for Rigid Piles window will open up. Fill in the blank data and check off
End-Bearing Tip. Click Mesh to continue.
8. The Mesh Generation window will open up. Check off Rectangular Mesh and
X-direction below the header Generation Direction. Fill in the blank data spots and click
APPLY to proceed. Return to the Pile Foundation menu.
9. From the Pile Foundation menu, select Pile-Soil to open the Data for Pile-Soil
System window. Fill in the blank spaces on the Pile Elements data set and click Soil
Elements.
10. Fill in the Soil Elements data set and click OK to return to the Pile
Foundation window.
11. The Properties of Soil Below the Foundation window will open. Fill in all of
the missing blanks. Click OK to return to the Pile Foundation window. Click Output to
continue.
12. The Output Options for Pile Foundation window will open up. Check off Unit
Displacements, Actual Loads, Print in Output File and Print Distribution. Check off
OK to proceed.
13. The Loading Type Choice window will open. Select Harmonic and click OK
to proceed.
14. The Harmonic Load Data window will open up. Check off out of phase below
the Applied Loads Phase header and fill in the missing data beneath the Frequency and
Amplitude of Forces/Moments headers. Click Phase Shift to continue.
15. The Form Load Harmonic Phase Shift window should open up. Fill in the
Phase Shift data in the empty blanks and click OK to proceed.
16. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Using the random vibration approach, calculate the expected peak response of the
rigid silo shown in Figure 8.1 to horizontal earthquake excitation given the spectrum m2 S g
1 (Figure 8.5) applied in x-direction. The silo rests on half-space. (A smoothed spectrum
is preferable.)
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window should be open. Check off Other beneath the
Units header. Check off Other beneath the Mass Units and Force Units headers and
fill in the blanks. Also, check off Meter below the Length Units header. Finally, check of
Echo Input Data and Output Draft Plots below the Output Flags headers. Click OK.
6. The Footing Base Dimensions and Data window will open up. Fill in the
missing data and click OK to return to the Half-Space Foundation window.
7. Click Below to open the Properties of Soil Below the Foundation window. Fill
in the blank data and click OK till you reach the Project Properties Menu. Click Edit
Load below the Project Control header.
8. The Loading Type Choice window will open up. Click Random and OK to
proceed.
9. The Form Load Random window will open. Fill in all of the blank data
including the chart with the Frequencies and corresponding Values. Click OK to return
to the Project Properties menu.
10. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 9
Using the random vibration approach, calculate the expected peak response of the
silo shown in Figure 8.1 to the horizontal earthquake excitation given by Figure 8.5; in this
case the silo is supported by a stratum of limited depth of 20 m (Figure 8.6). (For shallow
strata nonzero damping should be used.)
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will open up. Check off Other beneath the
Units, Force Units and Mass Units headers. Fill in the empty blanks. Check off Echo
Input Data beneath the Output Flags header. Click OK to return to the Project
Properties window.
4. Under the Project Control header, click Edit Foundation. The Foundation
Type Choice window will open up. Click Stratum and check off Surface, Circle, and
Rigid. Click OK.
6. The Footing Base Dimensions and Data window will open. Fill in the blank
data spots and click OK to return to the Stratum Foundation window.
7. Click Stratum to open the Form Footing Stratum Data window. Fill in the
blanks with the appropriate data and click OK to return to the Project Properties window.
8. Under the Project Control header click Edit Load to open the Loading Type
Choice window. Click Random and OK to continue.
9. The Form Load Random window will open up. Input the missing data and fill
out the Frequencies and Values in the empty table. Click OK to return to the Project
Properties window.
10. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 10: One Pile Foundation With Parabolic Soil Shear Modulus
Distribution and Harmonic Load
(a)
(b)
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. Check off Imperial Units under the Units header. Also check off Echo Input
Data and Stiffness/Damping Matrices under the Output Flags header. Click OK to
return to the Project Properties window.
4. Under the Project Control header click Edit Foundation to proceed to the
Foundation Type Choice window. Click Pile and check off Surface, Circle, and Rigid.
Click OK.
5. The Pile Foundation window will open up. First click Pile Cap to continue.
6. The Footing Base Dimensions and Data window will open up. Fill in the blanks
and click OK to return to the Pile Foundation window. Click Piles to continue.
7. The Pile Head Arrangements for Rigid Piles window will open up. Fill in the
blank spaces and then click Individual to proceed.
8. The Coordinates for Individual Piles window will open up. Fill in the X-Coord
and and Y-Coord data and click OK to continue.
9. The Data for Pile Soil System window will open up. Check off Parabolic and
Weak Zone below the Soil Medium Type and Other below Calculating Pile
Geometrical Data. Fill in the data and click Soil Elements at the bottom of the window.
10. Fill in the missing data for the Soil Elements and click OK to return to the Pile
Foundation window.
11. Fill in the missing data. Click OK to return to the Project Properties window.
Click Soil Below to open up the Properties of Soil Below the Foundation window. Fill in
the data and click OK to return to the Project Properties window.
12. Below the Project Control header on the Project Properties window, click
Edit Load. The Loading Type Choice window should open up. Click Harmonic and OK
to continue.
13. The Harmonic Load Data window will open up. Fill in the data below all of
the headers. Click OK to return to the Project Properties window.
14. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will open up. Check off Other under the Units,
Force Units, and Mass Units headers and fill in the empty blanks with the correct units.
Ensure that Meter is checked off beneath Length Units. Below Output Flags, check off
Echo Input Data, Stiffness/Damping Matrices, and Output Draft Plots. Click OK
4. In the Project Properties window, click Edit Foundation below the Project
Control header to open the Foundation Type Choice window. Select Rigid foundation
type and ensure Surface and Circle are checked off. Click OK.
5. The Rigid Body Foundation window will open up. Click Structure to proceed.
6. The Structure and Foundation Data window will open up. Input the appropriate
data in the spaces and click OK to return to the Rigid Body Foundation window. Next,
click Rigid Body.
7. The Rigid Body Foundation Data window should open up. Input the rigid-body
data and check off Floating below the End Conditions header. Click OK to return to the
Rigid Body Foundation menu. Click Elements.
8. The Rigid Body Foundation/Soil Elements Data window should open up.
Check off Weak Zone below the Side Layers Properties header. Fill in the table and click
Elements for Side Soil Layers to proceed.
9. The Elements for Side Soil Layers section should be open. Fill in the table with
the appropriate data. Click Elements for Soil Layers Below.
10. The Elements for Soil Layers Below section should be open. Enter the data
into table and click OK to return to the Rigid Body Foundation menu. Now, click Below.
11. The Properties of Soil Below the Foundation window will open up. Fill in the
empty boxes and click OK to return to the Project Properties menu.
12. Click Edit Load below the Project Control header. The Loading Type
Choice window will open up. Click Transient and click OK to proceed.
13. The Transient Load Data window will open up. Input the number of Data
Points and time interval. Check off Graph beneath the headers labeled Output of FFT
for Load and Output of FFT for Response. Finally, click Time History beneath the
Amplitude of Forces header.
14. The Time History window will open up. Input the Values for each Time. Click
OK to return to the Project Properties menu.
15. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
a 26 ft thic
3
, Vs
, Vs
(The contents of output file DYNAPLT are used to produce the response
curves sown in Figure 8.10).
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window should open up. Check off Imperial Units
beneath the Units header. Also, check off Echo Input Data and Output Draft Plots.
Click OK to return to the Project Properties menu.
4. Click Edit Foundation and the Foundation Type Choice window will open up.
Click Composite and check off Rectangle below the Footing Shape header. Click OK.
6. The Footing Base Dimensions and Data window will open up. Input the
Footing data, Dimensions, and Coordinate data. Click OK to return to the CompositeMedium Foundation menu.
7. Click Below to open the Properties of Soil Layer Below Foundation. Input the
Layer data and click OK to return to the Project Properties menu.
8. Click Edit Load below the Project Control header to open the Loading Type
Choice window. Click Harmonic in the Loading Type Choice window and click OK.
9. The Harmonic Load Data window will open. Fill the blank boxes and check off
out of phase below the Applied Loads Phase header. Click Phase Shift.
10. The Form Load Harmonic Phase Shift window will open. Fill the boxes in
with the given data. Click OK to return to the Project Properties window.
11. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 13: Embedded Pile Foundation With Battered Piles and Harmonic
Load
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will open up. Check off SI Units below the
Units window. Also, check off Echo Input Data and Output Draft Plots below the
Output Flags header. Click OK to return to the Project Properties window.
4. In the Project Properties window click Edit Foundation below the Project
Control header. The Foundation Type Choice window will open up. Click Pile and check
off Embedded, Rectangle and Rigid. Click OK.
5. The Pile Foundation window will open up. Click Pile Cap to proceed.
6. The Footing Base Dimensions and Data window will open up. Input the
Footing, Dimensions, and Coordinates data. Click OK to return to the Pile Foundation
menu. Click Side Layers.
7. The Data for Cap/Footing Side Layers window will open up. Fill in the data
missing in the table and click OK to return to the Pile Foundation window. Click Piles.
8. The Pile Head Arrangements for Rigid Piles window will open up. Check off
End-Bearing Tip below the Pile Tip Condition header. Input the Pile Properties and
Pile Material Properties. Click Mesh and the Mesh Generation window will open up.
Check off Rectangular Mesh and X-direction below Generation Direction. Input all of
the missing data and click Apply. Then click Batter.
9. The Data for Battered Piles window will open up. Input the Alpha and Phi
values and click OK twice to return to the Pile Foundations menu. Click Pile-Soil to
proceed.
10. The Data for Pile-Soil system window will open up. Check off Circular Solid
below the Calculating Pile Geometrical Data header. Input the data in the table and click
Soil Elements at the bottom of the window.
11. The Soil Elements table will open up. Fill in the table and click OK to return
to the Pile Foundation window. Click Soil Below to continue.
12. The Properties of Soil Below the Foundation window will open up. Input the
Shear Wave Velocity, Unit Weight, Poissons Ratio, and Material Damping. Click OK
to return to the Project Properties window.
13. Click Edit Load below the Project Control header. The Loading Type
Choice window will open up. Click Harmonic and OK to proceed.
14. The Harmonic Load Data window will open up. Input the Frequency and
Amplitude of Forces/Moments data. Check off out of phase below the Applied Loads
Phase header. Click Phase Shift to proceed.
15. The Form Load Harmonic Phase Shift window will open up. Fill in the empty
boxes and click OK to return to the Project Properties window.
16. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 14:
3
s
an amplitude equal to 900. KN acting at the mat centre and a couple with an
amplitude of 6.0 kNm acting in a vertical plane passing through the mat centre and
parallel to the plate side. The mass of the pedestal and the machine resting on it are
represented by three lumped masses, m1, m2 and m3, shown in the Figure.
The Procedure:
a)
b)
The vertical couple is replaced by two equal but opposite vertical loads,
placed 1.0 m apart at nodes 24 and 26,
P = couple/distance = 6.0/1.0 = 6.0 KN.
c)
The mass of the mat is accounted for in the program while the keyword
MASS in the input enters the additional lumped masses only. (If there are
no additional lumped masses the input data should read MASS=0.)
d)
The keyword DYNFLEX places the soil stiffness matrix referring to the
nodes, 49 x 49 in this example, into a file of the same name for further use.
This file is not a part of the output file but it can be printed and is shown here
for illustration. (For a large number of nodes, this file can be very large and
impractical to print.)
If the printout of the stiffness matrix of the soil is desired it can be printed
within the output file using the command FLEXIBLE. (This is demonstrated in the
Example on p. 8.63.)
8.128
The elements of the soil stiffness matrix are, as in other cases,
K(I,J) = K1 + iK2
where K1 = real stiffness, and K2 = imaginary part of the stiffness yielding
the equivalent viscous damping constant
c(I,J) = K2
NOTES:
a)
In this example, both the mat and the finite elements used are square.
Rectangular mats with any aspect ratio (length to width ratio) can be
analyzed by DYNA3, however. The program would operate for any aspect
ratio for the finite elements but it is recommended that this ratio be kept
between 1.0 and 2.0 (0.5 and 1.0) for better accuracy.
b)
The MATF option may require a long computing time, especially when it is
run for many elements and frequencies on a PC. Example 14 takes 25 to 60
minutes on a PC depending on the system. Therefore, if the user wishes to
run this example just for a check, the number of frequencies may be reduced.
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will open up. Click SI Units below the Units
header. Check off Echo Input Data and Output Draft Plots below the Output Flags. Click
OK to return to the Project Properties menu. Click Edit Foundation below the Project
Control window.
4. The Foundation Type Choice window will open up. Click Mat and click OK.
5. The Mat Foundation window will open up. Click Data to proceed.
6. The Flexible Mat Data window will open up. Input the Mat Properties,
Response, and Number of Nodes data. Check off DYNAFLEX. Click OK to return to the
Mat Foundation window. Click Masses.
7. The Flexible Cap Lumped Masses window will open up. Input the Node and
Mass data. Click OK to return to the Mat Foundation menu. Click Soil.
8. The Properties of Soil Below the Foundation window will open up. Input the
Shear Wave Velocity, Unit Weight and Poissons Ratio. Click OK to return to the
Project Properties window. Click Edit Load below the Project Control header.
9. The Loading Type Choice window will open up. Click Harmonic and click
OK to proceed.
10. The Mat Harmonic Load Data window will open up. On the drop down menu,
select Non Quadratic. Input the Frequency data and check off the Output Displacement
Phase Shift option. Finally, click Loads.
11. The Vertical Harmonic Loads window will open up. Input the Node and Load
data and click OK to return to the Project Properties window.
12. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.
Example 15: Piles Battered in a General Plane with Printing of Group Stiffness
Matrix in Vertical and Lateral Directions
Analyze the stiffness and damping properties of a group of 6 steel piles used to
support an offshore structure. The piles are battered in a general plane (i.e. not contained in
a plane parallel to one of the two major planes, XZ and YZ). The required input properties
are given below and in Figure 24. Consider complete pile-soil separation for the top 0.5 m.
Print the vertical and coupled horizontal (lateral) group stiffness matrices referred
to the pile heads whose level is defined at 1.0 m above the mud line. The frequencies to
consider are 0.2 Hertz and 1.0 Hertz. The units to use are kN, Mg = 1000 kg, m.
Pile properties: Density = 77 kN/m3; Pile length (includes 1 m free length) = 31.0 m
Outside radius = 0.5 m, Inside radius = 0.46 m (A = 0.12 m2, I = 0.0139 m4)
Youngs modulus Ep = 2x108 kN/m2; Static load on pile = 4000 kN
Soil properties: Poissons ratio = 0.45; Damping, D=0.05. For soil below tips, Va =
300 m/s and = 21 kN/m3. Other properties are given in the table below.
1. First, open the program DYNA and click on the icon labeled New. A window
will pop up prompting you to enter the name of the project and specific case. Click Okay
to continue.
2. The Project Properties window should open up. Click Edit Settings under
Project Control to proceed.
3. The Project Settings window will open up. Check off Other beneath the
Units, Force Units, and Mass Units headers. Fill in the appropriate boxes with the
type of unit. Check off Hertz below the Frequency Unit header. Also, check off Echo
Input Data below the Output Flags window and click OK to return to the Project
Properties window.
4. Click Edit Foundation below the Project Control header. The Foundation
Type Choice window will open up. Click Pile. Check off Rectangle below the Footing
Shape header. Click OK to proceed.
5. The Pile Foundation window will open up. Click Pile Cap to continue.
6. The Footing Base Dimensions and Data window will open up. Input the Footing,
Dimensions, and Coordinates. Click OK to return to the Pile Foundation window. Click
Piles.
7. The Pile Head Arrangements for Rigid Piles window will open up. Fill in the
Pile and Material Properties. Click Individual.
8. The Coordinates for Individual Piles window will open up. Input the X-Coord
and Y-Coord data and click OK. Click Batter.
9. The Data for Battered Piles window will open up. Input the Alpha and Phi
values for the piles and click OK twice to return to the Pile Foundation window. Click
Pile-Soil.
10. The Data for Pile-Soil System window will open up. Check off Other under
Calculating Pile Geometrical Data and Weak Zone under Soil Medium Type. Fill in
the Pile Elements data in the table. Click Soil Elements.
11. The Soil Elements section should now be open. Input the required data to
complete the table. Click OK to return to the Pile Foundation window and click Soil
Below.
12. The Properties of Soil Below the Foundation window is now open. Input the
required data. Click OK to return to the Pile Foundation window. Click Output.
13. The Output Options for Pile Foundation window should be open. Check off
Print in Output file below both the Piles Vertical Stiffness header and Piles Lateral
Stiffness. Click OK to return to the Project Properties window.
14. Click Edit Load below the Project Control header. The Loading Type
Choice window should open up. Click Harmonic and click OK to proceed.
15. The Harmonic Load Data window should open up. In the drop down menu,
select Non-Quadratic. Input the data for Frequency, Amplitude of Forces/Moments, and
make sure in phase is checked off below the Applied Loads Phase header. Click OK to
return to the Project Properties window.
16. Finally, click RUN PROJECT to run the analysis. Click the Enter button
several times to advance the windows. The Output File should open up, but if it does not,
select Output File from beneath the Project Control heading.