## Abstract

### Purpose

Transfer relations represent analytical solutions of the linear theory of circular arches, relating each one of the kinematic and static variables at an arbitrary cross-section to the kinematic and static variables at the initial cross-section. The purpose of this paper is to demonstrate the significance of the transfer relations for structural analysis by means of three examples taken from civil engineering.

### Design/methodology/approach

The first example refers to an arch bridge, the second one to the vault of a metro station and the third one to a real-scale test of a segmental tunnel ring.

### Findings

The main conclusions drawn from these three examples are as follows: increasing the number of hangers/columns of the investigated arch bridge entails a reduction of the maximum bending moment of the arch, allowing it to approach, as much as possible, the desired thrust-line behavior; compared to the conventional *in situ* cast method, a combined precast and *in situ* cast method results in a decrease of the maximum bending moment of an element of the vault of the studied underground station by 46%; and the local behavior of the joints governs both the structural convergences and the bearing capacity of the tested segmental tunnel ring.

### Originality/value

The three examples underline that the transfer relations significantly facilitate computer-aided engineering of circular arch structures, including arch bridges, vaults of metro stations and segmental tunnel rings.

## Keywords

## Citation

Zhang, J.-L., Liu, X., Yuan, Y., Mang, H.A. and Pichler, B.L.A. (2021), "Transfer relations: useful basis for computer-aided engineering of circular arch structures", *Engineering Computations*, Vol. 38 No. 3, pp. 1287-1302. https://doi.org/10.1108/EC-06-2020-0313

## Publisher

:Emerald Publishing Limited

Copyright © 2020, Jiao-Long Zhang, Xian Liu, Yong Yuan, Herbert A. Mang and Bernhard L.A. Pichler.

## License

Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/legalcode

## 1. Introduction

The curvature of arches renders their structural analysis more challenging and expensive than that of straight beams. Transfer relations, representing analytical solutions of the linear theory of slender circular arches, were developed by Zhang *et al.* (2017). They are arranged, in a matrix-vector form, as follows: an 1 × 7 vector on the left-hand side, containing the kinematic and static variables at an arbitrary cross-section of the arch, is obtained as the product of a 7 × 7 matrix and an 1 × 7 vector, containing the kinematic and static variables at the initial cross-section of the arch. They serve as a powerful tool for structural analysis of arch structures.

The present paper is focused on structural analysis of three exemplary arches, based on the transfer relations.

Structural analysis of an arch bridge is the topic of the first example. In case of arches with low degrees of static indeterminacy and specific types of external loads, bridge designers may prefer analytical methods for structural analysis to numerical method. Concerning analytical methods, Fox (2000) reviewed the one of consistent deformations for pinned arches, which are statically indeterminate to the first degree. Melbourne (2008) presented a method based on Castigliano’s second theorem. Dym and Williams (2011) derived the governing differential equations of arches based on Sanders’ shell theory (Sanders, 1961) and variational calculus (Dym and Shames, 1973). Solutions of these equations for pinned or clamped arches, under radial and gravitational line loads, and for shallow arches, subjected to a concentrated axial load, acting at the end of the arch, were presented in Dym and Williams (2011) and Dym (2011). In most cases, arch bridges are complex in terms of boundary conditions and types of external loads, necessitating the use of numerical methods of structural analysis (Wei *et al.*, 2019). However, significant pre-processing efforts may be involved, particularly in case of sensitivity analysis of arch bridges, frequently requiring several changes of the geometric dimensions of the arch and the number and arrangement of the hangers or the columns (Hemp, 1974; Cheng, 2010). To improve this situation, transfer relations are considered in this work. Sensitivity analysis of the arch bridge “Zollamtssteg” in Vienna, with regard to the number of hangers, was carried out.

Structural analysis of the vault element of an underground station is the topic of the second example. The motivation for this example is the development of underground space in large cities. It refers to both the increasing number of underground infrastructures (Chen *et al.*, 2018) and the increasing scale of underground space (Sadaghiani and Dadizadeh, 2010). It has resulted in innovative underground structures and corresponding construction techniques, see e.g. Kivi *et al.* (2012), who developed a pre-supporting system consisting of arches and central columns, using it in the construction of a metro station in Tehran, Fang *et al.* (2012) who proposed a method for construction of metro stations by means of shallow tunneling, and Ding *et al.* (2019) who reported on a prefabricated large-span metro station in Changchun. The investigated underground station represents the first column-free metro station in Shanghai (Zhang *et al.*, 2019c). It consists of a large-span vault without columns. The vault was constructed by a combined precast and *in situ* cast method (Zhang *et al.*, 2020). It required subdivision of the vault into a number of identical elements. A typical vault element consisted of a precast part and an *in situ* cast part. Compared to the conventional *in situ* cast method, the period of construction, based on the combined precast and *in situ* cast method, was shortened significantly, Structural analysis of the vault element was carried out with the help of transfer relations. The discussion focused on the advantage of the combined precast and *in situ* cast method over the *in situ* cast method in the form of reducing the maximum bending moment of the vault element.

Structural analysis of a real-scale test of a segmental tunnel ring is the topic of the third example. Segmental tunnel rings typically consist of reinforced concrete segments and the joints connecting the segments (Ninić and Meschke, 2017). The structural behavior of the individual segmental tunnel rings is governed by the local behavior of the segment-to-segment joints (Zhang *et al.*, 2019b). Modeling of the joints is a non-trivial task, because their complex contact mechanical behavior depends on many factors (Liu *et al.*, 2017; Li *et al.*, 2015; Zhang *et al.*, 2019d). Among them are the geometric layout of the joints, the potential use of prestressed connecting bolts, the stresses acting in the contact zone and potential dislocations and relative rotation angles (Majdi *et al.*, 2016). To investigate the structural behavior of segmental tunnel rings, many real-scale tests have been carried out, see e.g. Luttikholt (2007) for railway tunnels, Liu *et al.* (2016) for metro tunnels, Qiu *et al.* (2019) for gas transmission tunnels and Huang *et al.* (2019) for water conveyance tunnels. As regards the present work, access to the recorded experimental data of the test by Liu *et al.* (2016) provided the motivation to develop a hybrid method for structural analysis of segmental tunnel rings. The measured displacement discontinuities at the joints were used, together with prescribed external loads, as input for hybrid structural analysis. It was based on transfer relations, allowing to consider these discontinuities in a straightforward manner. This eliminates the interfacial nonlinearities that may result from:

bending-induced partial segment-from-segment separation (Gladwell, 1980; Janßen, 1983);

the material behavior, including crushing of concrete and yielding of the steel bolts (Tvede-Jensen

*et al.*, 2017); andthe time-dependent viscoelastic material behavior of concrete (Schlappal

*et al.*, 2017, 2019).

Thus, the developed method significantly simplifies structural analysis.

The three described applications from bridge and subsurface engineering were selected to underline the versatility of the transfer relations. They allow for computing analytical solutions of the theory of thin circular arches for highly statically indeterminate arch bridges, for statically determined vaults representing three-hinged arches and for three-times kinematic segmental tunnel rings.

The paper is organized as follows. The transfer relations are presented in Section 2. Sections 3–5 are dedicated to the aforementioned three examples. Section 6 contains the main conclusions drawn from the present study.

## 2. Transfer relations

The transfer relations read as (Zhang *et al.*, 2017) follows:

*θ*; the radial and the tangential component of the displacement of the axis of the arch,

*u*and

*v*; the normal force

*N*; the shear force

*V*; and the bending moment

*M*. The matrix on the right-hand side is the so-called transfer matrix. Its top-left six-by-six submatrix refers to the solution for an unloaded part of the arch. The mathematical expressions for some of the nonzero elements of this submatrix are given in equation (2), where

*R*,

*EA*, and

*EI*denote the radius, the extensional stiffness and the bending stiffness of the arch, respectively. The first six elements in the last column of the transfer matrix in equation (1) refer to the superposition of analytical solutions for dead load, interfacial discontinuities of kinematic variables and point loads, given in Zhang

*et al.*(2017), as well as for a uniform temperature change, given in Zhang

*et al.*(2018), ground pressure, given in Zhang

*et al.*(2019b) and the overload on the ground surface, given in Zhang

*et al.*(2020). In equation (1), the vector on the right-hand side contains the kinematic and static variables at the initial cross-section of the arch (index“

*i*”), i.e. for φ = 0.

*u*,

_{i}*v*,

_{i}*θ*,

_{i}*M*,

_{i}*N*, and

_{i}*V*are obtained from the boundary conditions (Zhang

_{i}*et al.*, 2017). Once they are known, the kinematic and static variables at an arbitrary cross-section can be obtained in form of a simple matrix-vector product, see equation (1). These relations serve as the vehicle for structural analysis of the three previously mentioned arch structures. The algorithm which steers structural analysis of arch structures is listed in Table 1.

The displacement field **u** is evaluated with the help of the computed radial and tangential displacement components, *u* and *v*. **u** is given as follows:

**e**

*and*

_{r}*e*

_{φ}represent base vectors in the radial and the tangential direction, respectively, and

*r*denotes the radial coordinate. The normal stresses

*σ*and the shear stresses

*τ*are obtained from:

*N*and

*M*stand for the normal force and the bending moment, respectively;

*I*and

*S*denote the second and the first moment of area, respectively, and

*b*stands for the width of the cross-section.

## 3. Structural analysis of an arch bridge

The investigated arch bridge is the “Zollamtssteg” in Vienna [Figure 1(a)]. It consists of two parallel circular arches, resting on clamped supports, the bridge deck and two pairs of equally spaced vertical hangers/columns, connecting the deck and the arches [Figure 1(b)]. The number of each pair of hangers/columns is 21. The span *L* and the radius *R* of the arches are equal to 30.40 m and 16.73 m, respectively. The boundary conditions read as
*f*” represents the index of the final cross-section. The arches are made of steel with Young’s modulus amounting to 210 GPa. Their cross-sections are of “I” shape. The thickness of the flange and the web is equal to 2.6 mm and 4 mm, respectively. Both the width and the height of the cross-section amount to 50 mm. The extensional stiffness *EA*, the bending stiffness *EI*, and the dead load *q _{a}* of the arches are given as follows:

Herein, one half of the bridge, consisting of one arch, one half of the bridge deck, and one pair of hangers/columns, is investigated. The dead load, *q _{d}* = 0.0070 MN/m, refers to one half of the bridge deck. It is carried by the arch via the hangers, above the bridge deck, and the columns, below the bridge deck.

The sensitivity of the load-carrying behavior of the arch with respect to the number of equally spaced hangers/columns, *n*, was investigated, considering dead load of the arch and the bridge deck. Dead load of the hangers/columns was not taken into account, because it is much smaller than that of the arches and the bridge deck. The sensitivity analysis followed an approach similar to the one described in Zhang *et al.* (2018). At first, classical structural analysis of the deck, representing a continuous beam with *n* intermediate supports (Figure 2), was carried out. The forces of the hangers/columns were equal to the support forces. Finally, structural analysis of the arch, subjected to *n* point loads, was accomplished by means of the transfer relations (1). *n* ranges from 1 to 30. Hence, altogether 30 analyses were performed. The maximum bending moments of the arch,
*n* increases from 1 to 10, the maximum bending moment decreases by 60%. A further increase of *n* results in an insignificant decrease of the maximum bending moment of the arch. Therefore, with regards to the analyzed arch bridge, 10 equally spaced hangers/columns are sufficient for a good distribution of the deck weight.

Knowledge of the internal forces allows for stress analysis. With the help of equations (4) and (5), the von Mises stresses are obtained from:

The transfer relations are appealing to structural analysis of arch bridges, because of the following reasons:

Compared to mesh-based numerical simulation methods, the transfer relations reduce the amount of preprocessing in the framework of sensitivity analysis of arch bridges with respect to number of hangers/columns. The former typically require re-meshing of the arch when either changing the number of hangers or the locations of their connection to the arch, because each of such connections requires node-to-node connectivity. However, the transfer relations allow to account for the forces and locations of hangers/columns by means of load integrals, i.e. in a straightforward manner. This advantage is non-trivial, because preprocessing (including mesh generation) generally dominates the time needed for performing structural analysis.

Compared to discretization-based numerical solutions, the transfer relations, representing

*analytical*solutions, do not suffer from discretization errors. Therefore, there is no need for quantification of this error in form of a convergence study.Transfer relations are computationally inexpensive, because the solution of a given problem only involves evaluations of the product of a 7 × 7 transfer matrix and a 7 × 1 vector.

## 4. Structural analysis of a vault element of a metro station

The metro station considered is the Wuzhong Road Metro Station in Shanghai, see the details described in Zhang *et al.* (2020). Its vault, which is part of the station, is a segment of a cylindrical shell. It consists of 56 identical arch elements, arranged in a parallel fashion along the axis of the vault. The width *B*, the span *L* and the radius *R* of the elements are equal to 2.95 m, 19.8 m and 15.1 m, respectively. The vault was constructed by a combined precast and *in situ* cast method (Zhang *et al.*, 2020). Thus, a vault element consists of a precast part and an *in situ* cast part. The precast part serves as the permanent formwork of the *in situ* cast concrete, see Figure 5 for part of the permanent formwork. Each vault element consists of two precast arch segments. The shape of the cross-section of these segments is an upside-down *π* [Figure 6(a)]. The thickness of the flange amounts to 100 mm. The width and the height of the ribs are equal to 25 mm and 50 mm, respectively. The width *B* and the height *H* of the rectangular cross-section of the integral arch elements amounted to 2.95 m and 1.0 m, respectively. The strength class of the used concrete is C35. Its Young’s modulus is equal to 31.5 GPa.

The loads depend on the construction program. The precast segments were produced in a factory. After curing, they were transported to the construction site. Two segments each are assembled to a three-hinged arch, serving as a permanent formwork. The arch carries the dead load of the integral arch, *q* = 0.059 MN/m [Figure 6(a)]. After hardening of the *in situ* cast concrete, the integral arch becomes a clamped arch, subjected to equivalent ground pressures in the vertical and the horizontal directions, denoted as *σ _{v}* and

*σ*, respectively [Figure 6(b)]. These pressures are given as follows (Zhang

_{h}*et al.*, 2020):

*β*denotes the inclination angle of the initial cross-section relative to the

*x*-axis,

*h*stands for the height of the covering soil at the crown of the vault,

_{s}*p*

_{0}stands for the overload on the ground surface and

*γ*and

*K*denote the unit weight and the lateral pressure coefficient of the soil, respectively. As for the investigated vault element,

*β*,

*h*,

_{s}*p*

_{0},

*γ*and

*K*are equal to 0.1408 rad, 2.65 m, 20 kPa, 20 kN/m

^{3}and 0.65, respectively.

Two modes of structural analysis were carried out with the help of the transfer relations. The first one simulates the combined precast and *in situ* cast construction method. The results were obtained by superposition of results from structural analysis of the three-hinged arch, subjected to dead load, and results from structural analysis of the clamped arch, subjected to the ground pressure, see the magenta curves in Figure 7. The second mode simulates the conventional in-situ-cast construction method. The results were obtained from structural analysis of the clamped arch, subjected to both dead load and the ground pressure, see the green curves in Figure 7. The effect of the construction method on the deformed configuration and the distribution of the bending moments of the vault element is significant. This is underlined by the fact that the deflection at the crown of the vault element, constructed by the combined precast and *in situ* cast method, is by 308% larger than that obtained by the *in situ* cast method [Figure 7(a)], but that the maximum bending moment in case of the former is by 46% smaller than that in case of the latter, [Figure 7(b)]. However, the value of the deflection at the crown of the vault element, constructed by the combined precast and *in situ* cast method, is only 0.026% of the span of the vault. Hence, the percentage increase of this deflection compared to the one obtained by the *in situ* cast method is irrelevant. The effect of the construction method on the normal and the shear forces is insignificant [Figures 7(c) and 7(d)].

An advantage of using transfer relations for structural analysis of the vault element is the way of consideration of the boundary conditions. They are satisfied by means of determination of the six kinematic and static variables at the initial cross-section. The three boundary conditions at this cross-section directly provide the solutions for three of the six variables. The three boundary conditions at the final cross-section provide three linear equations for the remaining three kinematic and static variables at the initial cross-section. A change of the boundary conditions results only in a change of the six kinematic and static variables at the initial cross-section, whereas the transfer matrix remains unchanged. This is particularly useful for a comparative analysis of vault elements with different boundary conditions.

## 5. Hybrid structural analysis of a real-scale test of a segmental tunnel ring

The re-analyzed real-scale test of a segmental tunnel ring was carried out at Tongji University (Liu *et al.*, 2016; Figure 8). The radius *R* of the centerline of the tested ring was equal to 2.925 m. The rings consisted of six reinforced concrete segments. Young’s modulus *E*, the uniaxial compressive strength *f _{c}*, and the tensile strength

*f*of the concrete were determined from multiscale models (Königsberger

_{t}*et al.*, 2018; Hlobil, 2016), as 43.6 GPa, 62 MPa and 3.17 MPa, respectively. The width

*B*and the height

*H*of the cross-section of the segments amounted to 1.2 m and 35 cm, respectively. The ring was loaded all way up to its bearing capacity. This was done by means of three groups of altogether 24 equally spaced hydraulic jacks, acting on the outer surface of the ring. The experimental monitoring involved the jack forces and the displacement/rotation discontinuities, developing at the segment-to-segment joints, as well as the vertical and the horizontal convergences.

During the test, cracking of the reinforced concrete segments, a reduction of the initial contact area of the segment-to-segment joints, and crushing of concrete in compression at the joints was observed. Therefore, the analysis of the bearing capacity of the tested segmental tunnel ring requires consideration of the nonlinearities resulting from both cracking of the segments and the complex behavior of the joints, including segment-from-segment separation and the nonlinear behavior of the materials of the joints. Herein, a hybrid method (Zhang *et al.*, 2017) was used for structural analysis of the tested ring. This way, the complex behavior of the joints was considered in a straightforward manner, because the measured relative rotations at the joints, together with the external loading, were used as input. Bending-induced tensile cracking of the segments requires the extension of the linear transfer relations to consideration of the nonlinearities resulting from such cracking.

The hybrid analysis consists of two load cases (Zhang *et al.*, 2019a). This approach follows scientific work by Blom (2002) and by El Naggar and Hinchberger (2008). They assumed that the relative rotation angles at the joints result in rigid-body displacements of the segments. Load Case I refers to these rotations. They are estimated from monitored data, based on the Bernoulli–Euler hypothesis, and then post-processed such as to refer to rigid-body displacements of the segments, see Zhang *et al.* (2019a) for details. This post-processing is carried out with the help of transfer relations including the load integrals for the relative rotation angles at the joints. Based on the analysis of Load Case I, rigid-body displacements of the segments are obtained. Load Case II refers to the point loads, while the relative rotation angles are set equal to zero. In this case, the corresponding structural analysis is one of a closed ring without joints, subjected to point loads. To consider bending-induced tensile cracking of the six segments, they are subdivided into 34 elements. The length of the elements in the tangential direction is approximately equal to the distance of neighboring cracks. The elements contain a crack band in the middle and an undamaged zone on each one of the two sides. The extensional stiffness and the bending stiffness of the crack band is quantified by a multiscale model for tensile softening of concrete (Hlobil, 2016). Then, the effective bending and extensional stiffnesses of the damaged element are quantified, using the Voigt–Reuss–Hill estimate (Hill, 1952). These stiffnesses are constant within the elements. They are used for simulations of the segmental tunnel ring, based on the linear transfer relations. Superposition of the two load cases completes the nonlinear hybrid analysis of the segmental tunnel ring. This is admissible, even though the analysis of Load Case II is nonlinear, because Load Case I only results in rigid-body displacements and the equilibrium conditions of the structure are formulated in the undeformed configuration.

The model-predicted convergences are compared with corresponding experimental measurements (Figure 9). Their agreement underlines the usefulness of the chosen analysis approach. It is notable that rigid-body displacements of the segments (see the green curves in Figure 9) are responsible for approximately 95% of the convergences, whereas the deformations of the segments (see the black curves in Figure 9) are only responsible for the remaining 5%. Therefore, the relative rotation angles at the joints govern the convergences of the analyzed segmental tunnel ring.

The normal stresses of the concrete are evaluated from the obtained internal forces, the effective stiffnesses of the elements and the constitutive model of the concrete (Zhang *et al.*, 2019a). Then, the level of loading of the concrete is quantified by the stress levels *ζ*, defined as (Hellmich *et al.*, 2001):

*ω*denotes the damage variable, which is equal to the crack density parameter by Budiansky and O’Connell (Budiansky and O’Connell, 1976).

*ζ*= 1 indicates that concrete cracks in tension and

*ζ*= –1 that it crushes in compression. When the bearing capacity was reached,

*ζ*was equal to 1 at the inner surface in the crown and the bottom regions as well as at the outer surface in the two lateral regions, see Figure 10 for the distribution of the stress levels in the deformed configuration. However, the blue color, associated with

*ζ*= –1, does not appear. Because the development of plastic hinges requires the compressive stress of concrete to reach the compressive strength, i.e.

*ζ*= –1, there are no plastic hinges inside the segments. Further analysis of the joints has shown that the bearing capacity of the segmental tunnel ring is a consequence of the development of plastic hinges at the two top joints and the two lateral joints.

Herein, the relative rotation angles at the joints are post-processed such that they refer to a *symmetric* mode of rigid-body displacements of the segments. A generalized mode of rigid-body displacements was recently considered in Jiang *et al.* (2020). In addition, the described hybrid method of analysis of segmental tunnels requires measured relative rotations at the segment-to-segment interfaces. If such measured quantities are not available, it is necessary to resort to a classical mode of analysis of such tunnels. To this end, linear transfer relations are combined with a nonlinear interface law (Zhang *et al.*, 2019b).

The hybrid approach eliminates the need to consider the interfacial nonlinearities and, thus, significantly simplifies structural analysis (Zhang *et al.*, 2017). The transfer relations are very valuable for such hybrid structural analysis, because they consider, in a *straightforward* manner, the relative rotation angles at the joints. As for finite element methods, prescribing these relative rotations, however, is a challenging task.

## 6. Conclusions

Transfer relations, representing analytical solutions of the linear theory of circular arches, were used for computer-aided engineering analysis of an arch bridge, a vault of a metro station, and a segmental tunnel ring. The following conclusions are drawn from these three examples:

Increasing the number of hangers/columns of arch bridges results in less concentrated loading of the arch, rendering a reduction of the maximum bending moment of the arch. For the investigated arch bridge, ten equally spaced hangers/columns are sufficient for a good distribution of the deck weight. Further increase of this number results in an insignificant reduction of the bending moments.

As regards the vault element of the metro station at Wuzhong Road, Shanghai, a combined precast and

*in situ*cast method of construction results in a decrease of the maximum bending moment by 46%, compared to the one obtained by the conventional*in situ*cast method.The hybrid method allows for re-analysis of the bearing capacity of the tested segmental tunnel ring. The analysis results show that the local behavior of the joints governs both the structural convergences and the bearing capacity.

It is interesting to compare the used approach for structural analysis based on transfer relations to the alternative and very popular approach of numerical analysis based on the finite element method (FEM) for curved beams. Notably, results from FE simulations converge, with increasing discretization efforts (i.e. with increasing fineness of the used FE meshes), toward analytical solutions. Thus, every application of the FEM requires a convergence analysis, involving at least one coarser FE mesh and typically two consistently refined FE meshes. Thereby, mesh generation takes a significant amount of time, particularly so in the presently discussed context of linear structural analysis. This underlines the usefulness of the transfer relations, because they provide direct access to analytical solutions of the linear theory of thin circular arches. The resulting advantages, demonstrated by the three presented examples, are summarized as follows:

Transfer relations reduce the complexity of structural analysis of circular arches to that of structural systems consisting of straight beams. These relations are particularly attractive for sensitivity analysis regarding, e.g. the number of hangers/columns of arch bridges, because they can be easily specialized for different numbers of these structural elements.

Transfer relations are beneficial to a comparative analysis of vault elements with different boundary conditions, to which they can be easily specialized by means of the integration constants.

Transfer relations are very valuable for hybrid structural analysis of segmental tunnel rings, because they can consider, in a

*straightforward*manner, the relative rotation angles at the joints by means of load integrals.

## Figures

### Figure 8.

(a) Photo and (b) geometric dimensions of the analyzed segmental tunnel ring according to Liu *et al.* (2016) and Zhang *et al.* (2017), respectively

Flowchart illustrating the use of the transfer relations for structural analysis of arch structures

Input | Radius, extensional stiffness, and bending stiffness of the arch, R, EA and EI, and, if applicable, the values and the acting positions of external loads, e.g. point loads, interfacial discontinuities of kinematic variables, dead load, uniform temperature change, ground pressure, etc. |

(i) | Determine three of the integration constants (u, _{i}v, _{i}θ, _{i}M, _{i}N and _{i}V) from the boundary conditions at the initial cross-section_{i} |

(ii) | Determine the remaining three integration constants from the boundary conditions at the final cross-section |

(iii) | Evaluate the transfer relations (1) for any cross-sections of interest |

Output | Kinematic and static variables at any cross-section of interest |

## References

Blom, C.B.M. (2002), “Design philosophy of concrete linings for tunnels in soft soils”, Ph.D. thesis, Delft University of Technology.

Budiansky, B. and O’Connell, R.J. (1976), “Elastic moduli of a cracked solid”, International Journal of Solids and Structures, Vol. 12 No. 2, pp. 81-97.

Chen, Z.L., Chen, J.Y., Liu, H. and Zhang, Z.F. (2018), “Present status and development trends of underground space in Chinese cities: evaluation and analysis”, Tunnelling and Underground Space Technology, Vol. 71, pp. 253-270.

Cheng, J. (2010), “Optimum design of steel truss arch bridges using a hybrid genetic algorithm”, Journal of Constructional Steel Research, Vol. 66 Nos 8/9, pp. 1011-1017.

Ding, P., Tao, L., Yang, X., Zhao, J. and Shi, C. (2019), “Three-dimensional dynamic response analysis of a single-ring structure in a prefabricated subway station”, Sustainable Cities and Society, Vol. 45, pp. 271-286.

Dym, C.L. (2011), “End-loaded shallow curved beams”, Journal of Structural Engineering, Vol. 137 No. 7, pp. 782-784.

Dym, C.L. and Shames, I.H. (1973), Solid Mechanics: Variational Approach, McGraw-Hill, New York, NY.

Dym, C.L. and Williams, H.E. (2011), “Stress and displacement estimates for arches”, Journal of Structural Engineering, Vol. 137 No. 1, pp. 49-58.

El Naggar, H. and Hinchberger, S.D. (2008), “An analytical solution for jointed tunnel linings in elastic soil or rock”, Canadian Geotechnical Journal, Vol. 45 No. 11, pp. 1572-1593.

Fang, Q., Zhang, D. and Wong, L.N.Y. (2012), “Shallow tunnelling method (STM) for subway station construction in soft ground”, Tunnelling and Underground Space Technology, Vol. 29, pp. 10-30.

Fox, G.F. (2000), “Arch bridges”, in Chen, W.F. and Duan, L. (Eds), Bridge Engineering Handbook, 1st ed. Boca Raton, CRC Press, Chapter 17, pp. 17/1–17/12.

Gladwell, G.M. (1980), Contact Problems in the Classical Theory of Elasticity, Springer Science and Business Media.

Hellmich, C., Mang, H.A. and Ulm, F.J. (2001), “Hybrid method for quantification of stress states in shotcrete tunnel shells: combination of 3D in situ displacement measurements and thermochemoplastic material law”, Computers and Structures, Vol. 79, pp. 2103-2115.

Hemp, W. (1974), “Michell framework for uniform load between fixed supports”, Engineering Optimization, Vol. 1 No. 1, pp. 61-69.

Hill, R. (1952), “The elastic behaviour of a crystalline aggregate”, Proceedings of the Physical Society. Section A 65, pp. 349-354.

Hlobil, M. (2016), “Micromechanical analysis of blended cement-based composites”, Ph.D. thesis. Vienna University of Technology, Austria, and Czech Technical University in Prague, Czech Republic.

Huang, X., Liu, W., Zhang, Z., Wang, Q., Wang, S., Zhuang, Q., Zhu, Y. and Zhang, C. (2019), “Exploring the three-dimensional response of a water storage and sewage tunnel based on full-scale loading tests”, Tunnelling and Underground Space Technology, Vol. 88, pp. 156-168.

Janßen, P. (1983), “Tragverhalten von tunnelausbauten mit gelenktübbings [load-carrying behavior of segmented tunnel linings]”, Ph.D. thesis, Technische Universität Carolo-Wilhelmina zu Braunschweig, In German.

Jiang, Z., Liu, X., Schlappal, T., Zhang, J.L., Mang, H. and Pichler, B.L.A. (2020), “Asymmetric serviceability limit states of symmetrically loaded segmental tunnel rings: hybrid analysis of real-scale tests”, Tunnelling and Underground Space Technology, Under review.

Kivi, A.V., Sadaghiani, M. and Ahmadi, M. (2012), “Numerical modeling of ground settlement control of large span underground metro station in Tehran metro using central beam column (CBC) structure”, Tunnelling and Underground Space Technology, Vol. 28, pp. 1-9.

Königsberger, M., Hlobil, M., Delsaute, B., Staquet, S., Hellmich, C. and Pichler, B. (2018), “Hydrate failure in ITZ governs concrete strength: a micro-to-macro validated engineering mechanics model”, Cement and Concrete Research, Vol. 103, pp. 77-94.

Li, X., Yan, Z., Wang, Z. and Zhu, H. (2015), “A progressive model to simulate the full mechanical behavior of concrete segmental lining longitudinal joints”, Engineering Structures, Vol. 93, pp. 97-113.

Liu, X., Bai, Y., Yuan, Y. and Mang, H.A. (2016), “Experimental investigation of the ultimate bearing capacity of continuously jointed segmental tunnel linings”, Structure and Infrastructure Engineering, Vol. 12 No. 10, pp. 1364-1379.

Liu, X., Zhang, C., Zhang, C. and Yuan, Y. (2017), “Ultimate load-carrying capacity of the longitudinal joints in segmental tunnel linings: illustrated with the typical segment joints of the Shanghai metro rapid transit system”, Structural Concrete, Vol. 18 No. 5, pp. 693-709.

Luttikholt, A. (2007), “Ultimate limit state analysis of a segmented tunnel lining”, Master’s thesis, Delft University of Technology.

Majdi, A., Ajamzadeh, H. and Nadimi, S. (2016), “Investigation of moment-rotation relation in different joint types and evaluation of their effects on segmental tunnel lining”, Arabian Journal of Geosciences, Vol. 9 No. 7, pp. 1-15.

Melbourne, C. (2008), “Design of arch bridges”, in Parke, G. and Hewson, N. (Eds), ICE Manual of Bridge Engineering, 2nd ed. Thomas Telford, London, pp. 305-344.

Ninić, J. and Meschke, G. (2017), “Simulation based evaluation of time-variant loadings acting on tunnel linings during mechanized tunnel construction”, Engineering Structures, Vol. 135, pp. 21-40.

Qiu, Y., Feng, K., He, C., Zhang, L. and Wang, C. (2019), “Investigation of the ultimate bearing capacity of a staggered assembly segmental lining for an urban gas transmission tunnel”, Sustainable Cities and Society, Vol. 48, p. 101551.

Sadaghiani, M.H. and Dadizadeh, S. (2010), “Study on the effect of a new construction method for a large span metro underground station in Tabriz-Iran”, Tunnelling and Underground Space Technology, Vol. 25 No. 1, pp. 63-69.

Sanders, J.L. (1961), “Nonlinear theories for thin shells”, Quarterly of Applied Mathematics, Vol. 21, p. 35.

Schlappal, T., Schweigler, M., Gmainer, S., Peyerl, M. and Pichler, B. (2017), “Creep and cracking of concrete hinges: insight from centric and eccentric compression experiments”, Materials and Structures, Vol. 50 No. 6, p. 244.

Schlappal, T., Kalliauer, J., Vill, M., Gmainer, S., Mang, H.A., Eberhardsteiner, J. and Pichler, B.L. (2019), “Serviceability limits of reinforced concrete hinges”,* Engineering structures*, p. 109861.

Tvede-Jensen, B., Faurschou, M. and Kasper, T. (2017), “A modelling approach for joint rotations of segmental concrete tunnel linings”, Tunnelling and Underground Space Technology, Vol. 67, pp. 61-67.

Wei, J., Guan, M., Cao, Q. and Wang, R. (2019), “Spatial combined cable element for cable-supported bridges”, Engineering Computations, Vol. 36 No. 1, pp. 204-225.

Zhang, J.L., Liu, X., Zhao, J. and Mang, H. (2020), “Application of a combined precast and in-situ-cast construction method for large-span underground vaults”, Tunnelling and Underground Space Technology, Under revision.

Zhang, J.L., Hellmich, C., Mang, H.A., Yuan, Y. and Pichler, B. (2018), “Application of transfer relations to structural analysis of arch bridges”, Computer Assisted Methods in Engineering and Science, Vol. 24, pp. 199-215.

Zhang, J.L., Mang, H.A., Liu, X., Yuan, Y. and Pichler, B. (2019a), “On a nonlinear hybrid method for multiscale analysis of a bearing-capacity test of a real-scale segmental tunnel ring”, International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 43 No. 7, pp. 1343-1372.

Zhang, J.L., Schlappal, T., Yuan, Y., Mang, H.A. and Pichler, B. (2019b), “The influence of interfacial joints on the structural behavior of segmental tunnel rings subjected to ground pressure”, Tunnelling and Underground Space Technology, Vol. 84, pp. 538-556.

Zhang, L., Feng, K., Li, M., Sun, Y., He, C. and Xiao, M. (2019d), “Analytical method regarding compression-bending capacity of segmental joints: theoretical model and verification”, Tunnelling and Underground Space Technology, Vol. 93, p. 103083.

Zhang, J.L., Vida, C., Yuan, Y., Hellmich, C., Mang, H.A. and Pichler, B. (2017), “A hybrid analysis method for displacement-monitored segmented circular tunnel rings”, Engineering Structures, Vol. 148, pp. 839-856.

Zhang, J.L., Zhao, J.B., Liu, X., Jiang, S.C., Yang, Z.H. and Jiang, H.C. (2019c), “Experimental investigation on structural behavior of a stacked reinforced-concrete arch”, in Derkowski, W., Gwoździewicz, P., Hojdys, Ł., Krajewski, P. and Pańtak, M. (Eds), Proceedings of the Fib Symposium 2019: Concrete – Innovations in Materials, Design and Structures, International Federation for Structural Concrete, Kraków, Poland, pp. 2270-2277.

## Acknowledgements

Financial support by the Austrian Science Fund (FWF), provided for the project P 281 31-N32 “Bridging the Gap by Means of Multiscale Structural Analyses”, is gratefully acknowledged. In addition, the first author is indebted to the National Natural Science Foundation of China (Grant No. 51908424), the Shanghai Pujiang Program (Grant No. 19PJ1409700) and to the Austrian Federal Ministry of Education, Science and Research (BMBWF) (Grant No. ICM-2019-14045) for financial support of this work. The second author is also indebted to the National Natural Science Foundation of China (Grant No. 51578409), for financial support of this work.