Open access peer-reviewed chapter - ONLINE FIRST

Perspective Chapter: Simplified Matrix Calculation for Analysis of Girder-Deck Bridge Systems

Written By

Alvaro Gaute-Alonso and David Garcia-Sanchez

Submitted: December 17th, 2021 Reviewed: December 22nd, 2021 Published: February 10th, 2022

DOI: 10.5772/intechopen.102362

From the Edited Volume

Applied Methods in Design and Construction of Bridges, Highways and Roads - Theory and Practice [Working Title]

Dr. Khaled Ghaedi

Chapter metrics overview

View Full Metrics

Abstract

In the design of girder-deck bridge systems, it is necessary to determine the cross-sectional distribution of live loads between the different beams that make up the cross section of the deck. This article introduces a novel method that allows calculating the cross-sectional distribution of live loads on beam decks by applying a matrix formulation that reduces the structural problem to 2 degrees of freedom for each beam: the deflection and the rotation of the deck slab at the center of the beam’s span. To demonstrate the proposed method, the procedures are given through three different examples by applying loads to a bridge model. Deflection, bending moment, and shear force of the bridge girders are calculated and discussed through the given examples. The use of the proposed novel method of analysis will result in significant savings in material resources and computing time and contributes in the minimization of total costs, and it contributes in the smart modeling process for girder bridge behavior analysis allowing to feed a bridge digital twin (DT) model based on Inverse Modeling holding the latest updated information provided by distributed sensors. The presented methodology contributes also to speed up real-time decision support system (DSS) demands.

Keywords

• girder bridge decks
• optimized matrix method
• structural grillage models

1. Introduction

Girder bridge decks are a structural typology commonly used in the design of road and railway bridges, and therefore, any optimization in their calculation has a depth impact on the project phase. As explained in [1], the design process relies solely on the designers’ experience, intuition, and ingenuity resulting in a depth cost material, time, and human effort. It is very common to use structural grillage models [2, 3, 4, 5, 6, 7, 8, 9, 10] to calculate the cross-sectional distribution of live loads between the different beams that make up the cross section of the deck. Another way to deal with the design and calculation of such decks would be to apply different formulations contained in the bridge design standards that allow approximating the cross-sectional distribution of the bending moment and shear stress caused by live loads through what is known as the load distribution factor “LDF” [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. The “LDF” associated with each case study is conditioned by the type and number of beams, their spacing and length, as well as the existence or not of transverse diaphragms that bring transverse rigidity to the deck.

However, the LDF does not allow determining the distribution of bending stresses on all beams; therefore, the design is oversized. The need for a method to calculate the cross-sectional distribution across all types of beams, without resorting to complex structural grillage models or finite element models in specific structure calculation programs, or to approximate parametric methods based on the “LDF,” is one of the authors’ motivations for the development of the research work that has given rise to this article.

2. Traditional methods for girder bridge deck analysis

Structural grillage models began to be used for the analysis of cross-sectional distribution on beams in the 1960s. These models divide the beam deck into longitudinal and transverse beams (Figure 1). Longitudinal beams are responsible for providing the longitudinal bending stiffness of the deck, considering as many longitudinal beams as beams conform to the analyzed beam deck. The structural section of each of the longitudinal beams shall be the result of the section composed of the beam analyzed and the effective depth of the contributing deck with that beam [24, 25]. The cross-sectional distribution of the structural model is provided by the cross-beams and the torsional stiffness of the longitudinal beams. The structural section of the cross-beams corresponds to a rectangular section with an equivalent depth to the slab thickness and a depth according to the discretization used in the grillage models. The analysis of such structural models involves the use of specific structure calculation programs that provide the computing power necessary for the resolution of the proposed matrix problem.

The use of structural grillage models allows obtaining the structural response of girder bridge decks to live loads, adequately manifesting the cross-sectional distortion of the girder bridge deck and the distribution of stresses on each of the longitudinal beams that make up the deck. However, such models involve complex, time-consuming analysis that involves the need to use specific structure calculation programs, making it necessary to use simplified methods for the start of the design process.

The concept of LDF was first introduced using empirical formulas at the American Association of State Depthway Officials (AASHO) in 1931 [26]; these methods propose calculating the cross-sectional distribution on beam bridge decks roughly. The LDF is calculated from a series of formulations that parametrically treat the calculation of the percentage of bending moment and shear stress supported by the most requested longitudinal beam. The parameters that condition the calculation of the LDF are the depth of the beam, the span length, the spacing, the number of beams, the position of the load, and the beam position. This method provides an approximate value of the maximum bending stresses on the beams but is not able to reproduce the cross-sectional distribution of longitudinal bending between all the beams that make up the deck.

3. Proposed method for the study of the cross-sectional load distribution on a girder bridge deck

The authors propose the use of a method that allows obtaining the cross-sectional distribution of live loads on girder bridge decks without using empirical formulas of LDFs or complex structural models involving the use of specific structure calculation programs. The proposed method is based on using a virtual model that reflects the transverse stiffness of the slab deck, supported on a series of springs that provide the flexural stiffness Eq. (1) and torsional stiffness Eq. (2) of the longitudinal beams that make up the girder bridge deck (Figure 2).

Kv,i=48·EIiL3E1
Kt,i=2·GJiLE2

where EIi = longitudinal bending stiffness of beam “i”; GJi = longitudinal torsional stiffness of beam “i”; L = distance between bridge supports.

The proposed method considers 2 degrees of freedom for each longitudinal beam on the bridge deck: (1) the deflection and (2) the rotation of the deck slab at the center of the beam’s span. Figure 2 represents the structural model scheme for a girder bridge deck composed of three longitudinal beams. The matrix approach that solves the structural problem of the cross-sectional distribution of live loads between the different beams that make up the deck is raised in matrix Eqs. (3), (4), (5). The stiffness of the springs that represents the longitudinal bending and the torsional stiffness of the longitudinal beams that make up the model corresponds to the stiffness equivalent to the center of the span.

Ke=EIeLe3126Le6L4Le2126Le6Le2Le2126Le6Le2Le2126Le6Le4Le2=K11,eK12,eK21,eK22,eE3
KBi=kv,i00kt,iE4
P1Mt,1P2Mt,2P3Mt,3=K11,1+KB1K12,10000K21,1K22,1+K11,2+KB2K12,20000K21,2K22,2+KB3×fv,1θ1fv,2θ2fv,3θ3E5

Where EIe = transverse bending stiffness of element “e” of the upper deck slab; Le = length of element “e” of the upper deck slab; fv,i = vertical displacement experienced by longitudinal beam “i”; Ɵi = transverse rotation of the bridge deck over the longitudinal beam “i.”

3.1 Loads applied in the center of the span of the longitudinal beams

The maximum longitudinal bending stress in each of the longitudinal beams that make up the bridge deck corresponds to the application of a point load in the center of the beam span. The distribution of the maximum bending moment and maximum shear stress in the different longitudinal beams is obtained by Eqs. (6) and (7), respectively.

Mfmáx,n=Q·L4·fv,ni=1Nfv,iE6
Qmáx,n=Q·L2·fv,ni=1Nfv,iE7

where Q = point load value; L = girder bridge span length.

Example 3.1.1.

In the study of the structural behavior of the bridge deck represented in Figure 3, it is desired to know the distribution of the bending moment and the shear force in each of the longitudinal beams generated by the application of the following load states: (a) a point load of 300 KN in the center of the span of beam 1, and (b) a point load of 300 KN in the center of the span of beam 2. The girder bridge is made up of four longitudinal beams 1.05 meters deep (two end beams and two central beams) and an upper slab 0.25 meters thick and 16.83 meters wide. Considering the effective width of the upper slab in each of the longitudinal beams, the inertia to the longitudinal bending of the end and central beams takes the value of 0.1445 m4 and 0.1904 m4, respectively. The longitudinal torsional inertia of the longitudinal beams takes the value of 4.3·10−3 m4. The spacing between longitudinal beams is 5.13 meters, while the distance between the support devices of each longitudinal beam is 25 meters. Both the upper slab and the longitudinal beams are made of structural concrete whose modulus of elasticity reaches 35,000 MPa.

The stiffness of the springs on which the upper slab of the bridge deck rests, and which simulates the bending and torsional stiffness of the longitudinal beams (Figure 4), is calculated as follows:

End beams

Kv,end=48·EIendL3=48·3.5·107KNm2·0.1445m425m3=15,537KNm
Kt,end=2·GJendL=2·3.5·107KNm22·1+0.2·4.3·103m425m=5,017KN·m
KB,end=kv,end00kt,end=15,537KNm005,017KN·m

Central beams

Kv,c=48·EIendL3=48·3.5·107KNm2·0.1904m425m3=20,472KNm
Kt,c=2·GJendL=2·3.5·107KNm22·1+0.2·4.3·103m425m=5,017KN·m
KB,c=kv,end00kt,end=20,472KNm005,017KN·m

The mechanical characteristics of the beam elements that simulate the transverse distribution provided by the upper slab (Figure 4) are calculated as follows:

I=112·20m·0.25m3=2.6·102m4
L=5,13m
K=E·IL3126L6L4L2126L6L2L2126L6L2L2126L6L4L2
K=81,020KNm207,800KN207,800KN710,690KN·m81,020KNm207,800KN207,800KN355,340KN·m81,020KNm207,800KN207,800KN355,340KN·m81,020KNm207,800KN207,800KN710,690KN·m
K=K11K12K21K22

The global stiffness matrix used in the proposed method for the analysis of the transverse distribution of live loads in girder bridge decks is obtained as follows:

KG=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end

The transversal distribution of each load state is obtained by planting the compatibility between loads and displacements:

300KN0000000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end×xfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Deflections (Figure 5) and rotations

Maximum bending moments

Mfmáx,n=Q·L4·fv,ni=1Nfv,i=1,588KN·m499KN·m25KN·m187KN·m

Maximum shear forces

Qmáx,n=Q·L2·fv,ni=1Nfv,i=127KN40KN2KN15KN

00300KN00000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end×xfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Deflections (Figure 6) and rotations

Maximum bending moments

Mfmáx,n=Q·L4·fv,ni=1Nfv,i=566KN·m859KN·m478KN·m28KN·m

Maximum shear forces

Qmáx,n=Q·L2·fv,ni=1Nfv,i=45KN69KN38KN2KN

3.2 Loads applied at a distance “x” from one of the supports of the longitudinal beams

The proposed method is applicable to calculate the structural response of a girder bridge deck to the application of a vertical load in any cross section. Using the Maxwell-Betti reciprocity theorem [27], for the application of a point load of value “Q” at a distance “x” from one of the two support points of a longitudinal beam, the deflection at the center of the span is obtained by applying the formulation reflected in Eq. (8). Likewise, the distribution of the maximum bending moment and maximum shear stress in the different longitudinal beams is obtained by Eqs. (9) and (10), respectively.

fcl,Qx,i=fv,i·sinπ·xLE8
Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i/xL/2E9
Qmáx,n=Q·LxL·fv,ni=1Nfv,i/xL/2E10

Example 3.2.1.

In the study of the structural behavior of the bridge deck analyzed in Example 3.1.1., it is intended to know the deflection in the center of the span of the longitudinal beams, as well as the distribution of the maximum bending moment and the maximum shear force in each of the longitudinal beams generated by the application of the following load states: (a) a point load of 300 KN at a distance equivalent to L/3 from one of the supports of beam 1, and (b) a point load of 300 KN at a distance equivalent to L/3 from one of the supports of beam 2 (Figure 7).

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=13.12mm4.12mm0.20mm1,55mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=1,412KN·m443KN·m22KN·m167KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=169KN53KN3KN20KN

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=4.12mm6.25mm3.48mm0.20mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=503KN·m764KN·m425KN·m25KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=60KN92KN51KN3KN

3.3 Loads applied at any point on the bridge deck

The proposed method is applicable to calculate the transverse response of girder bridge decks to the application of a vertical load at any point of the cross section. As in any matrix calculation, if the vertical load acts on a section of slab between longitudinal beams, the degrees of freedom of the structural model are locked and the reactions in the locked degrees of freedom are calculated (rigid step). Subsequently, the degrees of freedom are released and loaded with the reactions obtained in the previous step to obtain the vertical displacement of each of the longitudinal beams that make up the bridge deck (flexible step).

Example 3.3.1.

In the study of the structural behavior of the bridge deck analyzed in Example 3.1.1., it is intended to know the deflection in the center of the span of the longitudinal beams, as well as the distribution of the maximum bending moment and the maximum shear force in each of the longitudinal beams generated by the application of the load states described in Figure 8.

Rigid step (Figure 9).

The load generated by the action of the truck axles on longitudinal beams 1 and 2 is obtained considering the compatibility between loads and movements in the upper slab between longitudinal beams 1 and 2.

VLB1MLB1Qi0Qi0VLB2MLB2=K11,1K12,1K21,1K22,1+K11,200000000K12,200000000K21,200000000K22,2+K11,3K12,3K21,3K22,3·00v2Ɵ2v3Ɵ300

Load on longitudinal beams 1 and 2

Q=1.46·QKN1.09·QKN·m0.54·QKN0.70·QKN·m

Flexible step (Figure 10).

1.46Q1.09Q0.54Q0.70Q0000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,endxfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Truck axle 1

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=2.39mm1.20mm0.21mm0.30mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=268KN·m135KN·m24KN·m34KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=31KN16KN3KN4KN

Truck axle 2

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=4.27mm2.15mm0.38mm0.53mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=468KN·m236KN·m41KN·m59KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=39KN20KN3KN5KN

Truck axle 3

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=3.69mm1.86mm0.33mm0.46mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=409KN·m206KN·m36KN·m51KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=41KN21KN4KN5KN

Truck axle 4

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=3.43mm1.73mm0.30mm0.43mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=385KN·m194KN·m34KN·m48KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=45KN22KN4KN6KN

Deflection and efforts generated by the action of the four-axle truck

fv,i=13.78mm6.94mm1.22mm1.72mm
Mfmáx,i=1,323KN·m667KN·m117KN·m165KN·m
Qmáx,i=138KN70KN12KN17KN

4. Conclusions

The proposed method for the analysis of the cross-sectional distribution of live loads on girder bridge decks allows determining the cross-sectional distribution in different configurations of girder bridge decks without the need to resort to complex calculation models that involve depth computing power and excessive analysis time. The method is also applicable for modern synthetic materials, such plastic composites, self-repair. The simplicity of the method allows an easy integration into optimal bridge design strategies [28] or more heuristic approaches [29, 30, 31, 32, 33] to challenge today’s competitive world intelligently.

Acknowledgments

This work has received funding from the European’s Union Horizon 2020 research and innovation program under the grant agreement No 769373 (FORESEE project). This paper reflects only the author’s views. The European Commission and INEA are not responsible for any use that may be made of the information contained therein.

References

1. 1. Ahsan R, Rana S, Sayeed Nurul G. Cost optimum design of posttensioned I-girder bridge using global optimization algorithm. Journal of Structural Engineering. 2012;138(2):273-284. DOI: 10.1061/(ASCE)ST.1943-541X.0000458
2. 2. Douthe C, Caron J, Baverel O. Gridshell structures in glass fibre reinforced polymers. Construction and Building Materials. 2010;24(9):1580-1589. DOI: 10.1016/j.conbuildmat.2010.02.037
3. 3. Peloux LD, Tayeb F, Lefevre B, Baverel O, Caron J-F. Formulation of a 4-DoF torsion/bending element for the formfinding of elastic gridshells. Amsterdam: Proceedings of the International Association for Shell and Spatial Structures (IASS) Symposium; 2015
4. 4. Rombouts J, Lombaert G, Laet LD, Schevenels M. A novel shape optimization approach for strained gridshells: Design and construction of a simply supported gridshell. Engineering Structures. 2019;192:166-180. DOI: 10.1016/j.engstruct.2019.04.101
5. 5. D. Veenendaal and. P. Block. An overview and comparison of structural form finding methods for general networks. International Journal of Solids and Structures. 2012;49:3741-3753. DOI: 10.1016/j.ijsolstr.2012.08.008
6. 6. Samartin Quiroga A. Notas al cálculo de esfuerzos en tableros de puentes. Hormigón y Acero. 1971;22(98):115-135
7. 7. Manterola J. Cálculo de tableros por el método del emparrillado. Hormigón y Acero. 1977;28(122):93-149
8. 8. Jing-xian S. Application of T-beam Grillage model in reconstruction design of dangerous bridge. 4th International Conference on Mechanical, Control and Computer Engineering (ICMCCE); 2019. DOI: 10.1109/ICMCCE48743.2019.00190
9. 9. Zhou Y, Ji Y. Comparison and Analysis of the Results of Grillage Method and Single Beam Method to Continuous Box Girder with Variable Width. Banff, Canada: 4th International Conference on Transportation Information and Safety (ICTIS); 2017
10. 10. Feng JP, Zhang XL, Zhu Z, Huang PM, Wang D, Niu YW. Application on Grillage Method in Box Girder Bridges. Applied Mechanics and Materials. 2013;1:886-890. DOI: 10.4028/www.scientific.net/amm.438-439.886
11. 11. Dicleli M, Erhan S. Live Load Distribution Formulas for Single-Span Prestressed Concrete Integral Abutment Bridge Girders. Journal of Bridge Engineering. 2009;14:6. DOI: 10.1061/(ASCE)BE.1943-5592.0000007
12. 12. Fu CC, Elhelbawey M, Sahin MA, Schelling DR. Lateral distribution factor from bridge field testing. Journal of Structural Engineering. 1996;122(9):1106-1109. DOI: 10.1061/(ASCE)0733-9445(1996)122:9(1106)
13. 13. Gheitasi A, Harris DK. Overload flexural distribution behavior of composite steel girder bridges. Journal of Bridge Engineering. 2015;20(5) 04014076:1-15. DOI: 10.1061/(ASCE)BE.1943-5592.0000671
14. 14. Hess S, Filosa F, Ross BE, Cousins TE. Live Load Testing of NEXT-D Bridges to Determine Distribution Factors for Moment. Journal of Performance of Constructed Facilities. 2020;34(4) 04020063:1-9. DOI: 10.1061/(ASCE)CF.1943-5509.0001452
15. 15. J. Huang and. J. Davis. Live load distribution factors for moment in NEXT beam bridges. Journal of Bridge Engineering. 2018;23(3). 06017010:1-7. DOI: 10.1061/(ASCE)BE.1943-5592.0001202
16. 16. Huo XS, Wasserman EP, Zhu P. Simplified method of lateral distribution of live load moment. Journal of Bridge Engineering. 2004;9(4):382-390. DOI: 10.1061/(ASCE)1084-0702(2004)9:4(382)
17. 17. Idriss RL, Liang Z. In-service shear and moment girder distribution factors in simple-span prestressed concrete girder bridge. Journal of the Transportation Research Board. 2010;2172:142-150. DOI: 10.3141/2172-16
18. 18. Kim YJ, Tanovic R, Wight G. Load configuration and lateral distribution of NATO wheeled military trucks for steel I-girder bridges. Journal of Bridge Engineering. 2010;15(6):740-748. DOI: 10.1061/(ASCE)BE.1943-5592.0000113
19. 19. Semendary AA, Steinberg EP, Walsh KK, Barnard E. Live-load moment-distribution factors for an adjacent precast prestressed concrete box beam bridge with reinforced UHPC shear key connections. Journal of Bridge Engineering. 2017;22(11):04017088/1-04017088/18
20. 20. Yalcin OF, Diclel M. Comparative study on the effect of number of girders on live load distribution in integral abutment and simply supported bridge girders. Advances in Structural Engineering. 2013;16(6):1011-1034. DOI: 10.1260/1369-4332.16.6.1011
21. 21. Harris DK. Assessment of flexural lateral load distribution methodologies for stringer bridges. Engineering Structures. 2010;32(11):3443-3451. DOI: 10.1016/j.engstruct.2010.06.008
22. 22. Kong S, Zhuang L, Tao M, Fan J. Load distribution factor for moment of composite bridges with multi-boxgirders. Engineering Structures.215.110716. 2020;215(110716):1-19. DOI: 10.1016/j.engstruct.2020.110716
23. 23. Terzioglu T, Hueste MBD, Mande JB. Live Load Distribution Factors for Spread Slab Beam Bridges. Journal of Bridge Engineering. 2017;22(10) 04017067:1-15. DOI: 10.1061/(ASCE)BE.1943-5592.0001100
24. 24. G. o. S. Ministry of Public Works. 18.2.1. Ancho eficaz del ala en piezas lineales. EHE - 08. Instrucción de Hormigón Estructural; 2011
25. 25. G. o. S. Ministry of Public Works. 4.5. Anchura eficaz elástica. RPX-95. Recomendaciones para el proyecto de puentes mixtos para carreteras; 1996
26. 26. American Association of State Highway and Transportation Officials. Standard specifications for highway bridges. 1st ed. Washington, DC; 1931
27. 27. Baker WF, Beghini LL, Mazurek A, Carrion J, Beghini A. Maxwell’s reciprocal diagrams and discrete Michell frames. Structural and Multidisciplinary Optimization. 2013;48(2):267-277. DOI: 10.1007/s00158-013-0910-0
28. 28. Kuang Y, Ou J. Self-repairing performance of concrete beams strengthened using superelastic SMA wires in combination with adhesives released from hollow fibers. Smart Materials and Structures. 2008;17. DOI: 10.1088/0964-1726/17/2/025020
29. 29. Ghani SN. Performance of global optimization algorithm EOP for non-linear non-differentiable constrained objective functios. New York: Proceedings of IEEE International Conference on Evolutionary Computation; 1995
30. 30. Ghani S. A versatile algorithm for optimization of a nonlinear non-differentiable constrained objective function. UKAEA Harwell. R-13714. HMSO Publications Centre; 1989
31. 31. Hassanain MA, Loov RE. Cost optimization of concrete bridge infrastructure. Canadian Journal of Civil Engineering. 2003;30(5):841-843. DOI: 10.1139/L03-045
32. 32. Jones H. Minimum cost prestressed concrete beam design. Journal of Structural Engineering. 1985;111:11. DOI: 10.1061/(ASCE)0733-9445(1985)111:11(2464)
33. 33. Lounis Z, Cohn MZ. Optimization of precast prestressed concrete bridge beam systems. Precast/Prestressed Concret Institute Journal. 1993;38(4):60-78. DOI: 10.15554/pcij.07011993.60.78

Written By

Alvaro Gaute-Alonso and David Garcia-Sanchez

Submitted: December 17th, 2021 Reviewed: December 22nd, 2021 Published: February 10th, 2022