FB-MultiPier Help Manual PDF
FB-MultiPier Help Manual PDF
FB-MultiPier Features
entire bridge. The 2D Bridge window shows the bridge layout in plan and the 3D Bridge window shows the 3D
visualization of the bridge.
Pier rotation
Each pier can have a rotation about the vertical (z) axis. This is ideal for modeling skew bridges and radial piers
on curved alignments.
Bridge superstructure
The bridge superstructure is incorporated into the model using an equivalent beam that connects the centerline
of two piers. The bearing connections at the pier supports can be released, constrained, or user-defined using a
custom load-displacement curve.
Two rows of bearing locations
Two independent lines of bearings accommodate the transfer of load from the bridge superstructure to the piers.
Because the bearings are offset from the center of the pier cap, any pier cap torque induced from unequal spans
is automatically included.
Wind Load Generation
Wind loads can be automatically generated for the bridge superstructure. The resulting loads are applied to the
bearings at each pier, more specifically to the bearing locations of the transfer beam. These are based upon
tributary areas and are redistributed during analysis if the superstructure is continuous.
Dynamic Pier Analysis
Time step integration
Time history load functions and ground acceleration records can be applied to the model. Different time step
integration methods are available as well as a variety of analysis control parameters. Concentrated masses and
dampers can be added to the model to simulate added mass and energy dissipation effects.
Modal analysis
The modal analysis option performs a frequency analysis of the model. Both frequencies and mode shapes are
provided as output results.
Dynamic soil modeling
Soil gap modeling is available to model energy dissipation due to hysteretic damping. Cyclic degradation
parameters are also available to modify the lateral soil response during dynamic loading.
Animated results
The 3D model displacement results can be animated for a time step integration analysis. Animation results can
be played and paused and a slider bar is provided for selectively viewing individual time step results.
Time-Displacement plots
The displacement results for any model node can be plotted over time.
Seismic database
Ground acceleration records and response spectrums are provided for notable earthquakes.
1 Program Menus
1. File Menu
2. View Menu
3. Control Menu
4. Wizard Menu
5. Help Menu
6. About FB-MultiPier
l
l
l
l
l
l
l
Version number
Expiration date
License type
Number of seats
License file path where the license file is located
Physical RAM
Drive Space
2 Model Data
1. New Project/Problem Page
2. Analysis Page
3. AASHTOPage
4. DynamicsPage
5. Pushover Page
6. Pile Cap Page
7. Pile Page
8. Soil Page
9. Pier Page
10. Wall Structure Page
11. Extra Member Page
12. Load Page
13. Spring Page
14. Concentrated Mass/Damper Page
15. Retained Soil Page
16. Bridge Page
Select the unit type (English or Metric) in the "Select a New Problem Type" window.
Note: With this option, the Pier page becomes the Wall Structure page.
For a complete list of problem options go to the Problem Page.
Note: With this option, the Pier page becomes the Wall Structure page.
For a complete list of problem options go to the Problem Page.
Pile/Pier Behavior
2.
Section Properties
3.
Iteration Control
4.
Analysis Type
5.
Soil Behavior
6.
7.
Design Options
8.
Stiffness Options
9.
Print Control
Deflections do not cause secondary moments; no P-delta moments (moments of the axial force times the
displacements of one end of element to another).
Nonlinear Behavior:
Uses input or default stress strain curves which are integrated over the cross-section of the piles/pier
components.
Uses P-delta moments (moments of the axial force times the displacements of one end of element to
another). Furthermore, since the user subdivides the pile into a number of sub-elements, the P-y moments
(moments of axial force times internal displacements within members due to bending) are also modeled.
Return to the Analysis Page.
Static
2.
Dynamic
Check "User-defined phi" to enter a custom phi factor, or leave the option unchecked if you want to use the
default value.
Return to the Analysis Page.
1.
Pile Displacements
2.
3.
Pile Properties
4.
Pier Displacements
5.
6.
Pier Properties
7.
8.
9.
10.
11.
Unbalanced Forces
12.
13.
14.
15.
16.
17.
18.
Interaction Data
19.
Coordinates
20.
XML Data Printing Creates XML output file that can be used to extract FB-MultiPier data. See FBMultiPier XML Specification documentation.
Note: The AASHTO Combination option in the Design Options section of the Analysis Page must be selected
for this Page to appear.
2.
3.
4.
5.
6.
7.
8.
9.
10.
The AASHTO Load Case Manager manages the type and number of load cases in your model. Changes made
with this manager apply to every pier.
Note: Certain load case types are grouped together. Example, "Wind on Structure" and "Wind on Live
Load". Changing the number of cases for one of these types will automatically change the number of cases for
the other type.
Return to the AASHTO Page.
Figure: 2.3.e AASHTO Load Fatctor Table using Load Factored Design (LFD)
AASHTO Load Fatctor Table using Load Factored Design (LFD). When using LFD, Gamma (load factor) is
multiplied by each coefficient in the corresponding group. These products can be seen in the AASHTO Load
Preview Combination dialog. For example, Dead Load (D) has a default value of 1.00. It is multiplied by the the
GROUP-I Gamma value 1.3. Thus, the AASHTO Load Preview Combination dialog displays the value 1.3 (1.00 *
1.3 = 1.3); This is the value used in the analysis.
Return to the AASHTO Page.
Shows the load combination that will be run. Color changes indicate limit states.
Table 2
Table 3
From AASHTO LRFD Bridge Specifications, 2010, by the American Association of State Highway and
Transportation Officials, Washington, D.C. Used by permission.
Return to the AASHTO Page.
Section 3.8.2 of the AASHTO Bridge Design Specification requires a vertical wind force to be considered in Limit
States Strength III and Service IV.
To facilitate this requirement, beginning with FB-MultiPier Version 4.16, the FB-MultiPier program enables users
to include two separate load cases, i.e., VP (vertical pressure) and VR (vertical pressure for a reversed load
case), only for Strength III and Service IV.
Note: VR appears if and only if Reversible Loads is selected in the Limit States to Check found on the
AASHTO page.
2.3.10 Updated Load Factors from FB-MultiPier Version 4.15 and Earlier
Versions
The default load factors of DD, EH, EV, ES, CR, SH, TG and TU have been updated since FB-MultiPier Version
4.15. Detailed description of changes are summarized in the table.
Analysis Type
Two dynamic analysis types are available.
1.
Time Step Integration - Uses implicit integration to solve for results at every time step.
2.
Modal Response Analysis Applies static loads and then performs a response spectrum analysis using
the equilibrium (deformed) position. Performs a CQC of the modal analysis results.
Modal analysis requires the user to select the number of modes to use in the analysis. For a Modal
Response Spectrum analysis, enter the number of modes to consider and the damping ratio used for the
response spectrum. Then, check the modal contribution factors in the printed output file to ensure that at
least 90% of the structure mass is participating in the response.
2.
3.
Wilson-Theta.
Enter a constant value for the time step and the number of time steps to consider in the analysis.
Dynamic Relaxation
The reported analysis results do not include the effect of static loads (i.e. self weight). Adding the static results
and response spectrum results may not be conservative and is left to engineering judgment.
Dynamic Relaxation feature is utilized to initialize the system in equilibrium with permanent (static) loading, while
minimizing non-physical system vibrations due to the presence of permanent loads.
Global Damping
Two types of damping input are available.
1.
For Time Step Integration type, Rayleigh damping can be specified. The damping is proportional to the
mass and stiffness. Factors can be entered for the pier, piles, and soil. Note that if the Damping
checkbox is checked, but the mass and stiffness damping coefficients for the Soil are both defined as
0.0, then radiation damping will be applied to embedded nodes in the model. See "Dynamic Soil
Properties" for additional details about radiation damping
2.
For Modal Response Type, modal analysis damping can be specified. Enter the damping ratio.
function. Rows can be inserted or deleted as needed. The "Update Table" button sorts the values according to
increasing time. You can drag and drop a range of data points from a spreadsheet directly into the table.
Click on the Run Pushover Analysis checkbox to activate the pushover analysis module.
There must be 2 load cases. The first load case is used to apply permanent loads that will not be incremented
(i.e. self weight). The second load case is used to specify the load that will be incremented.
Enter the number of pushover steps and the load increment factor. The load increment factor multiplies the loads
in the second load case to create an accumulating load that is applied until convergence cannot be achieved.
For example, a load increment factor of 1.0 would add 100% of the original load to each incremental load case. If
the original load increment was 10 kips, the second load increment would be 20 kip load, the third increment 30
kips, and so forth for the number of load steps. The failure load is printed to the output file when a load is reached
that can not converge to a solution.
Check the checkbox to activate the pile cap "Bearing Resistance" feature
Check the checkbox to activate the pile cap Bearing Resistance feature. When activated, FB-MultiPier assigns
vertical soil springs to each of the nodes in the pile cap. The Bearing Resistance feature can only be utilized if
the mid-surface of the pile cap is at or below the elevation at the top of soil layer 1 (i.e. for submerged pile cap).
Note that any single soil set number can be applied to the cap elements for the purposes of using the Bearing
Resistance feature.
Based on bearing pressure vs displacement curves for each vertical soil spring, vertical reaction forces are
calculated for a given vertical displacement. Note that vertical springs are active at all cap nodes that do not
correspond to pile head nodes. In particular, this process is driven by the current vertical pressure acting at a
given pile cap node, and the corresponding tributary bearing area associated with that pile cap node:
Eqn: 2.6.A
where,
F is the vertical pressure at a pile cap node for a given vertical displacement
Z is the vertical displacement inducing the pressure at the pile cap node
K is the subgrade modulus of soil at the mid-surface of pile cap
Qu is the ultimate bearing resistance of the soil at the mid-surface of pile cap
Note that the vertical pile cap bearing springs only act in compression.
Subgrade modulus (K) provides the initial slope of the bearing pressure versus displacement curve, whereas the
ultimate bearing resistance (Qu ) provides the curve asymptote. Subgrade modulus (K) and ultimate bearing
resistance (Qu ) values are calculated as described below.For Sand:
a) Subgrade modulus (K)
For sand, subgrade modulus is directly input from within the Soil Page.
b) Bearing capacity (Qu ) The ultimate bearing capacity for pile caps in Sand ( 0 and c = 0) is calculated by the general bearing
capacity equation, Eqn: 2.6.B , by Meyerhof (1963)
Eqn: 2.6.B
For zero cohesion (c = 0) Qu becomes:
Eqn: 2.6.C
where,
Nc, Nq and N are bearing capacity factors (Table 1)
Sc, Sq and S are shape factors (Table 1)
Dc, Dq and D are depth factors (Table 1)
Table 1
Factors ( > 0)
Values
Nc
Nq
N
Sc
Sq and S
Dc
Dq and D
is angle of internal friction (rad), and is equal to zero for cohesive soils.
is moist unit weight of sand
v is effective vertical stress at the mid-surface of the pile cap.
B is width of pile cap
L is length of pile cap
z is vertical distance of mid-surface of pile cap from the ground surface.
For Clay:
a) Subgrade modulus (K)
For clay, subgrade modulus is selected from the table (Reese, 1984) below, based on the directly input value of
undrained shear strength (Cu):
Cohesive Soil
Properties
Shear Strength
Unit
Soft
Medium
Stiff
Very Stiff
Hard
(Cu,
250-500
500-1000
1000-
2000-
4000-
psf)
Subgrade Modulus
(Static Loading)
Subgrade Modulus
(cyclic Loading)
2000
4000
8000
(k, pci)
30
100
500
1000
2000
(k, pci)
200
400
800
b) Bearing capacity (Qu ) The ultimate bearing capacity for pile caps in Clay ( = 0 and c 0) is calculated by the general bearing capacity
equation, Eqn: 2.6.D , by Meyerhof (1963)
Eqn: 2.6.D
For zero internal friction angle ( = 0) bearing factor N is zero. Accordingly, Qu becomes:
Eqn: 2.6.E
where,
Nc, Nq and N are bearing capacity factors (Table 2)
Sc, Sq and S are shape factors (Table 2)
Dc, Dq and D are depth factors (Table 2)
Table 2
Factors ( = 0)
Values
Nc
Nc = 5.14
Nq
Nq = 1
N = 0
Sc
Sc = 1 + 0.2(B/L)
Sq and S
Sq = S = 1
Dc
Dc = 1 + 0.2(z/B)
Dq and D
Dq = D = 1
A different edit window appears depending upon the type of section description selected.
If "Gross Properties" is selected then the Gross Pile Properties window will appear if you click "Edit" from the "Edit
Cross Sections" in Section Description. Note that gross section properties can only be specified for linear pile
behavior.
Otherwise, if "Full Cross Section" is selected, then the Full Cross-Section Pile Properties window will appear.
Return to the Pile Page.
Use Pile/Shaft Segment List to Add and remove pile/shaft segments of the selected pile type.
Database Section Selection
If the "Use Database Section" option is selected, the user can select from a predefined set of cross-sections.
In the Gross Pile Properties page, there is only one option (Linear Pile) when you click on the "Retrieve Section"
button.
If the "Customize Current Section" option is selected, the user can customize the current cross section.
Furthermore, the user can also save custom cross sections by clicking the "Add To Database" button. The user
can also delete the Database section by clicking the "Delete Section" button
Section Type
Select a cross section type from the following:
1.Circular Pile
2.Square Pile
3.H-Pile
Note: this option is only available if the "Customize Current Section" option is selected.
Segment Dimensions
Enter the following data for the dimensions of the segment:
1.Length
2.Area
3.DiameterOnly available for a circular pile
4.WidthOnly available for a square pile
5.DepthOnly available for a square pile
6.[Unit] Weight
Note: Except Length, Other options are only available if the "Customize Current Section" option is selected.
Section Properties
Enter the following data for the dimensions of the segment:
1.Inertia 2 AxisThe moment of inertia about the 2-axis
2.Inertia 3 AxisThe moment of inertia about the 3-axis
3.Torsional Inertia
4.Youngs Modulus
5.Shear Modulus
Note: this option is only available if the "Customize Current Section" option is selected.
Return to Pile Section Description page.
1. Length
2. Diameter (d)
3. Unit Weight
Rectangular Section:
1. Length
2. Width (w)
3. Depth (h)
4. Unit Weight
H-Pile:
1. Length
2. Unit Weight
Pipe Pile / Pipe Pile (Concrete Filled):
1. Length
2. Diameter (d)
3. Thickness
4. Unit Weight
Detailed Cross Section
By clicking the "Section Details" button, one can edit the bar groups and material properties of the cross section
in a spreadsheet format. Select a segment from the "Section List" and a pile set from the "Pile Type" list to edit.
Choose between Mild steel or Prestress for the type of steel in the bar groups. If Prestressed is chosen, then
enter the prestress after losses.
For the type of Shear Reinforcement, the user can select either spiral or tied.
Casing
Casing can be added by checking "Casing" Checkbox. Click "Edit Casing" to enter the Shell Thickness.
casing diameter. Those portions of the cross section that are used for demand and capacity calculations are
displayed below. Refer to the "Confined Concrete Model" section for additional details concerning the
constitutive relationships used for each cross section component.
Choose between Mild steel or Prestress for the type of steel in the bar groups. If Prestressed is chosen, then
enter the prestress after losses.
Hollow Section
Hollow Section can be added by checking "Hollow Section" Checkbox. Click "Edit Void" to enter the diameter for
a circular void, or the length and width for a rectangular void of the Hollow Section.
H-Pile Properties
Enter the depth, width, web thickness, and flange thickness of the H-pile in the text boxes.
Section Orientation
Select the orientation of the H-pile (Web horizontal or web vertical).
Return to the Full Cross Section Pile Properties page.
Enter a section length, diameter, shell thickness, and unit weight in the "Section Dimension" Panel.
Concrete is not included in this cross section. (fc and Ec are set to zero.)
Return to the Full Cross Section Pile Properties page.
Enter a section length, diameter, shell thickness, and unit weight in the "Section Dimensions" Panel.
Concrete is included in this cross section. To set the section constitutive properties for concrete and steel, click
the 'Edit Properties' button.
For more information on inputting data for the pipe pile (concrete filled) section, see the "circular section" type.
Return to the Full Cross Section Pile Properties page.
Note: Linear interpolation is used to calculate stresses for strain levels beyond those provided in the stressstrain curve.
Return to the Full Cross-Section Pile Properties or Full Cross Section Pier Properties page.
3. Torsional
A. Hyperbolic
B. Custom T-
4. Tip
A. Driven Pile
B. Drilled Shaft Sand
C. Driven Pile Sand (API)
D. Drilled Shaft Clay
E. Driven Pile Clay (API)
F. Drilled Shaft IGM
G. Custom Q-Z
For a selected soil, one can enter/edit the properties by selecting one of the option(Lateral, Axial, Torsional or
Tip) and then clicking the button explained below:
The "Edit" button will activate "Additional Soil Properties" dialog. Clicking the "Dynamic Properties " button will
open the Soil Dynamics Dialog which will allow the user to input additional soil properties that pertain only to
dynamic type analysis.
The"Plot" button will activate the "Printable Soil Graph" dialog which allows the user to view and print the various
soil curves (load vs. deflection graph) based upon the selected options (Lateral, Axial, Torsional or Tip).
The "Group" button will activate the "Advanced Soil Data" dialog to specify for the soil model. The user can
specify the P-multipliers.
The "Table" button will activate the "Soil Table" dialog for entering soil properties. This feature allows the user to
view/modify multiple soil sets and layers and the same time and quickly enter properties for each. See "Soil-Pile
Interaction" for details on the soil properties.
When using a custom soil curve, one can enter/edit the properties by selecting one of the option(Lateral, Axial,
Torsional or Tip) and then clicking the "Edit" button.
The "Specify Top and Bottom Layer Props" checkbox allows you to enter different soil properties at the top and
bottom of each layer. The values will be interpolated across the layer.
2.8.5 Elevations
Enter the elevation of the water table, and the top and bottom elevation of each layer.
Note: Changes made to soil layer elevations AFTER extra members have been applied to the model may
cause extra members to shift position. The explanation 'Extra Member Positioning with Respect to Model
Changes' describes this further.
User can rotate the Pier and/or Bearing Pads on the Bridge Page. The rotation of the bearing is independent of
that of the Pier rotation and is only present for the Bridge Models.
Modify the properties of a linear pier cross section in the following fields:
Pier Component
Select the pier component to edit, or add and remove components.
2. Depth (h)
3. Unit Weight
H-Section:
1. Unit Weight
Bullet Section:
1. Diameter (d)
2. Width (w)
3. Unit Weight
Parabolic Taper Cantilever Properties
Enter the depths for a cantilever with a parabolic taper.
Casing
Casing can be added by checking "Casing" Checkbox. Click "Edit Casing" to enter the Shell Thickness.
casing diameter. Those portions of the cross section that are used for demand and capacity calculations are
displayed below. Refer to the "Confined Concrete Model" section for additional details concerning the
constitutive relationships used for each cross section component.
Choose between Mild steel or Prestress for the type of steel in the bar groups. If Prestressed is chosen, then
enter the prestress after losses.
Hollow Section
Hollow Section can be added by checking "Hollow Section" Checkbox. Click "Edit Void" to enter the diameter for
a circular void, or the length and width for a rectangular void of the Hollow Section.
H-Section Properties
Enter the depth, width, web thickness, and flange thickness of the H-Section in the text boxes.
Section Orientation
Select the orientation of the H-Section (Web horizontal or web vertical).
Return to the Full Cross Section Pier Properties page.
First choose a material type on the left, and then enter the properties for that material in the text boxes.
Note: Linear interpolation is used to calculate stresses for strain levels beyond those provided in the stressstrain curve.
Return to the Full Cross-Section Pile Properties or Full Cross Section Pier Properties page.
3.To taper a cross section, check the appropriate 'Taper' checkbox. This will create another column in the table,
so that cross section data can be entered for each end of the pier component. Example: column bottom, column
top.
4. To create custom Stress/Strain Curves, select "Custom Stress/Strain" in the 'Material Properties' field. Then
click the 'Custom Curves' button in the 'Custom Curves' field.
5.For quick access to directions on using the table, click the "Help >>"button.
6.When printing the table, to make the table more easily fit on a single page, hide the 'Help' section on the right
side of the table, by clicking the 'Help' button, so that the arrows point to the right (Help >>).
Note: Extra members are not available in the following models: Pile and Cap Only, Stiffness, Single Pile,
Column Analysis.
1. Load Case
2. Buoyancy
3. Applied Nodal Loads
4. Loads
5. Load Table
6. AASHTOLoad Table
7. Distributed Load
8. Vessel Collision
Check "Prescribed Displacement" to apply a displacement (rather than a load) to a node. Loads and
Displacements can not be applied at the same node in the same load case.
In AASHTO load mode, the Load Page looks as follows:
2.12.2 Buoyancy
The buoyant force on the bridge substructure that is submerged, i.e., below the water table, is automatically
computed if a buoyancy factor greater than 0 is selected in non-AASHTO mode or if buoyancy is activated
(checked on) in AASHTO mode. The computation includes piles, pile cap, pier columns. Partial buoyancy of the
pile cap is accounted for, where the volume of the pile cap that is submerged will be used to determine the
buoyant force acting on the pile cap. A convenient way to check buoyancy and self-weight calculations is to
include only these loads, run the program, and then view the "Sum of Total Soil Spring Loads", Z direction in the
output file.
Return to the Load Page.
2.12.4 Loads
Select whether or not pre-loading conditions (i.e. thermal stresses, construction loads, shoring, etc.) are present.
For the pre-loading situation, the equilibrium loads are found from the pre-loading. Then, after equilibrium is
established, the analysis uses the equilibrium conditions to calculate the solution for the load cases.
Note that preloading can be applied (or not) for any load case of interest.
Enter point loads in the Xp, Yp, and Zp directions, and moments about the Xp, Yp, and Zp axes.
Also, enter factors for self-weight and buoyancy (for non-AASHTO loads). Note that self-weight corrections are
made for overlaps between:
1) Pile portions that lie within the physical pile cap thickness;
2) Pier column portions that lie within the physical pile cap thickness; and,
3) Pier column portions that lie within the physical pier cap thickness.
Check Applied Displacement to specify a displacement rather than a load for a node.
Return to the Load Page.
The "Load Table" is used to define nodal loads in a spreadsheet-style format. Static and dynamic loads are
separated into two separate tables that can be toggled using the "Table Format" options.
Following are the Load Tables for Static and Dynamic analysis types.
Static Loads
Enter the load case, node, and load values.
For a selected Load Case, the Distributed Load dialog helps to convert constant, uniform, and linear load
distribution types along an element length into nodal loads. Figure: 2.12.i shows the Distributed Load dialog.
The 3D View window (right side in the main program window) displays the nodes for the selected structural
member.
The user can also highlight an individual pile or pier column of interest using the Highlight Member drop-down
menu. For example, when the user selects Pile 1 from the Highlight Member drop-down menu, Pile 1 in the
3D View window will then be highlighted. This feature facilitates location of a given pile or pier column of
interest.
Select one of the following load distribution options from the Load Types drop-down menu.
1. Constant
2. Uniform
3. Linear
The Constant load distribution option generates identical loads. The Uniform load distribution option
generates loading based on tributary length. The Linear load distribution option generates loading based on
the interpolation between the starting and ending nodes load magnitudes and tributary length.
In the table, the starting and ending node numbers can be input for the element lengths over which the
distributed load is to be applied. The node numbers can be directly input in the table, or the node numbers can
be graphically selected by clicking nodes in the 3D View window. After selecting the nodes over which the
distributed load will act, input the Xp, Yp, and Zp components of the distributed load or moment in the table.
Note: Nodes can only be selected from among the nodes displayed in the 3D View window, and the
selected nodes should be located along a straight line. For example, if pile nodes are selected, both selected
nodes must be on the same pile; if column nodes are selected, both selected nodes must be on the same
column; if pile cap nodes are selected, both nodes must be in the same pile cap grid row or grid column.
The Clear Load Generation Table button will clear all contents from the table.
Clicking the Generate Load On The Selected Member button will generate nodal loading on selected nodes
using the distribution option, load magnitude entered, and the tributary length.
The Generate Load On Other Members button is active only for Pile and Column members. Clicking the
Generate Load On Other Members button for the selected pile member will generate equivalent loads on all
Other Member piles having the same pile type, batter conditions and soil set. Similarly, clicking the Generate
Load On Other Members button for the selected column member will generate equivalent loads on all Other
member columns having similar height. For example, a pier model consist of two pile types (Pile Type 1, and
Pile Type 2) with similar batter condition and soil set. If nodes on the pile consisting of Pile Type 1 is selected,
then the Other Members are all other pile consisting of Pile Type 1.
Note: The Generate Load On Other Members button will NOT generate loads on the selected piles nodes,
but rather on the Other Members of the selected pile/column (if present) in the model.
Each time the Generate Load On The Selected Member or Generate Load On Other Members button is
clicked, the generated load is displayed in the Generated Load table. The Generated Load table is viewable
after clicking the Show Table button (check Figure: 2.12.j ) located at the bottom right of the dialog. The
Generated Load History drop-down menu can be used to check the previously generated load history. Users
can clear the generated load history by clicking the Clear Generated Load History button. Clearing the load
history does NOT delete the loading from the model.
Figure: 2.12.l Flow chart with empirical equations for calculating barge bow force-deformation (Getter and
Consolazio 2011)
User can input 'User-Defined Elastic-Plastic Curve', Or 'User-Defined General Curve'. For 'User-Defined ElasticPlastic Curve', input the Yield Load and Crush Depth. For 'User-Defined General Curve' enter the Crush Depth,
the Loading and Unloading Curves.
Figure: 2.14.b Concentrated Mass and Damper in 3D View (Thin Element Mode)
2.15.2 Surcharge
Depending upon the type of surcharge selected, different parameters will be required.
The Pier Rotation Angle specifies a pier rotation about the vertical z-axis. The pier rotation is specified as
clockwise positive in the FB-MultiPier coordinate system and is typically used to model skew or radial piers on a
curved alignment. The Bearing Angle (rotation of the bearings) is independent of that of the Pier Rotation and is
measured about the vertical z-axis.
Select the number of Bearing Rows and specify if the span should be continuous. Specific boundary conditions
can be selected and customized by clicking the Edit Supports button.
Superstructure
Select a Span to edit from the span combo box. The "C/C Length" indicates the span length from the center
bearing line of one pier to the center bearing line of the next pier. Click the "Edit Span" button to edit the span
section properties.
Custom bearing connections can be specified by selecting a boundary condition from the combo box. Boundary
conditions can be Released (free to move), Constrained (prevented from movement), or Custom (user-defined
load-displacement curve).
There are two versions of this dialog that are displayed based on the number of bearing rows requested.
Single Row: Only a single option is available
Click the "Edit Custom Bearings" button to define custom bearings using a load-displacement curve. This curve
can be applied to any of the six degrees of freedom for a bearing connection. A maximum of 20 values can be
used to define a custom bearing load-displacement relationship. Values should be entered from smallest to
largest displacement. Click the "Add" button to add a new load-displacement curve. Click the "Del" button to
remove an existing load displacement curve. Click the "Update Plot" button to refresh the load-displacement plot.
Enter the Transformed Section Properties for the bridge superstructure. There are two methods to do this.
Method 1 is to type them directly into the edit boxes in the Transformed frame on this dialog. Method 2 is to click
the Section Properties and Temperatures button, and have the program compute these properties based on
explicitly input girder and slab section properties on the Section Properties dialog. The program uses these
properties to model an equivalent beam that connects the centerline of two pier caps. The Begin Height and End
Height parameters are used to offset the beam from the center of gravity of the pier cap to the center of gravity of
the span.
Begin Height and End Height are measured from the c.g. of the pier cap to c.g. of the bridge
span. Live Load Height is measured from the c.g. of the pier cap to the c.g. of the Live Load (i.e. at
6 ft above the roadway per AASHTO).
Span End Conditions are set independently for each side of the span. Different end conditions may exist
based on the construction; FB-MultiPier can simulate these conditions by assigning various properties to
the Transfer Beam.
l
Custom User assigned custom properties (Must be selected to enable Custom Properties
Button)
Span end conditions can either be automatically calculated by the program, or alternatively, the end conditions
can be specified on the Bridge Span Properties dialog. In particular, each of the transfer beam, vertical link, and
(if applicable) continuity link properties can be input.
Clicking on "Transfer Beam" will open the Transfer Beam properties dialog and one can select stiff (default), soft,
or custom properties.
The span end spring stiffnesses and span mass values can be automatically generated when OPTS
model is generated from a multiple-pier Bridge Model. Refer Section 6.14 for details.
2.
User can manually calculate span end spring stiffnesses and mass values depending on the bridge
model. For more information on how to calculate the span end spring stiffnesses and mass values click
here.
3.
The user can modify the span end spring stiffness and mass values by using Import Spring/Mass Data
button in the OPTS Data dialog. The spring stiffnesses and mass values can be imported from an .opt
file after clicking on the Import Spring/Mass Data button. The imported values will then display in the
Span End Boundary Spring table and the Span End Lumped Mass text box.
The .opt file can be created using either of the following 2 scenarios:
a.
The .opt file is created when the analysis to convert a bridge model to an OPTS model occurs
(Refer Section 6.13.2 ). The .opt file is automatically placed in a special folder called OPTS_
Formation when the conversion from a bridge to an OPTS model occurs. This OPTS_
Formation folder is located in the same directory as the bridge file that was converted to an
OPTS model. NOTE: In this scenario, all creation and processing of the .opt file is automatically
handled by FB-MultiPier.
b.
It is also possible to manually create the .opt file without running the analysis. To do so, follow
the .opt file format (as described in the spring and mass headers in the Help Manual) and
manually type the values.
Note: Clicking the mouse scroll wheel button will toggle between the Picking mode and the 2D Mouse
Control.
*This action will change the currently selected pile. The currently selected pile is the pile that is displayed on the
'Pile Page' and in the 'Soil Edit' window. In the 'Soil Edit' window, the currently selected pile is displayed along
with its assigned soil set. For more information on the integration between the selected pile and the selected soil
set, click here.
Additional grid points can be clicked to apply this Pile Type. To view the various cross sections in a Pile Type that
has multiple
For Grid Spacing Table, refer to: See "Grid Spacing Table"
Return to the Pile Edit Window page.
1. 3D Mouse Control
a. Hold the left mouse button down and drag to rotate the view
b. Hold the left mouse button and the shift key down and drag to pan the view
c. Hold the Control key and the left mouse button down to enable zooming
i. With the key and button pressed down move forward to zoom in
ii. With the key and button pressed down move backward to zoom out Zoom Feature Tutorial
2. Picking Node Mouse Control
3. Picking Element Mouse Control Allows the user to select items in the view to edit in certain dialogs
a. Pick end nodes in the extra members dialog
b. Pick the loaded nodes in the load dialog
c. Pick the node to apply springs to in the spring dialog
d. View the coordinates of the node in most dialogs
4. Add/Remove Pier Cap Element
a. Click on pier cap element to remove or to add it back to model after it has been removed.
5. Piles show or hide the piles in this window
6. Caps show or hide the pile cap
7. Nodes Hide or show the various types of model nodes
a. All Nodes
b. Pier Nodes
c. Span Nodes
d. All Shell Nodes
e. Bearing Locations
f. Pile Cap Nodes
g. Pile Nodes
h. Column Nodes
i. Cantilever Nodes
j. Pier Cap Nodes
8. Pier show or hide the pier structure
9. Girders this feature is not yet available
10. Loads show or hide the loads
11. Springs show or hide the springs
12. Dampers show or hide the dampers (for models using dynamic analysis)
13. Masses show or hide the masses (for models using dynamic analysis)
14. Soil show or hide the soil
15. Retained Soil show or hide the retained soil (for retaining wall models)
16. Numbering show or hide the node numbering
17. Axes (Local) show or hide the local axis
18. Axes (Global) show or hide the global axis
19. Pier Data show or hide the pier global coordinates
20. Thin Elements display the model in thin element mode
21. XZ Plane View show the model as rotated in the XZ plane
4 Program Results
1. Pile Results
2. Pier/Bent Results
3. Pile Interaction Diagrams
4. Pier Interaction Diagrams
5. 3D Results
6. XML Report Generator
Click the Redraw Curves button to update the drawing in the plot windows. This button should be clicked after
changing which forces are selected among the various force checkboxes.
For a time step analysis, select a member force combo box item to display the maximum member force, location,
and corresponding time step.
Click on one of the plot windows to display the maximum and minimum values.
Return to the Pile Results or the Pier/Bent Results page.
4.1.3 Graphs
Graphs are plotted corresponding to the colored piles on the Pile Selection view.
This dialog is reached by right clicking in any plot window that contains data on the Pile Results Page.
- Print: Print the graph. Clicking this option will open the graph as a bitmap in your computer's picture viewer (for
example, "Window Picture and Fax Viewer"). >From here, the graph can be printed. This will allow the graph to
be printed without the print dialog displaying over the graph. (To print the entire dialog, click the 'Print' button at
the bottom of this dialog).
Table Options:
- Print: prints the graph and its contents. (To print the entire dialog, click the 'Print' button at the bottom of this
dialog).
- Save Data: saves the table data to a text file.
*Ifmore than one pile/column if graphed, the member with the maximum value will be displayed in the graph
title.For example, (max at Pile 1).
**The colors used to plot the curves are the identical colors used on the Pile/Pier Forces Dialog.
Printable Forces Tutorial
Return to the Pile Results page.
Click the Redraw Curves button to update the drawing in the plot windows. This button should be clicked after
changing which forces are selected among the various force checkboxes.
For a time step analysis, select a member force combo box item to display the maximum member force, location,
and corresponding time step.
Click on one of the plot windows to display the maximum and minimum values.
Return to the Pile Results or the Pier/Bent Results page.
4.2.3 Graphs
Graphs are plotted corresponding to the colored piles on the Pier Selection view.
- Print: prints the graph and its contents.(To print the entire dialog, click the 'Print' button at the bottom of this
dialog).
- Save Data: saves the table data to a text file.
*Ifmore than one pile/column if graphed, the member with the maximum value will be displayed in the graph
title.For example, (max at Pile 1).
**The colors used to plot the curves are the identical colors used on the Pile/Pier Forces Dialog.
Printable Forces Tutorial
Return to the Pier/Bent Results page.
The Data Table shows a list of all points along an Pile Interaction Diagram curve. To display the table, right click
in the Interaction Diagram window. This will launch a popup menu. Then select the 'Data Table' menu item.
Clicking the "Pier Interaction" on the menu opens a drop down list where the user can select the type of
interaction diagram to view AND selection mode (Thick/Thin)
Type of Pier Interaction Diagram:
1. Biaxial Moment Interaction
2. Uniaxial Moment..2 Axis
3. Uniaxial Moment..3 Axis
Please check out "Interaction Diagrams" for more information.
4.5 3D Results
View the three-dimensional results of the analysis using the following windows:
1. 3D Display Control
2. 3D Results Window
e.Rotation About Y
f.Shear 2
g.Shear 3
h.Moment 2
i.Moment 3
j.Axial
k.D/C Ratio
l.Soil Axial
m.Soil Torsional
n.Soil Lateral X
o.Soil Lateral Y
p.Allow Multiple Forces see Max Min Forces Dialog
Return to the 3D Display Control page.
Figure: 4.5.b Max and Min Forces Dialog for Piles (Left) and Piers (Right)
Displacement vs. Time. This includes graphs for X Translation, Y Translation, and Z Translation.
2.
Select the node you want to plot the results for. Then use the "Results Plotting" combo box to select the desired
graph. Then click the "Plot" button.
default speed is 24. To view the model for a specific time step, ensure the animation is paused (by clicking the
Pause button). Then, use the "Select Time Step" slider to select a time step. The model will draw in the 3D
Results window for the selected time step, and the displacements for the selected node will be shown in the 3D
Display Control dialog.
View the results of the analysis. Elements that have a demand/capacity (D/C) ratio exceeding 1.0 are shown with
a red (highlight) marker.
"Right Click" in the 3D Result window to bring up the view edit menu with the following options:
1. 3D Mouse Control
a. Hold the left mouse button down and drag to rotate the view
b. Hold the left mouse button and the shift key down and drag to pan the view
c. Hold the Control key and the left mouse button down to enable zooming
i. With the key and button pressed down move forward to zoom in
ii. With the key and button pressed down move backward to zoom out Zoom Feature Tutorial
2. Picking Mouse Control Allows the user to select items in the view to edit in certain dialogs
a. Pick end nodes in the extra members dialog
b. Pick the loaded nodes in the load dialog
c. Pick the node to apply springs to in the spring dialog
d. View the coordinates of the node in most dialogs
3. Picking Forces Control Allows the user to view result forces for selected element.
4. Nodal Plot
5. Nodal Results Data
6. Pile Cap Forces
7. Bridge Span Force Plots
8. Nodes
a. All Nodes
b. Transfer Beam Nodes
c. Bearing Nodes
d. Pier Cap / Column Nodes
e. Cantilever
f. Pile Cap Nodes
g. Pile Nodes
9. Node Numbering
10.
c. Pile Elements
d. Column Elements
e. Pile Cap Elements
f. Pier Cap / Bent Cap
13. Axes (Local)
14. Axes (Global)
15. Pier Data
16. Plastic Hinge Zones
17. Undisplaced Model
18. Bridge View Display the full bridge structure
19. Reset View Return the view back to the default setting
3. Show Point Values - This is a tool to help view the plotted data. When selected (checked), the displacement
(or rotation) values display next to their corresponding point.
Node119I_Axial txt,
Node119I_DCRatio txt,
Node119I_SoilLatX txt,
Node119I_SoilLatY txt,
Node119I_SoilAxial txt,
Node119I_SoilTors txt.
By default, these files are saved at the same location as the input file. The contents of each file can easily be
copied into another application for further use, such as EXCEL. The format of the text file is shown below. Each
line in the file contains the load case and its corresponding value (displacement, rotation, force, etc). The number
of lines in the file equals the number of load cases / load combinations / time steps (this particular example has 5
load cases):
1 0.0225
2 0.0667
3 0.1385
4 0.2604
5 0.5097
Other features on this dialog include the following:
1. Node Type - Many nodes in the model occupy both an I node position on one element, and a J node position
on an adjacent element. The I node and J node buttons are used to specify which of these node types is desired
for the currently selected node. (Some exceptions to a node being both an I node and a J node would be pile tip
nodes and cantilever tip nodes).
2. Element Type - Pier Cap nodes that are positioned at the top of a column exist as both a pier cap element
node and a column element node. Thus a choice must be made as to which element type the selected nodes
forces are desired. The same is true for pile head nodes in Pile Bent models, which exist as both a pile element
node and a bent cap element node.
Additional Notes:
1. This feature is not currently available for Pile Cap nodes. Pile Head nodes are NOT considered pile cap
nodes; they are pile nodes, and thus can be used with this dialog.
2. The Soil Axial forces saved by this dialog do NOT include the Pile Tip Spring Force, as shown in the output
(.out) file in the Near Field Soil Resistance section.
Return to the 3D Results Window page.
1) Open and run an analysis with a model that contains a pile grid size of 2x2 or greater (e.g., the default General
Pier model).
2) Click the 3D Results button on the top toolbar, and right-click within the 3D Results window. Then select the
Pile Cap Forces option from the popup menu ( Figure: 4.5.k ).
3) Left-click on a node of interest in the 3D Results window to select it. Press and hold the Ctrl key to select
any number of additional nodes along a pile cap gridline. Alternatively, to automatically select all nodes along a
gridline, check the End Point Nodal Selection checkbox on the Pile Cap Forces dialog. Then, left-click on the
first node of interest and (while holding the Ctrl key) left-click on the last node on the gridline. To ensure
meaningful results, please refrain from skipping any nodes along the gridline while making multiple-node
selections.
4) Select the desired moment (about Xp or Yp) (e.g., as shown for the pier model in Figure: 4.5.l ). The
maximum and minimum moments in the desired direction can be found across all load cases, or for a selected
individual load case. The current load case can be changed using the Load Case edit box on the program
toolbar.
5) Click the Generate button to display the maximum and minimum moment plots. Absolute maximum shear
parallel to Zp is also calculated along the axis on which the moments (Xp or Yp) are acting ( Figure: 4.5.l ).
Figure: 4.5.l Maximum and minimum moment (about Xp) plots along the gridline at the base of columns of the
default general pier.
Return to the 3D Results Window page.
You will be prompted to Allow blocked content. Click the Allow blocked content button to allow the blocked
content.
You will then need to browse to the location were the input file (.in) is saved. The .XML file, which will have the
same name as the input file (.in), should be in the same location as the input file (.in).
Click Yes to allow the ActiveX control interaction. If prompted, click the Allow blocked content button to allow
the blocked content.
There are two web pages from which to customize the data to be viewed. The first page is the Model Data
Report, where the data you have selected is displayed.
The second page is called the Model Data Menu, where you select the data you want to display by checking
the checkboxes (for example, the Print Text Report checkox). Then click the Click to Update Display button
which updates the data on the Model Data Report page.
5 Design Tables
1. Overview
2. GUIOptions
3. Condensed Tables
4. Extended Tables
5. Piercap Envelopes
6. User Walkthrough
5.1 Overview
The design tables represent a summary of the pertinent force results required for the design of the bridge pier(s)
modeled. Various tables and a full report are available for convenient review and are suitable for reporting
results to clients. Controlling forces and load cases are identified, as well as structural element identification and
location in the structure. For the pier cap, design force envelopes are plotted for convenient review and checking.
Figure: 5.3.b Maximum Demand Capacity Ratio for All Pile Sections
.
Figure: 5.4.b Shear Envelope
6 General Modeling
1. Column Connection to the Pile Cap
2. Taper Modeling
3. Node Numbering
4. Span Length Calculation
5. Preliminary Soil Values
6. Buried Pile Cap
7. Bearing Positioning Troubleshooting
8. Extra Member Positioning with Respect to Model Changes
9. Soil and Pile Integration in Program Windows
10. Step-By-Step Guide (Distributed Load Dialog)
11. One Pier Two Span (OPTS) Modeling
12. Coupled Vessel Impact Analysis (CVIA)
In FB-MultiPier loads are transferred between pier columns and foundation members (e.g., piles) via the pile cap.
However, the pier column bases in FB-Multipier are modeled so as to prevent the generation of unrealistic stress
concentrations, as load is transferred from the pier column base into the pile cap nodes. Namely, internal forces
that develop in pier column bases are spread across four adjacent pile cap nodes, as shown in the figure below.
Accordingly, the propensity for numerical, fictitious stress concentrations is reduced.
The load-spreading process is carried out internally by the program: based on the coordinates of the column
base node, the four nodes adjacent to the base node are automatically identified, and "connector" elements are
defined between the four pile cap nodes and the column base node. In this way, the load is distributed to the pile
cap in the form of a column "footprint". The connector elements that are created are "rigid" elements, where the
connector element properties are automatically assigned relatively large constitutive and geometric properties.
In particular, the connector element stiffnesses are proportional to the properties of the columns, where the
connector element stiffnesses are amplified (relatively) by a factor of 10,000. Importantly, the weight of the
connectors are maintained at nominal levels so as not to unduly influence the system configuration. For pileand-cap models that contain stiffness nodes (i.e., where no pier columns are present), the connector element
stiffnesses are assigned in proportion to the pile cap, where again, a stiffness amplification factor of 10,000 is
used. The end conditions at the connectors are defined so that there is no moment transfer at the ends of the
connectors. Therefore, the connectors facilitate manifestation of column-end bending behavior without
generating localized moments at the connector-ends. Salient properties associated with the nodes that the
"connector" elements frame into are listed in the output file in the paragraph "PROPERTIES FOR CONNECTOR
ELEMENTS".
Figure: 6.2.b Engine Model Discretization (Solid View) based on taper input data
Figure: 6.4.a Bearing Offset Measured from the Centerline of the Pier Pile Cap to the Centerline of the Bearings
In Figure: 6.4.a , y1, y2, y3 are local pier origins. Y Coordinate can be specified by the user on the "Bridge
Page".
Span length between two piers each with single row of bearings:
L = (y2 - y1) - (Pier1's pile cap width / 2) + (Pier2's pile cap width / 2)
Span length between two piers (eg. Pier 2 and Pier 3) each with double row of bearings:
L = (y3 - y2) - (Pier2's pile cap width / 2 + Pier2's bearing pad offset (right row)) + (Pier3's pile cap width / 2 Pier3's bearing pad offset (left row))
Example y1 = 0ft; y2 = 60ft; y3 = 120ft; Pile caps = 16ft x 16ft; Bearing offset = 0.75ft.
Single row of bearings pier 1, double row bearing pier 2 and 3. So Span 1 Length is given as:
L1 = (y2 - y1) - (Pier1's pile cap width / 2) + (Pier2's pile cap width / 2 - Pier2's bearing pad offset (left row))
= (60ft - 0ft) - (16ft /2) + (16ft / 2 - 0.75ft) = 59.3 ft
Span 2 Length is given as:
L2 = (Pier2's pile cap width / 2 + Pier2's bearing pad offset (right row)) + (Pier3's pile cap width / 2 - Pier3's
bearing pad offset (left row))
When the user-inputted bearing spacings cannot be accommodated by FB-Multipier, the following error
message displays: "Error Generating Nodes for the Bearings." The specific cause of this error is most often the
result of the close proximity of bearings in one row versus the other row, such that an unacceptably small
element would need to be created between the two bearings (see Figure 2). So the question arises What is the
smallest allowable element length? This is a model-specific value. It is determined as follows: the maximum
number of pier cap (bent cap) nodes for any model is 256. Thus, the maximum number of pier cap (bent cap)
elements is 255. The minimum allowable bearing spacing is the pier cap length divided by 255. Positioning
bearings closer than this length results in an element length that is unacceptably short.
The figures below help to better explain these concepts. The figure below illustrates Bearing Spacing, a term
specifically defined as the user-inputted spacing between bearings within the same row (ie, bearings within the
left row, or bearings within the right row).The figure below also illustrates a case where the column offset is zero
in both rows, as the first bearing in both the left and right row are aligned with the centerline of column 1.
Given the following input data,the figure below illustrates a case where the Error Generating Nodes for the
Bearings error message is caused by too small a distance between bearings, thus creating an element length
that is unacceptably small.
Pier Cap Length = 50.000 feet
Number of Left Row Bearings = 3
Number of Right Row Bearings = 4
Left Row Bearing Spacing, Uniform = 20.000 feet
Right Row Bearing Spacing, Uniform = 13.333 feet
Distance From First Bearing in Left Row to Leftmost column = 0.000 feet
Distance From First Bearing in Right Row to Leftmost column = 0.000 feet
Bearing Xp coordinates in the Left row = 0.000 feet, 20.000 feet, and 40.000 feet
Bearing Xp coordinates in the Right row = 0.000 feet, 13.333 feet, 26.666 feet, and 39.999 feet
When modeling changes are made after an extra member has been applied, it is possible that the position of the
extra member could be automatically adjusted by FB-MultiPier. If and how the extra member is repositioned
depends upon the type of modeling change that has been made. There are two types of modeling changes that
are relevant to extra member positioning: a) changes to the total number of nodes in the model, and b) changes
that require nodal positioning adjustments but do not affect the total number of nodes in the model.
Changes to the total number of nodes in the model: There are numerous modeling changes of this type. Some
changes are obvious, such as changing the number of column nodes, beam nodes, pile nodes, or free length
nodes. Other changes that affect the number of nodes are less explicit, such as changing the number of piles,
eliminating a free length, or converting from one bearing row to two bearing rows. When changes of this type are
made, FB-MultiPier attempts to retain the original position of the extra member, by changing the extra member
nodal assignments (ie, the node numbers used for the extra members i-node and j-node). Specifically, the extra
members previously assigned i-node and j-node are replaced with the nodes that are now closest to the extra
member elements original position. How much the extra member is shifted away from its original position
depends on how far away the closest nodes are to the extra members original position. Figure: 6.8.a and
Figure: 6.8.b illustrate an example of this (all pile cap and pile nodes have been hidden in these figures to
simplify the viewing of the model). A user has applied an extra member to the structure between columns 1 and
2. The user has assigned the extra member an i-node of 85 and a j-node of 95 (see Figure: 6.8.a ).
Figure: 6.8.b Column Node count has been increased, and i-node and j-node are reassigned.
Changes that require nodal positioning adjustments but do not affect the total number of nodes in the model:
There are numerous model changes of this type as well. For example, changing the pile free length or cantilever
length, changing the column height, changing the pile length, etc. These types of changes do not create or
delete nodes. FB-MultiPier simply redistributes the existing nodes, thereby keeping the total number of nodes
constant. When this type of modeling change is completed, the extra member maintains its nodal assignments.
This means that the extra members position in the model could change, because the extra members i-node and
j-node could change position. Figure: 6.8.c and Figure: 6.8.d illustrates an example of this. In this model, the
pile free length is 10. The pile free length is the distance from the pile cap centerline to the top of the soil. An
extra member has been applied at the bottom of the free length, at an elevation of -10 (see Figure: 6.8.c ).
There are 6 free length nodes and 5 nodes in the soil.
The model is then modified by a user. The user modifies the model by removing the top soil layer. This soil layer
had a top of layer elevation of -10 and a bottom elevation of -30. Thus, removing the soil layer extends the pile
free length down to an elevation of -30 (see Figure: 6.8.d ). The free length nodes are then redistributed over
the now-larger free length. Note that the number of free length nodes has NOT changed, nor has the number of
nodes in the soil. There are still 6 free length nodes 5 nodes in the soil. However, the element lengths have
increased in the free length because these nodes must now evenly divide a 30 free length instead of the original
10 free length. Conversely, the element lengths have decreased on the portion of the pile that is in the soil,
because the 5 soil nodes are now divided over a shorter portion of the pile. The extra members i-node and jnode shift downward to an elevation of -30 (the new elevation at the bottom of the free length). Because the
extra member maintains its i-node and j-node assignments, the extra member also moves down to an elevation
of -30.
When this action is taken, if the currently selected pile uses the selected soil set, then the pile will also display in
the Soil Edit window, along with its nodes. Clicking on these nodes will display the soil curve in the Printable Soil
Graph Dialog. However, if the currently selected pile does NOT use the selected soil set, then nodes will not be
drawn on the pile. This indicates the pile does not use the displayed soil set.
Using this mode, soil sets are assigned to pile. After each soil set is assigned, it will display in the Soil Edit
window, along with the pile that is located within the assigned soil set in the pile grid. Note in Figure: 6.9.f that
soil set 2 has just been applied to pile 5. The Soil Edit window now displays Pile 5 and its assigned soil set,
which is Soil Set 2.
This piles nodes will also display. Clicking on these nodes will display the soil curve in the Printable Soil Graph
Dialog.
When the selected pile is changed, the Soil Edit window will display it, along with the piles assigned soil set.
The following is list of ways to set/change the currently selected pile:
a) On the Pile Page, make a selection in the Pile edit box. See Figure: 6.9.d
When this action is taken, the currently selected pile will display in the Soil Edit window, along with its nodes, as
well as this piles assigned soil set. The assigned soil set for this pile also is displayed on the Pile Page, in the
Soil Set combo box.
b) In the Pile Edit window, use the Add/Remove Pile menu option.
A convenient feature is available that facilitates generation and application of distributed loads to FB-MultiPier
models. Namely, the Distributed Load feature allows for the generation of constant value loads, uniform loads, or
linearly varying loads, which can then be applied to piles, pile caps, columns, and pier caps.
As a step-by-step demonstration of this feature, consider Figure: 6.10.a , which is a model view of a single pile
located in a stream bed. For this scenario, the horizontal load associated with the stream is taken into
consideration, and the distributed load generator is used in creating and applying the loads to the model. For the
remainder of this article, the series of steps taken within the FB-MultiPier graphical user interface are listed to
showcase how such loads are generated and applied. Please refer to the Distributed Load page for additional
details about the graphical user interface components in the Distributed Load dialog.
Step-by-step preliminaries: Consider a pier with nine driven piles, which are partially embedded in the ground,
and where free standing water is present above the ground surface. Now, temporarily focus upon a single pile
within the pier ( Figure: 6.10.a ), and assume that the height of the water head above the ground surface is 30 ft
and that the pile diameter is 3 ft. The water is subject to a current, and therefore, horizontal loading acts upon the
pile as a function of height below the waterline (linearly varying from 0 kip/ft to 5.62 kips/ft at the ground surface).
The Distributed Load dialog greatly facilitates creation and application of the linearly varying load on the pile
nodes positioned above the ground surface. Follow the steps below to apply the load using Distributed Load
dialog.
Figure: 6.10.a A single pile excerpted from the pier scenario, partially embedded in soil and subjected to
horizontal loading from water currents
Step 1: Open a model in FB-MultiPier and model the pier, ground surface, and waterline.
Step 2: Click the Distributed Load button on Load page. This will open the dialog shown in Figure: 6.10.b .
Step 7: The Generate Load On Other Members can then be selected , which for this case, automatically assigns
the same load group to the respective nodes along each of the other piles in the pier model ( Figure: 6.10.d ).
Step 8: Click OK to apply the loads that have been generated.
Figure: 6.10.d Generation of load along all other piles of the same type
Figure: 6.11.c Message warning about the loss of all load data except the load case one.
4. Dynamic analysis parameters are considered as shown in Figure: 6.11.d . On the Dynamics Page, select
the Vessel Collision dynamic loading option.
Figure: 6.11.e Message warning about the loss of existing load function data
6. Navigate to the Load Page ( Figure: 6.11.f ), to apply the dynamic vessel collision load
Figure: 6.11.i Flow chart with empirical equations for calculating barge bow force-deformation (Getter and
Consolazio 2011)
9. Click OK at the bottom of the Vessel Collision dialog to return to the Load Page. Within the 3D View located
to the right of the Load Page:
1. Select the node on the pier that corresponds to the impact location (Node 17);
2. Click the Add button on the Load Page to add the impact location to the load table;
3. Click on the letter 'S' just to the left of the newly added node number to toggle the load type from static (S) to
dynamic (D).
Note that the direction in which the vessel collision will be generated and applied is automatically determined,
given the vessel velocity data previously defined above in Figure: 6.11.h . For more details, please refer to
Section 6.12
Also, note that for vessel collision analysis, all other applied nodal loads are treated as static (S) loads. Along
these lines, note that the Dynamic Relaxation feature is typically activated and is utilized to initialize the system
in equilibrium with permanent (static) loading, while minimizing non-physical system vibrations due to the
presence of permanent loads.
Figure: 6.11.j Applying the vessel collision load at the point vessel impact
10. The CVIA model is ready to run. After carrying out the analysis, the user can check structural
demand/capacity ratios on pier components of interest.
11. The vessel impact load-history is generated by the analytical engine and stored in the .VES file ( Figure:
6.11.k ) and is saved in the same folder as the input file. An excerpt of the .VES file is shown in Figure: 6.11.l .
Figure: 6.11.k .VES file with impact load-history (for first 6 time steps only)
6.12 Defining the Coupled Vessel Impact Analysis (CVIA) Collision Scenario
The following discussion assumes that the Coupled Vessel Impact Analysis has already been activated. If this is
not the case, then please first refer to the step-by-step guide for activating CVIA within FB-MultiPier models.
Navigate to the Load page to define the vessel collision scenario.
1. Select the node on the pier that corresponds to the impact location (e.g., Node 17 as shown in Figure: 6.12.a
);
2. Click the Add button on the Load Page to add the impact location to the load table;
3. Click on the letter 'S' just to the left of the newly added node number to toggle the load type from static (S) to
dynamic (D).
way, emphasis can remain on characterizing bridge response to vessel collision loading, rather than on forming
an appropriate impact load-history. However, in case it is of interest, the impact load-history is printed to an ACSII
file with extension .VES, where this file is located in the same directory as all other analysis output.
As a further convenience to engineers wishing to make use of the CVIA feature, only the initial vessel X and Y
(global direction) velocity components ( Figure: 6.12.c ) need to be specified to define the vessel orientation.
Other orientation-related inputs, such as activation of the respective degree(s) of freedom (DOF) dynamic
loading flags, are automatically handled by FB-MultiPier. For example, as shown in Figure: 6.12.c , a non-zero
Vessel X Velocity is defined, and accordingly, FB-MultiPier automatically activates the X dir DOF dynamic
loading flag ( Figure: 6.12.d ).
Figure: 6.12.d Flags to indicate loaded degree(s) of freedom (handled automatically by FB-MultiPier)
5. As further demonstration, consider a scenario where both (global) X and Y velocity components are needed to
define the vessel collision scenario. Accordingly, click on the Vessel Collision button again, and within the
Vessel Collision dialog, enter a non-zero Vessel Y Velocity as shown in Figure: 6.12.e and Click OK.
Figure: 6.12.f Flags to indicate loaded degree(s) of freedom (handled automatically by FB-MultiPier)
The One Pier Two Span (OPTS) model can be modeled using two scenarios as follows:
A. Default One Pier Two Span Problem Type
B. Automatic Generation, Using Pre-Existing Bridge Models
The following steps can be taken to form the OPTS model using the New Problem Type, One Pier Two Span:
1. Open FB-MultiPier, and click the New icon
to open the Select New Problem Type dialog. The default
OPTS model can be selected by clicking on "One Pier Two Span" radio button, and then clicking OK.
Figure: 6.13.c Bridge Page with Edit button to open "OPTS Data dialog
3. Figure: 6.13.d shows the "OPTS Data dialog. The user can see the default spring stiffnesses placed at the
ends of Span 1 and Span 2.
The default spring stiffnesses located at each end of the default OPTS model are generated so as to mimic Span
1, Span 2, and Pier 2 of the default FB-MultiPier bridge model. For more information on how the spring
stiffnesses and mass values were calculated, refer section 6.14 .
In addition to opening the default One Pier Two Span (OPTS) model in FB-MultiPier, one can also generate an
OPTS model from an existing multiple-pier bridge model. Note that multiple-pier bridge models must contain
three or more piers in order to be converted to OPTS models.
The following steps can be taken to form the OPTS model from an existing multiple-pier bridge model:
1. Open any existing multiple-pier bridge model. Open FB-MultiPier program and Click on icon
browsing window. Select and open the multiple-pier bridge model of interest.
to open
For this particular case we will open a default bridge model in FB-MultiPier. Open FB-MultiPier, and click the
New icon
to open the Select New Problem Type dialog. The default bridge model can be selected by
clicking on "Bridge (Multiple Piers)" radio button ( Figure: 6.13.e ). Click OK.
Figure: 6.13.g Selection of "One Pier Two Span" (OPTS) Problem Type
3. The OPTS conversion dialog presents a listing of all piers that are eligible to be used in generating an OPTS
model. Note that the end-most piers are not permitted for conversion to OPTS models. For this demonstration,
Pier 2 is selected from the list.
Figure: 6.13.i Analysis window pop-up while program automatically calculate span end spring stiffnesses and
mass values
5. The OPTS model generated, after the analysis, is as shown below in Figure: 6.13.j .
Figure: 6.13.j OPTS model extracted from default bridge (multiple-pier) model
6. For the OPTS model, the Bridge Page has an additional option where you can view/edit the OPTS Data.
Click the Edit button to open "OPTS Data dialog.
6.14 Automatic Span-End Spring Stiffness and Mass Data Calculation (For One
Pier Two Span Model)
Calculation of span end spring stiffness and mass data for OPTS problem type in FB-MultiPier:
When a One Pier Two Span (OPTS) model is extracted from a full bridge model, the spring stiffness and mass
data at the ends of span 1 and span 2, relative to the selected pier, are calculated automatically.
The following is the process used by the analytical engine to form the span end spring stiffnesses in OPTS
model:
Consider a five pier bridge model as shown below. Select the pier of interest to be converted to an OPTS model
(In this case, Pier 3). As shown below, the program splits the bridge into three distinct portions ( Figure: 6.14.a ).
Isolate the extreme bridge portion (e.g., Left portion of bridge as shown in Figure: 6.14.b )
2.
The node at which the bridge is split is the selected node of interest ( Figure: 6.14.b ). The program
calculates the flexibility matrix of the structure at the particular equilibrium state, following general
principles. To do that the program internally applies unit forces at the selected node of interest. The unit
forces are applied successively in all six possible directions (Fx, Fy, Fz, Mx, My, Mz).
3.
The displacements from each solution at the selected node comprise the columns of the flexibility matrix
(e.g., the displacements from the solution under the application of the Fx unit load comprise the first
column of the flexibility matrix).
4.
After the flexibility matrix is obtained, the program calculates the matrix inverse, which constitutes the
equivalent stiffness at the node of interest.
This calculated equivalent stiffness matrices, for both left and right portion of the bridge, is applied to the
respective span ends of the OPTS model as shown in Figure: 6.14.b .
Figure: 6.14.b Applying equivalent spring stiffnesses at the two span ends of OPTS model.
The following is the process used by the analytical engine to form the span end mass data in OPTS model:
As another example, consider a five pier bridge model as shown below. Select the pier of interest to be
converted to an OPTS model (In this case, Pier 3). Lumped masses are calculated for half of the span to the left
and right of the pier of interest as illustrated in Figure: 6.14.c .
Figure: 6.14.c Applying the lumped mass at the two span ends of OPTS model.
6.15 Coupled Vessel Impact Analysis - One Pier Two Span (CVIA-OPTS) Case
Study
Application of the Coupled Vessel Impact Analysis (CVIA) module and One Pier Two Span (OPTS) Modeling
Technique for Dynamic Finite Element Analysis of Vessel-Bridge Collisions
In the following, an example is developed in a step-by-step manner using the Coupled Vessel Impact Analysis
(CVIA) module and the One Pier Two Span (OPTS) modeling technique for an in-service bridge configuration.
Collision between an aberrant barge and a bridge pier is investigated, for which FB-MultiPier is showcased as a
state-of-the-art design-oriented bridge modeling tool for analyzing dynamic bridge response to vessel collision
loads. Namely, the SR-20 at Blountstown Bridge located in northwest Florida is selected for this example. The
SR-20 at Blountstown Bridge is 1.6 miles long and spans the Apalachicola River in northwest Florida. Vessel
navigation is permitted within a 150 ft channel. A six-span portion of the bridge is modeled in FB-MultiPier (
Figure: 6.15.a ), with Pier 59 as the pier of interest (i.e., the impacted pier).
Figure: 6.15.d Select the Vessel Collision option on the Dynamics page
3. Select Pier #3 and navigate to the Load page ( Figure: 6.15.e ) to apply the dynamic vessel collision load.
of 2,550 kips on Pier 59 (as listed in the structural drawings). Using this static vessel collision load and the
empirical load determination equations given in the American Association of State Highway Transportation
Officials (AASHTO) Guide Specifications and Commentary for Vessel Collision Design of Highway Bridges
(AASHTO 2009), the collision conditions shown in Figure: 6.15.g are back-calculated.
Figure: 6.15.h Flow chart with empirical equations for calculating barge bow force-deformation (Getter and
Consolazio 2011)
6. From the Load page, select the impact location (node), using 3D View and the Load page:
1. Select the node that best represents the impact location on Pier 59 (Node 12 in this case).
2. Click the Add button to add the impact node to the Node Applied list.
3. Click on the letter 'S' just to the left of the newly added node number to toggle the loading type at the
impact location from static (S) to dynamic (D).
Note that the direction in which the vessel collision load will be generated and applied is automatically
determined, given the global X and Y components of initial vessel velocity, which were previously input within
the Vessel Collision dialog. For more details, please refer to Section 6.12
Also, note that for vessel collision analysis, only a single impact location can be defined throughout the
model. All other applied nodal loads are considered to be static (S) loads.
Figure: 6.15.i Applying the vessel collision load at the point vessel impact
7. The model input necessary to make use of the CVIA feature is now complete, and the dynamic vessel collision
analysis is now ready to be carried out. Upon completion of the analysis, the structural demand/capacity ratios
on pier components of interest can be checked in a manner consistent with any other set of FB-MultiPier analysis
results.
8. The vessel impact load-history is generated automatically by FB-MultiPier during the analysis and stored in
the .VES file ( Figure: 6.15.j ). This file (the .VES file) is saved in the same folder as the input file. The full impact
load-history, excerpted from the .VES file generated as part of the current example, is shown in Figure: 6.15.k .
The following is a step-by-step guide for creating an OPTS model, which makes use of the previously discussed
SR-20 bridge model.
1. Open the SR-20 full bridge model in FB-MultiPier and navigate to the Problem page ( Figure: 6.15.l ).
Figure: 6.15.m Select the of "One Pier Two Span" (OPTS) Problem Type
3. When the OPTS Conversion dialog opens, select Pier 3 (which corresponds to Pier 59 of the SR-20 bridge),
and click the OK button ( Figure: 6.15.n ).
Figure: 6.15.o OPTS model (with Pier 59), created from the SR-20 full bridge model
5. For the OPTS model, the Bridge Page has an additional option where engineers can modify the stiffness and
mass data at the model extents, if so desired. Such modifications are made by navigating to the Bridge Page,
and then clicking the Edit button to open the "OPTS Data dialog. Recall that the spring stiffness and mass data
at the ends of span 1 and span 2, relative to the selected pier, were calculated automatically, and so, care should
be taken if it is elected to modify these values. For more information on how to calculate the spring stiffnesses
and mass values please refer section 6.14 .
The Span (deck) element is used to simulate the behavior of the bridge deck. This is modeled as a
series of linear frame elements (default = 10) with properties specified by the user.
2.
The Vertical Link transfers loads from the span (deck) to a transfer beam. The vertical link is used to
account for the eccentricity of the centerline of the bridge deck relative to the centerline of the pier cap.
Also, the vertical link properties can be user defined, or automatically computed based on the span
(deck) end properties.
3.
The Transfer Beam transfers load between the bearings and the vertical link. The Transfer Beam
properties are either based on custom user input or are generated by the program based on the Span
end conditions.
4.
The Bearings are modeled as six springs to represent the response of the bearing in all degrees of
freedom. The spring properties are either based on custom user input or are generated by the program
based on the Span end conditions.
5.
The Offset Rigid Links are used to model the offsets between the bearings to the pier cap centerline.
For piers that are modeled with two rows of bearings, the offset distance is user defined.
6.
The Pier Cap is modeled as a nonlinear frame element projected along the centerline of the physical
pier cap.
The Temperatures button launches the "Temperatures" dialog. This button is only enabled when the "Elastic
Thermo-Mechanical Analysis" check box is checked (selected).
Figure: 7.6.a Cross-section Dimensions (for clarity, only a single girder is displayed)
The series of dialogs that are accessed in making use of the transformed section properties feature are
described as follows. For any bridge model, navigate to the Bridge page and click Edit Span. Then, click
Section Properties and Temperatures ( Figure: 7.6.b ). Next, click the Compute Transformed Section
Properties checkbox to enable access to the Section Properties dialog ( Figure: 7.6.c ). Please see "Section
Properties" dialog for instructions on inputting the cross section data into FB-MultiPier (FBMP). Shown in Figure:
7.6.d is the Section Properties dialog, which houses detailed span element cross section data, where all of the
dimensions denoted in Figure: 7.6.a can be input. Subsequent to populating the detailed span section
properties, and clicking OK on the Section Properties and Temperatures dialog, a notification will appear to
confirm that the transformed section properties for the span will indeed be overwritten.
As demonstration of this feature, assume that the steel girder components (flanges, webs) from "Cross-section
Dimensions (for clarity, only a single girder is displayed)" are transformed to concrete components based on the
modular ratio n, where n = Egirder / Eslab = 29,000/3,625 = 8. Accordingly, the transformed section (in this
scenario) utilizes an elastic modulus value of 3,625 ksi. For the remainder of this discussion, automated
calculations of transformed section properties are compared to corresponding manual calculations.
The manually calculated transformed section properties are:
Transformed Cross Section Area:
Area of Single Steel Girder, As = (3.25 in*26 in+0.69 in*108 in + 3.25in *26in) = 243.52 in^2
Area of Concrete Slab, Ac = 9 in *656 in= 5904 in^2
Transformed Area, A = Ac + 4* As*n = 5904 in^2+ 4*243.52 in^2 *8 = 13696.64 in^2
The manual calculation of transformed area A matches the corresponding value generated using FBMP (
Figure: 7.6.b ).
The "Girder Spacing Typical Section" button launches the dialog shown in Figure: 7.7.b . The "Girder Spacing
Typical Section" dialog displays the spacing between the girders. Spacings are measured from the center-line of
each girder. The girder spacings should be set via this dialog.
Note: The girder spacing provided in the "Girder spacing Typical Section" is independent of the bearing
spacings, which are set via the "Bearing Locations" dialog launched from the "Pier" page.
Figure: 7.7.d Girder Dimensions (Expanded View) Non Uniform Girder Configuration
When girders have uniform shape and variable height, one set of girder dimensional data is input per element.
Note there are ten (10) rows of input data, as this model contains 10 span elements ( Figure: 7.7.e ).
The "Table PDF" button creates a .PDF file of the data displayed in the "Girders" table.
The "Plot" button draws the 2-dimensional image of the slab and girders in the "Plot" window ( Figure: 7.7.g ). It
is recommended to click the "Plot" button after making changes to the dimensional input, to ensure the image is
based on the most current input.
Properties of the Vertical Link match the axes between the Span (S) and the Vertical Link (VL). The VL properties
can be user defined. Alternatively, if the VL properties are to be automatically calculated by the program, then the
goal is to select properties that are proportional to the ends of the span.
Eqn: 7.8.A
Eqn: 7.8.B
Eqn: 7.8.C
Eqn: 7.8.D
Eqn: 7.8.E
Eqn: 7.8.F
Where,
LVL Length of the vertical link element.
AVL Cross sectional area of the vertical link element.
JVL Torsional constant of the vertical link element.
I3VL Moment of inertia (I3) of the vertical link element.
I2VL Moment of inertia (I2) of the vertical link element.
EVL Youngs modulus of the vertical link element.
GVL Shear modulus of the vertical link element.
LS Length of the attached total bridge span.
AS Cross sectional area of the nearest horizontal bridge span element.
JS Torsional constant of the nearest horizontal bridge span element.
I3S Moment of inertia (I3) of the nearest horizontal bridge span element.
I2S Moment of inertia (I2) of the nearest horizontal bridge span element.
ES Youngs modulus of the nearest horizontal bridge span element.
GS Shear modulus of the nearest horizontal bridge span element.
node separating each element (nodes 3 through 11); node 12 is located at the top of the right elevation beam;
node 13 is located at the base of the right elevation beam. The number and location of the remaining
superstructure nodes depend on the number of bearing locations. Node 14 is the first bearing location on the left
rigid transfer beam. There is one transfer beam node per bearing location (nodes 14 through 19 as depicted in
figure TR-1). The right transfer beam nodes then follow (nodes 20 through 25 as depicted in figure TR-1).
The transfer beam is used to model the end conditions of the span and to transfer the load to the bearings. The
transfer beams properties are dependent on the properties of the vertical link.
The user can select from one of the three descriptions for the transfer beam:
1. Stiff (default)
2. Soft
3. Custom
Figure: 7.9.b Axis for Transfer Beam and Vertical Rigid Beam
Equation used to generate the transfer beam properties for the stiff (default) and soft options are described
below.
Stiff (default) properties:
Eqn: 7.9.A
Eqn: 7.9.B
Eqn: 7.9.C
Eqn: 7.9.D
Eqn: 7.9.E
Eqn: 7.9.F
Soft properties:
Eqn: 7.9.G
Eqn: 7.9.H
Eqn: 7.9.I
Eqn: 7.9.J
Eqn: 7.9.K
Eqn: 7.9.L
Where,
LTB Length of the transfer beam element.
ATB Cross sectional area of the transfer beam element.
JTB Torsional constant of the transfer beam element.
I3TB Moment of inertia (I3) of the transfer beam element.
I2TB Moment of inertia (I2) of the transfer beam element.
ETB Youngs modulus of the transfer beam element.
GTB Shear modulus of the transfer beam element.
LVL Length of the nearest vertical link element.
AVL Cross sectional area of the nearest vertical link element.
JVL Torsional constant of the nearest vertical link element.
I3VL Moment of inertia (I3) of the nearest vertical link element.
I2VL Moment of inertia (I2) of the nearest vertical link element.
EVL Youngs modulus of the nearest vertical link element.
GVL Shear modulus of the nearest vertical link element.
Eqn: 7.10.A
Eqn: 7.10.B
Eqn: 7.10.C
Eqn: 7.10.D
Eqn: 7.10.E
Eqn: 7.10.F
Eqn: 7.10.G
Eqn: 7.10.H
Eqn: 7.10.I
Eqn: 7.10.J
Eqn: 7.10.K
Eqn: 7.10.L
Figure: 7.11.a Axis for Offset Rigid Link and Transfer Beam
Eqn: 7.11.A
Eqn: 7.11.B
Eqn: 7.11.C
Eqn: 7.11.D
Where,
LRL User input pad offset from pier cap centerline; 1 inch for single row; 1 inch if no offset is specified.
ARL Cross sectional area of the rigid link element.
JRL Torsional constant of the rigid link element.
I3RL Moment of inertia (I3) of the rigid link element.
I2RL Moment of inertia (I2) of the rigid link element.
E Maximum Youngs modulus among all the pier cap elements (excluding cantilever).
G Maximum Shear modulus among all the pier cap elements (excluding cantilever).
LPC Distance from the first column centerline to last column centerline.
APC Maximum cross sectional area among all the pier cap elements (excluding cantilever).
JPC Maximum Torsional constant among all the pier cap elements (excluding cantilever).
I3PC Maximum Moment of inertia (I3) among all the pier cap elements (excluding cantilever).
I2PC Maximum Moment of inertia (I2) among all the pier cap elements (excluding cantilever).
Note: For the Bridge model, it is up to the user to apply the loads that are associated with wind on substructure.
type is set on the "Section Properties" dialog (see the "Type" section in the "Section Properties" dialog topic).
Figure: 7.14.b Figure: 7.14.c , Figure: 7.14.d show the schematics for each girder type.
The "Girder Temperatures and Heights" table is used to input girder temperature data. The schematic in the
"Schematics" section (described above) is useful in understanding the locations on the girder at which the
temperature inputs are applied.
The "Temperatures (Expanded View)" is launched via the "Expand Table" button. This dialog is used to expand
the viewable size of the temperature and heights table, making viewing large sets of data with many table rows
easy. This expanded table also helps illustrate the input data required for each girder configuration, as follows:
When girders are uniform with constant height, one set of girder temperature data is input, and this set of data is
applied to all elements in all girders. Note there is one (1) row of input data ( Figure: 7.14.e ).
The "Table PDF" button creates a .PDF file of the data displayed in the "Girder Temperatures and Heights" table.
The "Plot" button launches the "Plot" dialog.
The "Plot" dialog displays the plot of the temperature data for the selected girder on the selected span ( Figure:
7.14.i ).
The "Girder" combo box specifies the girder for which the thermal gradient plot is displayed. This combo box is
not active when the "Uniform (across all girders)" property is selected on the "Section Properties" dialog is
selected (checked), because with this property selected, there is only one (1) set of temperature data per span.
The "Show Slab Temperatures" check box specifies whether to display the slab temperatures on the plot. This
check box does not change the curve line; it merely controls whether the slab temperatures are visible or hidden.
It is useful when slab temperatures and the uppermost girder temperatures are very close together such that they
overlap. Hiding the slab temperatures in this instance makes viewing the girder temperature data easier to see.
The "Show Girder Temperatures" check box specifies whether to display the girder temperatures on the plot. This
check box does not change the curve line; it merely controls whether the girder temperatures are visible or
hidden. It is useful when slab temperatures and the uppermost girder temperatures are very close together such
that they overlap. Hiding the girder temperatures in this instance makes viewing the slab temperature data easier
to see.
The "Print" button creates a printout of the thermal gradient plot.
8 Setup Options
1. Expanding Memory
2. Program Settings
Any installed version can have its permissions changed by entering encrypted numbers into the license file.
This is done by choosing the Control->Update license option from the main menu. The update can be done by
hand or automatically through the Internet.
E-mail/Fax/Phone License Update
This option is for installations that do not have an Internet connection. To do this installation, call the BSI support
number (check the web for the phone number) and you will be stepped through the process. Numbers from your
computer need to be given to the BSI representative and we can Fax or E-mail the encoded numbers you will
need to type into the program.
Internet License Update
This option requires the computer on which you are installing the license file be connected to the Internet. Then,
all numbers are communicated through the Internet and the license updated automatically. The computer can
either be a stand-alone system or the network server for a multiple seat license.
Transfer License
There is a built in function that allows you to transfer you license to another machine. This allows you to move
the license file from your current server or workstation to a new machine.
Troubleshooting
The license file (both for servers and individual workstations) is locked to a machine based on hardware
components contained in the machine. If you change or modify your hardware (drives, motherboards etc) your
installation may not function. To do this, you should first transfer the license, then modify your hardware, and
then re-install the license on the machine.
Novell systems: Be sure that the directory where the license file is saved is accessible to any user. The user
must have read, write, modify, erase and create rights for that directory.
To update the software license at any time, select Update Software License from the Control menu while viewing
the intro Logo window. Doing so brings up the License Configuration Wizard.
The initial License Configuration Wizard screen shows four options for updating the software license. The
options are shown below:
This option is used for a single installation of the software that does not rely on network to run the program. A
license of this type is individually purchased per machine.
Click the Next button to continue. The next screen allows the user to update the license by Email
([email protected]).Click Next. Please contact the BSI for assistance if needed.
The Update Complete screen will then be shown after successfully entering the numerical codes. In order to
apply the changes to the program configuration, the FB-MultiPier program needs to be restarted. Clicking the
Finish button will update and automatically close the program. The program will now run in an unlocked state.
License Update Tutorial
This option is used for a single installation of the software on a network server. This license update is identical
to stand alone workstation update, except that the license is configured on the network server. This option
would be used to run the program directly on the server to take advantage of the server hardware configuration
(i.e. more memory, hard disk space, etc.). A license of this type is individually purchases per machine.
Select Update a License on a Network Server from the initial screen and follow the steps outline for Updating a
License on a Stand Alone Workstation.
Because this process can not be reversed, the user must check the box to confirm the remove the license from
the current computer before proceeding. Doing so will enable the Next button. Click the Next button to remove
the license.
License Update Tutorial
Click the Next button to continue. The Update Complete page is now shown. In order to apply the changes to
the program configuration, the FB-MultiPier program needs to be restarted. Clicking the Finish button will
update and automatically close the program. The program will now run in an unlocked state.
10 Toolbar Icons
1. Description of Toolbar Icons
2. General Pier Wizard
The General Pier Wizard creates a general pier problem using detail specified information.
11 Batch Analysis
1. Batch Mode
2. RunningFB-MultiPier Engine in Batch Mode
2.
Save the file as run.bat (run is an arbitrary name, the extension must be .BAT).
4.
As an example, here is a batch file that will run the program for three input files.
12 Soil-Pile Interaction
Soil-Pile Interaction section characterizes both the axial and lateral soil-pile interaction. The axial soil-pile
interaction is modeled with hyperbolic t-z curves. The tip resistance is modeled with q-z curves given in axial
soil-pile interaction section. The lateral soil-pile interaction is modeled with nonlinear p-y curves. Four of the p-y
models are the same as those given in FHWA's COM624P manual (See "Wang and Reese").
1. Group Interaction
2. Soil Resistance Due to Pile Rotation
3. Soil Properties
4. Lateral Soil-Pile Interaction
5. Axial Soil-Pile Interaction
6. Torsional Soil-Pile Interaction
3x3
4x3
5x3
6x3
7x3
Average
Lead Row
2nd Row
3rd Row
4th Row
5th Row
6th Row
7th Row
245
178
142
294
205
151
142
294
222
160
151
142
302
205
178
142
142
142
285
222
178
151
142
142
142
284
206
167
148
142
142
142
Group
(Measured)
1664
2375
2909
3336
3790
Group
(Predicted)
1898
2398
2843
3270
3697
Error (%)
14
2.3
2.5
Note that the individual row contributions, with the exception of the trail row, appear to be only a function of row
position. Also, using the average for the row (with exception of trail row) does a good job of predicting the
measured group response. Consequently, the approach recommended by Hannigan et. al. (2006) with Pmultipliers listed in AASHTO 2012 has been implemented in the program.
The following P-multipliers are recommended for lateral loading at 3D pile spacing:
0.8, 0.4, 0.3, 0.3, ..0.3 where 0.8 is the lead row and 0.3 is the trail row value
For 5D pile spacing the following P-multipliers are recommended:
1.0, 0.85, 0.7, 0.7, , 0.7 where 1.0 is the lead row and 0.7 is the trail row value.
These multipliers generally represent group efficiencies of 70-75% for 3D spacings and 95% for 5D pile spaced
groups. Also, the multipliers were found to be independent of soil density (sands).
NOTE: The program will apply the P-multipliers to the correct pile rows (lead to trail) based on the direction the
piles move. The P-multipliers are always given in trail to lead order. This does NOT depend on the direction of
the applied load.
In the case of battered piles (A frame) as shown below:
Modeling of pile inclination (batter) is available in FB-MultiPier, as indicated in Figure: 12.2.a for a pile with
horizontal (H) versus vertical (L) slope. Instructions for making use of the FB-MultiPier UI to incorporate pile
batter into piles are listed here. In the following discussion, soil resistance of battered piles is focused upon,
where considerations made with respect to processing of soil resistance forces and stiffness are detailed.
12.3.2 Buoyancy
The buoyant force on the bridge substructure that is submerged, i.e., below the water table, is automatically
computed if a buoyancy factor greater than 0 is selected in non-AASHTO mode or if buoyancy is activated
(checked on) in AASHTO mode. The computation includes piles, pile cap, pier columns. Partial buoyancy of the
pile cap is accounted for, where the volume of the pile cap that is submerged will be used to determine the
buoyant force acting on the pile cap. A convenient way to check buoyancy and self-weight calculations is to
include only these loads, run the program, and then view the "Sum of Total Soil Spring Loads", Z direction in the
output file.
(psf)
Eqn: 12.3.A
where
= 5 for sands with fines
10 for clean normally consolidated sand
15 for clean overconsolidated sand
= atmospheric pressure ( 2000 psf)
= corrected SPT blow-count (blows/ft)
(psf)
Eqn: 12.3.B
where
= subgrade modulus (pcf)
= width of pile (ft)
= poissons ratio
(psf)
Eqn: 12.3.C
where
= subgrade modulus (pcf)
= depth below ground surface (ft)
For Clay
(psf)
Eqn: 12.3.D
where
= range of beta is shown in the table below
= undrained shear strength (psf)
The following typical values may be used for the Poisson's ratio for soils:
= 0.2 to 0.45 for sand
= 0.4 to 0.5 for clay
or a spatial average, for the values of over depth may be used for soils consisting of both sand and clay.
Eqn: 12.3.E
For Sand
(ksf)
Eqn: 12.3.F
where
= maximum shear modulus (ksf)
= corrected SPT blow-count (blows/ft)
use the values of Young's Modulus, E from equations Eqn: 12.3.A Eqn: 12.3.B and Eqn: 12.3.C into Eqn:
12.3.E to calculate shear modulus for sand.
(psf)
Eqn: 12.3.G
where
= 5 for sand fines
10 for clean normally consolidated sand
15 for clean overconsolidated sand
= atmospheric pressure( 2000 psf)
= corrected SPT blow-count (blows/ft)
= subgrade modulus (pcf)
= width of pile (ft)
= poisson's ratio
= depth below ground surface (ft)
For Clay
Use the values of, Young's Modulus, E from equation Eqn: 12.3.D into Eqn: 12.3.E to calculate shear modulus
for clay.
(psf)
Eqn: 12.3.H
where
= range of beta shown in Eqn: 12.3.D
= undrained shear strength (psf)
Eqn: 12.3.I
where
C
Eqn: 12.3.J
valid only for v 0.25 tsf (24 kPa) (Bowles, 1977 )
Normalizing for atmospheric pressure (pa): (1 atm = 101.3 kPa = 1.06 tsf )
Eqn: 12.3.K
Larger values should be used for granular material with 5% or less of fine sand and silt.
For numerical implementation, the average correlation can be expressed as
Eqn: 12.3.L
where
Eqn: 12.3.M
where
qu = unconfined compressive strength
Figure: 12.3.a Correlations between SPT N-value and Unconfined Compressive Strength
Figure: 12.3.b SPT Blow Count vs. Friction Angle and Relative Density
1. O'Neill's Sand
2. Sand of Reese, Cox, and Koop
3. Sand (API)
4. O'Neill's Clay
5. Matlock's Soft Clay Below Water Table
6. Reese's Stiff Clay Below Water Table
7. Reese and Welch's Stiff Clay Above Water Table
8. Clay (API)
9. Limestone (McVay)
10. Limestone (McVay) 2-3 Option
11. Custom P-Y
Eqn: 12.4.A
where
= a factor used to describe pile shape;
= 1.0 for circular piles;
A= 0.9 for cyclic loading;
= 3-0.8 z/D 0.9 for static loading;
D= diameter of pile;
pu = ultimate soil resistance per unit of depth;
Eqn: 12.4.B
Eqn: 12.4.C
wherez = depth in soil from ground surface;
= effective unit weight of soil;
Ka = Rankine active coefficient;
= (1 - sin )/(1 + sin )
Kp = Rankine passive coefficient;
= 1/Ka ;
Ko = at-rest earth pressure coefficient;
= 1 - sin ;
= angle of internal friction;
= 45o + /2 .
The p-y relationship given in equation Eqn: 12.4.A depends on the soil parameters k (lb/in3 or N/m3 ) and
(deg), which may be obtained from insitu SPT data. For sand, use SPT to find ( Figure: 12.4.b ) and to find k
(F/L) ( Figure: 12.4.c ).
A comparison between O'Neill's p-y curve for sand and Reese et. al. (1974) curve (SOIL=2) is shown in Figure:
12.4.a for =35 k=150 lb/in3 , and buoyant =52.6 lb/ft3 at a depth of 25 ft. Evident from the figure, O'Neill's curve
fits Reese's initially, but differs for Pu (generally the case).
Figure: 12.4.a Comparison of ONeills and Reese, Cox, and Koops P-Y Curves
Figure: 12.4.b SPT Blow Count vs. Friction Angle and Relative Density
Figure: 12.4.d P-Y Curves for Static and Cyclic Loading of Sand (after Reese, et al, 1974)
Eqn: 12.4.D
Eqn: 12.4.E
where z
p'o
D
C1,C2,C3
=
=
=
=
=
The variation of coefficients of C1,C2,and C3 with are given in the graph below.
Eqn: 12.4.F
where pu = ultimate bearing capacity, which is defined as smaller value from the Eqn: 12.4.D nd Eqn: 12.4.E
subgrade modulus, force per volume units, is determined from Figure: 12.3.c , which is a function of
k =
Figure: 12.4.f O'Neill's Integrated Method for Clay (b) Cyclic Loading Case
Figure: 12.4.g ONeills Integrated Method for Clay (b) Static Loading Case
Figure: 12.4.h P-Y Curve for Soft Clay Below Water Surface (Static Loading)
Figure: 12.4.i P-Y Curve for Soft Clay Below Water Surface (Cyclic Loading)
Figure: 12.4.j Reese et al (1975) Cyclic P-Y Curve for Stiff Clay Located Below the Water Level
Figure: 12.4.k Reese et al (1975) Static P-Y Curve for Stiff Clay Located Below the Water Table
Figure: 12.4.l Welch and Reese (1972) Static P-Y Curve for Stiff Clay Above Water Table
Figure: 12.4.m Welch and Reese (1972) Cyclic P-Y Curve for Stiff Clay Above Water Table
Eqn: 12.4.G
Eqn: 12.4.H
where c
p'o
z
D
J
XR
=
=
=
=
=
=
Eqn: 12.4.I
where ' = effective unit weight of soil in weight density units
In FB-MultiPier analysis, the value of J is set equal to 0.5, which is recommended for Gulf of Mexico clays. The
data used in piecewise linear p-y curves of API clay are given in the table:
=
=
=
=
Eqn: 12.4.J
where c = strain occurring at one-half the maximum stress on laboratory undrained compression tests of
undisturbed soil samples
The data for the PY curves presented below is based on the report "Development of Modified T-Z curves for large
diameter piles/drilled shafts in limestone for FBPIER"(McVay et. Al. (2004)). The data for the back computed
curves were obtained from 12 lateral load tests performed in the centrifuge with diameters of 6 and 9 ft,
embedment (L/D) of 2, 3, and 4 and rock strengths of 10 and 20 tsf. (The report recommends that full scale field
tests be employed to validate the curves presented). Each lateral load test gave multiple P-Y curves, which were
averaged to obtain a representative curve.
Presented in Figure: 12.4.p are back-adjusted P-Y curves for all twelve-centrifuge tests with side shear
considerations; i.e., two shaft diameters (6 and 9), three embedment lengths (L/D = 2, 3, and 4) and two rock
strengths (10 tsf and 20 tsf). Also shown in the figure are the predicted P-Y curves for soft and stiff clay models.
Figure: 12.4.p P-Y curves from 12 lateral tests corrected for side shear.
Evident from the figure, even though the lateral resistance is normalized with rock strength and diameter, there is
quite a bit of variability in the P-Y curves. Therefore the curves were normalized even further to be represented
by a single trend-line. The P values are normalized with q 0.15 D0.85. Figure: 12.4.q shows the normalized Pu
Y curves for Florida Limestone corrected for side friction. Note that the curves are valid for all the experimental
results (i.e., 6 and 9 diameter shafts, different rock strengths, etc.). Note also that the P-Y curves are unit
dependent. That is for the English system, the rock unconfined compressive strength (q ), the shaft diameter and
u
rocks lateral resistance, P must be in ksf, feet and kips/ft, respectively. For the Metric system, the rock unconfined
compressive strength (q ), the shaft diameter and rocks lateral resistance, P must be in KN/m2, m and KN/m,
u
respectively. The normalized curves can be obtained by the following equations:
Eqn: 12.4.K
Eqn: 12.4.L
where:
D = Pile diameter
qu = Unconfined Compressive strength
y = Pile Displacement
P-Y curves for Florida Limestone with/without the inclusion of 2-3 rotation are given below. The data for the P-Y
curves presented below are derived from the technical research report "Development of Modified T-Z curves for
large diameter piles/drilled shafts in limestone for FBPIER", McVay et al. (2004). Note that the P-Y curves are unit
dependent: for the English system, the rock unconfined compressive strength (qu), the shaft diameter and rocks
lateral resistance, P must be in ksf, feet and kips/ft, respectively. For the Metric system, the rock unconfined
compressive strength (qu), the shaft diameter and rocks lateral resistance, P must be in KN/m2 , m and KN/m,
respectively. The curves can be obtained by the following equations:
Eqn: 12.4.M
Eqn: 12.4.N
where:
D = Pile diameter
qu = Unconfined Compressive strength
y = Pile Displacement
Note that p is held constant for values of y/D exceeding 0.1
FB-MultiPier can generate two types of P-Y curves for the Florida limestone to allow the user to either include or
exclude the effect of the side shear contribution during the analysis. This option is activated when choosing "Soil
Resistance due to Pile Rotation about 2 and 3 axes" from the soil page, under Soil Layer Models Lateral. More
information on "Soil Resistance Due to Pile Rotation". When the analysis is requested to include the side shear
contribution to the lateral response mechanism then the program calculates the additional term dP. If on the
other hand the option is not selected then the effect of side shear is not calculated. Based on the discussion
above the user should use this feature with the necessary caution and only where the use is justified. That is
when this option is chosen then care must be taken so that the appropriate P-Y curve is used.
Eqn: 12.5.A
where
Eqn: 12.5.B
At a particular location on the pile, is the shear stress being transferred to the soil for a given z displacement,
o
where r is the radius of the pile and r is the radius out from the pile where axial loading effects on soil are
o
m
negligible, assumed to be initially equal to the product of the pile length, (1- soil's Poisson's ratio), and the ratio of
the soil's shear modulus at the pile's center to the value at its tip. The user must supply Gi, the initial shear
modulus of soil, v, Poisson's ratio of soil, and f, the maximum shear stress between the pile and soil at the depth
in question. Evident from the equation above, the side springs are highly nonlinear.
The following data is based on tests which were performed on 6 diameter shafts embedded 18 (L/D = 3) into
the rock and are described in the report "Development of Modified T-Z curves for large diameter piles/drilled
shafts in limestone for FBPIER"(McVay et. Al. (2004) ). All of the plots, Figure: 12.5.b Figure: 12.5.d show
the loaddis place ment data which mobilize significant axial resistance with small displacements (i.e., 80%
capacities at 0.5% of diameter). Axial load tests in lower strength 5 tsf rock, proved unattainable, because the
rock mass fractured from the shaft to the boundaries of the bucket.
The three normalized curves are quite similar and can be represented by a single curve (shown in bold line), with the
following equations:
0 R 0.5
Eqn: 12.5.C
0.5 R 3
Eqn: 12.5.D
3 R
Eqn: 12.5.E
where: R = z/D*100.
fs = skin friction
fsmax = ultimate unit skin friction
Kim (2001) analyzed data from 33 axial load tests (Osterberg) from various bridge sites throughout Florida and
recommended the normalized T-Z curve for the natural Florida Limestone given in Figure: 12.5.g . A
comparison of Kims normalized T-Z curve with the synthetic rock curve, Figure: 12.5.f is also shown in Figure:
12.5.g . Evident from the figure there is a very good agreement between the normalized T-Z behavior of the
natural limestone and the synthetic rock.
Eqn: 12.5.F
where K = dimensionless coefficient of lateral earth pressure (ratio of horizontal to vertical normal effective stress
(for unplugged K=0.8 and for plugged K=1.0)
p'
= effective overburden pressure in stress units
o
= friction angle between the soil and pile wall, which is defined as
Eqn: 12.5.G
where = internal friction angle
It is recommended that the ultimate (limiting) values of unit skin friction, fult, be considered, which are given in
Table G.4.3-1 API RP2A LRFD In FB-MultiPier analysis, the t-z curves are generated using piecewise linear
function based on the data given in the table:
Eqn: 12.5.H
where c = undrained shear strength of the soil in stress units
= a dimensionless factor, which is defined as
Eqn: 12.5.I
Eqn: 12.5.J
Eqn: 12.5.K
where p'o = effective overburden pressure in stress units
The data used in piecewise linear t-z curves of API clay are given in the table:
=
=
=
=
NOTE: In FB-MultiPier analysis, a ratio of residual stress to ultimate stress is set equal to 0.9
The T-Z curves used for drilled and cast insitu piles/shafts are based on the recommendations found in Wang
and Reese (1993). In particular, the curves are based on the trend lines and are computed for each node. Trend
lines of stress transfer for axial end bearing and side resistance are provided for the following materials:
1.
2.
3.
4.
Sand
Clay
Intermediate Geomaterial
User Defined
Sand
Valid for 30
Eqn: 12.5.L
Eqn: 12.5.M
Eqn: 12.5.N
valid for depths ranging from 5 to 87.5 ft (1.5 to 26.7 m)
The immediate settlements are computed using non-linear t-z springs, with the shape presented in following
Figure: 12.5.j . The equations are provided but it should be referred that there is a considerable scatter around
the trend line.
for R 0.908333
Eqn: 12.5.O
where
Eqn: 12.5.Q
Figure: 12.5.j Trend Lines for Drilled Shaft Side Friction in Sand
Clay
The immediate settlements are computed using non-linear t-z springs, with the shape presented in following
Figure: 12.5.k . The equations are provided but it should be referred that there is a considerable scatter around
these trend lines.
Eqn: 12.5.S
Figure: 12.5.k Trend Lines for Drilled Shaft Side Friction in Clay
Intermediate Geomaterial
The design of drilled shafts founded in intermediate Geomaterials is directly from FHWA's Load Transfer for
Drilled Shafts in Intermediate Geomaterials .
Intermediate Geomaterials are characterized as one of the following 3 Types:
1. (Type 1) Argillaceous geomaterials: Heavily overconsolidated clay, clay shale, saprolite and mudstone.
2. (Type 2) Calcareous Rock: Limestone and Limerock
3. (Type 3) Very Dense Granular Geomaterials: residual, completely decomposed rock, and glacial till.
l
Note:
Types 1 and 2 are considered to be cohesive materials with an undrained strength, qu in the range of 0.5
to 5.0 Mpa.
Type 3 is primarily cohesionless and has Nspt from 50 to 100
Method 1 proposed by FHWA's Load Transfer for Drilled Shafts in Intermediate Geomaterials, for Type 1 and 2
materials has been coded herein.
l
Valid for IGM Type 1 and 2; 0.5 < qu < 5.0 Mpa; Recovery > 50 %;
Appropriate for very short sockets (L/D <2) or very long sockets (L/D>20);
Where there is strong layering in the formation, or where part of the socket is artificially roughened and
part is smooth
Required Data:
l
Number of Layers
qu (Mpa)
Mass Modulus - Em
Thickness
Em/Ei
Em/Ei
(percent)
(closed joints)
(open joints)
100
1.00
0.60
70
0.70
0.10
50
0.15
0.10
20
0.05
0.05
Note: Values of Em/Ei for RQD values between those shown can be estimated by linear interpolation on RQD.
User can select the "Custom T-Z" curve option for applying user-defined T-Z curves to subsurface pile nodes that
fall within a given soil layer.
The user defined curves are specified by a set of TEN points. Please enter the data set in two columns.
User defined T-Z data
Z1 T1
Z2 T2
Z3 T3
... ...
Z10 T10
Where
Ziis the ith Z (axial pile deflection) value on the user specified T-Z curve.
Tiis the ith T (shear stress) value on the user specified T-Z curve.
Eqn: 12.5.A
where Qf is the ultimate tip resistance (force), Gi and v are the initial shear modulus and Poisson's ratio of the soil
at the pile tip. r0 is again the radius of the pile, and Qb is the mobilized tip resistance.
The unit end bearing, in stress units, can be calculated by the equation
Eqn: 12.5.B
where p'o = effective overburden pressure, in stress units, at the pile tip
Nq = dimensionless bearing capacity factor, which is defined as
Eqn: 12.5.C
where ' = effective internal friction angle
Based on a total unit weight of soil, an effective unit weight, i.e.soil = soil water is calculated for p'o that is
subsequently used to compute the unit end bearing in Eqn: 12.5.B
It is recommended that the ultimate (limiting) values of unit end bearing qult be considered, which are given in
Table G.4.3-1 API RP2A LRFD.
The ultimate end bearing capacity, in force units, may be calculated by the equation
Eqn: 12.5.D
where A = sectional area at tip of pile and is calculated, based on pile end condition
Plugged a gross sectional area is used to compute ultimate end bearing capacity
Unplugged a cross sectional area is used to compute ultimate end bearing capacity
The q-z curves are generated using piecewise linear function based on the data given in the table:
=
=
=
=
Eqn: 12.5.E
where c
undrained shear strength of the soil, in stress units, at the pile tip
Eqn: 12.5.F
where A = sectional area at the tip of pile, which is based on the pile end condition
Plugged a gross sectional area is used to compute ultimate end bearing capacity
Unplugged a cross sectional area is used to compute ultimate end bearing capacity
In FB-MultiPier analysis, q-z curves of API clay are defined as a piecewise linear function based on the data
given in the table:
The Q-Z curves used for drilled and cast insitu piles/shafts are based on the recommendations found in Wang
and Reese (1993). In particular, the curves are based on the trend lines and are computed for each node. Trend
lines of stress transfer for axial end bearing and side resistance are provided for the following materials:
1.
2.
3.
4.
Sand
Clay
Intermediate Geomaterial
User Defined
Sand
Eqn: 12.5.G
The immediate settlements are computed using non-linear Q-z springs, with the shape presented in figure
shown below. The equation is provided but is should be referred that there is a considerable scatter around the
trend line.
Eqn: 12.5.H
Figure: 12.5.d Trend Lines for Drilled Shaft End Bearings in Sand
Clay
Eqn: 12.5.J
where cub = average undrained shear strength of the clay (computed 1 to 2 diameters below the shaft)
for Bb > 75 in (1.90 m)
Eqn: 12.5.K
Eqn: 12.5.L
Eqn: 12.5.M
for R 6.5
Eqn: 12.5.N
Figure: 12.5.e Trend Lines for Drilled Shaft End Bearings in Clay
Intermediate Geomaterial
The design of drilled shafts founded in intermediate Geomaterials is directly from FHWA's Load Transfer for
Drilled Shafts in Intermediate Geomaterials.
Note:
Types 1 and 2 are considered to be cohesive materials with an undrained strength, qu in the range of 0.5 to
5.0 Mpa.
Type 3 is primarily cohesionless and has Nspt from 50 to 100
Method 1 proposed by FHWA's Load Transfer for Drilled Shafts in Intermediate Geomaterials, for Type 1 and 2
materials has been coded herein.
l
Valid for IGM Type 1 and 2; 0.5 < qu < 5.0 Mpa; Recovery > 50 %;
Appropriate for very short sockets (L/D <2) or very long sockets (L/D>20);
Where there is strong layering in the formation, or where part of the socket is artificially roughened and
part is smooth
Required Data:
l
Number of Layers
qu (Mpa)
Mass Modulus - Em
Thickness
Em/Ei
Em/Ei
(percent)
100
70
(closed joints)
1.00
0.70
(open joints)
0.60
0.10
50
20
0.15
0.05
0.10
0.05
Note: Values of Em/Ei for RQD values between those shown can be estimated by linear interpolation on RQD.
User can select the "Custom Q-Z" curve option for applying user-defined Q-Z curves to subsurface pile nodes that
fall within a given soil layer.
The user defined curves are specified by a set of TEN points. Please enter the data set in two columns.
User defined Q-Z data
Z1 Q1
Z2 Q2
Z3 Q3
... ...
Z10 Q10
Where
Ziis the ith Z (axial deflection of pile tip) value on the user specified Q-Z curve.
Qiis the ith Q (end bearing) value on the user specified Q-Z curve.
The torsional stiffness of a pile embedded in soil is modeled using T- springs, where T is the torque applied to
the pile and is the angle of twist, in radians. The springs are located at the nodal points. T- springs can be
modeled as follows:
1. Hyperbolic Curve
2. Custom T-
The non-linear T- behavior of the soil is modeled using an hyperbolic curve, with initial slope as a function of
the shear modulus G. The ultimate value is based on the ultimate shear stress at the contact pile/soil.
Eqn: 12.6.A
where:
r = radius of the pile
o
= shear stress along L
o
For a long rigid pile embedded in a soil with shear modulus G, Randolph (1981) deduced the expression for the
torque per unit length
Eqn: 12.6.B
This expression does not consider the pile tip stiffness. For a long pile the tip contribution may be considered
negligible.
Eqn: 12.6.C
Eqn: 12.6.D
Eqn: 12.6.E
The ultimate shear stress can be obtained with the same procedures as for axial skin friction. As for the initial
shear modulus, it should be determined from in-situ tests.
12.6.2 Custom T-
User can select the "Custom T-" curve option for applying user-defined T- curves to subsurface pile nodes that
fall within a given soil layer.
The user defined curves are specified by a set of TEN points. Please enter the data set in two columns.
User defined T- data
1 T1
2 T2
3 T3
... ...
10 T10
Where
iis the ith (pile rotation) value on the user specified T- curve.
Tiis the ith T (axial stress) value on the user specified T- curve.
Soil Resistance Due to Pile Rotation is only applicable to "Limestone (McVay) 2-3 Option" soil type. This option
is used for the program to calculate and apply rotational springs to the pile nodes in the ground. These springs
are based on the axial resistance of the piles (skin friction) as well as the rotation of the piles. It is particularly
important in soil layers where the piles can develop large values of skin friction.
Eqn: 12.7.A
where: 1 and 2 are the values of the strain on the opposite sides of the shaft.
Figure: 12.7.a shows in detail how the bending strains are obtained from the measured strains.
Soils Lateral Resistance P(F/L) Form Bending Moments and Skin Friction
The difference in the moment at two different elevations is caused by soils lateral (P force/length) and axial force
(T force/length) resistance at the soil-shaft interface. The contribution to moment in the case of the latter is a
function of shaft diameter, and the soils T-Z curve as well as the rotation of the shaft. Shown in the figure below
are the forces acting on a buried element of length dz.
Eqn: 12.7.B
Consequently, from lateral force equilibrium, Figure: 12.7.b , the soil lateral P (force/length) is found as
Eqn: 12.7.C
If the side shear, T ( Figure: 12.7.b ), is taken into account, then moment equilibrium results:
Eqn: 12.7.D
OR
Eqn: 12.7.E
Eqn: 12.7.F
Substituting Eqn: 12.7.E into Eqn: 12.7.F , then the soil lateral resistance, P, is obtained:
Eqn: 12.7.G
Evident from Eqn: 12.7.G vs. Eqn: 12.7.C , the side shear on the buried element will act to shed load away from
the lateral resistance, P. The moment/unit length, Ms, of the side shear is obtained from the T-Z curve for the soil.
The value of T requires the displacement, Z, at a point on the shaft.
Eqn: 12.7.H
The moment per unit shaft length about O, M , is found by multiplying T by the distance to the cross-section
s1
1
centroid, R , as
1
Eqn: 12.7.I
The total moment per unit length may be found by summing the moments acting on all the slices:
Eqn: 12.7.J
where: n = number of slices
1. Membrane Element
2. Plate Element
3. Flat Shell Element
4. Mindlin Theory
5. Generalized Stress and Strain
6. Special Element for FB-MultiPier
7. Mesh Correctness and Convergence
There are two common versions of plate theory used in finite elements: Kirchoff and Mindlin. Kirchoff plate
bending theory is derived in a similar fashion to beam bending but includes bending in both directions. The
derivation assumes that the normal displacement, vertical displacement w, controls. In Kirchoff theory the
rotation, , in the plate is the derivative of w. This is the same as beam theory. This means that shear
deformations are ignored. In Mindlin theory, shear is included and the rotation is the sum of the derivative of w
and the shear angle. FB-MultiPier uses a Mindlin formulation.
The results from all plate elements consist of moments. Some plate elements also give the transverse shear, Q,
as a result. It is important to note that the moments and shear results are per unit length of plate. The following
figure ( Figure: 13.2.a ) gives the sign convention for moment and shear results.
Figure: 13.3.a Flat Shell Element as a Combination of Membrane and Plate Element
The flat shell element can be used to model structures where both bending and stretching effects need to be
considered. Many small flat shell elements can be used to form curved surfaces. The modeling of bridge decks,
wide flange beams and curved shell structures are three such structures where flat shell elements are commonly
used.
Eqn: 13.5.A
Where E* is a modified constitutive matrix. Notice that this is just like the equation for stress and strain except we
have moments and curvature. In plates, the displacement unknowns are the normal displacement and the two
rotations. Following the analogy of generalized stress, moments are equivalent to stress and curvature is
equivalent to strain. This means when using these elements in modeling, we treat the moment gradient like we
would stress to determine the level of shape function and number of elements required for an accurate analysis.
In addition, the difference in moment at a common node between two elements indicates the adequateness of
the mesh.
Notice that for the four - four node elements, the difference between the elements is 28%. This large percentage
error indicates a poor mesh (or not enough elements). Looking at the two - nine node model we see a closer
difference. Here the error is 14.0%. This indicates that the mesh is marginal but probably sufficient. Finally we
look at the 40 element model. Here the error is much better and only 3%. The 40 element model is very good.
The difference in element stresses at a node is an important measure of model correctness. In general, we
do not have the exact displacements in order to check our model. Hence, the stress check is necessary to
verify convergence of our model. If the difference in stresses between elements is small the finite element
mesh is good.
14 Nonlinear Behavior
A discrete element formulation is used to model the nonlinear behavior of piles, pier columns, and pier cap
members in FB-MultiPier. The discrete element accounts for both constitutive (i.e., material) and kinematic (i.e.,
geometric) nonlinear behavior. The nonlinear material behavior is modeled by integrating stress-strain curves
over the cross-section that is applicable to each element. The nonlinear geometric behavior is modeled using Pdelta moments (the product of the axial force and relative displacements that occur across a given element).
Additionally, because pile, pier column, and pier cap members are discretized into a number of sub-elements, Py moments (the product of axial force and bending-induced internal displacements within a given element) are
also accounted for in the discrete element formulation.
Salient aspects of the discrete element model are discussed in the links provided below.
Eqn: 14.1.A
Eqn: 14.1.B
The elongation of the center section of the element is calculated as follows:
Eqn: 14.1.C
The angle changes for the center section about the z and y-axes are then defined below:
Eqn: 14.1.D
Eqn: 14.1.E
The discretized vertical and horizontal angle changes at the two universal joints are then:
;
Eqn: 14.1.F
;
Eqn: 14.1.G
and the twist in the center part of the element is defined as:
Eqn: 14.1.H
Thus, the internal deformations of the discrete element model are uniquely defined for any combination of
element end displacements.
The curvature for small displacements at the left and right universal joints about the y and the z axes are defined
as follow:
At the left joint,
;
Eqn: 14.1.I
At the right joint,
;
Eqn: 14.1.J
The axial strain at the center of the section is given by:
Eqn: 14.1.K
Eqn: 14.1.L
And, to satisfy equilibrium:
Eqn: 14.1.M
Eqn: 14.1.N
Eqn: 14.1.O
The relationship for strain at any point in the cross-section is:
Eqn: 14.1.P
The stress at any location in the section is found using the appropriate material stress-strain curve described
subsequently.
Numerical integration of equations is done using Gaussian Quadrature. To use the method of Gaussian
Quadrature, the function being integrated must be evaluated at those points specified by the position factors.
These values are then multiplied by the appropriate weighting factors and the products accumulated. Figure:
14.1.c shows a square section with 25 integration points (a 5x5 mesh). The number of defaults integration points
for square pile is set at 49 (a 7 by 7 mesh). For circular sections, the section is divided into circular sections (12
radial divisions and 5 circumferential divisions as shown in Figure: 14.1.d ). The sections are integrated at the
centroid of each sector using weighting factors of 1.0. The stress in all steel bars is evaluated at the centroid and
a weighting factor of 1 is used for each bar.
Figure: 14.1.c Section Integration Divisions - Cross Section of square pile showing integration points
Figure: 14.1.d Section Integration Divisions - Circular pile cross section showing steel rebars
When a circular void is encountered in a square section, the force is first computed on the un-voided section and
then the force that would be acting on the voided circular area is computed and subtracted from the force
computed for the non-void section. Circular sections with voids are divided into sectors omitting the voided
portion.
Even for nonlinear material analysis, the torsional moment M5 is assumed to be a linear function of the angle of
twist, 5 , and the torsional stiffness GJ, where J is the torsional constant and G is the shear modulus as shown
next
Eqn: 14.1.Q
Eqn: 14.1.R
Eqn: 14.1.S
And from equilibrium of the end bars:
;
Eqn: 14.1.T
;
Eqn: 14.1.U
;
Eqn: 14.1.V
;
Eqn: 14.1.W
where f1- f3 and f7 - f9 are the acting end forces, and f4 - f6 and f10 - f 12 are the end moments.
Note: Linear interpolation is used to calculate stresses for strain levels beyond those provided in the stressstrain curve.
14.2.1 Concrete
The figure below shows the default stress-strain curve used by FB-MultiPier in determining internal demands for
concrete portions of member sections in nonlinear analyses. The function values are dependent upon f'c and Ec,
as input by the user. Note that this stress-strain model incorporates tension stiffening, and is applicable for
unconfined normal weight concrete with compressive strengths of approximately 41 MPa (6,000 psi) or less
(MacGregor and Wight 2005).
The tensile portion of the curve is assumed as linear up to a stress of fr. Piece-wise linear tension softening is
modeled out to a tensile strain of 0.002, and no tensile stress contribution is modeled for tensile strains greater
than 0.002. The tension softening portion of the stress-strain curve is intended to account for the uncracked
portions of in-tension concrete sections, which are capable of carrying some non-zero level of tensile stress. The
value of fr is based on the fixed value of r shown in the figure below, and the modulus of elasticity Ec input by the
user. For English units this will give a value of fr of 7.5 f'c.
The compression portion of the concrete curve is highly non-linear, and is defined by the Modified Hognestad
parabola. For compressive strains bounded between the strain at peak compressive stress (0 ) and the ultimate
compressive strain (cu ), a linear descending curve segment is employed. For any members that generate
compressive strain levels exceeding cu , the stress level is estimated using linear extrapolation. For instances
where such large compressive strain levels are reached (or for other instances where the default curve may not
be applicable), users are recommended to supply a user-defined stress-strain curve. Additionally, confined
concrete model is available.
Eqn: 14.2.A
The default relations for the mild steel stress-strain curve are given by,
Eqn: 14.2.B
Eqn: 14.2.C
Eqn: 14.2.D
The default stress-strain curve generated for steel with f' =60 ksi and E =29600 ksi is shown in the figure below.
y
c
Figure: 14.2.c Prestressing Steel Stress-strain Curve for fsu = 270 ksi.
Eqn: 14.3.A
where
f
cc
x is given by:
Eqn: 14.3.B
The expression suggested for
cc
cc
Eqn: 14.3.C
where
f
co
co
Eqn: 14.3.D
is the tangent modulus of elasticity for unconfined concrete and is given by:
Eqn: 14.3.E
E
sec
cc
and
cc
Eqn: 14.3.F
For f , the confined concrete strength, Mander used the five-parameter failure criterion proposed by William
cc
and Warnke and the tri-axial test data of Schickert and Winkler. In the case of circular columns confined by
circular hoops or spirals, the confined concrete compressive stress has been shown to be:
Eqn: 14.3.G
where
f is effective confining pressure, and may be obtained from the equilibrium of internal forces acting on the
l
dissected sections shown in Figure: 14.3.b
For the cover concrete in columns, assuming uniform yield of the jacket, the equilibrium of forces requires:
Eqn: 14.3.H
where
f
Figure: 14.3.b Confining Action of Steel Jacket and Internal Hoops [4]
Eqn: 14.3.I
Substituting into Eqn: 14.3.H we obtain
Eqn: 14.3.J
By using f = f in Eqn: 14.3.G , the compressive strength of the cover concrete confined by the steel jacket
l
lj
can be determined.
Additional confinement is provided to the concrete core by the transverse reinforcement. The additional lateral
pressure, f , may also be determined from the equilibrium of forces. Assuming uniform yield of the transverse
lh
steel yields the following equation:
Eqn: 14.3.K
where
d
is the diameter of the concrete core defined along the center line of the confining steel
Eqn: 14.3.L
where
A is the area of an effectively confined concrete core
e
Eqn: 14.3.M
where
A
cc
Eqn: 14.3.N
where
A is the total longitudinal steel area .
s
By assuming an arching action between circular hoops in the form of a second -degree parabola with an initial
tangent slope of 45E, the confinement effectiveness ratio has been shown to be:
Eqn: 14.3.O
where
s is the clear distance between the hoop.
Similarly, the confinement effectiveness coefficient for a circular spiral has been shown to be:
Eqn: 14.3.P
By introducing as the ratio of the volume of transverse confining steel to the volume of confined concrete i.e.:
s
Eqn: 14.3.Q
Eqn: 14.3.R
The lateral confining pressure due to transverse steel in Eqn: 14.3.J may be written as:
Eqn: 14.3.S
Thus using f =f + f in Eqn: 14.3.G will allow the enhanced compressive strength of the concrete core to be
l
lj
lh
determined.
Scott et al (1989) proposed an expression for the ultimate compressive strain, , which is given by:
cu
Eqn: 14.3.T
where
Eqn: 14.3.U
Eqn: 14.3.V
Eqn: 14.3.W
Eqn: 14.3.X
It is assumed that the stress-strain curve for unconfined concrete follows Eqn: 14.3.Y during the earlier stages of
loading up to 2 . For compressive strains larger than 2 , the strains are assumed to decrease linearly with
co
co
strains up to the spalling strain . A value of 0.005 has been adopted for . The longitudinal compressive
sp
sp
stress for unconfined concrete may be written as:
For <= 2 ,
c
co
Eqn: 14.3.Y
For 2
co
< <= ,
c
sp
Eqn: 14.3.Z
For >
c sp
Eqn: 14.3.AA
14.3.3 Reinforcement
To avoid congestion of reinforcement, earlier design practices tended to use large diameter bars, up to #14 or
#18, however, such practice may lead to potential bond problems in cases where the column main reinforcement
were lapped at insufficient length with starter bars in the plastic hinge regions. Consequently, such columns are
characterized by very rapid flexural strength degradation under the design seismic loads. The current Caltrans
(1981) approach has been to avoid lap splicing of the main reinforcement in the potential plastic hinge region of
bridge columns. The analytical model developed here assumes full yield of the main reinforcement including
strain hardening.
The monotonic uniaxial stress-strain curve of a typical reinforcing steel is shown by an elastic region, a yield
plateau, a strain hardening region, followed by a falling branch after peak stress up to the strain at which fracture
occurs. A typical stress-strain curve for the reinforcing steel is shown in Figure: 14.3.d .
The monotonic uniaxial stress strain curve for reinforcing steel is defined by the following equations:
For the elastic range, i.e. <=
s
y
Eqn: 14.3.AB
Where
is the axial strain in the reinforcing steel
s
f is the stress in the reinforcing steel
s
E is the modulus of elasticity of the reinforcing steel
s
For the yield plateau, i.e. < < ,
y s sh
Eqn: 14.3.AC
where
Eqn: 14.3.AD
where
is the ultimate strain in the reinforcing steel
su
f
is the ultimate stress in the reinforcing steel and
su
Eqn: 14.3.AE
Eqn: 14.3.AF
It has been shown by Mizra and MacGregor (1979) that the ratio of ultimate to yield strength was f /f = 1.55.
su y
The steel model adopted for the program assumes a modulus of elasticity of 29000 ksi and a slightly lower
ultimate to yield strength ratio of 1.50. The other mechanical properties assumed for the stress strain model are:
Eqn: 14.3.AG
Eqn: 14.3.AH
where
f is equal to 40 for ksi units. This would be converted to any other consistent set of units
yl
The above equations are non-dimensional, allowing the model to be used with any grade steel. They were
obtained by interpolating from the values given by Priestly for 40 and 60 ksi steel.
It should be noted that the tangent modulus at the onset of strain hardening may be obtained by taking the
derivative of Eqn: 14.3.AD with respect to steel strain, and operated at the strain-hardening strain, :
s
sh
Eqn: 14.3.AI
60"
30'
4"
5.2 ksi
25 #14
Yield Strength fy
Transverse Steel
68.9 ksi
#5 Spiral at 3.5"
71.5 ksi
1000 kips
The force deformation curve for the full-scale column is given in the figure above. As can be seen from the figure,
the data from FB-MultiPier program is generally close to both the COLRET values and the experimental data. For
the majority of the curve the FB-MultiPier values are less than the COLRET values. Also, it is noted that the initial
stiffness of the response is higher from the CORLET than obtained from the FB-MultiPier Analysis and the
measured response.
24"
12'
4"
5.2 ksi
26 #6
Yield Strength fy
Transverse Steel
45.7 ksi
#2 hoops at 5 in.
51.0 ksi
48"
.188"
47 ksi
400 kips
The force-deformation curve for the jacketed column is given in the figure above. Looking at the results, we can
see that FB-MultiPier provides a close estimation of the experimental and COLRET curves until the post yield
region of curve where we see a reduction in the lateral load capacity predicted by FB-MultiPier in comparison to
the experimental and COLRET values. It is also noted that the CORLET program show slightly greater strengths
than that for the test.
14.3.9 Conclusions
A model for the prediction of the non-linear response of circular concrete piles with confinement has been
presented. More details on the model are available in Stone and Cheok (1989) .
This model has been incorporated into the FB-MultiPier computer program that is used specifically for analyzing
bridge pier structures consisting of pier columns and cap supported on piles or shafts. This allows the user of the
program to model the behavior of concrete piles confined by hoops, spirals and/or a steel jacket subjected to a
broad variety of loadings.
In the comparative studies conducted, the FB-MultiPier results show generally less of an increase in strength and
ductility than those given by the COLRET program. This is due to the following differences between the FBMultiPier program and the CORLET program.
First, Eqn: 14.3.T is used to compute the maximum concrete strain, gives less strain than the procedure
cu
using COLRET. COLRET uses a more complex procedure that was only documented for grade 40 and grade 60
steel. The FB-MultiPier program is written to handle a wide variety of inputs and thus it used the more
conservative Eqn: 14.3.T which is applicable for any grade of steel.
Second, the COLRET program assumes the entire area contained within a diameter ds is confined in integrating
the stresses over the column area, whereas FB-MultiPier conservatively uses only the effectively confined area
of the core.
Finally, in the case of an external steel jacket, FB-MultiPier neglects the longitudinal stiffness of the jacket, when
using the confined model and the COLRET program takes this stiffness into account.
These differences tend to give a somewhat conservative solution, which is probably best for a general purpose
program to be used for a wide variety of applications. The program can of course be modified to accommodate
more detailed models in the future.
diagrams and linear analysis of the piles and structure is provided for preliminary design phases and must not
be used for a final design.
Also it should be noted that the nonlinear analysis procedures while very thorough may not reflect all of the
criteria required for design of slender structures in appropriate codes. A list of factors that may affect the final
design that are not considered in the analysis are as follow:
1) Creep in concrete.
2) Initial imperfections or out of straightness of members.
3) Residual stresses.
4) A separate analysis phi factor to account for the possibility of undersized members.
All of the factors could potentially increase the defections and thus the P and Py moments above those given by
the nonlinear analysis routines and should be accounted for by the designer as appropriate.
The routine operates by computing numerous horizontal slices of the P , M , M (nominal strength) interaction
n nx ny
surface. The result is a series of M , M interaction curves for various magnitudes of axial load, P . Next, each
nx ny
n
interaction curve is represented in the form.
Eqn: 14.4.A
The moments M and M represent the nominal moment strength at axial load, P , for uni-axial bending about
0z
0y
n
the z and y axes respectively. The exponents, and are computed in by the least squares method. They
enable the above expression to fit the computed interaction curve, and vary with axial load. The actual DC is
computed, by interpolation for the axial force P / using the stored values of axial force and exponents, and .
u
The program then computes the parameter, DC, as; the ratio of the length of the vector for the actual forces (in 3d space M , M and P) divided by the length of the vector with the same direction as the actual results but
ux uy
with that vector touching the 3-D failure surface. This has the effect of assuming that the moments and axial load
will increase proportionally until failure.
Note that the routines handle square, rectangular or circular sections with prestressed or non-prestressed steel
and H-piles that themselves may be encased in concrete. The capabilities and limitations of sections were
discussed in detail earlier in the User Guide. The interaction diagrams for all sections that contain concrete are
handled in a similar manner and will be discussed next. Then the case of the H-pile section without concrete will
be discussed last.
Sections with reinforced or prestressed concrete.
The routine assumes a planar strain distribution across the section. The criterion for section failure is that the
concrete reaches the crushing strain cu. ( = -0.003 in/in) at one corner of the section. This conservatively
cu
ignores any effect of the ties or spiral reinforcement on the compressive strength or crushing strain of the
concrete.
All tensile stresses in the concrete are neglected. This includes both tension in uncracked regions and tension
stiffening in cracked regions.
AASHTO doesnt permit the design of a perfectly axially loaded column. A certain minimum eccentricity of load
must always be included. This is accomplished by limiting the applied factored axial force, P , to a factor times
u
P , where P is the nominal capacity of the section in an axially loaded column. For a tied column, this
0
0
maximum load is 0.8P , while for a column with spiral reinforcement the maximum load is 0.85P . When a
0
0
factored load, Pu , larger than these limits is input to the routine, the routine responds that the section is
inadequate.
The routine also includes a maximum axial tension for a section based on all the mild steel bars attaining the
yield stress fy for mild steel and all the prestressing strands attaining the ultimate stress fsu . If the factored Pu
exceeds this maximum tension force the routine responds that the section is inadequate.
Solutions for and for approximately 30 points between the maximum compressive and tensile forces are
attempted. The values of axial force are obtained for strains in the extreme bar farthest away from the corner of
the section where limiting (tensile and compressive) reinforcement (mild and prestressed) strains are taken in
accordance with those values given in Sec. 5.7.2.1 of the AASHTO LRFD code
The strength reduction factor, , is determined according to the unified requirements for prestressed and
nonprestressed concrete of AASHTO - LRFD. In these requirements, the magnitude of is based on the net
tensile strain occurring in the most heavily strained steel bar or strand when the nominal strength of the section is
attained (when the concrete crushes). The net tensile strain is that portion of the steel strain associated with the
development of tensile strain in the concrete adjacent to the steel bar or strand. For nonprestressed steel, the net
tensile strain is exactly the total strain in the steel. For prestressed steel, the net tensile strain is the total strain in
the steel minus the sum of the effective prestress strain in the steel and the effective prestress strain in the
concrete adjacent to the steel. The latter term can be thought of as the decompression strain.
The factor can be obtained for net tensile strain in the extreme tension steel in nonprestressed reinforcement
by linearly interpolating between 0.75 (at compression controlled reinforcement strain limit ) to 0.90 (at tension
cl
controlled strain limit ). Whereas the factor for prestressed reinforcement can be obtained for net tensile
tl
strain in extreme steel by linearly interpolating between 0.75 (at compression controlled strain limit ) to 1.00 (at
cl
tension controlled strain limit ). (Refer to Sec. 5.5.4.2.1 in the AASHTO LRFD code).
tl
H-piles embedded in concrete are treated as if fully bonded to the concrete and are thus treated just as if they
were an equivalent group of a large number of small rectangular reinforcing bars.
H-piles without Concrete
H-Piles not embedded in concrete have their interaction diagrams computed in a similar manner to the sections
with concrete except that steel does not have a small limiting compression strain due to its ductility. Thus the 30
points defining the range of axial forces are obtained by locating the neutral axis in successive positions across
the depth of the section and the steel is assume fully yielded on both sides of the neutral axis. Solutions are
attempted for and to represent the interaction curves for 30 Points with varying axial load from a peak in
tension to compression.
Steel Only Sections
The interaction diagram of steel sections such as H-pile and pipe pile sections is based on development of full
plastic section. Steel code requirements must be checked for design purposes.
A novice user may then be tempted to say that one should stick to linear programs and avoid such difficulty.
However, the counter argument can be made that a linear analysis will almost always find a solution even if the
user puts in a totally unreasonable loading.
For the sake of discussion, assume that a relatively simple structure is being modeled by FB-MULTIPIER,
perhaps even a single pile cap with one or two piles with some vertical load applied which is held constant and
then a lateral load is applied gradually. Several different scenarios of lateral load versus lateral displacement are
possible as shown in the following figure.
Suppose the designer wants to demonstrate that the behavior is indeed type 1 versus type 2. A push over
analysis could be done and this requires a displacement-controlled solution. A large spring would be placed at
the node where the lateral load is applied and then a series of large loads would be applied. The spring would
take the larger amount of the load but by properly choosing the spring stiffness and load, the displacements
could be controlled and the load absorbed by the structure could be found and the pushover results plotted.
In rare instances the response of a structure may resemble that shown for case 3. Here at a load near d the curve
flattens and may even decrease. However, for increasingly large displacements the load may start to rise again.
It will be very difficult to obtain converged solutions for loads near d. However, if a much larger load is applied a
solution may be found on the curve well above d. This type of behavior generally occurs when some type of local
failure occurs. If the structure has sufficient ductility it may then be able to find a new path to distribute the forces
and carry some additional load, albeit with a considerable reduction in stiffness. An example of this type of
behavior is when the gravity loading is small and because of a large lateral load a pull out occurs on one of the
piles. The question then arises, should the design based on the post pull out behavior be used?
Clearly the use of nonlinear analysis programs do not remove the responsibility of the designer to monitor the
local responses of the structure. Fortunately the program outputs detailed information about the behavior of the
soil and pile that can and must be reviewed before a structure can be said to be adequate.
Finally, case 4 in which the structure appears to move against the loads must be considered. For very slender
structures with very large gravity loading, the stiffness of the structure will go negative when the elastic buckling
loading of the structure is exceeded. Again this is a rare case and would almost never happen for a designer
evaluating a real structure. However, someone trying the program out with arbitrary dimensions and loads might
create such a condition and then be disturbed that the program is giving obvious unreasonable results. A linear
analysis program would of course produce even more possibly dangerous results; it would indicate a positive
displacement, which would then not give any indication that something was wrong with the structure.
15 Dynamic Behavior
The dynamic behavior used in FB-MultiPier is explained in detail in the following topics:
1. Analysis Type
2. Time Stepping Parameters
3. Dynamic Relaxation
4. Damping
5. Global Mass
6. Time Functions
7. Dynamic Soil Properties
8. Structural Material Behavior
Eqn: 15.1.A
where
velocity vector;
vector.
is the nodal
For a given structural system, the solution to Eqn: 15.1.A is necessary to determine structural response to
external loading
. However, for systems that possess kinematic nonlinearity; constitutive nonlinearity;
many degrees-of-freedom (DOF); or any combination thereof, the obtainment of a direct analytical solution to
Eqn: 1.1.A is not feasible. In place of a direct analytical solution, FB-MultiPier employs implicit time integration
algorithms to obtain a numerical solution of Eqn: 15.1.A . Specifically, in FB-MultiPier, both the Newmark and
Wilson-Theta methods are available for finding the numerical solution to Eqn: 1.1.A at user-specified, discrete
time increments (i.e., time-steps, t). See Fernandez (1999) for details regarding the theory and implementation
of these methods in FB-MultiPier.
Determination of time-step t:
An important property of the implicit dynamic time-integration methods in FB-MultiPier is that these methods are
implemented so as to be unconditionally stable: for a given structural system, the numerical solution will not
diverge solely due to the selection of a very large time-step, t, for the analysis. It is important to note, however,
that other factors (e.g., loading a nonlinear system beyond the system capacity) can still result in a lack of
convergence for a given analysis. Furthermore, numerical stability in the implemented time-integration methods
does not guarantee a numerically accurate solution. The following expression is recommended in Tedesco et al.
(1999) for determining an appropriate time-step for a given implicit time-integration analysis, with consideration
of accuracy:
Eqn: 15.1.B
where
Transient loading
Transient loading can be specified as an arbitrary set of nodal load-histories,
, in FB-MultiPier. Such
loading is readily handled by the time-integration methods implemented in FB-MultiPier, given that a proper
time-step t has been selected.
The user may define their own stress strain curves for concrete and steel or use the default values described
below.
1.
2.
Concrete (Hysteretic)
3.
4.
Note: Linear interpolation is used to calculate stresses for strain levels beyond those provided in the stressstrain curve.
Overview
For nonlinear dynamic (time-history) analysis, frame element cross sections are divided into fibers for the
purposes of integrating stresses and calculating element internal forces. The default constitutive behavior for
concrete portions of cross-sections in nonlinear dynamic analyses is that of nonlinear elastic behavior.
Alternatively stated, the default approach in mapping concrete strains to stresses is such that no deviation occurs
from the backbone curve. While this approach increases the likelihood of convergence for a model, no
hysteresis occurs within concrete fibers.
Stress-strain curve description
The default stress-strain curves for concrete portions of structural members are described here (for unconfined
concrete cross-sections) and here (for confined concrete sections ). Alternatively, FB-MultiPier permits userdefined input of stress-strain curves for concrete portions of frame element cross sections. In this latter case, the
user-defined stress-strain curve is treated as the backbone curve (but with no deviation from the backbone
curve).
Overview
For nonlinear dynamic (time-history) analysis, frame element cross sections are divided into fibers for the
purposes of integrating stresses and calculating element internal forces. Hysteretic behavior can be incorporated
into the constitutive behavior of concrete portions of cross-sections in nonlinear dynamic analyses. Alternatively
stated, an option is available to introduce behaviors such as unloading, reloading, and softening of concrete
fibers while conducting time-history analyses. Incorporating hysteresis increases model fidelity with respect to
that of physical systems containing concrete members. However, the incorporation of concrete hysteresis
substantially constrains the circumstances under which convergence can be reached. Specifically, because FBMultiPier utilizes a tangent stiffness formulation for time-history analysis, convergence is not likely to be reached
in the event that a given system is force-driven AND, for said system, several fibers within any one frame element
cross-section develop strains corresponding to (or exceeding) a state of material softening. Note that if a
displacement-driven process is modeled, then it is possible to reach convergence despite reaching a state of
material softening among the computed response data.
Stress-strain curve description
The default stress-strain backbone curves for concrete portions of structural members are described here (for
unconfined concrete cross-sections) and here (for confined concrete sections ). Alternatively, FB-MultiPier
permits user-defined input of stress-strain curves for concrete portions of frame element cross sections. In this
latter case, the user-defined stress-strain curve is treated as the backbone curve. The algorithm governing
deviations from the backbone stress-strain curve is adopted from Yassin (1994) .
Accessing the feature in the User Interface (UI)
The concrete hysteresis option may be activated by navigating to the Advanced dialog, located on the
Dynamics page. From within the Advanced dialog, check the Concrete Hysteresis option. Note that this
option, when activated, applies to all nonlinear frame elements in the model.
Overview
For nonlinear dynamic (time-history) analysis, frame element cross sections are divided into fibers for the
purposes of integrating stresses and calculating element internal forces. The default constitutive behavior for
steel portions of cross-sections in nonlinear dynamic analyses is that of a bilinear approach to modeling
hysteresis. Alternatively stated, the default approach in mapping steel strains to stresses is such that deviations
from the input stress-strain curve (i.e., the backbone curve) adopt a slope that is equal to the initial slope of the
backbone curve. While this approach mimics increases the likelihood of convergence for a model, no hysteresis
occurs for concrete fibers.
Stress-strain curve description
The default stress-strain curves for longitudinal reinforcement within structural members are described here (for
mild steel bars) and here (for prestressed tendons). Alternatively, FB-MultiPier permits user-defined input of
stress-strain curves for steel portions of frame element cross sections. In this latter case, the user-defined stressstrain curve is treated as the backbone curve (where all deviations adhere to the initial tangent stiffness of the
backbone curve).
Overview
For nonlinear dynamic (time-history) analysis, frame element cross sections are divided into fibers for the
purposes of integrating stresses and calculating element internal forces. The option to incorporate relatively
more advanced hysteresis behavior into steel portions of cross-sections is available for nonlinear dynamic
analyses. Specifically, in contrast to the (default) Steel-Bilinear hysteresis modeling approach, phenomena such
as strain-hardening and the Bauschinger effect can be included among the simulated steel unloading, reloading,
and material evolution behaviors.
Stress-strain curve description
The stress-strain backbone curves for steel portions of structural members are described in Filippou et al.
(1983) , where only the Youngs modulus and yield stress need to be specified. Alternatively, FB-MultiPier
permits user-defined input of stress-strain curves for steel portions of frame element cross sections. In this latter
case, parameters that satisfy the curve given in Filippou et al. (1983) are derived from the user-defined stressstrain curve.
Accessing the feature in the User Interface (UI)
This option may be activated by navigating to the Advanced dialog, located on the Dynamics page. From
within the Advanced dialog, check the Steel Hysteresis (Strain-Hardening and Bauschinger Effect) option.
Note that this option, when activated, applies to all nonlinear frame elements in the model.
2.
3.
Wilson-Theta.
See Tedesco et al. (1999) for details regarding each technique. For a selected time step integration, the time
step size (sec) and the number of time-steps must additionally be input.
analysis, the structure achieves static equilibrium between the product of stiffness
response to the application of
and displacements
in
permanent loads,
, is necessary.
Dynamically, an abrupt application of self weight loading is unrealistic, and furthermore, can lead to undesired
dynamic response. In Figure: 15.4.a , an example of improper dynamic self-weight load (w) application is given
to demonstrate that instantaneous application of permanent loads can lead to exaggerated structural response.
Inertial forces caused by motion of the single-degree of freedom (SDF) mass push the system far beyond the
static displacement, ust. By applying the load w abruptly, the maximum dynamic displacement, umax, is 100%
larger than that predicted by a static analysis. Additionally, the excessive displacement leads to amplified
internal forces,
), vertical
dynamic structural response. The erratic pier behavior caused by instantaneous application of
renders any
assessment of structural responsesuch as determining maximum internal forces or load-moment interaction
unreliable.
The dynamic relaxation feature is used in FB-MultiPier to facilitate staged loading in transient dynamic analysis
by initializing the dynamic system such that it is in equilibrium with an arbitrary set of permanent, or static, loads
(e.g., gravity loads),
initialized as:
. When dynamic relaxation is employed, the equation of motion (recall Eqn: 15.1.A ) is
Eqn: 15.4.A
where
is the stiffness
pre-analysis; and,
is the initial transient load vector. It is important to note that
to be constant throughout the analysis.
and
are assumed
As expressed in Eqn: 15.4.A and shown in Figure: 15.4.b , for dynamic relaxation, structural response to
is initialized by means of a static pre-analysis. Meaning, two distinct analyses are conductedone static analysis
(with only
and
quantities are then used to define the initial state for a time-domain analysis. In the dynamic analysis,
is
and
are initialized to
and
, and
quantities are present in the equation of motion after the first numerical solution iteration.
and displacements
Figure: 15.4.f Total soil reaction (at analysis termination time) per model and analysis type.
15.5 Damping
Note: For models utilizing the English unit system, and for instances where the RX, RY, or RZ degrees of
freedom are activated, then the units of the Load Function for these degrees of freedom is kip-in/sec.
Cycles defines the number of loading cycles needed to reduce the soil strength by 50%.
The number of loading cycles is dependent upon the soil type and system natural frequencies. For details, see
the figures from "Sudip Basack. 2007".
2. Load rate
This effect occurs with an increase in the apparent soil (axial or lateral) resistance caused by rapid rates of
loading that occur in most extreme events, which is approximated as:
Eqn: 15.8.A
where,
pi = instantaneous soil resistance,
pc = resistance considering only cyclic loading,
tr = actual rate of loading (hz)
ts = corresponding rate of loading for standard slow cyclic loading (typically 0.01 to 0.1 hz)
f2 = soil factor
sands = 0.01 - 0.03, silts = 0.02 - 0.07, clays = 0.02 - 0.12, calcareous = 0.01 - 0.03
3. Shear Wave Velocity
Shear Wave velocity is used in the calculation of "Radiation Damping".
4. Degradation Soil Factor:
The soil degradation factor is used to scale soil p-y curves to mimic degraded soil behavior. Soil degradation
factors can typically range from 0.5 to 0.75.
5. Soil Damping Factor
The Soil Damping Factor is used to directly impose damping behavior into subsurface nodes. At the end of any
time step numbering two or greater, the global damping matrix entries corresponding to degrees of freedom of
subsurface nodes are updated by:
A.
B.
C.
D.
6. Soil Gap
Users can incorporate soil gapping into dynamic analyses. Refer to Chapter 9 of "Consolazio et al. 2006" for
details about how this feature incorporates hysterisis into the soil behavior.
Cyclic degradation can be activated on the "Soil Page". During analysis, cyclic degradation is carried out
usings:
Eqn: 15.8.B
where,
pc = soil resistance for the current cycle of loading,
pp = value of the resistance corresponding to the present value of deflection, i.e. y of p-y curve, on the previous
cycle of loading,
pd = fully degraded value of p at the present value of y specified by the user in py-curve
c1 = soil degradation parameter
.
Eqn: 15.8.C
where,
n50 = number of cycles that would be required to degrade the soil by 50%.
Radiation Damping:
Radiation damping of embedded nodes is utilized within the FB-MultiPier engine when global damping is active
(by checking the Damping checkbox on the "Dynamics Page"), but the mass and stiffness (Rayleigh) damping
coefficients for the Soil are both defined as 0.0. Note that radiation damping is only imposed on horizontal
translations using:
C = 4d( / GF)SW
Eqn: 15.8.D
where,
d = Effective diameter of the pile section.
= Unit weight of the soil.
GF = Gravity factor= g (Scale Factor); g = 9.81 m/sec2 or 386.2 in/sec2 ; Scale Factor = user input in the
Dynamic Page.
SW = Shear wave velocity of soil.
The model (pile, pile-and-cap, single pier, or pile bent) is analyzed based on the applied loading.
2.
After the solution is obtained for the applied loading, the program calculates the flexibility matrix of the
structure at the particular equilibrium state, following general principles. To do that the program internally
applies unit forces at the selected node of interest. The unit forces are applied successively in all six
possible directions (Fx, Fy, Fz, Mx, My, Mz).
3.
The displacements from each solution at the selected node comprise the columns of the flexibility matrix
(e.g., the displacements from the solution under the application of the Fx unit load comprise the first
column of the flexibility matrix).
4.
After the flexibility matrix is obtained, the program calculates the matrix inverse, which constitutes the
equivalent stiffness at the node of interest.
5.
The program reports both the flexibility and equivalent stiffness matrices.
The calculated stiffness (or flexibility) matrix is calculated after the equilibrium state of the structure is obtained.
This is necessary since models generally contain nonlinear elements (including the soil springs). Therefore, the
snapshot in time (i.e., the equilibrium state) at which the stiffness is calculated is very important.
If the program was not following the particular sequence (not obtaining the equilibrium solution prior to
calculating the equivalent stiffness), then the calculation of the equivalent system stiffness would be incorrect
since it would be obtained using information for a state of the structure other than the loading-induced
equilibrium state.
The stiffness can be thought of as being the tangent stiffness (instead of secant) for the simple reason that it is
calculated for a particular instance in time (equilibrium state).
Given below is a step by step guide to obtain a stiffness matrix at any nodes selected in pile, pile and cap, single
pier, and pile bent models.
Step 1: Select new problem type and click on pile and cap model. This will open a default pile and cap model in
FB-MultiPier (See Figure: 16.1.a )
The following explanation shows how to convert a 6x6 stiffness matrix from the FB-MultiPier global coordinate
system to a standard coordinate system defined below.
A 3x3 transformation matrix (T) is first defined to show how the two coordinate systems are related.
Eqn: 16.2.A
This requires 2 matrix multiplications to obtain the transformed stiffness matrix. This can be easily done using
either Excel or MathCad.
As a result, to convert the FB-MultiPier stiffness to a standard coordinate system, use the following.
Note: Both the locations and signs change for some of the stiffness terms.
Example
The FB-MultiPier stiffness matrix is given by
Then the stiffness matrix in the standard coordinate system would be.
To transform the 6x6 stiffness matrix generated by FB-MultiPier from the 2D x-z coordinates system to a
standard 2D X-Y coordinate system.
The following transformation is used to transform the stiffness matrix from the FB-MultiPier coordinate system to
the standard 2D coordinate system as follows.
Eqn: 16.2.B
As a result, to convert the FB-MultiPier stiffness to a standard coordinate system, use the following.
The engine input parameters are divided into two categories. The Global headers category describes control
data for the entire model. The Pier Specific headers category describes data that is specific to a particular pier in
the model.
Pier specific data is separated in the input file by a special pier header as follows:
__PIER#x
Where
xPier number
The following two lines specify the data to be printed to the output file. These lines are always required.
PRINT
L=L1 M=M1 D=D1 O=O1 S=S1 P=P1 T=T1 F=F1 C=C1 B=B1 J=J1 H=H1 I=I1 R=R1 N=N1 X=X1 G=G1 A=A1
Z=Z1 E=E1 K=K1 Q=Q1 V=V1 W=W1
Any of the values: L1,M1,D1... etc. can be either 0 or 1. Setting a value to 1 enables its printing. Setting the value
to 0 turns off the printing of that data block. The default is =0 (NO print). Only the options desired (=1) are
required. A SUMMARY OUTPUT TABLE WILL ALWAYS BE PRINTED.
Where
L1is the flag for printing of the Pile and Structural coordinates.
M1is the flag for printing the pier material properties.
D1is the flag for printing the pile displacement.
O1is the flag for the out of balance forces.
S1is the flag for the soil response forces.
P1is the flag for the pile element forces.
T1is the flag for the pier columns and pier cap displacement.
F1is the flag for the pier columns and pier cap force output.
C1is the flag for the pile cap stress/moment output.
B1is the flag for the bridge simulation spring force output.
J1 is the flag for the Bride Span Element Displacement.
H1 is the flag for the Bridge element forces.
I1is the flag for printing the output of the interaction diagram.
R1is the flag for printing section stress-strain data.
N1is the flag for printing missing pile information.
X1is the flag for XML data printing. This is used in conjunction with the Model Data Report Generator for data
extraction.
G1 is the flag for the Pile Material Properties
A1 is the flag for the Soil data per layer
Z1is the dynamics printing option.
E1 is the flag for the Soil data per pile node;
K1 is the flag for the soil graph per pile node
0 = print displacement, velocity, and acceleration results to a binary file
1 = print displacement, velocity, and acceleration results to an ASCII file
Q1is the flag for the printing span properties.
V1 is the flag for the to Stop Analysis after Eigenvectors are printed.
W1 is the flag for printing span temperature input data and loading.
The following lines specify the control parameters for the FB-MULTIPIER program. There are four lines of input
for the general control section. These lines are always required.
CONTROL line 1
NUMLC U=U1 D=DESIGN S=S1 V=VER N=NPLNOD Z=RSOIL E=DYTPE H=DYRELAX P=PRELOAD
R=NPUSH S=STIFF K=NSTIF F=PROTL=RIGSTIFF Y=TRANSSTIFF B=RIGMOD G=TRANSMOD A=TORSION
W=SPANTEMPERATURE line 2
Where
NUMLCis the number of load cases (INTEGER)
U1= 0 is mixed for English units (ld-f, kips, inches, ft)
= 1 is for SI units (kilonewtons and meters)
= 2 is for metric (kilonewtons and millimeters)
DESIGN is the flag for the preliminary design option(currently limited to off, or '0')
S1= 0 for no stiffness creation
= 1 for stiffness creation
VER=Version for English Units
=0 is for English units are Kip & Inches.
=1 is for metric units of newtons and meters.
=2 is for metric units of newtons and millimeters.
NPLNOD= number of nodes per pile
RSOIL is the flag for the rotated soil option (not currently available in GUI)
=0 do not use rotated soil option
=1 use rotated soil option
DYTYPE is the flag for Dynamic Analysis Type (which is selected on the Dynamics Page)
=0 Time Step Integration
=1 Modal Response
DYRELAX is the flag for Dynamic Relaxation (which is selected on the Dynamics Page)
=0 do not apply dynamic relaxation
=1 apply dynamic relaxation
PRELOADis the preload option
=0 for no preload
=1 for preload
NPUSHis the pushover option
=0 for normal analysis
=1 for pushover analysis
STIFF is the stiffness creation option
=0 no stiffness creation
=1 stiffness creation - one column, or (foundation stiffness option and modal analysis)
=2 stiffness creation - with multiple columns (modal analysis)
NSTIFis the type of stiffness.
= 0 secant stiffness (default).
= 1 tangent stiffness. (required for preload and time step dynamics, optional for a static analysis)
PROTis the soil resistance due to pile rotation about the 2 and 3 axis
= 0 do not include soil resistance
= 1 include soil resistance (default)
RIGSTIFF is the flag for Rigid Link Stiffness (not currently available in GUI)
TRANSSTIFF is the flag for Transfer Beam Stiffness (not currently available in GUI)
RIGMOD is the flag for Rigid Link Modulus (not currently available in GUI)
TRANSMOD is the flag for Transfer Beam Modulus (not currently available in GUI)
TORSION is the flag for Torsion PG(not currently available in GUI)
SPANTEMPERATURE is the flag to indicate whether or not span temperature loading is utilized, and if so, the
scope of activity.
= -1 is for use when temperature loading is utilized while in AASHTO analysis mode.
= 0 is for when span temperature loading is not utilized.
= 1 is for when span temperature loading is utilized in non-AASHTO analyses.
NUMLC is the number of different load cases for analysis within a single input file. Each load case is a separate
analysis of the same structure with a different set of loads. It is intended to reduce the number of input files that
need to be created in order to analyze a bridge pier structure. An input of 0 will execute the data check mode.
This halts the execution of FB-MultiPier after all the structure data has been generated and writes to the plot
database files for viewing with post processor. The latter is useful for data checking.
U1 identifies the standard FDOT English or metric pile sections.
S1 identifies that the analysis is to create an equivalent foundation stiffness. There can be no structure, only a
cap. All six loads must be given (three forces and three moments). These will be applied at the center of the pile
cap.
S=IFLEX T=ITIP, TSTIF P=NLOPT F=PHI PHITENSION PHI FLEXURE line 3
Where
IFLEXcontrols how the soil is to be modeled (INTEGER)
=0 user supplied P-multipliers must be given
=1 all user supplied P-multipliers are set to 1 internally in FB-MULTIPIER
=2 pile restraint only occurs through tip springs (i.e. no soil); soil information may be supplied, but is ignored.
ITIPis for the linear tip spring option (IFLEX=2) (INTEGER)
=0 for no linear tip springs on piles
=1 for axial tip springs on piles of stiffness TSTIF
=2 all d.o.f. at tip have springs with stiffness TSTIF
TSTIFis the stiffness of linear tips springs (REAL)
PHI is the user defined Phi factor over-ride used for creating the interaction diagrams. (REAL)
PHI TENSION is the user-defined tension Phi factor (not currently available in GUI)
PHI FLEXURE is the user-defined flexure Phi factor (not currently available in GUI)
NLOPT chooses linear or nonlinear piles
=1 for linear piles
=2 for nonlinear piles (cracked concrete, steel yielding and P-).
=3 for linear piles where interaction diagrams are generated
The no soil model (IFLEX=2) can be useful in testing the model and comparing its results to other solutions. In
this case, the user must make sure the structure is stable through the proper use of tip springs (ITIP) and pile cap
fixity (KFIX). The tip spring model allows the user to add either linear springs to the axial (ITIP=1) or to all (ITIP=2)
degrees of freedom at the bottom of each pile. In the case of IFLEX=0 or 1, ITIP or TSTIF are still active in addition
to any soil tip properties specified through the use of soil tip modeling.
I=MAXITER T=TOLER M=MEM X=TRANS V=VERL=LETTER line 4
Where
MAXITERis the maximum # of iterations for the nonlinear soil analysis (INTEGER)
TOLERis the tolerance on the maximum out-of-balance force for any node in the system in the nonlinear
analysis (REAL)
MEMis the amount of memory used during analysis, always in Megabytes (MB)
TRANSis the option to use the transformed section properties
= 0 do not use the transformed section
= 1 do use the transformed section
VERVersion number used to generate input file
LETTERis the version letter. This value, combined with the version number, fully describes the version of the
program
The out of balance forces are obtained in the following manner. The stiffness matrix is multiplied times the
current set of displacements to obtain a force vector. This force vector is then compared with the applied forces
on the structure. If the structure is in static equilibrium then the two force vectors would be identical. The
difference between the two sets of forces are the out of balance forces.
The following default values are be used for the maximum # of iterations for nonlinear analysis (MAXITER) and
the tolerance on the out-of balance forces (TOLER) for convergence:
MAXITER= 50
TOLER= 1.0
FB-MULTIPIER offers the option to use linear or nonlinear piles and piers. Linear piles will converge more
quickly and should be used for preliminary design and when nonlinear sections are not significant. NLOPT (on
previous line) chooses which type of pile behavior will be used.
The following information is used by with multiple pier generation. The information under the SUBSTR header
describes the pier to pier geometry.
SUBSTR
NPIERS
PIERNUM X= XDIST Y= YDIST R= ZROT A= LBROT B= RBROT
Where
NPIERS is the number of piers (or pile bents)
XDIST is the x-distance from the global origin (located at the first pier)
YDISTis the y-distance from the global origin
ZROTis the pier rotation angle about the global z-axis (used for curved alignments)
LBROTis the left bearing row rotation angle about the global z-axis (used for curved alignments)
RBROTis the right bearing row rotation angle about the global z-axis (used for curved alignments)
This section must end with a blank line.
The pier rotation option allows each pier to be rotated about the global z-axis (clockwise positive) to
accommodate skew and curved bridge alignments. The rotation is performed as follows based on the bottom left
corner of the pile cap
The following information is used by with multiple pier generation. The information under the SUPPROP header
describes the superstructure properties. These properties are used in establishing the pier to pier connectivity.
SUPPROP
A=AREA, AREASUP I=I3, I2 J= TOR E=EMOD G=GMOD S=WEIGHT F=BEGIN X=AREA, I3, I2, TOR, EMOD,
GMOD B=END Y=AREA, I3, I2, TOR, EMOD, GMOD C= CONTIN Z= AREA, I3 , I2, TOR, EMOD, GMOD T=
BEGINVERT K= AREA, I3, I2, TOR, EMOD, GMOD P= ENDVERT Q= AREA, I3, I2, TOR, EMOD, GMOD H=H1,
H2 L=LH V=VARSP N=NODES
(One line for each span)
If variable span properties are present, each property line that follows applies to an element along the span.
There are a total of 10 elements per span. The above line contains the properties for Element #1.
A=AREA I=I3, I2 J= TOR E=EMOD G=GMOD(One line for each span element for variable span properties)
Where
AREAis the transverse cross-sectional area of the corresponding superstructure component
AREASUPis the span profile area used for wind load generation
I3is the moment of inertia of the corresponding superstructure component (strong axis bending)
I2is the moment of inertia of the corresponding superstructure component (weak axis bending)
TORis the torsional moment of inertia of the corresponding superstructure component
EMODis the elastic modulus of the corresponding superstructure component
GMODis the shear modulus of the corresponding superstructure component
WEIGHTis the unit weight of the superstructure
BEGINis the flag for the begin (left) span end condition (0 = Diaphragm; 1 = No Diaphragm; 2 = Custom)
ENDis the flag for the end (right) span end condition (0 = Diaphragm; 1 = No Diaphragm; 2 = Custom)
XProperty list for begin (left) span end condition
YProperty list for end (right) span end condition
CONTIN is the flag to indicate custom input of continuity link properties (0=default; 1=custom)
ZProperty list for continuity link element
BEGINVERTis the flag to indicate custom input of begin vertical link properties (0=default; 1=custom)
ENDVERTis the flag to indicate custom input of end vertical link properties (0=default; 1=custom)
KProperty list for begin vertical link element
PProperty list for end vertical link element
H1is the vertical distance from the center of the pier cap to the center of gravity of the superstructure (beginning
of span)
H2is the vertical distance from the center of the pier cap to the center of gravity of the superstructure (end of
span)
LHis the live load height from the center of the pier cap to the center of the live load
VARSPindicates whether variable span properties are present (0;no, 1;yes).
NODES is the total number of nodes on the bridge span (not currently available in GUI, and is set to 11, creating
10 span elements per span)
:
This section must end with a blank line.
Each span of the superstructure is modeling with a single beam (divided into sub-elements) that spans from the
center of the back span pier to the center of the forward span pier. The superstructure beam is connected to rigid
beams at the back span and forward span, which distribute the load to the bearing locations.
The following information is used by with multiple pier generation. The information under the PADPROP header
describes the load-displacement behavior for the bearing locations. This information is only provided for userdefined substructure to superstructure connectivity.
PADPROP
NPROP
X1,X2,X3,,X20
F1,F2,F3,,F20
T=TYPE1,TYPE2.....TYPE100
Where
NPROPis the number of custom load-displacement curve definitions
F1F20are the load values in the load-displacement relationship (20 points max.)
X1X20are the displacement values in the load-displacement relationship (20 points max.)
TYPE1....TYPE100 are the bearing curve types. A value of 1 is written for translational, a value of 2 for rotational.
(Repeat curve pair values for each user-defined curve definition)
:
This section must end with a blank line.
A maximum of 20 points can be used to describe the load-displacement relationship for the bearing. Load values
should be entered for both positive and negative displacements. Zero force values can be entered intentionally.
For example, to model no vertical reaction due to girder uplift, enter a load-displacement curve for the vertical
displacement with positive load and displacement values (for loads acting downward) and zero force values for
negative displacement (loads acting upward). Displacement values must be entered in order from the largest
negative displacement to the largest positive displacement.
Examples:
17.2.9 Pushover
For static pushover analysis,
PUSH
N=NPSTEP I=PINCR
Where
NPSTEPis the number of load steps
PINCR is load increment (e.g. 1.0 = increase initial load by 100% for each load step)
2 load cases are required
Use load case #1 for permanently applied loads
Use load case #2 for initial load to be incremented
Under the CONTROL header,
R=NPUSH
Where
NPUSH = 0 for normal analysis
= 1 for pushover analysis
line1
For Load Case 1, store the load case number (LCNUM) and indicate if the load case will utilize the PreLoad
feature (LCFLAG), where a value of 1 indicates that the PreLoad feature is used for the corresponding load
case number, and a value of 0 indicates that the PreLoad is not being used for the corresponding load case
number.
L= LCNUM P = LCFLAG
line2
Repeat input for each load case (excluding the PRELOAD header).
= 1 for STRENGTH-I
= 1 for STRENGTH-II
= 1 for STRENGTH-III
= 1 for STRENGTH-IV
= 1 for STRENGTH-V
= 1 for EXTREME-I
= 1 for EXTREME-II
= 1 for SERVICE-I
= 1 for SERVICE-II
= 1 for SERVICE-III
= 1 for SERVICE-IV
= 1 for FATIGUE-I
= 1 for FATIGUE-II
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(LFD)
L1
L2
L3
L4
L5
L6
L7
L8
L9
L10
L11
= 1 for GROUP-I
= 1 for GROUP-IA
= 1 for GROUP-III
= 1 for GROUP-II
= 1 for GROUP-IV
= 1 for GROUP-V
= 1 for GROUP-VI
= 1 for GROUP-VII
= 1 for GROUP-VIII
= 1 for GROUP-IX
= 1 for GROUP-VESSEL
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
(0 otherwise)
Where
NDYNSis the type of dynamics solution.
= 0 Step by step integration (default).
= 1 Spectrum analysis (for the structure only).
NDAMPis the damping option.
= 0 no damping (default).
= 1 damping.
ALPHA1,BETA1coefficients for Rayleigh damping for the structure (C=ALPHA*M+BETA*K).
ALPHA2,BETA2coefficients for Rayleigh damping for the piles.
ALPHA3,BETA3coefficients for Rayleigh damping for the soil.
NPMAXis the maximum number of time steps for the analysis.
DMPis the mass density for the piles. This a default global value used if a weight density if not given for a pile
cross-section.
DMSis the mass density for the pier (excluding pile cap). This a default global value used if a weight density if
not given for a pier cross-section.
NPRTis the output option for time step analysis.
= 0 maximum displacements and maximum forces caused by the maximum displacements (default).
= 1 all displacements and maximum forces.
= 2 maximum displacements and all forces.
= 3 all displacements and all forces.
Note that NPRT = 2 or 3 options only allow the program to compute the element forces for the options above
(because this may take some time for a large structure), to print them out, you still have to use these in addition to
the print out option. For example, if you want the pile forces for every time step use O=4 and set P=1 under the
PRINT label. If you want only the structure forces use T=1. The maximum forces are the forces caused by the
maximum displacements, note that these can be smaller than the maximum forces for the structure. For options 2
and 3 a summary of the maximum forces and the time step when it occurred will be printed out at the end.
SMASS is the concentrated mass adopted for the soil. This mass is applied to all the translational DOF X, Y and
Z to represent the attached soil mass.
NSHM is the option for the mass matrix for the cap.
= 0 consistent mass matrix (default).
= 1 lumped mass matrix.
NBMMis the option for the mass matrix for the structure and piles.
= 0 consistent mass matrix (default).
= 1 lumped mass matrix.
NDYSOL is the option for the type of numerical solution. This option is only valid for step by step solution
(NDYNS=0)
= 0 Newmark's method - average acceleration (default).
= 1 Newmark's method - linear acceleration. This method is conditionally stable and is not guaranteed to
converge unless dt/T < 0.551 for all vibration modes.
= 2 Wilson-Theta method.
NMSE is the option for multiple support excitation.
= 0 standard analysis.
= 1 multiple support excitation.
EPP is the strain rate for concrete (default = 1e-5)
NCMODis the option for the concrete model in nonlinear analysis.
= 0 ; nonlinear elastic (default)
= 1 ; hysteresis
NSTEELis the steel model
= 0 ; bilinear (default)
= 1 ; strain - hardening and Bauschinger effect
D1...D7this is an option if the user wants to save a specific NODE displacement to a file for possible later plotting
or checking.
D1 - is the number of NODES that will me saved (maximum = 6).
D2...D7 - is the number of the NODE to be saved. For example:
P=6,1,2,3,4,5,6
D1 = 6 - six NODES will be saved to the file.
1...6 - NODES 1 to 6 will be saved to the file. They do not have to be in any specific order.
Note that the displacements will be saved in a text file with the name: 'inputname'.DSn, the velocities in
'inputname'.VSn, and the accelerations in 'inputname'.ASn, where n is the file number for the chosen node. In the
example above if the input file is called test.in, the displacements for node 1 are saved in test.DS1, for node 2 in
test.DS2 and so on.
NFREQ is the option for computing the period.
= 0 (default) - initial period for the structure will not be computed.
= 1- computes the initial periods for the structure.
CFRQ is the frequency of loading (used for cyclic degradation, in Hz)
NSCS is the option for subtracting the stress in the concrete at the level of the steel bar for dynamic analysis.
= 0 (default) - does not subtract.
= 1 - subtracts.
JPIL is the option to track the forces in a specific pile element.
JPIL(1) = pile number
JPIL(2) = element number, must be between 1 and 16. 1 is top element, 16 is bottom element. If not specified is
set to 1.
JPIL(3) = element node, must be 1 or 2. 1 is bottom element node, 2 is top element node. If not specified is set to
1.
This information is saved to the file INPUT.DFO, where INPUT is the input file name.
This section MUST end with a blank line.
Options shown with a strikethrough font will be implemented in the future.
Where
Ti,Fiare the time and force (or acceleration) values for the point being specified.
This section MUST end with a blank line.
Z= LFN
Where
LFNis the name of the load function. This can be up to 20 characters in length.
C=PB
Where
PB is the Peer Berkeley decryption string, to describe the Load Function.
SPECTRUM DEFINITION LINES (Repeat for as many lines as necessary)
The next lines specify the spectrum function values. The lines MUST contain FOUR pairs of numbers each. The
number of lines is dictated by the number of points used to specify the function (SP). The points do NOT need to
be at even spacing.
F1,A1 F2,A2 F3,A3 F4,A4
Where
Fi,Aiare the frequency (rad/sec) and acceleration (length/sec^2) values for the point being specified.
The acceleration values should represent the maximum absolute acceleration (spectral
acceleration), not the relative acceleration (between the structure and moving support).
This section MUST end with a blank line.
For steel girders (SUPXSECTYPE= 1) that are prismatic and identical in cross-section (GIRGEOM= 0), the
following parameters are read once. For steel girders (SUPXSECTYPE= 1) that are prismatic and mutually
varying in cross-section (GIRGEOM= 1), the following parameters are read once for each girder. For steel girders
(SUPXSECTYPE= 1) that are non-prismatic but identical (GIRGEOM= 0), the following parameters are read once
for each span element. For steel girders (SUPXSECTYPE= 1) that are non-prismatic and mutually varying in
cross-section (GIRGEOM= 1), the following parameters are read for each girder of each span element.
Specify bottom flange depth (BFLDEPTH) and bottom flange width (BFLDEPTH).
D= BFLDEPTH W= BFLWIDTH line 4
Specify web depth (WEBDEPTH) and web width (WEBWIDTH).
E= WEBDEPTH X= WEBWIDTH line 5
Specify top flange depth (TFLDEPTH) and top flange width (TFLWIDTH).
F= TFLDEPTH Y= TFLWIDTH line 6
For concrete girders (SUPXSECTYPE= 2) that are prismatic and identical in cross-section (GIRGEOM= 0), the
following parameters are read once. For concrete girders (SUPXSECTYPE= 2) that are prismatic and mutually
varying in cross-section (GIRGEOM= 1), the following parameters are read once for each girder. For concrete
girders (SUPXSECTYPE = 2) that are non-prismatic but identical (GIRGEOM = 0), the following parameters are
read once for each span element. For concrete girders (SUPXSECTYPE = 2) that are non-prismatic and mutually
varying in cross-section (GIRGEOM = 1), the following parameters are read for each girder of each span
element.
Specify bottom flange depth (BFLDEPTH) and bottom flange width (BFLWIDTH).
D= BFLDEPTH W= BFLWIDTH line 4
Specify bottom taper depth (ie, depth of taper between bottom flange and web) (BTDEPTH).
E= BTDEPTH line 5
Specify web depth (WEBDEPTH) and web width (WEBWIDTH).
F= WEBDEPTH X= WEBWIDTH line 6
Specify top taper depth (ie, depth of taper between web and top flange) (TTDEPTH).
G= TTDEPTH line 7
Specify top flange depth (TFLDEPTH) and top flange width (TFLWIDTH).
H= TFLDEPTH Y= TFLWIDTH line 8
For box girders (SUPXSECTYPE = 3) that are prismatic and identical in cross-section (GIRGEOM = 0), the
following parameters are read once. For box girders (SUPXSECTYPE = 3) that are prismatic and mutually
varying in cross-section (GIRGEOM = 1), the following parameters are read once for each girder. For box girders
(SUPXSECTYPE = 3) that are non-prismatic but identical (GIRGEOM = 0), the following parameters are read for
on girder of each span element. For box girders (SUPXSECTYPE = 3) that are non-prismatic and mutually
varying in cross-section (GIRGEOM = 1), the following parameters are read for each girder of each span
element.
Specify girders bottom flange depth (BSLDEPTH) and the girders bottom slab width (BSLWIDTH).
D= BSLDEPTH W= BSLWIDTH line 4
Specify number of web-walls (NWEBW) and the width at top of web walls (ie, the horizontal distance from the
top, outside of the outermost web walls) (WEBWDIST).
N= NWEBW O= WEBWDIST line 5
Specify the web wall depth (WEBWDEPTH) and web wall width (ie, the average web-wall width)
(AWEBWWIDTH).
E= WEBWDEPTH X= AWEBWWIDTH line 6
Specify the top flange depth (TFLDEPTH) and top flange width (TFLWIDTH).
line 1
For Span 1, define the reference temperature (REFT), and if using the AASHTO Combinations Design Options
feature, the uniform span temperature (UNIFT).
R= REFT U=UNIFT line 2
For Span 1, skip the next line (the perimeter SPANNO is only read by the GUI).
S= SPANNO line 3
For Span 1, define the coefficient of thermal expansion (GALPHA) for the girders.
A= GALPHA line 4
For Span 1, define the coefficient of thermal expansion (SALPHA) for the slab.
B= SALPHA line 5
For Span 1, specify the slab top (SLABTOPT) temperature, the temperature 4 in. below the slab top (SLABFINT),
and the slab bottom temperature (SLABBOTT). Note that these temperatures are read only once for each span.
S= SLABTOPT T=SLABFINT U= SLABBOTT line 6
For steel girders that are prismatic and identical in cross-section (as determined in the SPANDIM header), the
following parameters are read once. For steel girders that are prismatic and mutually varying in cross-section,
the following parameters are read once for each girder. For steel girders that are non-prismatic but identical in
cross-section, the following parameters are read once for each element. For steel girders that are non-prismatic
and mutually varying in cross-section, the following parameters are read for each girder of each span element.
Specify the girder temperatures at the bottom of the girder (T1), top of the bottom flange (T2), first arbitrary web
depth (T3), second arbitrary web depth (T4), third arbitrary web depth (T5), top of the web (T6), and top of the top
flange (T7).
T= T1 T2 T3 T4 T5 T6 T7 line 7
Specify the arbitrary vertical offsets (D3, D4, and D5) that correspond to the arbitrary web depth temperatures
(T3, T4, and T5). Note that the offsets are cumulative, where D3 is taken relative to the bottom of the web, for the
current element being defined.
D= D3 D4 D5 line 8
For concrete girders that are prismatic and identical in cross-section (as determined in the SPANDIM header),
the following parameters are read once. For concrete girders that are prismatic and mutually varying in crosssection, the following parameters are read once for each girder. For concrete girders that are non-prismatic but
identical in cross-section, the following parameters are read once for each element. For steel girders that are
non-prismatic and mutually varying in cross-section, the following parameters are read for each girder of each
span element.
Specify the girder and slab temperatures at the bottom of the girder (T1), top of the bottom flange (T2), top of the
bottom flange to web taper (T3), first arbitrary web depth (T4), second arbitrary web depth (T5), third arbitrary
web depth (T6), bottom of the web to top flange taper (T7), bottom of the top flange (T8), and top of the top flange
(T9).
T= T1 T2 T3 T4 T5 T6 T7 T8 T9 line 7
Specify the arbitrary web depths (D4, D5, and D6) that correspond to the arbitrary web depth temperatures (T4,
T5, and T6). Note that the offsets are cumulative, where D4 is taken relative to the bottom of the web, for the
current element being defined.
D= D4 D5 D6 line 8
For box girders that are prismatic and identical in cross-section (as determined in the SPANDIM header), the
following parameters are read once. For box girders that are prismatic and mutually varying in cross-section, the
following parameters are read once for each girder. For box girders that are non-prismatic but identical in crosssection, the following parameters are read once for each element. For box girders that are non-prismatic and
mutually varying in cross-section, the following parameters are read for each girder of each span element.
Specify the box-girder and slab temperatures at the bottom slab of the box-girder (T1), top of the box-girder
bottom slab (T2), first arbitrary web-wall depth (T3), second arbitrary web-wall depth (T4), third arbitrary web-wall
depth (T5), top of the web-wall (T6), and top of the top flange (T7).
T= T1 T2 T3 T4 T5 T6 T7 line 7
Specify the arbitrary web depths (D3, D4, and D5) that correspond to the arbitrary web depth temperatures (T3,
T4, and T5). Note that the offsets are cumulative, where D3 is taken relative to the bottom of the web, for the
current element being defined.
D= D3 D4 D5 line 8
Repeat input for each span (excluding the SPANTHERM header and the reference temperature, REFT).
Note: 1. In the English system, units are feet for span offsets.
2. In the SI system, units are meters for span offsets.
line 2
line 3
line 4
line 5
line 6
line 7
line 8
line 9
line 10
line 11
line 12
line 13
where
X, Y, Z, P, Q, R correspond to the rows of the 6x6 stiffness matrices, and where lines 2-7 pertain to the left extent
of the OPTS model, while lines 8-13 pertain to the right extent of the OPTS model.
S11, S12, S13 are the individual stiffness terms. For each of the six data rows allotted for the left and right
model extents, the rows correspond to X-translation, Y-translation, Z-translation, X-rotation, Y-rotation, and Zrotation degrees-of-freedom, respectively. For example, entry S44 of line 5 corresponds to the diagonal value of
X-rotational stiffness for the left extent of the OPTS model.
Note: For Spring Stiffness, the English units are kip/in and kip-ft/rad for translational and rotational terms,
respectively. SI units are kN/m and kN-m/rad, respectively.
:
where
NODE is the span node number to which the Mass will be applied
MTX is the additional mass in the x-axis translational DOF
MTY is the additional mass in the y-axis translational DOF
MTZ is the additional mass in the z-axis translational DOF
MRX is the additional mass in the x-axis rotational DOF
MRY is the additional mass in the y-axis rotational DOF
MRZ is the additional mass in the z-axis rotational DOF
SPAN is the span to which the mass(es) will be applied
Note: For Span mass, the English units are kip-sec^2/in and kip-sec^2/rad for translational and rotational
DOF, respectively.
N=NABPBU X=ABMAXU
NABPBU is the number of points in current unloading curve [limited to 25 points]
ABMAXU is the crush level at which the loading and unloading curves intersect
ABU(1) PBU(1)
ABU(2) PBU(2)
ABU is the abscissa for the current user-defined vessel bow unloading curve
PBU is the ordinate for the user-defined vessel bow unloading curve
STYPE is the segment type code, where HPILE = 0, PRECAST = 1, DRILLED SECTION = 2, and PIPE= 3.
DSEG is the segment data type, where UNKNOWN = 0, LINEAR = 1, and NONLINEAR= 2.
SUNIT is the units code, where ENGLISH = 0, and METRIC = 1.
SHAPE is the shape description. These include H-Pile, PreCast, Drilled Shaft, Pipe Pile, and Generic.
LNR is the linear / non linear decription. These include nonlinear and linear.
Finally, for both linear and nonlinear piles, six additional lines are needed to define the pile geometry.
or
For user specified Stress Strain Curves (MATOPT=2) plus up to 5 sets of stress strain points for user defined
curve
K=KTYPE L=XPL M=MATOPT S=KSTEEL(1), KSTEEL(2), KSTEEL(3), KSTEEL(4)
Where
XPLis the pile length for this segment for plumb and battered piles (REAL)
MATOPT is the material input option (INTEGER)
MATOPT=1 means input FPC, FY or FSU,ES,EC and KSTEEL on this line and default stress strain curves
will be generated.
MATOPT=2 means describe stress strain curves for steel and concrete in INPUT #6B and #6C
No FPC, (FY or FSU), ES and EC values to be entered for MATOPT=2.
FPC is the compression stress, f'c, for concrete (REAL)
FPC = 0 for tubular steel sections
ECis the modulus of elasticity of concrete (REAL)
FY(1) is the yield stress, Fy , for mild steel (REAL)
FSU(2)is the ultimate stress for prestressed strands (REAL)
FY(3)is the yield stress for H-pile section (REAL)
FY(4)is the yield stress for tubular steel section (REAL)
(only for circular pile)
ES(1)is the modulus of elasticity of mild steel (REAL)
ES(2)is the modulus of elasticity of prestressing strand (REAL)
ES(3)is the modulus of elasticity of H-pile section (REAL)
ES(4)is the modulus of elasticity of tubular steel section (REAL)
KSTEEL(I)is the steel type option (INTEGER)
KSTEEL(I) = 1 includes steel type
KSTEEL(I) = 0 does not include steel type
KSTEEL(1) for mild steel reinforcement
KSTEEL(2) for prestressing steel strands
KSTEEL(3) for H-pile section
KSTEEL(4) for tubular steel section (only for circular piles)
KTYPEspecifies the cross sectional shape of the pile
KTYPE=1 for a round pile.
KTYPE=2 for a square pile.
KTYPE=3 for an H-pile
Tubular and H-pile steel sections can be input by negating concrete as described above or in the non-linear user
defined stress strain curves and inputting the section properties described in the sections for the input of circular
piles and H-piles.
Stress-Strain Curve for Concrete, used with NLOPT=2 or 3 and MATOPT=2
NC=NPCC, SIGC(1), SIGC(2),,, line 1
EPSC(1), EPSC(2),,,line 2
Where
NPCCis the number of points on concrete curve (INTEGER)
line 1
Stress-Strain Curve for Tubular Steel, used with NLOPT=2 and MATOPT=2 and KSTEEL(4) = 1
S4=NPSC, SIGS(1), SIGS(2),,,line 1
EPSS(1), EPSS(2),,, y= y
line 2
where
NPSCis the number of points on the tubular steel curve (INTEGER)
SIGS(1) is the first stress point on the tubular steel curve (REAL)
SIGS(2) is the second stress point on the tubular steel curve (REAL)
EPSS(1) is the first strain point on the tubular steel curve (REAL)
EPSS(2) is the second strain point on the tubular steel curve (REAL)
yis the steel yield strain
For Nonlinear Analysis of Square/Rectangular Piles, used with NLOPT=2 or 3 and KTYPE=2
W=WIDTH D=DEPTH V=DV B=BV P=PREST N=ISTNOPT S=SW
Where
WIDTH is the width of square pile parallel to the local Z axis (REAL)
DEPTHis the depth of rectangular pile parallel to the local Y axis (REAL)
DVis the diameter of the void (DV=0 for no void) (REAL) (It can also be the width of a rectangular void IF BV>0)
BVis the depth of a rectangular void. If BV>0, then the void is rectangular and the DV value is used for the width
of the void.
PRESTis the prestressing stresses after release & all losses for standard sections only(AASHTO 9.15.1, 9.16.2)
(REAL)
PREST=0 for non-prestressed (i.e. reinforced concrete)
ISTNOPT is the standard section option (INTEGER)
ISTNOPT=1 means use FDOT standard reinforcement for input width as shown below (INTEGER)
Note: WIDTH MUST BE one of the values from a) through f) from Figure: 17.3.b or Figure: 17.3.c
STNOPT=2 means describe the reinforcement in the section for the nonlinear analysis of nonstandard
rectangular piles. (Use next lines)
SWis the unit weight of the pile, used for self-weight calculations. If SW>0, self-weight is included in the
analysis.
Figure: 17.3.c Standard FDOT Prestressed Concrete Pile Sections (meters, kN)
Figure: 17.3.d Standard FDOT Prestressed Concrete Pile Sections (millimeters, kN)
For nonlinear Analysis of Nonstandard Square/Rectangular Piles, used with NLOPT=2, KTYPE=2, and ISTNOPT=
2
NG=NGRPS HPI= IHPILE BM=BMETH X=MINSPACE Z=TYPE
AS, Y, Z, PREST N=N1 D=D1 repeat NGRPS times
Where
NGRPSis the # of groups of bars/strands(INTEGER)
IHPILEis the H-pile option.
IHPILE = 1 for H-pile embedded in the concrete,
else IHPILE = 0
BMETH0=Custom; 1=Percentage
MINSPACEMinumum spacing between two bars
TYPEBar Type Number
ASis the bar or strand areas (REAL)
Yis the local Y coordinate for bar or strand (REAL)
Zis the local Z coordinate for bar or strand (REAL)
PRESTis the prestressing stress in the strands after all losses (REAL)
N=N1 D=D1is code to generate multiple bars in (INTEGER)
N=N1 D=2means generate N1 bars/strands in the local Y direction as follows:
First bar is at coordinates Y,Z
if N1 = 2, second bar is at coordinate -Y,Z
if N1 > 2, then second bar is at coordinate -Y,Z
and remaining N1-2 bars/strand are equally spaced between first two bars/strands
EXAMPLE INPUT
For steel H-piles used with KTYPE=3 or HP=1 in either circular or square sections
Two lines are required:
OR=ORIENTline 1
[D=DEPTH U=WEIGHT]line 2, for standard H-pile sections
or
[D=DEPTH TW=WEB B=WIDTH TF=FLANGE] line 2, for user defined sections
Where
ORIENTis the orientation of the H-pile.
ORIENT=2 for web parallel to the Local Y axis, or 3 for web parallel to thelocal Z-axis. (INTEGER)
DEPTH is the depth of the H-pile in inches (REAL). (Use the nominal depth for standard sections)
WEIGHT is the standard unit weight of the H-pile in lb/ft3 (REAL)
WEB is the web thickness in inches (REAL)
WIDTH is the flange width in inches (REAL)
FLANGEis the flange thickness in inches (REAL)
Note: For metric examples H-pile dimensions will be soft converted to metric units.
After the cross section data is input, SIX additional lines defining the pile system are required
For Nonlinear Analysis of Oblong Piers, used with NLOPT=2 and KTYPE=4 NOTE: This type is ONLY
available for pier elements NOT for piles.
R= OR RW=RWIDTH D=DIAM T=VT V=DV B=WV S=WDEN
Where
ORis 2 or 3 and defines orientation, see Figure: 17.3.k (INTEGER)
RWIDTH is the width of rectangular portion (REAL)
DIAMis the diameter of semi-circular ends (REAL)
VTis void type, see Figure: 17.3.k VT may be 1 or 2 for OR = 2 or 3 (INTEGER)
DVis the diameter of the void for VT=1 (DV=0 for no void)
DVis the depth of the void parallel to DIAM for VT=2 (DV=0 for no void) (REAL)
WVis the width of the void parallel to RWIDTH for VT=2
(WV= 0 for no void) (REAL)
WDENis the self weight of the concrete
Reinforcement specification (Rectangular middle is similar to steel generation for rectangular sections)
NG=NGRPS >M=BMETH X=MINSPACE Z=TYPE
AS, Y, Z, PREST N=N1 D=D1 repeat NGRPS times
Where
NGRPSis the # of groups of bars/strands (INTEGER)
BMETH 0=Custom, 1=Percentage
MINSPACE Minumum spacing between two bars
TYPE Bar Type Number
ASis the bar or strand areas (REAL)
Yis the local Y coordinate for bar or strand (REAL)
Zis the local Z coordinate for bar or strand (REAL)
PRESTis the prestressing stress in the strands after all losses (REAL)
N=N1 D=D1is code to generate multiple bars in (INTEGER)
N=N1 D=2means generate N1 bars/strands in the local Y direction as follows:
First bar is at coordinates Y,Z
if N1 = 2, second bar is at coordinate -Y,Z
if N1 > 2, then second bar is at coordinate -Y,Z
and remaining N1-2 bars/strand are equally spaced between first two bars/strands
Reinforcement specification: (Semi-circular ends are similar to steel generation for circular sections.)
NL=NLAY
PREST, NBS, D=DSL, A=ASIrepeat NLAY times
Where
NLAYis the number of circumferential steel layers (INTEGER)
NBSis the number of bars in the layer (INTEGER)
PRESTis the effective prestressing stress in the strands for the layer
Note: If mild steel is present along with prestressing strands, the prestressing stress on the concrete is reduced
due to the area of mild steel, and the strain in the concrete due to the prestressing is assumed to be shared with
the mild steel.
Required for all types: Input for free length, number of sub-elements, axial efficiency and pile head fixity
F=FLNG H=KFIX S=NSUB A=AXEFF G=GAP C=KBCAPpile configuration line 1
or
E=ECAP H=KFIX S=NSUB A=AXEFF G=GAP C=KBCAP
Where
FLNGis the length of pile between the pile cap and the ground surface, the free length. It can be zero. If < 0, the
cap is analyzed as a buried cap.(not currently written to GUI)
ECAPis the elevation of the pile cap. This is assumed at the top of the pile heads, which is the same as the
centroid of the pile cap. Since the pile cap is modeled using a shell element, the pier column base, the pile
heads and the neutral axis of the pile cap all meet in the same location. This modeling does NOT account for the
thickness of he pile cap in the geometry of the system (it is included in the behavior).
KFIXis for the pile head fixity into the cap (INTEGER)
Where
DX1is the spacing between the first and second row of piles in the X direction. (REAL)
DX2is the spacing between the second and third row of piles in the X direction. (REAL)
Pile Spacing in the Y-direction
The pile system may have even or uneven spacing in the Y direction. If only one value is given (DY1), then the
spacing is uniform. Otherwise, values MUST be given for each spacing value between every row of piles. There
must be NPY-1 values given for uneven spacing.
DY1, DY2,...pile configuration line 4
Where
DY1is the spacing between the first and second row of piles in the Y direction. (REAL)
DY2is the spacing between the second and third row of piles in the Y direction. (REAL)
Input for P-multipliers in the X-direction
P-multipliers used for the X direction given in order from trail to lead row of piles ( Figure: 17.3.m ). Multipliers
have to be specified for existing rows only. The program assigns the values in the correct order depending upon
the resultant loads in the X direction.
PYMX1, PYMX2, ...pile configuration line 5
Where
PYMX1is the multiplier for the trail row (REAL)
PYMX2is the multiplier for the second row (REAL)
P-multipliers used for the Y direction given in order from trail to lead row of piles. Multipliers have to be specified
for existing rows only. The program assigns the values in the correct order depending upon the resultant loads in
the Y direction.
PYMY1, PYMY2,...pile configuration line 6
Where
PYMY1is the multiplier for the trail row (REAL)
PYMY2is the multiplier for the second row (REAL)
PILESET
The next line can be repeated for as many pile types as need to be specified.
PILEx PSETx (repeat for each pile type greater than 1)
Where
PILEis the pile number to which the cross section type is applied.
PSETis the pile cross section set number to apply to this pile.
Example:
PILESET
12(pile # 1 of pile type #2)
22(pile # 2 of pile type #2)
32(pile # 3 of pile type #2)
63(pile # 6 of pile type #3)
This data is used to specify any removed piles. If none are removed, skip this section.
MISSING
NMPIL
Where
NMPILis the number of missing piles from the pile group (INTEGER). This value may be zero.
Specify missing piles by x-row, y-row pile coordinate system. The coordinate system of the pile rows is shown in
Figure 6. One line is used for each missing pile. Repeat the following lines NMPIL times.
IXORD, IYORD
Where
IXORDis the x row location of missing pile (INTEGER)
IYORDis the y row location of missing pile (INTEGER)
D=FM, G=FS, E=ETOP, EBOT, B=PBOT, S=STYPE, A=TANDB, T=PTOP N=N50 L=SLR V=SWVS F=SFDF
M=NSMOD; J=RSDAMP line one
, RK, , Cu, 50 or qu (Limestone), 100 or Cavg, or K, G, ,f S=Stype, A=TANDB or Fsmax line two
Where
is the angle of internal friction (REAL)
RKis the soil modulus k (REAL)
is the total unit weight of the soil (REAL)
Cuis the undrained shear strength (REAL)
JS is the joint spacing
50 is the major principal strain @ 50% maximum deviator stress in a UU triaxial compression test (REAL)
or
qu(Limestone)is the unconfined compressive strength of limestone
100 is the major principal strain @ failure in a UU triaxial compression test (SOIL=3) (REAL)
or
Cavg is the average undrained shear strength for the soil layer (REAL) (SOIL=5 & 6)
or
K is the Dimensionless Coefficient of Lateral Earth Pressure (REAL) (SOIL=10)
Gis the shear Modulus of the soil (REAL)
is Poisson's ratio of the soil (REAL)
f is the vertical failure shear stress on pile-soil interface (REAL)
or
Fsmax is the Ultimate Side Friction (REAL) (SOIL=10)
THICKNESSis the thickness of the soil layer (REAL)
LSMis the Lateral Soil Model. It selects one of seven different lateral P-Y curves (INTEGER)
1 =Sand (O'Neill, 1984) requires , RK,
2 =Sand (Reese,Cox,Koop, 1974) requires , RK,
3 =Clay (O'Neill) requires Cu, 50, 100,
4 =Clay - Soft clay below water table; (Matlock, 1970) requires , Cu, 50
5 =Clay - Stiff clay below water table; (Reese, 1975) requires RK, , Cu, 50, Cavg
6 =Clay - Stiff clay above water table; (Reese, 1975) requires , Cu, 50, Cavg
7 =user defined P-Y curve for lateral soil response. Requires four additional lines of input (2 for top and 2 for
bottom of layer).
8 =Limestone (McVay)
9 =Limestone (McVay) No 2-3 Rotation
10 =Sand (API)
11 =Clay (API)
ASMis the axial soil model. There are 5 allowable axial soil models.
1 =Driven Pile (McVay et al, 1989) requires G, , f
2 =Drilled Shaft on Sand (ONeill et al, 1996) requires
3 =Drilled Shaft on Clay (ONeill et al, 1996) requires Cu
4 =Drilled Shaft on Intermediate Geo Material IGM (O'Neil) requires Surface Type, qu,Core Recovery, Em,
Em/Ei
5 =user defined T-Z curve. Requires four additional lines of input (2 for top and 2 for bottom of layer)
ASM type 4
1 =Rough surface
2 =Smooth surface
qu(IGM)is the unconfined compressive strength for intermediate geomaterialASM type 4
or
fsmaxis the ultimate unit skin friction
Core Recoveryis the IGM core recovery in percentageASM type 4
Em is the IGM mass modulusASM type 4
Em/Ei is the ratio of IGM mass modulus to intact material modulusASM type 4
qtis the split tensile strength (used only rough surface and Florida Limestone) ASM type 4
JZ is the joint spacing (REAL)
FM is the fracture constant M
FS is the fracture constant S
ETOPis the elevation at the top of this soil layer
EBOTis the elevation at the bottom of this soil layer
PTOPis the elevation of the piezometric head at the top of the layer
PBOTis the elevation of the piezometric head at the bottom of the layer
STYPEis the soil layer type (used for graphical interface only)
0 =Cohesionless
1 =Cohesive
2 =Rock
N50is the number of cycles necessary to degrade the soil by 50%.
SLRis the rate of loading for slow cyclic loading.
SWVSis the shear wave velocity for each soil layer.
SFDFis the fully degraded soil factor.
NSMOD0 = (default) - no soil gap, soil loads and unloads on the same curve.
1 =gap model, soil forms a gap when unloading parallel to initial stiffness in either tension or compression.
RSDAMPis the force proportional soil damping factor (lateral only)
(e.g. 0.01 applies 1% of the lateral soil force as a damping force)
TANDBspecify both Top and Bottom soil layer properties for the select layer.
0 =Use one set of properties per layer
1 =Specify top and bottom properties
SPT Sounding values (values are provided per soil set)
This input lines consists of two input lines per soil set. Contained in the first input line are the elevations for which
SPT data are available, and the SPT values are supplied in the second input line. In the event that no SPT values
are supplied for a given soil set, then two lines with 0 values will be supplied.
SPTELEV1 SPTELEV2
SPTVALUE1 SPTVALUE2
:
Where
SPTELEV1 is the highest elevation value for which an SPT blow count value is available.
SPTELEV2 is the second highest elevation value for which an SPT blow count value is available.
SPTVALUE1 is SPT blow count value corresponding to SPTELEV1.
SPTVALUE2 is SPT blow count value corresponding to SPTELEV2.
User defined P-Y data - ONLY FOR LSM=7
User defined soils require TWO additional lines of input.
(Two lines define the P-Y curve for the soil layer.)
Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8, Y9, Y10
P1, P2, P3, P4, P5, P6, P7, P8, P9, P10
Where
Yiis the ith Y value on the user specified P-Y curve.
Piis the ith P value on the user specified P-Y curve.
The user defined curves are specified by a set of TEN points.
User defined T-Z data - ONLY FOR ASM=5
User defined axial soil model requires TWO additional lines of input.
(Two lines define the T-Z curve for the soil layer.)
Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8, Z9, Z10
T1, T2, T3, T4, T5, T6, T7, T8, T9, T10
Where
Ziis the ith Z value on the user specified T-Z curve.
Tiis the ith T (axial stress) value on the user specified T-Z curve.
The user defined curves are specified by a set of TEN points.
User defined T- data - ONLY FOR TSM=2
User defined torsional soil model requires TWO additional lines of input.
(Two lines define the T- curve for the soil layer)
1, 2, 3, 4, 5, 6, 7, 8, 9, 10
T1,T2, T3, T4, T5, T6, T7, T8, T9, T10
Where
iis the ith value on the user specified T- curve.
Tiis the ith T (torque) value on the user specified T- curve.
When using the BENT header, a minimum of two material property lines are required. The first is for the pier cap
and the second is for the center section of the pier cap.
The MAST heading is used for high mast lighting/sign type structures.
MAST
N=N1 S=S1 H=H1O=O1 C=C1 B=B1, B2 W=W1 X=CC F=KFLOOD \
E=SPELEV D=CONT A=NUMLM, NUMPR T=TC, TCANT, TBEAM R=RH1, RH2, RH3 J=NLOPT
(all one line)
When using the MAST header, a minimum of two material property lines are required. The first is for the column,
the second is for the mast/sign portion. Next comes any additional properties (NUMPR). Next comes two lines for
high mast line loading, as follows:
LC= MLL
LA= ALL
Where
MLL is the line loading applied to the mast in the positive Y direction. One value is written on this line per load
case.
ALL is the line loading applied to the arm in the positive Y direction. One value is written on this line per load
case.
Next, the data for the additional members (NUMLM) is written.
The SOUND header is for use when sound walls are required.
SOUND
N=N1 S=S1 H=H1 O=O1 C=C1 B=B1, B2 W=W1 X=CC F=KFLOOD \
E=SPELEV D=CONT A=NUMLM, NUMPR T=TC, CANT, TBEAM J=NLOPT
The following line defines the soil layers behind the wall
O=IOPTI S=ISURG L=NLAYE
line 1
Where
IOPTIis equal to 1 for pressure at rest
is equal to 2 for active case computed with Coulomb expression
ISURG is 0 for no surcharge
is 1 for uniform surcharge
is 2 for line load
is 3 for strip load
NLAYEis the number of layers
This line defines the basic soil geometry
A=THETA S=BETA H= HWATE G=GWATE Q=Q1, Q2, Q3 line 2
Where
THETAis the inclination of the back of wall measured clockwise from horizontal plane (degrees)
BETAis the inclination of ground slope behind wall measured counterclockwise from the horizontal plane
(degrees)
HWATE is the Zcoordinate of ground water level (reference iscenter of pile cap)
GWATEis the unit weight of water
Q1, Q2, Q3are parameters for surcharge definition
If ISURG = 0Q1, Q2, Q3 are not used
If ISURG = 1Q1 = uniform surcharge
If ISURG = 2Q1 = line load intensity
Q2 = Horizontal distance of line load from back of wall
If ISURG = 3Q1 = Intensity of load
Q2 = Horizontal distance of load from back of wall
Q3 = Width of strip load
Soil Layer Property Lines (one line for each layer, NLAY)
T=THICK S=NSLAY P=COHES, PHI, DELTA G=GAMMA, GASAT
(one line per layer, the bottom layer being layer #1)
Where
THICKis the layer thickness
NSLAYis the number of sub-layers in which the layer will be divided
COHES is the cohesion of the soil
PHIis the friction angle of soil (degrees)
DELTAis the angle of friction soil/wall (degrees)
GAMMAis the unit weight of the soil
GASATis the saturated unit weight of the soil
One material property line is required when using the RETAIN header. Then any additional properties and extra
members.
The definition of the parameters for all structures are given below.
Where
N1is # of columns of the bridge bent supported on the pile group (INTEGER)
S1, S2, S3is spacing of the pier columns. For retaining walls and sound walls, S1 is the wall width. (REAL)
H1is height of the pier columns (REAL,)
O1is offset of the pile cap from the column (REAL)
C1is # of column nodes (INTEGER)
B1is # of pier cap nodes (Figure F21) (INTEGER)
B2is # of pier cap cantilever nodes (Figure F21) (INTEGER)
NPAD(L)is the number of bearing locations (left row of bearings)
POFF(L)is the offset from the first bearing location (left row of bearings)
PSPCx(L)is the bearing location spacing value (left row of bearings)
NPAD(R)is the number of bearing locations (right row of bearings)
POFF(R)is the offset from the first bearing location (right row of bearings)
PSPCx(R)is the bearing location spacing value (right row of bearings)
PUNF(L)is the uniform bearing location spacing (left row of bearings)
PUNF(R)is the uniform bearing location spacing (right row of bearings)
PUNFis uniform spacing between bearing locations, same for all locations (REAL)
For a single row of bearing locations, the left and right row parameters should be the same.
W1is cantilever length of top of bent (REAL)
CC is the cantilever flag, where a value of 1 means the model has cantilevers, and a value of 0 means the model
does not have cantilevers.
NUMLMis number of extra beam elements (Figure F21) (INTEGER)
NUMPRis number of extra beam properties (INTEGER)
TCis # of segments for tapered column (INTEGER), equal to zero for no tapered columns. This overrides C1
TCANTis # of segments for tapered cantilevers (INTEGER), equal to zero for no tapered cantilevers. This
overrides B2
TBEAM is the number of segments for tapered beams (INTEGER), equal to zero for no tapered beams
NLOPTselects the non-linear option for the pier structure analysis
NLOPT=1 for linear material
NLOPT=2 for nonlinear material
NLOPT=3 for linear material where interaction diagram are generated
RH1 is the depth of the pier cap at the cantilever base
RH2 is the depth of the pier cap at the center of the pier cap
RH3is the depth of the pier cap at the cantilever tip
KFLOODis the flag to tell if the column (if under the water table) is flooded or not. If flooded, the buoyancy will
use the net area. If not flooded, it will use the gross area (net=area-void).
SPELEVis the pier cap superelevation slope (+ or -) beginning at leftmost pier column. Expressed as a decimal
(not a percent).
CONTis the bridge span continuity option (over the pier)
0 for discontinous spans (does not transfer moment)(default)
1 for continous spans (transfers moment)
line 1
EPSC(1), EPSC(2),,,line 2
Stress-Strain Curve for Mild Steel, used with NLOPT=2 and MATOPT=2 and KSTEEL(1) = 1
S1=NPSC, SIGS(1), SIGS(2),,,line 1
EPSS(1), EPSS(2),,, y= yline 2
Stress-Strain Curve for Prestressing Steel, used with NLOPT=2 and MATOPT=2 and KSTEEL(2) = 1
S2=NPSC, SIGS(1), SIGS(2),,,line 1
EPSS(1), EPSS(2),,,line 2
Stress-Strain Curve for H-pile Steel, used with NLOPT=2 and MATOPT=2 and
KSTEEL(3) = 1
S3=NPSC, SIGS(1), SIGS(2),,,
EPSS(1), EPSS(2),,, y=y
line 1
line 2
Stress-Strain Curve for Tubular Steel, used with NLOPT=2 and MATOPT=2 and
KSTEEL(4) = 1
S4=NPSC, SIGS(1), SIGS(2),,,line 1
EPSS(1), EPSS(2),,, y=y
line 2
For Nonlinear Analysis of Square/Rectangular Piers, used with NLOPT=2 or 3 and KTYPE=2
W=WIDTH D=DEPTH V=DV B=BV N=ISTNOPT P=PREST S=SW
For nonlinear Analysis of Nonstandard Square/Rectangular Piers used with NLOPT=2, KTYPE=2, and
ISTNOPT= 2
NG=NGRPS HPI= IHPILE M=BMETH X=MINSPACE Z=TYPE
AS, Y, Z, PREST N=N1 D=D1
For Nonlinear Analysis of Round Piers, used with NLOPT=2 and KTYPE=1
NL=NLAY NG=NGRPS D=DP V=DV S=SW TH=DS IC=ICON, T=TR HPI=IHPILE \
BM=BMETH X=MINSPACE Z=TYPE
[PREST, NBS, D=DSI, A=ASI]
One of the next four lines is necessary for ICON ? 1(hoop or spiral steel is present)
FYH=FYHOOP HS=HOOPS N=NHOOP
or
FYH=FYHOOP HS=HOOPS D=DHOOP
or
FYS=FYSPI SP=SPIRS N=NSPI
or
FYS=FYSPI SP=SPIRS D=DSPI
For steel H-piles used with KTYPE=3 or HP=1 in either circular or square sections
Two lines are required:
OR=ORIENTline 1
[D=DEPTH U=WEIGHT] line 2, for standard H-pile sections
or
[D=DEPTH TW=WEB B=WIDTH TF=FLANGE] line 2, for user defined sections
For Nonlinear Analysis of Bullet Piers, used with NLOPT=2 or 3 and KTYPE=2
R=ORIENT W=WIDTH D=DEPTH T=VTYPE V=DV B=BV S=SW N=ISTNOPT P=PREST (line 1)
NG=NGRPS BM=BMETH X=MINSPACE Z=TYPE (line 2)
LOC is the location of the extra member in the model. A value of 1 means
the extra member connects column nodes, a value of 2
means the extra member connects pier cap nodes, and a value of 3 means the
extra member connects pile nodes.
TYPE is the cross section type code. A value of 0 means a pier type, and a
value of 1 means a pile type.
TAPERED COLUMN AND CANTILEVER SECTIONS
Columns and Cantilever Pier Cap sections can be set to tapered (non-prismatic) by setting TC and/or TCANT to
values greater than 0. When material properties, linear or non-linear, are set for tapered sections, 2 sets of
properties [base and top (tip)] are required instead of the one set required for prismatic sections. Figure: 17.3.y
and Figure: 17.3.z and sample inputs, below, illustrate the addition of tapered column and cantilever properties
to the input file.
top of column
These are load input lines. As many lines as needed can be used. One line must be supplied for each loaded
joint and each load condition. This can be skipped if no concentrated nodal loads are applied. This can happen
in the case of mast or sound walls where wind load is applied or in retaining walls where soil pressure is
applied.
Note, torsion in the pile cap can only be applied where piles are located.
In the input file after the LOAD header,
LOAD
NF, NL, NI, L=LC, F=FX, FY, FZ, MX, MY, MZ, D=DX, DY, DZ, RX, RY, RZ, T=TYPE (one line per nodal load)
Where
NFis the starting node number
NLis the ending node number(not available in GUI)
NIis the node numbering increment(not available in GUI)
LCis the load case number
FXis the force in the global X-direction
FYis the force in the global Y-direction
FZis the force in the global Z-direction
MXis the moment about the global X-axis
MYis the moment about the global Y-axis
MZis the moment about the global Z-axis
DX is the force in the global X-direction (only written for applied displacement)
DY is the force in the global Y-direction (only written for applied displacement)
DZ is the force in the global Z-direction (only written for applied displacement)
RX is the moment about the global X-axis (only written for applied displacement)
RY is the moment about the global Y-axis (only written for applied displacement)
RZ is the moment about the global Z-axis (only written for applied displacement)
TYPEis the load type specified in AASHTO (ignore for non-AASHTO loads)
LRFD Loads:
TYPE =
DCDead load of components
DDDowndrag
DWDead load of wearing surfaces and utilities
EHHorizontal earth pressure load
EVVertical earth pressure load
ESEarth surcharge load
EL Locked-in Construction Stresses
PS Post Tensioning
CR Creep
SHShrinkage
LLLive load
IMImpact
CEVehicular centrifugal force
BRVehicular braking force
PLPedestrian live load
LSLive load surcharge
WAWater Load
WSWind load on structure
WL Wind load on live load
VP Vertical Wind Pressure
VR Vertical Wind Press. - Reverse Loads
FRFriction
TUUniform temperature
TGTemperature gradient
SESettlement
EQEarthquake
ICIce load
CTVehicular collision force
CVVessel collision force
LFD Loads:
TYPE =
DDead load
LLLive load (AASHTO Type "L")
IMImpact (AASHTO Type "I")
EEarth pressure
BBuoyancy
WSWind load on structure (ASSHTO Type "W")
WLWind load on live load
LFLongitudinal force from live load
CFCentrifugal force
RRib shortening
SShrinkage
TTemperature
EQEarthquake
SFStream flow pressure
ICEIce pressure
The following information is used by piers with bearings locations. The information under the LOADBP header
describes the concentrated loads applied to the bearing locations.
LOADBP
PADNUM L= LC F= FX, FY, FZ, MX, MY, MZ T=TYPE B=DIR (one line per nodal load)
Where
PADNUM is the bearing number
LCis the load case number
FXis the force in the global X-direction
FYis the force in the global Y-direction
FZis the force in the global Z-direction
MXis the moment about the global X-axis
MYis the moment about the global Y-axis
MZis the moment about the global Z-axis
TYPEis the load type specified in AASHTO (ignore for non-AASHTO loads)
DIRis the bearing row ("L" for left or "R" for right)
:
This section must end with a blank line.
SSWINDL is the longitudinal wind pressure on the superstructure (at each angle)
PIERWINDT is the transverse wind pressure on the pier (at each angle)
PIERWINDL is the longitudinal wind pressure on the pier (at each angle)
LTWIND is the transverse wind line load on the live load (at each angle)
LLWIND is the longitudinal wind line load on the live load (at each angle)
Parameters with a strikethrough font are not currently used. These parameters were used by the previous wind
load generator (based on the AASHTO-LRFD 1997 Interim Revisions).
Note: This section must end with a blank line.
The wind load generator calculations are as follows:
Wind Load on Structure (WS)
l
Vertical loads at the bearing locations are determined using a rigid beam and spring model
Note: Since the wind load on the column is applied at the centroid (and not the pier cap), ratio of
clwarm/colheight is used to reduce the wind load in order to apply it at the level of the pier cap.
Vertical loads at the bearing locations are determined using a rigid beam and spring model
The following information is used by with multiple pier generation. The information under the PADBC header
describes the bearing location to superstructure connectivity. This information is provided per pier.
PADBC
L=LEFTPAD S= FX, FY, FZ, FRX, FRY, FRZ O=OFFSET
Or
R=RIGHTPAD S= FX, FY, FZ, FRX, FRY, FRZ O=OFFSET
Where
LEFTPADis the bearing location index number in left row of bearing locations
RIGHTPADis the bearing location index number in right row of bearing locations
FXis the fixity for the local x-direction
FYis the fixity for the local y-direction
FZis the fixity for the local z-direction
FRXis the fixity for rotation about the local x-axis
FRYis the fixity for rotation about the local y-axis
FRZis the fixity for rotation about the local z-axis
For all six directions: 0 for released (free), 1 for constrained
Values greater than 1 indicate the custom connection material property number. This custom connection is
described by a load-displacement relationship. See PADPROP header.
OFFSETis the bearing offset (measured from the centerline of the pier cap to the center of the bearing). This
value must be greater than zero when two rows of bearings are used.
This section must end with a blank line.
For a single row, the left and right bearing parameters should be the same.
Figure: 17.3.ad Bearing Connection Layout for One and Two Rows
The next line specifies the mass to be added to a node for each of the six global directions. There is one header
per pier.
NS,NF,NI M=MX,MY,MZ,MRX,MRY,MRZ
Where
NSis the starting node to add the mass to.
NFis the final node to add the mass to.(not written by GUI)
NIis the increment to generate additional node numbers at between NS and NF at which to add mass. (not
written by GUI)
MX\
MY
MZ
} are the mass values for the translational and rotational X,Y,Z directions
MRX |
MRY |
MRZ/
This section must end with a blank line.
MZ
} are the dampers values for the translational and rotational X,Y,Z directions
MRX |
MRY |
MRZ/
This section must end with a blank line.
:
where
The line is repeated for each load case.
L is load case
R = 0 corresponds to off
R = 1 corresponds to on
numPiers
Where
numPiersis the number of bridge piers
pierCoordX, pierCoordY, pierRot
Where
pierCoordXis the nodal x-coordinate for the pile cap origin (for that pier)
pierCoordYis the nodal y-coordinate for the pile cap origin (for that pier)
pierRotis the rotation angle about global z-axis (for that pier)
File: name.PLS
This file contains information for the bearing row to bridge span connectivity (per pier).
nodesLeft, nodesRight, spanNodeLeft, spanNodeRight, spanNodeLeftHeight, spanNodeRightHeight
Where
nodesLeftis the number of connection nodes for the left bearing row
nodesRightis the number of connection nodes for the right bearing row
spanNodeLeftis the connector node number for the begin of bridge span
spanNodeRightis the connector node number for the end of bridge span
spanNodeLeftHeightis the elevation above the pier cap (c.g.) for the begin of bridge span
spanNodeRightHeightis the elevation above the pier cap (c.g.) for the end of bridge span
(If there is a left row of bearings i.e. nodesLeft > 0)
(nodesLeft number of lines)
padLeftCoordX, padLeftCoordY, padLeftCoordZ
Where
padLeftCoordXis the nodal x-coordinate for the bearing connection node
padLeftCoordYis the nodal y-coordinate for the bearing connection node
padLeftCoordZis the nodal z-coordinate for the bearing connection node
(If there is a right row of bearings i.e. nodesRight > 0)
(nodesRight number of lines)
padRightCoordX, padRightCoordY, padRightCoordZ
Where
padRightCoordXis the nodal x-coordinate for the bearing connection node
padRightCoordYis the nodal y-coordinate for the bearing connection node
lmsdy,lpsdy,lmsto,lpsto,lmdiz,lpdiz,lmdix,lpdix,
lmdiy,lpdiy
lmsh2is the load case with the max pile shear-2
lpsh2is the pile number with the max pile shear-2
lmsp is the load case with the max pile shear-3
lpsp is the pile number with the max pile shear3
lmrm2 is the load case with the max pile moment-2
lprm2 is the pile number with the max pile moment-2
lmrm3 is the load case with the max pile moment-3
lprm3 is the pile number with the max pile moment-3
lmaxl is the load case with the max axial force
lpaxl is the pile number with the max axial force
lmtoris the load case with the max torsion
lptor is the pile number with the max torsion
lmsaxis the load case with the max soil axial force
lpsaxis the pile number with the max soil axial force
lmsdx is the load case with the max soil lateral-x force
lpsdx is the pile number with the max soil lateral-x force
lmsdy is the load case with the max soil lateral-y force
lpsdy is the pile number with the max soil lateral-y force
lmsto is the load case with the max soil torsion
lpsto is the pile number with the max soil torsion
lmdizis the load case with the max pile axial-displacement
lpdiz is the pile number with the max pile axial-displacement
lmdix is the load case with the max pile x-displacement
lpdix is the pile number with the max pile x-displacement
lmdiy is the load case with the max pile y-displacement
lpdiy is the pile number with the max pile y-displacement
MINIMUMS
lmsh2,lpsh2,lmsp,lpsp,lmrm2,lprm2,lmrm3,lprm3,
lmaxl,lpaxl,lmtor,lptor,lmsax,lpsax,lmsdx,lpsdx,
lmsdy,lpsdy,lmsto,lpsto,lmdiz,lpdiz,lmdix,lpdix,
lmdiy,lpdiy
lmsh2is the load case with the min pile shear-2
lpsh2is the pile number with the min pile shear-2
lmsp is the load case with the min pile shear-3
lpsp is the pile number with the min pile shear3
lmrm2 is the load case with the min pile moment-2
lprm2 is the pile number with the min pile moment-2
lmrm3 is the load case with the min pile moment-3
lprm3 is the pile number with the min pile moment-3
lmaxl is the load case with the min axial force
The next two sections are repeated twice and both are repeated NUMLC times, for each load case.
Mtype, nume
Mtypeis the element type (=3 for structure, =2 for piles)
Numeis the number of elements
Axial
Axialis the axial force for the member for the appropriate load case.
The next two sections are repeated twice and both are repeated NUMLC times, for each load case.
Mtype, nume
Mtypeis the element type (=3 for structure, =2 for piles)
Numeis the number of elements
Rmom
Rmomis the maximum moment in the member for the appropriate load case.
Do I=1,numlc
Do j=1,8 (the eight sets of results)
Read() (stress(k), k=1, 8* #elements)
Enddo
Enddo
The next set of records is repeated for each cross section for which capacity information exists.
Nlcvis the number of contour slices for this cross section
PTUV, YPC, ZPC, YPT. ZPT
Ptuvis the ultimate axial tension strength
Ypcis the y shift for the plastic centroid (compression)
Zpcis the z shift for the plastic centroid (compression)
Ypt is the y shift for the plastic centroid (tension)
Zpt is the z shift for the plastic centroid (tension)
(PNC(J)J=1, 13) (repeated nlcv times)
pncis the table of capacity results. Where the values are:
pnc(1) = * Compression capacity
pnc(2) = * moment capacity about local 3 axis (M1)
pnc(3) = * moment capacity about negative local 2 axis (M2)
pnc(4) = * moment capacity about negative local 3 axis (M3)
pnc(5) = * moment capacity about local 2 axis (M4)
pnc(6) = 1
pnc(7) = 1
pnc(8) = 2
pnc(9) = 2
pnc(10)=3
pnc(11)=3
pnc(12)=4
pnc(13)=4
The and s are used as a pair for the following capacity equation:
If the compression is in the 1st quadrant (+2,+3) then use M1, M2, 1, 1
If the compression is in the 2nd quadrant (-2,+3) then use M3, M2, 2, 2
If the compression is in the 3rd quadrant (-2,-3) then use M3, M4, 3, 3
If the compression is in the 4th quadrant (+2,-3) then use M1, M4, 4, 4
XMLis the distance from the I end where the maximum midspan moment exists.
FRATIis the capacity ratio at the I end.
FRATJis the capacity ratio at the J end.
AXLIis the axial force at the I end of the member.
AXLJis the axial force at the J end of the member.
NOTE: All values are single precision real numbers (4 bytes). Also, the pile elements come first, then the
structure elements.
NPlt
Where
NPltis the results version number (currently version 1)
Ndynam
Where
Ndynamis the type of analysis (0 for static, 1 for dynamic, 2 for response spectrum analysis)
NTimeStep
Where
NTimeStepis the time step used (for time domain dynamic analysis, otherwise 0)
NVEC
Where
NVECis the number of eigenvectors (for response spectrum analysis, otherwise 0)
File: name.EIG
This file contains eigenvalues (frequencies) and eigenvectors (mode shapes) used in the response spectrum
analysis.
NVEC
Where
NVECis the number of number of eigenvectors
NNODE
Where
NNODEis the number of nodes in the model
FREQ1, FREQ2, FREQ3, .. FREQN
Where
FREQxis the vibration frequency for mode x.
Loop over the number of eigenvectors and over each node in the model
NODE, PHIXX, PHIYY, PHIZZ, PHIRX, PHIRY, PHIRZ
Where
NODEis the model node number (integer)
PHIXXis the eigenvector in the x-direction (double)
PHIYYis the eigenvector in the y-direction (double)
PHIZZis the eigenvector in the z-direction (double)
PHIRXis the eigenvector about the x-axis (double)
PHIRYis the eigenvector about the y-axis (double)
PHIRZis the eigenvector about the z-axis (double)
File: name.ASH
This file contains design code and limit state information.
nCodeType
Where
nCodeTypeis the AASHTO design code used for load combinations
(0 for LRFD, 1 for LFD)
nGroup1, nGroup2, , nGroup11
Where
nGroup1are the limit states that were analyzed
(0 for analyzed, 1 for not analyzed)
(CritPl(J),J=1,11)
Where
CritPlis load combination number with the maximum pile demand/capacity ratio for each analyzed limit state (0
if not analyzed)
(CritCol(J),J=1,11)
Where
CritPlis load combination number with the maximum pier column demand/capacity ratio for each analyzed limit
state (0 if not analyzed)
(CritPierCap (J),J=1,11)
Where
CritPlis load combination number with the maximum pier cap demand/capacity ratio for each analyzed limit
state (0 if not analyzed)
References
Andrade, P., "Materially and Geometrically Non-Linear Analysis of Laterally Loaded Piles Using a Discrete
Element Technique," Masters Report, Universtiy of Florida, Gainesville, 1994.
Bamford, W. E. (1969). Anisotropy, and the natural variability of rock properties, Proc. Symp. Rock Mech.,
Sydney, 1-10.
Bowles, J. E. (1987). Elastic Foundation Settlement on Sand Deposits, Journal of Geotechnical Engineering,
ASCE, Vol. 113, No. 8, pp. 846-860.
Bowles, J. E., "Foundation Analysis and Design", McGraw-Hill, New York, 1977.
Brown, D., Morrison, C., and Reese, L. (1988), "Lateral Load Behavior of a Pile Group in Sand," ASCE Journal of
Geotechnical Engineering, Vol. 114, No. 11, pp. 1261-1276.
Brown, E. T. (1970). Strength of models of rocks with intermittent joints, J. Soil Mech. Fdn. Eng., ASCE, 96,
1935-9.
California Department of Transportation (CALTRANS), Bridge Design Practice Manual, Sacrament, CA, 1981.
Chai, Y., H., Preistly, M., J., N., and Seible, F. (1991), "Flexural Retrofit of Circular Reinforced Bridge Columns by
Steel Jacketing,", University of California, San Diego, La Jolla.
Consolazio et. al. (2006), "Barge Impact Testing of the St. George Island Causeway Bridge, Phase III", University
of Florida, Final report, FDOT Contract No. BC-354 RPWO 76.
Das, B. M. (1980). A Procedure for estimation of Ultimate Uplift Capacity of Foundations in Clay, Soils and
Foundations, Vol . 20, No. 1 pp. 77-82.
Das, B. M., and Jones, A. D. (1982). Uplift Capacity of Rectangular Foundations in Sand, Transportation
Research Record 884, National Research Council, Washington, D. C., pp. 54-58.
De Beer, E. E. (1970). Experimental Determination of the Shape Factors and Bearing Capacity Factors of Sand,
Geotechnique, Vol. 20, No. 4, pp. 387-411.
FHWA, PCSTABLS-Users Manual, Federal Highway Administration, Springfield, VA 1988.
Filippou, F. C., Popov, E. P., Bertero, V. V. (1983). "Effects of Bond Deterioration on Hysteretic Behavior of
Reinforced Concrete Joints". Report EERC 83-19, Earthquake Engineering Research Center, University of
California, Berkeley.
Gazioglu, S. M., and ONeill, M. W. "Evaluation of P-Y Relationships in Cohesive Soils," from Analysis and
Design of Pile Foundations, proceedings of a symposium sponsored by the ASCE Geotechnical Engineering
Division, ASCE National Convention, San Francisco, CA, pp. 192-213.
Georgiadis, M. "Development of P-Y curves for Layered Soils," Proceedings, Geotechnical Practice in Offshore
Engineering, American Society of Civil Engineers, pp. 536-545.
Getter, D. J., and Consolazio, G. R., (2011). Relationships of barge bow force-deformation for bridge design:
Probabilistic consideration of oblique impact scenarios. Transportation Research Record 2251, Transportation
Research Board, Washington, DC, 3-15, 2011.
Hanna, A. M., and Meyerhof, G. G. (1981). Experimental Evaluation of Bearing Capacity of Footings Subjected
to Inclined Loads, Canadian Geotechnical Journal, Vol. 18, No. 4, pp. 599-603.
Hannigan, P. J., Goble, G. G., Likins, G. E., and Rausche, F., (2006). "Design and Construction of Driven Pile
Foundations Volume I and II," FHWA Publication No. FHWA-NHI-05-042. U.S. Department of Transportation,
Federal Highway Administration, Washington D.C.
Hoek, E. and Brown, E. T. (1980). Empirical strength criterion for rock masses, J. Geotech. Engng. Div., ASCE,
106 (GT 9), p. 1013-1035.
Hoek, E. and Brown, E. T. (1988). The Hoek-Brown failure criterion a 1988 update , Rock Engineering for
Underground Excavations, Proc. 15th Canadian Rock Mech. Symp. J. C. Curran, Ed. Toronto: Dept. Civil
Engineering, University of Toronto, p. 31-38.
Hoek, E. and Brown, E. T. (1997). Practical estimates or rock mass strength , Int. J. Rock Mech. & Mining Sci. &
Geomechanics Abstracts. 34(8), p. 1165-1186.
Itani, A. M. "Future Use of Composite Steel-Concrete Columns in Highway Bridges." AISC Engineering
Journal.33, No.3 pp. 110-115 1996.
Janbu, N., Bjerrum, L., and Kjaernsli, B. (1956). Veiledning vedlosning av fundamenteringsoppgaver,
Publication No. 18, Norwegian Geotechnical Institute, pp. 30-32.
Kim, Myoung-Ho, "Analysis of Osterburg and Stanamic Axial Load Testing and Conventional Lateral Load
Testing", Masters Thesis, University of Florida, Gainsvelli, Florida, 2001.
Kulhawy, F. and Mayne, P. "Manual for Estimating Soil Properties for Foundation Design." Electric Power
Research Institute (EPRI) Report. EPRI EL-6800. Project 1493-6. Aug. 1990, pp. 5-17.
Kumbhojkar, A. S. (1993). Numerical Evaluation of Terzaghis Ny, Journal of Geotechnical Engineering,
American Society of Civil Engineers, Vol. 119, No. 3, pp. 598-607.
Mandel, J., and Salencon, J. (1972). Force protante dun sol sur une assise rigide (tude thorique),
Geotechnique, Vol. 22, No. 1, pp. 79-93.
Mander, J. B., Priestley, M. J. N.,Park R., "Theoretical Stress Strain Model for Confined Concrete" ASCE Journal
of Structural Engineering 114,pp. 1804-1826, 1988.
Mander, J. B., Priestley, M. J. N.,Park R., "Observed Stress-Strain Behavior of Confined Concrete" ASCE Journal
of Structural Engineering 114,pp. 1827-1849, 1988.
Matlock, H. "Correlations for Design of Laterally Loaded Piles in Soft Clay," Paper No. OTC 1204, Proceedings,
Second Annual Offshore Technology Conference, Houston, Texas, Vol. 1, 1970, pp. 577-594.
McVay,M. C., O'Brien, M., Townsend, F. C., Bloomquist, D. G., and Caliendo, J. A. "Numerical Analysis of
Vertically Loaded Pile Groups," ASCE, Foundation Engineering Congress, Northwestern University, Illinois, July,
1989, pp. 675-690.
McVay,M. C., M., Niraula, L. "Development of Modified T-Z Curves for Large Diameter Piles/Drilled Shafts in
Limestone for FB-Pier," Report Number 4910-4504-878-12, National Technical Information Service, Springfield,
VA June 2004.
Meyerhoff, G. G. (1953). The Bearing Capacity of Foundations Under Eccentric and Inclined Loads,
Proceedings, Third International Conference on Soil Mechanics and Foundation Engineering, Zrich, Vol. 1, pp.
440-445.
Meyerhoff, G. G. (1963). Some Recent Research on the Bearing Capacity of Foundations, Canadian
Geotechnical Journal, Vol. 1, No. 1, pp. 16-26.
Meyerhoff, G. G. (1974). Ultimate Bearing Capacity of Footings on Sand Layer Overlying Clay, Canadian
Geotechnical Journal, Vol. 11, No. 2, pp. 224-229.
Meyerhoff, G. G. and Adams, J. I. (1968). The Ultimate Uplift Capacity of Foundations, Canadian Geotechnical
Journal, Vol. 5, No. 4, pp. 225-244.
Meyerhoff, G. G. and Chaplin, T. K. (1953). The Compression and Bearing Capacity of Cohesive Soils, British
Journal of Applied Physics, Vol. 4, pp. 20-26.
Meyerhoff, G. G. and Hanna, A. M. (1978). Ultimate Bearing Capacity of Foundations on Layered Soil under
Inclined Load, Canadian Geotechnical Journal, Vol. 15, No. 4, pp. 565-572.
Meyerhoff, G. G. (1965). Shallow Foundations, Journal of Soil Mechanics and Foundations Division, American
Society of Civil Engineers, Vol. 91, No. SM2, pp. 21-31
Mitchell, J. S., "A Nonlinear Analysis of Bi-Axially Loaded Beam-Columns Using a Discrete Element Model," Ph.
D. Dissertation, The University of Texas at Austin, August, 1973.
Mirza, S., A., and MacGregor, J., G., "Variability of mechanical properties of reinforcing bars,", ASCE Journal of
Structural Engineering, 105(ST5):921-937, May 1979.
Murchison, J. M. and ONeill, M. W. "Evaluation of P-Y Relationships in Cohesionless Soils," from Analysis and
Design of Pile Foundations, proceedings of a symposium sponsored by the ASCE Geotechnical Engineering
Division, ASCE National Convention, San Francisco, CA, pp. 174-191.
ONeill, M. W. and Dunnavant, T. W., "A Study of Effect of Scale, Velocity, and Cyclic Degradability on Laterally
Loaded Single Piles in Overconsolidated Clay.", Rep. No. UHCE 84-7, Dept. of Civil Engineering, Univ. of
Houston, Houston, TX, 1984.
ONeill, M.W., Townsend, F. C., Hassan, K. M., Buller, A., and Chan, P. S., "Load Transfer for Drilled Shafts in
Intermediate Geomaterials," Report No. FHWA-RD-95-172, November 1996.
Peck, R. B., W. E. Hanson, and T. H. Thornburn (1953). Foundation Engineering, John Wiley and Sons, New
York.
Peck, R.B., Hanson, W. E. and Thornburn, T. H., "Foundation Engineering", John Wiley & Sons, 1974.
Pinder, Terrence, "A Model for Concrete Under The Effect of Transverse Confinement", Report presented to the
graduate committee of the Department of Civil Engineering, University of Florida, Summer 1997.
Precast Conrete Institute, PCI Design Handbook, 4th Ed., Chicago IL, 1992.
Randolph, M.F., "Piles Subjected to Torsion," Journal of the Geotechnical Division, ASCE, Vol. 107, No. GT8,
August, 1981, pp. 1095-1111.
Reese, L. C., Cox, W. R. and Koop, F. D "Analysis of Laterally Loaded Piles in Sand," Paper No. OTC 2080,
Proceedings, Fifth Annual Offshore Technology Conference, Houston, Texas, 1974 (GESA Report No. D-75-9).
Reese, L. C., Cox, W. R. and Koop, F. D. "Field Testing and Analysis of Laterally Loaded Piles in Stiff Clay,"
Paper No. OTC 2312, Proceedings, Seventh Offshore Technology Conference, Houston, Texas, 1975.
Reese, L. C. and Welch, R. C. "Lateral Loading of Deep Foundations in Stiff Clay," Journal of the Geotechnical
Engineering Division, American Society of Civil Engineers, Vol. 101, No. GT7, Proceedings Paper 11456, 1975,
pp. 633-649 (GESA Report No. D-74-10).
Sayed, S. M. and Bakeer, R. M. "Efficiency Formula For Pile Groups," Journal of the Geotechnical Engineering,
American Society of Civil Engineers, Vol. 118, No. 2, Paper No. 26553, 1992, pp. 278-299.
Scott, B., D., Park R., and Priestly, M., J., N., "Stress-strain Behavior of concrete confined by overlapping hoops at
low and high strain rates,", ACI Journal , 79(1): 13-27, Jan./Feb 1982.
Stoll, U.W., "Torque Shear Test of Cylindrical Friction Piles," Civil Engineering, ASCE, Vol. 42, No. 4, April., 1972,
pp.63-64.
Stone, W., C. and Cheok, G., S., "Inelastic Behavior of Full-Scale Bridge Columns Subjected to Cyclic Loading".
Report No. NIST/BSS-166, National Institute of Standards and Technology, U.S. Department of Commerce,
Gaithersberg, MD 20899, Jan. 1989.
Sudip Basack, (2007) "Experimental Investigation of Pile Group under Lateral Cyclic Load in Layered Soil", Dept.
of Applied Mechanics, Bengal Engineering and Science University
Terzaghi, K. (1943). Theoretical Soil Mechanics, Wiley, New York.
University of Florida, "Development of Modified T-Z Curves for Large Diameter Piles/Drilled Shafts in Limestone
for FB-Pier", UF Project #4910-4504-878-12.
Vesic, A. S. (1963) Bearing Capacity of Deep Foundation in Sand. Highway Research Record Mo. 39, National
Academy of Sciences, pp. 112-153.
Vesic, A. S. (1973) Analysis of Ultimate Loads of Shallow Foundations. Journal of Soil Mechanics and
Foundations Division, American Society of Civil Engineers, Vol. 99, No. SM1, pp. 45-73.
Wang, S. T., and Reese, L. C., "COM624P Laterally loaded pile analysis for the microcomputer, ver. 2.0"
FHWA-SA-91-048, Springfield, VA, 1993.
Welch, R. C., and Reese, L.C., "Lateral Load Behavior of Drilled Shafts", Research Rep. No. 3-5-65-99,
conducted for Texas Highway Department and U.S. Department of Transportation, Federal Highway
Administration, Bureau of Public Roads, Center for Highway Research, Univ. of Texas at Austin, TX, 1972.
Yassin, M., "Nonlinear analysis of prestressed concrete structures under monotonic and cyclic loads", PhD.
dissertation, University of California, Berkeley, 1994.
Tutorials
The tutorials below require an Internet connection and are located on the BSI server. (Click Link to download)
Element Forces
Covers the Element Forces Dialog.
Help About
The Help About displays important program information.
License Update
Explains how to update a license in FB-MultiPier.
Max Min
Covers the Max Min Dialog.
Percentage Steel
Explains how to use the Percentage Steel method to enter reinforcing steel.
Pile Sets
Explains how to create Pile Sets.
Printable Forces
Covers the Printable Forces Dialog.
Soil Plot
Covers the Soil Plot Dialog.
Soil Table
Demonstrates how to enter data in the Soil Table.
Zoom Feature
Covers the Zoom features in FB-MultiPier.