FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007 DIAGONAL STRESS-STRAIN MODEL FOR FRP-CONFINED RECTANGULAR CONCRETE COLUMNS Domingo A. MORAN 1 1 2 Chris P. PANTELIDES 2 Ph.D. Candidate, University of Utah, Salt Lake City, Utah, USA Professor, University of Utah, Salt Lake City, Utah, USA Keywords: concrete columns; composite materials; ductility; rehabilitation; retrofit; seismic hazard; stress-strain relations. 1 INTRODUCTION Despite the development of the ACI440.2R-02 [1] design guidelines and the successful application of Fiber Reinforced Polymer (FRP) jacketing systems for circular, square and rectangular columns, research into the constitutive relationships governing the dilation and compressive behavior of rectangular FRP-confined concrete is continuing [2,3,4,5]. For circular FRP-confined concrete sections subjected to uniform axial compression the confining FRP jacket provides a uniform kinematic restraint at the surface of the concrete core; and subjects the concrete to a uniform triaxial compressive state of stress. However, for rectangular and square FRP-confined concrete sections, the effectively confined concrete core is in a complex triaxial compressive state of stress. The thin rectangular FRP jacket provides a high axial kinematic restraint along the faces of the rectangular section as a result of its high axial stiffness, and a very low out-of-plane kinematic restraint as a result of its low flexural stiffness. As a result of the geometry of the rectangular FRP jacket, and the difference between the axial and flexural stiffness of the jacket, the concrete core is subjected to a non-uniform triaxial compressive state of stress. The resultant passive confining pressure provided by the FRP jacket decreases significantly as the aspect ratio of the rectangular FRP-confined (FCC) or concrete filled FRP tube (CFFT) section increases, as it is concentrated along the diagonal axis of the rectangular FCC or CFFT section. Despite the development of stress-strain models for circular and square FRPconfined concrete sections there is a lack of models that can accurately predict the dilation and compressive behavior of rectangular FRP-confined concrete sections of varying aspect ratios and jacket corner radii. In this paper, a new and unique stress-strain model for FRP-confined concrete sections is introduced in which the concept of diagonal dilation and equilibrium of the FRP-confined concrete is considered. A variable Poisson’s ratio model, developed by the authors [6], in combination with analytical strain compatibility and diagonal equilibrium relationships developed herein for rectangular FRP-confined concrete sections is combined with a Popovics [7] type uniaxial stressstrain model for FRP-confined concrete sections, using the following concepts: (1) The increase in compressive strength and ductility of FRP-confined concrete are expressed in terms of the internal damage resulting from diagonal dilation of concrete and lateral restraint provided by the FRP jacket; (2) The dilation behavior of FRP-confined concrete sections in compression depends on the mechanical properties of the unconfined concrete, the lateral restraint provided by the FRP jacket, the geometry of the concrete section, and the extent of internal damage in the concrete core. Modeling of both the compressive and dilation behavior of FRP-confined sections is achieved using a minimal number of experimentally obtained coefficients, which were found to affect both the compressive and dilation behavior. 2 CONFINEMENT EFFECTIVENESS OF FRP-CONFINED CONCRETE The confinement effectiveness of the confining element k cc , is measured by how much the strength, f cc , and strain, ε cc , of the confined concrete is increased over the unconfined concrete strength, f co and strain ε co using the following empirical relationships [8]: 1 ε f f k cc = cc = 1 + k1k r ; kεc = cc = 1 + k 2 k r ; k r = r f co ε co f co (1) where k1 and k 2 are experimentally determined confinement and strain effectiveness coefficients, respectively; k εc = strain effectiveness; k r = confinement ratio; and f r = average confining pressure. Various investigators have proposed nonlinear-failure criteria for concrete, most of which are empirical in nature. Considering that the structure of concrete material within the confined concrete core is heterogeneous, consisting of granular aggregates (i.e. crushed stone or gravel and sand), a binding material (i.e. cement paste) and pores. Treating the concrete core as a Mohr-Coulomb frictional cohesive material with some inter-particle attraction or apparent cohesion cc ; the concrete’s resistance to deformation is developed from a combination of particle sliding, rolling, crushing and void compaction that develops during the deformation of the concrete core, which depend on the average angle of internal friction φc of the confined concrete. Fig. 1 Mohr-Coulomb failure criteria for frictional-cohesive materials (a) Mohr’s failure stress circles including terminology (b) Mohr-Coulomb failure envelope versus normal stress. τ c , of the concrete material as shown in Fig. 1, can be approximated by: (σ − σ 3 ) φ (2) τc = 1 cos(2θ c ) = cc + σ n tan φ c ; θ c = 45° + c The shear strength, 2 2 σn = where σ1 = σ1 + σ 3 2 + σ1 − σ 3 2 major principal compressive stress; cos(2θ c ) σ 3 = minor (3) principal compressive stress; σ n = θ c = average angle of inclination of the failure shear plane. For the cases of concrete in uniaxial compression for which σ 3 = 0 and σ 1 = f cc = f co , and normal stress on the failure shear plane; σ 1 = f cc = 0 , and solving for the major compressive stress σ 1 using Eqs. (2) and (3), setting σ 1 = f cc and cc = cu , yields: uniaxial tension for which σ 3 = − f to ⎛ 1 + sin φc f cc = σ 3 ⎜⎜ ⎝ 1 − sin φc and ⎞ ⎛ cos φc ⎟⎟ + 2cu ⎜⎜ ⎠ ⎝ 1 − sin φc 2 ⎞ ⎟⎟ = σ 3 tan 2 (θ c ) + 2cu tan (θ c ) ⎠ principal (4) FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007 cu = f to f = co 2 2 kt kt ; kt = f to ; f to = 0.32( f co )2 3 f co (5) where cu is the apparent uniaxial cohesion intercept of the concrete core, as shown in Fig. 1; k t = ratio of the uniaxial tensile strength f to to the uniaxial compressive strength f co of the concrete core. Using Eqs. (4) and (5), setting k 3 = σ 3 f co , k 3e = k e k 3 , and setting k 3 = k 3e , the confinement effectiveness of the confined concrete core k cc of Eq. (1), can be found in terms of its frictional-cohesive properties and the effective minor principal compressive stress ratio k 3e in the rectangular FRP confined concrete section, as follows: f 1 + sin φ c k cc = cc = k t k1 + k 3e k1 ; k1 = = tan 2 θ c f co 1 − sin φ c k 3e = α uψ 3 k de ; α u = (6) σ 1 ; ψ 3 = 3e = [1 + cos(2θ d − 90°) + sin (2θ d − 90°)] (7) f de 2 (α sh )2 1 α sh = Rj Hc ε H = = tan (θ d ) ; α j = Bc Hc εB (8) ⎛ 2t j k e ⎞ f ⎟E j k de = de = E je ψ sh ε d = K jeγ sh ε d ; E je = ⎜⎜ ⎟ B B f co ⎝ Hc ⎠ C je 2ψ sh α sh 1 ; γ sh = ; K je = ; C je = C sh E je ψ sh = B C sh f co 2 2 β d 1 − α jα sh ( ) ( ) ) 2 ⎤ 2 ⎡ (1 − 2α j α sh ) ⎥ ke = 1 − ⎢ 3 ⎢1 − α sh (4 − π )(α j )2 ⎥ ⎦ ⎣ ( ) ( β d = β sh β eff ; β sh = 1 + where σ 3e = α u = uniformity 1 αu ; C sh = ; β eff = (1 + α sh ) − (4 − π )α jα sh 1 − α sh (4 − π )(α j )2 1 2 − eeff ; eeff = 1 − α u (9) (10) (11) (12) effective minor principal compressive stress provided by the FRP jacket; coefficient; ψ3 = minor principal stress coefficient; compressive stress in the confined concrete core; (E je )B = σ3 = minor principal effective FRP jacket transverse Bc of the section; ψ sh = diagonal equilibrium coefficient; K je = normalized effective FRP jacket confinement stiffness; γ sh = jacket equilibrium coefficient; k e = confinement efficiency of the FRP jacket; C je = effective FRP jacket stiffness along the minor dimension confinement stiffness; α sh = section aspect ratio; αj = jacket corner aspect ratio; εd = ε c ; C sh = jacket β sh = aspect ratio diagonal strain in the FRP-confined concrete section at any given axial strain, confinement ratio coefficient; βd = diagonal strain compatibility coefficient; β eff = jacket coefficient; eeff = eccentricity of the effectively confined elliptical concrete core. The effectively confined concrete core of the section of Fig. 2, is modeled as an equivalent elliptical core having dimensions H eff , Beff and a distance between the foci of C eff , coefficient; as shown in Fig 2(a). The concrete section is confined by an FRP jacket having a normalized 3 effective stiffness K je , subjected to an axial compressive strain diagonal strain εd , εc , along the diagonal axis D j , of the section. inclination of the main diagonal D j ; εB and εH and a resultant expansive In addition, θd = angle of are the average FRP jacket strains along the minor and major dimensions of the section, respectively. Fig. 2 Typical Rectangular FRP-confined concrete section (a) Typical section geometry (b) Diagonal equilibrium of section. Note that the k t k1 and k 3e k1 terms of Eq. (6), can be considered to be the apparent cohesive intercept parameter and apparent frictional slope parameter, respectively, of the confinement effectiveness k cc , of confined concrete. Also, note that k cc of Eq. (6) yields k cc of Eq. (1) only k1 = (k1 )u = 1 k t . As a result, the confinement effectiveness k cc of Eq. (6), can be considered to be the lower bound value of k cc at all levels of confinement and k cc of Eq. (1), can when be considered to be its upper bound value. 2.1 Nonlinear Failure Envelopes of FRP-confined concrete For a Mohr-Coulomb type maximum strength failure envelope, such as that introduced in Eqs. (2)-(12) to properly model the behavior of concrete in a triaxial compression state of stress, it must also include the cases of uniaxial tension, tension-compression and tension-tension stress states and shall satisfy the following conditions: (1) It should pass through the point of uniaxial tension, i.e. when σ 3 = − f to , σ 1 = 0 ; (2) It should pass through the point of uniaxial compression, i.e. when σ 3 = 0 , σ 1 = f co ; (3) It should fit the failure envelope of concrete in a triaxial compression state of stress, and (4) It should approach the hydrostatic state when the confinement is very large (i.e. as σ 2 = σ 3 → ∞ , σ 1 → σ 2 = σ 3 ), where σ 2 is an intermediate principal compressive stress. The following Mohr-Coulomb failure envelope developed herein for confined concrete satisfies these conditions. This model considers that the inherent curvature of the failure envelopes of confined concrete, as shown in Fig. 1(b), can be attributed to the remolding of the internal structure of the concrete core as damage in the concrete progresses which is attributed to the degradation of the angle of internal friction φ c [9] of the confined concrete core, as shown in Fig. 3. As the 4 FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007 internal structure of the concrete core degrades to its constituent materials (from internal crack growth due to loss of adhesion provided by the cement paste), and is found to depend on the state of stress ( σ n , σ 1 and σ 3 ) in the confined concrete core as follows: Fig. 3 Plot of friction angle versus the normal stress on the failure shear plane of the confined concrete core. ⎛ (k1 )be − 1 ⎞ ⎟ ; ψ dil = φ o − φbe ⎟ kn ⎝ (k1 )be + 1 ⎠ 1+ (k n )m φ c = φ be + ψ dil ; φbe = sin −1 ⎜⎜ (k1 )be = neff (13) φ + φbe 1 + sin φ b + (1 − neff )(k1 )b ; (k1 )b = (14) ; φm = o 1 − sin φb 2 φ k φ o = φu + σ (φ u − φbe )[1 + cos(2θ u )] ; φ u = 90° − 2 tan −1 k t ; θ u = 45° + u 2 kt kσ = cu (σ n )m = k t (φ t − φ u ) 2k t (φ t − φ be ) + (φ u − φ be )[1 + cos(2θ u )] ; (k n )m = where the basic angle of internal friction of the concrete φb (σ n )m f co = (15) (φu − φbe )[1 + cos(2θ u )] (16) 2(φ o − φ u ) is a material constant of the concrete core that depends on the aggregate type (crushed stone or gravel), aggregate gradation and sandaggregate ratio of the concrete mix; for normal-weight normal strength concrete an average value of φ b ≅ 36° is used herein, which was determined from analysis of dry triaxial compression concrete cylinder tests [10]. In addition, ψ dil = angle of dilatancy of the concrete which depends on the state of stress of the concrete core, and is governed by Eq. (13). In reference to Fig. 3, φ o = initial angle of internal friction of the concrete core corresponding to the initial normal stress (σ n )o = 0 , on the shear failure plane at the instant when inclination of the failure shear plane corresponding to coefficient; (14); (σ n )m = (k n ) m (σ n )o ; σ n = 0 ; θo = initial angle of kσ = median angle normal stress median-angle normal stress corresponding to the median angle φm of Eq. = normalized median-angle normal stress; k n normalized normal stress in the concrete core. In addition, as shown in Figs. 1 and 3, the apparent uniaxial cohesion cu , and the uniaxial normal stress (σ n )u and the uniaxial angle of internal friction 5 φ u correspond to the uniaxial compression case for which k 3e = σ 3e = 0 and σ 1 = f co . Excess pore water effects [10,11] are accounted for in the angle of friction model of Eqs. (13)-(16) by utilizing the effective basic confinement effectiveness coefficient (k1 )be of Eq. (14), the corresponding effective basic angle of internal friction φbe of Eq. (13), and the effective porosity, neff , of the granular concrete core with 0 < neff ≤ 1.0 . An average value of an average value of neff = 2.5% is used herein for FCC sections and neff = 55.5% is recommended for CFFT sections. 2.2 Elastic and Plastic confinement effectiveness of FRP-confined concrete Typical plots of the confinement effectiveness k c , versus the effective minor principal stress k 3e , are shown in Fig. 4. Using Hooke’s law prior to cracking of the concrete core (i.e. the elastic region), the elastic confinement effectiveness k ce , as shown in Figs. 4(b) and 4(d), can be ratio found as follows: f E Eco 1 − µ do k ce = ce = −ω E ε d ; ω E = k1e K o = do ; k1e = ; Ko = (1 − µ do ) f co f co f co µ do µ do = − [ ] ∂ε d E = co = β d (µ B )o ; Eco = Eci 1 + α Eν ci (µ B )o (β sh )2 ; α E = ∂ε c E do (µ B )o =− ∂ε B = ∂ε c ν ci [ 1 + α E (1 + ν ci )1 − ν ci (β sh )2 (17) (E je )B (18) Eci ; Eci ≈ 3320 f co + 6900 (19) ] where f ce = average “elastic” stress in the concrete core; ω E = effective elastic confinement index, which represents the average initial “elastic” slope of the confinement effectiveness versus the effective minor principal stress ratio curves of Figs. 4(b) and 4(d); k1e = elastic confinement effectiveness coefficient, and K o = normalized initial stiffness of the FRP-confined concrete. Also, E co and E do are the initial diagonal and axial elastic modulus of the FRP-confined concrete section, respectively; µ do and minor εB (µ B )o are the initial dilation rates along the diagonal εd , and strain directions, respectively. In addition, ν ci = initial Poisson’s ratio of the unconfined concrete core, where 0.15 ≤ ν ci ≤ 0.25 ; the value ν ci = 0.18 is used herein; αE = modular ratio; E ci = initial tangent modulus of the unconfined concrete core. As shown in Figs. 4(a) and 4(c), in the plastic compressive region ( i.e. when ε co ≤ ε cp ≤ ε cu , and when ε dp where ε do ≤ ε dp ≤ ε du FRP jacket ε c = ε cp where ), the plastic confinement effectiveness, k cp , of the FRP-confined concrete section can be found by rewriting k cc and k εc of Eq. (1) as follows: k cp = ϕ1 p = f cp f co k1 p k1 = 1 + k1 p k 3e ; kεp = ; ϕ2 p = k2 p (R p )i k1 ε cp ε co ; ϕR = = 1 + k 2 p k 3e ; R p = ϕ1 p ϕ2 p = 6 (R p )i Rp k2 p k1 p ; ϕµ = ⎛ k2 p ; Rp = ⎜ i ⎜k ⎝ 1p ( ) µd µ do (20) ⎞ ⎟ = 3 1 + α u (21) ⎟ ⎠i ( ) FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007 Fig. 4 Typical dilation behavior of moderate stiffness FRP-confined concrete sections (a) Plot of normalized dilation and confinement coefficients and (b) Plot of confinement effectiveness-versus minor principal compressive stress-ratio where f cp = axial plastic compressive stress at a given plastic compressive strain corresponding plastic diagonal expansive strain ε dp ; k1 p = ε cp , and plastic confinement effectiveness coefficient; k 2 p = plastic strain coefficient; R p = plastic strain ductility ratio; k εp = plastic strain ( )i of effectiveness [13,14,15]. The ideal plastic strain ductility ratio of FRP-confined concrete R p ( )i = 6.0 for concrete subjected to a triaxial compression state Eq. (21), is considered equal to R p of stress [10], as occurs when the section becomes more circular as For high aspect ratio sections for which ratio, α sh >> 2.0 and α sh → 1.0 and α j → 0.50 . α j → 0 , the ideal plastic strain ductility (R p )i , approaches that of concrete in a biaxial compression state of stress for which (R p )i = 3.0 [12]. ( ) In Figs 4(a) and 4(c) it can be observed that the ideal plastic confinement effectiveness k cp i represents the tangent of the confinement effectiveness k c , versus the effective minor principal stress ratio k 3e , curve of the FRP-confined section when ( )i can be given found by using by setting k cc = k cp yields: 7 ϕ R = 1.0 (i.e. when ϕ1 p = ϕ 2 p ); and in Eq. (6) and using Eqs. (7) and (9) which (k cp )i = ( f cp )i (22) = k t k1 − ω p ε d ; ω p = α uψ 3γ shω je ; ω je = K je k1 f co where ω p = ideal plastic confinement index of the rectangular FRP jacket, ω je = effective confinement index of the FRP jacket. The plastic confinement index ωp of Eq. (22) of the rectangular FRP jacket (i.e. the tangent plastic slope shown in Figs. 4(b) and 4(d)), indicates the upper bound tangent slope of the plastic confinement effectiveness k cp versus the effective minor principal stress ratio k 3e , curve of the FRP-confined concrete when ϕ R = 1.0 or (i.e. when ϕ1 p = ϕ 2 p ). The softening and flattening of the ϕ1 p , ϕ 2 p and ϕ R curves in Figs. 4(a) and 4(c), can be attributed to the softening of the mechanical response of the FRP-confined concrete that results from remolding of the internal structure of the passively confined concrete core due to transverse dilation of the confined concrete at high axial deformations as damage builds up. As the lateral stiffness of the confining FRP jacket increases (i.e. higher the kinematic restraint) the confinement effectiveness k c and the normalized plastic confinement coefficient ϕ1 p curves shift from softening behavior to essentially bilinear hardening behavior. 3 UNIAXIAL STRESS-STRAIN MODEL FOR FRP-CONFINED CONCRETE The typical axial compressive stress-strain behavior of rectangular FCC and CFFT members can be modeled using a Popovics [7] type fractional stress-strain model introduced by the authors [6] with the use of the ideal plastic strain ductility ratio R p of Eq. (21) and the ideal plastic i confinement effectiveness k cp of Eq. (22). In addition, the diagonal dilation of the rectangular i FRP confined concrete section can be modeled using a damage based dilation model introduced by the authors [6] which considers that the dilation behavior of the FRP confined section is dependent on the stiffness of the FRP composite jacket, the geometry of the section, the mechanical properties of the unconfined concrete and FRP jacket, and the extent of internal damage as measured by the expansive diagonal strain, ε d , in the confined concrete core. In the ( ) ( ) analytical expressions developed herein, the only coefficients that were experimentally obtained and used in the uniaxial stress-strain model are the average value of neff = 55.5% recommended for CFFT sections, the average value of the basic angle of friction concrete core, and the diagonal plastic dilation rate, φb = 36° of the confined µ dp , of the dilation model introduced by the authors [6]. 4 COMPARISON OF MODEL WITH EXPERIMENTAL RESULTS The stress-strain curves obtained using the damage based model presented in Eqs. (2)-(22) in addition to the dilation and stress-strain model introduced by the authors [6] are compared to the compressive tests of square and rectangular FCC tests [16] in Fig. 5(a), and [2] in Figs. 6(a) and 6(b), and to circular FCC tests [16] in Fig. 5(b). From analysis of Figs. 5 and 6, it can be observed that the proposed damage-based model can accurately capture both the strain softening behavior of low FRP jacket stiffness rectangular and circular FRP-confined concrete sections, and the strain hardening behavior of moderate to high FRP jacket stiffness rectangular and circular FCC concrete sections, as well as rectangular and circular CFFT sections (not included herin dure to space limitations). In addition, the model can accurately capture the effects of an increase in the aspect ratio of the rectangular FRP-confined concrete section. 8 FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007 Fig. 5 Comparison of analytical model predictions versus tests of (a) square and rectangular and (b) circular FCC sections [16]. Fig. 6 Comparison of analytical model predictions versus tests of (a) square and (b) rectangular FCC sections [2]. 5 CONCLUSION A uniaxial stress-strain model using basic principles of mechanics was introduced for describing the compressive behavior and dilation behavior of circular and rectangular FRPconfined (FCC) concrete sections confined by bonded FRP jackets with rounded corners, and nonbonded (CFFT) sections. The distinguishing features of the proposed model are: (a) the use of a damage-based variable diagonal secant Poisson’s ratio and dilation rate model for FRP-confined concrete, which a function of the stiffness of the confining FRP composite jacket, the geometry of the FRP-confined concrete section, the mechanical properties of the unconfined concrete core and FRP jacket, and the extent of internal damage as measured by the expansive diagonal strain in the FRP-confined concrete; and (b) the introduction of a Mohr-Coulomb type maximum strength failure envelope for confined concrete, in which the confinement effectiveness of the confined concrete was found to depend on the uniaxial tensile and compressive properties of the concrete material, and the internal stress-dependent angle of internal friction of the confined concrete core. 6 ACKNOWLEDGEMENT The authors would like to acknowledge financial support from the National Science Foundation under contract No. CMS 0099792. The opinions expressed in this article are those of the writers, and do not necessarily reflect the opinions of the sponsoring organization. 9 REFERENCES [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] ACI440.2 R-02 (2002) “Guide for the Design and Construction of Externally Bonded FRP systems for Strengthening Concrete Structures” ACI Committee 440, October 2002. Lam, L., and Teng, J.-G. (2003) “Design-oriented Stress-Strain Model for FRP-Confined Concrete in Rectangular Columns”, J. of Reinforced Plastics and Composites, 22 (13), 2003, 1149-1186. Marques, S.P.C., Marques, D.C.D.C., daSilva, J.L., and Cavalcante, M.A.A. (2004) “Model for Analysis of Concrete Confined by Fiber-Reinforced Polymer” J. of Composites for Construction, 8 (4), 332-340. Carey, S. A. and Harries, K.A. (2003) “The effects of shape, “gap”, and scale on the behavior and modeling of variably confined concrete” Rep. No. ST03-05, Department of Civil and Envr. Engrg., U. of South Carolina, pp. 189 Mirmiran, A., Zaggers, K. and Yuan, W. (2000). “Nonlinear finite element modeling of concrete confined by fiber composites” Finite Elements in Analysis and Design, 35, 7996.Pantazopoulou, S. J., and Mills, R.H. (1995). “Microstructural aspects of the mechanical response of plain concrete.” ACI Mater. J., 92(6), 605-616. Moran, D.A., and Pantelides, C.P. "Stress-Strain Model for Rectangular FRP-Confined Concrete: A Damage Approach." Civil Engineering Infrastructure Systems (CEIS 2006), Harajli, M.H., Harik, I.E., and Naaman, A. E., (eds.), American University of Beirut, Beirut, Lebanon, June 12 - 14, 2006, Paper 49, 12 pp. Popovics, S. (1973),” A Numerical Approach to the Complete Stress-Strain Curve of Concrete” Cement and Concrete Research, 1972, Vol. 3. Pp. 583-599. Richart, F.E., Brandtzaeg, A., Brown, R.L., (1928), “A study of Failure of Concrete under Combined Compressive Stresses”, Engrg. Experiment Station Bulletin No. 185, Urban, Ill. Maksimovic, M. (1996) “A Family of Nonlinear Failure Envelopes for Non-Cemented Soils and Rock Discontinuities” Electronic J. of Geotechnical Engrng, pp. 169. Imran, I., (1994) “Applications of Non-associated Plasticity in Modelling the Mechanical Response of Concrete”, PhD. Dissert., Dept. of Civil Engrg., Univ. Toronto, 1994, 326 pp. Imran I., and Pantazopoulou, S.J. (1996) “Experimental Study of Plain Concrete Under Triaxial Stress” ACI Materials J. 93 (6), 589-601. Darwin, D. and Pecknold, D. (1977) “ Nonlinear Biaxial Stress-Strain Law for Concrete”, Journal of the Engineering Mechanic Division, ASCE, Vol. 103, No. EM2., April 1977, pp. 229-241. Moran, D.A., and Pantelides, C.P. (2002a) “Stress-Strain Model for Fiber-Reinforced Polymer-Confined Concrete”, J. of Composites for Construction, 6 (4), 2002, 233-240. Moran, D.A., and Pantelides, C.P. (2002b) “Variable Strain Ductility Ratio for FiberReinforced Polymer-Confined Concrete”, J. of Composites for Construction, 6 (4), 224-232. Moran, D.A., and Pantelides, C.P. (2005) “Damage Based Stress-Strain Model for FiberReinforced Polymer Confined Concrete.” ACI Structural J. 102(1), 54-61. Rochette, P., and Labossière, P. (2000). “Axial testing of rectangular column models confined with composites.” J. Compos. Constr., 4(3), 129-136. 10
© Copyright 2026 Paperzz