ONIGBGA PRJ Master 4
ONIGBGA PRJ Master 4
INTRODUCTION
1.1Preamble
Sensitivity analysis is emerging as a fruitful area of engineering research. The reason for this
interest is the recognition of the variety of uses for sensitivity derivatives in its early stage,
sensitivity analysis found its predominant use in assessing the effect of varying parameters in
mathematical model system. More recently there has been strong interest in promoting
systematic structural optimization as a useful tool for the practicing structural design
engineer on large problems a process still underway. Early attempt to use formal
optimization for large structural systems result in excessively long and expensive computer
runs. Examined of the optimization procedure indicated that the predominant contributor to
the cost and time was calculation derivative. As a consequence emerging interest in
sensitivity analysis has emphasized efficient computational procedure.
In addition researchers have developed and applied sensitivity analysis for approximate
analysis, analytical model improvement and assessment of design trends so that structural
sensitivity analysis become more than a utility for optimization and is a versatile design tool
in its own method.
The design of structures for buildings and bridges is mainly concerned with the provision and
support of load-bearing horizontal surfaces. Except in long-span bridges, these floors or
decks are usually made of reinforced concrete, for no other material has a better combination
of low cost, high strength, and resistance to corrosion, abrasion, and fire. The economical
span for a reinforced concrete slab is little more than that at which its thickness becomes just
sufficient to resist the point loads to which it may be subjected or, in buildings, to provide the
sound insulation required. For spans of more than a few meters it is cheaper to support the
slab on beams or walls than to thicken it.
When the beams are also of concrete, the monolithic nature of the construction makes it
possible for a substantial breadth of slab to act as the top flange of the beam that supports it.
1
At spans of more than about 10 m, and particularly where the susceptibility of steel to
damage by fire is not a problem, as for example in bridges and multi-storey car parks, steel
beams become cheaper than concrete beams. It used to be customary to design the steelwork
to carry the whole weight of the concrete slab and its loading; but by about 1950 the
development of shear connectors had made it practicable to connect the slab to the beam, and
so to obtain the T-beam action that had long been used in concrete construction.
(Johnson.2005)
1.2 Limitation
The project is limited to:
Composite steel-concrete structural members have become increasingly popular in the design
and construction of floor systems, structural frames and bridges in the country. In civil
engineering practice, reference to standard technical codes is essential in the design of
structural members. Performing analysis on some member to determine the economic and
safety effects of the change in partial factors of safety in the new Eurocode. This is however,
very inefficient from a computational point of view. Some techniques that take advantage of
the properties of the simple solution are developed to make a sensitivity analysis.
2
1.4 Aim
The aim of this project is to determined the sensitivity analysis of composite beam using EC4
1.5 Objectives
Present a program for optimization design of composite beam using MATLAB
Perform a sensitivity analysis on selected members
Present a design to minimize weight and maximize the performance of the composite
beam
3
CHAPTER TWO
LITERACTURE REVIEW
2.1 Historical Background
Composite construction as we know it today was first used in both building and bridge in
U.S. over a century ago. The first forms of composite structures incorporated the use of steel
and concrete for flexural members, and the issue of longitudinal slip between these elements
was soon identified.
Composite steel–concrete beams are the earliest form of the composite construction method.
In U.S. a patent by an American engineer was developed for the shear connectors at the top
flange of a universal steel section to prevent longitudinal slip. This was the beginning of the
development of fully composite systems in steel and concrete. Concrete-encased steel
sections were initially developed in order to overcome the problem of fire resistance and to
ensure that the stability of the steel section was maintained throughout loading.
The steel section and concrete act compositely to resist axial force and bending moments. A
composite tubular column was developed because they provided permanent and integral
formwork for a compression member and were instrumental in reducing construction time
and consequently costs. They reduce the requirement of lateral reinforcement and costly
tying, as well as provide easier connection to steel beams of a framed structure. (Mosley and
Bungey.2007)
4
Figure 2.1 Component of Composite Construction
Composite slabs have been introduced recently to consider the increase in strength that can
be achieved if the profiled steel sheeting is taken into account in strength calculations.
Composite slabs provide permanent and integral reinforcement, which eliminates the need
for placing and stripping of plywood and timber formwork. More recently, composite slab
and beam systems have been developed for reinforced concrete framed construction; this
provides advantages similar to those attributed to composite slabs for reinforced concrete
slab and beam systems. ( Ray.S.1998)
5
complete structure reduces by about 25 percent. The advantages can be fully utilized as
summarized below:
Faster construction for maximum utilization of rolled and or fabricated components
(structural steel members) and hence quick return of the invested capital.
Reductions in overall weight of structure and thereby reduction in foundation
Reinforced cement concrete (RCC) slab is in compression and steel joist is in tension.
Hence, most effective utilization of the materials can be achieved.
Composite sections have higher stiffness than the corresponding steel sections (in a steel
structure) and thus bending stress as well as deflection are lesser.
Reduced depth allows provision of lower cost for fire proofing of beam’s exposed faces.
Cost of formwork is lower compared to RCC construction.
Cost of handling and transportation is minimized for using major part of the structure
fabricated in the workshop.
Easy structural repair modification maintenance. ( Johnson.2005)
In conventional composite construction, concrete slabs rest over steel beams and are
supported by them. Under load these two components act independently and a relative slip
occurs at the interface if there is no connection between them. With the help of a deliberate
and appropriate connection provided between the beam and the concrete slab, the slip
between them can be eliminated. In this case the steel beam and the slab act as a composite
6
beam and their action is similar to that of a monolithic T- beam. Concrete is stronger in
compression than in tension, and steel is susceptible to buckling in compression. By the
composite action between the two, we can utilize their respective advantages to the fullest
extent. Generally, in steel-concrete composite beams, steel beams are integrally connected to
prefabricated or cast-in-situ reinforced concrete slabs.
8
Therefore, shear connectors are to be designed to cater for integral action of the composite
structure at all load conditions on the following basis:
a) Transmission of longitudinal shear along the contact surface without slip.
b) Prevention of vertical separation of the in-situ RC slab from the pre-fabricated structural
beam.
9
2.7.1 Structural Eurocode Standards
EN 1990: Eurocode : Basis of structural design
EN 1991: Eurocode 1: Actions on structures
EN 1992: Eurocode 2: Design of concrete structures
EN 1993: Eurocode 3: Design of steel structures
EN 1994: Eurocode 4: Design of composite steel and concrete structures
EN 1995: Eurocode 5: Design of timber structures
EN 1996: Eurocode 6: Design of masonry structures
EN 1997: Eurocode 7: Geotechnical design
EN 1998: Eurocode 8: Design of structures for earthquake resistance
EN 1999: Eurocode 9: Design of aluminum structures. (BSEN.2004)
11
There are nine parts of BS5950:
Part 1: Code of practice for design in simple and continuous construction:
Hot rolled sections;
Part 2: Specification for materials, fabrication and erection: hot rolled sections;
Part 3: Design in composite construction;
Section 3.1: Code of practice for design of simple and continuous composite
beams
Section 3.2: Code of practice for design of composite columns and framess
Part 4: Code of practice for design of floors with profiled steel sheeting;
Part 5: Code of practice for design of cold formed sections;
Part 6: Code of practice for design of light gauge sheeting, decking and cladding;
Part 7: Specification for materials and workmanship: cold formed sections;
Part 8: Code of practice for fire resistant design;
Part 9: Code of practice for stressed skin design.
12
Sensitivity analysis was performed to determine the effect of the difference in partial factor
of safety considered in the Eurocode (EN 1990, EN 1992-1-1, EN 1993-1-1 and EN 1994-
11). The ratio of live load to dead load (ranging from 0.2 to 1.0) had been used to vary the
ultimate design load given that correspond to optimization design of composite section
performed for a simply supported composite beam for different combination of span and
loadings. A graph of optimum design against load ratio was plotted.
The model is design to considered all relevant decision variables is that may have an impact
on the cost optimisation of composite beams. These include the following for the steel
section:-
The depth = hw
Steel depth = d
13
The web thickness = tw
Flange = tfl,
tc
tf1
tf2
b2
The moment resistance, shear resistance and the concentrated load bearing resistance of beam
whose plate element are slender, may be significantly influenced by local buckling
constraints. Therefore, the beam cross section can be classified as class 1, 2, 3 or 4 depending
on the ability of the element to resist local buckling. The section are classified by comparing
the slenderness of each compression element with the proper limit from table 5.2 of EC3 (EN
1993-1-1). This is express as
γ=
√
c 235
t fy
for the purpose of this project, the cross section of class 1 will be considered.
The local web buckling constraint can be expressed as t w ≥ Bhw
Where the limiting web slenderness for plastic design has the value
14
1
=83 ԑ ε
B
ε=
√ 235
fy
The limiting plate slenderness for plastic design at the compression flange
1
=10 ԑ
δ
The shear resistance of the section
Av f
V c , RD= , Y mo=1.0 is the material partial factor of safety
y
Y mo √ 3
15 tc
d1 Naf
PNA
d2
d2
Figure 3.1 Plastic stress blocks when PNA lies in flange of steel beam
Mpi, Rd = Nc (d1 + d2 + tc/2) + Naf (d2 + d1/2)
The distance between the bottom of the concrete slab and d PNA,( d 1) is found by equating
the tension force which the compression force
Na−Nc
Naf = 2
Na−Nc
Therefore, d1 = 2b (fy / ymo)
1
The distance d2 between the PNA of the composite section and plastic centroid is
Aa – b 1 d1 b1 ( tf 1−d1 )
d2 = – +d 1−tf 1
2tw tw
beff
tc Nc
16
Naf
d1 Naw
Figure 3.2 Plastic stress blocks when PNA lies within the web
d 1 – tf
Mpe, Rd = Nc (d1 + d2 + tc) + Naf (d2 + d1/2) + Naw ( + d2)
2
Aa – b 1 d1 d 1−tf
d2 = −
2tw 2
C) Plastic neutral axis (PNA) in concrete slab: Na > Nc
beff
tc a Nc
PNA d2
Na
d
tw
Figure 3.3 Plastic stress blocks when PNA lies in concrete slab
a
Mpc, Rd = Nc (d2 + 2 )
17
Where a= depth of the concrete equivalent rectangular stress block
Nc
( )
a= 0 .85 Fck beff . with Yc=1.5
Yc
The distance, d2 between the PNA of the composite section and the plastic centroid is
d a
d2 = 2 + t c – 2
Composite beam
X
Sb
X
Composite beam
Sb
Composite beam
beff
hc
18
hs hp
Floor Details
Beam Span = L
Beam Spacing = Sb
Slap depth = hs
hp
Effective depth of slab= he = (hs 2 )
L
Effective width (beff is lesser of 4 and Sb)
Loads
19
Variable action, = qk
WL
VRD = Design Shear force = 2
Deflection
where w = gk +q k
In order to carry-out the study, four categories of universal sections have been considered
with each of the categories subjected to different loading conditions. The following shows
categories of sections considered:
20
Member Slab Beam Beam Extra Live Design Characteristi
Categor depth length spacing dead load steel c strength
y M M M load KN /m
2
strength of concrete,
KN/ N /mm
2 f cu
2
m
2
N /mm
21
Member Depth Width Thickness Dept Second Plastic Area of
Designation of of h of moment Modulus Section
Sectio Sectio We Flang Web of Area Wply Aa
n n b e d Ixx /Iy
ha b tw tf mm mm3 mm4 mm2
mm mm mm mm
22
20313325 203.2 133.2 5.7 7.8 172.4 2.34107 2.58105 3.20103
3.5 MATLAB
Matlab is a high performance language for technical computing, visualized and program in
an easy to use environment where express in familiar mathematical notation. Matlab requires
the user to provide the input parameters and the variables present in the determined part of
the program.
CHAPTER FOUR
RESULTS AND DISCUSSIONS
23
4.1 Design Moment and Moment Resistance
The design moment, design shear force, actual deflections and their corresponding resistance
and deflection limit as specified by Eurocode are presented in the table below for the
following loading conditions and trial section been considered.
Table 4.1 Applied Bending Moment and Capacity
24
KNm KNm
533 X 210 X 122 1550.2 42489
533 X 210 X 109 1550.2 34131
533 X 210 X 101 1550.2 30212
533 X 210 X 92 1550.2 24106
533 X 210 X 82 1550.2 20673
457 X 152 X 82 695.50 20638
457 X 152 X 74 695.50 16648
457 X 152 X 67 695.50 13801
457 X 152 X 60 695.50 10911
457 X 152 X 52 695.50 8575.1
305 X 127 X 48 251.98 5412.0
305 X 127 X 42 251.98 4116.7
305 X 127 X 37 251.98 3199.0
203 X 133 X 30 108.81 1838.1
203 X 133 X 25 108.81 1312.7
26
The result obtain from the analysis shows that the design moment, shear force and the
actual deflection are less compared to resistance moment, shear resistance and the limit
deflection as specified in EC4
4.2 Optimisation
The optimisation of composite beam of the trial members in term of minimum weight is
presented
Table 4.2 Optimization for Minimizing Weight
27
Table 4.2 show that there is a reduction in the minimum weight of optimisation
calculated for different section of the trial members.
Table 4.3 Load Sensitivity Result for Member size 533 210 92
Weight (KN)
Load Ratio EC4
0.2 913.71
0.4 777.41
0.6 646.11
0.8 526.92
1.0 418. 52
1000
900
f(x) = − 620.425 x + 1028.791
800
Optimum Weight (KN)
700
600
500
400
300
200
100
0
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
Load Ratio
28
Fig 4.1 Sensitivity analysis for 53321092 UB
Table 4.4 Load Sensitivity Result for Member size 305 127 42
0.2 762.72
0.4 613.08
0.6 474.29
0.8 349.97
1.0 246.31
900
800
Optimum Weight (KN)
Load Ratio
29
Sensitivity analysis showed that at little variation of live loads, the composite beam
designed to EC4 attained the target optimum design up to the worst loading scenarios
where the live load equals the dead load.
From the graph a linear equation of optimum weight was gotten so as to calculate the actual
amount optimum weight.
Table 4.5 Actual Amount Optimum Weight
Load Optimum
Ratio Weight (KN)
0.2 903.92
0.4 779.81
0.6 655.76
0.8 531.68
1.0 407.60
CHAPTER FIVE
CONCLUSIONS AND RECOMMENDATIONS
5.1 CONCLUSIONS
The following conclusion can be drawn from the result obtained in this research
1. Findings from the research depicts that composite beam designed to EC4 is less
susceptible to failure as it possess higher safety in bending.
2. The analysis carryout help us to minimize weight and maximized the performance of
the composite beam using MATLAB
3. The sensitivity analysis showed that live load to dead load ratio is inversely
proportional to the optimisation design of the system.
5.2 RECOMMENDATIONS
30
For improvement in the research the following recommendation were made
1. Other methods of Sensitivity analysis should be used to compare the results obtained in
this research.
2. In structures where high risk of failure is probable, the EC4 should be used.
3. The Nigerian design engineers should embrace the use of EC4 due to the increased
safety index and the cost reduction offered by new code.
REFERENCES
British Standards Institution (2004). Design of Composite Steel and Concrete Structures-
General rules and rules for buildings. London, BS EN 1994-1-1: 2004.
British Standards Eurocode (2002). Basis of Structural Design, BS EN 1990:2002.
Chung K.F and Lawson R.M (1994). Composite Beam Design to Eurocode. The Steel
construction Institute, London.
Chung K.F and Lawson, R.M (2000). Simplified Design of composite Beams, Eurocode 4.
Kenny T.C (2009). Comparison of composite beam design between British Standard part 3.1
and Eurocode 4 part 1.1. Unpublished B.Eng. Thesis. University Technology, Malaysia.
Nowak A.S and Collins, K.R (2000). Reliability of Structures.
Tahir M.M (2008). Experimental test on composite and non-composite connection using
trapezoidal web profiled steel section. Journal of steel structure,8: 43-58.
Singiresu S.RAO. Engineering Optimization Theory and Practice 4th Edition.
31
Jasbri S. Arora (Elsevier, 2004), Introduction to Optimum Design 2nd Edition.
R.P. JOHNSON, (2005). Volume1, Composite Structures of Steel and Concrete 2nd Edition.
Daniel A. Tortorelli, Design Sensitivity Analysis. Overview and Review, University of
Illinois at Champaign-Urbana.
William R. Spillers and Keith M. MacBain Structural Optimization.
B.Davison and W.Owens, (2003), Steel Designerss Manual,6th Edition.
B.Mosley and J.Bungey (2007), Reinforced Concrete Design.
www.iitg.ernet.in
www.shodhganga.inflibnet.ac.in
www.design12.weebly.com
APPENDIX
32
% Y = unit weight of concrete
% b = effective breadth
% gk = dead load
% qk = live load
% N = design load
% Ea = extra dead load
% output:
% Med = design moment
***************************************************************************
***************************************************************************
**************************************************************************
>> Med
length of beam in m:15
33
Unit weight of concrete:23.6
Slab thickness in m:0.15
Decking profile height in m:0.06
Extra dead load in m:2.5
Ea =
2.5000
gk =
5.3320
qk =
N = 55.1183
34
>> Med
Length of beam in m:12
Unit weight of concrete:23.6
Slab thickness in m:0.15
Decking profile height in m:0.06
Extra dead load in m:2.5
Ea =
2.5000
gk =
5.3320
qk =
N=
39.5946
35
design moment in (KNm):
712.7028
***************************************************************************
***************************************************************************
**************************************************************************
>> Med
length of beam in m:9
Unit weight of concrete:23.6
Slab thickness in m:0.13
Decking profile height in m:0.06
Extra dead load in m:2.5
Ea =
2.5000
gk =
4.8600
qk =
36
N=
24.8873
Ea =
2.5000
gk =
4.7420
37
qk =
2.5000
N=
17.7655
38
% b = effective breadth
% gk = dead load
% qk = live load
% N = design load
% Ea = extra dead load
% output:
% Vrd = design shear force
***************************************************************************
***************************************************************************
***************************************************************************
>> Vrd
length of beam in m:15
39
Unit weight of concrete:23.6
Slab thickness in m:0.15
Decking profile height in m:0.06
Extra dead load in m:2.5
Live load in kN/m^2:5
>> Vrd
length of beam in m:12
Unit weight of concrete:23.6
Slab thickness in m:0.14
Decking profile height in m:0.06
Extra dead load in m:2.5
Live load in kN/m^2:4
40
Decking profile height in m:0.06
Extra dead load in m:2.5
Live load in kN/m^2:3
41
function [moment capacity in the beam flange] = Mc
(Nc,d1,d2,tc,Na,Naf,Fy,Ymo,b1,A,tf1,tw,b,fck)
% Mc calculate the moment capacity in the beam flange
% Detailed explanation goes here
% input:
% Na = resistance of steel section
% Nc = resistance of concrete section
% Naf = resistance of the steel flange above the neutral axis
% tc = slab thickness
% fy = characteristic strenght of steel
% fck = characteristic strenght of concrete
% tf1 = flange thickness
% tw = web thickness
% b = effective breadth
% A = area of steel section
% Ymo = factor of safety of steel
% b1 = width of steel section
% d1 = distance between bottom of the concrete slab and PNA
% d2 = distance between PNA of the composite section and plastic centroid
% output:
% Mc = moment capacity in the beam flange
42
Ymo = input('factor of safety of steel:');
b1 = input('width of steel section in mm:');
d1 = (Na - Nc)/(2*b1*(fy/Ymo));
h = input ('slab thickness in mm:');
tf1 = input('flange thickness in mm:');
tw = input('web thickness in mm:');
d2 = A+d1*(b1+2*tw)+2*tf1*(tw-b1);
disp(' ');
disp('moment capacity in the beam flange in (Nmm):');
disp(Nc*(d1+d2+h/2)+Naf*(d2+d1/2));
end
***************************************************************************
***************************************************************************
***************************************
>> Mc
Characteristic strength of steel in N/mm^2:355
Area of steel section in mm^2:1.55*10^4
Characteristic strength of concrete in N/mm^2:35
Effective breadth in mm:3750
Factor of safety of steel:1
Width of steel section in mm:211.9
Slab thickness in mm:150
Flange thickness in mm:21.3
Web thickness in mm:12.7
43
***************************************************************************
***************************************************************************
***************************************************************************
>> Mc
Characteristic strength of steel in N/mm^2:355
Area of steel section in mm^2:1.39*10^4
Characteristic strength of concrete in N/mm^2:35
Effective breadth in mm:3750
Factor of safety of steel:1
Width of steel section in mm:210.8
Slab thickness in mm:150
Flange thickness in mm:18.8
Web thickness in mm:11.6
>> Mc
Characteristic strength of steel in N/mm^2:355
Area of steel section in mm^2:1.29*10^4
Characteristic strength of concrete in N/mm^2:35
Effective breadth in mm:3750
Factor of safety of steel:1
Width of steel section in mm:210
Slab thickness in mm:150
Flange thickness in mm:17.4
44
Web thickness in mm:10.8
45
Width of steel section in mm:208.8
Slab thickness in mm:150
Flange thickness in mm:13.2
Web thickness in mm:9.6
>> Mc
Characteristic strength of steel in N/mm^2:355
Area of steel section in mm^2:1.05*10^4
Characteristic strength of concrete in N/mm^2:30
Effective breadth in mm:3000
Factor of safety of steel:1
Width of steel section in mm:155.3
Slab thickness in mm:140
Flange thickness in mm:18.9
Web thickness in mm:10.5
46
Area of steel section in mm^2:9.45*10^3
Characteristic strength of concrete in N/mm^2:30
Effective breadth in mm:3000
Factor of safety of steel:1
Width of steel section in mm:154.4
Slab thickness in mm:140
Flange thickness in mm:17
Web thickness in mm:9.6
47
***************************************************************************
***************************************************************************
***************************************************************************
>> Mc
Characteristic strength of steel in N/mm^2:355
Area of steel section in mm^2:7.62*10^3
Characteristic strength of concrete in N/mm^2:30
Effective breadth in mm:3000
Factor of safety of steel:1
Width of steel section in mm:152.9
Slab thickness in mm:140
Flange thickness in mm:13.3
Web thickness in mm:8.1
48
Moment capacity in the beam flange in (Nmm):
8.5751e+09
***************************************************************************
***************************************************************************
***************************************************************************
>> Mc
Characteristic strength of steel in N/mm^2:275
Area of steel section in mm^2:6.12*10^3
Characteristic strenght of concrete in N/mm^2:25
Effective breadth in mm:2250
Factor of safety of steel:1
Width of steel section in mm:125.3
Slab thickness in mm:130
Flange thickness in mm:14
Web thickness in mm:9
49
Width of steel section in mm:124.3
Slab thickness in mm:130
Flange thickness in mm:12.1
Web thickness in mm:8
50
Characteristic strenght of concrete in N/mm^2:20
Effective breadth in mm:1750
Factor of safety of steel:1
Width of steel section in mm:133.9
Slab thickness in mm:125
Flange thickness in mm:9.6
Web thickness in mm:6.4
51
Program for Shear Resistance
***************************************************************************
***************************************************************************
***************************************************************************
52
Program for Actual Deflection
function [ actual deflection] = Def( W,gk,qk,Y,he,hs,hp,l,b,Ia,Ea,E)
% Def calculate the actual deflection
% Detailed explanation goes here
% input:
% l = lenght of beam
% hs = slab thickness
% hp = decking profile height
% he = effective depth of slab
% Y = unit weight of concrete
% b = effective breadth
% gk = dead load
% qk = live load
% W = design load at serviceability limit state
% Ea = extra dead load
% Ia = second moment of area
% E = Elastic modulus ofb steel
% output:
% Def = actual deflection
53
W = (gk + qk)*b
Ia = input('second moment of area in m^4:');
E = 205
disp(' ')
disp('actual deflection in (m):')
disp(5*W*l^4/384*E*Ia)
end
55