Original article Biomath 1 (2012), 1212107, 1–5 B f Volume ░, Number ░, 20░░ BIOMATH ISSN 1314-684X Editor–in–Chief: Roumen Anguelov B f BIOMATH h t t p : / / w w w . b i o m a t h f o r u m . o r g / b i o m a t h / i n d e x . p h p / b i o m a t h / Biomath Forum A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures Environmentally Damaged and Strengthened by Cable-bracings Dedicated to Academician Blagovest Sendov’s 80th anniversary Konstantinos Liolios∗, Stefan Radev†, Asterios Liolios‡, Ivan Georgiev§, Krassimir Georgiev¶ ∗Department of Environmental Engineering, Democritus University of Thrace, Xanthi, Greece kliolios@env.duth.gr †Institute of Mechanics, Bulgarian Academy of Sciences, Sofia, Bulgaria stradev@imbm.bas.bg ‡Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece aliolios@civil.duth.gr §Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Sofia, Bulgaria john@parallel.bas.bg ¶Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, Sofia, Bulgaria georgiev@parallel.bas.bg Received: 1 August 2012, accepted: 10 December 2012, published: 29 December 2012 Abstract—A computational treatment is presented for the mathematically rigorous analysis of civil engineering structures, which have been environmentally damaged and subsequently strengthened by cable-elements. The problem is treated as an inequality one, where the governing conditions are equalities as well as inequalities. The cable behavior is considered as nonconvex and nonmonotone one and is described by generalized subdifferential relations including loosening, elastoplastic - fracturing and other ef- fects. Using piece-wise linearization for the cable behavior, a linear complementarity problem, with a reduced number of unknowns, is solved by optimization algorithms. Finally, an example from Civil and Environmental Engineering praxis is presented. Keywords-Civil and Environmental Engineering; Non- Convex Analysis; Computational Mechanics; Cable-braced structures; Optimization Algorithms; I. INTRODUCTION Environmental actions can often cause significant damages to civil engineering structures, see e.g. [1], [2]. Main such defect is the strength degradation, causing a reduction of the load bearing capacity. To handle such defects, sometimes cable-like members are used as a first strengthening and repairing procedure. These cable-like members can undertake tension but buckle and become slack and structurally ineffective when subjected to a sufficiently large compressive force. Thus the governing conditions take an equality as well as an inequality form and the problem becomes nonlinear. So, the problem of structures containing as above cable-like members be- longs to the so-called Inequality Problems of Mechanics, as their governing conditions are of both, equality and inequality type [3]–[6] . A realistic numerical treatment of such problems Citation: K Liolios, S Radev, A Liolios, I Georgiev, K Georgiev, A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures Environmentally Damaged and Strengthened by Cable-bracings, Biomath 1 (2012), 1212107, http://dx.doi.org/10.11145/j.biomath.2012.12.107 Page 1 of 5 http://www.biomathforum.org/biomath/index.php/biomath http://dx.doi.org/10.11145/j.biomath.2012.12.107 K Liolios at al., A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures ... can be obtained by mathematical programming methods (optimization algorithms) [3]–[5]. The early numerical realizations of these approaches were based mainly upon the principle of minimum complementary energy. Thus, an equivalence principle for the analysis of statically undetermined structures with unilateral constraints has been proposed and proven by G. Nitsiotas [3]. Further, for the rigorous mathematical investigation of the problem, convex analysis and the variational or hemivariational inequality concept have been introduced and used, see Panagiotopoulos [6], [7] and [11]. The aim of this paper is to deal with the develop- ment of a simple numerical procedure for the rigorous analysis of civil engineering structures containing cable- like members by using a version of the direct stiffness (displacement) method of structural analysis. The present procedure is based on the finite element method and the equivalence principle, proposed by G. Nitsiotas in [3]. Using this principle, the analysis of such structures can be reduced to a Linear Complementarity Problem (LCP), which can be solved by various effective quadratic pro- gramming algorithms [5]. A numerical example shows the direct applicability on the computer and the effec- tiveness of the procedure presented herein. II. METHOD OF ANALYSIS A. Problem formulation A frame structure containing N cable-like members is considered. The structure is discretized according to the natural finite element method [4], [6]. For the cables, pin-jointed bar elements with unilateral behavior are used. Following Panagiotopoulos [6], [7]. the behavior of the cables, including loosening, elastoplastic or/and elastoplastic-softening-fracturing and unloading - reload- ing effects, can be expressed mathematically by the subdifferential relation: si(di) ∈ ∂̂Si(di) (1) Here si and di are the (tensile) force and the de- formation (elongation), respectively, of the i-th cable element, ∂̂ is the generalized gradient and Si is the superpotential function [6], [7], [10], [11]. By definition, relation (1) is equivalent to the following hemivariational inequality, expressing the Virtual Work Principle for inequality problems: S ↑ i (di, ei −di) ≥ si(di).(ei −di) (2) Here S↑i denotes the subderivative of Si , and ei, di are kinematically admissible (virtual) deformations. For numerical treatments of practical inequality problems, a piece-wise linearization is usually applied to relations (1) and (2), see e.g. [4], [8]. So, the unilateral behavior for the i-th cable-element (i = 1, . . . , N) is expressed by the following relations [3], [4]: ei = F0isi + ei0 −vi (3) si ≥ 0, vi ≥ 0, si.vi = 0. (4) Here ei , F0i, si, ei0, and vi denote the strain (elonga- tion), natural flexibility constant, stress (tension), initial strain and slackness, respectively. From (3) it is clear that the slackness vi can be considered as an unknown initial strain which constitutes a reversible negative elongation [3]. Further, relations (4) express that either a non- negative stress-tension or a non- negative slackness exists on any cable. For the remaining structure (besides the cables), the usual frame finite element models, which exhibit a bilateral behavior, are used. B. Numerical solution approach Now the equivalence principle, proposed by G. Nitsio- tas in [3], is applied for the whole structure. According to this principle, the structure under consideration be- haves as an equivalent, linearly elastic structure, under the condition that in each cable-element either a non- negative stress or a fictitious, unknown, non- negative slackness appears-see relations (3). Thus, collecting in (N×1) vectors t and v the stress and slackness behavior of all the N cable-elements, corresponding, the following Linear Complementarity Conditions hold: t ≥ 0, v ≥ 0, tT v = 0. (5) Further, following the Stiffness (Displacement) Method of Structural Analysis, we consider the cable- element as solidified rods and we assume that the so-modified structure is a statically stable one with bilateral rod-elements. So, the tension vector t is decomposed as follows [9]: t = Cv + t0 (6) Here t0 is the stress vector of the solidified cable- elements, now acting as normal bilateral rods, due to external actions and C is the natural influence matrix of v on t. For both it is assumed a linearly elastic, bilateral behavior for the stable structure, where. the cables are considered as already solidified bars. So, the natural Biomath 1 (2012), 1212107, http://dx.doi.org/10.11145/j.biomath.2012.12.107 Page 2 of 5 http://dx.doi.org/10.11145/j.biomath.2012.12.107 K Liolios at al., A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures ... stiffness matrix C is symmetric and in general positive semi-definite. Thus, if t0 is known, then vectors t and v can be determined by solving the Linear Complementarity Problem (LCP) formed by relations (5) and (6). For the solution of this problem, various effective algorithms are available [3]–[10]. Most of these algorithms reduce the above linear complementarity problem to a quadratic programming one of the form: Min { 1 2 vT Cv + vT t0 : v ≥ 0 } (7) After the previous preparation we can now formulate the following numerical procedure for the analysis of structures containing cable like members: a) Considering the cables as having been solidified (normal bilateral bars), the vector t0 due to external actions is determined by the Finite Element Method. b) Under the same assumption and by the same method as in a), the influence matrix C is determined. In this matrix, Cij is the stress (axial force) in the solidified cable-element i caused by a unit-shortening vj = 1 imposed in the solidified cable-member j, (i, j = 1, . . . , N). c) The Linear Complementarity Problem of relations (5) and (6) is solved to provide the sought vector v. So it is computed which cable-elements are activated (under tension) and which are not (under non-zero slackness). d) The final stress state of the structure is determined by taking into account the external actions and the computed forces t of the active cable-elements. Thus, the whole procedure requires the linear elastic analysis of the modified (with solidified cable-elements) structure (N + 2) times, where N is the number of the cables, and the solution of a quadratic programming problem or a LCP. Alternatively, after having computed t, the structure is analyzed due to external actions by omitting the slack cables for which the step c) has given zero tension values. III. NUMERICAL EXAMPLE The reinforced concrete plane frame structure of Fig. 1 had been initially analyzed, designed and constructed to bear the shown loads, without the shown cable- elements. The concrete class is C40/50 and the elasticity modulus Eb = 3.5×107kN/m2. The shown dimensions width/height of rectangular sections are in centimeters [cm]. The analysis of the above frame without cable- elements is obtained by using any available finite element method code, e.g. SAP2000 code [13]. Due to environmental actions [1], [2], corrosion and cracking had been taking place. This had caused a reduc- tion for the section inertia moments, which, according to [14], is estimated to be 10% for the columns and 50% for the beams. So it was necessary for the system to be strengthened. As a first repairing and strengthening procedure, ten (N = 10) cable members with cross- sectional area Fr = 8 cm2 have been added as shown in Fig. 1. The steel class is S1400/1600 with an elasticity modulus Es = 210GPa. These cables are placed as counter diagonals. As it is not known in advance which of them are activated or not by the given loads, the purpose here is to compute what happens. The application of the presented numerical procedure gives first the values of the slackness of the not activated cable-elements: v1 = 0.848×10−3m, v3 = 10.321×10−3m, v5 = 1.082×10−3m, v8 = 9.564×10−3m, v10 = 1.652×10−3m. Further, the elements of vector t, where t = [S1, S2, . . . , S10]T , are computed to have the following values (in kN) for the non-active cables: S1 = S3 = S5 = S8 = S10 = 0.0, whereas for the active cables it holds: S2 = 10.17kN, S4 = 346.04kN, S6 = 18.84kN, S7 = 342.08kN, S9 = 25.81kN. Thus, cables 2,4,6,7 and 9 are the only ones active, having zero slackness. The other cables 1,3,5,8 and 10 cannot contribute to the system resistance under the given loads in Fig. 1. Using the previous results, the final stress state is computed. In Fig. 2 is shown indicatively the final bending moments diagram for the strengthened frame containing the active cable-elements only. Comparing the diagram in Fig. 2 with the corresponding one for the initial frame without cables, the efficiency of the strengthening can be checked. IV. CONCLUSION The inequality problem of the cable-braced civil en- gineering structures can be treated numerically by the herein presented approach. This approach takes into ac- count the unilateral behavior of cable elements, uses the equivalence principle of G. Nitsiotas [3] and so leads to a linear complementarity problem, with a reduced number of problem unknowns. Thus, the numerical realization of Biomath 1 (2012), 1212107, http://dx.doi.org/10.11145/j.biomath.2012.12.107 Page 3 of 5 http://dx.doi.org/10.11145/j.biomath.2012.12.107 K Liolios at al., A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures ... 4,6 m 3,4 m 3,1 m 4,9 m 18 KN 180 KN 4,1 m 3,0 m 24 KN 32KN/m 52 KN/m 220 KN 8 KN/m 30/40 30/50 30/60 30/60 30/60 30/45 30/50 30/65 30/55 30/70 30/65 30/50 30/30 25/25 25/30 30/35 10 9 1 2 3 4 5 6 7 8 Fig. 1. The cable-braced structural system of the numerical example. Fig. 2. Bending moments diagram (in kNm) for the frame with the 5 active cable-elements No 2,4,6,7 and 9. Biomath 1 (2012), 1212107, http://dx.doi.org/10.11145/j.biomath.2012.12.107 Page 4 of 5 http://dx.doi.org/10.11145/j.biomath.2012.12.107 K Liolios at al., A Linear Complementarity Numerical Approach to the Non-Convex Problem of Structures ... the proposed approach is obtained by available computer codes of the finite element method and of mathematical programming algorithms. Moreover, as it has been ver- ified in an example, the herein developed approach can treat in a realistic way the general problem of cable- braced structures in civil and environmental engineering praxis. An extension of the presented approach for the case of the earthquake response of such structures has been recently presented [12]. ACKNOWLEDGMENT This work was partially supported by the Bulgarian Science Fund under the project DMU 03-62. REFERENCES [1] H. S. Peavy, D. R. Rowe, and G. Tchobanoglous, Environmental Engineering, McGraw-Hill, New York, 1985. [2] A. Moncmanova, Environmental Deterioration of Materials, WIT Press, Southampton, 2007. http://dx.doi.org/10.2495/978-1-84564-032-3 [3] G. Nitsiotas, Die Berechnung statisch unbestimmter Tragwerke mit einseitigen Bindungen, Ingenieur-Archiv, vol. 41, pp. 46-60, 1971. http://dx.doi.org/10.1007/BF00536162 [4] G. Maier, A Quadratic Programming Approach for Certain Classes of Non-linear Structural Problems, Meccanica , vol. 3, pp. 121-130, 1968. http://dx.doi.org/10.1007/BF02129011 [5] R. W. Cottle, Fundamentals of Quadratic Programming and Lin- ear Complementarity, in M.Z. Cohn, D.E. Grierson and G. Maier (Eds.), ”Engineering Plasticity by Mathematical Programming”, Pergamon Press, Oxford, 1979. [6] P. D. Panagiotopoulos, Inequality Problems in Mechanics and Applications. Convex and Nonconvex Energy Functions, Birkhuser Verlag, Basel, Boston, 1985. http://dx.doi.org/10.1007/978-1-4612-5152-1 [7] P. D. Panagiotopoulos, Hemivariational Inequalities. Applica- tions in Mechanics and Engineering, Springer Verlag, Berlin, 1993. http://dx.doi.org/10.1007/978-3-642-51677-1 [8] A. A. Liolios, A Linear Complementarity Approach for the Non- convex Seismic Frictional Interaction between Adjacent Struc- tures under Instabilizing Effects. Journal of Global Optimization, vol. 17, no. 1-4, pp. 259-266, 2000. http://dx.doi.org/10.1023/A:1026789817828 [9] Ang. Liolios, Numerical Simulation of Structures with Unilateral Constraints and Solution by the Finite Element Method, Graduate Degree Thesis, Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece, 2007. [10] E. S. Mistakidis and G. E. Stavroulakis, Nonconvex optimization in mechanics. Smooth and nonsmooth algorithmes, heuristic and engineering applications, Kluwer, London, 1998. [11] R. T. Rockafellar, Convex Analysis, Princeton Univ. Press, Princeton, 1970. [12] Ang. Liolios, K. Chalioris, Ast. Liolios, St. Radev, and Kon. Li- olios, A Computational Approach for the Earthquake Response of Cable-braced Reinforced Concrete Structures under Environ- mental Actions. In: Lirkov, I., Margenov, S. and Wasniewski, J. (Eds): ”Large-Scale Scientific Computing”, LNCS 7116, pp. 590-597, Springer-Verlag, Berlin Heidelberg, 2012. http://dx.doi.org/10.1007/978-3-642-29843-1 [13] Computers and Structures, Inc. SAP2000, Integrated Finite Element Analysis and Design of Structures, Berkeley, 2005. [14] T. Pauley and M. J. N. Priestley, Seismic design of reinforced concrete and masonry buildings, Wiley, New York, 1992. http://dx.doi.org/10.1002/9780470172841 Biomath 1 (2012), 1212107, http://dx.doi.org/10.11145/j.biomath.2012.12.107 Page 5 of 5 http://dx.doi.org/10.2495/978-1-84564-032-3 http://dx.doi.org/10.1007/BF00536162 http://dx.doi.org/10.1007/BF02129011 http://dx.doi.org/10.1007/978-1-4612-5152-1 http://dx.doi.org/10.1007/978-3-642-51677-1 http://dx.doi.org/10.1023/A:1026789817828 http://dx.doi.org/10.1007/978-3-642-29843-1 http://dx.doi.org/10.1002/9780470172841 http://dx.doi.org/10.11145/j.biomath.2012.12.107 Introduction Method of Analysis Problem formulation Numerical solution approach Numerical Example Conclusion References