Next Article in Journal
Luminescence and Electrochemical Activity of New Unsymmetrical 3-Imino-1,8-naphthalimide Derivatives
Next Article in Special Issue
A Nonlinear Elastic Model for Compressible Aluminum Alloys with Finite Element Implementation
Previous Article in Journal
Halloysite in Different Ceramic Products: A Review
Previous Article in Special Issue
Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables

Department of Civil and Industrial Engineering, University of Pisa, Largo Lazzarino, 1-56122 Pisa, Italy
Materials 2021, 14(19), 5502; https://doi.org/10.3390/ma14195502
Submission received: 18 August 2021 / Revised: 17 September 2021 / Accepted: 21 September 2021 / Published: 23 September 2021

Abstract

:
Cables are typically used in engineering applications as tensile members. Relevant examples are the main cables of suspension bridges, the stays of cable-stayed bridges, the load-bearing and stabilizing cables of tensile structures, the anchor cables of floating mooring structures, the guy-ropes for ship masts, towers, and wind turbines, the copper cables of electrical power lines. Since cables are characterized by non-linear behavior, analysis of cable structures often requires advanced techniques, like non-linear FEM, able to consider geometric non-linearity. Nevertheless, a traditional simplified approach consists in replacing the cable with an equivalent tie rod, characterized by a suitable non-linear constitutive law. Currently used equivalent constitutive laws have been derived by Dischinger, Ernst and Irvine. Since the equivalence is restricted to taut cables, characterized by small sag to chord ratios, these traditional formulae are not appropriate for uniformly loaded sagging cables: the main cables of suspension bridges are a particularly emblematic case. Despite some recent attempts to find more refined solutions, the problem is still open, since closed form solutions of general validity are not available. In the paper, general analytical formulae of the non-linear constitutive law of the equivalent tie rod are proposed, distinguishing two relevant cases, according as the length of the cable can vary or not. The expressions, derived by applying the general form of the theorem of virtual work, can be applied independently on the material, on the sag to chord ratio, on the load intensity and on the stress level, so allowing the replacement of the whole cable with a single equivalent tie rod. The expressions are critically discussed referring to a wide parametric study also in comparison with the existing formulae, stressing the influence of the most relevant parameters.

1. Introduction

Cables are widely used in engineering fields as typical load bearing tensile members: relevant examples are the main cables of suspension bridges, the stays of cable-stayed bridges, the load-bearing and stabilizing cables of tensile structures, the anchor cables of floating mooring structures, the guy-ropes for masts, ship masts, towers, and wind turbines, the copper cables of electrical power lines, and so on. Nowadays, depending on the application, cables can be made resorting to a large variety of different materials. Suitable materials span from the historical ones, like natural fibers (cotton, flax, jute, silk), the use of which is lost in mists of time, to the modern ones, like high-strength steel, up to the most advanced ones, like aramid, glass, polyester, and carbon fiber. Innovative structural materials are more and more proposed: for suspension and cable-stayed bridges [1,2], for other general structural applications [3,4], as well as for strengthening, restoration and repair of historical and heritage buildings [5], even considering hybrid solutions [6].
Since cables are characterized by geometric non-linear behavior, analysis of cable structures requires suitable advanced techniques, like non-linear FEM, able to consider this kind of non-linearity, while the material is assumed linear elastic. A traditional and very effective approach to simplify the structural model consists in replacing the cable, whose behavior is governed by geometric nonlinearity, with an equivalent tie rod connecting the ends of the cable, characterized by a suitable non-linear constitutive law. The rationale of this approach is emulating the apparent along the chord stiffness of the cable with by the axial stiffness of the equivalent tie rod: it is crystal clear that in this way, the complexity of structural models can be strongly reduced, as it occurs in finite element analyses, where geometric non-linearities can be captured by representing the cable with a unique equivalent element characterized by material non-linearity.
Under the basic assumptions that bending and shear stiffnesses of the cable are negligible in comparison with its axial stiffness, and that cables’ configurations are characterized by small sag to chord ratios, the assessment of the non-linear constitutive law of the equivalent tie rod has been the subject of several studies [7,8,9,10], which provide suitable approximate expressions for the tangent or the secant elastic modulus of the equivalent straight tie. Widely used classical expressions have been derived by Dischinger [7,8] and Ernst [9], considering a horizontal cable. Said a the chord length, A 0 the area of the cross section, and p the unit weight of the cable, acting in the vertical direction, Dischinger [7,8] expressed the equivalent elastic tangent modulus, E t , e q , as ratio between the variation, d σ 0 , of the horizontal component of the normal stress in the cable, σ 0 , and the variation of the deformation, d ε , along the chord:
E t , e q σ 0 = d σ 0 d ε = E 1 + p   a 2 12   A 0 2   σ 0 3   E ,
being E is the elastic modulus of the cable material. Of course, from the equilibrium equations it derives that σ 0 is independent on the considered cross section. Considering a variation of the normal stress between σ 0 i and σ 0 f , Ernst [9], following the Dischinger’s approach, derived the equivalent secant modulus as:
E s , e q = Δ σ 0 Δ ε = E 1 + p   a 2 24   A 0 2   σ 0 i 3   1 + σ ¯ σ ¯ 2 E ,
where
σ ¯ = σ 0 f σ 0 i .
As these studies are based on simplified assumptions, the equivalence is strictly restricted to stretched cables subjected to the self-weight, so characterized by catenary shape, even if Equations (1) and (2) have been often adopted also for uniformly loaded cables, in case the sag to chord ratio is small enough to approximate the actual deformed shape, which is a parabola, with a catenary. In fact, let be the uniformly distributed load q acting perpendicularly to the chord A B of the cable, and assume the axes, x and y , of a 2D Cartesian coordinate system, whose origin is in A , given by A B , and q , respectively, the configuration y x of the cable under q is given by the parabolic funicular curve (Figure 1).
y x = q   2   N 0 x a x = 4   f   a x 1 x a ,
where f is the sag of the cable and N 0 ,
N 0 = q   a 2   8   f ,
is the x -component of the normal force N in the cable, which is obtained from the equilibrium condition of the portion A C of the cable (Figure 2). For the sake of simplicity, it can be assumed that x is horizontal, and y is vertical.
In principle, Equation (4) is derived assuming the cable inextensible, i.e., having infinite axial rigidity, but it can be properly modified to take into account the Hooke’s law for the cable [10], so making unnecessary the inextensibility hypothesis. Equation (4) thus applies independently on the cable material properties the material is elastic, and the initial length of the cable, L , is greater than the chord, a : L > a , provided it is written referring to the final, stretched, configuration. The condition L > a ensures that if the external load is nil, the cable is unstressed. In case L < a , since the installation of the cable is conditional on suitable pre-strain, ε ¯ p , of the cable itself
ε ¯ p = a L 1 ,
elastic deformations need must be explicitly considered.
To improve the equivalent model, widening its field of application, Irvine [10] suggested to express the tangent elastic modulus in the form
E t , e q 1 = d σ d ε = E 1 + λ 2 12 ,
where the deformed configuration of the cable is considered by means of the so-called characteristic parameter of the cable, λ 2
λ 2 = q   a N 0 2 a   E σ 0   L e ,
where L e is the virtual length of the cable, depending on the arc length generalized coordinate s ,
L e = 0 a d s d x 3 d x .
Recalling Equation (4), after some elementary manipulations, Equation (8) gives
L e = 0 a 1 + y 2 3 2 d x = a 8   1 + 4   f a 2   5 + 2   4   f a 2 + 3 4   a f   asinh 4   f a ,
which, in case a f , can be satisfactorily approximated, following Irvine, by
L e a   1 + 8   f a 2 .
For the sake of completeness, we recall that actual length of the cable L is
L = 0 a 1 + y 2 1 2 d x = a 2   1 + 4   f a 2 + 1 4   a f   asinh 4   f a ,
which can be estimated, when a > 0.1   f , by
L a   1 + 8 3   f a 2 ,
When applicable, Irvine expression, Equation (7), is generally more reliable than the Dischinger formula, Equation (1); anyhow, not fully considering the influence of the modifications of the cable configurations on the position of the external loads, these classical expressions are both dependent on important limitations, which limit their field of application to taut cables. It must be remarked that relevant structures can be characterized by uniformly loaded sagging cables, where the cited approaches are not appropriate: the main cables of suspension bridges are particularly emblematic cases. For this reason, several theoretical and numerical studies have been carried out [11,12,13,14,15,16,17,18,19,20,21], with the aim of achieving more general solutions, valid also for relaxed cables. Moreover, in the last years, further improvements of the classical expressions have been proposed, in view of the implementation of appropriate non-linear cable elements in finite element analysis [22,23,24,25], also discussing the influence of bending stiffness and shear stiffness [26,27]. Anyhow, when significant, the effects of bending stiffness and shear stiffness on the cable configuration can be inferred as local perturbations of the solution derived disregarding them.
It must be remarked that uniformly loaded sagging cables characterize relevant structures, where the cited approaches are not appropriate: the main cables of suspension bridges are particularly emblematic cases. Despite the already cited relevant literature [22,23,24,25] illustrating some attempts to find more refined solutions, the problem is still open, since closed form solutions of general validity are not yet available.
In the paper, a general analytical expression of the non-linear constitutive law of an equivalent tie rod simulating uniformly loaded cables is proposed, distinguishing two relevant cases, according to whether the length of the cable can vary or not:
  • the first case, which is rather new, corresponds to the case when an extremity of the cable is fixed and the final section on the other end can move over a fixed pulley, so that the chord length cannot vary;
  • the second case, which is the one traditionally considered, corresponds to the case when the cable length is fixed.
The expressions, derived by applying the general form of the theorem of virtual work, can be applied to uniformly loaded cables, independently on the material, on the sag to chord ratio, and on the load intensity, so allowing the replacement of the whole cable with a single equivalent tie rod. The theoretical approach is a general implementation of a procedure already successfully proposed to derive the equivalent along the chord stiffness of horizontal [28] and inclined stays [29], loaded by the self-weight.
It must be remarked that the proposed approach can be also a powerful tool to obtain reliable prior information, in view of the multifaceted implementation of advanced Bayesian updating techniques for the identification of mechanical parameters [30,31,32,33,34].
The obtained expression is thus critically discussed in some relevant case studies, also in comparison with the existing formulae, as well as with the results of non-linear finite element analysis, stressing the influence of the most relevant parameters as well as its advantages.

2. The Non-Linear Constitutive Laws of the Equivalent Tie-Rod

The non-Linear constitutive law of the equivalent tie-rod can be derived by applying the theorem of virtual works to the cable, considering a horizontal virtual relative displacement between the cable ends, d x [28,29]. Without loss of generality, we assume the end A is fixed, so that the relative displacement is at the end B , d x B .
As anticipated, two relevant cases can be envisaged, depending on the boundary conditions in B , according as whether the distance between A and B , i.e., the chord length a , is fixed or not:
(a)
if the chord length is fixed, and the cable runs on a fixed pulley in B , the sole effect of the virtual displacement d x B is a variation of the cable length L : d L = d x B , and d a = 0 (Figure 3);
(b)
if cable ends, A and B , are fixed, the effect of the relative displacement is a variation of the cable chord a : d a = d x B (Figure 4).
The former case typically describes the behavior of cables during the tightening phase, the latter the behavior of structural cables in usual operational conditions.

2.1. Theorem of Virtual Work for the Cable

Referring to both previously mentioned cases, the virtual work equation for the cable can be expressed in the form,
N B + d N B · d x B + 0 a q · d y   d x = V σ d ε   d V ,
where d y is the variation of the cable ordinate due the modification of the cable configuration, σ is the normal stress, d ε the variation of the longitudinal strain, and V the volume. Evidently, if the cable is assumed inextensible, d ε = 0 , and the virtual work of normal stresses is nil.
Since q and y are parallel and d N B · d x B can be disregarded in comparison with N B · d x B , Equation (14) becomes
N B · d x B + 0 a q   d y d x = V σ d ε   d V .
The virtual displacement d x B in B causes a variation of the cable configuration, d y , which is the effect of the variations of the chord length, d a , as well as of the normal force, d N . The effect of the previously occurring elastic deformations is implicitly considered, setting the length of the “inextensible” cable equal to its final, deformed, length.
The total variation of the ordinate d y of a point of the cable, whose abscissa is x , depends on the variation of the chord length d a , and on the variation of the horizontal component of the normal force d N 0 . From another perspective, the total variation of the ordinate d y can be seen as sum of two contributions, the one associated with the variation of the configuration of the cable, assumed inextensible, d y i n , the other associated with the elastic stretching of the cable, d y e l . It must be remarked that, while the variation of the horizontal component of the normal force, d N 0 , influences both the inelastic contribution, d y i n , and the elastic contribution, d y e l , the variation of the chord length, d a , affects only d y i n , being negligible its effect on d y e l .
The contribution due to the modification of the inextensible cable configuration, can be obtained differentiating the cable equation (Equation (4)):
d y i n = y a d a + y N 0 d N 0 ,
which reduces to
d y i n = y N 0 d N 0 ,
in the previously cited case (a), when the chord length is fixed: d a = 0 .
The variation of the cable ordinate, d y e l , due to the elastic deformation is a function of the vertical component of the normal force only. A very good estimate of the contribution of the variation of the elastic deformations can be obtained approximating the parabola with a catenary, whose total weight is equal to total applied load. The unit self-weight of the equivalent catenary q * is thus
q * = p a L c ,
being L c is the length of the catenary,
L c = 2   N 0 p *   sinh q *   a 2   N 0 ,
from which,
q * = 2   N 0 a   arcsinh q   a 2   N 0 ,
which is independent on the Lagrangian coordinate s .
As suggested by Irvine [10], y e l s can be expressed as
Δ y e l s = N 0 E   A 0 y 0 s q * s 2 2   E   A 0 .
Differentiating Equation (21) with respect to N 0 and disregarding the effects of the normal force variation on p q * , we find
d y e l = y 0   E   A 0 + N 0 E   A 0 d y 0 d N 0 s   d N 0 .
Integrating on the cable length, we obtain
0 a q y e l N 0 d x 0 L c q * y e l N 0 d s   d N 0 = q * y 0   E   A 0 + N 0 d N 0 d y 0 E   A 0 0 L c s   d s = q *   L c 2 2   E   A 0   sinh q *   a 2   N 0 q *   a 2   N 0   cosh q *   a 2   N 0 = q *   L c 2 2   E   A 0   q   a 2   N 0 1 + q   a 2   N 0 2   arcsinh q   a 2   N 0 = q 2   a 3 4   N 0   E   A 0   q   a 2   N 0 arcsinh q   a 2   N 0 1 + q   a 2   N 0 2  
Obviously, the present study duly considers effects of elastic deformations: in fact, in the following, only the variation of the integration limits associated with the elastic elongation of the cable are disregarded, according to common practice.

2.2. The Non-Linear Constitutive Law for a Cable Running on a Fixed Pulley

In case the cable runs on a fixed pulley at its end B (Figure 3), it results in d a = 0 . Recalling Equations (4), (5), and (17), Equation (15) becomes:
N B   d l + 0 a q d y e l   d x + 0 a q y N 0 d N 0     d x = 0 L N A 0 d N E   A 0 A 0   d s ,
which reduces, after some elementary passages, to
N 0   1 + y a 2   d x + 0 a q y e l N 0 d N 0     d x + 0 a q y N 0 d N 0     d x = 0 L N 0 cos θ E   A 0   d N 0 cos θ   d s , .
where θ is the angle between the deformed configuration and the x -axis (Figure 2).
Further manipulations of Equation (25) lead to:
N 0   1 + y a 2   d x + 0 a q y e l N 0 d N 0     d x + 0 a q y N 0 d N 0     d x = N 0   d N 0 E   A 0 0 a   1 + y 2 3 2   d x .
Evaluating separately the individual contributions of each relevant term of Equation (26), we obtain
N 0   1 + y a 2   d x = N 0   1 + q   a 2   N 0 2   d a ,
0 a q y N 0 d N 0   d x = 0 a q 2 x 2   N 0 2 x a   d x   d N 0 = q 2 a 3 12   N 0 2   d N 0 ,
and
0 a   1 + y 2 3 2   d x =   1 8   a   1 + q   a 2   N 0 2   5 + 2 q   a 2   N 0 2 + 6   N 0 q   asinh q   a 2   N 0 ,
which obviously corresponds to Equation (10), and finally, recalling Equation (29):
1 + q   a 2   N 0 2   d x = q 2 a 3 12   N 0 3   d N 0 +   d N 0 8   E   A 0   { a   1 + q   a 2   N 0 2   5 + 2 q   a 2   N 0 2 +   6 N 0 q   asinh q   a 2   N 0 2   q 2 a 3 N 0 2 q   a 2   N 0 asinh q   a 2   N 0 1 + q   a 2   N 0 2   } .
By dividing both members of Equation (30) by a , and recalling that
d ε = d x a ,
is the apparent along the chord deformation of the cable, and of the equivalent tie-rod, it results
1 + ϱ   a 2   σ 0 2   d ε = ϱ 2 a 2 12   σ 0 3   d σ 0 +   d σ 0 8   E     {   1 + ϱ   a 2   σ 0 2   5 + 2 ϱ   a 2   σ 0 2 + 6   σ 0 ρ   a   asinh ϱ   a 2   σ 0 2   ϱ 2 a 2 σ 0 2 ϱ   a 2   σ 0 asinh ϱ   a 2   σ 0 1 + ϱ   a 2   σ 0 2   } ,
where ϱ is the specific load, i.e., the ratio between the uniformly distributed load and the area of the cable
ϱ = q A 0 .
The non-linear constitutive law of the equivalent tie-rod, expressed in terms of tangent elastic modulus, is thus:
E t , e q , 0 = d σ 0 d ε =   1 + ϱ   a 2   σ 0 2   E 1 3 ϱ   a 2   σ 0 2 E σ 0 + 1 8   1 + ϱ   a 2   σ 0 2   5 + 2 ϱ   a 2   σ 0 2 + 3 4   σ 0 ρ   a   asinh ϱ   a 2   σ 0 ϱ   a 2   σ 0 2 ϱ   a 2   σ 0 2   σ 0   asinh ϱ   a 2   σ 0 1 + ϱ   a 2   σ 0 2   ,
which, introducing the non-dimensional parameter ξ ¯ ,
ξ ¯ = ϱ   a 2   σ 0 ,
can be written as
E t , e q , 0 =   1 + ξ ¯ 2   E ξ ¯ 2 3   E σ 0 + 1 8   1 + ξ ¯ 2   5 + 2   ξ ¯ 2 + 3 8   ξ ¯     asinh ξ ¯ ξ ¯ 2 ξ ¯   asinh ξ ¯ 1 + ξ ¯ 2   ,
or, equivalently, stressing the dependency on the cable sag f
E t , e q , 0 = 1 + 4   f a 2   E 1 3 4   f a 2 E σ 0 + 1 8 1 + 4   f a 2   5 + 2 4   f a 2 + 3 8   a 4   f   asinh 4   f a 4   f a 2 4   f a   asinh 4   f a 1 + 4   f a 2   .

2.3. The Non-Linear Constitutive Law for a Cable with Fixed Ends

In case both ends of the cable are fixed (Figure 4), i.e., the inextensible length of the cable is not varying, the virtual displacement d x entirely results in a variation of the chord length: d a = d x . In this case, recalling Equation (16) and, again, Equations (4), (5) and (15) becomes
N 0   d a + 0 a q d y e l   d x + 0 a q y N 0 d N 0 + y a d a d x = 0 L N A 0 d N E   A 0 A 0   d s ,
from which one obtains the analogous of Equation (25)
N 0     d a + 0 a q y e l N 0 d N 0     d x + 0 a q y N 0 d N 0 + y a d a d x = 0 L N 0 cos θ E   A 0   d N 0 cos θ   d s .
Remembering the previously calculated contributions, Equations (23) and (28), and given that it is
0 a q y a d a   d x = q   a 2   N 0 2 d a ,
it results:
1 + q   a 2   N 0 2   d a = q 2 a 3 12   N 0 3   d N 0 +   d N 0 8   E   A 0   { a   1 + q   a 2   N 0 2   5 + 2 q   a 2   N 0 2 + 6   N 0 q   asinh q   a 2   N 0 2   q 2 a 3 N 0 2 q   a 2   N 0 asinh q   a 2   N 0 1 + q   a 2   N 0 2   } .
Dividing both members of Equation (41) by a , the variation of the chord deformation
d ε = d a a ,
can be derived
1 + ϱ   a 2   σ 0 2   d ε = ϱ 2 a 2 12   σ 0 3   d σ 0 +   d σ 0 8   E     {   1 + ϱ   a 2   σ 0 2   5 + 2 ϱ   a 2   σ 0 2 + 6   σ 0 ρ   a   asinh ϱ   a 2   σ 0 2   ϱ 2 a 2 σ 0 2 ϱ   a 2   σ 0 asinh ϱ   a 2   σ 0 1 + ϱ   a 2   σ 0 2   } ,
where ϱ is the specific load (Equation (33)). In this case, the tangent elastic modulus expressing the non-linear constitutive law of the equivalent tie-rod is
E t , e q , a = d σ 0 d ε =   1 + ϱ   a 2   σ 0 2   E 1 3 ϱ   a 2   σ 0 2 E σ 0 + 1 8   1 + ϱ   a 2   σ 0 2   5 + 2 ϱ   a 2   σ 0 2 + 3 4   σ 0 ρ   a   asinh ϱ   a 2   σ 0 ϱ   a 2   σ 0 2 ϱ   a 2   σ 0 2   σ 0   asinh ϱ   a 2   σ 0 1 + ϱ   a 2   σ 0 2   ,
or, equivalently,
E t , e q , a =   1 + ξ ¯ 2   E ξ ¯ 2 3   E σ 0 + 1 8   1 + ξ ¯ 2   5 + 2   ξ ¯ 2 + 3 8   ξ ¯     asinh ξ ¯ ξ ¯ 2 ξ ¯   asinh ξ ¯ 1 + ξ ¯ 2   ,
or, as a function of the cable sag f ,
E t , e q , a = 1 + 4   f a 2   E 1 3 4   f a 2 E σ 0 + 1 8 1 + 4   f a 2   5 + 2 4   f a 2 + 3 8   a 4   f   a s i n h 4   f a 4   f a 2 4   f a   a s i n h 4   f a 1 + 4   f a 2   .

2.4. Preliminary Remarks

Looking at the expressions of the equivalent elastic moduli, Equations (34) and (44), or Equations (37) and (46), it can be remarked that their limits exactly match the expectations. In fact, considering that the equilibrium conditions imply that, for q 0 ,
lim     f a   0   σ 0 = ,
and
lim       f a   σ 0 = 0 ,
the limits of E t , e q , 0 and E t , e q , a are
lim     f a   0   E t , e q , 0 = E ; lim     f a   0   E t , e q , a = E ,   and
lim       f a   E t , e q , 0 = 0 ;   lim       f a   E t , e q , a = 0 .
In addition, on a given set of conditions, the rigidity of the cable with fixed ends is higher than the rigidity of the cable with one end running on a fixed pulley. This aspect was already noticed for cables loaded by self-weight in [28,29].
It must be also underlined that the contribution to the equivalent axial stiffness due to the elastic increment of the cable ordinate, d y e l , is significant only for high values of the sag to chord ratio, f / a > 0.1 .
A particularly relevant feature of uniformly loaded cables is that the constitutive law of the equivalent tie-rod is a function of σ 0 / ϱ , and E / σ 0 : as the upper limit of σ 0 is controlled by the tensile strength, f t , of the cable material, and σ 0 depends on ϱ via Equation (5), the performances of the cable thus pivot on the apparent specific strength, f t / ρ ¯ , and on the apparent specific stiffness, E / ρ ¯ , of the cable material, being ρ ¯
ρ ¯ = ϱ g ,  
a sort of equivalent density, obtained regarding the specific load as an equivalent specific weight. Expressed in this way, governing parameters are analogous to the specific strength, f t / ρ , and to the specific stiffness, E / ρ , respectively, characterizing the behavior of a cable, made with a material whose density is ρ , transversally loaded by the self-weight only, as extensively discussed in [28]. Anyhow, it must be stressed that a substantial difference exists between the two cases: in fact, while density ρ is a material property, the equivalent density ρ ¯ is not; consequently, cables made by different materials can exhibit similar behaviors, provided their apparent specific strength, f t / ρ ¯ , and their apparent specific stiffness, E / ρ ¯ , are the same.
The expressions are valid also for creep sensitive materials, for which long term phenomena are relevant, by replacing the Young modulus of the material with the appropriate age-adjusted effective elastic modulus, taking into account the time dependency by means of a suitable creep coefficient [35].

3. Parametric Studies

The formulae previously derived allow to directly obtain the tangent elastic modulus of the equivalent tie-rod once the mechanical properties of the cable material, the cable geometry, and the applied load are known. However, in the author’s opinion, the most effective graphical representation of the ratio between the equivalent elastic modulus and the elastic modulus of the cable material is as a function of the sag to chord ratio, f / a , and of the normal stress to elastic modulus ratio, σ 0 / E . For this reason, in the following, the constitutive laws of the equivalent tie-rod are plotted as a function of the previously cited relevant parameters.
From the operational point of view, the sag to chord ratio is commonly in the range 0–0.125, but in the following we explore the range 0–1.0, independently of its practical feasibility. Regarding the ratio σ 0 / E , some additional consideration is necessary to fix its upper limit. Although currently adopted materials for structural cables are steel, aramid fiber, and carbon fiber [28,36], for the sake of the present study, we can explore a much more general set of natural and artificial materials, covering in principle every material suitable to build ropes and cables, as summarized in Table 1. In the Table, density, elastic modulus, ultimate strength, and ratios between the ultimate strength and the elastic modulus are given for relevant fibers. The mechanical properties of the cable materials reported in Table 1 have been derived summarizing data coming from several sources: the interested reader can specifically refer not only to the relevant literature [30,37], but also to some web sources [38,39], which provide links to cable material producers, too.
Looking at the Table, it clearly results that for most materials f t / E < 1 % ; for polybenzoxazole (PBO), and carbon fibers is f t / E = 1.5 2 % ; for aramid fibers, and liquid crystal aromatic polyester (LCP) fibers f t / E = 2.5 4.5 % ; while nylon and polyester fibers exhibit higher values. Moreover, copper and aluminum, which are very important in terms of applications, are characterized by very small values of f t / E , being f t / E = 0.15 % , and f t / E = 0.22 % , respectively. As the design limit of σ 0 is a fraction of the ultimate strength of the material, in the following, we consider that σ 0 / E varies in the interval 0 ,   0.02 .

3.1. The Equivalent Elastic Modulus of Cables Running on Fixed Pulleys

Considering the case of cable running on a fixed pulley, the ratio, E t , e q , 0 / E , between the tangent elastic modulus of the equivalent tie-rod, E t , e q , 0 , and the elastic modulus of the cable material, derived from Equations (34), (36), and (37), is diagrammatically represented in the 3D-graph of Figure 5. To facilitate the examination of the diagram, it is plotted in terms of contour lines in Figure 6, Figure 7 and Figure 8. More precisely, the E t , e q , 0 / E f / a diagrams in Figure 6 and Figure 7 illustrate the contour lines obtained for σ 0 / E values belonging to the intervals 0.01 % ,   0.1 % , and 0.2 % ,   2.0 % , respectively, while the E t , e q , 0 / E σ 0 / E diagrams in Figure 8a,b show the contour lines for sag to chord ratios, f / a , varying in the range 0.0025 ,   0.1 , and 0.125 ,   1.0 , respectively.

3.2. The Equivalent Elastic Modulus of Cables with Fixed Ends

Adopting criteria similar to those illustrated in Section 3.1, the tangent elastic modulus of the equivalent tie-rod, E t , e q , a , can be calculated also in case the ends of the cable are fixed. The ratio, E t , e q , a / E , is represented in the 3D-graph of Figure 9.
Again, to facilitate the examination, the diagram is plotted in terms of contour lines in Figure 10, Figure 11 and Figure 12. The E t , e q , a / E f / a diagrams in Figure 10 and Figure 11 illustrate the contour lines obtained for σ 0 / E values belonging to the intervals 0.01 % ,   0.1 % , and 0.2 % ,   2.0 % , respectively, while the E t , e q , a / E σ 0 / E diagrams in Figure 12a,b show the contour lines for sag to chord ratios, f / a , varying in the range 0.0025 ,   0.1 , and 0.125 ,   1.0 , respectively.

4. Discussion

The examination of the diagrams demonstrates that, as expected, the ratio between the tangent modulus of the equivalent tie-rod and the elastic modulus of the material, E t , e q / E :
  • decreases as the sag to chord ratio, f / a , increases and that the reduction rate is particularly marked in the region where f / a and σ 0 / E , are low;
  • it is a quasi-linear function of the stress when f / a 0.08 ;
  • rises as the load, and then the ratio σ 0 / E , increases; again, the increment rate is particularly marked in the region where f / a and σ 0 / E , are low;
  • for this reason, for creep sensitive materials, when the age-adjusted moduli E t , e q t , and E t are taken into account, it results:
    E t , e q t E   t > E t , e q E ,
    since the effect of the creep is an increase of the ratio σ 0 / E t , as better illustrated in the diagrams of Figure 13 and Figure 14, where the E t , e q , 0 / E f / a diagrams, and the E t , e q , a / E f / a diagrams, respectively, are parameterized for high value of σ 0 / E , in the interval 4.0 % ,   20 % . Moreover, these diagrams, in combination with the ones previously given in Section 3, allow to consider all possible practical cases, including consideration of long terms effects.
It must be stressed that the influence of the boundary conditions, or, in other words, of the the end restraints, is relevant only for f / a 0.08 , and becomes very significant for flabby cables: in fact, the ratio E t , e q , a / E t , e q , 0 is:
E t , e q , a E t , e q , 0 = 1 + 4 f a 2 .
Considering how they have been derived, the proposed expressions allow one to assess the constitutive laws of the equivalent tie-rod for uniformly loaded cables, whatever the stress and the sag. Recalling that, said γ the specific weight of the cable material and ξ ,
ξ = γ   a 2   σ 0 ,
the non-dimensional parameter of the catenary, analogous to ξ ¯ (see Equation (29)), the constitutive laws for the equivalent tie-rod of horizontal cables loaded by the self-weight [28] are given by:
E t , e q , 0 E = cosh ξ 1 12   ξ   9   sin h ξ + sinh 3   ξ 24 sinh 3 ξ + 12 ξ   sinh ξ sinh 2   ξ + E 8   σ 0 cos h   2   ξ 2 1 ξ sinh 2   ξ ,
for the cable running on a fixed pulley at one end, and by
E t , e q , a E = cos h 2 ξ 1 12   ξ   9   sin h ξ + sinh 3   ξ 24 sinh 3 ξ + 12 ξ   sinh ξ sinh 2   ξ + E 8   σ 0 cos h   2   ξ 2 1 ξ sinh 2   ξ ,
for the cable with fixed ends, it is thus possible to assess parabolic or catenary cables, whichever the sag to chord ratio and the stress level.

4.1. Comparison with “Historical” Formulae

To further enrich the discussion, in the present Subsection the results obtained with the proposed formulae are compared with those obtained using the most used “historical” formulae, namely Dischinger’s formula [1,2], Equation (1), and Irvine’s formula [10], Equation (7), deriving the horizontal component of the tensile force on the cable from Equation (5), and assuming q p . Regarding Irvine’s formula, they are considered both expressions of the virtual length of the cable, L e : the simplified version, Equation (11), and the accurate version, Equation (10). In the following, for the sake of completeness, also the case of cable with one end running on a fixed pulley is considered, even if it is not contemplated by the traditional formulae.
Although these historical expressions are normally used also for catenary cables, we will discuss here only uniformly loaded cables. Further information about catenary cables is given in [28,29]; moreover, the assessment of the equivalent uniformly distributed load and of the equivalent parabolic shape of sagging cables transversally loaded by the self-weight only is not an easy task, involving subtle equilibrium considerations.

4.1.1. Comparison with the Dischinger’s Formula Results

The ratios between the tangent elastic modulus of the equivalent tie-rod, evaluated by means of the Dischinger’s formula Equation (1), E t , e q , D , and by using the expressions derived before, E t , e q , 0 (Equation (34)), and E t , e q , a (Equation (44)) are illustrated in Figure 15 and Figure 16, respectively, as a function of the sag to chord ratio f / a , parameterized with respect to four relevant values of the σ 0 / E ratio, namely, 0.1%, 2%, 10%, and 20%.
The diagrams clearly show that Dischinger’s formula provides acceptable results only when the sag to chord ratio is very small and the stress is relatively high. In any case, Dischinger’s formula systematically underestimates the correct value.

4.1.2. Comparison with the Irvine’s Formula Results

The comparisons with Irvine’s formula results are summarized in Figure 17 and Figure 18, assuming L e expressed by the approximate Equation (11), and in Figure 19 and Figure 20, assuming L e expressed by the accurate Equation (10), analogously to what just done for Dischinger’s formula in Section 4.1.1. In the diagrams, E t , e q , I r v s indicates the equivalent tangent modulus obtained using simplified Equation (11), and E t , e q , I r v e the equivalent tangent modulus obtained using Equation (10).
Analyzing the graphs, in the common case of cables with fixed ends,
(a)
when the virtual length L e is calculated by means of the simplified expression, Equation (11), Irvine’s formula always underestimates the equivalent elastic modulus, leading to satisfactory results only for very small values of the sag to chord ratio f / a : in this field, it is more accurate than the Dischinger’s formula, being the ratio E t , e q , I r v s / E t , e q , a practically insensitive to the stress level;
(b)
when the virtual length L e is calculated by means of the accurate expression, Equation (10), Irvine’s formula can be adopted when the f / a is small. In its field of application, it generally leads to more precise evaluations of the equivalent modulus, both in comparison with Dischinger’s formula and with the simplified approach, anyhow underestimating the correct value, and resulting as much more sensitive to the stress level. In addition, as soon as f / a increases, the equivalent modulus is sensibly overestimated.
In the special case of a cable with the end running on a pulley, equivalent values calculated using the Irvine’s formula are always overestimated, especially when the accurate formula is adopted for L e .

4.2. Comparison with Non-Linear FEA

Of course, for a sounder validation of the proposed approach, the analytical results should be analyzed in the light of relevant experimental studies results, especially concerning the context of the cables characterized by high values of the sag to chord ratio, where the proposed method provides estimates much more accurate than the traditional methods. Unfortunately, in the author’s knowledge, there are no available experimental results concerning the along the chord equivalent stiffness of such kind of cables. Nevertheless, a first validation of the method can be obtained comparing the theoretical results with those obtained with some refined non-linear finite element software. With this aim, cables characterized by different sag to chord ratios and by uniformly distributed loads of varying intensity have been studied using the Cosmos/M finite element software package.
To satisfactorily explore the range of the previously considered cases, four different values have been considered for the sag to span ratio: f / a = 0.04 ; 0.1 ; 0.25 ; 1.0 and six different values have been considered for the ratio between the horizontal component of the normal stress and the elastic modulus of the cable material: σ 0 / E = 0.05 % ; 0.1 % ; 0.25 % ; 0.5 % ;   1.0 % ; 2.0 % .
The cable has been discretized using 46 non-linear 2D truss elements, in such a way that the x-distance between two consecutive nodes was constant and equal to a / 46 .
During each non-linear run, the finite element model was first loaded by the uniformly distributed loads, represented by vertical forces applied on each node; subsequently, an imposed horizontal displacement of the 0.5 %   a was imposed to the right end of the cable. The FE models are summarized in Figure 21.
The loading phase and the imposed displacement phase were both subdivided in 500 steps of constant amplitude, adopting an iterative Newton–Raphson algorithm to seek the convergence. The slope of the horizontal reaction–displacement curves in the initial stage of the imposed displacement phase allowed to derive the apparent stiffness. A typical reaction–displacement curve for flabby cables is illustrated in Figure 22, where the loading phase correspond to the first fictitious unitary time interval on the abscissa, and the imposed displacement phase to the second unit time interval. In the diagram is quite evident the hardening behavior of the cable in the displacement phase ( t 1 ) .
The numerical results are compared with the pertinent theoretical curves in the diagrams in Figure 23a, referring to the cases f / a = 0.04 , and f / a = 0.10 , and in Figure 23b, referring to f / a = 0.25 , and f / a = 1.0 , where the numerical results are represented by shaded squares. The comparison seems satisfactory, even if some small discrepancy appears, especially for small values of the σ 0 / E ratios. In any case, further comparisons are necessary to fully validate the procedure.
The loading phase and the imposed displacement phase were both subdivided in 500 steps of constant amplitude, adopting an iterative Newton–Raphson algorithm to seek the convergence. The slope of the horizontal reaction–displacement curves in the initial stage of the imposed displacement phase allowed to derive the apparent stiffness. A typical reaction–displacement curve for flabby cables is illustrated in Figure 22, where the loading phase correspond to the first fictitious unit time interval on the abscissa, and the imposed displacement phase to the second unit time interval. In the diagram is quite evident that the hardening behavior of the cable in the displacement phase.
The numerical results are compared with the pertinent theoretical curves in the diagrams in Figure 23a referring to the cases f / a = 0.04 , and f / a = 0.10 , and in Figure 23b, referring to f / a = 0.25 , and f / a = 1.0 , where the numerical results are represented by shaded squares.
The comparison seems satisfactory, even if some small discrepancy appears, especially for small values of the σ 0 / E ratios. In any case, further comparisons are necessary to fully validate the procedure.

5. Conclusions

A traditional and very effective approach to assess the non-linear behavior of cable structures consists in replacing the cables with equivalent tie rods characterized by suitable non-linear constitutive laws, which have been the subject of several research works. Aiming to enlarge the field of application of the classical formulae, based on drastic simplifications, overcoming the severe limitations of the traditional approaches, analytical expressions of the non-linear constitutive laws of equivalent tie rod are proposed for horizontal cables transversally loaded by uniformly distributed loads, distinguishing two relevant cases, according as the length of the cable can vary or not:
  • in the first case, rather new, an extremity of the cable is fixed, and the other end moves on a fixed pulley, so leaving unchanged the chord length; while
  • in the second case, the one traditionally considered, both ends of the cable are fixed, so that the unstrained length of the cable is unaffected.
The expressions, derived by applying the general form of the theorem of virtual work, can be applied independently on the cable material, on the sag to chord ratio, on the load intensity, and on the stress level, so allowing the replacement of the whole cable with an appropriate single equivalent non-linear tie rod. Moreover, the expressions can be easily modified to include long-term effects, by introducing an appropriate age-adjusted elastic modulus, for example accounting for the creep coefficient.
Starting from a wide parametric study, exhaustive sets of curves are provided and commented for each relevant operational case, underlying the capabilities of the method, also in comparison with the traditional approaches. Finally, a first validation is obtained comparing the theoretical predictions with the numerical results obtained with a refined nonlinear finite element analysis.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Zhang, X. Mechanics feasibility of using CFRP cables in super long-span cable-stayed bridges. Struct. Eng. Mech. 2008, 29, 567–579. [Google Scholar] [CrossRef]
  2. Xiong, W.; Cai, C.S.; Xiao, C. The use of carbon fiber-reinforced polymer (CFRP) composites for cable-stayed bridges. In Advanced Composites in Bridge Construction and Repair; Kim, J.L., Ed.; Woodhead Publishing: Cambridge, UK, 2017; pp. 210–264. [Google Scholar] [CrossRef]
  3. Meier, U. Carbon fiber reinforced polymer cables: Why? Why Not? What If? Arab J. Sci. Eng. 2012, 37, 399–411. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, Y.; Zwingmann, B.; Schlaich, M. Carbon fiber reinforced polymer for cable structures—A Review. Polymers 2015, 7, 2078–2099. [Google Scholar] [CrossRef] [Green Version]
  5. Micelli, A.; Cascardi, A.; Marsano, M. Seismic strengthening of a theatre masonry building by using active FRP wires. In Proceedings of the 16th International Brick and Block Masonry Conference, Padova, Italy, 26–30 June 2016; Modena, C., da Porto, F., Valluzzi, M.R., Eds.; CRC Press/Balkema: Leiden, NL, USA, 2016; pp. 753–761. [Google Scholar]
  6. Cai, H.; Aref, A.J. On the design and optimization of hybrid carbon fiber reinforced polymer-steel cable system for cable-stayed bridges. Compos. Part B Eng. 2015, 68, 146–152. [Google Scholar] [CrossRef]
  7. Dischinger, F. Hängebrücken für schwertse Verkehrslasten, I. Bauingenieur 1949, 24, 65–75. (In German) [Google Scholar]
  8. Dischinger, F. Hängebrücken für schwertse Verkehrslasten, II. Bauingenieur 1949, 24, 107–113. (In German) [Google Scholar]
  9. Ernst, J.H. Der E-modul von Seilen under Berücksichtigung des Durchhanges. Bauingenieur 1965, 40, 52–55. (In German) [Google Scholar]
  10. Irvine, H.M. Cable Structures; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  11. Como, A.; Grimaldi, A.; Maceri, F. Statical behaviour of long-span cable-stayed bridges. Int. J. Solids Struct. 1985, 21, 831–850. [Google Scholar] [CrossRef]
  12. Daniell, W.E.; Macdonald, J.H.G. Improved finite element modeling of a cable-stayed bridge through systematic manual tuning. Eng. Struct. 2007, 29, 358–371. [Google Scholar] [CrossRef]
  13. de Sá Caetano, E. Cable Vibrations in Cable-Stayed Bridges; IABSE: Zürich, Switzerland, 2007. [Google Scholar]
  14. Fleming, J.F. Nonlinear static analysis of cable-stayed bridge structures. Comp. Struc. 1979, 10, 621–635. [Google Scholar] [CrossRef]
  15. Freire, A.M.S.; Negrão, J.H.O.; Lopes, A.V. Geometrical nonlinearities on the static analysis of highly flexible steel cable-stayed bridges. Comp. Struc. 2006, 84, 2082–2140. [Google Scholar] [CrossRef] [Green Version]
  16. Gimsing, N.J. Cable Supported Bridges: Concept and Design, 2nd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 1997. [Google Scholar]
  17. Hajdin, N.; Michaltsos, G.T.; Konstantakopulus, T.G. About the equivalent elastic modulus of elasticity of cables of cable-stayed bridges. Facta Univ. Niš Univ. 1998, 1, 569–575. [Google Scholar]
  18. Li-zhong, W.; Zhen, G.; Feng, Y. Quasi-static three-dimensional analysis of suction anchor mooring system. Ocean Eng. 2010, 37, 1127–1138. [Google Scholar] [CrossRef]
  19. Pao-Hsii, W.; Hung-Ta, L.; Tzu-Yang, T. Study on nonlinear analysis of a highly redundant cable-stayed bridge. Comp. Struc. 2002, 80, 165–182. [Google Scholar] [CrossRef]
  20. Troitsky, M.S. Cable Stayed Bridges: Theory and Design, 2nd ed.; BSP Professional Books: Oxford, UK, 1988. [Google Scholar]
  21. Tschemmernegg, F.; Obholzer, A. Einfach abgespannte Seile bei Schrägseilbrücken. Bauingenieur 1981, 56, 325–330. (In German) [Google Scholar]
  22. Crusells-Girona, M.; Filippou, F.C.; Taylor, R.L. A mixed formulation for nonlinear analysis of cable structures. Comp. Struct. 2017, 186, 50–61. [Google Scholar] [CrossRef]
  23. Li, C.; He, J.; Zang, Z.; Liu, Y.; Ke, H.; Dong, C.; Li, H. An improved analytical algorithm on main cable system of suspension bridge. Appl. Sci. 2018, 8, 1358. [Google Scholar] [CrossRef] [Green Version]
  24. Tan, H.; Zeng, Y.; Zhang, X. Cable spring effect and its longitudinal restraint stiffness on towers. Proc. ICE Bridge Eng. 2020, 173, 78–85. [Google Scholar] [CrossRef]
  25. Bertrand, C.; Acary, V.; Lamarque, C.-H.; Ture Savadkoohi, A. A robust and efficient numerical finite element method for cables. Int. J. Num. Methods Eng. 2020, 121, 4157–4186. [Google Scholar] [CrossRef]
  26. Hussein, H.; Gouttefarde, M.; Pierrot, F. Static modeling of Sagging Cables with flexural rigidity and shear forces. In Advances in Robot Kinematics 2018; Lenarcic, J., Parenti-Castelli, V., Eds.; Springer International Publishing: Cham, Switzerland, 2019; pp. 310–318. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Wang, J.; Ye, G.; Xu, R. Bending stiffness of parallel wire cables including interfacial slips among wires. ASCE J. Struct. Eng. 2018, 144, 04018164. [Google Scholar] [CrossRef]
  28. Croce, P. Equivalent axial stiffness of horizontal stays. Appl. Sci. 2020, 10, 6263. [Google Scholar] [CrossRef]
  29. Croce, P. Non-linear behavior of heavy stays. Int. J. Solids Struct. 2013, 50, 1093–1107. [Google Scholar] [CrossRef]
  30. Noii, N.; Khodadadian, A.; Ulloa, J.; Aldakheel, F.; Wick, T.; Stijn, F.; Wriggers, T. Bayesian inversion for unified ductile phase-field fracture. Comput. Mech. 2021, 68, 943–980. [Google Scholar] [CrossRef]
  31. Khodadadian, A.; Noii, N.; Parvizi, M.; Abbaszadeh, M.; Wick, T.; Heitzinger, C. A Bayesian estimation method for variational phase-field fracture problems. Comput. Mech. 2020, 66, 827–849. [Google Scholar] [CrossRef] [PubMed]
  32. Croce, P.; Beconcini, M.L.; Formichi, P.; Landi, F.; Puccini, B.; Zotti, V. Bayesian Methodology for Probabilistic Description of Mechanical Parameters of Masonry Walls. ASCE-ASME J. Risk Uncertain. Eng. Syst. Part A Civ. Eng. 2021, 7, 04021008. [Google Scholar] [CrossRef]
  33. Marsili, F.; Croce, P.; Friedman, N.; Formichi, P.; Landi, F. Seismic reliability assessment of a concrete water tank based on the Bayesian updating of the finite element model. ASCE-ASME J. Risk Uncertain. Eng. Syst. Part B Mech. Eng. 2017, 3, 021004. [Google Scholar] [CrossRef]
  34. Beconcini, M.L.; Croce, P.; Marsili, F.; Muzzi, M.; Rosso, E. Probabilistic reliability assessment of a heritage structure under horizontal loads. Probabilistic Eng. Mech. 2016, 45, 198–211. [Google Scholar] [CrossRef]
  35. Yang, D.; Zhang, J.; Song, S.; Zhou, F.; Wang, C. Experimental investigation on the creep property of carbon fiber reinforced polymer tendons under high stress levels. Materials 2018, 11, 2273. [Google Scholar] [CrossRef] [Green Version]
  36. Croce, P. New frontiers of composites applications in heritage buildings: Repair of exposed masonry of St. Nicola Church in Pisa. J. Compos. Sci. 2021, 5, 218. [Google Scholar] [CrossRef]
  37. Pejis, T. 1.5 High Performance Polyethylene Fibers. In Comprehensive Composite Materials II; Beaumont, P.W.R., Zweben, C.H., Eds.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 86–126. [Google Scholar] [CrossRef]
  38. Specific Modulus. Available online: https://en.wikipedia.org/wiki/Specific_modulus (accessed on 14 August 2021).
  39. Specific Strength. Available online: https://en.wikipedia.org/wiki/Specific_strength (accessed on 14 August 2021).
Figure 1. Deformed (parabolic) configuration of a cable subject to a uniformly distributed load q .
Figure 1. Deformed (parabolic) configuration of a cable subject to a uniformly distributed load q .
Materials 14 05502 g001
Figure 2. System of forces acting on the cable portion A C .
Figure 2. System of forces acting on the cable portion A C .
Materials 14 05502 g002
Figure 3. Virtual deformations of a cable running on a fixed pulley in B : d a = 0 .
Figure 3. Virtual deformations of a cable running on a fixed pulley in B : d a = 0 .
Materials 14 05502 g003
Figure 4. Virtual deformations of a cable with ends fixed in A and B .
Figure 4. Virtual deformations of a cable with ends fixed in A and B .
Materials 14 05502 g004
Figure 5. E t , e q , 0 / E f / a σ 0 / E diagram for a cable running on a fixed pulley: d a = 0 .
Figure 5. E t , e q , 0 / E f / a σ 0 / E diagram for a cable running on a fixed pulley: d a = 0 .
Materials 14 05502 g005
Figure 6. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ,   0.1 % ) : (a) 0 < f / a     0.03 ;   (b) 0.03 < f / a     0.08 ; (c) 0.08 < f / a     0.2 ; (d) 0.2 < f / a     1.0 .
Figure 6. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ,   0.1 % ) : (a) 0 < f / a     0.03 ;   (b) 0.03 < f / a     0.08 ; (c) 0.08 < f / a     0.2 ; (d) 0.2 < f / a     1.0 .
Materials 14 05502 g006
Figure 7. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.2 % ,   2.0 % ) :   (a) 0 < f / a     0.2 ;   (b) 0.2 < f / a     1.0 .
Figure 7. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.2 % ,   2.0 % ) :   (a) 0 < f / a     0.2 ;   (b) 0.2 < f / a     1.0 .
Materials 14 05502 g007
Figure 8. E t , e q , 0 / E σ 0 / E diagrams parameterized in terms of f / a :   (a) 0.0025     f / a     0.1 ;   (b) 0.125     f / a     1.0 .
Figure 8. E t , e q , 0 / E σ 0 / E diagrams parameterized in terms of f / a :   (a) 0.0025     f / a     0.1 ;   (b) 0.125     f / a     1.0 .
Materials 14 05502 g008
Figure 9. E t , e q , 0 / E f / a σ 0 / E diagram for a cable with fixed ends.
Figure 9. E t , e q , 0 / E f / a σ 0 / E diagram for a cable with fixed ends.
Materials 14 05502 g009
Figure 10. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ,   0.1 % ) :   (a) 0 < f / a     0.025 ;   (b) 0.025 < f / a     0.05 ; (c) 0.05 < f / a     0.1 ; (d) 0.1 < f / a     1.0 .
Figure 10. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ,   0.1 % ) :   (a) 0 < f / a     0.025 ;   (b) 0.025 < f / a     0.05 ; (c) 0.05 < f / a     0.1 ; (d) 0.1 < f / a     1.0 .
Materials 14 05502 g010
Figure 11. E t , e q , a / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.2 % ,   2.0 % ) :   (a) 0 < f / a     0.2 ;   (b) 0.2 < f / a     1.0 .
Figure 11. E t , e q , a / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.2 % ,   2.0 % ) :   (a) 0 < f / a     0.2 ;   (b) 0.2 < f / a     1.0 .
Materials 14 05502 g011
Figure 12. E t , e q , a / E σ 0 / E diagrams parameterized in terms of f / a :   (a) 0.0025     f / a     0.1 ;   (b) 0.125     f / a     1.0 .
Figure 12. E t , e q , a / E σ 0 / E diagrams parameterized in terms of f / a :   (a) 0.0025     f / a     0.1 ;   (b) 0.125     f / a     1.0 .
Materials 14 05502 g012
Figure 13. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 4.0 % ,   20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 13. E t , e q , 0 / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 4.0 % ,   20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g013
Figure 14. E t , e q , a / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 4.0 % ,   20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 14. E t , e q , a / E f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 4.0 % ,   20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g014
Figure 15. E t , e q , D / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 15. E t , e q , D / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g015
Figure 16. E t , e q , D / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 16. E t , e q , D / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g016
Figure 17. E t , e q , I r v s / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 17. E t , e q , I r v s / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g017
Figure 18. E t , e q , I r v s / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1 .
Figure 18. E t , e q , I r v s / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1 .
Materials 14 05502 g018
Figure 19. E t , e q , I r v e / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 19. E t , e q , I r v e / E t , e q , 0 f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g019
Figure 20. E t , e q , I r v e / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 20. E t , e q , I r v e / E t , e q , a f / a diagrams parameterized in terms of σ 0 / E ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g020
Figure 21. Finite element model of the cable: (a) f / a = 0.04 ;   (b) f / a = 0.10 ; (c) f / a = 0.25 ;   (d) f / a = 1.0 .
Figure 21. Finite element model of the cable: (a) f / a = 0.04 ;   (b) f / a = 0.10 ; (c) f / a = 0.25 ;   (d) f / a = 1.0 .
Materials 14 05502 g021
Figure 22. R t diagram for the cable: loading phase 0 t 1 ; imposed displacement phase 1 < t 2 .
Figure 22. R t diagram for the cable: loading phase 0 t 1 ; imposed displacement phase 1 < t 2 .
Materials 14 05502 g022
Figure 23. Comparison between E t , e q , I r v e / E t , e q , a f / a theoretical curves parameterized in terms of σ 0 / E and nonlinear finite element analysis results ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Figure 23. Comparison between E t , e q , I r v e / E t , e q , a f / a theoretical curves parameterized in terms of σ 0 / E and nonlinear finite element analysis results ( σ 0 / E 0.01 % ; 2.0 % ; 10.0 % ; 20.0 % ) :   (a) 0 < f / a     0.3 ;   (b) 0.3 < f / a     1.0 .
Materials 14 05502 g023
Table 1. Relevant mechanical properties of natural and artificial cable and rope materials.
Table 1. Relevant mechanical properties of natural and artificial cable and rope materials.
Cable MaterialDensity
[kg/m3]
E
[GPa]
f t
[MPa]
f t E
Aramid fiber (high modulus)144011230002.68%
Aramid fiber (normal modulus)144070.529004.11%
Liquid crystal aromatic polyester (LCP) fiber14106628304.29%
Polybenzoxazole (PBO) fiber156027039501.46%
Carbon fiber 1156017025001.47%
Carbon fiber 2 1180024048502.02%
Steel strand785018017700.98%
Steel wire785020619000.92%
Copper89401102400.22%
Aluminum2700691050.15%
Nylon fiber11404.5661013.38%
Polyester fiber138013.87905.72%
Cotton rope15407.92252.85%
Hemp rope1490323000.94%
Flax rope1540273401.26%
Jute rope150025.82300.89%
Abaca (Manila hemp) rope1320303001.00%
Sisal1320302500.83%
Silk (silkworm)1320106506.50%
Silk (spider)1100129007.50%
1 Carbon fiber for structural applications [31].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Croce, P. Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables. Materials 2021, 14, 5502. https://doi.org/10.3390/ma14195502

AMA Style

Croce P. Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables. Materials. 2021; 14(19):5502. https://doi.org/10.3390/ma14195502

Chicago/Turabian Style

Croce, Pietro. 2021. "Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables" Materials 14, no. 19: 5502. https://doi.org/10.3390/ma14195502

APA Style

Croce, P. (2021). Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables. Materials, 14(19), 5502. https://doi.org/10.3390/ma14195502

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop