The Modified Compression Field Theory Then and Now
The Modified Compression Field Theory Then and Now
The Modified Compression Field Theory Then and Now
SYMPOSIUM VOLUME
Editors:
Denis Mitchell and Abdeldjelil Belarbi
SP-328
Shear in Structural Concrete
Sponsored by ACI/ASCE Committee 445
Editors:
Denis Mitchell and
Abdeldjelil Belarbi
SP-328
First printing, September 2018
Discussion is welcomed for all materials published in this issue and will appear ten months from this journal’s date
if the discussion is received within four months of the paper’s print publication. Discussion of material received
after specified dates will be considered individually for publication or private response. ACI Standards published
in ACI Journals for public comment have discussion due dates printed with the Standard.
The Institute is not responsible for the statements or opinions expressed in its publications. Institute publications
are not able to, nor intended to, supplant individual training, responsibility, or judgment of the user, or the supplier,
of the information presented.
The papers in this volume have been reviewed under Institute publication procedures by individuals expert in
the subject areas of the papers.
Copyright © 2018
AMERICAN CONCRETE INSTITUTE
38800 Country Club Dr.
Farmington Hills, Michigan 48331
All rights reserved, including rights of reproduction and use in any form or by any means, including the making of
copies by any photo process, or by any electronic or mechanical device, printed or written or oral, or recording
for sound or visual reproduction or for use in any knowledge or retrieval system or device, unless permission in
writing is obtained from the copyright proprietors.
ISBN-13: 978-1-64195-025-1
Preface
Editors:
Denis Mitchell (ACI 445)
Abdeldjelil Belarbi (Chair of ACI 445)
SP-328—3
Abstract: The Modified Compression Field Theory (MCFT) was introduced almost 40 years ago as a simple and
effective model for calculating the response of reinforced concrete elements under general loading conditions with
particular focus on shear. The model was based on a smeared rotating crack concept, and treated cracked reinforced
concrete as a new orthotropic material with unique constitutive relationships. An extension of MCFT, known as the
Disturbed Stress Field Model (DSFM), was later developed which removed some restrictions and increased the
accuracy of the method. The MCFT has been adapted to various types of finite element analysis programs as well as
structural design codes. In recent years, the application of the method has been extended to more advanced research
areas including extreme loading conditions, stochastic analysis, fiber-reinforced concrete, repaired structures, multi-
scale analysis, and hybrid simulation. This paper presents a brief overview of the original formulation and its
evolvement over the last three decades. In addition, the adaptation of the method to advanced research areas are
discussed. It is concluded that the MCFT remains a viable and effective model, whose value lies in its simple yet
versatile formulation which enables it to serve as a foundation for accurately solving many diverse and complex
problems pertaining to reinforced concrete structures.
3.1
Vahid Sadeghian and Frank Vecchio
Biography:
Vahid Sadeghian is a PhD graduate at the Department of Civil Engineering at the University of Toronto. His research
interests include hybrid (experimental-numerical) simulation, multi-platform analysis of reinforced concrete
structures, and performance assessment of repaired structures.
Frank J. Vecchio, FACI, is a Professor at the Department of Civil Engineering at the University of Toronto. He is a
member of Joint ACI-ASCE Committees 441, Reinforced Concrete Columns, and 447, Finite Element Analysis of
Reinforced Concrete Structures. His research interests include nonlinear analysis and design of reinforced concrete
structures, constitutive modeling, performance assessment and forensic investigation, and repair and rehabilitation of
structures.
INTRODUCTION
Lessons Learned from Shear Failures
Predicting the shear behavior of reinforced concrete (RC) structures is a crucial and challenging task which has not
been fully solved yet. Over the last few decades, underestimating the complexity of the problem and relying on limited
knowledge of shear behavior have led to several catastrophic failures. Some of the more recognized shear failures,
and the resulting experimental investigations that greatly improved our understanding of the behavior of concrete in
shear, are summarized in the following.
The basis for shear design provisions of most current codes is a truss model developed by Mörsch (1909) more than
a century ago. In this model, the compression component of the shear load is carried by concrete struts, while the
tension component is carried by longitudinal and transverse reinforcement. In 1955, the shear failure of a US Air Force
warehouse raised questions about the accuracy of such shear design provisions. The experimental investigations that
followed demonstrated the detrimental effect of axial tensile force on the shear strength of lightly reinforced members
(Elstner and Hognestad, 1957). Later, Bhide and Collins (1989) tested 24 panel elements under combined shear and
tension and concluded that the reduction in shear capacity was influenced by the amount and distribution of the
longitudinal reinforcement. For members containing sufficient longitudinal reinforcement, this reduction was much
less than in members with no tension reinforcement. Based on a review of previous experimental studies, Collins and
Kuchma (1999) demonstrated that the failure of the US Air Force warehouse was mainly due to the size effect in
shear, rather than the influence of axial tensile force. One of the most well-known test programs to investigate shear
behavior was conducted by Kani et al. (1979) which included examining the response of more than 300 specimens.
The experimental findings resulted in development of new shear design equations that were included in the ACI code
(1971) and mostly remain unchanged to date. In 1991, the collapse of the Sleipner A offshore platform, the most
expensive shear failure in history with a total economic loss of about $700 million, and the resulting experimental
investigation conducted by Gupta (2001) indicated that the beneficial effects of the axial compression force on the
shear strength had been seriously overestimated and was among the factors contributing to the failure. It was concluded
that the shear design provisions of ACI code for structures with high axial load or prestressed structures can be unsafe.
While the problem of shear behavior was assumed to be nearly solved in recent years, the collapse of the Laval bridge
in 2006, caused by shear failure of a 1.25 m (49.21 in.) thick reinforced concrete slab, showed that it is still too early
to draw such a conclusion. The experimental tests conducted by Sherwood et al. (2007) on shear-critical slab-strip
specimens indicated the importance of accurately considering the size effect and aggregate interlock in predicting
shear behavior. Collins et al. (2008) compared the observed failure loads of these specimens against the failure loads
calculated using Canadian, American, European, and British design codes. There were considerable inconsistencies
between predictions of the four design codes particularly with variation in the specimen depth or addition of stirrups.
Collins et al. (2015) further investigated the size effect by testing a very deep slab specimen under an off-center point
load. To assess the ability of existing analysis and design procedures, prior to the test, 66 engineers were invited to
predict the failure loads of the two shear spans. It was concluded that the failure load of the shear span with no shear
reinforcement was dangerously overestimated by many engineers. The ‘size effect’, and the divergent behavior
between narrow deep elements (beams) versus wide shallow elements (slabs), remain contentious issues.
Analytical Models for Shear Behavior
The aforementioned catastrophic shear failures, and many others, have inspired engineers to develop progressively
more refined analytical models over the years. The Mörsch truss model assumed that the inclination of the concrete
compression struts was constant and equal to 45 degrees, a simplification that was later found to be incorrect. Hence,
modified truss models were proposed with variable angle of inclinations. These models were further improved by
considering the contributions of the concrete to the shear strength in the form of residual tensile strength and aggregate
3.2
The Modified Compression Field Theory: Then and Now
interlock at the crack surface. The application of truss models was extended to non-yielding domains by using the
theory of plasticity. In addition, for beams with large span-to-depth ratio where the direct transfer of force by the
compression strut is not fully valid, a “tooth model” originally proposed by Kani (1964) was found to be a better
representation of the shear behavior. The aforementioned improvements are discussed in detail by the ACI-ASCE
Committee 445 (1999). In recent years, more advanced truss models have been developed capable of analyzing other
types of shear-critical RC members including walls under cyclic loads (Panagiotou, et al., 2012) and columns with the
arch action (Pan and Li, 2013).
The brittle nature of the typical shear failure encouraged many researchers to apply fracture mechanics principles in
modeling the shear strength of RC members. Hillerborg et al. (1976) proposed the well-recognized fictitious crack
model which represented cracks in a discrete manner and was based on a concrete tensile strength versus crack width
relationship. However, the discrete representation of the cracks could not capture the influence of micro-cracking on
the shear behavior. Bazant and Gambarova (1984) proposed the crack band model in which the fracture surface
consisted of a zone of micro-cracks and the size effect was considered in the energy equation. Van Mier (1996) showed
that residual tensile stresses exist within the crack band as a result of a phenomena known as “crack bridging”. In
recent years, the application of methods based on fracture mechanics has been extended to other types of loading
conditions (e.g., seismic and impact) and other concrete materials (e.g., fiber-reinforced concrete). These
advancements are summarized by IA-FraMCoS (2016). However, most fracture mechanics models require
assumptions for the crack path and the governing failure mechanism prior to analysis.
In an alternative approach to calculating shear strength, several empirical equations were derived by curve fitting to a
series of experimental data. Empirical methods are the basis for the shear design provisions of the American and most
European design codes. However, as discussed in detail by Collins et al. (2008), because of the complexity of the
shear behavior and the narrow range of test parameters, the application range of an empirical method is typically
limited. For a comprehensive representation of situations encountered in practice, models with a more rigorous
theoretical basis are required.
One effective rational approach which has attracted considerable interest in recent decades is the Modified
Compression Field Theory (MCFT) (Vecchio and Collins, 1986). The MCFT is the refined version of the Compression
Field Theory (CFT) (Mitchell and Collins, 1974), a smeared rotating crack model originally developed for analysis of
beams in shear and torsion. In both models, the compatibility and equilibrium conditions as well as the stress-strain
relationships are formulated in terms of average values of stresses and strains. The CFT assumes that the cracked
concrete cannot resist any tension and the shear force is carried by a field of diagonal compression. The inclination of
the compression field is determined using an approach similar to that used in the Tension Field Theory (Wagner, 1929)
developed for the post-buckling shear resistance of thin metal girders. The MCFT improved the original model by
accounting for compressive strength reduction in concrete due transverse cracking, concrete tensile strength after the
cracking, and local stress conditions at the crack surface. The MCFT is the basis for the shear design provisions of the
Canadian design code (2014) and the AASHTO LRFD Bridge Design Specifications (2014). Vecchio (2000) extended
the MCFT and proposed the Disturbed Stress Field Model (DSFM) which removed some of the restrictions and
improved the accuracy of the model. Fig. 1 compares the equilibrium, compatibility and constitutive relationships in
the CFT, MCFT, and DSFM models. Bentz et al. (2006) proposed a simplified version of the MCFT that can predict
the shear strength of RC elements in a form suitable for inclusion in code applications. In recent years, the application
of the MCFT has been extended to more advanced research areas including extreme loading conditions, stochastic
analysis, fiber-reinforced concrete, repaired structures, multi-scale analysis, and hybrid simulation.
This paper presents a summary of the MCFT formulation and the adaptation of the model to sectional and finite
element analysis methods. The extension of the model to the DSFM and the resulting improvements, as well as the
findings of three comprehensive verification studies performed on the MCFT and DSFM models are reviewed. In
addition, a comprehensive discussion regarding the application of the model to advanced research areas and related
references is provided.
3.3
Vahid Sadeghian and Frank Vecchio Modified Compression Field Theory (MCFT) –
Compression Field Theory (CFT) Disturbed Stress Field Model (DSFM)
Original Formulation
Average stresses Average stresses Average stresses Local stresses Average stresses Average stresses Local stresses
perpendicular to crack parallel to crack parallel to crack along crack surface perpendicular to crack parallel to crack along crack surface
Equilibrium
(h) (i)
vci ≤ vcimax ; Ɵɛ = Ɵ𝜎
Average strain conditions are similar to CFT Average strains Slip strains Total strains
Slip deformations are not considered (j)
Slip deformations are not considered Ɵɛ ≠ Ɵ𝜎
Concrete in tension Concrete in compression Tension stiffening Compression softening Tension softening Tension stiffening
Constitutive Relationships
Reinforcement response Shear stress on crack Reinforcement response Compression softening Reinforcement response
3.4
The Modified Compression Field Theory: Then and Now
3.5
Vahid Sadeghian and Frank Vecchio
of average stresses and assuming no dowel action from the reinforcement, the equilibrium equations can be written as
follows:
σx = ρx × fsx + fc1 − vcxy cot Ɵc (1)
σy = ρy × fsy + fc1 − vcxy tan Ɵc (2)
fc1 − fc2
τxy = vcxy = (3)
tan Ɵc + cot Ɵc
where ρx and ρy are the smeared reinforcement ratios in the X and Y directions, fc1 and fc2 are the principal tensile and
principal compressive stresses in concrete, and Ɵc is the crack inclination. Fig. 1(a) and Fig. 1(c) show the free body
diagrams for the equilibrium equations.
Local stress conditions — in cracked reinforced concrete, the tension stiffening mechanism allows the concrete
to carry tensile stresses after cracking resulting from bond with the reinforcement. Since the tensile strength of concrete
is zero at crack locations, the average tensile stress must be carried across the cracks by local increases in the
reinforcement stresses. From equilibrium equations, as shown in Fig. 1(d), the local stresses in the reinforcement (fsxcr,
fsycr) at the crack locations can be determined as:
1
fsxcr = (f + fci + vci cot Ɵc + ρsx fsx ) ≤ fyx (4)
ρsx c1
1
fsycr = (f + fci − vci tan Ɵc + ρsy fsy ) ≤ fyy (5)
ρsy c1
where vci and fci are the shear and compressive stresses on the crack surface (see Fig. 1(i)). The shear stress is
transferred across the crack surface based on the friction mechanism, also known as aggregate interlock. The interface
shear stress is calculated from the following equation based on an experimental study conducted by Walraven (1981):
fci2
vci = 0.18vcimax + 1.64fci − 0.82 ≤ vcimax (6)
vcimax
where
√−fc′
vcimax = (7)
24w
0.31 +
a + 16
where a and w are the maximum aggregate size and the average crack width in millimeters (1.0 mm = 25.4 in.) and
fc′ is the concrete cylinder compressive strength in MPa (1.0 MPa = 145.0 psi). The average crack width (w) is
estimated as the product of the principal tensile strain (ɛ1) and the average crack spacing (sƟ):
w = ɛ1 sƟ (8)
The average crack spacing (sƟ) is determined from the nominal crack spacings in the X and Y directions (sx, sy) as
follows:
1
sƟ = (9)
|sin Ɵ𝑐 | |cos Ɵ𝑐 |
+
sx sy
In addition, the concrete average tensile stress (fc1) should be limited by the reserve capacity of the reinforcement. In
the solution procedure of the MCFT, this limit is checked using a relationship between vcimax and fc1 derived by
substituting Eq. (7) into the local equilibrium equations (Eq. (4) and Eq. (5)). For instance, assuming the reinforcement
in the Y direction is the weaker reinforcement, then the following condition should be satisfied:
3.6
The Modified Compression Field Theory: Then and Now
fc1 ≤ vcimax [0.18 + 0.3(1.64 − cot Ɵc )2 ] tan Ɵc + ρsy (fyy − fsy ) (10)
Compatibility Relations
Since perfect bond is assumed between the concrete and reinforcement, compatibility conditions require that the
deformation of the concrete and reinforcement be identical. Therefore, the average strain in the reinforcement should
be equal to that in the concrete in the corresponding direction:
ɛx = ɛcx = ɛsx ; ɛy = ɛcy = ɛsy (11)
where εx and εy are the element total strains, εcx and εcy are the concrete total strains, and εsx and εsy are the
reinforcement total strains in the X and Y directions. If the total strains in the X and Y directions are known, with the
aid of a Mohr’s circle, the total strains in the principal directions (εc1, εc2) and the crack inclination (Ɵc) can be found:
(ɛcx + ɛcy ) 1
ɛc1 , ɛc2 = ± √(ɛcx − ɛcy )2 + γ2cxy (12)
2 2
1 γcxy
Ɵc = tan−1 [ ] (13)
2 ɛcy − ɛcx
It must be remembered that the inclinations of the principal strains and principal stresses are assumed equal in the
MCFT formulation.
Constitutive Relations
Stress-strain relations for concrete — the compressive response of concrete is shown in Fig. 1(o). The
compression softening effect is taken into account by representing the principal compressive stress (fc2) as a function
of both the principal compressive strain (ɛc2) and the principal tensile strain (ɛc1):
ɛc2 ɛc2 2
fc2 = fc2max [2 ( ′ ) − ( ′ ) ] (14)
ɛc ɛc
fc′
fc2max = ɛ ≤ fc′ (15)
0.8 − 0.34 ( c1 ′ )
ɛc
where fc′ and ɛ′c are the concrete cylinder compressive strength and the corresponding strain.
The tensile response of concrete is shown in Fig. 1(n). Prior to cracking, the response is assumed to be linear:
fc1 = Ec × ɛc1 (ɛc1 ≤ ɛcr ) (16)
where Ec is the concrete elastic modulus and ɛcr is the concrete cracking strain.
After cracking, experimental results showed that even at large tensile strain values, significant levels of tensile stress
were carried by the concrete between the crack locations. This observation led to the development of a preliminary
tension stiffening relationship:
fcr
fc1 = (ɛc1 ≥ ɛcr ) (17)
1 + √200ɛc1
where fcr is the concrete cracking stress. Collins et al. (1996) suggested that for larger elements using “500” instead
of the “200” in the denominator is more suitable.
Stress-strain relations for reinforcement — the reinforcement response in compression and tension is
represented by a bilinear relationship (see Fig. 1(q)):
fsx = Esx × εsx ≤ fyx (18)
fsy = Esy × εsy ≤ fyy (19)
where Esx and Esy are the elastic modulus of the reinforcement in the reference directions.
APPLICATION TO BEAMS
Vecchio and Collins (1988) developed an iterative procedure based on the MCFT for analysis of reinforced concrete
beams. In this procedure, the cross section is divided into a number of concrete layers and longitudinal reinforcing bar
layers. Each layer is defined based on its geometry and material properties. The procedure assumes “plane sections
remain plane” which allows the calculation of the longitudinal strain throughout the cross section as follows:
(εb + εt )
εxi = εt − yi (20)
h
where ɛxi is the longitudinal strain of the i layer, ɛt and ɛb are the strains of the top and bottom layers, yi is the distance
th
from top of the section to the ith layer, and h is the section height.
Two types of methods were proposed to determine the shear stress distribution: 1) a rigorous method in which the
shear stresses are calculated from the analysis of two sections of the beam located close to each other, and 2)
approximate methods based on a uniform shear stress distribution or a parabolic shear strain distribution assumption.
3.7
Vahid Sadeghian and Frank Vecchio
Fig. 4 shows the layered beam model with the linear longitudinal strain distribution and the parabolic shear strain
distribution.
Knowing the longitudinal strain and shear stress distributions, and with the aid of the MCFT formulation, the
longitudinal stresses for the concrete layers (fcx) and the reinforcement layers (fsx) can be determined. To satisfy
equilibrium requirements, the resultant forces from the internal stresses must equal to the applied sectional forces.
Therefore:
m n
3.8
The Modified Compression Field Theory: Then and Now
From the compatibility conditions at the crack surface, the local strain of the ith reinforcement component (ɛscri) can
be determined using the average reinforcement strain (ɛsi) and the local incremental strain (Δɛcri).
εscri = εsi + Δεcri × cos 2 (Ɵ − αi ) (25)
The local incremental strains are computed such that the resultant local reinforcement stresses satisfy the equilibrium
equation expressed in Eq. (24). Also, the requirement for the principal tensile stress can be written as:
n
3.9
Vahid Sadeghian and Frank Vecchio
3.10
The Modified Compression Field Theory: Then and Now
where s is the average crack spacing as defined in the MCFT and δs is the slip along the crack surface determined
based on a slip model. The default slip model currently used is the one proposed by Walraven (1981) as follows:
vci
δs = (38)
1.8w −0.8 + (0.234w −0.707 − 0.20)fcc
where fcc is the concrete compressive cube strength, vci is the shear stress along the crack surface determined from Eq.
(27), and w is the average crack width estimated from Eq. (8).
Constructing a Mohr’s circle for the average shear slip strain (γs ), yields the following relations for the strain
components relative to the X and Y directions:
γs γs
ɛsx = − sin 2θ ; ɛsy = sin 2θ ; γsxy = γs cos 2θ (39)
2 2
Fig. 7 shows the average deformations, local deformations due to slip along the crack surface, and total deformations
of an RC element.
Because of the assumption of perfect bond between the concrete and reinforcement, the average strain in the
reinforcement should be computed from the total strains. At the crack surface, the local stresses and strains in the
reinforcement are calculated using a similar procedure to that described for the MCFT in Eq. (24) and Eq. (25).
Constitutive Relations
Stress-strain relations for concrete — to consider the compression softening response of concrete, a reduction
factor (βd) was introduced as the following:
1
βd = ≤ 1.0 (40)
1 + Cs Cd
The Cd factor is a function of the principal compressive strain (ɛc2) and the coexisting principal tensile strain (ɛc1).
Further analysis of the data obtained from panel element tests revealed that considering the crack slip deformations in
the compatibility relations results in a more softened element response, but less softened concrete response, relative
to the MCFT; therefore, the rate of the strength degradation due to the transverse cracking must be reduced. This
reduction is taken into account by the Cs factor. Using Cs = 0.55 produced the best fit to the test data when applying
the DSFM; Cs = 1.0 is implied when using the MCFT.
3.11
Vahid Sadeghian and Frank Vecchio
The βd factor is used to reduce both the peak strength (fp) and the corresponding strain (ɛp) in the concrete compressive
stress-strain relationship as follows (see Fig. 1(t)):
fp = βd fc′ ; ɛp = βd ɛ′c (41)
Similar to the MCFT, the tensile response of concrete before cracking is considered using a linear relation presented
in Eq. (16). After cracking, the DSFM considers two independent mechanisms (tension softening and tension
stiffening) to determine the residual concrete tensile stress. The post-cracking tensile stress resulting from the tension
softening mechanism (fc1a) is determined using a fracture mechanics relation as follows:
a
(ɛc1 − ɛcr )
fc1 = fcr [1 − ] (42)
(ɛts − ɛcr )
where ɛts is the terminal strain calculated from the concrete fracture energy. The tension softening mechanism is
particularly significant in lightly reinforced concrete members.
The post-cracking tensile stress produced from the tension stiffening mechanism (fc1b) is calculated based on a model
proposed by Bentz (2000) which was later modified by Vecchio (2000) to account for reinforcement with arbitrary
directions:
b
fcr
fc1 = (43)
1 + √ct ɛc1
The MCFT uses a constant number for the ct factor (200 for small elements and 500 for large elements). Bentz (2000)
showed that the tension stiffening is also a function of the reinforcement ratio (ρ) and the reinforcing bar diameter
(db), and defined the ct factor accordingly.
Fig. 1(r) and Fig. 1(s) show the tension softening and the tension stiffening mechanisms in the DSFM, respectively.
The average principal tensile stress (fc1) is defined as the larger of the two mechanisms:
a b
fc1 = max(fc1 , fc1 ) (44)
Stress-strain relations for reinforcement — the DSFM uses a stress-strain relationship for the reinforcement
similar to that of the MCFT except it considers strain hardening effects. Prior to strain hardening, the response can be
expressed using the relationship in Eq. (18). After the strain hardening, the following equation is used to determine
the response of the ith reinforcement component in both tension and compression:
fsi = fyi + Eshi (ɛsi − ɛshi ) (ɛshi < ɛsi < ɛui ) (45)
where Eshi is the strain hardening modulus, ɛshi is the strain hardening strain, and ɛui is the ultimate strain. The
reinforcement response is shown in Fig. 1(u).
VERIFICATION STUDIES
Over the last several years, the accuracy of the MCFT and DSFM models have been examined through the analyses
of various types of test specimens (e.g., Vecchio et al., 2001; Palermo and Vecchio, 2007; Hunter 2016) as well as
real-world structures (e.g., Vecchio, 2002; Vecchio et al., 2004). In this section, summaries of three comprehensive
verification studies that were performed on panel element, shear wall, and beam specimens using the VecTor2 program
are presented as examples.
Verification against Panel Element Tests
Vecchio et al. (2001) modeled three series of panel element specimens, a total of 40 specimens, using both the MCFT
and DSFM models and compared the analytical results against the experimental behavior. The specimens varied in
terms of the concrete compressive strength, the longitudinal and transverse reinforcement ratios, and the loading
condition. In general, the analytical and experimental load-deflection responses agreed reasonably well. With the
MCFT model, the computed shear strength had a mean analytical to experimental value of 1.05 with a coefficient of
variation (COV) of 10.3%. With the DSFM model, the load-deflection responses were estimated more accurately than
the MCFT, and resulted in a mean analytical to experimental value of 1.02 with a COV of 9.6% for the shear strength.
The study also concluded that the analytical responses for uniaxially reinforced panels were highly sensitive to the
concrete cracking strength. Since this is an extremely variable parameter, regardless of the analysis procedure,
modeling of such members with a high level of accuracy may not be consistently achieved.
Verification against Shear Wall Tests
Palermo and Vecchio (2007) examined the accuracy of the DSFM model by performing a parametric study on 21
large-scale shear walls ranging from slender to squat with various types of cross sections tested under reversed cyclic
loading conditions. The analytical procedure fairly accurately captured the experimental behavior in terms of the peak
load capacity, ductility, energy dissipation, and failure mode. The ratio of the computed-to-observed shear strength
obtained using VecTor2 had a mean of 1.01 and a COV of 6.1%; for the lateral displacement at ultimate strength, the
corresponding numbers were 1.09 and 36.4%. Clearly, accurately calculating deformation capacity consistently
remains a challenge. Vecchio et al. (2001) conducted a similar analytical study on two series of large-scale shear walls,
3.12
The Modified Compression Field Theory: Then and Now
a total of 13 specimens, tested under monotonic loading conditions. The specimens varied in terms of the level of axial
load, height-to-width ratio, and amount of reinforcement. The results computed by the MCFT and DSFM models were
almost identical. The ratio of the calculated-to-measured strength for all the shear walls had a mean of 1.01 with a
COV of 5.3%. The computed failure mode consisted of concrete crushing in the toe region combined with a sliding
shear failure along the base, agreeing well with the experimentally observed behavior. Overall, the ductility of the
shear walls was computed with reasonable accuracy as well.
Verification against Beam Tests
Hunter (2016) assessed the professional factor (a random variable that captures the uncertainty of a numerical model)
of the VecTor2 program for analysis of shear-critical beams with no transverse reinforcement. The study included
modeling a total of 376 shear-critical beams, which had sufficient testing information, selected from a database
reported by Reineck et al. (2013). The beams had a rectangular cross section with span-to-depth ratios greater than
2.41, and were subjected to a one-point or two-point bending test. Based on the analysis results, the professional factor,
defined as the ratio of the calculated-to-measured peak load, had a mean of 0.94 and a COV of 18.7%. Similar analyses
were performed with the CSA A23.3 and ACI-318 design codes. The CSA A23.3 code whose shear design provisions
are based on a simplified version of the MCFT resulted in acceptable and safe predications for the shear strength of
the majority of the beams. However, the results obtained from the ACI-318 code were highly conservative for beams
with effective depths less than 500 mm and, more importantly, were unconservative for beams with effective depths
greater than 750 mm. This was mainly attributed to neglecting the size effect in the shear design provisions of the
ACI-318 code. Another parametric study on the behavior of shear-critical beams was undertaken by Vecchio (2001).
In this study, three series of beams with little or no transverse reinforcement (a total of 54 specimens) were modeled
with the MCFT and DSFM. The ratio of the calculated-to-observed shear capacity obtained from the MCFT had a
mean of 1.00 with a COV 20.3%. The DSFM resulted in improved analytical responses with the mean of 1.07 and
COV of 16.9%. For shear-critical beams containing no shear reinforcement, these should be considered good results.
Fig. 8 shows the calculated-to-measured load capacity for all the specimens of the aforementioned verification studies.
vcal / vexp
RECENT ADVANCES
In recent years, the application of the MCFT has been extended to more advanced research areas. In addition, pre- and
post-processors have been developed for the VecTor suite of software contributing to the programs’ utility. A summary
of these advancements and the related references are presented in the following.
Pre- and Post-Processors
A generalized user-friendly pre-processor, FormWorks-Plus (Sadeghian and Vecchio, 2015), has been developed to
facilitate the process of creating finite element models for the VecTor suite of programs. FormWorks-Plus is an
extended version of FormWorks (Wong and Vecchio, 2002) which was a graphical user-interface specialized for the
VecTor2 program. The pre-processors provide easy-to-use features for the user to create the finite element mesh and
define materials, loads, and support conditions. They also perform a series of checks to ensure the finite element model
is defined properly and in accordance with the format required for the VecTor programs inputs. FormWorks is
enhanced with a two-dimensional auto-meshing feature, facilitating the mesh generation process. FormWorks-Plus,
3.13
Vahid Sadeghian and Frank Vecchio
which is compatible with all the VecTor programs, provides a wide range of viewing options, enabling modeling
structures with complex geometries.
Most academic analysis software, including the VecTor programs, store the analysis results in a series of text files.
However due to the substantial volume of output data, accessing the relevant analysis information is often a time-
consuming and error-prone process. A graphical post-processor, named Janus (Chak, 2013), has been developed that
allows the user to comprehensively recall and manipulate structural analysis results both at the system-level and
element-level. Janus is capable of displaying different response characteristics of the structure including deformed
shape, stress and strain variations in the concrete and reinforcing bar elements, and various types of graphs. The
program has been recently improved to display the crack propagation and the element sectional view for frame
structures (Loya et al., 2017). Janus is compatible with all the VecTor programs. For VecTor2, the post-processor
program Augustus (Bentz, 2000) is the preeminent tool for visualization of the analysis results.
Application to Advanced Loading Conditions
Dynamic analysis — in recent years, with elevated terror threats around the world, the need for accurate analysis
of RC structures subjected to blast and impact loads has been significantly increased. Test results have shown that RC
structures typically experience a shear-dominant behavior under impact loads. The MCFT, which was successfully
employed for the analysis of shear-critical members under static loads, was further extended to consider dynamic
loads (Saatci and Vecchio, 2009). For this purpose, the original tangent-stiffness-based formulation of Newmark’s
method (Newmark, 1959) was presented in the form of total-load secant-stiffness formulation similar to that used in
the MCFT. The procedure was implemented in VecTor2 and verified through the analysis of a series of beam
specimens subjected to impact loads. The analytical and experimental results agreed well in terms of the failure modes,
displacements, and damage levels. The minor discrepancies were attributed to deficiencies regarding modeling the
strain rate effects. Guner and Vecchio (2012) derived a simplified version of the method for sectional analysis of RC
frames under dynamic loads. They incorporated proper concrete and reinforcement strain rate models into the dynamic
formulation of the MCFT, improving the accuracy of the method. Trommels (2013) performed a numerical study to
assess the analysis capabilities of VecTor2 and VecTor3 for impact and impulse loading conditions. The study
demonstrated the importance of proper modeling of the support conditions and the shape of the pressure-time profile.
Hrynyk (2013) extended the method to three-dimensional analysis of reinforced concrete and fiber-reinforced concrete
shells subjected to high-mass low-velocity impact loads. Lulec (2017) further developed the method for analysis of
three-dimensional and axisymmetric RC solid structures under blast and hard missile impact loading conditions.
Missile perforation velocities were estimated using a semi-analytical relationship based on the MCFT. A series of
verification tests was performed to examine the accuracy of the analytical procedure. Overall, the analytical and
experimental results were in good agreement in terms of the time-displacement response and damage mode.
Fire analysis — advanced modeling of RC structures subjected to fire has been the focus of researchers for many
decades. However, current analytical tools typically have major deficiencies in considering simultaneously both
thermal- and structural-related mechanisms. More importantly, they are typically unable to accurately compute the
shear behavior of RC structures under elevated temperatures. ElMohandes and Vecchio (2016) proposed a time
stepping procedure, developed based on the MCFT, which considers the thermal- and structural-related mechanisms
in a coupled manner. The procedure consists of two levels of analysis. At the first level, a nonlinear heat transfer
analysis is performed at certain time steps to determine the thermal expansion strains and reduction factors for the
temperature-dependent mechanical properties of the concrete and reinforcement. The heat transfer analysis can be
paired with a moisture transfer analysis that enables calculating the additional tensile stresses resulted from the vapor
pressure developed in the concrete pores. At the second level, a nonlinear structural analysis is performed taking into
account the first level analysis results as well as the externally applied loads. The proposed procedure was
implemented in VecTor3 and was successfully verified through coupled thermal and structural analysis of a series of
RC column specimens reported in the literature.
Fatigue analysis — recently, Isojeh et al. (2017a; 2017b; and 2017c) experimentally investigated the behavior of
reinforced concrete and steel-fiber concrete elements subjected to fatigue loading and proposed a set of damage and
constitutive models for the concrete composites (including steel fibers) and a crack-growth model for the
reinforcement. These models were incorporated into the DSFM formulation enabling analysis of concrete structures
under a fatigue loading condition. With this analysis procedure, first a set of material damage values for a certain
number of fatigue loading cycles is evaluated. Then, the structure is analyzed under the externally applied loads taking
into account the fatigue damage effects. The analytical procedure was corroborated by examining the behavior of a
series of small-scale reinforced concrete and steel-fiber concrete deep beams tested under fatigue loads. In general,
the analytical load-deflection responses agreed well with the experimentally observed behavior.
3.14
The Modified Compression Field Theory: Then and Now
3.15
Vahid Sadeghian and Frank Vecchio
Unreinforced masonry — masonry is a composite orthotropic material comprised of blocks connected by mortar
joints oriented in two perpendicular directions (i.e., head and bed joints). Facconi et al. (2014) developed a unique
macro-modeling approach based on the DSFM for the analysis of unreinforced masonry structures. Unlike
conventional smeared crack models, the proposed method is capable of considering the global behavior of the
composite material in an average sense as well as the local shear slip response of the mortar joints. The equilibrium
and compatibility relations are formulated similar to the DSFM. For the equilibrium equations, unlike with reinforced
concrete, the shear stresses along the crack surface are assumed to be zero. For the compatibility equations, the total
slip strain is calculated from the summation of the slip strains of the head and bed joints defined with respect to the X
and Y directions. Also, the average crack spacing is defined as a function of the head and bed joint spacings. The
compressive stress-strain relationship is represented with a concrete model modified according to the masonry
characteristics. A failure criterion for masonry in compression, which takes into account the failures of the components
and the composite material, is adopted. The masonry tensile response is assumed to be linear elastic prior to cracking
and decay gradually after cracking (i.e., tension softening mechanism). The shear slip behavior of the masonry joints
is represent with an elastic-plastic model based on the Mohr-Coulomb yield criterion. Details of the aforementioned
models are provided by Facconi et al. (2014). The analytical procedure was verified against a full-scale masonry wall
and building façade tests reported in the literature. The procedure was able to capture the experimental behavior for
both ductile and brittle failure modes reasonably well.
Application to RC Structures Repaired with FRP Sheets
Fiber-reinforced polymer (FRP) composites have been shown to be an effective repair method to improve the
performance of damaged RC structures. Various types of analytical models have been proposed to predict the behavior
of RC structures repaired with FRP sheets. However, most of the studies have focused on the flexure-critical structures
and the shear behavior of such structures has not been fully understood. Montoya et al. (2004) adopted the MCFT
formulation and developed a constitutive relationship for concrete under the triaxial compression. The model was
corroborated by analysis of several steel- and FRP-confined concrete column specimens. Sato and Vecchio (2003)
developed a set of analytical models to represent the bond characteristics between the concrete and FRP sheets, and
incorporated them into the DSFM formulation. The contributions of the FRP sheets to the crack formation and tension
stiffening effects were formulated according to the equilibrium and compatibility conditions. A parametric study was
performed to reduce the formulations into simplified models, enabling their implementation into the DSFM. The
models were validated through analyses of a series of RC beam specimens containing externally bonded FRP sheets.
Kim and Vecchio (2008) successfully employed the analytical procedure proposed by Sato and Vecchio (2003) for
modeling a shear-critical RC frame repaired with FRP wrap. The bond between the concrete and the FRP wrap was
modeled with two-noded non-dimensional link elements. The damage effects prior to repair were taken into account
by recording the stress and strain history of the elements. Sadeghian and Vecchio (2016) extended the application of
the method to multi-scale modeling, enabling analysis of large repaired RC structures. Details of the proposed multi-
scale modeling approach is presented in the following sections.
Application to Stochastic Analysis
Stochastic simulation is often used to evaluate reliability indices for structural design codes. A reliability index defines
the probability of structural failure in the limit state design method. Hunter (2016) implemented stochastic simulation
techniques into VecTor2 for reliability analysis of shear-critical RC structures. Various types of statistic models for
the concrete and reinforcement were adopted from the literature. The program was extended to perform Monte Carlo
(MC) simulation, Latin Hypercube simulation (LHS), random field generation (both MC and LHS), and correlated
sampling (both MC and LHS). Based on the selected stochastic simulation method, random samples are generated for
the concrete and reinforcement material properties. The random sampled values are then used to modify the user
defined material properties and generate random inputs for the analysis. By repeating the analysis with randomly
generated inputs, the stochastic simulation derives probabilistic information regarding the structure. The application
of the implemented stochastic simulation techniques was demonstrated by reliability analysis of shear-critical RC
beams without transverse reinforcement for the CSA A23.3 (2014) and ACI-318 (2014) design codes. Habibi (2017)
employed the stochastic tools of VecTor2 to investigate the sensitivity of the analysis results to the input parameters
for corroded RC beam specimens. It was concluded that for RC beams damaged by uniform corrosion, modeling the
rate of corrosion as a single variable rather than a random field results in more accurate statistical responses.
Application to Multi-Platform Analysis and Hybrid Simulation
Due to the complex behavior of cracked reinforced concrete, modeling an entire structure in fine detail using a
nonlinear finite element analysis program may not be practical. In addition, situations often arises where a multi-
disciplinary analysis may be warranted (e.g., soil-structure interaction analysis). Sadeghian (2017) developed a
simulation framework, named Cyrus, for multi-platform analysis of RC structures. In this approach, based on the
mechanical characteristics of the members and the anticipated locations of critical regions, the structural system is
3.16
The Modified Compression Field Theory: Then and Now
divided into several substructures. Each substructure can be modeled using a wide range of analysis programs that are
compatible with Cyrus including the VecTor suite of software. Cyrus coordinates the analysis and takes into account
the interaction between the substructures. This allows a user to integrate different types of MCFT and DSFM
formulations with each other (e.g., sectional analysis and finite element analysis) or with other types of analysis
procedures available in structural software. The application of the proposed simulation framework was demonstrated
through multi-platform analysis of several RC structures. In addition, a new interface element, named F2M element
(Sadeghian et al., 2017a), was developed to connect layered beam elements to membrane elements in RC frames. The
F2M element, which is based on the DSFM formulation, computes realistic stress distributions and allows for
transverse expansion at the connection section between the two elements.
The simulation framework was further extended to combine numerical modules with experimental modules to
accommodate hybrid testing (Sadeghian et al., 2017b). The physical modules are incorporated into the framework
using a generalized controller interface program compatible with a wide range of laboratory equipment and testing
configurations. The hybrid simulation capabilities of the framework were verified by conducting a small-scale
experimental program using a six degrees-of-freedom hydraulic testing facility. The hybrid simulations included
integration of the VecTor2 and VecTor5 programs with a small-scale test specimen.
CONCLUSIONS
The MCFT, developed almost 40 years ago, remains an effective general model for analysis of reinforced concrete
structures with detailed consideration of shear-related mechanisms. The MCFT shares the same principles as its
predecessor model, CFT, as well as with its refined version, the DSFM. All three models are based on the smeared
rotating crack concept, treat stresses and strains in an average sense, satisfy equilibrium and compatibility conditions,
and represent cracked reinforced concrete as an orthotropic material with unique constitutive relationships. The
models can be distinguished primarily in two aspects: 1) advancements in the constitutive relationships and 2)
improvements in local considerations at the crack surface. Relative to the CFT, the MCFT takes into account the
compression softening and tension stiffening mechanisms in the concrete constitutive relationships, and introduces
the need to examine local stress conditions at crack locations. The DSFM incorporates further refinements in the
constitutive models also takes into account shear slip along crack surfaces within the element compatibility
formulations, thus allowing for the divergence of principal stress and principal strain directions.
Over the last three decades, the MCFT has been adapted to various types of sectional and finite element analysis
programs as well as structural design codes. In addition, the formulations have been extended to several state-of-the-
art research topics including extreme loading conditions, advanced material types, stochastic analysis, repaired
structures, multi-scale analysis, and hybrid simulation. Many advancements were incorporated into the MCFT over
the years; however, the fundamentals of the model remain unchanged.
The factors that make the MCFT a versatile and enduring model, which can be successfully applied to many diverse
and complex structural problems with a high level of accuracy, can be summarized as follows:
1) The MCFT incorporates a rational and transparent approach based on simple and well-defined equations derived
from fundamental principles of structural analysis and solid mechanics.
2) The nature of its formulation makes it easily adaptable to finite element analysis algorithms and other simple or
advanced analysis procedures.
3) The MCFT provides a framework whereby a diverse range of constitutive models, behavior mechanisms, and
material types can be included, enabling the consideration of a broad range of structures and loading conditions.
4) The MCFT does not rely on assumptions, material characterizations, or analysis calibrations that may vary from
one type of structural element or load condition to another.
5) The smeared crack concept enables a macro-modeling approach, reducing modeling and computational effort.
6) Over the past four decades, the MCFT and its related analysis tools has benefitted from constant evolution and
continuous verification using results from comprehensive well-defined testing programs.
ACKNOWLEDGEMENT
The concept and preliminary formulation of the Modified Compression Field Theory were conceived by Professor
Michael P. Collins, University of Toronto. The vision, guidance and encouragement he provided in its further
development and application were essential to its success.
REFERENCES
AASHTO LRFD Bridge Design Specifications 7th edition, American Association of State Highway and Transportation
Officials, Washington, U.S.A., 2014, 2150 pp.
3.17
Vahid Sadeghian and Frank Vecchio
ACI-ASCE Committee 445, “Recent Approaches to Shear Design of Structural Concrete,” American Concrete
Institute, 1999, 50 pp.
ACI Committee 318, “Building Code Requirements for Reinforced Concrete and Commentary,” American Concrete
Institute, Detroit, U.S.A., 2014, 520 pp.
ACI Committee 318, “Building Code Requirements for Reinforced Concrete,” American Concrete Institute, Detroit,
U.S.A., 1971, 78 pp.
Bazant, Z. P., and Gambarova, P., “Crack Shear in Concrete: Crack Band Microplane Model,” ASCE Journal of
Engineering Mechanics, V. 110, 1984, pp. 2015-2035.
Bentz, E. C., “Sectional Analysis of Reinforced Concrete Members,” PhD Thesis, Department of Civil Engineering,
University of Toronto, 2000, 310 pp.
Bentz, E. C., Vecchio, F. J., and Collins, M. P., “The Simplified MCFT for Calculating the Shear Strength of
Reinforced Concrete Elements,” ACI Structural Journal, V.103, No. 4, 2006, pp. 614-624.
Bhide, S., and Collins, M. P., “Influence of Axial Tension on the Shear Capacity of Reinforced Concrete Members,”
ACI Structural Journal, V. 86, No. 5, 1989, pp. 570-581.
Chak, I. N., “Janus: A Post-Processor for VecTor Analysis Software,” M.A.Sc. Thesis, Department of Civil
Engineering, University of Toronto, Canada, 2013, 193 pp.
Collins, M. P., Bentz, E. C., Quach, P., Proestos, G., “The Challenge of Predicting the Shear Strength of Very Thick
Slabs,” Concrete International, V. 37, No. 11, 2015, pp. 29-37.
Collins, M. P., and Kuchma, D., “How Safe Are Our Large, Lightly Reinforced Concrete Beams, Slabs and Footings?”
ACI Structural Journal, V. 96, No. 4, 1999, pp. 482-490.
Collins, M. P., Mitchell, D., Adebar, P., and Vecchio, F. J., “A General Shear Design Method,” ACI Structural
Journal, V. 93, No. 1, 1996, pp. 36-45.
Collins, M. P., Mitchell, D., and Bentz, E. C., “Shear Design of Concrete Structures,” The Structural Engineer, V. 86,
No. 10, 2008, pp. 32-39.
Collins, M. P., Vecchio, F. J., and Mehlhorn, G., “An International Competition to Predict the Response of Reinforced
Concrete Panels,” Canadian Journal of Civil Engineering, V. 12, No. 3, 1985, pp. 624-644.
CSA Committee A23.3, “Design of Concrete Structures,” Canadian Standard Association, Rexdale, Ontario, Canada,
2014, 297 pp.
Deluce, J. R., Lee, S. -C., and Vecchio, F. J., “Crack Model for Steel Fiber-Reinforced Concrete Members Containing
Conventional Reinforcement,” ACI Structural Journal, V. 111, No. 1, 2014, pp. 93-102.
ElMohandes, F., and Vecchio, F. J., “Reliability of Temperature-Dependent Models for Analysis of Reinforced
Concrete Members Subjected to Fire,” ACI Structural Journal, V. 113, No. 3, 2016, pp. 481-490.
ElMohandes, F., and Vecchio, F. J., “VecTor3 User’s Manual,” Department of Civil Engineering, University of
Toronto, Toronto, Canada, 2013, 153 pp.
Elstner, R. C., and Hognestad, E., “Laboratory Investigation of Rigid Frame Failure,” ACI Journal, V. 53, No. 1, 1957,
pp. 637-668.
Facconi, L., Plizzari, G., and Vecchio, F. J., “Disturbed Stress Field Model for Unreinforced Masonry,” ASCE Journal
of Structural Engineering, V. 140, No. 4, 2014.
Ferche, A. C., Panesar, D. K., Sheikh, S. A., and Vecchio, F. J., “Toward Macro-Modeling of ASR-Affected
Structures,” ACI Structural Journal, V. 114, No. 5, 2017, pp. 1121-1129.
Guner, S., and Vecchio, F. J., “Simplified Method for Nonlinear Dynamic Analysis of Shear-Critical Frames,” ACI
Structural Journal, V. 109, No. 5, 2012, pp. 727-737.
Guner, S., and Vecchio, F. J., “Pushover Analysis of Shear-Critical Frames: Formulations,” ACI Structural Journal,
V. 107, No. 1, 2010, pp. 63-71.
Gupta, P. R., and Collins, M. P., “Evaluation of Shear Design Procedures for Reinforced Concrete Members under
Axial Compression,” ACI Structural Journal, V. 98, No. 4, 2001, pp. 537-547.
Habibi, S., “Finite Element Modeling of Corrosion Damaged Reinforced Concrete Structures,” M.A.Sc. Thesis,
Department of Civil Engineering, University of Toronto, Canada, 2017, 139 pp.
Hillerborg, A., Modéer, M., and Petersson, P-E., “Analysis of Crack Formation and Crack Growth in Concrete by
Means of Fracture Mechanics and Finite Elements,” Cement and Concrete Research, V. 6, No. 3, 1976, pp. 773-781.
Hrynyk, T. D., “Behavior and Modeling of Reinforced Concrete Slabs and Shells under Static and Dynamic Loads,”
Ph.D. Thesis, Department of Civil Engineering, University of Toronto, Canada, 2013, 455 pp.
Hrynyk, T. D., and Vecchio, F. J., “Modeling of Steel-Concrete Composite Elements under In-Plane and Out-of-Plane
Loads,” ASCE Journal of Structural Engineering, V. 142, No. 10, 2016.
Hrynyk, T. D., and Vecchio, F. J., “Capturing Out-of-Plane Shear Failures in the Analysis of Reinforced Concrete
Shells,” ASCE Journal of Structural Engineering, V. 141, No. 12, 2015.
3.18
The Modified Compression Field Theory: Then and Now
Hunter, M. D., “Towards Stochastic Finite Element Analysis of Reinforced Concrete Structures,” M.A.Sc. Thesis,
Department of Civil Engineering, University of Toronto, Canada, 2016, 262 pp.
IA-FraMCoS Conference Proceedings, Fracture Mechanics of Concrete and Concrete Structures International
Conference, University of California, Berkeley, U.S.A., 2016.
Isojeh, B., El-Zeghayar, M., and Vecchio, F. J., “Fatigue Behavior of Steel Fiber Concrete in Direct Tension,” ASCE
Journal of Materials in Civil Engineering, V. 29, No. 9, 2017a.
Isojeh, B., El-Zeghayar, M., and Vecchio, F. J., “Concrete Damage under Fatigue Loading in Uniaxial Compression,”
ACI Materials Journal, V. 114, No. 2, 2017b, pp. 225-235.
Isojeh, B., El-Zeghayar, M., and Vecchio, F. J., “Simplified Constitutive Model for Fatigue Behavior of Concrete in
Compression,” ASCE Journal of Materials in Civil Engineering, V. 29, No. 7, 2017c.
Kani, G. N. J., “The Riddle of Shear Failure and Its Solution,” ACI Journal, V. 61, No. 28, 1964, pp. 441-467.
Kani, M. W., Huggins, M. W., and Wittkopp, R. R., “Kani on Shear in Reinforced Concrete,” Department of Civil
Engineering, University of Toronto, Canada, 1979, 225 pp.
Kim, S. W., and Vecchio, F. J., “Modeling of Shear-Critical Reinforced Concrete Structures Repaired with Fiber-
Reinforced Polymer Composites,” ASCE Journal of Structural Engineering, V. 134, No. 8, 2008, pp. 1288-1299.
Lee, S. -C., Cho, J. -Y., and Vecchio, F. J., “Analysis of Steel Fiber-Reinforced Concrete Elements Subjected to
Shear,” ACI Structural Journal, V. 113, No. 2, 2016, pp. 275-285.
Lee, S. -C., Cho, J. -Y., and Vecchio, F. J., “Simplified Diverse Embedment Model for Steel Fiber-Reinforced
Concrete Elements in Tension,” ACI Materials Journal, V. 110, No. 4, 2013a, pp. 403-412.
Lee, S. -C., Cho, J. -Y., and Vecchio, F. J., “Tension-Stiffening Model for Steel Fiber-Reinforced Concrete Containing
Conventional Reinforcement,” ACI Structural Journal, V. 110, No. 4, 2013b, pp. 639-648.
Lee, S. -C., Cho, J. -Y., and Vecchio, F. J., “Diverse Embedment Model for Steel Fiber Reinforced Concrete in
Tension: Model Development,” ACI Materials Journal, V. 108, No. 5, 2011, pp. 516-525.
Loya, A. S., Lourenco, D. D. S., Guner, S., and Vecchio, F. J., “User's Manual of Janus for VecTor5,” Department of
Civil and Environmental Engineering, University of Toledo, Ohio, U.S.A., 2017, 28 pp.
Lulec, A., “Simplified Analytical Tools for Impact and Impulsive Loading Analysis of Reinforced Concrete
Structures,” Ph.D. Thesis, Department of Civil Engineering, University of Toronto, Canada, 2017, 196 pp.
Mitchell, D., and Collins, M. P., “Diagonal Compression Field Theory-A Rational Model for Structural Concrete in
Pure Torsion,” ACI Journal Proceedings, V. 71, No. 8, 1974, pp. 396-408.
Montoya, E., “Behavior and Analysis of Confined Concrete,” Ph.D. Thesis, Department of Civil Engineering,
University of Toronto, Canada, 2003, 321 pp.
Montoya, E., Vecchio, F. J., and Sheikh, S. A., “Compression Field Modeling of Confined Concrete,” Structural
Engineering and Mechanics, V. 12, No. 3, 2006, pp. 231-248.
Montoya, E., Vecchio, F. J., and Sheikh, S. A., “Numerical Evaluation of the Behavior of Steel- and FRP-Confined
Concrete Using Compression Field Modeling,” J. of Engineering Structures, V. 26, No. 11, 2004, pp. 1535-1545.
Mörsch, E., “Concrete-Steel Construction,” McGraw-Hill Book Company, New York, U.S.A., 1909, 368 pp.
Newmark, N. M., “A Method of Computation for Structural Dynamics,” ASCE Journal of the Engineering Mechanics
Division, V. 85, No. 3, 1959, pp. 67-94.
Palermo, D., and Vecchio, F. J., “Simulation of Cyclically Loaded Concrete Structures Based on the Finite Element
Method,” ASCE Journal of Structural Engineering, V. 133, No. 5, 2007, pp. 728-738.
Pan, Z., and Li, B., “Truss-Arch Model for Shear Strength of Shear-Critical Reinforced Concrete Columns,” ASCE
Journal of Structural Engineering, V. 139, No. 4, 2013, pp. 548-560.
Panagiotou, M., Restrepo, J. I., Schoettler, M., and Kim, G., “Nonlinear Cyclic Truss Model for Reinforced Concrete
Walls,” ACI Structural Journal, V. 109, No. 2, 2012, pp. 205-214.
Polak, M. A., and Vecchio, F. J., “Nonlinear Analysis of Reinforced Concrete Shells,” ASCE Journal of Structural
Engineering, V. 119, No. 12, 1993, pp. 3439-3462.
Reineck, K., Bentz, E. C., Fitik, B. F., Kuchma, D. A., and Bayrak, O., “ACI-DAfStb Database of Shear Tests on
Slender Reinforced Concrete Beams without Stirrups,” ACI Structural Journal, V. 110, No. 5, 2013, pp. 867-876.
Saatci, S., and Vecchio, F. J., “Nonlinear Finite Element Modeling of Reinforced Concrete Structures under Impact
Loads,” ACI Structural Journal, V. 106, No. 5, 2009, pp. 717-725.
Sadeghian, V., “A Framework for Multi-Platform Analytical and Experimental Simulations of Reinforced Concrete
Structures,” Ph.D. Thesis, Department of Civil Engineering, University of Toronto, Canada, 2017, 321 pp.
Sadeghian, V., Kwon, O., and Vecchio, F. J., “Modeling Beam-Membrane Interface in Reinforced Concrete Frames,”
ACI Structural Journal, Accepted for publication, 2017a.
Sadeghian, V., Kwon, O., and Vecchio, F. J., “Small-Scale Multi-axial Hybrid Simulation of a Shear-Critical RC
Frame,” Earthquake Engineering and Engineering Vibration, V. 16, No. 4, 2017b, pp. 727-743.
3.19
Vahid Sadeghian and Frank Vecchio
Sadeghian, V., and Vecchio, F. J., “Application of Multi-Scale Modeling on Large Shear-Critical Reinforced Concrete
Structural Systems Repaired with FRP Sheets,” Innovations in Corrosion and Materials Science, V. 6, No. 2, 2016,
pp. 106-114.
Sadeghian, V., and Vecchio, F. J., “A Graphical User Interface for Stand-Alone and Mixed-Type Modeling of
Reinforced Concrete Structures,” Computers and Concrete, V. 16, No. 2, 2015, pp. 287-309.
Sato, Y., and Vecchio, F. J., “Tension Stiffening and Crack Formation in RC Members with FRP Sheets,” ASCE
Journal of Structural Engineering, V. 129, No. 6, 2003, pp.717-724.
Sherwood, E. G., Bentz, E. C., and Collins, M. P., “Effect of Aggregate Size on Beam-Shear Strength of Thick Slabs,”
ACI Structural Journal, V. 104, No. 2, 2007, pp. 180-190.
Susetyo, J., Gauvreau, P., and Vecchio, F. J., “Steel Fiber-Reinforced Concrete Panels in Shear: Analysis and
Modeling,” ACI Structural Journal, V. 110, No. 2, 2013, pp. 285-295.
Trommels, H., “Towards Simplified Tools for Analysis of Reinforced Concrete Structures Subjected to Impact and
Impulsive Loading: A Preliminary Investigation,” M.A.Sc. Thesis, Department of Civil Engineering, University of
Toronto, Canada, 2013, 196 pp.
Van Mier, J. G. M., “Fracture Processes in Concrete,” CRC Press, U.S.A., 1996, 464 pp.
Vecchio, F. J., “Contribution of NLFEA to the Evaluation of Two Structural Concrete Failures,” ASCE Journal of
Performance of Constructed Facilities, V. 16, No. 3, 2002, pp. 110-115.
Vecchio, F. J., “Disturbed Stress Field Model for Reinforced Concrete: Formulation”, ASCE Journal of Structural
Engineering, V. 126, No. 8, 2000, pp. 1070-1077.
Vecchio, F. J., “Finite Element Modeling of Concrete Expansion and Confinement,” ASCE Journal of Structural
Engineering, V. 118, No. 9, 1992, pp. 46-56.
Vecchio, F. J., “Nonlinear Finite Element Analysis of Reinforced Concrete Membranes,” ACI Structural Journal, V.
86, No. 1, 1989, pp. 26-35.
Vecchio, F. J., Bentz, E. C., and Collins, M. P., “Tools for Forensic Analysis of Concrete Structures,” Computer and
Concrete, V. 1, No. 1, 2004, pp 1-14.
Vecchio, F. J., and Collins, M. P., “Predicting the Response of Reinforced Concrete Beams Subjected to Shear using
the Modified Compression Field Theory,” ACI Structural Journal, V. 85, No. 3, 1988, pp. 258-268.
Vecchio, F. J., and Collins, M. P., “The Modified Compression-Field Theory for Reinforced Concrete Elements
Subjected to Shear,” ACI Structural Journal, V. 83, No. 2, 1986, pp. 219-231.
Vecchio, F. J., and Collins, M. P., “Response of Reinforced Concrete to In-Plane Shear and Normal Stresses,”
Department of Civil Engineering, University of Toronto, Canada, 1982, 332 pp.
Vecchio, F. J., Lai, D., Shim, W., and Ng, J., “Disturbed Stress Field Model for Reinforced Concrete: Validation,”
ASCE Journal of Structural Engineering, V. 127, No. 4, 2001, pp. 350-357.
Vecchio, F. J., and McQuade, I., “Towards Improved Modeling of Steel-Concrete Composite Wall Elements,” Journal
of Nuclear Engineering and Design, V. 241, No. 8, 2011, pp. 2629-2642.
Vecchio, F. J., and Selby, R. G., “Towards Compression Field Analysis of Reinforced Concrete Solids,” ASCE
Journal of Structural Engineering, V. 117, No. 6, 1991, pp. 1740-1758.
Wagner, H., “Ebene Blechwandtra¨ger mit sehr du¨nnem Stegblech,” Z. Flugtech. Motorluftschiffahrt, Berlin,
Germany, V. 20, No. 8-12, 1929.
Walraven, J. C., “Fundamental Analysis of Aggregate Interlock,” ASCE Journal of the Structural Division, V. 107,
1981, pp. 2245-2270.
Wong, P. S., Trommels, H., and Vecchio, F. J., “VecTor2 and FormWorks User’s Manual, 2nd Edition,” Department
of Civil Engineering, University of Toronto, Canada, 2013, 311 pp.
Wong, P. S., and Vecchio, F. J., “VecTor2 and FormWorks User’s Manual,” Department of Civil Engineering,
University of Toronto, Canada, 2002, 217 pp.
3.20