Natural Periods of an Elastic Column
Note
To view this project in FLAC3D, use the menu command . Choose “Dynamic/ NaturalPeriodsElasticColumn” and select “NaturalPeriodsElasticColumn.f3dprj” to load. The project’s main data files are shown at the end of this example.
A column of elastic material resting on a rigid base has natural periods of vibration, depending on the mode of oscillation and the confining conditions. Three cases are examined: an unconfined column; a confined column in compression; and a column in shear.
The column is loaded by applying gravity in either the \(x\)- or \(y\)-direction and observing the oscillations with zero damping. The case of confined compression is modeled by inhibiting lateral displacement along the vertical boundaries, which prevents lateral deformation of the mesh. For unconfined compression, lateral displacement is not inhibited. For the column in shear, vertical motion is inhibited to eliminate bending modes; the loading is applied laterally.
The theoretical value for the natural period of oscillation, T, is given by Equation (1),
where \(E^*\) is the appropriate modulus selected from Table 1.
Confined Compression | Unconfined Compression | Shear |
---|---|---|
\(K + (4/3)\ G\) | \({{9KG} \over {3K + G}}\) | \(G\) |
(Young’s modulus) | ||
2.5714 × 104 | 1.2 × 104 | 1.0 × 104 |
Material properties are given in Table 2.
Properties | Symbol | Value | Comment |
---|---|---|---|
bulk modulus | K | 2.0 × 104 | for compression tests |
shear modulus | G | 0.428562 × 104 | |
Poisson’s ratio | \(\nu\) | 0.4 | |
bulk modulus | K | 1.0 × 104 | for shear tests |
shear modulus | G | 1.0 × 104 | |
Poisson’s ratio | \(\nu\) | 0.125 | |
density | \(\rho\) | 1.0 | |
applied gravity | \(g_z\) | -10.0 | for compression tests |
\(g_x\) | 0.1 | for shear tests | |
column height | L | 800 | |
column width | W | 100 |
The theoretical periods and calculated (FLAC3D) natural periods of oscillation averaged over several periods by the FISH function average-period are compared in Table 3.
Confined Compression | Unconfined Compression | Shear | |
---|---|---|---|
Theoretical | 19.96 | 29.21 | 32.00 |
FLAC3D | 19.96 | 29.26 | 31.98 |
Data Files
ConfinedCompression.f3dat
;-----------------------------------------------------------------------
; Script file for dynamic problem 'Natural periods of an Elastic Column'
; case 1) Confined compression
;-----------------------------------------------------------------------
model new
model largestrain off
model title "Natural periods of an Elastic Column - Confined Compression"
model config dynamic
;; Create zones
zone create brick point 0 (0,0,0) point 1 (100,0,0) ...
point 2 (0,100,0) point 3 (0,0,800) ...
size 1,1,8
;; Assign material model and properties
zone cmodel assign elastic
zone property bulk 2e4 shear 0.428562e4 density 1
;; Boundary Conditions
zone gridpoint fix velocity-x
zone gridpoint fix velocity-y
zone gridpoint fix velocity-z range position-z 0
;; Remaining setup
model gravity 0 0 -10
history interval=1
zone dynamic damping local 0
zone history velocity-z position (0,0,800)
model history dynamic time-total
;; Solve to 200 seconds
model solve time-total=200
model save 'ConfinedCompression'
UnconfinedCompression.f3dat
;-----------------------------------------------------------------------
; Script file for dynamic problem 'Natural periods of an Elastic Column'
; case 2) unconfined compression
;-----------------------------------------------------------------------
model new
model largestrain off
model title "Natural periods of an Elastic Column - Unconfined Compression"
model config dynamic
;; Create zones
zone create brick point 0(0,0,0) point 1(100,0,0) ...
point 2(0,100,0) point 3(0,0,800) ...
size 1,1,8
;; Assign model and properties
zone cmodel assign elastic
zone property bulk 2e4 shear 0.428562e4 density 1
;; Boundary Conditions
zone gridpoint fix velocity-z range position-z -1,1
;; Remaining setup
model gravity 0 0 -10
history interval=1
zone dynamic damping local 0
zone history velocity-z position (0,0,800)
model history dynamic time-total
;; Solve to 200 seconds
model solve time-total=200
model save 'UnconfinedCompression'
Shear.f3dat
;-----------------------------------------------------------------------
; Script file for dynamic problem 'Natural periods of an Elastic Column'
; case 3) Shear
;-----------------------------------------------------------------------
model new
model largestrain off
model title "Natural periods of an Elastic Column -- Shear"
model config dynamic
;; Create zones
zone create brick point 0(0,0,0) point 1(100,0,0) ...
point 2(0,100,0) point 3(0,0,800) ...
size 1,1,8
;; Assign model and properties
zone cmodel assign elastic
zone property bulk 1e4 shear 1e4 density 1
;; Boundary Conditions
zone gridpoint fix velocity range position-z -1,1
zone gridpoint fix velocity-z
;; Remaining setup
model gravity 0.1 0 0
history interval=1
zone dynamic damping local 0
zone history velocity-x position (0,0,800)
model history dynamic time-total
;; Solve to 200 seconds
model solve time-total=200
model save 'Shear'
⇐ Verification Problems | Comparison of FLAC3D to SHAKE for a Layered, Linear-Elastic Soil Deposit ⇒
Was this helpful? ... | PFC 6.0 © 2019, Itasca | Updated: Nov 19, 2021 |