Example 9 Eng
Example 9 Eng
Example 9. Analysis of structure on soil foundation with the use of SOIL system
Description:
Reinforced concrete floor slab with dimensions in column axes 4 x 6 m, thickness 150 mm.
Reinforced concrete base slab with dimensions in column axes 4 x 6 m, thickness 500 mm.
Reinforced concrete columns of rectangular section with dimensions 400 x 800 mm.
Height of the framework – 3 m.
Analysis is carried out for floor slab with dimensions of finite elements 0.2 x 0.2 m and base slab with
dimensions of FE 0.4 x 0.4 m and 0.4 x 0.2 m (see Fig.9.1).
Loads:
– load case 1 – dead weight;
2
– load case 2 – dead uniformly distributed load g1 = 0.5 t/m applied to floor slab;
2
dead uniformly distributed load g2 = 1 t/m and dead concentrated vertical load P = 100t applied to base
slab (see Fig.9.2, load case 2);
– load case 3 – concentrated horizontal load F = 2t applied to floor slab (see Fig.9.2, load case 3).
On the taskbar, click the Start button, and then point to All Programs. Point to the folder that
contains LIRA SAPR / LIRA-SAPR 2015 and then click LIRA-SAPR 2015.
Click OK .
It is also possible to open the Model type dialog box with a pre-defined type of model. To do this,
on the LIRA-SAPR menu (Application menu), point to New and click Model type 5 (Six DOF per
node) command . One more way to do the same: on the Quick Access Toolbar, click New and
in the drop-down menu select Model type 5 (Six DOF per node) command . Then you should
define only problem name.
To save all output data files for the problem in certain directory, select appropriate check box. The
directory name will coincide with the name of the problem. This directory will appear in the directory for
files with analysis results. This is helpful if you have to find output data files for certain problem, then
transfer these files or review and evaluate them with the help of Windows Explorer or other file
managers.
Click Apply .
Click Apply .
button .
In the Display dialog box (see Fig.9.5), select the Node numbers check box on the Nodes tab.
Click Redraw .
To add columns:
On the Select toolbar (by default, it is displayed at the bottom of the screen), click PolyFilter .
In the PolyFilter dialog box, select the Filter for nodes tab (the first tab) (see Fig.9.6).
Select By numbers of nodes check box and specify the following node numbers (separated either
by comma or by space character): 43, 58, 263, 278, 398, 428, 1138, 1168.
Click Apply .
Fragmentation .
On the Create and edit ribbon tab, on the Create panel, point to Add element drop-down list and
click Add bar . The Add element dialog box is displayed with the Add bar tab open (see
Fig.9.7).
To add bar elements between nodes, specify with the pointer the following pairs of nodes in
sequence: No. 43 and 398, 58 and 428, 263 and 1138, 278 and 1168 (in this case the rubber-band
line is automatically stretched between the nodes that you select).
To restore design model in initial view after fragmentation, on the Select toolbar, click Restore model
fragmentation button .
On the Create and edit ribbon tab, on the Stiffness and restraints panel, click Perfectly rigid body
button .
On the Select toolbar, point to Select nodes and click Select nodes button .
Select with the pointer nodes where base slab is connected to perfectly rigid body 1 (column body)
No. 22-24, 42-44 and 62-64 (nodes will be coloured red).
You can select elements either with a single click or by dragging selection window around
appropriate elements.
If required, it is possible to enlarge the fragment of the model with the Zoom command (located on the
Select toolbar) or with the mouth wheel button.
In the Perfectly rigid body dialog box (see Fig.9.8), in the Assign principal node of current PRB
area, select the Specify principal node with pointer check box.
Then, on the model of base slab, specify with the pointer node No.43 (node will be coloured black).
To generate the first perfectly rigid body, in the Define perfectly rigid body area, click Add PRB
button .
All operations with rigid body are carried out only with their principal nodes, e.g. restraints, local
(nodal) coordinate system, initial displacement, coupled degrees of freedom.
Make sure that the Specify principal node with pointer check box is selected and specify with the
pointer node No.58 on the base slab.
On the Select toolbar, point to Select nodes and click Select nodes button .
Select with the pointer nodes where base slab is connected to perfectly rigid body 2 - nodes No.37-
39, 57-59 and 77-79.
To generate the second perfectly rigid body, in the Define perfectly rigid body area, click Add PRB
button .
When the Select nodes command is active, select with the pointer nodes where base slab is
connected to perfectly rigid body 3 - nodes No.242-244, 262-264 and 282-284.
In the Perfectly rigid body dialog box, in the Assign principal node of current PRB area, select
the Specify principal node with pointer check box.
Then, on the model of base slab, specify with the pointer node No.263.
To generate the third perfectly rigid body, in the Define perfectly rigid body area, click Add PRB
button .
When the Specify principal node with pointer check box is selected, on the model of base slab,
specify with the pointer node No.278.
On the Select toolbar, point to Select nodes and click Select nodes button .
Select with the pointer nodes where base slab is connected to perfectly rigid body 4 - nodes No.257-
259, 277-279 and 297-299.
To generate the fourth perfectly rigid body, in the Define perfectly rigid body area, click Add PRB
button .
To restore design model in initial view after fragmentation, on the Select toolbar, click Restore model
fragmentation button .
On the Select toolbar, point to Select nodes and click Select nodes button .
Select with the pointer nodes where floor slab is connected to perfectly rigid body 5 - nodes No.359-
363, 396-400 and 433-437.
In the Perfectly rigid body dialog box, in the Assign principal node of current PRB area, select
the Specify principal node with pointer check box.
Then, on the model of floor slab, specify with the pointer node No.398.
To generate the fifth perfectly rigid body, in the Define perfectly rigid body area, click Add PRB
button .
To restore design model in initial view after fragmentation, on the Select toolbar, click Restore model
On the Select toolbar, point to Select nodes and click Select nodes button in order to make this
command not active.
button .
In the Stiffness and materials dialog box (see Fig.9.9a), click Add. The list of standard section types
will be presented in the Add stiffness dialog box (see Fig.9.9b).
a b
Figure 9.9 Dialog boxes: a – Stiffness and materials, b – Add stiffness
Double-click the Rectangular bar icon in the list. The Define standard section dialog box opens. In
this dialog box you can define material properties for selected type of the section.
In the Define standard section dialog box specify the following parameters for Rectangular bar
(see Fig.9.10):
modulus of elasticity – E = 3e6 t/m2 (for the U.S. keyboard layout);
geometric properties – B = 40 cm; H = 80 cm;
unit weight of material – Ro = 2.75 t/m3.
To preview schematic presentation, click Draw.
In the Stiffness and materials dialog box, select the third tab Plates, solids, numerical and double-
click the Plates icon in the list.
In the Stiffness for plates dialog box (see Fig.9.11), specify the following parameters for Plate (floor
slab):
modulus of elasticity – E = 3e6 t/m2 (for the U.S. keyboard layout);
Poisson's ratio – ν = 0.2;
thickness – H = 15 cm;
unit weight of material – Ro = 2.75 t/m3.
In the Stiffness and materials dialog box, in the List of stiffness types, select ‘2.Plate H 15’.
Click Copy.
In the Stiffness and materials dialog box, in the List of stiffness types, select ‘3.Plate H 15’.
Click Edit.
In another Specify stiffness for plates dialog box specify parameter for base slab:
thickness - H = 50 cm.
Click OK .
To hide the library of stiffness parameters, click Add in the Stiffness and materials dialog box.
You can select elements either with a single click or by dragging selection window around
appropriate elements. When the Select vertical bars command is active, you can drag selection
window around the whole model and only vertical bars will be selected.
In the Stiffness and materials dialog box, click Apply . The elements become unselected. It
indicates that the current combination of stiffness type and material is assigned to selected elements.
In the Stiffness and materials dialog box, on the Properties tab, in the List of stiffness types, click
the stiffness type 2. Plate H15.
Click Set as current type.
In the Stiffness and materials dialog box, on the Properties tab, in the List of stiffness types, click
the stiffness type 3. Plate H50.
Click Set as current type.
To avoid geometric unstability in the XOY-plane, additional boundary conditions are defined for
base slab (with the help of one-node elements FE 56).These elements are defined at nodes of the
base slab. The stiffness for all FE 56 will be defined as equal to the following value: 70% of stiffness of
elastic foundation C1z multiplied to load area associated with one FE 56.
As the stiffness of elastic foundation is unknown, we accept initial values of stiffness for FE 56 as
equal to 200 t/m. When moduli of subgrade reaction are calculated, the stiffness value for FE 56 may
be clarified.
button .
In the Stiffness and materials dialog box, click Add. The dialog box expands to display the library of
stiffness parameters. In the Add stiffness dialog box, select the Plates, solids, numerical tab (the
third tab).
Double-click the Numerical for FE 56 icon in the list. The Numerical description for FE 56 dialog
box opens.
In the Numerical description for FE 56 dialog box (see Fig.9.12), specify the following parameters
for the section:
stiffness of FE per unit length in tension-compression along the global X-axis – Rx = 200 t/m;
stiffness of FE per unit length in tension-compression along the global Y-axis – Ry = 200 t/m.
To confirm the specified data, click OK.
To add FE 56:
When the Select block button is active, select with the pointer any node or element of the base
slab.
On the Create and edit ribbon tab, on the Create panel, point to Add element drop-down list and
Click Apply .
When FE 56 are added to design model, stiffness for these FE was assigned as the current one
in the Stiffness and materials dialog box. This stiffness was automatically assigned to these added
FE.
Select with the pointer any node or element of the base slab.
On the Create and edit ribbon tab, on the Stiffness and restraints panel, click Moduli of subgrade
reaction button .
In the Define moduli C1 and C2 dialog box (see Fig.9.14), make sure that the Plates check box is
selected.
To define moduli of subgrade reaction, select the From soil model option and define:
uniformly distributed load on soil Рz = 12 t/m2.
Click Apply .
To unselect all nodes, on the Select toolbar, click Unselect all button .
On the calculate C1, C2 tab, define the following parameters that will be applied for calculation of
moduli of subgrade reaction:
in the Data for calculation area, define coefficient for depth of compressible stratum as equal
to 0.5;
in the Merge loads area, clear the Replace adjacent loads or loads of neighbouring values
with one load equal to average value check box.
Then in the Soil model dialog box, click Attach soil model button.
You could also open Sol model dialog box by clicking the Soil model button (Advanced
edit options ribbon tab, SOIL panel).
In the Open file with soil model dialog box (see Fig.9.16), when the file name Example9 is
mentioned, click Open.
In the SOIL system window, on the MODEL menu, click Grids (button on the toolbar).
In the Grids dialog box (see Fig.9.17), select the First point check box and define coordinates (m):
X = -1; Y = -1.
Then define the step value for the grid:
Along X-axis Along Y-axis
Step value Steps Step value Steps
1 18 1 14.
Click Apply .
Close the Grids dialog box.
To calculate moduli of subgrade reaction we will use soil properties that were defined according
to results of geological engineering survey and are located in the Soil properties table. To activate
this table, on the MODEL menu, click Soil properties (button on the toolbar).
Click Apply .
Click Apply .
Click Apply .
Close the Boreholes dialog box.
In the Loads dialog box (see Fig.9.19), select the type of load – Arbitrary load (the second tab).
Select the Elevation check box and specify the value equal to 97 (m).
For the specified intensity of soil pressure as 20 t/m2, click the Specify load (arbitrary polygon) on
plan button .
Then specify (with the pointer) load contours on the grid of soil as presented in Fig.9.20 (beginning
from the point with coordinate X = 16, Y = 12 and move the pointer 8m left, then 4m back, then 4m
right, 5m back, 4m right and come back to initial point).
Click Apply .
To save the data about design model, on the FILE menu, click Save (button on the toolbar).
Close the SOIL system window and return to the main window of LIRA-SAPR.
To change the number of the current load case, click the Next load case button located on the
Status bar or on the toolbar.
On the Select toolbar, point to Select nodes drop-down list and click Select nodes button .
Select nodes No.106, 115, 206 and 215 with the pointer.
On the Create and edit ribbon tab, select the Loads panel, then select Load on nodes command
Click OK .
In the Define loads dialog box, select the fourth tab Load on plates .
Click OK .
In the Define loads dialog box, click Uniformly distributed load button once again.
2
In the Load parameters dialog box specify P = 0.5 t/m .
Click OK .
To change the number of the current load case, click the Next load case button located on the
Status bar or on the toolbar.
On the Select toolbar, point to Select nodes drop-down list and click Select nodes button .
Then select node No.321 (the nearest left node on the floor slab).
In the Define loads dialog box, select the second tab Load on nodes and define direction along
the X-axis.
Click OK .
On the Create and edit ribbon tab, on the Loads panel, click Edit load cases . The Edit load
cases dialog box is displayed on the screen (see Fig.9.27).
For load case 1 – in the Edit selected load case area, in the Type box, select Dead and click Apply
.
For load case 2 – in the Edit selected load case area, in the Type box, select Dead and click Apply
.
For load case 3 – in the Edit selected load case area, in the Type box, select Short-term and click
Apply .
On the Analysis ribbon tab, select the More calculations panel and click DCL button .
As the type of load cases was defined in the Edit load cases dialog box (see Fig.9.27), the DCL
table is generated automatically with parameters accepted by default for every load case. Now you
have to modify parameters for the third load case and define combinations.
In the Design combinations of loads dialog box (see Fig.9.28), select building code SNIP 2.01.07-
85* and for load case 3, in the Load factor cell, define factor as equal to 1.4.
To carry out static analysis, on the Analysis ribbon tab, select the Analysis panel and in the
When analysis procedure is complete, to review and evaluate analysis results, select the Results
ribbon tab.
On the Select toolbar, click Flags of drawing button . In the Display dialog box, clear the
Loads check box on the General tab.
Click Redraw .
In the mode of analysis results visualization, by default design model is presented with account of
nodal displacements (see Fig.9.30).
To display diagram Qz , on the Results tab, select Forces in bars panel and click Shear force
panel, select the Stress mosaic plot command in the Stress mosaic/contour plots drop-down
list.
To present stress mosaic plot for Nx, click Stress Nx button on the same panel.
To present stress mosaic plot for Rz (soil pressure), click Stress Rz button on the same panel.
To present the full picture of stress mosaic plots for Rz in base slab, select the slab with the Select
block command and perform fragmentation.
To restore design model in initial view, on the Select toolbar, click Restore model .
Elements located inside the PRB are not considered in analysis. That's why, forces in such
elements are absent with the exception of soil pressure.
Elements of floor slab that are located inside the PRB may be deleted from design model because
they have no influence on analysis of the framework.
Mosaic plot Pz .
click Apply .
To visualize analysis results by DCL, click Results by DCL button on the Status bar.
To present mosaic plot Rz (soil pressure), on the Results ribbon tab, on the Stress in plates and
Make sure that Stress mosaic plot Rz button is active for the second DCL. Then on the
Results ribbon tab, on the Tools panel, click Apply soil pressure button .
In the Apply soil pressure dialog box (see Fig.9.32), when All elements option is selected, click
OK.
The system automatically goes back to the mode of creating design model and analysis results
become unavailable.
To carry out analysis, on the Analysis ribbon tab, select the Analysis panel and in the Analyse
When analysis is complete, to present stress mosaic plot Rz (soil pressure), on the Results ribbon
tab, on the Stress in plates and solids panel, click Stress mosaic plot Rz . Stress mosaic plot
should be presented for the second DCL.
To compare obtained values of soil pressure Rz with the load on soil Pz, on the WINDOW menu, click
New window and then on the WINDOW menu, click Arrange icons.
In the new window with design model, to present mosaic plot of uniformly distributed load on soil Pz,
on the Advanced results ribbon tab, on the Subgrade moduli panel, click Mosaic plot Pz .
For the next refinement of subgrade moduli, make sure that Stress mosaic plot Rz option is
active for the second DCL and then on the Results ribbon tab, on the Tools panel, click Apply soil
pressure button .
In the Apply soil pressure dialog box (see Fig.9.32), when All elements option is selected, click
OK.
In the new warning box (see Fig.9.33), click OK.
Recalculation of the problem and comparison of obtained values of soil pressure Rz to the load on
soil Pz should be made in a manner identical to the one described above.
The structural engineer decides in each case whether refinement of subgrade moduli is required.
It is recommended that iterations for refinement of subgrade moduli are made up to the state when
difference between values of soil pressure Rz and load on soil Pz does not exceed 5% (but not more
than 5 iterations).
One iteration is taken to mean applying soil pressure and recalculation of the problem with new value
of load on soil.