Abstract
This article presents the modeling of a complete bolted connection based on a model with one bolt connecting two or three plates. Initially, the behavior of this model with one bolt is analyzed by comparing it with existing bibliography for 3 different types of applied load: tension, shear and a combination of these two. This model includes all necessary considerations: contacts between the plates and the nut, head and shank of the bolt; contact between the plates, as well as friction between them; and pre-load on the bolt. The model also responds properly to loads parallel and perpendicular to the contact surface between the plates, including prying action effects. These calibrated models are then introduced as super-elements in empty spaces left on the full connection, through a relatively simple process using the finite element software ANSYS®. Upon filling these spaces, two complete connection models are evaluated: one with a single plate and one between two T-stubs. The results obtained with these models are compared with standard forecasts. These two connection types have a practical application in the way they will be analyzed, and also as part of more complex connections: bolted girder splices (in the region of web beam), beam-to-column connections with splice plates or end plate connections, beam splices with end-plates, etc.
Keywords
Modeling of bolted connections; Simple shear connections; T-stub connections; Super-elements
1 INTRODUCTION
During the last years, in the discussions involving structural analysis and global performance of steel structures, the semi-rigid behavior of the connections has become a topic of extreme importance, and was introduced in the calculation procedures of AISC (1980)American Institute of Steel Construction (AISC), Steel Construction Manual, 8th Edition, American Institute of Steel Construction, Chicago, IL, 1980. and, later, in the design methodology proposed by Eurocode 3 (1993)Eurocode 3 (1993). Design steel of structures: Part 1.1 - General rules and rules for buildings - Revised Annex J: Joints in building frames..
The semi-rigid behavior of the connections has been incorporated into the practice of structural analysis, not only because it better represents the actual behavior of the structure relative to its overall response, but also because it contributes to the design of more economical structural elements (Maggi et al., 2005Maggi, Y. I., Gonçalves, R. M., Leon, R. T., Ribeiro, L. F. L. (2005). Parametric analysis of steel bolted end plate connections using finite element modelling, Journal of Constructional Steel Research 61: 689-708.).
Among the connecting means employed in the practice, the bolted connections distinguish themselves for the simplicity and rapidity in the assembly process of steel structures. However, the large number of variables in this type of connection - the thickness of the plates involved in the connection, the diameter and positioning of the fasteners, etc. - make it difficult to analyze their behavior (Queiroz et al., 2013Queiroz, G., Carvalho, H., Vilela, P. M. L., Garcia, L. B. L. (2013). Sobre as regiões nodais das estruturas de aço formadas por barras, Revista da Estrutura de Aço, v.2, No. 3, p. 206-217.).
There are numerous references in literature dealing with bolted connections through numerical analysis, and most of them used one-dimensional and two-dimensional elements in the numerical models development. With the advancements of finite element programs and with a greater processing capacity of the computers, a refinement in the models is noticed, and these bolted connections started to be simulated with the use of solid elements, making it possible to consider the nodal region stiffness more realistically.
Yilmaz and Bekiroglu (2014)Yilmaz, O. e Bekiroglu, S. (2014). Numerical Simulations of Bolted Steel Connections in Single and Double Shear under Pretension Effect. In: 11th International Congress on Advances in Civil Engineering. Istanbul, Turkey., for example, simulated bolts subjected to single and double shear loadings, through ANSYS® software, evaluating the initial pretension influence on the behavior of these connections. The models were developed with solid elements, the efforts between the parts were transmitted using contact elements, however the thread geometry of the bolts was not considered. Moreover, Fransplass et al. (2015)Fransplass, H., Langseth, M., Hopperstad, O. S. (2015). Experimental and numerical study of threaded steel fasteners under combined tension and shear at elevated loading rates, International Journal of Impact Engineering 76: 118-125. investigated the behavior of threaded steel fasteners subjected to combined tension and shear loading at elevated loading rates, also using three-dimensional finite element model, but in this case, the thread geometry was simulated with high levels of detail. This numerical simulation captured the overall physical behavior of the threaded assembly, yet the values of ultimate load and initial elastic stiffness for all the load cases investigated were higher than the experimental results obtained in Fransplass et al. (2015)Fransplass, H., Langseth, M., Hopperstad, O. S. (2015). Experimental and numerical study of threaded steel fasteners under combined tension and shear at elevated loading rates, International Journal of Impact Engineering 76: 118-125. research.
Since nodal region is composed of several elements with different behaviors, and the simulation of this region, as well as the interaction between its components, makes bolted connections modeling complex, a finite element model with one bolt connecting two or three square plates was developed in this work, in such a way that all the necessary considerations were included: contact between the plates and the nut, the head and the shank of the bolt; contact and friction between the plates; pretension in the bolt. This unitary model, that adequately responded to forces acting parallel and perpendicular to the contact surface between the plates, capturing the effects of prying action, was finally inserted in some complete connections, with different amount and position of the bolts, which shows the versatility of this form of numerical modelling.
2 METHODOLOGY
The methodology adopted for this study consists in:
-
validating a finite element model with only one bolt through comparisons with experimental results available in the literature, considering bolts subjected to tension, shear and a combination of tension and shear;
-
inserting this model with one bolt in a complete connection, modeled with finite elements too, comparing the numerical results of these connections with actual standards.
Two types of connections were chosen: with the bolt group loaded eccentrically and subjected, essentially, to shear; and a T-stub type connection, where the bolts are subjected to tension. All models were developed using the APDL language of software ANSYS®.
3 UNITARY NUMERICAL MODEL
3.1 Numerical model description
The element type used on the modelling of the components is SOLID186. The plates of the unitary model are square and present a central hole with diameter that corresponds to the diameter of the standard hole indicated by the AISC (2010)American Institute of Steel Construction (AISC), ANSI/AISC 360-10, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2010. standard, bigger than the nominal diameter of the bolt. The geometry of the bolt - Figure 1 - presents some simplifications: the nut and the head of the bolt have the same diameter, D, and the same thickness, t, with values based on American specifications for high-strength bolts, ASTM A325 (ASTM, 1985)American Society for Testing and Materials (ASTM), High-Strength Bolts for Structural and Steel Joints, ASTM Designation A325M-84a, Philadelphia, 1985.; the thread region, whose effective area is equivalent to, approximately, 75% of the shank area, was modelled with the nominal diameter of the bolt, however with a less resistant material.
Equations 1-4 show how to obtain the mentioned diameter, D, and thickness, t:
where d nut is the considered diameter of the nut, calculated using the average of the diameters inscribed and circumscribed to the hexagon, represented, respectively, by W and W/cos 30° - see Figure 2;
dhead is the considered diameter of the bolt head, calculated using the average of the diameters inscribed and circumscribed to the hexagon, represented, respectively, by F and F/cos 30° - see Figure 2,
H1, H2 is the height of the bolt head and of the nut, respectively.
The trilinear stress-strain relationship used to define the material property of the plates and of the thread and shank portions of the bolt is shown in Figure 3. The first line segment represents the elastic range of the material, with slope equivalent to the Young’s modulus, E, equal to 200 GPa (ABNT, 2008Associação Brasileira de Normas Técnicas (ABNT), NBR 8800:2008, Projeto de estrutura de aço e de estrutura mista de aço e concreto de edifícios, Rio de Janeiro, 2008.), and yield stress equal to f y. The second and third line segments represent the plastic range, where it is assumed that the tensile strength, f u, is reached when the deformation values exceed 2%. The ultimate deformation, ε u, corresponds to the deformation attained at the rupture of the material. Yield stress and tensile strength for the thread material are equal to these values for the shank material of the bolt multiplied by 0.75 (Vilela, 2016Vilela, P. M. L. (2016). Modelagem de ligações parafusadas pelo método dos elementos finitos. Dissertation (in Portuguese), Federal University of Minas Gerais, Brazil.).
For an adequate simulation of the contact between the parts, surface-to-surface contact elements TARGE170 and CONTA174 were used. In the determination of the contact pairs, the ANSYS® default values were adopted for most of the parameters, except the FKN value for the contact pair established between the hole and the bolt, adjusted to 0.01 - recommended value for regions with bending deformations (ANSYS - 17.0, 2016ANSYS - 17.0. Release 1.0 Documentation for Ansys. Canonsburg, United States, 2016.). The tangent stiffness factor, FKT, is calculated automatically by the software and is proportional to the normal stiffness and to the friction coefficient, μ, assumed to be 0.3 - based on AISC (2010).
In the definition of the contact elements CONTA174, the KEYOPT(4) and KEYOPT(10) parameters were considered with a value of 2, which mean, respectively, that the contact detection occurs at the nodes of the contact elements (and not at the Gauss points), and that the contact stiffness is updated at each iteration (not at each load step).
It is important to check, through the command CNCHECK,DETAIL, the status of all contact pairs present in the model before starting the simulation. The parts that should have been in contact since the beginning of simulation must have “closed” contact status, which does not occur only for the region between the bolt and the plates hole, that has “open” contact status due to the gap between these elements.
The bolt preload was applied through the PRETS179 element, like in Wang et al. (2013)Wang, Y. Q., Zong, L., Shi, Y. J. (2013). Bending behavior and design model of bolted flange-plate connection, Journal of Constructional Steel Research 84: 1-16..
According to Salamet and Garlock (2010)Salamet, S. and Garlock, M. (2010). Guidelines for Modeling Three Dimensional Structural Connection Models using Finite Element Methods. International Symposium “Steel Structures: Culture & Sustainability”. Istanbul, Turkey, September 21-23., for regions with high stress concentrations, such as on bolts and plate holes, mesh convergence studies suggest 20 to 24 elements around the circumference of a typical bolt diameter (22.225 mm to 31.75 mm). Figure 4 presents the mesh of the unitary model components, with 24 elements around the circumference that define the diameter of the plate holes and bolt shank. The mesh configuration of the rest of the model is based on this value, which resulted in a unitary model with 3240 elements for the plates discretization and 3888 elements for the discretization of the bolt.
Figure 5 shows an example of a full connection with empty spaces, where the unitary model is inserted. The mesh of the matching edges must be the same so that the parts can be joined together.
3.2 Results
Several studies have shown that there are many variables governing the behavior of threaded steel fasteners subjected to tension, shear or combined tension and shear loading conditions (Fransplass et al., 2015Fransplass, H., Langseth, M., Hopperstad, O. S. (2015). Experimental and numerical study of threaded steel fasteners under combined tension and shear at elevated loading rates, International Journal of Impact Engineering 76: 118-125.).
Rumpf and Fisher (1963)Rumpf, J. L. and Fisher, J. W. (1963). Calibration of A325 bolts, Lehigh University, Fritz Engineering Laboratory, Report No. 288.5, December 1963. studied the behavior and performance of individual A325 high-strength bolts under direct tension and torqued tension, Wallaert and Fisher (1964)Wallaert, J. J. and Fisher, J. W. (1964). Shear Strength of High-Strength Bolts, Lehigh University, Fritz Engineering Laboratory, Report No. 288.20, July 1964. determined the behavior of single high-strength bolts under static shear loadings. Chesson et al. (1964)Chesson, E., Faustino, N. L., Munse, W. H. (1964). Static Strength of High-Strength Bolts under Combined Shear and Tension, University of Illinois, March 1964. aimed at defining the strength and behavior characteristics of single, high-strength bolts under static loadings of tension and shear. All these works were carried out as part of the Large Bolted Joints Project, under the technical guidance of the Research Council on Riveted and Bolted Structural Joints, and the unitary model was validated according to the experimental results presented by these three works, considering bolts subjected to tension, shear and a combination of these efforts. The main aspects of the unitary model simulations (number of analysis, bolt diameters, plate thicknesses, etc.) are presented in Table 1.
Figure 6 and Table 2 presents numerical and experimental results of two bolt series, B (ASTM A325) and C (ASTM A354), for each T-S ratio evaluated in Chesson et al. (1964)Chesson, E., Faustino, N. L., Munse, W. H. (1964). Static Strength of High-Strength Bolts under Combined Shear and Tension, University of Illinois, March 1964.. T and S are, respectively, the tension and the shear forces applied on the bolts while T u is their ultimate strength in tension. The experimental interaction curve proposed by the authors is also represented in Figure 6 (green curve).
In the last column of Table 2, it can be observed that the numerical results (blue points in Figure 6) showed good agreement with the experimental results (red points in Figure 6) and the maximum difference between them is of 5%. Further numerical results for the unitary model, which also showed good agreement with the experimental ones, are not presented in this work. More details can be found in Vilela (2016)Vilela, P. M. L. (2016). Modelagem de ligações parafusadas pelo método dos elementos finitos. Dissertation (in Portuguese), Federal University of Minas Gerais, Brazil..
4 COMPLETE BOLTED CONNECTIONS MODELS - ECCENTRICALLY LOADED CONNECTIONS
4.1 General
When the line of action of an applied load does not pass through the center of gravity of a bolt group, the moment resulting from this eccentricity must be considered in the connection design. The eccentricity produces both rotation about the centroid of the bolt group and translation of one connected element with respect to the other. The combined effect of this rotation and translation is equivalent to a rotation about a point, defined as instantaneous center of rotation, and its location depends on the geometry of the bolt group as well as the direction and point of the load application (AISC, 2005American Institute of Steel Construction (AISC), ANSI/AISC 360-05, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2005.).
The AISC Steel Construction Manual presents some tables that employ the instantaneous center of rotation method for different bolt patterns and eccentric conditions. The available strength of the bolt group for any of the arrangement of the tables is defined by ϕR n, which is equivalent to the coefficient C multiplied by ϕr n, the available strength of a single bolt.
4.2 Numerical model description
This work analyzes some configurations of eccentrically loaded bolt groups provided in the AISC manual, varying the values of the horizontal component of the load eccentricity, e x, number of bolts in one vertical row, n, and bolt spacing, s. All these variables are represented in Figure 7, that illustrates three of the complete eccentrically loaded models that were simulated: the first has 3 bolts, spaced vertically by 76.2 mm; the second, 3 bolts spaced by 152.4 mm; and the third, 5 bolts spaced by 76.2 mm. A typical mesh employed in these numerical models is also displayed in Figure 7.
On the evaluated numerical models, a thickness of 25.4 mm was adopted for plate 1 and a thickness of 12.7 mm was used for plate 2, which length varies according to the value of e x. The friction coefficient between the plates was considered equal to 0.3. All bolts are of type ASTM A325, have a diameter of 19.05 mm, thread length of 3.18 mm to guarantee shear plane out of the thread region, and an initial preload equal to, approximately, 70% of the minimum specified tensile strength of the bolt. The material of the plates is ASTM A36. Table 3 describes the values adopted for the definition of the trilinear curves of these materials.
Values of the trilinear curve stress-strain for the materials of the complete model with eccentrically loaded bolt group.
4.3 Results
Table 4 presents a comparison between the values of the coefficients C recommended by the AISC manual and the coefficients obtained through the numerical simulations for the indicated values of n, e x and s. This coefficient was determined dividing the resistant shear force of the complete connection (last converged step) by the resistant shear force of the unitary model.
A good agreement was observed between the numerical results and the AISC (2005)American Institute of Steel Construction (AISC), ANSI/AISC 360-05, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2005. recommended values. Plates with different thicknesses were not evaluated, once numerically obtained coefficients C were calculated based on the ultimate bolt resistance, that is, considering the attendance of the other ultimate limit-states applicable to the plates. Figure 8 presents a distribution of von Mises stresses for the last converged step of the models with n equal to 3 and 5, s equal to 76.2 mm and e x equal to 152.4 mm. When a group of bolts is eccentrically loaded, the bolts of the extremities are the most requested. Such behavior can be seen in Figure 8, where the external bolts failed, while the internal bolts did not develop all of their strength.
Distribution of von Mises stresses on eccentrically loaded bolt group models: a) n=3; b) n=5.
Figure 9, Figure 10 and Figure 11 present the load-displacement curves for all analyzed cases, on which the ultimate strengths of the numerical models is defined. The node selected to plot all the graphs is localized at the top right corner of the plate 2, and its total displacement value is the vector sum of the displacements along x-axis and y-axis (there is no displacement out of the plane, along z-axis).
5 COMPLETE BOLTED CONNECTIONS MODELS - T-STUB CONNECTIONS
5.1 General
The T-stub connection is one of the simplest connections where the bolts are subjected to tension. Depending on the geometric characteristic and the material employed on the connection, additional forces near the flange extremities can appear. This effect, known as prying action, increases the effort in the bolt, and should be taken into account in the analysis of the connection (Latour et al., 2014Latour, M., Rizzano, G., Santiago, A., da Silva, L. S. (2014). Experimental analysis and mechanical modeling of T-stubs with four bolts per row, Journal of Constructional Steel Research 101: 158-174.).
The American code provides a design procedure that includes the prying action, where bolt diameters and connected flange thicknesses are determined in such a way that bolt resistance, as well as flange stiffness and strength, are above the required.
Figure 12 presents the variables involved on consideration of prying action. The force per bolt can be determined by T+q, where T represents the tension force per bolt due to the external acting force, 2T, and q is the additional tension force per bolt due to prying action.
When the connection geometry is known, the available tensile strength including the effects of prying action, T avail, can be obtained multiplying the available tensile strength per bolt, B, by Q:
where Q is the factor that represents the presence or not of sufficient strength and stiffness of the flanges to develop all the available tensile strength of the bolts. The expressions used to calculate factor Q, showed in Table 5, depends on the value of á’, a parameter that either maximizes the bolt available tensile strength for a given thickness or minimizes the thickness required for a given bolt available tensile strength (AISC, 2005).
5.2 Numerical model description
ASTM A325 bolts, with 19.05 mm diameter, and ASTM A572 Gr50 tee shapes were considered for evaluation of the behavior of T-stub connection subjected to tension. Since bolts of these connections are essentially subjected to tension, the thread material was used to define the shank and thread of these model bolts and the stress-strain values are presented in Table 6.
The variables in Figure 12 have the following values: distance between the holes, g, equal to 120 mm; distance from the bolt centerline to the face tee stem, a, equal to 40 mm; tributary length per pair of bolts, perpendicular to the page, p, equal to 120 mm. The other variable values are shown in Table 7, as well as the tee flange and web thicknesses of the evaluated cases. Figure 13 illustrates one of the complete simulated T-stub model and a typical mesh employed in this type of connection.
5.3 Results
Knowing the geometry of the connection, it is possible to calculate the required flange thickness, t c, to develop the available tensile strength of the bolt, B, without the effect of prying action. Besides these values, Table 8 shows the available strength of the evaluated connections, T avail, considering the effect of prying action, if it occurs.
For flange thickness equal to 12.5 mm, the last step considered in the numerical model occurred for a load equal to 294 kN, when the stress level at the interface between the web and the flange reached the material strength of the plates, indicating connection failure. Figure 14a presents the von Mises stresses for the last step of the mentioned tee, where it is possible to observe an accentuated bend of the plates and, also, the presence of extension regions with stress above steel yield stress at the centerline of the holes.
For flange thickness equal to 25.0 mm, the last step considered in the numerical model occurred for a load, T numerical, equal to 646 kN, close to the value indicated by the formulation of AISC (2005)American Institute of Steel Construction (AISC), ANSI/AISC 360-05, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2005.. Figure 14b shows the von Mises stresses for this tee shape, where it is possible to observe a zone of intense plastification next to the web. Like in the previous analyzed model, the material strength of the plates is reached at the interface between the web and the flange.
For flange thickness equal to 37.5 mm, the last converged step occurred for a load, T numerical, equal to 702 kN, which is equivalent to the available strength of the 4 bolts group. The high stiffness of flange tee shapes avoids the occurrence of prying action effects, what is in accordance with the AISC (2005)American Institute of Steel Construction (AISC), ANSI/AISC 360-05, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2005. standard. Figure 15 exhibits the von Mises stresses diagram for the connection ultimate load, indicating the presence of few points with stress above the yield limit.
Von Mises stresses for complete T-stub model with flange thickness equal to: a) 12.5 mm; b) 25.0 mm
Figure 16 presents the load-displacement curves for all analyzed cases. The node selected to plot all the graphs is localized at the intersection of the web with the flange of T-stub, and its total displacement value is the displacement along z-axis.
6 CONCLUSIONS
This work developed a numerical methodology to evaluate the behavior of one bolt subjected to tension, shear and a combination of these efforts, using software ANSYS®. This unitary model was inserted, subsequently, in two types of complete connections, and the results were compared with standard predictions.
In the connections with an eccentrically loaded bolt group, it was possible to determine the coefficients C for different geometric configurations and all of them presented good agreement with the values suggested by AISC (2005). The maximum difference between the results was of 2.3%.
In the numerical T-stub connections, distinguished failure modes were observed for different evaluated flange thicknesses. For flange thicknesses smaller than 19.0 mm, the failure was characterized by the rupture of the plates due to bending. For flange thickness equal to 25.0 mm, the failure occurred simultaneously at the flange (interception with the web) and at the bolts. When flange thickness is equal to 37.5 mm, there were no prying action effects, occurring bolt ruptures due to direct acting of the external load through the tee web. The available strength of all the evaluated models differ up to 9% in relation with that determined by AISC (2005).
Then, the proposed numerical methodology proves to be interesting, since it enables the evaluation of any bolted connection, including connections not provided by the reference standards.
References
- American Institute of Steel Construction (AISC), Steel Construction Manual, 8th Edition, American Institute of Steel Construction, Chicago, IL, 1980.
- American Institute of Steel Construction (AISC), ANSI/AISC 360-05, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2005.
- American Institute of Steel Construction (AISC), ANSI/AISC 360-10, Specification for Structural Steel Buildings, American Institute of Steel Construction, Chicago, IL, 2010.
- American Society for Testing and Materials (ASTM), High-Strength Bolts for Structural and Steel Joints, ASTM Designation A325M-84a, Philadelphia, 1985.
- ANSYS - 17.0. Release 1.0 Documentation for Ansys. Canonsburg, United States, 2016.
- Associação Brasileira de Normas Técnicas (ABNT), NBR 8800:2008, Projeto de estrutura de aço e de estrutura mista de aço e concreto de edifícios, Rio de Janeiro, 2008.
- Chesson, E., Faustino, N. L., Munse, W. H. (1964). Static Strength of High-Strength Bolts under Combined Shear and Tension, University of Illinois, March 1964.
- Eurocode 3 (1993). Design steel of structures: Part 1.1 - General rules and rules for buildings - Revised Annex J: Joints in building frames.
- Fransplass, H., Langseth, M., Hopperstad, O. S. (2015). Experimental and numerical study of threaded steel fasteners under combined tension and shear at elevated loading rates, International Journal of Impact Engineering 76: 118-125.
- Latour, M., Rizzano, G., Santiago, A., da Silva, L. S. (2014). Experimental analysis and mechanical modeling of T-stubs with four bolts per row, Journal of Constructional Steel Research 101: 158-174.
- Maggi, Y. I., Gonçalves, R. M., Leon, R. T., Ribeiro, L. F. L. (2005). Parametric analysis of steel bolted end plate connections using finite element modelling, Journal of Constructional Steel Research 61: 689-708.
- Queiroz, G., Carvalho, H., Vilela, P. M. L., Garcia, L. B. L. (2013). Sobre as regiões nodais das estruturas de aço formadas por barras, Revista da Estrutura de Aço, v.2, No. 3, p. 206-217.
- Rumpf, J. L. and Fisher, J. W. (1963). Calibration of A325 bolts, Lehigh University, Fritz Engineering Laboratory, Report No. 288.5, December 1963.
- Salamet, S. and Garlock, M. (2010). Guidelines for Modeling Three Dimensional Structural Connection Models using Finite Element Methods. International Symposium “Steel Structures: Culture & Sustainability”. Istanbul, Turkey, September 21-23.
- Vilela, P. M. L. (2016). Modelagem de ligações parafusadas pelo método dos elementos finitos. Dissertation (in Portuguese), Federal University of Minas Gerais, Brazil.
- Wallaert, J. J. and Fisher, J. W. (1964). Shear Strength of High-Strength Bolts, Lehigh University, Fritz Engineering Laboratory, Report No. 288.20, July 1964.
- Wang, Y. Q., Zong, L., Shi, Y. J. (2013). Bending behavior and design model of bolted flange-plate connection, Journal of Constructional Steel Research 84: 1-16.
- Yilmaz, O. e Bekiroglu, S. (2014). Numerical Simulations of Bolted Steel Connections in Single and Double Shear under Pretension Effect. In: 11th International Congress on Advances in Civil Engineering. Istanbul, Turkey.
-
Available Online: May 16, 2018
Publication Dates
-
Publication in this collection
2018
History
-
Received
31 July 2017 -
Reviewed
21 Apr 2018 -
Accepted
14 May 2018