Dimensional Structural Mass Optimization of A Welded I-Profile Bridge Crane Girder
Dimensional Structural Mass Optimization of A Welded I-Profile Bridge Crane Girder
https://doi.org/10.31803/tg-20191106234034
Nedim Pervan, Adis J. Muminovic, Elmedin Mesic, Mirsad Colic, Vahidin Hadziabdic
Abstract: This paper presents the methodology for the development of an optimization model for the optimization of the cross-section dimensions of a bridge crane girder designed
as a welded I-profile. To carry out this optimization, the CAD/CAE software package CATIA V5 was used. In order to develop an optimization model, a CAD geometrical model
and structural analysis model were developed. Optimization was carried out by the iterative method using a simulated hardening algorithm. Additionally, the optimization process
is carried out by using the PEO (Product Engineering Optimization) CATIA module that contains tools for setting the optimization criteria, design parameters, constraints, and
algorithms. The goal of the optimization is to achieve the minimal mass of the girder, while satisfying all functional and geometrical constraints. As a result of the optimization
process, minimal girder dimensions were obtained and due to that, a minimal amount of material can be used for the manufacturing of the girder.
function of some unknown variables. Optimization can be material of the girder. The maximum displacement should
defined as a mathematical problem of the minimization or not be bigger than L/1000. Additionally, it is necessary to
maximization of the function. check the possibility for lateral torsional buckling [15, 16].
There are many different methods for optimization. The The second step is to collect the data about the
most common method for engineers is the method of optimization problem. A static schematic diagram is shown
simulated hardening. This method is based on the process of in Fig. 2a. The length of the girders is 16 m and the maximal
metal cooling during hardening. It was developed by load of the bridge crane is 20 T.
Kirkpatrick, Gellat and Vecchi [13]. For the numerical finite element analysis and for the
Simulated hardening is the stochastic method for optimization analysis of the girders, the CAD/CAE system
searching for a good approximation of the global minimum CATIA is used. The material of the bridge crane is stainless
of the function. The simulated hardening method is based on steel S235.
the process of the heating and controlled cooling of metal In the third step, it is necessary to define the parameters
materials with the goal of exerting influence on the metal that describe the system. The mass of the girder is in the
structure and inside errors. This algorithm can find the global function of the volume (dimensions) of the girder and
minimum in a certain amount of time [14]. material density. The dimensions of the girder cross section
The algorithm begins with the initial design vector X1 are shown in Fig. 2b, and these dimensions are taken as
(iteration number i = 1) and a high value of the temperature design parameters.
T. By random generation of a new point Xi+1, the nearby first The length of the girder L and the material (density of the
point, the value of the goal function can be calculated by material ρ) are constant during the optimization process. Due
using the following equation: to that, they are removed from the optimization problem
formulation. The same applies to the rails. Rails have
ΔE = Ei +1 − Ei = Δf = fi +1 − fi = f (X i +1 ) − f (X i ) (1) standard dimensions (a = 50 mm and c = 40 mm). The rest of
the girder cross section dimensions (b, t1, h and t) are free
If fi+1 is smaller than fi (negative value of Δf), point Xi+1 parameters and they are taken as the design parameters for
is accepted as a new design vector. If Δf has a positive value, optimization.
point Xi+1 is accepted as a new design vector with the
probability of e−ΔE/kT. That means that if the value of a
randomly generated number is bigger than e−ΔE/kT, point Xi+1
is not accepted as a design vector. With this process, one
iteration of the optimization process is finished.
If the point Xi+1 is not accepted, the generation process of
the new point is based on the Metropolis criteria:
∆E
P (Ei +1 ) = exp (2)
kT
3 DEVELOPMENT OF CAD AND NUMERICAL MODELS where the geometrical and parameter variation are natural"
[17-19].
During the development of the 3D CAD model, the goal The parametric model of a bridge crane girder is
is not just to get the visualization of the design, but also to developed by using the CATIA software package. First, the
get a functional 3D CAD model. A functional 3D CAD 3D model of the bridge crane girder is developed by using
model can be developed by using parametric modelling. the geometrical and dimensional constraints. After the 3D
The term parametric design is associated with the model is developed, parametrization of design parameters,
parametric systems and it is defined as "a process of defined in Chapter 2 of this paper, can be carried out. An
designing with parametric models in a virtual surrounding example of the relation between the parameters is shown in
Fig. 3.
The finite element structural analysis via the Generative the first three or four modes are important in practical
Structural Analysis in the software package CATIA module applications [27]. For the bridge crane girder, only the first
is carried out by using the above-mentioned developed 3D mode is taken into account as a relevant one. This will reduce
CAD parametric model. The first step in the process of the the time of calculations (Fig. 5).
finite element structural analysis is to carry out the
discretization of the girder by using finite elements [20-26].
The second step is to apply the constraints. In this case,
the girder is supported at both ends. One end has all
translations fixed and all rotation allowed, the second end is
a movable support with all translations and rotations allowed
(Fig. 2a).
The third step in the process of the finite element
structural analysis is to apply the loads on the girder. All
loads from real working conditions must be applied to the
girder. The girder is loaded with their own weight, the weight
of the load which it is carrying and the weight of the winch.
The results of the finite element structural analysis,
carried out on the initial 3D CAD model, are the values of the
displacements and Von Mises stresses. These values will be
used as the constraints for optimizations (Fig. 4). Figure 4 Results of a static finite element analysis on the initial 3D model of a
As it has earlier been mentioned, one of the constraints bridge crane girder: a) vector of displacements, b) Von Mises stresses
should be the value for lateral torsional buckling. To get this
value, it is necessary to carry out a new structural analysis for Fig. 5a shows the deformation of the girder without
this case. The static analysis of lateral torsional buckling is lateral bracing during lateral torsion buckling loads. The
taken as a referent analysis. coefficient of deflection for this case is 0.324. Fig. 5b shows
The analysis of lateral torsional buckling in the CATIA the deformation of the girder with three lateral bracings. For
software package gives 10 modes of the torsional twist. Only this case, the coefficient of deflection is 4.45. It can be
m
1
, t1 , h, t , r ) f (b, t1 , h, t ) + r ∑
Φ (b= (8)
j =1 c j = (b, t1 , h, t )
Figure 5 Deformation of the girder; a) without lateral bracing, b) with three lateral
bracings
To solve the problem, a previously made computer
program in FORTRAN was adjusted [29].
4 DEVELOPMENT OF AN OPTIMIZATION MODEL
4.1 The Mathematical Model 4.2 The Numerical Model
For the development of an optimization model, the The same problem was solved numerically, and for that
nonlinear SUMT (Sequential Unconstrained Minimization purpose, the CATIA module Product Engineering Optimizer
Technique) method was used. According to this method, the was used. The module has iterative methods that can be used
so-called penalty functions are introduced. In these for any parameter [30]. In this case, the so-called simulated
processes, approximation is carried out by adding the penalty hardening algorithm is used. This algorithm was chosen
member to the goal function in the area where constraints are because the goal function is unknown, and because this
not satisfied [28]. For the goal function, a girders mass is algorithm can avoid local minimums and find the global
selected, which has the following form: maximum. For this algorithm, the speed of convergence can
be chosen [31-33].
After the definition of the goal function, design
f (b, t1 , h, t )= (2b ⋅ t1 + h ⋅ t + a ⋅ c) ⋅ L ⋅ ρ ⋅ g (4)
parameters and constraints. The results of the optimization
can be presented in the form of diagrams for different
where: b, t1, h, t, a and c - are shown in Fig. 2, ρ – the specific combinations of design parameters. (Fig. 6 and Fig. 7). Fig.
mass of steel, g - speed of gravity. 6a shows the changes of the goal function in the form of the
The constraints that the design should meet: value of mass for all iterations. In the same figure, a diagram
- flanges: of the best solution (the best value of the goal function) is
also shown. Fig. 6b shows the changes of the height of the
σ v − υ ⋅ σ u , p1 ≥ 0 web and the deviation of the values of the deflection from the
σ v − υ ⋅σ u, p2 ≥ 0 allowed values during the iteration process. It can be noticed
(5) that when increasing the values for the height of the web,
t1 − 1.2t ≥ 0 deviation is going close to zero, which means that the
2t − t1 ≥ 0 deflection has lover values. From Fig. 7a, it can be noticed
that the maximal deflection is 15.96 mm, which means that
- the web: the constraints for deflection are fulfilled.
Furthermore, from Fig. 7b, it can be noticed that
maximal stresses occur at the lover part of the girder and that
σ v − υ ⋅σ u, r ≥ 0
the maximal value is 82.4 MPa, which is lower than the
σ ux − υ ⋅ σ s , r ≥ 0 constraint of 160 MPa.
(6)
σ v − σ ux ≥ 0
t −5 ≥ 0
Figure 6 Parameter changes during the iteration process a) relation between the current and best mass of the girder b) relation between deflection and the height of the web
5 RESULTS
Figure 8 Mass of the girder; a) for different lengths of the girder and a different number of laterals bracing for the load of 200 kN, b) in correlation to the allowed deflection for a
length of 16 m, the free length of torsional buckling of 4 m and the ratio of h/t = 120
Fig. 9a shows the optimal values of the mass of the girder Optimizer are the same as the results obtained by the SUMT
for the lengths of 12 m and 16 m and the loads of 50, 100, (Sequential Unconstrained Minimization Technique)
160, 200, 250 and 320 kN. These values of mass result from method.
the structural optimization of the software CATIA. It is The h/t ratio also has an important effect on the mass of
important to notice that the results obtained by the structural the girder. For the same length of the girder, the same free
optimization in the CATIA module Product Engineering length of the torsional buckling and the same loads, a
different mass can be obtained by changing the height and Girders with a high value of the h/t ratio have a problem
thickness of the web. An increase of the h/t ratio has a because buckling can occur. It is necessary to check the
positive impact on the deflection of the girder, and the girder for buckling and if necessary, to add an additional
constraint for the allowed deflection was fulfilled faster. vertical and/or horizontal bracing on the web.
Fig. 9b shows the correlation between the mass of the
girder and the h/t ratio for the load of 200 kN, the allowed
deflection L/1000 and the free length of the buckling of 4 m.
Figure 9 Mass of the girder; a) in correlation to the load of the crane for ratio h/t = 120 and free length of buckling of 4 m, b) in correlation to the ratio h/t for load of 200 kN,
allowed deflection L/1000 and free length of buckling of 4 m
6 CONCLUSION tool which is easy to use and which gives a lot of useful data
about the behavior of the design for different input
The goal of the optimization process carried out in this parameters (dimensions, loads etc.). When using this module,
paper was to find the minimal value of the mass of the bridge it is possible to design lighter and more economical designs.
crane girder designed as a welded I-profile. The mass of the Additionally, an important advance of this module lies in its
girder depends on the input parameters such as the length of fast calculation properties, which reduce the time of the
the girder, the load of the bridge crane, the allowed deflection design or redesign processes. The task which needs to be
and the ratio between the height and the thickness of the carried out by engineers is to develop an optimization model,
girder web. Optimization is carried out for different run the calculation, and analyse the results. Engineers need
combinations of the above-mentioned input parameters. to have some real life experience to precisely develop the
For a numerical analysis and optimization, the CAD/CAE optimization model and to analyse the results in the right
system CATIA V5 was used. A structural numerical analysis manner.
is the first step that needed to be carried out to start with the A disadvantage of this module is the need to have
optimization process. Using the structural numerical analysis computers with height characteristics, the so-called
data for displacements, the stresses and buckling coefficients workstations. With workstations, it is possible to develop
can be obtained at all points of the design. With the help of advance optimization models and to carry out calculations
computers and software, the calculation process is fast. with a small amount of errors.
Additionally, it is easy to control the values of all parameters The optimization methodology developed in this
and to obtain certain conclusions about the behaviour of the research can be used for the optimization of similar designs.
design for the applied constraints. The mass of the girder can be additionally minimized by
It is important to notice that for a successful optimization analysing others shapes of the girder cross section. For the
of the design, it is not enough to know how to work with results presented in this paper, it can be noticed that the
computer software, but it is also very important to understand material at the end of the girder is not used properly. It is
the real constraints of the design during exploitation. This possible to remove some of the material at those places of the
knowledge is important if someone wants to develop an girder. A girder can have a different height at the ends and at
optimization model which is very similar to the real design. the middle.
It is important to know how design is supported, where the
loads are, how to implement some of the standards and 7 REFERENCES
recommendations based on experience. Every optimization
model and optimization process have some errors, but it is [1] Ognjanović, M. (1990). Design Methodology. Faculty of
important to minimize those errors and to develop such a Mechanical Engineering Belgrade, Belgrade.
model that is very close to the real one. [2] Singiresu, S. R. (2009). Engineering Optimization: Theory and
The Product Engineering Optimizer module in Practice. Fourth edition, John Wiley & Sons, Inc., New Jersey.
[3] Muminovic, A. J., Muminovic, A., Mesic, E., Saric, I., &
combination with the finite element method is a powerful Pervan, N. (2019). Spur Gear Tooth Topology Optimization:
Finding Optimal Shell Thickness for Spur Gear Tooth Analysis of the Car Jack. Advances in Science and Technology
produced using Additive Manufacturing. TEM Journal, 8(3), Research Journal, 13(3), 24-30,
788-794. https://doi.org/10.18421/TEM83-13 https://doi.org/10.12913/22998624/109791
[4] Swantner, A. & Campbell, M. I. (2012). Topological and [20] Mešić, E., Avdić, V., Pervan, N., & Repčić, N. (2015). Finite
parametric optimization of gear trains. Engineering element analysis and experimental testing of stiffness of the
Optimization, 44(11), 1351-1368. Sarafix external fixator. Procedia Engineering, 100, 1598-
https://doi.org/10.1080/0305215X.2011.646264 1607. https://doi.org/10.1016/j.proeng.2015.01.533
[5] Gardan, N. & Schneider, A. (2015). Topological optimization [21] Mesic, E., Muminovic, A. J., & Repcic, N. (2012). Geometrical
of internal patterns and support in additive manufacturing. Modelling and Structural Analysis of the Sarafix Fixator
Journal of Manufacturing Systems, 37(1), 417-425. Configurations. Annals of DAAAM for 2012 & Proceedings of
https://doi.org/10.1016/j.jmsy.2014.07.003 the 23rd International DAAAM Symposium, 0069-0074.
[6] Mirzendehdel, A. M. & Suresh, K. (2016). Support structure [22] Muminovic, A. J., Saric, I., & Repcic, N. (2011). Software for
constrained topology optimization for additive manufacturing. Deformity and Stress Condition of 2D Truss Beams Analysis.
Computer-Aided Design, 81, 1-13, Annals of DAAAM for 2011 & Proceedings of the 22nd
https://doi.org/10.1016/j.cad.2016.08.006 International DAAAM Symposium, 1225-1226.
[7] Mezzadri, F., Bouriakov, V., & Qian, X. (2018). Topology [23] Mešić, E., Pervan, N., Repčić, N., & Muminović, A. (2012).
optimization of self-supporting support structures for additive Research of influential constructional parameters on the
manufacturing. Additive Manufacturing, 21, 666-682. stability of the fixator Sarafix. Proceedings of the 23rd
https://doi.org/10.1016/j.addma.2018.04.016 International DAAAM Symposium, 561-564.
[8] Wu, J., Clausen, A., & Sigmund, O. (2017). Minimum [24] Mešić, E., Avdić, V., & Pervan, N. (2015). Numerical and
compliance topology optimization of shell–infill composites experimental stress analysis of an external fixation system.
for additive manufacturing. Computer Methods in Applied Folia Medica Facultatis Medicinae Universitatis Saraeviensis,
Mechanics and Engineering, 326, 358-375. 50(1), 74-80.
https://doi.org/10.1016/j.cma.2017.08.018 [25] Muminovic, A. J., Saric, I., & Repcic, N. (2014). Analysis of
[9] Muminovic, A. J., Mesic, E., Muminovic, A., Colic, M., & Stress Concentration Factors Using Different Computer
Braut, S. (2016). Optimisation of Central Disc Position on Software Solutions. Procedia Engineering, 69, 609-615.
Elastic Rotor Suported with Active Magnetic Bearings. https://doi.org/10.1016/j.proeng.2014.03.033
Journal of Multidisciplinary Engineering Science and [26] Muminovic, A. J., Saric, I., & Repcic, N. (2015). Numerical
Technology (JMEST), 3(9). Analysis of Stress Concentration Factors. Procedia
[10] Fancello, E. A. (2006). Topology optimization for minimum Engineering, 100, 707-713.
mass design considering local failure constraints and contact https://doi.org/10.1016/j.proeng.2015.01.423
boundary conditions, Structural and Multidisciplinary [27] Muminovic, A. J., Braut, S., Muminovic, A., & Saric, I. (2014).
Optimization, 32(3), 229-240. Numerical and Analytical Analysis of Elastic Rotor Natural
https://doi.org/10.1007/s00158-006-0019-9 Frequency. TEM Journal, 3(4), 323-328.
[11] Fancello, E., Haslinger, J., & Feijóo, R. (1995) Numerical [28] Štajer A. (2014). Rešavanje problema minimizacije sa
comparison between two cost functions in contact shape ograničenjima primenom kaznenih funkcija, University of
optimization. Structural Optimization, 9(1), 57-68. Novi Sad, Faculty of Sciences, Novi Sad.
https://doi.org/10.1007/BF01742646 [29] Colic, M., & Colic, M. (2014). Dimensioning and FEM
[12] Hilding, D., Klarbing, A., & Petersson, J. (1999). Optimization analysis of the lever tongs holding the by the friction force.
of structures in unilateral contact. Applied Mechanics Reviews, Proceedings of the 18th International TMT Symposium, 405-
52(4), 139-160. https://doi.org/10.1115/1.3098931 408.
[13] Singiresu, S. R. (2009). Engineering Optimization: Theory and [30] Pervan, N., Čolić, M., Šarić, I., & Hadžiabdić, V. (2016).
Practice. 4th edition, John Wiley & Sons, Inc., New Jersey. Analysis of the Haulage Ropes on Ropeways in Case of
[14] Jasbir, S. A. (2012). Introduction to optimum design. 3rd Accidental Loads. TEM Journal, 5(2), 171-174.
edition, Elsevier, Oxford. [31] Pervan, N., Mešić, E., Čolić, M., & Avdić, V. (2015). Stiffness
[15] Delić, M., Čolić, M., Mešić, E., & Pervan, N. (2017). Analysis of the Sarafix External Fixator based on Stainless
Analytical Calculation and FEM Analysis Main Girder Double Steel and Composite Material. TEM Journal, 4(4), 366-372.
Girder Bridge Crane. TEM Journal, 6, 48-52, [32] Pervan, N., Mešić, E., & Čolić, M. (2017). Stress analysis of
https://doi.org/10.18421/TEM61-07 external fixator based on stainless steel and composite material.
[16] Delić, M., Pervan, N., Čolić, M., & Mešić, E. (2019). International Journal of Mechanical Engineering &
Theoretical and experimental analysis of the main girder Technology, 8(1), 189-199.
double girder bridge cranes. International Journal of Advanced [33] Pervan, N., Mešić, E., Čolić, M., & Avdić, V. (2016). Stiffness
and Applied Sciences, 6(4), 75-80. analysis of the sarafix external fixator of composite materials.
https://doi.org/10.21833/ijaas.2019.04.009 International Journal of Engineering & Technology, 5(1), 20-
[17] Saric I., Pervan N., Muminovic A., & Colic, M. (2018). 26. https://doi.org/10.14419/ijet.v5i1.4912
Development of Integrated Intelligent Cad System for Design
of Shafts. Tehnički vjesnik, 25(Supplement 1), 99-104.
https://doi.org/10.17559/TV-20170521194820 Authors’ contacts:
[18] Saric, I., & Muminovic, A. J. (2018): Development of
integrated intelligent CAD system for synthesis and stress- Nedim Pervan, PhD, Assistant Professor
deformation analysis of pressure vessels. International Journal University of Sarajevo, Faculty of Mechanical Engineering,
of Engineering & Technology, 7(1). Vilsonovo setaliste no. 9,
71000 Sarajevo, Bosnia and Herzegovina
https://doi.org/10.14419/ijet.v7i1.8666
Tel.: +387 33 729 800
[19] Pervan, N., Muminovic, A. J, Muminovic, A., & Delic, M. [email protected]
(2019). Development of Parametric CAD Model and Structural