8.7 C
New York
Monday, February 24, 2025
Home Blog Page 51

Comparative Analysis of Triangular Pile Caps | Free PDF download

When the most economical number of piles that will satisfy design requirements is three, a triangular-shaped pile cap can be adopted. In this case, the piles are arranged in a triangular manner, and equidistant from each other (usually spaced at 3 times the diameter of the pile).

Following the same arrangement, a similarly shaped pile cap can be used to transfer the loads to the piles. Triangular pile caps can be analysed and designed using the strut-and-tie method, beam analogy, or finite element analysis.

In our post on Structural Aspects of Pile Foundation Design, we showed how we can determine the number of piles required to support a column load, and we went ahead to design a triangular pile cap based on the strut-and-tie method (truss analogy).

In this article, we are going to analyse the same triangular pile cap using bending theory and finite element analysis and compare the results obtained.

The layout of the pile cap is shown in Figure 1 and we are going to assume that it is subjected to an axial load of 2500 kN from a 300mm x 300mm column (neglecting the orientation of the column in Figure 1).

Structural layout of column supported by three piles
Figure 1: Layout of triangular pile cap

Bending Theory Analysis of Triangular Pile Cap

The load on each pile is computed based on its distance from the centroid of the column

Load on the single pile
N1 = N(Ay – Ay‘)/Ay
Where:
N = Axial load from column = 2500 kN
Ay = Vertical distance between the centriod of piles in the y-direction = 515 + 1045 = 1560 mm
Ay‘ = Vertical distance between the centroid of the column to the centroid of the single pile in the y-direction = 1045 mm

Hence N1 = 2500 x (1560 – 1045)/1560 = 825. 321 kN

Load on the twin piles
N2 = (N – N1 ) / 2
N2 = (2500 – 825.321)/2 = 837.3395 kN

Bending moment in the x-direction
Mx = N2 x (Ax – Cx)/2
Ax = Spacing of the twin piles = 1800 mm
Cx = Width of the column parallel to x-direction = 300 mm

Mx = 837.3395 x (1.8 – 0.3)/2 = 628 kNm

Bending moment in the y-direction
My = N1 x (Ay‘ – Cy/2)
Ay‘ = Spacing of the twin piles = 1045 mm
Cy = Width of the column parallel to y-direction = 300 mm
My = 825. 321 x (1.045 – 0.3/2) = 738.662 kNm

Finite Element Analysis of Triangular Pile Cap

The pile cap was modelled as a 1300 mm thick concrete plate element on Staad Pro, while the piles were modelled as 600 mm diameter columns with a length of 2m supported on a fixed base. The short column was adopted to avoid wide variation of result due to second-order effects and to avoid large displacement which is not very practicable and also not considered in manual analysis.

3D model of triangular pile cap
Figure 2: 3D rendering of the pile cap model


The results obtained from Staad Pro are given in Figures 3 – 6 below;

Axial load on Single pile
Figure 3: Axial load on Single Pile
Axial load on twin pile
Figure 4: Axial load on twin pile
Bending moment in the x direction
Figure 5: Bending moment in the x-direction
Bending moment in the y direction
Figure 6: Bending moment in the y-direction

The summary of the results between manual analysis (bending theory) and Staad Pro (finite element analysis) is shown in Table 1.

Action effectBending theoryFEA (Staad Pro)Percentage difference
Single pile load (kN)825.32825.770.054%
Twin pile load (kN)837.32837.120.024%
Mx (kNm)628756.0718.507%
My (kNm)738.66762.9043.22%
Table 1: Comparison of bending theory and finite element analysis results


From Table 1, it could be seen that the distribution of load to the piles was found to be approximately the same for both approaches, but a wide variation was observed for Mx bending moment between the two methods. The My bending moment was found to be close with a percentage difference of about 3.229%.

Therefore, Staad Pro can be used to model triangular pile caps and the result can be used for design purposes. The result obtained from this study is open to discussion, especially based on the difference obtained in bending moment values. We will be glad to hear from you!

To download the full structural analysis and design calculation sheet, click HERE

What is Site Instruction?

A site instruction is a formal written order given by the heads of a project to contractors or sub-contractors on specific projects issues such as delay in progress, defective work or materials, guidance on how to carry out a specific item of work, permission to proceed with an item of work, instruction on procurement and logistics, amendments to procedure, or instruction on general project challenges. Site instructions are formal documents that can be presented during disputes, claims, variations, and arbitration. As a result, they should be issued and treated with utmost care and caution by both parties.

Due to its formal nature, a site instruction must be specific, direct, and understandable. A good site instruction must contain the following information;

(1) The name of the firm and the individual issuing the instruction.
(2) The name of the firm and the individual receiving the instruction.
(3) The place, date, and time that the instruction was issued.
(4) The team that issued the instruction and the head of the team (if applicable).
(5) The current observation on site for the item of work for which instruction is being issued
(6) The instruction on what is to be done (with sketches if applicable).
(7) The signature of the person issuing instruction and the signature of the person receiving the instruction

It should be issued in a minimum of 3 duplicates where one copy is kept on site with the site manager, the consultant takes one copy, and another copy is filed in the project’s file at the company’s office. 

Read Also:
Propriety of bending schedules for construction purposes
Seven basic principles of floor and wall tiling

Who issues site instruction?
Site instruction can be issued by the project manager to the consultants or contractors. The consultants can issue instruction to the contractors, and contractors can issue instruction to the sub-contractors. Building regulatory agencies from the government can also issue site instruction to the contractors or consultants of a project. 

For example, if a client decides to make HVAC changes, which will lead to need for new ductworks (say for chiller pipes), the project manager (usually the architects) will have to instruct the HVAC consultants to modify their design, and perhaps instruct the structural engineers to make provision for passage of pipes through the beams. Note that this should be a technically coordinated decision. In as much as the project has started already, such alterations cannot be instructed verbally, hence the need for written site instruction to that effect. The structural engineers may deem the changes not too critical, and may opt to issue site instruction (with sketches) to the main contractors (usually civil works) on how to make adjustments for the pipes to pass through. 

Also, during site inspection, regulatory agencies can issue site instruction to contractors to increase set-backs if found contrary to the approved drawing. They can also order demolition of defective works, or rearrangement of wrongly placed bars etc.

When do you issue/request for site instruction?
(1) A project head should issue site instruction when proposing a change to what is in the approved construction drawing.

(2) If an error or challenge is discovered in the working drawing, the consultant must either issue a new drawing or give a site instruction on how to proceed.

(3) A project head should issue site instruction when proposing a solution to an unforeseen problem or challenge on site which can impact on time, procedure, and cost.

(4) A project head can issue notice of delay in form of a site instruction.

(5) A site manager/contractor should request for an instruction from consultants/regulatory agencies to proceed after an item of work has been satisfactorily completed. For instance, after checking arrangement and placements of reinforcements and formwork, the consultant/regulatory agency should issue an instruction to the contractor to proceed with concreting.

(6) A site manager/contractor should request for a written site instruction when asked to do something he is not comfortable or very familiar with.

(7) A site manager/contractor should request for site instruction to proceed and use a material or procedure he has not tried before.

Site meeting

What do you do after receiving site instruction?
Since site instructions are formal documents, the site manager should therefore take the following steps after receiving instruction:

(1) Contact the company’s head office using the approved means of communication and forward a duplicate copy of the instruction received for filing (an e-mail attachment or posting it to an online workplace platform can suffice).

(2) Discuss and analyse the impact of the instruction on delivery time, procedure, and cost of the project with the project team.

(3) The contractor should send a formal response to the instructor on the impact of the instruction if necessary. It is worthwhile to note that some instructions are very normal in construction and do not need deliberation for an experienced site manager. For instance, a site instruction from a consultant that formwork should not be removed until after 21 days should be seen as normal. However, if this will affect the project delivery time, the contractor can request for an instruction to remove the formwork at say 14 days. The project team can look at the situation technically and a new instruction can be issued to proceed or there can be an outright rebuttal by the consultant. However, it is now on record that there might be 7 days delay depending on the initially agreed programme of work.

Instructions to demolish, remove, or make good a major defective work must be approved from the general office with the input of the quantity surveyor of the project. Such decisions do not require a site manager’s unilateral action. 

(4) A good site manager should work hard and implement all approved corrections according to the instruction.

(5) After all the corrections have been effected, the consultant or agency should be invited to check that the instruction has been carried out properly, after which they will issue a new instruction to proceed with the next item of work.

Do you need good project managers, drawings, and construction advice for your construction works?Send an e-mail to structville@gmail.com or info@structville.com. Alternatively, send a Whatsapp message to +2347053638996

Structural Aspects of Pile Foundation Design: A Practical Example

In the design of a pile foundation, the geotechnical engineer is expected to hand over the soil investigation report to the structural engineer, who will proceed to provide the longitudinal reinforcements needed for the piles and also design the pile cap. The structural design of a pile cap is an important aspect of pile foundation design and the method of carrying out the design has been presented in this article.

The soil investigation report handed to the structural engineer for the design of pile foundation purpose should contain the embedment length of the piles, the recommended sizes of the pile, the safe working load of each pile size, and other information that may be necessary for the structural engineer to do his design properly. All aggressive materials in the soil should be stated as well so that adequate protection will be given to the pile material(s) for durability.

The first step in the structural design of a pile cap usually involves obtaining the number of piles required to support each column load. This is usually done using the service loads of the column, and relating them to the safe working load of the piles from the soil investigation report. In this article, we are going to show how the structural design of reinforced concrete pile foundations and pile caps can be done based on practical design and site experience.

Design Example
The frame of a 5 storey building is shown in Figure 1, and it is expected to be supported on piles with an embedment length of 20 m. The allowable working loads of the bored pile (CFA) is shown in Table 1. fy = 460 Mpa, fcu = 30 Mpa

Table 1: Pile safe working load

Pile diameter (mm)300450600750900
Safe working load (kN)246.74370.11493.48616.85740.22
Building frame modelled on Staad Pro
Fig 1: Frame of a 5-storey building
Column layout
Fig 2: Column load layout

Design for Column A1
Service axial load on column = 647 kN
Ultimate axial load on column = 885 kN
Size of column = 450 x 230 mm

Try 2 No of piles
Service load per pile = 647/2 = 323.5 kN
Let us adopt 600 mm diameter piles for uniformity and to have fewer pile boring points

Safe working load of φ600 mm piles = 493.48 kN > 323.5 kN Okay

Centre to centre spacing of piles = 3φ = 3 x 600 = 1800 mm
Overhang of pile cap edge from the pile = 150 mm

Total length of pile cap = 1800 + 600 + 2(150) = 2700 mm
Width of pile cap = 600 + 150 + 150 = 900 mm
Thickness of pile cap = 2φ + 100 = 2(600) + 100 = 1300 mm

The layout of the pile cap is therefore given as shown in Figure 3.

Structural arrangement of column supported by two piles
Fig 3: Pile cap type 1

Let us quickly carry out the structural design of pile cap Type 1 according to BS 8110-1:1997. You can also read Design of Pile Cap According to Eurocode 2.

From Table 3.61 of Reynolds et al. (2008), the tensile force to be resisted within the pile cap is given by;

Ft = N/(12ld)[3l2 – a2]

Where;
N = Column axial load at ultimate limit state
l = Centre to centre spacing of the piles
d = Effective depth of the pile cap 
a = dimension of the column side parallel to the length of the pile cap

Self weight of the pile cap (ULS) = 1.4 x Area x depth x 24 kN/m3 = 1.4 x 2.7m x 0.9m x 1.3m x 24 kN/m3 = 106.14 kN

N = 885 kN + 106.14 = 991.142 kN
l = 1.8 m
d = 1300 – 100 = 1200 mm = 1.2 m
a = 0.45 m

Ft = [991.142/(12 x 1.8 x 1.2)] x [3 x 1.82 – 0.452] = 364 kN
Ast = Ft/0.95fy = (364 x 1000)/(0.95 x 460) = 833 mm2
Asmin = 0.13bh/100 = 1690 mm2
Provide 6T20 @ 175 c/c (Asprov = 1974 mm2)

Check for shear
Critical position for shear on vertical section across full width of pile-cap occurs at distance from face of column given by:
av = 0.5(l – c) – 0.3φ = 0.5(1800 – 450) – (0.3 x 600) = 495 mm

The shear force carried by the piles V = 991.142/2 = 495.571 kN

The shear stress ν = V/bd = (495.571 x 1000)/(900 x 1200) = 0.458 MPa
Concrete resistance shear stress vc = 0.632(100As/bd)1/3(400/d)1/4

vc = 0.632 x [(100 x 1974)/(900 x 1200)]1/3 x (400/1200)1/4 = 0.632 x 0.557 x 0.759 = 0.275 MPa
For grade 30 concrete, vc = 0.275 x (30/25)1/3 = 0.292 Mpa
vc(2d/av) = 0.292 x [(2 x 1200)/495] = 1.415 MPa > 0.458 MPa  This is okay

Shear stress at column perimeter
ν = V/ud = (885 x 1000)/[(2 x 225 + 2 x 450) x 1200] = 0.546 MPa
This is less than 0.8√fcu = 4.38 Mpa. Therefore, this is okay.

Anti-burst bars of T12 @ 200 spacing should be provided
Main bars should be returned at least 900 mm into the sides to satisfy anchorage length requirements. You can take anchorage length to be conservatively 50 x diameter of reinforcement = 50 x 20 = 1000 mm

Design of Pile Cap Type 2
Service axial load on column = 1077 kN
Ultimate axial load on column = 1476 kN
Size of column = 450 x 225 mm
The number of φ600 mm piles required = 1077/493.48 = 2.184

Using 3 No of φ600 piles
Service load per pile = 1077/3 = 359 kN

Safe working load of φ600 mm piles = 493.48 kN > 359 kN This is okay

Let us adopt a triangular pile cap arranged in such a way that the column load will be equally distributed to the piles. This arrangement can be found in Table 3.16 of Reynolds et al (2008), and it is shown in Figure 4.

Layout of column supported by three piles for equal load distribution
Fig 4: Dimensions of triangular pile cap for equal load distribution (Reynolds et al, 2008)

hp = φ = diameter of pile = 600 mm
Spacing of piles = 3φ = 3 x 600 = 1800 mm
Overhang of pile cap edge from the pile = 150 mm
(α + 1)φ + 300 = (3 + 1)600 + 300 = 2700 mm
φ + 250 =  600 + 250 = 850 mm
φ + 300 =  600 + 300 = 900 mm
(6α/7 + 1)φ + 300 = 2442.857 mm (say = 2445 mm)
(2α/7 + 0.5)φ + 150 = 964.285 mm (say = 965 mm)
Thickness of pile cap = 2φ + 100 = 2(600) + 100 = 1300 mm

The layout of the pile cap is shown in Figure 5.

Structural layout of column supported by three piles
Fig 5: Structural arrangement of 3 pile caps

Self weight of the pile cap (ULS) = 1.4 x Area x depth x 24 kN/m3 = 1.4 x 5.166 m2 x 1.3m x 24 kN/m3 = 225.61 kN

Total load on pile cap at ULS = 1476 kN + 225.61 kN = 1701.61 kN
l = 1.8 m
a = 0.225 m
b = 0.45 m

Tensile force to be resisted by the reinforcement in the direction parallel to X-X;
Ft,x = N/(36ld)[4l2 + b2 – 3a2]
Ft,x = [1701.61/(36 x 1.8 x 1.2)] x [4 x 1.82 + 0.452 – 3 x 0.2252]  = 284 kN

Tensile force to be resisted by the reinforcement in the direction parallel to Y-Y;
Ft,y = N/(18ld)[2l2 – b2]
Ft,y = [1701.61/(18 x 1.8 x 1.2)] x [2 x 1.82 – 0.452]  = 275 kN

Let us use the highest value for design in anticipation that we will provide the same reinforcement in both directions
Ast = Ft/0.95fy = (284 x 1000)/(0.95 x 460) = 649 mm2
Asmin = 0.13bh/100 = 1690 mm2
Provide T20 @ 175 c/c in both directions (Asprov = 1974 mm2)

Shear resistance
The shear force carried by the piles V = 1701.61/3 = 567.2 kN
The shear stress ν = V/bd = (567.2 x 1000)/(1000 x 1200) = 0.472 MPa
vc(2d/av) = 0.292 x [(2 x 1200)/495] = 1.415 MPa > 0.472 MPa  This is okay
Shear will obviously not be a problem.

Design of Pile Cap Type 3
Service axial load on column = 1825 kN
Ultimate axial load on column = 2545 kN
Size of column = 400 x 400 mm
The number of φ600 mm piles required = 1825/493.48 = 3.69

Use 4 No of φ600 piles
Service load per pile = 1825/4 = 456.25 kN

Safe working load of φ600 mm piles = 493.48 kN > 456.25 kN This is okay

Let us adopt a square pile cap arranged in such a way that the column load will be equally distributed to the piles. This arrangement can be found in Table 3.16 of Reynolds et al (2008), and it is shown in Figure 6.

Structural layout of column supported by four piles
Fig 6: Pile Cap Type 3

Self weight of the pile cap (ULS) = 1.4 x Area x depth x 24 kN/m3 = 1.4 x 7.29 m2 x 1.3m x 24 kN/m3 = 318.43 kN
Total load on pile cap at ULS = 2545 kN + 318.43 kN = 2863.43 kN

Tensile force to be resisted by the reinforcement in both directions;
Ft = N/(24ld)[3l2 – a2]
Ft = [2863.43/(24 x 1.8 x 1.2)] x [3 x 1.82 – 0.402]  = 528 kN

Ast = Ft/0.95fy = (528 x 1000)/(0.95 x 460) = 1208 mm2
Asmin = 0.13bh/100 = 1690 mm2
Provide T20 @ 175 c/c in both directions (Asprov = 1974 mm2)

Shear resistance
The shear force carried by the piles V = 2863.43/4 = 715.9 kN
The shear stress ν = V/bd = (715.9 x 1000)/(1000 x 1200) = 0.595 MPa
vc(2d/av) = 0.292 x [(2 x 1200)/495] = 1.415 MPa > 0.472 MPa  This is okay

The structural engineer is expected to provide the following drawings;

(1) Setting out drawing showing the piling points and layout with a known reference point
(2) General column/pile cap layout/arrangement
(3) Pile cap/ground beam/ground floor slab layout
(3) Columns, piles and pile cap reinforcement drawings (detailing)
(5) Ground beams and ground floor slab reinforcement detailing
(6) Construction procedure sketches

Construction Considerations
(1) Ground beams are usually provided to chain the pile caps together, and to provide the needed support for the ground floor slab. There are construction scenarios where the ground floor slab is placed directly on the pile caps, but note that this concept is quite different from piled raft foundation. The ground beams are usually embedded into the pile caps, or may be allowed to sit directly on the pile caps depending on the site level. A typical construction drawing showing this interaction is given in Figure 7.

(2) Contractor should maintain a minimum concrete cover of 75 mm

(3) There may be need to be cast the pile cap in two stages to achieve the configuration shown in Figure 7.  The first casting will get to the bottom level of the ground beams (see Figure 8), then the ground beam reinforcements are laid (see Figure 9), before the final casting of the pile cap and ground beams to the required level (see Figure 10). Read about bonding of old and new concrete.

Pile cap and ground beam arrangement
Concreting of pile cap
Figure 8: Typical casting of pile cap to ground beam level.
Ground beam reinforcement
Fig 9: Typical ground beam reinforcement arrangement
Reinforced concrete beam and pile cap
Fig 10: Completed pile cap and ground beam

On getting to the stage shown in Figure 10, the bays are filled with sharp sand, and the ground floor slab is cast as appropriate.

Do you need help in design, consultancy, production of construction drawings, supervision, and project management, contact us today at Structville Integrated Services Limited. We are excellent at what we do, and we pride ourselves in professionalism and integrity. Send an e-mail now to info@structville.com copy structville@gmail.com or send a whatsapp message to +2347053638996.

References
[1] Reynolds C.E., Steedman J.C., Threlfall A.J. (2008): Reynolds’ Reinforced Concrete Designer’s Handbook 11th Edition. Taylor and Francis, New York

Calculation of Collapse Settlement of Collapsible Soils

Civil engineers are always faced with one challenge or another on-site, and one of the major challenges on-site has to do with the nature of the soil to support the structure to be constructed. Some soils exhibit excessive change in volume under constant load and changes in moisture content. Under field conditions, if the addition of water under constant load causes the void ratio to reduce drastically, the soil is said to have collapsed.

These soils are predominantly found in the arid regions and appear strong and stiff in their dry natural state, but lose strength and undergo high compression upon wetting.  According to Behzad (2013), almost all natural deposits of collapsible soils are either debris flow deposits or wind deposited soils (loess).

Soils which exhibit this behaviour are usually unsaturated granular soils which in a loose state are maintained by apparent cohesion. This apparent cohesion may be due to the presence of clays at the intergranular contact areas or the accumulation of soluble salts as binders (Murthy, 2012). In South-Western Nigeria, collapsible soils have been observed and studied by Owolabi and Ola (2014).

Collapsible Soils

For a soil sample to collapse, the soil structure must have an open, partially unstable, and unsaturated fabric which must be held in place by apparent cohesion. On wetting and application of high enough stress, the soil fabric will be expected to collapse.

For compacted cohesive soils, the collapse behaviour will depend on the percentage of fines, the moulding water content, the dry density, and the compaction energy employed Murthy (2012). Any soil compacted at the dry side of the optimum moisture content and at low dry density will also be expected to develop a collapsible structure.

Calculation of Collapse Settlement

A method of calculating the collapse settlement of a soil sample was suggested by Jennings and Knight (1975) and reported by Murthy (2012). This involves carrying out double oedometer tests of identical undisturbed samples in the laboratory. Both specimens are kept under a constant pressure of 1 kPa for a period of 24 hours at their natural moisture content after which one specimen is flooded and the other allowed to remain in its natural state.

After flooding for 24 hours, the consolidation test is allowed to progress as usual after which the e-LogP curve is plotted for both samples. The collapse settlement is calculated in two parts according to equations (1) and (2).

S1 = ΔenH /(1 + e0)  ————————- (1)
S2 = ΔecH /(1 + e0)  ————————- (2)

Where;
Δen =  Change in void ratio due to load ΔP as per the e-logP plot without change in moisture content
Δec =  Change in void ratio due to load ΔP as per the e-logP plot with increase in moisture content (settlement due to the collapse of the soil structure).
H = Thickness of the collapsible layer
e0 = Initial void ratio

Collapse settlement (S1) = S1 + S2

Solved Example

A building was constructed on collapsible soil and the double oedometer test carried out on the undisturbed sample gave the results shown in Table 1. The thickness of the collapsible layer is 4 m and the average overburden pressure was 60 kPa. Calculate the collapse settlement for an increase in pressure of 40 kPa.

Applied Pressure (kPa)102040100200400800
Void ratio at natural moisture content0.800.790.780.750.7250.680.61
Void ratio at soaked condition0.750.710.660.580.510.430.32
Table 1-


Solution

(1) The first step is to plot the e-log P curve for the two samples as shown below. This has been done using Microsoft Excel. 

e-log P plot for double oedometer test
Figure 1: e-log P plot for double oedometer test

(2) The virgin compression curve is plotted by drawing a line tangential to the soaked sample curve.

(3) After drawing the virgin compression line, the overburden pressure is identified, and for this problem, it is given at 60 kPa. This is traced up to meet the virgin compression line, and this marks the initial void ratio (e0).

(4) From point e0, an adjusted moisture curve is plotted to be perfectly parrallel to the natural moisture content curve.

(5) The increment pressure p + ΔP = 60 + 40 = 100 kPa. A vertical line is drawn from the pressure of 100 kPa to meet the soaked sample curve and adjusted moisture curve.

Steps (2) to (5) are shown in Figure 2.

Graph for calculation of collapse settlement
Fig 2: Full analytical plot for collapse settlement of the soil

(6) Trace down the points where the pressure increment makes contact with the soaked sample curve and adjusted moisture content curve and identify the associated void ratios e1 and e2.

For the problem at hand;
e0 = 0.64
e1 = 0.62
e2 = 0.58

Therefore;
Δen = e0 – e1 = 0.64 – 0.62 = 0.02
Δec = e1 – e2 = 0.62 – 0.58 = 0.04

S1 = ΔenH /(1 + e0) = (0.02 x 4)/(1 + 0.64) = 0.0487 m
S2 = ΔecH /(1 + e0) = (0.04 x 4)/(1 + 0.64) = 0.0975 m

The total collapse settlement is therefore given by S1 + S2 = 0.0487 + 0.0975 = 0.146 m = 146.2 mm

Methods of Foundation and Construction on Collapsible Soils

(1) Utilising deep foundations such as piles to send the load beyond the collapsible layer.

(2) Chaining the foundation using continuous reinforced ground beams and slab to minimise differential settlement.

(3) In a laboratory study carried out by AlShaba et al.. (2018) at Egypt, iron powder was used to reduce collapse settlement of collapsible soil by about 46% at a mix ratio of 6% of the dry weight of the soil. When the same mixture was compacted, a reduction ratio of about 86.66% was obtained.

(4) Prewetting and compaction of the collapsible soil.

(5) Ayeldeen et al (2016) obtained promising results on improving collapsible soils using biopolymers.

References
[1] Ayeldeen M., Nrgm A., El-Sawwaf M., Kitazume M. (2016): Enhancing the Behaviour of Collapsible Soils using Bioploymers. Journal of Rock Mechanics and Geotechnical Engineering. doi:10.1016/j/jrmge.2016.11.007
[2] AlShaba A.A., Abdelaziz T.M., Ragheb A.M. (2018): Treatment of Collapsible Soils by Mixing with Iron Powder. Elsevier – Alexandria Engineering Journal (57):3737-3745
[3] Behzad Kalantari (2013): Foundations of Collapsible Soils: A Review. Proceedings of the Institute of Civil Engineers Forensic Engineering 166(FE2):57-63
[4] Jennings J.E., and Knight K. (1975): An Additional Settlement of Foundation Due to Collapse Structure of Snady Soils on Wetting. Proceedings to thhe 4th International Conference on Soil Mechanics and Foundation Engineering, Vol 1. Paris
[5] Murthy V.N.S (2012): Textbook of Soil Mechanics and Foundation Engineering (First Edition). CBS Publishers and Distributors Pvt Ltd, India
[6] Owolabi F.A., Ola S.A. (2014): Geotechnical Properties of a Typical Collapsible Soil i South-Western Nigeria. Electronic Journal of Geotechnical Engineering (19):1721-1738



Structural Analysis of Portal Frames Subjected to Gravity Load

Gravity actions are normally used for the verification of portal frames for member resistances, lateral buckling, and torsional stability. Elastic design of portal frames is permitted by Eurocode 3, and the most significant load case (from experience) is the situation where the frame is subjected to gravity load from permanent actions and variable actions, taking into account the second-order effects and imperfection.

The load combination used for this is usually of the form shown in Equation (1);

p = 1.35gk + 1.5qk (1)


In the past, we have solved the problem of two hinged portal frames which are statically determinate due to the presence of internal hinge at the apex. However, it is not a very practical scenario as portal frames are required to be rigidly connected in order to achieve in-plane stability.

In this article, we are going to show how you can obtain the effects of actions (bending moment, shear force, and axial force) for a rigid portal frame subjected to an ultimate limit state load of 12 kN/m as shown in Figure 1. We will assume that the same universal beam section will be used for the column and the rafters (EI = Constant).

Portal Frame
Fig 1: Portal frame

A little consideration will show that the frame shown in Figure 1 is statically indeterminate to the first order. Therefore, using the force method of analysis, we will need to reduce the structure to a basic system.

This is a process by which we remove the redundants (excess reactive forces) in the structure, to make it statically determinate. The basic system adopted must also be stable. This is shown in Figure 2.

Basic System 1
Fig 2: Reduction of the frame to basic system

To proceed, we will now assign a unit force to the removed redundant horizontal reaction, and plot the bending moment diagram due to the unit force on the basic system as shown in Figure 3.

Bending Moment Diagram
Fig 3: Bending moment diagram due to horizontal unit force at point A

Since there is no other redundant support, we can apply the external load on the basic system. Due to the symmetry of the structure and the loading, we can verify that the vertical support reactions at A and E are given as shown in Figure 4;

Ay = Ey = (12 x 18)/2 = 108 kN

LOAD ON BASIC SYSTEM
Fig 4: External load and support reaction on the basic system

On observing Figure 4, we can see that there will be no bending moment in sections A-B and D-E, rather they will be subjected to compressive axial force of 108 kN. For member B-C, the bending moment will induced due to the reactive force and the externally applied uniformly distributed load. The free-body diagram of member B-C is shown in Figure 5.

We can verify that the geometrical properties of B-C are as follows;
The angle of inclination β = tan-1(1.5/9) = 9.462º
The length of the member z = √(1.52 + 92) = 9.124 m

section bc
Fig 5: Free body diagram of member of B-C

The equation for the bending moment of the section can be given by;
Mz = [108 x cos(β) x z] – [12 x cos(β) x z2]/2
Mz = 106.53z – 5.918z2

At z = 9.124 m
MCL = 106.53(9.124) – 5.918(9.124)2 = 479.32 kNm

The same thing is also happening at section C-D, and the bending moment diagram of the structure is shown in Figure 6.

BMD
Figure 6: Bending moment diagram of external load on basic system

The canonical equation for the structure is given by Equation (2);

δ11X1 + Δ1P = 0 ——– (2)

Where;
δ11 = Deflection at point 1 due to unit force at point 1
X1 = The actual external reaction at point 1 in the direction of the deflection due to externally applied load
Δ1P = Deflection at point 1 due to externally applied load

We are going to evaluate δ11 and Δ1P according to Vereshchagin’s rule which involves combining of the bending moment diagrams. For another example of how to apply Vereshchagin’s rule, click here.

Evaluation of δ11
In this case, the bending moment diagram of the unit force will combine with itself.

Combiner 1
Fig 7: Load Case combining with itself

δ11 = 2 x [1/3 x 8 x 8 x 8] + 2 x 1/3 x [(8 x 8) + (9.5 x 8) + (9.5 x 9.5)] x 9.124 = 1741.867/EI

Evaluation of Δ1P
Here, the bending moment diagram due to the unit force will combine with the bending moment diagram due to the externally applied load as shown in Figure 8.

Combiner 2
Fig 8: Load case 1 combined with externally applied load

Δ1P = 2 x (1/12) x 479.32 x [(5 x 9.5) + (3 x 8)] x 9.124 = 52115.345

Substituting back into the canonical equation;
X1 = 52115.345/1741.867 = 29.919 kN (This obviously represents the horizontal reactions Ax and Ex)

The final bending moment diagram can be obtained therefore as follows;
Mdef = M1X1 + M0
MA = ME = 0
MB = MD = (29.919 x -8) + 0 = -239.352 kNm
MC = (29.919 x -9.5) + 479.32 = 195.089 kNm

The bending moment diagram is shown in Figure 9, and the shear and axial forces can also be obtained accordingly.

Bending moment diagram of portal frame
Fig 9: Final bending moment diagram

Thank you for visiting Structville today… Please check back for more updates.

On the Consolidation Settlement of Pile Groups

When a group of piles are embedded in saturated clays, the pile group undergoes time-dependent consolidation settlement. This can be calculated using the consolidation settlement equations and generally involves evaluating the increase in stress when the soil beneath a pile group is subjected to pressure due to load from the superstructure/column Qg.

The consolidation settlement of pile groups can be estimated by assuming an approximate distribution method that is commonly referred to as the 2:1 method.

The calculation procedure involves the following steps as described in Das (2008);

1. Obtain the total length of embedment of the pile L, and the total load (service load) transmitted from the column to the pile cap (Qg).
2. The load Qg is assumed to act on a fictitious footing depth of 2L/3 from the top of the pile. The load Qg spreads out along 2 vertical:1 horizontal lines from this depth.
3. Calculate the effective stress increase caused at the middle of each soil layer by the load Qg:

∆σi‘ = Qg/[(Bg + zi) (Lg + zi)]

Where;
∆σi‘ = effective stress increase at the middle of layer i
Bg, Lg = Width and length of the plan of pile group, respectively
zi = distance from z = 0 to the middle of the clay layer, i

4. Calculate the settlement of each layer caused by the increased stress:

∆Sc(i) = [∆e/ (1 + e0(i)] × Hi

where;
∆Sc(i) = consolidation settlement of layer i
∆ei = change of void ratio caused by the stress increase in layer i
e0(i) = initial void ratio of layer i (before construction)
Hi = Thickness of layer i

5. Calculate the total consolidation settlement of the pile group by;
∆Sc(g) = ∑∆Sc(i)

Solved Example

Calculate the consolidation settlement of a group of 4 piles founded on layers of clay as shown below. The piles are 600mm in diameter and form a square plan of 1800mm x 1800mm. The service load from the column on the pile is 1350 kN.

consolidation settlement of pile group

Solution
Depth of fictitious footing = 2L/3 = 2(20)/3 = 13.33 m
By implication, it can be seen that the influence depth went beyond the first layer of clay. Therefore, consolidation settlement calculation commences from the second clay layer of the soil profile. Let us therefore call this first layer of influence, layer 1.

Layer 1
∆σ1‘ = Qg/[(Bg + z1) (Lg + z1)] = 1350 / [(1.8 + 4.585) × (1.8 + 4.585)] = 33.114 kPa

σ0(1)‘ = (2 × 17) + 12(18.5 – 9.81) + 5.415(19 – 9.81) = 188.04 kPa

∆Sc(1) = [(Cc(1)H1)/(1 + e0)] × log[(σ0(i)‘ + ∆σi‘ )/σ0(i)‘]
∆Sc(1) = [(0.28 × 9.17)/(1 + 0.73)] × log[(188.04 + 33.114)/188.04] = 0.1045m = 104.473 mm

Layer 2
∆σ2‘ = Qg/[(Bg + z2) (Lg + z2)] = 1350 / [(1.8 + 12.67) × (1.8 + 12.67)] = 6.4475 kPa

σ0(2)‘ = (2 × 17) + 12(18.5 – 9.81) + 10(19 – 9.81) + 3.5(18.3 – 9.81) = 259.895 kPa

∆Sc(2) = [(Cc(2)H2)/(1 + e0)] × log[(σ0(2)‘ + ∆σ2‘ )/σ0(2)‘]
∆Sc(2) = [(0.25 × 7)/(1 + 0.77)] × log[(259.895 + 6.4475)/259.895] = 0.0104m = 10.48 mm

Therefore total settlement ∆Sc(g) = ∆Sc(1) + ∆Sc(2) = 104.473 + 10.48 = 114.953 mm

Thank you so much for visiting Structville today, and God bless you.

References
Das B.M. (2008): Fundamentals of Geotechnical Engineering (3rd Edition). Cengage Learning, USA

Analysis of Beams on Elastic Foundation

When a beam lies on an elastic foundation under the action of externally applied loads, the reaction forces of the foundation are proportional at every point to the deflection of the beam. This assumption was introduced first by Winkler in 1867. There are cases in which beams are supported on foundations which develop essentially continuous reactions that are proportional at each position along the beam to the deflection of the beam at that position.

This is the reason for the name ‘elastic foundation’. There are many geotechnical engineering problems that can be idealized as beams on elastic foundations. This kind of modelling helps to understand the soil-structure interaction phenomenon and predict the contact pressure distribution and deformation within the medium (e.g. soil). The most common theory for a beam on elastic foundation modelling is the Winkler approach.

image 5
Figure 1: Schematic representation of beam on elastic foundation


For instance, the beam shown in Fig. 1 will deflect due to the externally applied load, and produce continuously distributed reaction forces in the supporting medium. The intensity of these reaction forces at any point is proportional to the deflection of the beam y(x) at that point via the constant ks:

R(x) = ks.y(x)  — Eq. (1)

Where kis the soil’s modulus of subgrade reaction which is the pressure per unit settlement of the soil (unit in kN/m2/m). The reactions R(x) act vertically and oppose the deflection of the beam. Hence, where the deflection is acting downward there will be a compression in the supporting medium. Where the deflection happens to be upward in the supporting medium tension will be produced which is not possible (for soils).

If we assume that the beam under consideration has a constant cross-section with constant width b which is supported by the foundation. A unit deflection of this beam will cause a reaction equal to ks.b in the foundation, therefore the intensity of distributed reaction (per unit length of the beam) will be:

R(x) = b.ks.y(x) = k.y(x) —-Eq. (2)

where k = k0.b is the constant of the foundation, known as Winkler’s constant, which includes the effect of the width of the beam, and has unit of kN/m/m.

The general 4th order differential equation for beam on elastic foundation is given by equation (3);

EI(d4y)/dx4 + k.y = q  —- Eq. (3)

The homogenous equation is given by;

EI(d4y)/dx4 + k.y = 0

(d4y)/dx4 + 4β4y = 0

Where β = ∜(k/4EI) = (k/EI)(1/4)

The general solution for the equation is available, which is given by;
y = eβx (C1sinβx + C2cosβx) + e-βx(C3sinβx + C4cosβx)

Warren and Richard (2002) published tables containing equations for analysing beams on elastic foundation subjected to different loads. The method described in the book has been employed in this article to analyse a beam on an elastic foundation and the results compared with the results from Staad Pro software.

Solved Example
A 600mm x 400mm rectangular beam is resting on a homogenous soil of modulus of subgrade reaction ks = 10000 kN/m2/m. The beam is 10m long and carries a concentrated load of W = 300 kN at 3m from the left-hand side. Neglecting the self-weight of the beam, and assuming freely supported ends, obtain the bending moment at the point of the concentrated load. Take the modulus of elasticity of concrete Ec = 21.7 x 106 kN/m2.

BEAM 2BON 2BELASTIC 2BFOUNDATION
Example of beam resting on an elastic foundation
beam 2Bsection

Solution 
Second moment of area of concrete beam IB = (bh3)/12 = (0.4 × 0.63)/12 = 7.2 x 10-3 m4
Flexural rigidity of the beam EcIB = 21.7 × 106 × 7.2 × 10-3 = 156240 kN.m2
β = (bks/4EI)(1/4) = [(0.4 × 10000)/(4 × 156240)](1/4) = 0.2828
βl = 0.2828 × 10 = 2.828 m; β(l – a) = 0.282(10 – 3) = 1.979

Where a is the distance of the concentrated load from the left end of the beam.
Since βl < 6.0, we can use Table 8.5 of Roark’s Table for Stress and Strain (Warren and Richard, 2002).

For a beam with both ends free;
RA = 0; MA = 0

The equations for bending moment and shear force along the beam is as given below;

Mx = MAF1 + RA/2βF2 – yA2EIβ2F3 – θAEIβF4 – W/2βFa2
Vx = RAF1 – yA2EIβ3F2 – θAEIβ2F3 – MAβF4WFa1

Where;
θA (rotation at point A) = [W/(2EIβ2)] × [(C2Ca2 – 2C3Ca1)/C11]
yA (vertical deflection at point A) = [W/(2EIβ3)] × [(C4Ca1 – C3Ca2)/C11]

We can therefore compute the constants as follows;
C2 = coshβl.sinβl + sinhβl.cosβl = [cosh(2.828) × sin(2.828)] + [sinh(2.828) × cos(2.828)] = (8.485 × 0.3084) + (8.426 × –0.9512 ) = –5.398

C3 = sinhβl.sinβl = (8.426 × 0.3084) = 2.5985

C4 = coshβl.sinβl – sinhβl.cosβl = [cosh(2.828) × sin(2.828)] + [sinh(2.828) × cos(2.828)] = (8.485 × 0.3084) – (8.426 × –0.9512) = 10.631

Ca1 = coshβ(l – a).cosβ(l – a) = cosh(1.974) × cos(1.974) = 3.669 × -0.392 = -1.438

Ca2 = coshβ(l – a).sinβ(l – a) + sinhβ(l – a).cosβ(l – a) = [cosh(1.974) × sin(1.974)] + [sinh(1.974) × cos(1.974)] = (3.669 × 0.9198) + (3.530 × -0.392) = 1.9909

Ca3 = sinhβ(l – a).sinβ(l – a) = (3.530 × 0.9198) = 3.247

C11 = sinh2βl – sin2βl = 8.4262 – 0.30842 = 70.902

θA = W/(2EIβ2) × [(C2Ca2 – 2C3Ca1)/C11]
θA = [300/(2 × 156240 × 0.28282 ) × [(–5.398 × 1.9909) – (2 × 2.5985 × –1.438)/70.902] = (0.012 × –0.0462) = -0.0005544 radians

yA (vertical deformation at point A) = [W/(2EIβ3)] × [(C4Ca1 – C3Ca2)/C11]
y = [300/(2 × 156240 × 0.28283)] × [(10.631 × –1.438) –2.5985 × 1.9909)/70.902] = (0.0424 × –0.2885) = –0.01223 m = -12.23 mm

Bending moment at point C
Substituting the values of deflection and slope into the equation for bending moment (note that the first and second terms of the equation goes to zero since RA = MA = 0);

Mx = – yAEIβ2F3  – θAEIβF4 – W/2βFa2
Mx = – (–0.01223 × 2 × 156240 × 0.28282)F3 – (– 0.0005544 × 156240 × 0.2828)F4 – 300/(2 × 0.2828)Fa2

Mx = 305.638F3 + 24.495F4 – 530.410Fa2

Substituting F3, F4, and Fa2 (see Table 8.5, Warren and Richard, 2002) into the equation;

Mx = 305.638(sinhβx.sinβx) + 24.495(coshβx.sinβx – sinhβx.cosβx) – 530.410[coshβ(x –a).sin β(x –a) + sinh β(x –a).cos β(x –a)]

The bending moment under the concentrated load (point C);
x = 3m; (x – a) = 3 – 3 = 0

βx = 0.8484
Mx = 305.638[sinh(0.8484) × sin(0.8484)] + 24.495[cosh(0.8484) × sin(0.8484) – sinh(0.8484) × cos(0.8484)]

Mc = 305.638 (0.953 × 0.750) + 24.495(1.382 × 0.750 – 0.9539 × 0.6611) = 218.45 + 9.942 = 228.392 kN.m


Verification
This manual calculation has been verified using Staad Pro software. The steps adopted were as follows:

(1) Modelling
The 10m beam was modelled as a one-dimensional line element connected by nodes at 1m length interval. This was to represent/attach soil springs at 1m interval.

node 2Bconnection

(2) Section Properties

The beam was modelled using concrete of modulus of elasticity = 21.7 x 106 kN/m2, with dimensions of 600mm x 400mm.

(3) Support 

Using the elastic mat foundation option will not work since the support will not form a closed loop (analysing this way will give a ‘colinear support error’). Therefore, the soil spring was modelled using the ‘FIXED BUT’ support option. The support must be released for moment since this is the ideal scenario for the structure we are trying to model.

The soil modulus of the subgrade reaction was multiplied by the width of the beam thus;

kb = 10000 kN/m2/m x 0.4m = 4000 kN/m2

Now, given that the nodes are spaced at 1m interval, the vertical spring constant was taken as 4000 kN/m. The general form of the foundation is given below;

soil 2Bspring 2Bmodel
Modelling of beam on elastic foundation

(4) Analysis and Results

When analysed using the static check option, the following results were obtained;

(a) Soil deformation

settlement

The vertical deflection at point A (node 1) was observed to be 12.195mm, against 12.23mm obtained using manual analysis.

(b) Bending Moment

Bending 2Bmoment
Staad Pro result for beam on elastic foundation

The bending moment obtained at point C was 228.146 kNm which is comparable with 228.392 kNm obtained using manual calculations. Therefore, it can be seen that both methods give approximately the same result.

References
Warren C.Y., Richard B. Y. (2002): Roark’s Formula for Stress and Strain (7th edition). McGraw Hill, USA


Analysis and Design of Cantilever Retaining Walls on Staad Pro

In this post, we are going to show how cantilever retaining walls can be analysed and designed on Staad Pro software, and also compare the answer obtained with classical solutions. We should know that retaining walls must satisfy geotechnical, equilibrium, structural, upheaval, seismic considerations, etc. As a result, the designer must ensure that by appropriate knowledge of materials, site conditions, etc, he/she will provide suitable dimensions of the retaining wall that will ensure resistance of the structure to overturning, sliding, bearing capacity failure, uplift, etc. After appropriate sizing of the retaining wall, the structural analysis and design will commence to determine the action effects (bending moments, shear forces, axial forces, deflection etc), and provision of proper reinforcements to resist the action effects.


In the past, Structville has published a 17 page document on geotechnical design of cantilever retaining walls subjected to earth load, pavement surcharge load, traffic load, etc. This loading situation can be found when retaining wall is used to support embankment carrying traffic road way.  It was interesting to see how Design Approach 1 (DA1) of Eurocode 7 was used to ensure the geotechnical stability of the wall. Just in case you missed it, kindly download the PDF from the link below;

Geotechnical Design of Cantilever Retaining Walls to Eurocode 7

In this post, let us consider the retaining wall sized and loaded as shown in Figure 2. This structure has been modelled on Staad Pro in order to determine the action effects due to the applied load.

Cantilever 2Bretaining 2Bwall
Fig 2: Cantilever retaining wall

The retaining wall is subjected to a 3m thick earthfill, and a variable surcharge pressure of 10 kPa. Given that the retained earth has an angle of internal friction of 30°, we can obtain the Rankine active earth pressure as follows;

ka = (1 – sin 30)/(1 + sin 30) = 0.333

front
[Get this publication on Design of Swimming Pools and Underground Water Tanks into your mailbox by clicking HERE or by clicking on the book cover]

Therefore the actions on the retaining wall that will be input into Staad Pro as are as follows;

Vertical actions
(1) Self weight (to be calculated automatically by Staad)
(2) Weight of earthfill (19 kN/m3 × 3m) = 57 kN/m2
(3) Surcharge load = 10 kN/m2

Horizontal Actions
(4) Triangular earth pressure = (0.333 × 19 × 3) = 18.98 kN/m2
(5) Uniform surcharge pressure = (0.333 × 10) = 3.333 kN/m2

The wall has been modelled per metre run on Staad, and plate mat foundation was utilised with coefficient of subgrade modulus of 100000 kN/m2/m.

Steps to adopt
(1) Model the retaining wall utilising plate element meshing, assign thickness of 0.4m to the base, and 0.3m to the wall. Also assign plate mat foundation of subgrade modulus 100000 kN/m2/m to the base in the y-direction.

RET 2B1
Fig 3: Modelling and meshing of the retaining wall

(2) Assign the following loads to the structure

Load Case 1 (LC1)
(a) Self weight to the whole structure
(b) Weight of earth fill to the heel of the retaining wall (57 kN/m2)
(c) Assume that the base is buried 1m into the ground, hence apply vertical pressure load of (19 kN/m2) to the toe but neglect all passive pressures.

a1
Fig 4: Permanent vertical actions on the retaining wall base
final 2Bfront 2Bcover
[Get this publication on Design of Residential Buildings Using Staad, Orion, and Manual Calculations into your mailbox by clicking HERE or by clicking on the book cover]

Load Case 2 (LC2)
(d) Triangular earth pressure to the wall (18.98 kN/m2)

a2
Fig 5: Horizontal active earth pressure on the wall

Load Case 3 (LC3)
(e) Uniform horizontal surcharge pressure to the walls (3.333 kN/m2)
(f) Uniform vertical surcharge pressure to the heel (10 kN/m2)

a3
Fig 6: Surcharge loads on the wall and on the base


Combination (Ultimate limit state)
pEd = 1.35LC1 + 1.35LC2 + 1.5LC3

A little consideration will show that the load cases 1 and 2 are treated as permanent actions, while load case 3 is treated as a variable action.

(3) Analyse the structure for the load cases
(N/B): You may need to increase the iteration limits for the load cases containing horizontal actions to converge

MX1
Fig 7: Main bending moment on the retaining wall
MXY
Fig 8: Twisting moment on the retaining wall
SQX
Fig 9: Shear stress on the retaining wall

On considering the bending moment diagram Mx (Figure 7), we can see that the maximum moment close to the base of the wall is 57 kNm/m. If we add the effect of torsion (Figure 8), the design moment can be taken as 57 + 2.99 = 59.99 kNm/m.

To carry out manual analysis, we will have to follow the steps given below to obtain the maximum moment at the base of the wall. The actions causing bending on the wall are the horizontal earth pressure and the horizontal surcharge pressure.

Moment from surcharge pressure = [(3.333 × 32)/2] = 14.998 kNm/m
Moment from horizontal earth pressure = [(18.98 × 3)/2] × (3/3) = 28.47 kNm/m

At ultimate limit state, MEd = 1.35(28.47) + 1.5(14.998) = 60.931 kNm/m

A little consideration will show that Staad Pro and Manual calculations gave almost the same value for wall bending moment. However, I expect the value of base moment from Staad Pro to be lower than the one from manual analysis. Kindly verify this at your private time.


Thank you for visiting Structville today and God bless you.

An Investigation on the Analysis of Beam and Raft Slab Using Staad Pro

Modern codes of practice are increasingly recognising the computational power of structural analysis  softwares. Staad Pro is a renowned structural analysis and design software that is used all over the world, and this post is aimed at investigating the analysis of beam and raft foundation using Staad Pro V8i. Winkler’s model has been used as a basis for the analysis. This post will help engineers who use Staad Pro to make decisions on whether raft slab can be confidently analysed on Staad Pro V8i environment, or there will be need to move to Staad Foundation, or any other foundation analysis software.

Winklers 2BModel
Fig 1: Winkler’s Model

Winkler’s model assumes that soil possesses stiffness which is considered to be the ratio of the contact pressure of the soil, and the vertical deformation associated with it. This relationship is assumed to be linear, and can be given by what is called the coefficient of subgrade reaction (ks).

ks = q/S ——————- (1)

where;
q is the contact pressure at a point in the footing, and
S is the settlement at the same point.

Equation (1) assumes that for granular soils, the value of ks is independent of the magnitude of pressure, and is the same at all points on the surface of the footing. These two assumptions are however not very accurate for obvious reasons, especially for flexible foundations. There are lots of literature available on the topic of ‘coefficient of subgrade reaction’, and you are advised to read them up on Google. However, it is accepted that the method gives realistic contact pressures especially when considering very low order of settlement. Coefficient of subgrade reaction is obtained in the field using in-situ plate load test. Plate load tests are usually carried out using (300 x 300) mm plate, and there is usually need to correct the obtained value for the actual width of the foundation (for a quick check of procedure, see page 696 of Murthy, 2012). I will not actually dabble into the accuracy of application of coefficient of subgrade reaction for foundation design in this post, but the reader is advised to consult wide range of literature since opinion on this seems to vary.

front
[Get this publication on Design of Swimming Pools and Underground Water Tanks into your mailbox by clicking HERE or by clicking on the book cover]

In Staad Pro software, the use of Winkler’s Model can easily be done by applying ‘foundation’ support option to the structure in question. Under this option, you can use ‘elastic mat’ or ‘plate mat’. Elastic mat is suitable for beam and plate combination (e.g beam and raft slab in a building) while plate mat is good when you have plate elements only forming the base (e.g base of an underground water tank). For training on how to use Staad Pro efficiently, contact the author on ubani@structville.com


In many parts of Nigeria, beam and raft slab is employed for low-medium rise buildings when the soil has low bearing capacity (usually between 100 kPa and 50kPa) or where providing separate bases will pose so many challenges. For manual analysis of such foundations, the rigid approach is usually used, which yields higher value of internal stresses, while on the other hand, finite element based foundation design softwares can be used, which usually yields lower internal stresses but higher settlement (flexible approach).

In the design of beam and raft slab foundation, it is usually assumed that the column load gets transferred/distributed to the ground beams, and tries to push them down into the soil. This action is resisted by the earth pressure intensity, which is mainly transferred to the ground floor slab. The schematic diagram of this action is shown below.

Schematic 2BDiagram 2Bof 2BBeam 2Band 2BRaft 2BFoundation
Fig 2: Load Path of Beam and Raft Foundation

Design Example
As an example, let us consider a two storey building with the foundation layout shown below.

GA
Fig 3: Typical Foundation Layout of Structure

For the beam and raft layout shown in Fig. 3, the necessary data are as follows;

Ground beam dimensions = 1200 x 250 mm
Thickness of ground floor slab = 150 mm
Dimension of all columns = 230 x 230 mm
First floor slab = 150 mm
First floor beams = 450 x 230 mm
Roof beams = 300 x 230 mm
Additional dead load on floors (finishes and partition allowance) = 2 kN/m2
Imposed load on building (variable action) = 2 kN/m2
Load at ultimate limit state = 1.35gk + 1.5qk
Coefficient of subgrade reaction ks =  of 20000 kN/m2/m (assumed for very soft clay, reference taken from CSC Orion Software)

The 3D view of the frame is shown in Figure 4 below;

3D 2BFRAME 2BOF 2BBUILDING
Fig 4: 3D Model of the Building

After analysis of the superstructure on Staad, the following ultimate column axial loads which are transferred to the foundation are shown below;

COLUMN 2BLOAD 2BON 2BRAFT 2BSLAB
Fig 5: Column Loads on the Foundation

From figure 5 above, I believe you can easily model the foundation in another software, or carry out manual analysis for checks of the procedure adopted here.

When the ground beams and elastic mat foundations are applied on Staad Pro, the results are as follows;

SUPPORTS 2BON 2BRAFT
Fig. 6: Typical View of the Structure with Elastic Mat Foundation




(1) Base Pressure

BASE 2BPRESSURE
Fig. 7: Base Pressure

A little consideration will show that we are obtaining a maximum base pressure of 73 kPa at the edges, which might be higher than the bearing capacity of the soil. The engineer is advised to review this properly. However, along grid lines 2 and 3, we have base pressure of 39 kPa which increases towards the edges. The minimum base pressure occurred at the mid-spans.

(2) Moment in the raft slab (x-direction)

MX
Fig 8: Mx Moment

The bending moment on the raft slab was discovered to be low when compared with the value that could be gotten using the rigid approach (the maximum moment can be seen to be within 10.9 and 13.4 kNm/m)

final 2Bfront 2Bcover
[Get this publication on Design of Residential Buildings Using Staad, Orion, and Manual Calculations into your mailbox by clicking HERE or by clicking on the book cover]

(3) Bending moment in the Y-direction

MY
Fig 9: My Moment


(4) Twisting Moment

MXY
Fig 10: Mxy Moment

(5) Bending Moment in the ground beams

BENDING 2BMOMENT 2BIN 2BGROUND 2BBEAMS
Fig 11: Bending Moment Diagram of the Ground Beams

For simplicity, the maximum design forces have been presented in the table below;

Maximum 2BMoments

The bending moment on ground beam of grid line (2) is given below;

BM2
Fig 12: Bending Moment for Ground Beam Grid line 2

What are your thoughts on the values obtained? At your convenience, you can model the raft slab on the software you normally use, and compare the results obtained. You can present your observations in the comment section for review.

Thank you for visiting Structville, and God bless you.

Citation
Murthy V.N.S (2012): Textbook of Soil Mechanics and Foundation Engineering (Geotechnical Engineering Series) 1st Edition. CBS Publishers and Distributors Pvt Ltd, India ISBN: 81-239-162-1


Analysis and Design of Box Culvert Using Staad Pro

In our last post, we were able to establish how we can load box culverts properly. If you missed the post, kindly follow the link below to read it;
Loading and Analysis of Box Culverts to Eurocode 2


In this post, we are going to describe how we can model, load, and analyse box culverts using Staad Pro software. Here is a quick a recap of the properties of the box culvert under consideration;

Geometry of the culvert
Total length of culvert = 8 m
Width of culvert c/c of side walls = 2.5 m
Height of culvert c/c of top and bottom slabs = 2.0 m
Length of wing walls = 2.12 m
Thickness of all elements = 300 mm
Thickness of asphalt layer = 70 mm
Materials property
Angle of internal friction of fill soil = 30°
Unit weight of water = 9.81 kN/m3
Unit weight of back fill soil = 19 kN/m3
Unit weight of concrete = 25 kN/m3
Unit weight of asphalt concrete = 22.5 kN/m3
fck = 30 Mpa
fyk = 500 Mpa
Concrete cover = 50 mm
box 2Bculvert 2Bdimensions
Fig 1: Section of the Box Culvert
The steps in modelling the structure on Staad Pro are as follows;

(1) Meshing
Here, the box culvert is idealised with dimensions based on centre to centre of the slabs and walls. This means that the width of the box culvert that will be input into Staad Pro is 2.5 m, while the depth will be 2.0 m. This can be started by forming the nodes in the global XY plane, and then copying and pasting for the length of 8 m in the Z global direction. The output of this operation is as given below;

initial 2Bmeshing
Fig 2: Initial Nodes for Commencement of Modelling

After forming this, the wing walls can also be formed, which is followed by meshing (rectangular or polygonal) to form the shell of the box culvert. The final output of the meshing operation is as shown below;

3D 2Bmeshing
Fig 3: Fully Meshed Box Culvert

The meshing process can be completed by adding plate thickness of 300 mm to all the elements.


(2) Assigning of support conditions/foundations
This is an important aspect of modelling structures. A purely rigid approach will involve using fixed supports, but note that employing 3D model for this purpose will not be very appropriate, but a simple 2D frame model will be better. There are many proposals on how culverts can be modelled as 2D frames, and the reader is advised to consult as many publications as possible. However, to incorporate the effects of soil-structure interaction (to a limited degree though) on our model, we can employ the use of ‘elastic plate‘ foundation option on our model.

If we assume that the supporting soil and the backfill are of the same material, then we can maintain the same angle of internal friction of 30°. Angle of internal friction of 30° can suggest a loose – medium dense sand in its undisturbed state, therefore we can take a modulus of subgrade reaction of 50,000 kN/m2/m for a well compacted sand. For a slightly compacted sand, you can take a value of 30,000 kN/m2/m.

So we can input this option into Staad Pro using ‘compression only’ option (see the dialog box below in Fig 4);

elastic 2Bmat 2Bfoundation 2Bon 2BStaad 2Bpro
Fig 4: Elastic Mat Foundation Dialog Box

When this is applied on the structure, we have the final model as shown in Fig 5;

3D 2BELASTIC 2BMAT 2BFOUNDATION
Fig 5: Application of Elastic Mat Foundation on the Model

What this model (Fig. 5)  means is that every node of the base slab is in contact with the soil, and the soil is represented by a spring of stiffness 50,000 kN/m2/m which is the subgrade modulus of the soil.

front
[Get this publication on Design of Swimming Pools and Underground Water Tanks into your mailbox by clicking HERE or by clicking on the book cover]
Do you wish to get trained in theory of structures, structural design, and use of Staad Pro for analysing and designing of complex structures? Then quickly contact us;
E-mail: info@structville.com or rankiesubani@gmail.com
WhatsApp: +2347053638996 

(3) Loading
This is another important aspect of modelling structures because analysis results would go very wrong if the loads are not applied properly. If you review our previous post, we considered two construction cases;

(A) where the culvert is buried under the soil, and
(B) where is there no earth fill on top of the culvert.

When there is no earth fill on the culvert, the traffic load is directly on the top slab of the culvert as tandem loads and as UDL, but when there is earth fill, traffic load is dispersed in the ratio of 2:1 as UDL on the culvert. We are going to consider 5 load cases in our analysis in this model;

(1) Self weight and other superimposed actions
(2) Vertical earth load
(3) Traffic load
(4) Surcharge load
(5) Horizontal earth pressure load

These  load cases are considered independently at first, and then combined with appropriate partial factors of safety to determine the design actions. Please note that effects of ground water and the pressure in the shell of the culvert when it is filled with water is an important load case too but was not considered in this post. We have determined the magnitude of these loads in our previous post, and we are going to apply them on the box culvert for Case A and Case B.

In our previous post, we were able to analyse the loads on the culvert for Case A as follows;

(1) Self weight: To be calculated automatically by Staad + 1.69 kN/m2 (self weight of asphalt wearing course)
NB: In some cases, the partial factor of safety for self weight of concrete elements and other superimposed dead loads like asphalt wearing course might be different, so in that case, it is very advisable to treat each of them as a separate load case on Staad.
(2) Vertical earth load on the culvert = 22.80 kN/m2
(3) Traffic load dispersed as UDL = 59.523 kN/m2
(4) Horizontal surcharge load = 5.0 kN/m2
(5) Horizontal earth pressure load = trapezoidal distribution with minimum earth pressure of 11.40 kN/m2 at the top of the culvert and 33.25 kN/m2 at the bottom of the culvert

Load cases 1-4 are very easy to apply on Staad by utilising uniform pressure action on plates, however while load cases 1-3 are applied in the global Y direction (gravity loads), load cases 4-5 are applied in the global X direction (horizontally). I will only demonstrate how load case 5 can be applied, since it will involve using the hydrostatic load command on Staad. This is shown on Fig. 6.
HORIZONTAL 2BEARTH 2BPRESSURE 2BLOAD 2BON 2BSTAAD 2BPRO
Fig 6: Application of Horizontal Earth Pressure as Hydrostatic Load on Staad
The steps involved in applying hydrostatic load can be described as follows:
(i) Launch the hydrostatic plate pressure load on Staad
(ii) Select/highlight all the plates that will receive the load
(iii) Input the minimum and maximum pressure loads based on the notation given on Staad (in this case you can see that our maximum pressure load served as W1), and also input the sign conventions properly in order to identify the direction of pressure.
(iii) Put the interpolation direction (in this case we interpolated the load in the Y-direction)
(iv) Add the load case
After all said are done, the output should be as given below;
PRESSURE 2BLOAD
Fig 7: Earth Pressure on the Culvert
At this point, we are going to define the load combination for ultimate limit state as succinctly as possible. The reader is advised to consult the relevant code of practice and standard textbooks on how to group and combine loads involving traffic actions, earth load, etc. In this case, EN 1990, EN 1991-1, EN 1991-2, and EN 1997-1 can be consulted. Note that culverts are also analysed for other loads such as temperature effects, collision on head walls, centrifugal actions, braking actions, etc.
The load combination principle adopted here is based on expressions (6.10a and 6.10b) of EN 1990.
From Table A2.4(B) of EN 1990:2002 + A1:2005, we are going to adopt the following partial factors:
  • All permanent actions including superimposed dead load and vertical earth pressure  γG = 1.35
  • Leading/main traffic action γQ,1 = 1.35
  • Traffic surcharge γQ,2 = 1.5
  • Horizontal earth pressure and ground water γQ,2  = 1.5
We will now apply this load combination on Staad Pro as shown in Fig. 8.

ULS
Fig 8: ULS Load Combination Dialog Box


(4) Analysis
On analysing the structure, we obtain the following results at ultimate limit state;

MX
Fig 9: Transverse Bending Moment
final 2Bfront 2Bcover

[Get this publication on Design of Residential Buildings Using Staad, Orion, and Manual Calculations into your mailbox by clicking HERE or by clicking on the book cover]

Do you wish to get trained in theory of structures, structural design, and use of Staad Pro for analysing and designing of complex structures? Then quickly contact us;
E-mail: info@structville.com or rankiesubani@gmail.com
WhatsApp: +2347053638996 
MY
Fig 9: Longitudinal Bending Moment
BASE 2BPRESSURE
Fig 10: Base Pressure
SQX
Fig 11: Transverse Shear Stress

A little consideration will show that the top slab is subjected to an ultimate design moment MEd of 56.5 kNm/m and an axial pull of NEd of 91.5 kN/m (check for membrane stresses in your own Staad Model).

NB: Shear and membrane forces in plates are expressed in Mpa in Staad, so will you have to multiply them with the thickness of the element to get the values in kN/m.

For the top slab of the culvert, the M-N interaction chart is given in Fig. 11 below for obtaining the design reinforcement.

Column 2Bdesign 2Bchart
Fig 12: M-N Interaction Chart on the Box Culvert

When designed, the area of steel for the axial compression and bending is 1223 mm2/m. Therefore, provision of H16@200 c/c (Asprov = 2010 mm2) on each face will be adequate.

Design Case B: No Earth Fill on Box Culvert
(1) Traffic Loads
When there is no earth fill on the box culvert, all we have to do is to remove the vertical earth fill load, apply direct traffic load on the top slab, and edit the horizontal earth load from trapezoidal to rectangular. As a reminder, the nature of Load Model 1 (which can be used for global and local verification) on the culvert is given in Fig 14 below;

Eurocode 2Btraffic 2Bload 2Bon 2Bculvert
Fig 14: Load Model 1 on the Culvert

The ideal thing is to apply a moving load on Staad after vehicle definition, such that the worst effect can be obtained. Note that you cannot apply a moving load directly on plate elements on Staad Pro, but you will need to create dummy beam members of negligible stiffness so that the axles can sit on them. In this post, we are not going to bother ourselves with the process, but we are going to treat the wheel load as static.

Influence line has shown that the most onerous bending moment is obtained when the front axle is 0.26 m beyond the mid-span of the culvert. Therefore, we are going to apply static wheel load at that location. Remember that it is always recommended to apply the full tandem system of LM1 whenever applicable. The critical location of wheel load on the box culvert for maximum moment is given in Fig 15 below;

most 2Bcritical 2Bwheel 2Bload 2Blocation 2Bon 2Bbox 2Bculvert
Fig 15: Most Critical Location of Wheel Load on the Culvert
When the static traffic load is applied on Staad and viewed longitudinally on the box culvert, we can see it as given in Fig 16 (note that the unloaded areas represent the wing walls).
Traffic 2Bload
Fig 16: Application of LM1 on Staad Pro

(2) Non-traffic Loads
For Design Case A, it is observed that the self weight of the structure and asphalt layer remains the same, the surcharge load also remains the same, but the horizontal earth pressure changes to triangular distribution with a maximum pressure of 21.85 kN/m2 at the bottom of the culvert.


(3) Analysis
When the structure is analysed, the results at ultimate limit state are as shown below;

MX2
Fig 17: Transverse Bending Moment
MY2
Fig 18: Longitudinal Bending Moment
BASE 2BPRESSURE 2B2
Fig 19: Base Pressure

Our analysis results have shown that when there is earth fill, the bending moment at ultimate limit state on top of the culvert is about 56.5 kNm/m, but when traffic is directly on top of the culvert, the bending moment is about 62 kNm/m. This is about 8.8% difference.

Thank you for visiting Structville today, and remember to share this post with your friends and colleagues.

Do you wish to get trained in theory of structures, structural design, and use of Staad Pro for analysing and designing of complex structures? Then quickly contact us;
E-mail: info@structville.com or rankiesubani@gmail.com
WhatsApp: +2347053638996