Diss Piscitelli LorenzoRuggero
Diss Piscitelli LorenzoRuggero
Diss Piscitelli LorenzoRuggero
Dissertation
submitted to the
and the
by
Lorenzo Ruggero Piscitelli
born 8 March 1986
from Bagno a Ripoli, Italy
2018
*) Either the German or the Italian form of the title may be used.
Dottorato di ricerca in Civil and Environmental Engineering
Ciclo XXX, novembre 2014 - Ottobre 2017
Settore Scientifico Disciplinare ICAR/09
Tutor:
Prof. Maurizio Orlando
Prof. Harald Kloft
Prof. Giovanna Ranocchiai
Coordinatore del dottorato:
Prof. Fabio Castelli
Acknowledgements:
These pages are the result of a long time combined effort of many individuals.
Perhaps, some of them are unaware of the significance and impact of their
contributions on my work and education as a whole; it is now a privilege for me to
release my feelings of profound gratitude to all of them.
To Prof. Orlando, “Maurizio”, as you invited me to call you since day one, for the
abiding support, advices and guidance. You’ve been the driver of my academic
experience since 2009, I could never have wished for a better mentor.
To Prof. Giovanna Ranocchiai, for all the essential contributions to this work, for
the time you always found to help me when I needed the most and for your
unfailing assistance.
To Prof. Harald Kloft for the interest you showed immediately on this topic and
your support during the many months spent at the ITE in Braunschweig.
To Marco Casucci for being the best co-worker a Ph. D student could ever dream
for: curious, proactive, exceptionally resourceful and totally dedicated. I hope our
paths may cross again very soon.
To Saverio Giordano, for sharing a lifetime of experience over the last… many
years! I will try to treasure and put to good use everything that you left me.
To Enzo Barlacchi and Franco Bruni, for always going the extra mile for me; I look
forward to working again side by side with the both of you.
To the research team in Milano: Prof. Luigi Biolzi for your precious guidance; to
Prof. Sara Cattaneo for having me constantly involved in absolutely anything
exciting taking place and for going above and beyond in supporting my work; to
Gabriele D’Incertopadre and Luca Marini, for the precious help throughout many
months of tests and data analysis; to Nicola Piazzolla and Daniele Spinelli for the
joyful moments and their priceless assistance in the everyday life of the laboratory.
Also, my work would not have been possible without the precious contributions
and support of many colleagues: Antonino “Ninni” Marra, Tommaso Massai,
Andrea Giachetti, Enzo Marino, Luca Salvatori, Alessandro Giusti and Luca
Pigolotti. There is an inestimable value in knowing that you can always find the
right answer at the right time, I am grateful for all your help and suggestions.
Thank you my dearest friends, you know who you are, I’m blessed to count you on
both hands. Thank you to my family, thank you “Nonna”, for teaching me the
respect for other people’s work, the value of collaboration over competition, the
importance of taking the time for a break when things get tense, for sitting together
to share a meal… and five more minutes for a coffee.
i
ii
“What do we mean by “understanding” something?
We can imagine that this complicated array of moving
things which constitutes “the world” is something like a
great chess game being played by the gods, and we are
observers of the game. We do not know what the rules of the
game are; all we are allowed to do is to watch the playing.
Of course, if we watch long enough, we may eventually
catch on to a few of the rules. The rules of the game are
what we mean by fundamental physics. Even if we knew
every rule, however, we might not be able to understand
why a particular move is made in the game, merely because
it is too complicated and our minds are limited. If you play
chess you must know that it is easy to learn all the rules,
and yet it is often very hard to select the best move or to
understand why a player moves as he does.
So it is in nature, only much more so.” 1
1
Basic Physics, volume I; lecture 2, Introduction published by Jeffrey Robbins
in "The Value of Science", part of the National Academy of Sciences.
2
Richard Phillips Feynman, New York on the 11 May 1918, Los Angeles 15
February 1988, American physicist, Nobel prize in physics 1965. To the
exceptional skills in many physics and science fields, he combined a sense of
humour out of the ordinary, the passion for music and arts. He defined himself
as a “Nobelist Physicist, teacher, storyteller, bongo player".
iii
INDEX
Index iv
Abstract vi
Workflow of the thesis vii
List of Acronyms ix
3 Mechanics 43
3.1 Elasticity 43
3.1.1 Linear elasticity 44
3.1.2 Non-linear elasticity (hyperelasticity) 44
3.1.3 Time-dependent linear elasticity (linear viscoelasticity) 50
3.2 Plasticity 58
3.2.1 Viscoplasticity 59
3.3 Adhesive bonds 60
4 Experimental tests 63
4.1 Short-term tensile properties of interlayers 63
4.1.1 Design and aim of tensile tests 64
4.1.2 Experimental results 65
4.2 Long-term viscoelastic properties of interlayers 72
4.2.1 Aim and design of tests 72
4.2.2 Preliminary tests and highlights 76
4.2.3 Performed tests 84
4.2.4 Experimental results and first data analysis 86
4.3 Progressive damage in Laminated Glass beams 93
4.3.1 Purpose of tests 94
4.3.2 Design and setup 95
4.3.3 Experimental results 103
4.3.4 Uniaxial tests on fully damaged LG specimens 113
iv
5 Critical analysis and discussion 117
5.1 Interlayer properties 117
5.1.1 Short-term properties of interlayers 117
5.1.2 Viscoelastic properties of interlayers in LG 131
5.1.3 Temperature-depentent adhesive properties 151
5.2 Post-failure properties of LG 154
5.2.1 Analysis of the tension-stiffening effect 154
5.2.2 Dynamic tests 155
5.2.3 Equivalent thermal expansion model 158
5.2.4 Fully damaged LG models 162
5.2.5 Delamination and TS effect 166
5.2.6 Synthesis of post failure LG properties 167
6 Applications 169
6.1 Cold-bent glazing 169
6.1.1 Design hypotheses and material properties 170
6.1.2 Design procedure and verifications 172
6.1.3 Results, design options and interlayers comparison 175
6.2 Post-failure analysis of a LG element 177
6.2.1 Data, constraints and parameters 178
6.2.2 Validation of the numerical model 180
6.2.3 Verifications and design 181
Appendix A 195
Appendix B 199
Appendix C 202
Appendix D 206
Appendix E 207
Appendix F 208
Appendix G 209
Appendix H 214
Appendix I 215
References 222
v
ABSTRACT
Structural laminated glass elements are being used ever more frequently in the
construction industry, following a growing architectural trend that looks for light and
transparency. Nevertheless, an analysis of both regulatory and scientific state-of-the-art
reveals several fields of inquiry which could benefit from deeper investigations. Namely,
properties of plastics used as interlayer materials within the glass plies are scarcely
investigated, professionals being far from unanimous on reliable techniques for
comparing different materials on the same grounds. Yet, such knowledge is needed for
reliable designs, especially in structural applications.
This manuscript presents the results of a multi-scale experimental research on the
mechanical response of three interlayers: PVB, SG and DG41. The former has been the
standard in glass lamination industry, while the latter two are more recent and
supposedly more performing from a mechanical point of view. The hyperelastic
behaviour is studied with simple tensile tests on interlayer specimens; in the end, a novel
generalized response model is proposed, which can be tuned to replicate the complex
short-term and finite-strain response of any thermoplastic using few coefficients. Also,
viscoelastic parameters of interlayers play an essential role in the global long-term
laminated glass elements response. The temperature-dependent viscoelastic problem is
investigated on a larger scale, using double-lap laminated glass joints under compressive
loadings. Tests were performed on specimens made of three glass plies under long-term
imposed actions in a temperature range between and . An existing procedure
was further developed, to provide insight on both creep and relaxation properties.
Finally, calibrated Prony series for viscoelastic models are provided together with
Williams Landel Ferry coefficients for time-temperature superposition, allowing to model
the viscoelastic responses of the three interlayer materials at arbitrary temperatures.
Limits and reliability of such models are discussed; simplified and ready-to use tables are
provided. An analysis on the correlation between mechanical actions and loss of
adhesion is performed.
The third-level of the experimental analyses investigates the mechanical behaviour of
progressively damaged, full-scale laminated glass beams. Risk analyses followed by fail-
safe designs are often mentioned by standards and technical documents, but few studies
have been carried out on the post-failure performance and effects of laminated beams
made with tempered glass plies. After partial failure, the load-bearing capacity depends
on the interlayers ability to generate coupling effects among fractured and undamaged
glass elements through adhesion and its own mechanical properties. Results from
dynamic and static tests are compared and the tension stiffening effect of interlayers in
partially damaged elements is investigated. An equivalent thermal expansion is proposed
to model the mutual effects which are generated among broken and intact plies in the
event of partial failure. Further tests have been carried out in fully damaged conditions,
to evaluate the residual load-bearing capacity and the effects of ageing such conditions.
Uniaxial tensile and compressive test on fully damaged laminated glass specimens are
performed. Results are used to model the response of fully damaged beams.
In the last part of the manuscript, examples of application of newly found results are
used in possible laminated glass structural designs: applications are provided for cold-
bending techniques and post-failure safety assessments.
vi
WORKFLOW OF THE THESIS
CHAPTER · I CHAPTER · II
Structural LG Introduction State-Of-The-Art
Discussion on fields of applications Regulatory and scientific review regarding
for LG with an introduction to properties of interlayers, models and
main components and features theories used for LG design and
advancements on assessing post-failure
performance
CHAPTER · IV
Experimental Tests
Short-term tensile properties of
1. interlayers and strain rate influence.
CHAPTER · V
Critical Analysis And Discussion
Results gathered with experimental tests are discussed in light of existing
mathematical formulations. The scope of empirical evidence is extended
with new or adapted models for interlayer materials and study is provided
on post-failure effects and performance of damaged LG elements.
vii
viii
LIST OF ACRONYMS
ix
x
A mio padre,
xi
Chapter 1 – Structural Glass Introduction
The use of glass in the construction industry has been a part of human
culture from centuries. Nonetheless, despite numerous advantages, glass and
architecture inherently have to deal with compromises. While glass is stiff, it
is brittle, although beautiful, it is costly. This first chapter introduces to the
topic of architectural glass and issues related to glass design, contextualizing
the aim of this manuscript and purpose of the presented researches.
“Don't build a glass house if you're worried about saving money on heating.”
~Philip Johnson
Throughout history of building and construction, glass has always been a material
fascinating both designers and the collective imagination. The research of light and
transparency seems to be inborn in human kind and is often associated to healthy
environments and high quality of life [1–3]. In fact, transparent or translucent
elements in buildings are key for connecting interiors to natural lighting, allowing
for natural environments heating through infrared radiation [4,5].
While in past centuries glass has been used in relatively small portions of buildings
facades, gradual technical evolutions have allowed transparent elements to
dominate most of the outer “skin” of buildings. This continuing transformation
eventually redefined our very perception of many recent large structures with a
feeling lightness, harmoniousness and environmental compatibility (Figure 1.1).
Even so, while detailed properties of glass will be discussed in §1.1.1, it is well-
known that brittleness and unpredictable failure are inherent features in the nature
of this material. To allow for a safe use of glass in buildings (that has to do both
with the perception of safety and structural safety itself), committing it to
increasingly more demanding structural tasks, designers have to deal this
fundamental Achilles' heel.
Figure 1.1 - a) North Rose of Notre Dame (Paris), 1270, still considered one of the largest
windows in the world with its 13.1 m diameter b) The Crystal Palace (London), 1850,
largest a cast-iron and plate-glass structure ever built c) International Forum glass canopy
(Tokyo), 1997, among the first structural applications of laminated glass.
The invention of Laminated Glass (LG) by the beginning of the last century, with a
French Patent filed in France in 1902 1 , was an essential milestone in glass
1
Patent 321.651, filed in May the 31st 1902 by the “société anonyme Le Carbone”, describing the
process of applying a layer of celluloid on a glass object, through immersion and drying, to reduce
the fragility of glass and shattering of fragments uprising from failure.
1/237
Chapter 1 – Structural Glass Introduction
architecture: a turning point in shifting the role of glass elements from carried to
self-sustained and eventually to structural ones [6]. In a nutshell, the lamination
process consists in coupling glass plies to polymers which are generally –but not
always- required transparency and must have prominent adhesive properties with
the surface of glass. The main task of the polymer is to secure glass fragments in
case of glass breakage, avoiding the scattering of potentially harmful shards.
Beyond to this original intent of mere fragments-retaining property, the concept
appeared to provide containing capacity to a glazing element, preventing to some
extent the perforation of colliding bodies [7]. Moreover, assembling many glass
plies in such arrays soon highlighted the structural coupling ability of interlayers,
namely that the resistance and stiffness of LG elements is greater than the sum of
the resistance and stiffness of the glass plies comprising this composite material.
While these features were observed ever since the invention of LG, nowadays,
more than a century later, the fundamental understanding and modelling of these
effects is still an ongoing research field.
Architectural LG usage has progressively extended from secondary elements (e.g.
windows and curtain walls) to structural elements with significant load-bearing
capacity. Nowadays applications of structural LG include floors, stairs, balconies,
canopies, point glazing systems, curtain walls, beams, pillars, sloped and overhead
glazing and essentially any glazing application in which glass is required to be safe
in case of failure and to be able to bear a substantial amount of load.
The architectural research for transparency has historically set the high road for
external glazing and see-through façades, focusing on connecting the interior to the
external environment. Figure 1.2 shows the evolution of the concept of seamless
transparent façades and environmental connection in twenty years, interpreted by
the design engineering firm Rice Francis Ritchie (RFR). In many recent buildings,
with increasingly higher knowledge and confidence on LG resources, planar glass
surfaces are progressively losing ground to curved, sinuous, shapes, aiming to
impart a sense of delicateness and harmony.
Figure 1.2 - a) & b) Cité des Sciences et de l'Industrie greenhouses and reception, Parc de
la Villette, Paris, RFR 1986 c) new TGV railway station in Strasburg, RFR-SNCF 2007
This increasing awareness and reliance on LG properties is helping architects and
engineers to bridge the gap from viewing glass merely as an external interface, to a
proper construction material. Numerous examples of practical LG applications can
nowadays be found in staircases, internal walls, balustrades and walkable areas.
Taking to extremes this general tendency are concepts like the one shown in Figure
1.3 from santambrogiomilano architects, where all traditional structural elements
have been replaced with LG counterparts, designed allow the occupant to be
entirely immersed in the natural scenery. Nonetheless, apart from these
architectural virtuosity exercises and touristic niches, arguably the housing market
not nearly close enough to wild transitions like the ones showcased in such
concepts. Besides, offices and public buildings often chose to actively rely on
2/237
Chapter 1 – Structural Glass Introduction
Figure 1.3 - “The Snow House” concept by Santambrogio & Arosio a) exterior b) interior
Over the past decades, an increase in the demand of structural load bearing glass
elements has occurred [8], stimulating the formation of an organic institutional and
regulatory framework. Transparent glass structures are required serviceability
conditions, granting that the failure of any glass element does not turn out into an
unexpected collapse of the structure and that sufficient load bearing capacity is
available until substitution or evacuation (fail-safe response). From a
standardization point of view, Europe is long time awaiting the approval of a
dedicated Eurocode 3 on glass products and design, which is expected to be a
cornerstone in glass industry. Aiming to fill this regulations void, provisional
norms have been issued, namely prEN13474 4 and later prEN16613 5 , ultimately
dividing the LG market in two macro-sectors: laminated glass and laminated
2
2015, Michael Abrams, Photographer and journalist for Stars and Stripes (www.stripes.com).
3
A series of 10 Standards on different topics, providing a common approach for the design of
buildings and other civil engineering works and construction products. For issues that are
contained within its regulatory framework, Eurocodes are the recommended means of giving a
presumption of conformity with the basic requirements of the Construction Products Regulation
for construction works and products that bear the CE Marking, as well as the preferred reference
for technical specifications in public contracts.
4
prEN13474 “Glass in buildings, Design of glass panes” was the first European attempt to supply
the glass professional and designers with general principles on the design of structural glass
elements. The document was divided in three parts:
prEN13474-1 “Basis of design for glass”, January 1999,
prEN13474-2 “Design for uniformly distributed loads”, February 2000,
prEN13474-3 “Calculation of strength and load resistance of glass”, 2005, updated 2009.
While the general philosophy of the document was in accordance with the Eurocodes, the norm
never saw a final approval and was ultimately replaced by prEN16612.
5
prEN16612 “Glass in buildings” resumed the earlier regulatory attempt of prEN13474 in recent
years, with the final aim to issue a dedicated Eurocode for glass design. The document is divided
in two sections:
prEN16612-2013 “Determination of the strength of glass panes by calculation and testing”,
prEN16613-2013 “Laminated glass and laminated safety glass”.
To this day, this document has not yet received final approval.
3/237
Chapter 1 – Structural Glass Introduction
safety glass: the first, doesn’t require safety properties and can be used for acoustic
isolation or decorative purposes, while the latter is required safety with respect of
the adhesion of glass fragments to the interlayers and the ability to withstand the
impact of colliding bodies according to standardized procedures. Lacking a strong
central coordination, the regulatory framework on glass production and design has
followed different - sometimes divergent - paths among European countries. Some
have interpreted and reused part of the provisions of American6 and Australian7
standards, which arguably are some of the most organic attempts to a
comprehensive regulation body. A list of applicable country-based regulations and
technical documents follows:
AS 1288-2006, Australia (2006),
DIN 18008, Germany (parts 1/2/3/4 2010-2013, parts 5/6/7 in preparation),
CSTB 3574, France (2011),
CNR-DT210/2013, Italy (2013),
ÖNORM B 3716, Austria (2016),
ASTM E 1300, US (2016).
Alongside those and other documents, which are openly attempting a wide-
ranging approach to issues related to glass design and production, lives a crowded
universe of standards (ISO, EN, VDA, DIN, ETSI, etc.), regulating different aspects
related to glazing. An extensive description of these norms, along with their Italian
UNI EN counterparts, was made in a recent work by M. Blascone [9].
An outlook on the LG market can be obtained using Google Trends8 for the query
“laminated glass” (Figure 1.4). Results show a decreasing trend which has been
stable for almost a decade, with higher peaks around 2004-2006. Small and rich
emerging countries are leading the sector interest along with developed markets 9.
100
80
60
40
20
Figure 1.4 – worldwide interest in time for the query “laminated glass”
6
ASTM E 1300 “Standard Practice for Determining Load Resistance of Glass in Buildings” for
the US, currently used throughout the United States and parts of Canada.
7
AS1288-2006 “Glass in buildings” prepared by the Joint Standards Australia/Standards New
Zealand Committee BD-007, Glazing and Fixing of Glass, to supersede AS 1288—1994.
8
A web facility of Google Inc., based on Google Search, that shows how often a particular search-
term is entered relative to the total search-volume across various regions of the world.
9
On a scale from 0 to 100, where 100 indicates the location with the highest frequency of search in
proportion to the total number of searches in that location: Singapore (100), Australia (63), Hong
Kong (63), U.A.E. (57), New Zeland (52), Malaysia (45), U.K. (42), South Africa (37), Philippines
(24), U.S. (20), Canada (19), India (17), etc. [updated November 2017]
4/237
Chapter 1 – Structural Glass Introduction
1.1 MATERIALS
LG products are assemblies of several glass plies, two at the very least, bound
together with adhesive polymers. These plastic sheets are usually one order of
magnitude thinner than the glass plies, and for their inner position among the
composite material are usually referred to as “interlayers” 10.
Lamination was introduced in the automotive and aeronautical industries for the
assembly of forward windshields, to prevent the scattering of sharp shards upon
glass failure from inducing likely lethal injuries to pilots and drivers. As the first
safety glass with cellulose interlayer turned yellow in a few years of light exposure,
the bonding layer was replaced in 1933 by cellulose acetate, made from the
reaction of cotton with acetic acid [10]. By 1939 this was again replaced by the
cheaper polyvinyl butyral (PVB), which has then dominated the marked up to this
day. Almost one century later, despite major improvements in manufacturing
processes and increasingly more options for polymers to be used as interlayers, the
basic concept of glass lamination has not yet been revolutionized: New LG
products now conquering the construction industry are conceptually comparable
to the earliest elements in cars and other means of transport.
Polymer/glass laminates are employed when safety and security preclude the use
of monolithic glass [11]. Contrary to monolithic glass elements, these composite
materials gain some level of impact resistance and ability to dissipate energy [12–
16], along with post failure load bearing capacity. These and other properties, most
of which are yet to be adequately modelled and understood, are always
accompanied by the ability to hold glass fragments: a property connected to the
adhesive properties of the glass/interlayers interface which is the first unalienable
requisite of LG products.
Intuitively, the properties of LG composites derive from the properties of the
comprising materials11. In the following paragraphs, a general overview of the key
features and properties of the individual materials is presented. Bearing in mind
the scope of this work, which is aimed primarily towards the use of transparent
glass elements in construction industry from a civil engineering point of view, the
reader should not be surprised not to find adequate mention of some among the
known and important properties of the materials. For instance, optical or thermal
transmittance properties of glass and interlayers will not be mentioned at all, while
only a general overview of the fracture mechanics is discussed. Should the reader
require more insight on those topics, references are provided to exhaustive
scientific literature works.
Modelling the mechanical properties of LG from a structural engineering point of
view is no easy task. To showcase some of the most relevant phenomena, the
following scheme provides a broad overview of models and key properties that
need to be mastered to perform a reliable and cost-effective LG design, divided for
the three key fields of inquiry:
10
from latin intĕr : in the middle, in the midst. As in: “situated between layers”.
11
In most applications, two materials are used: one type of glass and one of interlayer.
Nevertheless, specific designs may require or suggest the adoption of several glass types, diverse
interlayers or the embedding of instruments and various devices [309–311].
5/237
Chapter 1 – Structural Glass Introduction
6/237
Chapter 1 – Structural Glass Introduction
17
Sir Alastair Pilkington and Kenneth Bickerstaff, pioneered the technique in St Helens,
Merseyside in 1952. In this “float glass production process” (also known as the “Pilkington
process”), molten glass rests on top a bed of molten metal, entering as a liquid around 1000°C and
exits as a solid, at about 600°C. The cooling must be carefully controlled up to 450°C and can be
sped up successively. This method gives the sheet uniform thickness and very flat surfaces. Due to
this technique, glass plies have a “tin” side (because tin is often used as molten metal, but other
metals can be chosen) and an “air” side, with slightly different properties.
18
Process of slowly cooling hot glass after they have been formed.
19
Often named ‘‘Griffith flaws’’, in the order of a few micrometers, produced during handling by
the manufacturer as well as the user in ordinary ambient conditions.
7/237
Chapter 1 – Structural Glass Introduction
While the first three approaches aim to reduce the various coefficients in fracture
mechanics models [34–36], the introduction of surface compression remains the
easiest, more reliable and cost-effective method for glass strengthening. Tempering
and toughening processes were developed to induce compressive actions in the
outer surface of glass, with an equilibrating tensile stresses in the inner part of the
material, which is not affected by the same vulnerability to superficial defects
[32,37]. This has beneficial effects in the overall tensile resistance, as it is mainly a
function of stress concentrations around the imperfections. Tempering can be
obtained either with a thermal process or a chemical one20: while thermal tempered
products can be “fully tempered” or “thermally toughened”(essentially, a lighter
form of thermal tempering), building codes only contemplate one type of
chemically tempered product (Table 1.2). Other properties of interest for glass
design are collected in Table 1.3.
Table 1.2 - characteristic tensile strength for glass (min. required [38])
float (annealed) MPa
thermally toughened MPa
thermally tempered MPa
chemically tempered MPa
The release of the eigenstresses, induced to enhance the performance of glass
elements with toughening or tempering, influences crack patterns and distribution
upon failure. Figure 1.6 shows how the relief of internal stress distributions may
result in widespread shattering patterns upon glass failure.
Figure 1.5 - stress distribution in chemically or thermal tempered (or toughened) glass
Cracks propagate in glass at the about speed of sound in the material21. From a
20
The French Francois Barthelemy Alfred Royer de la Bastie is credited with first developing a
method of thermal glass tempering in 1874 [313]. Nonetheless, inadvertent effects of glass
tempering had been already observed in the past (cfr. footnote 16 p.6). Thermal tempering consists
in re-heating the glass above its softening point (620-675°C) and following rapid cooling via
compressed air jets. Because of the fast cooling the outer surface contracts while the core remains
partially molten. When the core finally cools down, it must deal with a rigid outer surface that is
therefore pulled toward the center and becomes tightly compressed.
In the chemical strengthening process, alkali-containing glass are immersed in an electrically
heated bath of alkali salt . During the bath, an ion exchange occurs between the host alkali
ions of glass and invading ions from the salt. Because the invading ions are larger in size than the
original ones, the resultant stuffing in the near-rigid atomic network of glass leads to the
development of high surface compression and balancing by interior tension, whose distribution
depends on the diffusion depth and glass thickness [33].
21
The assumption in fracture mechanics is that the propagation of cracks occurs when the rate of
energy return from fracture is greater than the energy necessary to create new surfaces. Once this
happens, the incremental fracture energy grows while the incremental surface energy stays the
same, and the system is unstable. Crack will then propagate "as fast as possible", which for an ideal
continuous media is the speed of sound. Other works [35] used high-speed photography to reveals
8/237
Chapter 1 – Structural Glass Introduction
Figure 1.6 – typical failure patterns for a) annealed, b) heat-strenghtened and c) fully
tempered glass (Overend, 2007 [39])
Figure 1.7 - tempered glass fragments expelled upon bending-induced failure (Laboratorio
Strutture, Università degli Studi di Firenze, Italy 2017)
that cracks can spread through glass at speeds of hundreds of meters per second, or roughly half the
speed of sound in glass.
9/237
Chapter 1 – Structural Glass Introduction
22
Terminal velocity where is the mass of the object, is gravity,
air density, is the external surface and the drag coefficient.
23
While most LG products assembled for architectural purposes have the interlayers positioned
within the glass plies (i.e. most of the external surface is glass), laminates may comprise a single
sheet of glass included among synthetic resin sheets. In this latter type of laminated safety glass
products, which constitutes a niche product and will not be further referenced or studied in this
manuscript, the surface of the synthetic resin sheet which is not laminated to the glass may
optionally be coated with or laminated to another synthetic polymeric material [56].
24
In most production lines, pre-assembled components are heated in a chamber filled with oil or
air to about and pressure of to to soften and fluidize the interlayer
sheets and cause them to adhere intimately to the glass panes. Nevertheless, vacuum-enhanced
[314] and non-autoclave [315] techniques have also been developed.
25
“Ethylene-vinyl acetate (EVA) is still the market leader and standard encapsulation material for
solar cells […]. For a long time there have been no alternatives on the horizon, as finished
modules produced since the Nineties have a service life of at least 20 years and all series products
have to be certified. It is therefore difficult for other plastics to gain a foothold on the market if the
points in their favour are rarely measurable.” (Kuraray, 2009). Nevertheless, the research in PV
laminates has literally boomed in the last decade [316] and numerous researchers are dedicated
now more than ever to finding cheaper and more efficient products [101,317–319].
10/237
Chapter 1 – Structural Glass Introduction
26
some manufacturers, resellers and professionals are used to call those “single thickness”,
“double thickness”, “safety thickness” and “heavy-duty thickness” respectively. Intermediate
thicknesses that multiples of are not the product of simple multiple folds are also found
in the market, but less frequently.
11/237
Chapter 1 – Structural Glass Introduction
27
The most widely used polymeric films for glass lamination are: Polyvinyl Butyral (PVB)
[168,320], Ethylene Vinyl Acetate (EVA) [321] Polycarbonate (PC, often used in combination
with other interlayers for its limited adhesive properties) [84], SentryGlas® (SG) [54] and DG41
(enhanced PVB). Examples of other interlayers and brands, sometimes with dedicated decorative
or insulation properties, are StormGuard® (enhanced PVB interlayer alleges improved
performance for hurricane-resistant applications), Vanceva® Color (decorative), Saflex®
SilentGlass (acoustic), S-LEC™ (used mostly in the automotive sector as an UV shield with
acoustic/thermal insulating properties). A fairly recent review on the research status, developments
and applications of different polymers in LG was performed by M. Teotia and R. Soni [48].
28
Just to name a few worldwide: 1. DuPont Co., Ltd (US). 2. Kuraray Co., Ltd. (Japan) 3. Solutia
Inc. (US) 4. Huakai Plastic Co., Ltd. (China).
12/237
Chapter 1 – Structural Glass Introduction
Figure 1.9 - thermo-mechanical curves for a) elastic modulus [67] and b) strain rate [68]
Working conditions of interlayers used in LG should never involve temperatures
close or higher than the melt temperature, which corresponds to the states of the
lamination process, one that should never be reached again during the lifetime of
architectural LG elements. Operational temperatures for interlayers commonly
correspond to the glassy, leathery and rubbery phases (I,II&III in Figure 1.9). For
this reason, the glass transition temperature is usually an important and insightful
parameter for LG design and is provided by manufacturers (Table 1.4), stressing
the difference between very different responses of the material above and below.
Note that, while SG and DG41 are patents from specific manufacturers, PVB and
EVA are a classes of materials produced by many companies; as a consequence,
results for those materials should be looked at as averages or likely ranges.
Table 1.4 – glass transition and (melt) temperature for some interlayer materials
SG DG41 PVB EVA
55.0a (94.0)a 46.0b (n/d) 19.5c (113.0)c -43.0d/-33.1e (65.0f/80.0e)
a b c
2014, Kuraray (www.kuraray.com); 2014, Saflex (www.saflex.com) and [69]; 2016, Kuraray, confirmed by
d e
some researches [70–72], values may vary for other manufacturers; 1983, Cuddihy et. al. [73]; 2012, Agroui
f
et. al. [74]; 2006, Chen et. al. [75].
Traditional methods and models used in civil engineering are scarcely able
accurately predict the mechanical behaviour of most plastic materials; besides, a
complete mathematical formulation able to describe such response is currently
missing. For most design purposes, a phenomenological methodology appears to
be the most reliable approach, with different empirical models which can be used
to solve specific problems, under well-specified boundary conditions.
29
Below the glass transition temperature, a reduced molecular mobility doesn’t allow the polymer
to reach a state of thermodynamic equilibrium [291].
13/237
Chapter 1 – Structural Glass Introduction
Haward and Thackray in 1968 [76] (HT in the following) proposed that the
response of thermoplastic polymers could be divided is the sum of two terms,
in turn describing different physical mechanisms. While cross-links among the
molecules, along with and other weaker chemical actions, are accountable of the
response for short-range deformations , those tend to be less relevant for
higher strains, where the response is driven from the tethering and entanglement
of the long and branched polymeric chains30.
Figure 1.10 - classical stress decomposition (Haward and Thackray, 1968 [76])
30While this model was issued for true stress-strain diagrams, an analogous decomposition can be
made in an engineering-stress over engineering-strain diagram, neglecting the cross-sectional
reduction due to Poisson effect [164]. In fact, multiple transitions between the two will be used for
the analysis of tensile tests in §5.1. While the substance of this decomposition, its original idea, will
be used within the upcoming analyses, the graphical interpretation brought by Haward and
Thackray in Figure 1.10 is believed to have been originally flawed, having considered the presence
of a strain softening region within true stress-strain diagrams. In the following §5.1.1.2, a
discussion is brought regarding the fact that softening may only be visible in engineering diagrams.
31The key element that has been considered for the design of the viscoelastic investigation
campaign (cfr. §3.3.7 and §4.1).
14/237
Chapter 1 – Structural Glass Introduction
poly-(vinyl butyral) and about 30% plasticizers and softeners to enhance plasticity
and toughness. Condensation leaves some free hydroxyl groups, which are then
able to tie to similar groups originating from the Silicon atoms of glass [79]; these
bonds are at the origin the origin of the good adhesive properties of PVB.
PVB is a recyclable material. In fact, no incompatibilities were observed in
blending various grades of commercial and recycled PVB, suggesting that blends
and recycled polymers themselves could be used as substitute materials in glass
lamination [71].
PVB has been either described as an amorphous [80,81] or semi-crystalline [82,83]
thermoplastic. The mechanical behaviour of PVB has been proven to be highly
nonlinear, time- and temperature-dependent and capable of major extension. The
scientific literature includes several experimental results on PVB, which will be
referenced and compared throughout this manuscript to the newly presented
results.
1.1.2.2 SG
SentryGlas® ionoplast Architectural safety glass interlayer (in short, SG) is the
commercial brand of a Ionoplast Polymer (IP) produced by Kuraray Inc. 32 ,
primarily composed of ethylene/methacrylic acid copolymers, with small amounts
of metal salts. Compared with PVB, SG allegedly has higher stiffness and strength,
so much that in an early announcing of the material in 1997 the company stated:
“SentryGlas® [Plus] interlayers are the building block for a new generation of
structural laminated glass: a high-strength polymer used for glass laminators and
window companies to make advanced safety glass for automobiles and buildings”.
Pioneering the present historic phase, witnessing a boom in interlayer material
development, the genesis of this material was arguably due to the growing
awareness among LG designers and architects that PVB had very limited
mechanical properties, especially concerning long actions at higher temperatures.
SG belongs to the IPs family. Ionoplast polymers and copolymers have been
proven useful in safety laminates intended for structures requiring a high degree
of penetration resistance. Examples include hurricane-resistant glazing and
structural elements such as glass staircases and glass balustrades. In particularly
demanding applications, the use of IP interlayers such as SG in safety laminates
have shown ballistic resistance [84]
In literature, SG has been described as a semi-crystalline thermoplastic material
[45,81]. Ionic polymers like SG are synthesized by connecting ionic groups to a
carbon chain backbone. If the ionic content within the non-polar polymer is at most
10%, the ionic polymer is defined as a ionomer [85,86]. While the production of SG
is a patented and preserved secret, we known that IPs are synthesized by
neutralizing the carboxylic acid groups of acid copolymers. Incomplete details on
those chemical procedures are available in patents and scientific papers [40,87].
The final fabrication processes of SG incorporate an extrusion blow molding,
resulting in a monolayer product of the desired thickness.
Arguably few details on the mechanical properties of SG are known, apart from
those willingly disclosed by the manufacturer. Preliminary independent analyses
32
This product was patented as SGP® by DuPont Glass Laminating Solutions/Vinyl, a segment of
the DuPont Packaging & Industrial Polymers company acquired by Japanese chemical
manufacturer Kuraray in 2014.
15/237
Chapter 1 – Structural Glass Introduction
with tensile tests on material specimens have shown that this material has an
elasto-plastic behaviour at room temperature (about ), with a noteworthy
initial modulus ( to ) [54,88,89] (cfr. Appendix D). From a thermo-
mechanical point of view, this polymer appears to have major advantages over
PVB for a large range of temperatures ( to 60°C) [54]. Analyses within this
manuscript aim to share more light on the properties of this material: as a material
itself and most of all as component of LG elements.
1.1.2.3 DG41
Saflex® DG, also known as “DG41”, is the most recent interlayer material among
those studied in this manuscript. It is produced from plasticized PVB, so much that
some [90], along with the manufacturer itself, have been referring to this new
material as “structural PVB”; which is arguably a poor choice of words, given that
PVB has been rightfully considered a structural material itself for over half a
century. DG41 was announced by Solutia Inc. in 2012, allegedly suited for:
“applications in which increased interlayer rigidity and high glass adhesion are
required relative to standard glazing interlayers” 33.
Where the part “standard glazing interlayers” should be intended referred to the
overwhelmingly popular use of PVB. DG41 is supposed to combine:
“the benefits of a rigid interlayer with the features of glass containment, UV
screening, edge stability and clarity” 33.
Upon release of this new product, J. Schimmelpenningh, manager for Solutia said:
“[…] Saflex DG interlayer tends to keep glass intact at high and low
temperatures, even after impact and under load. Saflex DG is unique in that
it combines the meaningful performance and processing characteristics of
standard PVB interlayers while allowing for balanced engineering practices,
which can lead to more efficient whole building design.”33
Such bold claims have not yet been adequately backed up by experimental
evidence, hence the need for unbiased analyses and studies (cfr. Appendix D).
DG41 is comprised by long and heavy molecules ( ) [69] and is
already difficult to process in extrusion, compared to other more highly plasticized
polymers (plasticizer level is about 20 parts of per hundred parts resin [91]). This
low plasticizer level in increases the melt state viscosity, resulting in a difficult
extrusion process. For the current state-of-the-art production lines, attempting to
increase further the glass transition temperature of DG41 interlayer by dropping
the amount of plasticizer will result in an unacceptably low flowability for
extrusion, therefore necessitating a casting procedure which is not cost-effective.
Being a PVB-modified material, the good chemical adhesion of DG41 to the glass
surface essentially relies on the same chemical bonds described in §1.1.2.1.
33
Solutia® press conference, St. Louis (US-MO), 9 February 2012.
16/237
Chapter 1 – Structural Glass Introduction
copolymers include a wide range of applications; also, extensive efforts are being
made to replace PVB films with alternate films.
Some polymers are not included in this treatise have found thriving and wide
success in specific applications. Namely, because of its UV-screening and ageing
properties [73], EVA is leader in the PV sector. EVA is a polymer which fluctuates
from partial crystalline and thermoplastic to amorphous and rubber-like, but an
amplified amount of vinyl acetate increases strength and ultimate elongation.
A fairly recent and extensive review of most interlayers used for glass lamination,
most of whom will not be dealt with in this manuscript, can be found in the work
of M. Teotia, R. K. Soni [48].
Lastly, as this topic will be present in part of the following discussion, it is worth
mentioning the existence of so-called “multilayer interlayers” [69], interlayers
comprising several coextruded or pre-laminated polymeric layers. Most
applications of this technique aim to include within the whole interlayer thickness
polymers with prominent mechanical properties but poor adhesion with glass. In
doing so, the adhesion is dealt with traditional and reliable materials such as PVB,
while the bulk of the mechanical resistance is assigned to the core layers (e.g.
polycarbonate).
1.1.3 ADHESION BETWEEN GLASS AND INTERLAYERS
Adhesion is the macroscopic measure of the chemical attraction molecules of one
kind exert on molecules of a different kind on an atomic level. According to the
properties of the surfaces, adhesion can rely upon several mechanisms [94–96]:
mechanical interlocking, chemical bonding, physical adsorption and electrostatic
effects [96,97].
Within LG elements, adhesion between glass and polyvinyl-based interlayers is
primarily due to the chemical bonds that form between alcohol groups comprised
within the polymeric chains and the silanol groups on the surface of the glass
(Figure 1.11). These combine to form a reversible but dense network of weak
hydrogen bonds and beyond this point to some extent true chemical bonds [79,98].
Figure 1.11 - PVB-glass chemical bonding mechanisms (Keller and Mortelmans, 1999 [79])
34
85°C and 85% R.H.
17/237
Chapter 1 – Structural Glass Introduction
techniques have been employed to measure the strength of the adhesive bond, no
single test has proved to be superior, more general or reliable than the others. In
the last decade, effective non-destructive adhesion monitoring techniques (NDT)
have been developed [100].
CST tests have been used to investigate the strength of the adhesion bond for PVB
and EVA blends [99], showing that the layers of EVA exhibit a hyperelastic
deformation until catastrophic delamination of the interface takes place and that
hydrogen bonds among glass and PVB are easily degraded by moisture leading to
a large decrease in adhesion after DH.
While the research on EVA and PVB (and DG41, as its adhesion is owed to the
same molecular interactions than PVB) allowed to draw guidelines for degradation
conditions, details on the chemical adhesion mechanisms of ionoplast interlayers
have not been disclosed by the manufacturer nor yet seem to have been sufficiently
investigated. Comparative observations on post-breakage effects of weathering on
LG elements with PVB and SG showed that the strength of the adhesion bond for
ionoplast-laminates is weaker than for polyvinyl-based polymers. Moreover,
experimental analyses at high temperatures ( ) showed that the strength
of the adhesion bond for ionoplast-laminates is weaker than the internal cohesion.
As adhesion is an essential requirement for LG (it is the reason for which glass
lamination was invented in the first place), Even though the analyses within this
work were not directly focused in investigating the causes of delamination nor
modelling its final effects, a description of the observations was deemed to be both
useful and necessary (cfr. §4.2.4.5, §5.1.3 and Appendix C).
18/237
Chapter 2 – State-of-the-art and Aim of the Research
2 STATE-OF-THE-ART
AND AIM OF THE RESEARCH
The expansion of the LG market in the construction industry, together with
the development of an ever-growing number of interlayer materials, comes
with new problems and questions. Most of these open topics are related to the
contribution of interlayers to the performance on laminates. The scientific
community and the standardization bodies are working to find reliable
methods for safety assessment, characterizing materials properties and
devising cost-effective production technique. In this Chapter, the scope of the
manuscript is put into perspective after an overview on these efforts and a
review of the state-of-the-art.
1
For instance, the performance of interlayers to be used in glass lamination should rather not be
evaluated for compressive or flexural actions nor for direct impact, as those properties are not
primarily tied to the role of the material within the sandwich structure of LG. On the other hand,
tensile or shear tests could yield useful results for practical applications.
2
In the LG field, it appears that the standardization bodies are guiding behind the scenes the
scientific research the toward increasing levels of detail. Namely, standards are providing
researchers with a framework of basic assumptions, ways of thinking and methodologies that can
be accepted by the housing market. A similar scientific evolution path is described by T. Kuhn in
[322] in terms of shifting paradigms: “the paradigm determines the kinds of experiments scientists
perform, the types of questions they ask, and the problems they consider important. A shift in the
paradigm alters the fundamental concepts underlying research and inspires new standards of
evidence, research techniques and pathways for theory and experiment.” Such paradigm-shifting
revolutions produce “a shift in the standards by which the profession determined what should
count as an admissible problem or as a legitimate problem-solution.”
19/237
Chapter 2 – State-of-the-art and Aim of the Research
While the topic of structural elements made with glass can be wide and
multifaceted, according to the intended purpose of this work, the discussion will
exclusively address issues related to LG and specific aspects connected to
interlayer mechanics. Therefore, the reader should not be surprised if other topics,
which are also debatably needing improvement or refinement, will not be
mentioned or thoroughly analysed (e.g. general design principles, properties of
glass as a structural component, issues related to sealants and other adhesives,
insulated glazing, etc).
2.1 SHORTCOMINGS IN THE REGULATORY FRAMEWORK
A brief but wide-ranging review of regulations on architectural LG has been
discussed when introducing the topic in Chapter 1. Nonetheless, it appears that
providing a complete overview on the regulatory framework is a complicated task
for at least three reasons: first, numerous European standards or guidelines live
side-by-side with former regulation, some active and others outdated or
withdrawn (but still used in practice) [7]; second, while the glass industry can rely
upon internationally acknowledged manufacturing standards, regulating the LG
industry involves the necessity to standardize an emerging interlayer market,
populated by a multitude of new products along with traditional ones; third, while
several simplified models exist3, a fundamental approach for LG design, one that
should be based on strong theoretical bases, is still missing.
While some standards do not officially attribute to interlayers any structural
capacity within LG elements4, even those who do are not providing nearly enough
information for actual model implementation. Italian technical document CNR-
DT210 [38] is arguably among of the latest and most complete attempts for
standardization of LG design. Nonetheless, the document states:
“Mechanical properties of the elastomers [used as interlayers] are strongly
dependent from time and temperature; in any case the stiffness of these polymeric
materials is much lower compared to glass, so even provided the adhesion is
effective, polymers are subject to high angular sliding.” 5
And,
“[…] the viscoelastic behaviour of polymeric interlayers is modelled as linear
elastic materials. In particular, the shear modulus G, which generally varies
according to the operating temperature, loads duration and aging, is considered
constant.” 5
Statements like these, considering the outcomes of the experimental campaigns
presented within this manuscript, will likely be seen as oversimplifying the
response of a diverse and complicated panorama of interlayer materials.
3
Examples of standard-adopted models for LG response are: Épaisseur équivalente des vitrages
feuilletés (CSTB 3488, 2011), Interlayer stiffness family model (prEN16612, 2013), Wölfel-
Bennison (sandwich model [323], 1987, also adopted by ASTM E1300), Effective Laminate
Thickness (2009, [89]), Enhanced Effective Thickness (EET, 2012 [113], adopted and integrated by
CNR-DT210).
4
On the other hand, such standards may propose an increase in the maximum design stress for
glass plies in LG elements, which means they are indirectly accounting for the collaboration
between glass slabs, due to the coupling effect provided by the interlayers (e.g. TRLV ,2006,
“Technische Regeln für die Verwendung von linienförmig gelagerten Verglasungen”).
5
Translated from Italian.
20/237
Chapter 2 – State-of-the-art and Aim of the Research
Throughout the past two decades, with a remarkable acceleration in the last five
years, the scientific community has been progressively moving towards a more
rigorous approach in dealing with LG load bearing capacity. This new paradigm is
based on first principles for laminates mechanics and on proper, time-dependent,
constitutive equations for the interlayers [64,107,108]. While this approach seems
promising, it also appears at this time excessively complicated and lengthy for
most practical LG structural applications, inevitably requiring complex numerical
models (provided the materials properties have been established in a reliable way).
Standards and technical documents on occasion mention to these numerical
analyses as level-2 and level-3 models, with 2-D shells and three-dimensional solid
elements respectively. Comments and examples come with a few prescriptions, but
mostly relate on geometric or numerical issues. So far, for instance, requirements
and suggestions have not included specific details on the constitutive equation to
be used for defining the interlayer elements response6, which is understandably
one of the key issues at stake.
For the sake of presenting features and voids in the existing regulatory framework,
along with perhaps some misleading instructions, the aforesaid Italian standard is
taken as case-study. From the point of view of a structural engineer attempting a
design over LG, it appears that two choices may be source of decisive dilemmas:
choice of the model for LG response among the many available,
selection of the most appropriate constitutive equations for interlayers: linear
elastic (LE), linear viscoelastic7 (LVE), non-linear elastic (hyperelastic HE), non-
linear viscoelastic7,8 (NLVE).
An often-underestimated aspect to bear in mind, is the general lack of
understanding of polymers and adhesives peculiarities by structural engineers and
architects. Training in material science guided those professions toward arguably
simpler – or at least different - materials, that are only marginally influenced from
variables like time and temperature. With that in mind, it is of paramount
importance that standards and guidelines help them in choosing the most
appropriate design path. Besides, the awareness is growing that polymer
particularities might require adapted assessment methods of corresponding
building products and applications, specifically regarding time-dependant
behaviour and ageing phenomena [109]. For instance, scientific evidence has
shown that viscoelastic models are capable to describe the evolution of interlayer
properties over time [37,64,110–112], while hyperelastic models may be more
suited for the interpretation of brief effects [14–16]. Simple distinctions and
directions like this one seem to be missing in the existing standards framework.
6
More precisely, there are references to elastic approaches à la Newmark (Newmark et al. [160]),
but no mention is made relatively to constitutive equation other than LE.
7
For linear viscoelastic materials, the isochronous response is linear, but the effective modulus
decreases with time, so that the stress-strain curves at different times are separated from one
another. When a viscoelastic material behaves nonlinearly, the isochronous stress-strain curves
begin to deviate from linearity at a certain stress level [324].
8
Non-linear viscoelastic models for interlayer materials have already been adopted in LG studies,
so it made sense to include them in this list, as they appears to be the frontier of the research in this
field [108,325]. Nonetheless, such models are not mentioned in standard and guidelines (e.g. CNR-
DT210), as neither are mentioned some more complex models such as the visco-elasto-plastic ones
(linear and non linear, that to my understanding have not yet been successfully adopted).
21/237
Chapter 2 – State-of-the-art and Aim of the Research
In dealing with time and temperature effects for the polymeric interlayers
response, standardization bodies seem to have at least partially acknowledged the
inherent difficulty of designers to deal with those matters with a rigorous
approach. Therefore, as previously hinted in the quote from the Italian standard,
linear elastic simplified models (e.g. EET model by L. Galuppi & G. F. Royer [113])
are being provided, which can account for constant values of the interlayer shear
modulus . To deal with a diminishing modulus over time and as a function of
temperature (due to the viscous deformations of the polymers), a designer ought to
perform multiple analyses for varying values of the shear modulus.
Provided a reliable formulation to calculate the effective thickness of laminated
glass can be identified, a significant obstacle arises when attempting to find reliable
values of for given time and temperature.
9
At an international level, the European Committee for Standardization (CEN) has issued the draft
standard prEN-16613 “Glass in building - Laminated glass and laminated safety glass -
Determination of interlayer mechanical properties” (first issued in 2013, updated in April 2017),
proposing a classification of polymers in “interlayer stiffness families”. This ranking system has
been said by many to be taken with great caution, because while being very simple and easy-to-
use, it is deemed to be excessively inaccurate and flawed [326] (for instance, it neglects the
influence of known important factors such as the loading conditions, typology and size of the LG
element, cross-sectional details on the laminated assembly, etc).
10
The term “full-scale” is not explicitly adopted in [38], but it seem implied by the context.
22/237
Chapter 2 – State-of-the-art and Aim of the Research
23/237
Chapter 2 – State-of-the-art and Aim of the Research
Even though the fundamental origin of all properties of interlayer materials could
be - in principle - reunited with a unified theory, this has so far proven to be out of
the reach for the current scientific understanding. Far from such research of a
fundamental truth, investigations on interlayer properties for architectural
purposes have mostly shifted towards a phenomenological level of understanding,
splitting up different properties (i.e. interface and bulk responses, chemical and
optical properties, etc). To introduce to the aspects which are most relevant for this
work, a breakdown in two categories seems useful:
investigations on the interface strength and effects of weathering and ageing,
works oriented towards assessing bulk properties of interlayers.
While a traditional classification system exists, dividing plastics in thermoplastics,
thermosets and elastomers according to the type and density of the cross-links [10]
(cfr. §1.1.2), such division will not be used for this treatise of the scientific research
state-of-the-art for three reasons: first, it is not always obvious which category a
product belongs, as pointed out by other researchers [7,59], therefore the
classification might emerge as equivocal or inconsistent, second, most of the
scientific research has been made on PVB, with far less investigations on other
materials, and third, this classification is rarely mentioned by the scientific works
which will be mentioned in the following.
As a general idea, good practice would be to design and perform experimental
analyses which are a good representation of the final design, rather than trying to
extrapolate mechanical properties from different – although standardized – testing
techniques. Joined to a cautious indication of safety factors, a comprehensive
overview of meaningful properties for well-defined applications is the most
reliable approach when designing an experimental campaign on thermoplastics.
11
Mainly during the lamination process, but also throughout the whole life span of the LG
elements, due to infiltration of water and moisture along the outer borders.
24/237
Chapter 2 – State-of-the-art and Aim of the Research
with the chemical ones) and that a mechanical action is ultimately needed to
detach the materials, even with a severely jeopardized adhesion.
The great influence of storage humidity for interlayers, prior to the lamination with
glass, which appears to be vastly exceeding the effect of moisture and weathering
for laminated elements [65]. An almost dedicated research in that direction was
performed by M. Froli and L. Lani from University of Pisa in 2011 [98], concluding
that the most important parameter in the lamination process is the humidity
storage of the PVB.
Among the latest researches on investigating adhesion in LG in 2012, Y. Sha et al.
[125] studied the correlation between the peak load and the mechanism of
debonding for PVB, resulting in a relationship which could be fed to FE models.
Using peel tests, and a different approach, a study on glass-PVB interface
properties was performed a few years after by J. Pelfrene et al. [126], highlighting
that there seems to be no clear relationship between the peeling force and the
Mode I activation of fracture for characterization of the interface strength [34–36].
Moreover, this work underlined that no certain correlation could be established
among the presented experimental results and other publications in this field,
because the large-strain behaviour of PVB interlayer is yet to be satisfactorily
documented and understood. Adhesive resilience of specimens with initial
delamination at different regions and ratios was recently carried out by E. Dural in
2016 [127]. Using a high-resolution camera and a numerical morphological
analysis, Dural was able to obtain a relationship between delamination ratio and
the residual strength of delaminated glass units. The analysis concluded that the
behaviour of LG units is highly nonlinear, and that the strength of laminated glass
is influenced by the location and extent of the delaminated surface.
Figure 2.1 - specimens used in the tensile tests (Santarsiero, 2017 [128])
C. Louter and other researchers have also shown interest in the topic of adhesive
properties of interlayers in a series of recent works between 2014 and 2017. First,
effects of environmental impacts on durability of laminated glass plates with
interlayers were investigated for three types of interlayers: SG, EVA and PVB [129].
This initial research opened to two relevant “twin” publications for shear [130] and
tensile [128] strain rate correlations with temperature for laminated connections
between glass and metallic connectors. In these two researches, the authors
investigated both properties of a Transparent Structural Silicon Adhesive (TSSA)
and SG interlayer, observing two interesting phenomena with the loss of adhesion:
the whitening in TSSA and the development of bubble within the SG adhesive,
hinting to the hypothesis of a dominant hydrostatic component of the stress tensor
as the consequence of the large confinement state of the adhesive. The
experimental investigations also showed that temperature and strain rate
variations have important effects on the mechanical response of the connections;
namely, that resistance is proportional to the logarithm of the strain rate, and
decreases at high temperature. More specifically, a logarithmic law was proposed
25/237
Chapter 2 – State-of-the-art and Aim of the Research
for the strain rate effects for both TSSA and SG connections, while linear and
inverse hyperbolic-tangent-based laws were instead proposed for the TSSA and SG
temperature effects, respectively.
Results of the last experimental analyses suggest that much has yet to be done for
understanding and reliably anticipate the adhesive properties of interlayers used
for LG. Besides, while the research on PVB can boast many researches, those are
admittedly not yet fully relatable for the numerous aspects of interlayer rheology
which come into play. Conversely, research on new interlayer and adhesives is not
yes sufficiently developed to yield wide-ranging conclusions.
Led by L. Andreozzi, within 2013 and 2015 a team of researchers from Universities
of Pisa and Firenze also delved into the mechanics of interlayer adhesion and
correlations to environmental effects [131–133]. In their project, the authors
presented studies on thermo-viscoelastic properties of PVB, for varying weathering
actions (humidity, thermal cycles and UV radiation). The analysis highlighted
distinct damage mechanisms, promoted in different extents by the different
exposure to environmental effects. Results proved for the first time that RH can
influence both the adhesive and bulk properties PVB, calling for further
investigations to separately analyse these features. The authors concluded that
degradation related to solar radiation can produce major changes to mechanical
properties of PVB interlayers, to an extent that originally reliable numerical
evaluations on the response of LG structures may become incorrect over time.
bulk properties
Among the first works to show interest in bulk properties of interlayers, a research
issued in 1981 by J. R. Huntsberger included a study of mechanical performance of
PVB [134]. As usual, boosted by the needs of the automotive department rather
than the architectural one, a model for windshield fracture event was described in
terms of interactions among mechanical responses of interlayers and fracture
characteristics of glass. Anyhow, this pioneering research still seems very much
oriented towards taking steps for increasing safety with respect of the expulsion of
glass fragments in the event of colliding bodies.
The first work to openly mention PVB properties for architectural glazing is
possibly the paper from C Vallabhan et al. published in the Journal of Materials in
Civil Engineering in 1992 [135]. Shear force-deformation characteristics of PVB were
studied via direct-shear tests (Figure 2.2), similarly to those performed by
geotechnical engineers for determining soil characteristics in a specified plane.
While failing to emphasize temperature effects and openly reference to viscoelastic
properties, the major conclusion of the research was that the shear modulus of PVB
interlayer is “very dependent on the amount shear strain, and not on its thickness”.
Figure 2.2 – scheme of a setup used for direct shear tests (Vallabhan, 1992 [135])
26/237
Chapter 2 – State-of-the-art and Aim of the Research
Thermal properties and the dependency on those of the mechanic response of PVB
were investigated a few years later by Egyptians N. M. S. El-Din and M. W. Sabaa
[136]. While the research was certainty more oriented towards a chemical and
optical analysis of degradation rather than dependency of mechanical properties,
the conclusions discussed how the intrinsic viscosity decreases with increasing
temperatures. According to the authors, this indicates that the thermal degradation
is mostly accompanied by a main-chain scission, rather than splitting of the cross-
links (cfr. §1.1.2).
Apart from similar trailblazing works, of which other examples can also be found,
as indicated in the first part on this paragraph, the largest part of the research on
interlayers properties is found within the last two decades. A detailed comparative
study, complete with a state-of-the-art of researches and methods used for
assessing of the mechanical interlayers, was recently issued from the same Italian
team from Pisa and Firenze mentioned for their analyses on adhesive and
durability properties of interlayers (see above). As the panorama of techniques and
hypotheses is admittedly extremely variegated, the authors took a first effort for
classifying test methods used for the characterization of interlayers, according to
the time load history. The convenient scheme proposed within the work of G.
Ranocchiai et al. [118] is re-proposed in Figure 2.3 with some integrations.
Figure 2.3 - schematic representation of the interlayer studies proposed by Ranocchiai et al.
in 2017 [118], integrated with a few additional and relevant works.
References for studies cited in Figure 2.3: a[43], b[137], c[11], d[138], e[139], f[140],
g[141], h[142], i[143], j[54], k[144], l[83], m[145], n[146], o[147], p[148], q[65], r[149], s[150].
27/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.4 - intermolecular and network interaction scenarios (Vegt, 2002 [59])
12
Visser in 2009 seems to address the same question with regard to the distinction between
reversible and irreversible deformations. [7,270]
28/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.4 shows a few qualitative examples of plausible decompositions for some
hypothetical thermoplastic materials on both engineering and true stress-strain
correlations (indeed, for a given material and testing conditions, there only is one
correct decomposition).
As an empirical rule, studies have confirmed that intermolecular effects, are
dominant for small and very small strain, while network effects, when present and
significant, govern the response for large deformations [72,151–154]. Likewise, the
emerging viscous properties of thermoplastics are believed to be related to these
two microscopic mechanisms [57,155]:
viscoelasticity due to the uncoiling of the molecules network,
viscoplasticity ascribable to breaking of intermolecular cross-links.
An attempt for a combined description of intermolecular and network response
was proposed by Klompen in 2005 [156] and recently re-discussed by Delincé in
2014 [7] (in the following, KD approach/decomposition). This analysis incorporates
both the thermo-mechanical history (represented by the parameter [157]) and
softening kinetics, attempting to decompose the total response into different terms,
each dependent only on some of the known influential parameters. The proposed
constitutive model was then used to study polycarbonate response [50,156],
distinguishing the contribution of secondary interactions due to molecular
interactions from the network effect 13. In turn, a further decomposition was
proposed for the driving stress in a strain-rate dependent rejuvenated yield stress
and a yield-drop , function of the thermomechanical history (Figure 2.5).
This decomposition assumes that both the yield-drop14 and “rejuvenated” stress-
hardening functions are independent from the strain rate 15.
13
This method is consistent with the original 1968 Haward and Thackray (HT) approach [76]
14
The softening part of the curves in Figure 4.4, well recognizable in the experimental analyses for
SG and to some extent also for DG41.
15
Both hypotheses that doesn’t seem to be accurate in our analyses (cfr. Figure 5.2 123 with
comments in §5.1.1.2). Moreover, the yield-drop function would be required to output initial
negative values to correctly describe the material response for very small strains, something that
doesn’t seem to have an intuitive physical meaning, if any.
29/237
Chapter 2 – State-of-the-art and Aim of the Research
16
Using these results, the effect of different interlayers and temperatures have been studied in
recent years by El-Shami et al. [327], comparing them with FE models and new experimental data.
30/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.6 - scheme of four-point bending tests used by Louter [171] and Briccoli Bati [145]
17
At that time known with the commercial name of SGP, where the “P” stands for “Plus”.
18
By the means of which an approximate solution to the elastic problem of the composite glass
beam subjected to four-point bending test was proposed.
19
Even if the authors never mentioned this term, the effect observed called “tension stiffening” by
other researchers in succeeding works.
31/237
Chapter 2 – State-of-the-art and Aim of the Research
20
Possibly for its user-friendly characteristics and adaptability to numerous practical design
situations, this is the approach is proposed in the Italian technical document CNR-DT210 [38].
21
The “effective laminate thickness” idea draws from the intuition of Bennison et al. [89], in turn
elaborating on the work of E. Wolfel n 1987 [323], to introduce an equivalent moment of inertia
for the LG cross section. Once the effective stiffness is established, deflection and curvature can be
calculated and, consequently, the maximum stress in the single glass plies.
22
[14,98,211,220,328–331,104,106,108,139,141,145,158,171], only to reference a few.
32/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.7 - 3D FE model by C. Bedon, 2017 [184], for the numerical interpretation of
2014 tests by C. Louter et al. on pre-stressed LG beams [186] (cfr §2.2.3)
2.2.3 REINFORCED LG
In literature, with the term “reinforced laminated glass” scientists and engineers are
generally referring to composite elements based on a LG core, with the addition of
stronger and/or more ductile materials which are either able to enhance the SLS or
ULS performance. The knowledge about synergistic features of LG reinforced
composites relates varies depending on materials and techniques used. Numerous
examples of such performance-enhanced LG can be found across the body of
international researches, once again with a decisive increase in the last decade.
The interest in reinforcements for LG structures start to come up by the end of the
last century. A concept for a steel-glass-polycarbonate beam presented in 1997 by
F. A. Veer et al. in the bi-annual Glass Performance Days conference was later on
updated in 2001 and once again in 2003 [187]: a three-ply LG with polycarbonate
interlayers and a stainless steel profile in the tension zone. In later years, more
examples and similar ideas have been proposed and tested; namely: plain or
perforated metallic plates and flanges (2003-2008) [188–190], steel-reinforced
concrete (2004) [191], plain or hollow steel bars (2005-2012) [45,192], stainless steel-
wire meshes, thin perforated metal sheets and fabrics of high strength fibres (2008)
[193], carbon fibre (2009-2016) [5,194], glass fibres [195], pultruded Glass Fibre
Reinforced Polymers (GFRP) [196,197] or timber (2010-2016) [198–200]. Some of
those concepts are presented in Figure 2.8.
23
Abaqus Unified FEA - SIMULIA™ by Dassault Systèmes®
33/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.8 - a) & b) embedded hollow/plain steel c) steel or GFRP profiles d) carbon-fibre
e) embedded perforated metal plate f) steel flanges g) timber h) reinforced concrete 24
For most of those studies, the whole purpose of the reinforcements was to enhance
the post-failure performance of LG elements, by avoiding sudden collapses.
Therefore, such reinforcements had been charged with minor – if any - structural
responsibilities for serviceability conditions.
Figure 2.9 – tenting the “Zappi1” beam specimen25 in 4-points bending (Veer, 1997 [201])
24
All of the examples proposed are shown in a vertically symmetric configuration. Nonetheless, it
is not uncommon in LG beams to find reinforcements only on the lower part of the cross-section,
which is affected by tensile actions in most SLS and whose post-failure behaviour can hardly
account any structural contribution from glass. On the other hand, it is known and it will be shown
again within this manuscript, that fractured glass is still able to withstand substantial compressive
actions (cfr. §4.3.4 and §5.2.4). For each of the presented solutions, the number of glass plies is
indicative and may in fact vary for different studies and applications.
25
A special method of laminating glass plates in an overlapping pattern, with the advantages that
failure is concentrated in the individual segments and specimens of any size can be made by
manufacturing them out of smaller components. While this techniques didn’t use steel
reinforcements, the same concept of segmented LG beams was used with stainless steel
reinforcements in later researches from the same team (Delft, 2003) [187], concluding that
reinforced LG is a safer structural material.
26
The origin of this group is ascribable to Prof. D. Dicke (1924 - 2003) in 1989, who asked to
describe “the ideal transparent structural material”. He proposed to design a structure with an
imaginary material having a motley set of characteristics, derived from those of well-known
structural materials. The quest for Zappi, this still unknown transparent wonder material, was
34/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.10 – five-points bending tests for SG-laminated beams reinforced with a) solid-
and b) hollow-steel profiles (martens, 2017 [205]; cfr. Figure 2.8)
35/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.11 - cross section of the GFRP-reinforced FL tested (Speranzini, 2015 [197])
Figure 2.12 - four-point bending results on P-FG (left) and P-TFG specimens (right)
For both solutions, the structural behaviour is similar; however, the mode of failure
is different: for P-FG the load increases after the first failure and diminishes
progressively with increasing number of cracks, until the final collapse, for of P-
TFG specimens, the initial failure is reached because of fragmentation of the whole
array of tempered glass plies. Overall, this solution proved that some kind of
ductile and hardening structural response can be obtained for reinforced LG
beams, with load-displacement curves which are similar to the desired load-
carrying behaviour described by K. Martens et al. [206] for pre-stressed LG (see
below).
2.2.4 PRE-STRESSED LG
The term “pre-stressed laminated glass” can be used to refer to LG elements whose
glass is subject by design to a – significant - tensile state prior to the application of
external loads. Two macro-categories of elements can be identified:
elements with axial self-balancing loads, designed to induce a compressive state
in the glass plies (“pre-stressed” or “post-tensioned” LG28),
out-of-plane “cold-bent” LG structures, whose internal stress distribution is
proportional to the extent of the curvature, thickness of the sandwich structure
and coupling abilities of the interlayer.
28
The term post-tensioned has been extensively used in scientific literature referring to such LG
elements. Anyhow, in my opinion this choice ay be slightly confusing and the term pre-stressed
should be always preferred. As a matter of fact, “post-tensioning” refers to the act of applying a
tensile force to some tensile-resistant component of a composite LG material after (i.e. post-) the
lamination process is concluded. Arguably, because the goal of the whole procedure is not to tear
those accessory elements, but rather to put the glass - or part of the glass - in a compressive state
before (i.e. pre-) the application of the design live loads, the term pre-stressed (or pre-compressed)
is more suited and less prone to misunderstandings.
36/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.13 - a) cross-section and lateral view of LG reinforced and pre-stressed using FRP
pultruded bars (Cagnacci, 2009 [211])
Limits of adhesive reinforcements bonding have been pointed out in 2014 with the
research of Louter et al., showing that pre-stress can be achieved by means of
mechanically anchored or adhesively bonded tendons. Due to the actions of the
tendons, initial fracture loads of pre-stressed LG elements can be increased up to
compared to reference beams. Also, with the bridging of the glass cracks, a
post-fracture load-carrying mechanism is created. Therefore, the residual load-
carrying capacity is enhanced up to – of the initial fracture load [186,207].
Figure 2.14 - method for bonding pre-tensioned tendons to the lower edge
of laminated glass beams (Louter, 2014 [186])
37/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.16 - desired load-carrying behaviour for LG beams (Martens, 2015 [206])
On the other hand, regarding cold-bent LG elements, the act of forcing a single-
curvature shape after the lamination process is concluded gives also birth to pre-
stressed LG elements: structures which need to cope with an internal state of stress
prior to the application of external live loads. The behaviour of such elements is
not easy to model, as the stress distribution in the glass plies varies with time after
the initial bending phase because of the viscoelastic properties of the interlayers,
which causes the decay of the shear-coupling effect [212].
The possibility to cold-bend LG has been used for less than twenty years, with the
earliest scientific reports dating back to the start of the century [213]. The
distinction between short-term and long-term response of this technique has been
clear since the early days: while at short term the glass panes “work together”, and
the sandwich structure will act approximately as a single layered pane of the
thickness of the whole sandwich, at long term the panes of the sandwich act mostly
independently from one another as separate panes (Figure 2.17).
Figure 2.17 – qualitative stress distribution in cold-bent LG at short (left) and long term
behaviour (right) respectively (Vakar and Gaar, 2004 [213])
38/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.18 – Roof over the IJsei bus station in Amsterdam in a) an early 2004 concept,
which was completed in 2011 and b) a picture in 2017.
29
Where curved glass plies are first obtained by heating them up to the material softening point
(cfr. §1.1.1) and curved individually into the desired shape using moulds. Hot-bent glass plies are
then successively laminated to form LG elements. While this process allows to shape virtually any
kind of curved glass elements, either single-curvature and double-curvature ones, it is very
expensive and doesn’t leave a perfectly reflective glass surface.
30Cold-bending a double curved shape is possible, but typically results in very high stresses in
glass, so that the economic advantage of cold-bending in manufacturing is lost. However, recent
advances in theoretical algorithms allow for the discretization of any surface using only single
curvature panels, thus allowing the construction of smooth double curvature glazing of any form
[332–334].
31
This technique can already count some practical application, like the Strasbourg Railway Station
(cfr. Figure 1.2.c) and the aforementioned bus station in Amsterdam (Figure 2.18).
39/237
Chapter 2 – State-of-the-art and Aim of the Research
concentrations can kick-start the uprising of cracks, thus influencing the maximum
load which can be carried by structural LG. As peak stress concentrations are
independent from the type of element and can affect any type of glass product (cfr.
§1.1.1), localized loads induced failures cannot be safely included into the design
process. Besides, what happens to the glass after it has failed is an important
design concern; and both designers and customers must accept that LG can fail, at
least partially, due to concentrated impacts and actions. Ultimately, a successful
and widespread strategy to avoid catastrophic effects is to include sacrificial plies
on the outer layers on LG [42,192,215]. Therefore, the understanding and
modelling of damaged LG elements must be seen as fundamental parts of the
design process. While studies on the behaviour of unreinforced LG beams are
available [62,211,216,217] (also cfr. §2.2.2), those rarely delve in performance
assessments after the failure of the glass plies [11,158,159,166].
In the last decades, researchers have picked up the challenge of understanding the
post-failure response and properties of LG. Once again, the research started from
the automotive department. One of the first research on post failure LG
performance can be traced back to the year 2000, in a work by Japanese J. Oda et al.
[218], studying variously damaged automotive windshields residual impact
resistance. Mechanisms of crack initiation and propagation were investigated by
Chinese team J. Xu, Y. Sun et al. between 2011 and 2017 with a couple of
noteworthy papers [219,220]. In the former, with preliminarily impact fracture
experiments, authors investigated crack propagation in PVB-LG; in the latter, an
explanation to the basic crack propagation mechanism on radial multi-crack was
proposed, providing guidance for retroactive impact information extraction and
improving cracking-resistant structural design and material selection.
Models have been proposed for in-plane and out-of-plane bending. The latter can
account for decisively more in-depth analyses, having yield in the last years both
an homogenized analytical model by L. Galuppi and G. Royer [104] and a Non-
linear discrete modelling proposal [215] by D. Baraldi, A. Cecchi and P. Foraboschi.
Both researches focused on tempered glass failure. Arguably, while influencing
dramatically the post-failure LG configuration, tempered glass failure produced a
more homogeneous material in the post-failure phase, compared to annealed glass
failure. Major research conclusions were that fully-damaged tempered LG
elements can still give rise to a resisting system, with a degree load-bearing
capacity and stiffness strongly dependent on the interlayer material. Results of the
analysis also showed that the interlayers are “stiffened” by the contact with the
adherent glass shards, but a degree of delamination is supposed to occur at the
glass-polymer interface. From a numerical point of view, it was shown that non-
linear Discrete Element Modelling (DEM) is a practical solution for analyzing the
behaviour of fully damaged LG elements made with tempered glass plies.
Compared to out-of-plane analyses, far fewer works seem to have been performed
on the post-failure behaviour of in-plane loaded LG. A recent experimental paper,
issued by Chinese X. Huang et al. [221], described an experimental campaign
which was carried out on a significant number of full-scale PVB- and SG-laminated
beams in short term three-point bending at room temperature. While PVB-
laminated were found to have very poor post-failure performance, the research
was able to produce a simple analytical model to predict the post-failure strength
of the SGP LG beams. Still, this work seems to be the only example regarding
studies in this area of research.
40/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.19 - a) homogenized model problem by Galuppi-Royer [104] with hint to interface
delamination and b) DEM example from Baraldi-Cecchi-Foraboschi [215]
41/237
Chapter 2 – State-of-the-art and Aim of the Research
Figure 2.20 - a) interlayer specimens for tensile tests and b) LG specimens tested for
assessing the viscoelastic properties of interlayers ( )
42/237
Chapter 3 – Mechanics
3 MECHANICS
Despite its intrinsic brittleness, glass itself can be included in the category of
“traditional” solid building materials. Most of increased complexity when dealing
with LG is due to the peculiar behaviour of thermoplastic interlayers. This
complexity is well introduced by A. Ali:
“Elastomers present a very complicated mechanical behaviour that exceed
linear elastic theory and contain large deformations, plastic and viscoelastic
properties and stress softening”. A. Ali et al., 2010 [222]
In the following, an overview of the fundamental mechanics needed to understand
the response of LG is provided. While some chapters are openly more dedicated
glass and others are dealing with interlayer mechanics, it is once again useful to
recall the words of K. Gilbert (cfr. p.11): observations and models are generally
dependent on the scale of human problems, so that mechanics can vary once the
scale is changed to larger perspectives.
This chapter deals with the following topics:
small strains linear elasticity,
large deformations non-linear elastic response,
time- and temperature-dependency of solids,
non-reversible processes and plastic deformations,
Despite being essential to the understanding of the response of brittle materials
such as glass, the topic of fracture mechanics is not dealt with in this manuscript,
as countless dedicated works on this topic can be found in literature [32,34,36].
3.1 ELASTICITY
The theory of the elastic solid is briefly introduced, deriving the relations and laws by
which an elastic material transforms because of externally applied actions.
Elastic materials share the ability to deform under external actions and the
property to restore the initial shape when freed from said actions, independently
from the conditions upon which those were applied. The rule according to which
loads are applied is called loading process [223].
An elastic material or body is described as a reversible system, where energy
spend throughout the deformation is internally stored and entirely released upon
removal of the external actions. The ability of such materials to naturally restore
their initial state allows to acknowledge the existence of an entity such as the elastic
energy or strain energy, capable to produce (at its expenses) negative work from a
43/237
Chapter 3 – Mechanics
Different materials may store elastic energy in different ways. From a macroscopic
point of view, a distinction of interest for this work may be performed according to
materials showing a linear elastic response, a non-linear elastic response and a
time-dependent elastic response [224].
The treatise in the following paragraphs will first briefly mention linear elastic
solids, which are a good model for glass in the pre-failure phase, while non-linear
elasticity and time-dependent elasticity will be functional for understanding the
response of thermoplastic interlayers. On occasion, when presenting the various
mathematical formulations, specific mechanics will be immediately referred to a
specific component of LG, to clarify and synthesize the discussion.
3.1.1 LINEAR ELASTICITY
Among elastic materials (see above), a Linear Elastic (LE) one is defined by linear
constitutive equations, meaning that the relations between stresses and strains
are linear and homogeneous [223]. For LE isotropic materials, the correlation
between the components of tensors and allows for the formulation of the so-
called generalized Hooke’s laws1 :
Where is the Young modulus2 and the the Poisson’s coefficient3. These relations
are valid for the small strain domain and will be used frequently throughout this
work when dealing with glass; as their use is common practice in engineering,
these will be usually implied without direct reference.
3.1.2 NON-LINEAR ELASTICITY (HYPERELASTICITY)
The demand of predictive models for mechanics of rubber-like materials has been
pushing the research since the beginning of the 20th century, overcome the limits of
linear elastic models that have guided many of the technical developments
throughout the industrial revolution (namely, the restriction to infinitesimal
strains). This progress has led to the development and critical assessment of
constitutive equations for the nonlinear elasticity of rubber materials known as
hyperelastic models [225]. Such mathematical formulations can be used to model
materials exhibiting a non-linear but reversible response (Figure 3.1). Primary
applications of these models have been dedicated to rubbery polymeric materials
1
Robert Hooke, Wight 1635 – London 1703
2
Thomas Young, Milverton 1773 - London 1829
3
Simeon Denis Poisson, Pithiviers 1781 - Paris 1840
44/237
Chapter 3 – Mechanics
and foams exposed to large, reversible shape changes (e.g. tires, seals, conveyor
belts, base isolations of buildings, etc) [226,227].
Hyperelastic models are suited for situations in which the non-linear response of a
material is mostly independent from strain rate and, more generally, from the
strain history. Alike elastic materials, the hyperelastic response only depends on
the strain (i.e. any work done deforming an hyperelastic material is stored as
potential energy in a thermo-dynamically reversible process).
for isotropic materials in finite strains, a function is usually defined using the
invariants , and of the Cauchy-Green tensor, and sometimes through the
elongation components in the three spatial dimensions:
Dealing with rubber elasticity in this work, namely, with the response of interlayer
polymers used in glass lamination, the large strain domain of isotropic materials
must be considered. To choose the most appropriate SED functions to work with,
we first need to assess how the compressibility of such materials come into play,
both from a mathematical point of view and from the perspective of practical
4
Some models postulate the existence of an arbitrary SED function, then derivate the formulation
the constitutive material equations; some other bypass the definition of a SED function and directly
provide empirical constitutive mathematical formulations instead.
45/237
Chapter 3 – Mechanics
Figure 3.2 - hyperelastic models used within this manuscript, highlighted among some
popular ones (references can be found in the work by Steinmann et al., 2012 [225])
Polymers used as interlayer materials in LG transition among various regimes of
mechanical behaviour depending on the temperature: glassy, leathery and rubbery
(cfr. §1.1.2). Densely cross-linked polymers are the most likely to show the “ideal”
rubbery behaviour [227], that hyperelastic expression are intended to approximate.
Within this manuscript, hyperelastic models will be used extensively in Chapter 5,
to elaborate on tests results of Chapter 4. Mathematical formulations of the used
phenomenological models (Figure 3.2) are presented in the following paragraphs.
Starting from the SED formulations proposed by various authors, a walkthrough
on the analytical developments needed to calculate the stress-strain relations in
uniaxial tensile tests is carried out in the following.
Neo Hookean model
Among phenomenological hyperelastic models, the Neo-Hookean (NH) is
significant for two reasons: it is the simplest and does not involve arbitrary
constants. The SED for is defined as follows:
Where is the first invariant of the isochoric part of the right Cauchy-Green tensor
from the deformation gradient , is correlated to the variation
of volume for compressible materials, and the constants and
are depending on the shear modulus and the bulk modulus respectively.
46/237
Chapter 3 – Mechanics
In an uniaxial stress test with the force applied in the direction, the transversal
elongations are equal . With this relation and the previous equations, one
finds as a function of the principal elongation component and :
Substituting the expression of in the formula of and then differentiating for the
three elongation components (proof in [234]) yields the expression of the
principal stress components for a compressible material in an uniaxial stress test:
If no action is applied on the sides, radial stress components and are null;
thus, it is possible to subtract the last two expressions and get rid of the term in :
The relation between and can then be found solving the nonlinear equation in
with numerical methods. To highlight the variation of for changing and ,
the Newton–Raphson method was used for results in Figure 3.3.
compressible
incompressible
with
Figure 3.3 shows the effects of varying values of the bulk and shear moduli on the
function and a qualitative representation of the NH models described in Table
3.1 for given values of moduli and .
47/237
Chapter 3 – Mechanics
Mooney-Rivlin model
The Mooney‐Rivlin (MR) hyperelastic model proposed in 1951, is a special case of
the generalized MR polynomial model [235]. Based on empirical descriptions, this
is today one of the most regularly used models [230]; the SED function is:
Where is the 2nd invariant of the isochoric part of the left Cauchy-Green tensor:
Closed form expressions for a compressible problem are scarcely used, involving
the solution of complex differential equations with arbitrary parameters. For an
incompressible material, using the same relations introduced for the NH model,
the principal stresses for a MR incompressible material are shown in Table 3.2.
Table 3.2 - uniaxial elongation stress for incompressible Mooney-Rivlin material
true stress engineering stress
with and
Note that for the MR model degenerates in the NH incompressible one, and
that for any other value of the MR model results in softer material behaviour for
finite strain, while still having the same initial tangent modulus.
48/237
Chapter 3 – Mechanics
Yeoh model
Like the MR one, the Yeoh model also provides a phenomenological description of
the material response. The model is based on R. Rivlin's observation that the elastic
properties of rubber may be described using a strain energy density function
which is a power series in the strain invariants [230]:
The constant is related to the bulk modulus, so the right sum is non-null only for
compressible materials. Alike to the MR model, despite compressibility being
included in the SED formulation, a closed form solution is unknown using the
complete expression of . Hence, the Yeoh model is normally applied for
incompressible materials only ( ) choosing :
with and
49/237
Chapter 3 – Mechanics
with
linear relaxation
5
Cited in over two thousand scientific works at the time this manuscript is written.
6
Ludwig Eduard Boltzmann, Vienna 1844 - Triest 1906.
50/237
Chapter 3 – Mechanics
given value, and the stress required to hold the length constant is recorded as a
function of time [240]. As a result, most of the data available for viscous properties
of materials is available as (creep) compliance vs. time. In a creep test, strain is
proportional to the applied constant stress , and depends from time according to
the properties of the material through a creep compliance7 function [241]:
The response to a creep or relaxation test for is the elastic response (cfr. §3.1.1
and §3.1.2), which depends from the instantaneous moduli and
. Alas, although both creep and relaxation ultimately issue from the
manifestation of the same physical effect 9 , there is no known closed form for
connecting the mathematical formulations of the two functions [114,237,240],
hence:
Relaxation moves toward its equilibrium quicker than creep , and the
shape of the function can be substantially different [114,241]. Therefore, if one has
available creep data (e.g. the coefficients for the discrete formulations) it is not
easy to get insight on the relaxation model (and vice-versa).
If strains are small, the partitioning hypothesis is applicable; being the elastic
strain and the viscoplastic strain [242]:
The analytically exact creep compliance function can be found by solving the
convolution integral [241]:
creep compliance
7
To the letter "the act of complying"; as creep itself is the measured effect, the compliance is
intended as the way the effect is measured and mathematically contextualized.
8
As one can deduct from the equations, the creep function is expressed in , while the
relaxation modulus is in .
9
For polymers, the reversible sliding process of macromolecules over time.
51/237
Chapter 3 – Mechanics
a) b)
spring, following Hooke's law dashpot, following Trouton's law
a) b)
Figure 3.8 - rheological representation of a) the Maxwell element and b) the generalized
Maxwell model, of which the Wiechert model is a special case if
“instantaneous” Wiechert
creep compliance
10
Sir Isaac Newton, Woolsthorpe 1642 - Kensington 1726.
11
Emil Johann Wiechert, Tilsit 1861 - Göttingen 1928
52/237
Chapter 3 – Mechanics
the instant actions are applied. In the sum, are the creep compliance coefficients,
and the retardation times associated with those coefficients 12 . For , an
infinite compliance modulus can be defined:
If all terms of the sum are known, this allows to find from and vice-versa:
Therefore, the model allows for a different and equivalent "infinite" formulation:
“infinite” Wiechert
creep compliance
12
Note that a coefficient is not defined, because the instantaneous compliance -as the name
suggests- occurs instantaneously.
13
Example: given , and , the
infinite creep compliance is . Any
given term of the sum has accumulated only a fraction its final viscous effect when .
Because of the mathematical formulation of the Wiechert compliance, this fraction is equal to
. From a practical point of view, the viscous effect of each term can be
considered complete when , which yields a completion ratio . Note that the time
also corresponds to the change in the curvature from positive to negative for each creep coefficient
contribution in the time domain. To better explain this property, a comparative figure is
shown with the compliance function (left) and the curves relative to each element of the sum
(right); points where each term has achieved 95% of the its final viscous contribution to the
compliance are highlighted.
53/237
Chapter 3 – Mechanics
relaxation modulus
“instantaneous” Wiechert
relaxation modulus
“infinite” Wiechert
relaxation modulus (Prony)
The correlation between the different coefficients of the modulus can be found by
calculating in an analogous way as in the creep tests, yielding:
The previous equation for the “infinite” Wiechert relaxation modulus, which
allows for a simple representation of the relaxation problem 14, is presented in a
formulation which is equivalent to a Prony series15 [243]:
“A linear viscoelastic model with a given Prony series is identical
to a multi-network Maxwell model”. J. S. Bergstrom, 2015 [237]
14
Example: given , and ,
we find the infinite creep compliance is
, using the same considerations discussed in footnote #13 above, the shape of the
compliance function is:
15
Gaspard Clair François Marie Riche de Prony, Chamelet 1755 - Asnières-sur-Seine 1839.
54/237
Chapter 3 – Mechanics
16
Based on the Marquardt-Levenberg Method [335].
55/237
Chapter 3 – Mechanics
Defining the storage modulus and loss modulus , the equation can be simplified:
(relaxation)
Storage modulus
(relaxation)
Loss modulus
The storage modulus, defined as the stress in phase with the strain in a sinusoidal
shear deformation divided by the strain, is a measure of the energy stored and
recovered per cycle. Since both and are measurements of the stored
elastic energy, a dynamic measurement at frequency is qualitatively equivalent
to a transient one at , so the graphs for and are
approximately mirror images with respect to a vertical axis (Figure 3.10). In
particular, when is changing very slowly, [114].
Nevertheless, it has been pointed out that relationships among LVE functions such
as storage and relaxation moduli have bases in the theory of linear differential and
integral equations [249]. While several approximate interconversion methods can
be found [248,250,251], Mun et al. [252] proposed an exact interconversion
procedure for LVE functions of asphalt concrete from frequency-domain to time-
domain (although details on calculations are not given).
56/237
Chapter 3 – Mechanics
Figure 3.9 – LG specimen and testing in a torsional rheometer (Andreozzi, 2014 [148])
Figure 3.10 - relaxation & storage moduli plotted against and axes (Ferry, 1980 [114])
17
Observations were made in the first decades of the 20th century by K. W. Wagner and applied to
dielectric measurements.
57/237
Chapter 3 – Mechanics
science between temperature and shift factor has been proposed in 1955 by M.
L. Williams, R. F. Landel and J. D. Ferry [254]:
“instantaneous” Wiechert
creep compliance
“instantaneous” Wiechert
relaxation modulus
3.2 PLASTICITY
While elastic materials will always internally store a form of energy, able to restore
their initial configuration when external actions are removed (cfr. §3.1). In
opposition, another class of materials can be defined for those which undergo non-
reversible changes under certain conditions. From a thermo-mechanical point of
view, this inability to restore the original shape is ascribed to the activation of
dissipative mechanisms within the material structure. The power dissipated
during plastic flow is called power of dissipation , and is essentially positive [255]:
In engineering, the study of plasticity has often been related to the study of metals,
for which it has been observed that plastic deformations consist of slip in crystal
dislocations [256]. Hence, it is generally assumed that the onset of yielding and
power dissipation occurs when a fixed value of the maximum stress is reached:
18
Most of the times, the least-squares differences method is used.
19
While in first applications of the WLF equation, averages for many polymers of coefficients
and had been used, it later became clear that the variation among different materials were too
large to allow for the use of such “universal” values. In a slightly better approximation, fixed
values and can be used, in conjunction with an adjustable reference
temperature , which usually fell about above the glass transition temperature .
58/237
Chapter 3 – Mechanics
Several criteria can be used to define yield conditions, Tresca20 and von-Mises21 ones
being the most commonly accepted for metals. A perfectly plastic solid can thus be
described as a solid which undergoes unlimited plastic strains under constant
reference stress, when subjected to a homogeneous stress state [255].
Nevertheless, the theory of perfect plasticity seems to be not particularly relevant
and helpful when dealing of LG structures, as neither glass or interlayer materials
fill the profile for perfectly plastic solids (cfr. §1.1.1 and §1.1.2).
3.2.1 VISCOPLASTICITY
The fundamental time-independence assumption of all theories of plasticity does
not allow for a simultaneous description of both plastic and rheologic properties of
materials. In fact, there is no yield function for viscoplastic solids, but a stress- or
strain-rate-dependent inelastic behaviour is observed instead, this means that
when permanent deformations are produced, those are dependent on the rate at
which external actions are applied [239,257].
Constitutive equations for viscoplasticity can be based on springs, dashpots, sliders
and other elements, to model different responses such as perfectly viscoplastic,
elastic perfectly viscoplastic and elasto-viscoplastic hardening solids (Figure 3.11).
A creep test yields linear plot, while a relaxation test has a vertical asymptote for
and an horizontal one at zero for .
Elastic perfectly plastic solids (EPP)
The rate of plastic strain in EPP materials is a function of the yield stress of the
slider element. A constant stress flows through the slider element when the elastic
limit is exceeded; otherwise, the plastic strain is nil. The Bingham-Maxwell22 model
is shown in Figure 3.11.b [255,259]:
20
Henri Édouard Tresca, Dunkirk 1814 - Paris 1885
21
Richard von Mises, Lemberg (now Lviv) 1883 - Boston, 1953; the yield surface is also known
under the name of Maxwell, Huber, Hencky and von-Mises yield condition.
22
Eugene Cook Bingham, Cornwall 1878 - Easton 1945, and James Clerk Maxwell, Edinburgh
1831 - Cambridge 1879.
59/237
Chapter 3 – Mechanics
This classical formulation is not the only one proposed for HEVP materials,
besides, the function f is not specified a priori, opening to an array of possible
mathematical interpretations. As this research field is still looking for definitive
models for finite strain response of polymers, a few studies have concurrently tried
to tackle this topic, either with a reinterpretation of classical models or introducing
new approaches23.
3.3 ADHESIVE BONDS
While cohesion refers to the strength of the bonds within homogeneous materials,
adhesion is a measure of the tendency of different particles (or surfaces) to “stick”
one another. In engineering practice, adhesive joints can be designed to endure a
certain set of loads or displacements. In LG industry, adhesion is a relevant factor
for assessing structural performance in three situations:
securing glass to interlayers,
securing LG to an external frame or supports (seals),
bounding of reinforcements.
In fact, the purpose and most of the times the behaviour of LG ultimately
originates from strength of the adhesive bond between glass and interlayer
materials (cfr. §1.1.3). While adhesive properties of different resins24 to a variety of
reinforcements have been investigated (e.g. steel profiles and tendons, glass or
carbon fibres, concrete, wood, GFRP, etc. cfr. §2.2.3 and §2.2.4 for more), these
properties are of limited interest for this research. Regarding external adhesive
joints for LG elements, traditional seals are made using silicone-based compounds
(silicone structural glazing) [38]. This field is also far from the aim of this work,
where the focus will be set to studying interactions among glass and interlayers, in
which loss of adhesion is relevant for two orders of reasons:
23
The affine network hypothesis [336–338], based on the assumption that the number of entangled
points in the molecular chains remains constant during deformation [339]; multi-scale models for
semi-crystalline polymers in finite- [340,341] and infinitesimal-strain regimes [342,343], using a
spherulitic morphology with a radial assembly of crystalline lamellae and amorphous layers [339];
generalized yield surfaces [339], based on the Drucker and Prager criterion, including invariants of
the stress tensor and deviatoric part of the stress tensor; hyperelastic strain-hardening springs in
HEVP one-dimensional arrays similar to the one in Figure 3.11.c [258], intended to reproduce the
finite strain elasto-viscoplastic non-linear response of polymeric materials.
24
Most of the times epoxy resigns [209].
60/237
Chapter 3 – Mechanics
Figure 3.12 - geometry of common adhesive joints for LG (Adams and Wake, 1984 [38,94])
A single-lap joint model shown in Figure 3.12.a can be used to understand the
mechanics of multi-layered LG elements connections with any number of glass
plies (e.g. three- or four-layered LG elements loaded transversally with flexural
actions can be imagined as series of single-lap joints among adjacent plies).
While several arrays may be considered for adhesion of LG elements to the
substructure (e.g. step, butt, butt strap and double butt strap [94]), for the purpose
of this work and understanding the experimental analyses, it is only of interest to
focus on some specific aspects related to double-lap adhesive joints of the kind
described in Figure 3.12b. As the figure suggests, loads in single- and double-lap
joints are parallel but not perfectly collinear with the interlayers. Therefore, the
bending moment generated perpendicularly to the plane must be resolved within
the joint with stresses normal to the direction of the applied loads. Tearing and
compressive stresses are higher at the end of the joint and smaller in the middle.
The relative influence of these normal stresses is a function of the geometry of the
joints: for a fixed thickness of the glass plies, the longer the joint, the smaller the
normal stresses compared to shear ones. For long joints compared to glass plies
thickness, normal stresses are very small and can be neglected [65].
Parametrical study at varying temperatures for the shear stress distribution in SG
used as an adhesive for metals on LG elements has been performed by M.
Santarsiero et al. [130] (Figure 3.13). The analysis shows that the shear stresses
25
For instance, the Italian technical document CNR-DT210 [38], one among those attempting a
wide-ranging approach to LG design (cfr. §1 and §2.1), doesn’t provide numerical specifications
regarding the strength of the adhesion between glass and interlayers. Instead, the document states
(translated from italian): “It is necessary to point out that adhesion between glass and interlayer
involves a chemical bond; this circumstance requires […] to perform mechanical tests on stratified
specimens for the evaluation of the overall behavior. Since there are currently no regulatory
indications in this regard, […] tests on plastic materials only are allowed”.
26
For high temperatures, LG elements are heated to , a temperature held constant for two
hours; for humidity tests, specimens are put for two weeks in saturation conditions (
at a temperature of .
61/237
Chapter 3 – Mechanics
Figure 3.13 - a) shear stress distribution in SG for varying temperatures b) shear stress
factor at different temperatures (Santarsiero, 2016 [130])
Knowing the stress distribution, to evaluate the risk of delamination and joint
safety margins, a tri-linear model like the one shown Figure 3.14, can be used,
where the horizontal segment corresponds to the loss of adhesion among surfaces.
This mechanical problem is talked using a zero-thickness interface concept, for
which two adherents come into contact through an interaction surface. The local
interface response in mode II [34] is adopted: a LE branch up to the adhesion
strength followed by softening to the final displacement corresponding to
the perfect fracture formation. The response for a monotonic displacement is [262]:
62/237
Chapter 4 – Experimental Tests
4 EXPERIMENTAL TESTS
LG composites bear an inherent complexity, for their non-homogeneous and
anisotropic structure, but also for the intrinsic properties of the thermoplastic
materials used as interlayers. Direct experimental tests are often used to
address the lack of reliable theoretical models. Within this research, extensive
and different experimental investigations have been carried out, regarding
some of the less known mechanisms which are significant in LG design.
1
Obtaining viscoelastic or hyperelastic parameters requires a number of various tensile testing
procedures and a very large number of specimens [118,137,152,344,345].
2
Details on the uniformity of stress and strain are discussed in §4.2.1.
3
The term “conventional” is used by numerous researchers, as such tests are often prescribed when
dealing about plastics materials and are detailed by ASTM and ISO standards [66,264].
63/237
Chapter 4 – Experimental Tests
Outcomes of tensile tests are presented in the following paragraphs, while a critical
analysis an discussion is provided in §5.1.1.
4.1.1 DESIGN AND AIM OF TENSILE TESTS
Short-term properties of interlayer materials have been investigated following the
standards guidelines applicable to uniaxial tensile tests of plastics [66,264].
Throughout tests, the end grips of flat, dog-bone4 shaped, material specimens are
pulled apart at constant velocity, inducing an increasing strain which is mostly
concentrated within the narrow section of the body (Figure 4.2). The average strain
is measured between two points which are at a distance on the un-deformed
specimens.
Tested specimens were extracted out of larger, thick, interlayer sheets: half
the thickness of layers used for laminating specimens for the campaign dedicated
to study the viscoelastic properties (cfr. §4.2). According to the standards [66],
“Type IV” specimens are suited when direct comparisons are required between
materials with different rigidity cases (that is, nonrigid and semirigid) [66], which
appears to be the case in our study (cfr. §1.1.2.1, §1.1.2.2 and §1.1.2.3). For
, specimen dimensions shown in Figure 4.2 are fixed 5 : ,
, , , . A conveniently shaped manual die-
cutter was used to provide identical specimens from the larger interlayer pieces.
Once extracted, specimens were checked for scratches and imperfections. Flawed
specimens were discarded.
4
A technical term often used in literature for such specimens, due to the visual analogy with a
stereotypical bone (also sometimes called “dumbbell” shape).
5
Radii of fillets, along with several other geometric specifications, are also provided in the
aforementioned ASTM standard [66], they are not reported here for brevity.
64/237
Chapter 4 – Experimental Tests
6
Instron Illinois Tool Works Inc. Canton, Massachusetts
65/237
Chapter 4 – Experimental Tests
66/237
Chapter 4 – Experimental Tests
The failure of specimens always originated in a point within the narrow area,
never in the grip parts nor at the interface with the machine clamps.
Tests duration of was in the order of twenty minutes for , two minutes for and
fifteen seconds for .
The engineering stress in the narrow area has been computed, knowing the
applied force and original geometry of the specimen: with a cross-sectional
nominal area (cfr. §4.1.1), . The mean strain in the narrow
area can also be calculated knowing the initial distance and the relative
displacement of the two measure points: .
Stress results of performed tests are presented in Figure 4.4. Results show that the
influence of the strain rate is more relevant for DG41 and PVB interlayers, while
SG has a more consistent response at different strain rates.
Maximum tearing strength is usually of little interest for architectural LG
applications, as it is reached for very high strains which are not compatible with
serviceability limit states imposed by design. Nonetheless, for completeness in the
discussion of experimental results, results in terms of ultimate strain and
engineering stress are given in Table 4.2.
Since the existence of an elastic limit for plastics has been debated [7,78,114], the
term “elastic modulus” is indicative, only to measure a quantitative value of the
rigidity. Standards suggest two supposedly equivalent approaches to evaluate the
elastic tensile modulus of polymers. The first method, also called chord slope method
[264], is illustrated in Figure 4.5.a: two points A and B are found on the stress-
strain curves, corresponding to fixed values of strain. The slope of the straight line
interpolating the points defines the tensile “chord” modulus :
67/237
Chapter 4 – Experimental Tests
A close-up view on the average experimental results for small strains is shown in
Figure 4.6. Dashed lines represent the slope of the tensile modulus calculated with
the chord slope method (the regression method yields very similar results).
Figure 4.6 – mean engineering stress-strain curves and resulting chord-slope moduli
according to ISO 527 [264] for a) SG b) DG41 and c) PVB
68/237
Chapter 4 – Experimental Tests
While curves obtained for SG, the hardest among the three materials, are arguably
fit to be studied with both methods proposed by the standards, the same can
hardly be said for the other materials (Figure 4.6.b and Figure 4.6.c). For very small
strains ( ) both PVB and DG41 initially exhibit a very stiff response,
qualitatively comparable to SG and even higher in the case of DG41 for high strain
rates (i.e. and ). Nonetheless, this very high initial stiffness is immediately
followed by a quasi-constant load over strain segment, containing the whole strain
interval , proposed by the standards to compute the tensile
elastic modulus. The right part of Figure 4.6 shows that the tensile chord modulus7
is an acceptable approximation of the initial tangent modulus for SG, but is neither
very telling or consistent for both DG41 and PVB. Mean values of moduli
calculated with both standard-proposed approaches are in Table 4.3; results show
that the difference among the two methods is minor and, for this case study, the
regression modulus is always slightly lower than the chord modulus. The
convenience of using standard methods to define the tensile properties will be
further discussed in Chapter 2.
Table 4.3 – mean chord- and (regression)-slope tensile modulus in MPa
claw displacement
SG DG41 PVB
velocity
5 mm/min 431 (430)* 5.209 (5.209) 0.725 (0.704)*
50 mm/min 435 (435) 16.45 (16.38) 1.671 (1.660)
500 mm/min 1060 (996) 16.50 (16.45) 1.489 (1.412)
* statistically non-fully reliable according to [66,264] due to limited number of specimens
7
Tensile regression moduli curves are not shown, though they yield similar conclusions.
8
The Poisson’s ratio is defined within the infinitesimal strain theory [164], it has been questioned
whether the concept and definition of may be extended to large deformations [346]. For an
incompressible isotropic material, the relative error between the infinitesimal and finite strain
domain can be estimated by .
69/237
Chapter 4 – Experimental Tests
the same straightforward approach used for SG would yield a conservative lower
bound value of the same diagram.
70/237
Chapter 4 – Experimental Tests
The time-lapse shooting frequency was limited by the data acquisition system to a
maximum of one shot every two seconds (Table 4.4). To allow for a practical use of
the time-lapse images, reference measurement lines have been manually drawn10
on specimens as shown in Figure 4.9.
Table 4.4 – SG tensile tests photographic time-lapse shooting details
shooting freq. duration shots
5 mm/min 1 shot every 5s 29’28” 355
50 mm/min 1 shot every 2s 2’35” 79
500 mm/min 1 shot every 2s 19” 11
9
Optical device used is a Tamron SP AF 90mm F/2.8 Di MACRO 1:1 on a Nikon D90 camera
with CMOS 12,3-megapixel sensor. Software for time-lapse control is open-source
digiCamControl © 2015 Duka Istvan.
10
A thin permanent marker was used, the kind of those used to write on CDs.
71/237
Chapter 4 – Experimental Tests
Figure 4.10 shows the complete shooting sequence for the higher test speed (500
mm/s), including a final after-failure picture.
Figure 4.10 – SG tensile test at 500 mm/min time lapse, all shots11
11
The time-lapse video can be seen at https://www.youtube.com/watch?v=O17Q_z9DZ7U
12
Prior to this campaign, a similar procedure was implemented with good results by the author and
other colleagues in a previously published paper [65] and other works still to be published.
72/237
Chapter 4 – Experimental Tests
have been achieved, has been highlighted by numerous researchers [7,263] and in a
previous campaign carried in 2013-2014 [65]. In the conclusions of his 2014
manuscript, Didier Delincé stated (cfr. [7] §VI.3):
“The determination of the contribution of the interlayer should be assessed by
considering it as a component rather than as a material, and test
configurations should be designed in accordance. In particular, separate
characterization of adhesive properties and bulk material properties of the
interlayer seems not really useful with purpose of assessing the end
performances of products.”
Rather than focusing on the interlayer alone, the design of tests, starting from the
setup itself and then defining the procedure, aims to investigate a scaled version of
the actual LG beam. The mindset is to carry out simple and well-controlled tests
which are a good representation of in-situ working conditions for interlayer
materials. LG specimens were assembled using three -thick plies,
wide, with two -thick interlayers (Figure 2.20.b p.42).
Drawing from the basic concept behind shear tests for bolted double-lap joints for
steel elements and/or other composites [269], tests for LG specimens were designed
to apply a vertical force on the central ply. For the geometric properties of the
cross-section and the large difference in stiffness between polymers and glass, this
procedure induces almost exclusively shear actions through the interlayers, with
minor tensile and compressive stresses at both ends of the specimen [118,226]. The
setup in Figure 4.11 is similar to the one used in 2014 (cfr. Biolzi et al. [65]), with
the addition of the base contrast element, which will prove essential for the
improved design of the procedure.
The three-layered LG specimens being loaded vertically through the central glass
ply, actions normal to the interlayer surface were deemed to be minor, and a
procedure to record or quantify these actions was not put in place. Besides, local
shear forces transferred through the interlayer can be confused with the mean
shear stress under some conditions13; this assumption was used to gather stress
information from the load recordings.
13
cfr. §2.2.1, the outcomes of the work by C. Vallabhan et al. [135] and also E. Viola [223]. This
topic is further discussed in §4.2.4.4.
73/237
Chapter 4 – Experimental Tests
Figure 4.12 - testing devices with specimens inside the climatic chamber before the tests
14
Testing devices were built by Officine Meccaniche Rinaldi Srl in Milan and delivered to the
laboratory in the Politecnico di Milano, were tests were performed.
15
The term “quasi” is used in the sense that strain is not properly consistent, as a proper relaxation
test would require, but varies slightly over time.
74/237
Chapter 4 – Experimental Tests
structure through the base contrast elements shown in Figure 4.1116, highly reliable
data could be gathered on the viscous recovery in a second phase, contrary to
previous testing methods [65]. The applied load in phase II is null, as the
dynamometric rings are physically disconnected from LG specimens, which are
left free to return over time to their lower-energy state (Figure 4.13.a). Using the
superposition principle [114,237,240,241], it is possible to look at the load removal
as to the application of an inverse, “virtual” load (upwards), which is equal to the
load recorded at the end of phase I. Conversely, the “virtual” strain
, essentially a re-scaled and flipped version of , can be
looked at as the creep compliance effect of a constant, “virtual” load step
, as shown in phase II’ in Figure 4.13.b.
Figure 4.13 - design two-phases test output for stress and strain a) expected data to be
recorded and b) conversion of a recovery phase II to a “virtual” creep phase II’
Such testing technique aims first and foremost to gather data on LVE parameters
for creep and relaxation of three interlayers: PVB, SG and DG41 (cfr. §1.1.2). In
addition, some permanent plastic strain is expected, especially at high
temperatures ( in Figure 4.13.a, whose scale may only be considered indicative),
if a significant pattern can be detected, further analyses can be made in that sense
using known models for PVP or HEVP (cfr. §3.2.1), always bearing in mind that
such tests are not specifically designed to investigate viscoplastic parameters.
The launch of tests is achieved by adjusting the loading bolt (Figure 4.11) to a fixed
design displacement , correlated to the design strain . More details are
discussed with the presentation of preliminary results in the next paragraph.
16
More details on the mounting and securing technique are shown in Figure 4.14 and Figure 4.15,
and further discussed in the following.
75/237
Chapter 4 – Experimental Tests
17
As properties of interlayer materials vary significantly with changing climatic conditions, for
higher temperatures the limit was usually defined by the maximum displacement, while at lower
temperatures the maximum load was reached for significantly lower values of the displacements
( mm). The number of steps to reach the target load/displacement was achieved
accordingly (e.g. first load step imposed is imposed up to 5kN or 0.25 mm, whichever comes first,
then a second load step up to 10kN or 0.50mm, whichever comes first, etc.). All load steps were
performed by alternating the fastening of the three testing devices, for the whole loading process to
be performed in about three to five minutes.
76/237
Chapter 4 – Experimental Tests
While a precise lower limit for the conditioning time was defined, to allow setup
and specimens to reach thermal equilibrium, the duration of phases I and II was
not defined a priori. Tests were carried out for the longest amount of time, while
simultaneously managing several factors: opening hours of the laboratory, rate of
load decay at a given temperature, upcoming tests, etc.
As tests were carried out from the highest temperatures to the lowest ones (Figure
4.24 p.85) the duration of phase I and phase II shown in Figure 4.13.a was adjusted
on-the-fly for each test. The aim of this way of proceeding was to provide data that
could be later joined in a general mastercurve (cfr. §3.1.3.5) with the lowest amount
and extent of “gaps”. This apparently confusing and arbitrary criterion will be
clarified with the data analysis performed in §5.1.2.1.
4.2.2.1 CALIBRATIONS
A verification of the calibration for all LVDT transducers was performed, both the
displacement ones18 and those connected to the dynamometric rings19 (Figure 4.11).
Effects of temperature on the measurements of displacement transducers were also
investigated in the full range of test temperatures by laying a still object at half the
stroke of the instrument ( ) and cyclically varying the temperature in the
chamber. Minor variations in the instrument readings were observed in a range of
, but most importantly for tests, at each temperature the conditioning
time interval for consistent readings ( ) was observed to be lower than
, allowing for reliable readings throughout tests (cfr. §4.2.4.3 and §5.1.2.6).
4.2.2.2 PREPARATIONS OF SPECIMENS
Due to the way LG glass is manufactured [271,272], the edge of the interlayer was
not always spot-on the edge of the glass plies (Figure 4.14). Protrusions of
interlayer outside the perimeter of the specimens was generally in the order of
, although not always present. Cutting out of the excess interlayer was
only needed in some areas, to allow for convenient contact among all surfaces:
between the glass plies and the base contrast elements, between the upper loading
edge and the central glass ply and between the specimens and the lateral fastening
angular shims. Excess of interlayer outside those areas was usually deemed to be
irrelevant to the outcome of tests and therefore not always completely removed.
This cleaning procedure was mostly needed for all the SG-laminated specimens, as
they have been individually laminated by superposition of three 100x100mm glass
plies and two interlayer foils of the same size. On the other hand, PVB- and DG41-
laminated specimens were produced by partitioning larger LG elements 20 .
Therefore, those specimens already presented adequate finishing and didn’t
require extensive manual finishing.
18
Calibration curves for displacement transducers were given by the HBM certification body,
supplier of Politecnico di Milano, for a 10mm range and 0.001mm accuracy, within an operating
temperature range of 10÷60 .
19
Dynamometric rings calibration was done by means of a 50kN Wheatstone bridge load cell, in
turn calibrated by the same certification body of the LVDT transducers. The resulting calibration
curve for the dynamometric ring itself showed good linearity, with maximum deviation of ±0.5%
in a 0÷25kN range, and an accuracy of 0.05kN.
20Partitioning was achieved by slicing of large LG elements through a diamond band saw,
to produce smaller pieces the size of the design specimens.
77/237
Chapter 4 – Experimental Tests
Figure 4.15 – four phases (clockwise) for the installation of angular shims
78/237
Chapter 4 – Experimental Tests
79/237
Chapter 4 – Experimental Tests
Results in Figure 4.17 are referred to an unloaded load measuring device. Being the
ambient temperature generally close to , one can observe that higher test
temperatures have a greater influence in creating a fictitious load measurement. A
thermodynamic model was calibrated to understand why the recordings have such
an odd shape, details are in Appendix A.
For some tests, multiple load measuring devices were not used to perform actual
tests on laminated glass specimens, but were nonetheless inserted in the climatic
chamber and subject to the same temperature variations. Compared results
showed good consistency among all non-loaded devices, implying that results
from one device can be successfully used to correct another, which was the
assumption adopted in the following data analysis.
Figure 4.18 – typical position and evolution of cracks in laminated glass specimens
80/237
Chapter 4 – Experimental Tests
Typical crack patterns and starting points of fractures are presented in Figure 4.18.
Most cracks were very minor: either affecting just the outer surface of the glass to a
depth of a few millimetres (Figure 4.18.a), either affecting a limited volume of the
whole ply (Figure 4.18.b). In both of those cases, the interlayer interface was not
severely affected by the failure of the glass plies, thus results regarding the
viscoelastic properties of the interlayer could have been accepted anyway. On the
other hand, a few cracks were much larger, affecting a vast portion of the plies and
the interface, possibly influencing the tests (Figure 4.18.c and Figure 4.18.d).
The starting point of cracks was always observed close to the edges, suggesting
that this phenomenon may be connected to defects in the alignment of plies or
minor tolerances in the finishing of surfaces: stresses which were meant to be
transferred through a large contact surface were in fact redirected and
concentrated into significantly smaller regions in proximity of the edges.
A noteworthy remark is that the finishing of the edges did not seem to be
correlated to the formations of cracks. The uprising of cracks was observed only for
some specimens and for loads higher than , translating in nominal
compressive stresses around for the central ply and for
the lateral ones. This was an odd result, as those stresses are even lower than the
tensile limit of glass [273–276].
To prevent the glass from cracking, a thin lead foil was inserted between the
bottom of the lateral plies and the base contrast elements (Figure 4.20). The
downside of this choice is that the recorded displacement at the base of the central
glass ply now must account for some plastic deformations of the lead interface. To
allow for correct readings, the plastic deformations in the lead foils were addressed
before the beginning of tests (i.e. before the conditioning phase). To this aim,
several very high load ramps of kN were imposed with the testing devices on
the lateral plies alone, with the temporary introduction of lead foils also on top of
them.
Figure 4.19 shows that increments in lead plastic deformations gradually reduce
with an increasing number of cycles. When the relative displacement between the
beginning and the end of a high-load cycle was smaller than m, a verification
cycle was performed up to kN, the maximum design force. Lead plastic strains
were deemed to be settled with a residual deformation lower than m.
81/237
Chapter 4 – Experimental Tests
Figure 4.20 – a) setup for adjustment of deformations in lead foils and b) test setup
Figure 4.21 – setup for calibration of rings and base contrast elements
21 The first remarks were made comparing preliminary PVB tests to available data [38,65].
82/237
Chapter 4 – Experimental Tests
Figure 4.22 – a) load and b) displacement curves for base elements calibration
Curves shown in Figure 4.23 are the result of this calibration process; details on the
construction of these mathematical models are discussed in Appendix B.
22Each time a test was performed, the load cell was removed and reinstalled, to ensure that
results were not influenced by some peculiar arrangement of the elements.
23To grant full consistency with the following tests, the dynamometric ring readings were
used to calibrate model. Redundant load recordings must be looked at as an internal
feedback.
83/237
Chapter 4 – Experimental Tests
24
In the last part, the applied load is null, as the dynamometric ring is physically disconnected
from any surrounding object and the temperature is constant.
25
is a progressive two-digits number.
84/237
Chapter 4 – Experimental Tests
The list of tests and specimens is presented in Table 4.5, while Figure 4.24 shows a
Gantt diagram of the campaign progression26.
26
The list is “missing” some specimens (e.g. tests for SG start with the S09 and S12), in part due to
those used in preliminary tests and in part to the rejection of specimens for assembly flaws (e.g.
glass plies not perfectly aligned).
85/237
Chapter 4 – Experimental Tests
On top of both Figure 4.25 and Figure 4.26, conditioning, loading and recovery
phases divide the time axis (cfr. §4.2.1). The duration of the phases matches what is
shown in Table 4.5. Significant points have been highlighted:
loading point end of the conditioning phase, prior to the chamber
opening and specimens loading by fastening of the
tightening bolts,
thermal stabilization after closing the climatic chamber, when the
after loading dynamometric rings are in thermal equilibrium (Figure
4.27),
unloading point end of the loading phase, before the opening of the
chamber for loosening testing devices and unloading the
specimens,
27
HBM QuantumX CX22B-W (www.hbm.com)
86/237
Chapter 4 – Experimental Tests
thermal stabilization after closing the climatic chamber, when thermal effects
after unloading on dynamometric rings are settled for the most part (cfr.
Figure 4.27),
end of the test end of data recording, before the final opening of the
chamber.
Thermal effects were introduced in §4.2.2.5 and are discussed in Appendix A. To
define the moment when such effects start to be negligible, points were identified
on the measured load curves: the thermal stabilization was assumed to be reached
in the loading and recovery phases when (Figure 4.27).
Figure 4.27 - load measured by the reference dynamometric ring upon chamber opening
Figure 4.25 shows that even starting from the so defined stabilization point, some
thermal effects must still be considered (line for device 041 after thermal
stabilization at loading/unloading); this has been taken care of in the following by
subtracting the load recorded by the reference dynamometric ring from the
readings of the instruments used for testing the LG specimens.
4.2.4.1 LOAD DRIFT CORRECTION AND INITIAL DISPLACEMENT
The problem was introduced in §4.2.2.8. Figure 4.28 shows what the load drift
looks like, and how it affects the loaded dynamometric rings and the reference one.
Dotted lines are superimposed to the recorded measurements to highlight the
linear trend of the drift, calculated by linear regression of the load curves in the
viscoelastic recovery phase28. Figure 4.29 shows the result of the correction, where
the linear drift effects are subtracted from the original load readings.
28The regression base for the linear drift curve is defined by all the points comprised in the
interval within the thermal stabilisation point after unloading and the end of the test point
shown in Figure 4.25. The choice was made for all dynamometric rings to be effectively
disconnected by all surrounding effects in that time interval and in thermal equilibrium.
87/237
Chapter 4 – Experimental Tests
The comparison Figure 4.28 and Figure 4.29 shows that the slope of viscous load
decay in the loading phase is more consistent after the drift correction.
Figure 4.30 shows the normalized recordings for the displacement curves, obtained
by subtracting the displacement at the loading point (cfr. Figure 4.26).
The imposed displacement during the loading phase is not perfectly constant, as it
should be instead by theoretical design of a stress relaxation test (cfr. §4.3.2).
Instead, the imposed displacement slowly rises after the initial imposed value, due
to the elastic deformations of the dynamometric ring. Nonetheless, it is important
to note that this example is willingly referred to SG (cfr. §1.1.2.2) to show the worst
case scenario: SG being the stiffest among tested interlayers, results for other
materials are more concordant to the hypotheses.
88/237
Chapter 4 – Experimental Tests
29
Not all specimens had identical geometry, some were slightly wider or shorter than design.
89/237
Chapter 4 – Experimental Tests
load and displacement (Figure 4.31 and Figure 4.32), assuming a linear
variation of the displacements30 [223].
Figure 4.33.a shows stress results for the case study: for the viscoelastic recovery
phase, the virtual – constant – load has been set equal to the last
recorded load in the relaxation phase just before unloading: the dash-dot line,
representing the virtual load, intersects the end of the corresponding curve in the
loading phase. The time axis for the creep curve has been shifted, the new origin is
the unloading point at time (cfr. Figure 4.25 p.86).
Figure 4.33 – a) stress and b) strain in loading and viscoelastic recovery phases
Shear strain curves shown in Figure 4.33.b for the unloading phase have been
shifted, normalized and turned upside down according to the procedure (cfr.
§4.2.1). The new time origin is once again the unloading point at time (black dot
in Figure 4.31 p.89).
Shear strain and stress data over time allow to compute the values two
experimental moduli and [114,241]:
30
While values of shear strain are way larger than those found for more conventional building
materials, a definition of shear strain like that shown in the equations can be successfully used to
represent the interaction of the interlayer in the LG element in models, aiming to compute shear
strain and the corresponding shear stress transferred among the glass plies [135,167].
90/237
Chapter 4 – Experimental Tests
Even so, looking at the slope of the strain curves in the loading phase in Figure
4.33.b and remembering that this case study is referred to the stiffest among tested
materials, it seems not completely unacceptable to confuse the experimentally
calculated modulus with the relaxation modulus. In the proceeding of the
dissertation, this approximation will be accepted, and the consequences will be
carried on knowingly, up to the final conclusions. Therefore, by concisely referring
to the relaxation modulus , in all the following discussions, we will in fact be
referring to the experimental approximation of the relaxation modulus , rather than to
its true definition.
The experimental creep compliance modulus , calculated over a constant
“virtual” load step , is not affected by the same approximation.
Figure 4.34 shows both the relaxation and creep compliance moduli over time; to
ease the comparison among them and dimensionally uniform the graph, the
inverse of the creep modulus is presented. Finally, note that the time axis has
now been converted to seconds, a choice which will also be carried on for all
following analyses (cfr. §5.1.2).
Figure 4.34 – relaxation and creep compliance moduli for the case-study
Final results for all tests are displayed in Figures 4.35 to 4.37. Values are displayed
in a log-log scale plot, helping to compare quantities ranging among several orders
of magnitude.
91/237
Chapter 4 – Experimental Tests
Figure 4.36 – a) relaxation and b) inverse compliance moduli for DG41 interlayer†
Figure 4.37 – a) relaxation and b) inverse compliance moduli for PVB interlayer†
______________________________________
† Dotted lines are the results for each different tested specimen (cfr. Table 4.5). For each
temperature, the mean value is presented as a thick plain line.
92/237
Chapter 4 – Experimental Tests
31 Typically, around 120-150°C and 10-15 atm, depending on the manufacturer and
interlayer material [93,347].
93/237
Chapter 4 – Experimental Tests
and §4.1 with further discussion in §5.1). While the increase in tensile strength
brought by such performance-enhancing processes is unquestionably beneficial for
ULG, the fragmentation after glass failure might not be advantageous from a
structural point of view. Depending on several factors that will be discussed in the
following, Partially Damaged- and Fully Damaged Laminated Glass (PDLG and
FDLG) can still withstand significant loads prior to the final collapse.
Studying properties of progressive damage among a LG element is significant for
assessing safety and performance for Ultimate Limit State (ULS) and Serviceability
Limit States (SLS) as well. As a matter of fact, it is common design practice to
include sacrificial glass plies within those LG elements that bear the most
responsibility or those which are more likely to be exposed to accidental failure.
Foraboschi in [42] provided a definition of the characteristics of sacrificial plies:
[The sacrificial ply] collects the live loads [and] is allowed to fracture in service. [The
sacrificial ply] is redundant to the serviceability and ultimate limit states of the LG
member and is included only to allow for the possibility of accidental glass breakage
due to contact damage. The sacrificial ply may be kept in service also in fractured
condition; it has to be substituted only when it has lost the capacity of shielding the
others layers.
4.3.1 PURPOSE OF TESTS
Interlayers are believed to have a negligible influence over the response of in-plane
loaded ULG elements, due to the vast mechanical differences between the two
materials (i.e. the glass is generally orders of magnitude stiffer compared to any
interlayer). Nonetheless, stiff interlayers are known to enhance significantly the
post-failure response [102]. The goal of this experimental campaign 32 , is to
investigate and quantify the effects of the tension-stiffening (TS) mechanism which
arises once one or more glass plies break. The insight that kick-started the
investigation is that glass fragments adherent to interlayers can still provide a
beneficial effect on the global properties of damaged LG.
Architectural LG elements can be designed to bear transversal loads (e.g. glazing,
roofing, etc.) or loads in their own plane (e.g. beams, columns, etc.). From a
structural point of view, elements which bear greater load-carrying responsibilities
are generally the latter. As a safe design of in-plane loaded LG could benefit from
better understood and quantified levels of resilience, a focus on such elements was
chosen for this research33.
Structural LG elements are often comprising several – if not only – fully tempered
glass plies. Tempered glass undergoes thermal or chemical treatments which
increases its strength and influences the breakage texture, producing a denser and
more widespread crack patterns compared to annealed of heat-strengthened glass
(cfr. §4.3). While residual properties of damaged LG elements ultimately rely upon
the strength of the adhesion bond and the intrinsic properties of interlayer
32
Presented studies are actually part of two separate campaigns with similar goals, carried out in
different time spans: a first campaign on SG-laminates was performed in summer 2015 and a
second one on DG41- and reinforced DG41-laminates in summer-fall 2016.
33
Research on post-failure response and progressive damage regarding transversally-loaded
elements is ongoing at the Structures Laboratory of the DICEA by the time this manuscript is
being completed. More information on the results of those tests and comparisons with the ones
presented within this manuscript can hopefully be expected soon, within dedicated research papers.
94/237
Chapter 4 – Experimental Tests
materials [104], as a general rule of thumb, the bigger size of the glass fragments,
the better the overall post-failure performance [77]. Performance analyses of
damaged LG made with tempered glass plies can therefore be regarded as a worst-
case scenario.
4.3.2 DESIGN AND SETUP
Tests were designed to assess the post-failure behaviour of simply supported SG-
laminated beams (type A), DG41-laminated beams (type B) and reinforced DG41-
laminated beams (type C). The “reinforcement” within type C specimens consisted
in 0.2 mm thick polycarbonate (PC) foils between two layers of DG41. A total of
eight specimens were assembled (Table 4.6), each using three 10 mm thick
tempered glass plies. Figure 4.39 shows details of specimens with nominal sizes
and thicknesses.
Figure 4.40 illustrates the design damage progression: first failure of the central
ply, then the lateral plies, one at a time. Rather than reproducing the most common
95/237
Chapter 4 – Experimental Tests
failure mode within real structures34, this modus operandi was chosen to safeguard
the instruments, allowing for the collection of a larger amount of data. Moreover,
this procedure allowed for studying a symmetrical PDLG configuration, which is
convenient for reasons that will be discussed throughout the critical analysis and
discussion of experimental results.
Figure 4.41. shows the two experimental setups which were designed around an
MTS hydraulic press, used to apply vertical loads in three-point bending tests. To
allow specimens to be simply supported, a HE300B steel beam was inserted within
the press columns, to serve as a base for steel rollers, in turn used as punctual
supports for the specimens. Vertical struts were needed to provide additional
vertical support, while lateral torsional restraints were installed for safety reasons:
preventing damage and injuries in the event of sudden lateral-torsional buckling35.
As the manufacturing process did not allow for a perfect glass plies alignment, two
solutions were adopted to rectify the contact surfaces: at both ends, aluminium
pads were placed on top of the steel rollers and a stiff epoxy resin was used to fill
the interface between the pads and the glass surface; at midspan, where the vertical
load is applied, a third aluminium element was inserted between the top of the
specimen and the hydraulic jack. A thin lead foil under the central aluminium pad
was used to compensate for small plies misalignments and help to evenly
distribute the applied load.
As illustrated in Figure 4.44 for type A and Figure 4.45 for type B&C specimens,
several instruments were used to monitor the LG response along the different
stages.
Piezoelectric accelerometers were installed on the specimens and the setup for
dynamic tests, to detect the response of impulse-induced vibrations.
For bending tests, vertical displacements were monitored at midspan (F1) and both
ends (F2 and F3) with LVDT sensors. For FDLG elements, the midspan transducer
was replaced with a cable potentiometer transducer, to cope with the higher sag
expected.
34 Which are likely to have the external plies fail first, due to local contact actions.
35Apart from asymmetrical failure configurations, when the lateral sag didn’t allow for the
beam to rest on a vertical position by itself, torsional restraints were intentionally distanced
from the specimen surface, allowing it to move freely.
96/237
Chapter 4 – Experimental Tests
Figure 4.41 – lateral setup scheme for a) type A and b) type B or C specimens
97/237
Chapter 4 – Experimental Tests
Figure 4.43 –pads used for the setup interfaces (on stage II PDLG specimen B2)
98/237
Chapter 4 – Experimental Tests
99/237
Chapter 4 – Experimental Tests
100/237
Chapter 4 – Experimental Tests
Table 4.7 lists all tests designed to investigate the effects of progressive failure on
LG beams36. Individual tests are described in the following.
I-1 Dynamic response test: multiple vertical impulses are applied at time
intervals at midspan and L/4 37. Accelerometers data is recorded.
36
As mentioned when introducing these experimental analyses, results discussed in this paragraph
and listed in Table 4.7 are referred to tests carried out during two distinct experimental campaigns.
Therefore, one should not be surprised observing that type B and C specimens are listing a larger
number of tests (and carry a larger measurement equipment), as those investigations could benefit
from the experience gathered during the analyses of type A specimens in the first campaign.
37At the beginning, impulses were given with a rubber hammer. However, this technique
proved to be sub-optimal during the first tests, for the formation of a feeble sonic wave
propagating through the beam and the setup, affecting the accelerometers accuracy. After
several attempts, impulses given with a swift palm hand strike showed to yield the best results, this
technique was then used for all other tests.
101/237
Chapter 4 – Experimental Tests
I-2 Three-point bending test: with hydraulic jack speed , a load ramp is
performed up to the design maximum load, held for two minutes, and an unload
ramp at the same speed back to zero load. Repeated twice. Design maximum load
is for type A and for type B and C. Strain data are collected from the
strain-gauges, applied load from the load cell, and vertical displacements at
midspan and both ends through LVDT transducers.
II-1 Failure of the central glass ply: failure of the central glass ply is induced with
a glass-breaker pick and a hammer. Strain gauges data are recorded during failure
and for 10 minutes after failure.
II-2 Long-time strain recording: while the specimen is left undisturbed, resting on
the supports at both ends, strain-gauges data is recorded continuously for 24 hrs
after the glass ply failure. In the same room, an additional strain-gauge is installed
on a separate glass surface in the same room, to allow for a posteriori compensation
of thermal effects.
II-4 Three-point bending test: similar procedure of test I-2, the design maximum
load is now for type A and for type B and C. Strain-gauges data are
now only collected from instruments installed on non-broken glass plies.
III-1 Failure of one lateral glass ply: the first lateral ply is broken with the same
procedure of test II-1. The first lateral ply designated to be broken is the one not
equipped with the strain gauges to .
III-2 Measure of the lateral sag: a caliper is used to measure the lateral sag at
midspan, induced by the shattered glass expansion in an asymmetrical failure
configuration. A thin spring steel cable is used to collimate the ends of the beam
during the measurement.
III-5 Three-point bending test: Same procedure than test I-2, the design maximum
load is now for both type A, B and C.
IV-1 Flexural failure: failure of the last glass ply is induced with a continued three-
point bending test. Lateral torsional restraints are placed close to midspan and
used to hold the cross-section in vertical position.
IV-2 Dynamic response test: similar procedure of test I-1, except vertical impulses
are applied at intervals.
IV-3 Three-point bending test: similar to the procedure of test I-2, the design
maximum load is not set a priori, while a maximum sag of is imposed to
safeguard the measure instruments. Strain-gauges data is replaced newly installed
“omega” transducers in a position close to midspan. Before testing, specimens
were supported on multiple points (two at least) with wooden wedges, to prevent
the onset of significant viscous deformations. Wedges have been removed
moments before the beginning of tests.
IV-4 Repeated three-point bending test: same procedure of test IV-3, performed 5
months after the previous test.
102/237
Chapter 4 – Experimental Tests
4.3.3.1 UNDAMAGED LG
The undamaged mechanical properties of glass have been drawn from
experimental analyses, specifically, by interpretation of tests I-2 (Table 4.7).
According to European Standards, the elastic modulus of glass is allowed to ranges
between and [38,278]. Experimental values of the Young modulus have
been obtained assuming that interlayers – in virtue of their very low modulus - do
not give a significant contribution to the response of in-plane loaded LG beams.
For a simply supported element of span , composed of glass plies, of thickness
with moment of inertia about the strongest axis, with a load at
midspan, the elastic modulus of glass can be obtained as:
Experimental values of have been evaluated on the second ramp load plateau of
tests I-2, when the design action is held constant for two minutes. Table 4.8
displays experimental values of the elastic modulus for glass. Tests on materials
coming from different manufacturers yielded comparable results.
While the interlayer modulus was neglected, the presence of the PC reinforcement
for type C specimens has been considered with an elastic modulus of
[66,279]. Figure 4.47 shows how the experimental modulus defined with this
procedure is consistent over the whole loading plateau; final listed in Table 4.8 for
all specimens are referred to the experimental point in the middle of the plateau. It
is worth pointing out that such values must be looked at as average moduli for the
three plies comprising each specimen.
Figure 4.47 - glass modulus evaluation over the 2nd load ramp for type A specimens
Data recorded by accelerometers in dynamic tests I-1 have been analysed using a
one-sided Fast Fourier Transformation (FFT), to evaluate the frequency content. As
an example, Figure 4.48 shows the result of this procedure for the undamaged
specimen A2.
The first natural frequency is easy to be found at , with a mutual effect
also measured by accelerometer placed on the setup. Similarly, the effect of the
natural frequency of the setup, and its effect on all three instruments installed on
the specimens, are found around . At higher frequencies, between and
103/237
Chapter 4 – Experimental Tests
After the first campaign, evidence showed that accelerometers and were
giving redundant results. Therefore, as shown by comparison of Figure 4.44 and
Figure 4.45, a slightly different instrument choice was made for the second
campaign on type B&C specimens: former instrument a3 was removed, while
was installed at , rather than . This choice was motivated by the fact that,
while the second modal form displacements are not massively different between
and , one can hope to glance at the third modal form. Unfortunately, this
approach did not give satisfactory results in practice.
Figure 4.49 shows the same analysis for a specimen of type B, investigated within
the second campaign, with the new instrument configuration. While insight on the
second modal form could still be gathered for most specimens, those results
seemed overall less reliable with respect of the first campaign results. Ultimately,
the data analysis focused on results for the first modal form only, which was
always very easily identifiable.
By analysing data from strain-gauges on the monitored cross-section close to
midspan (A-A in Figure 4.44 for type A specimens and Figure 4.45 for type B&C
specimens), information can be gathered on the likely position of the neutral axis
throughout bending tests.
104/237
Chapter 4 – Experimental Tests
Figure 4.50.a shows recorded data for the two load ramps, according to the
procedure for test I-2. Strain over the height of the cross-section has been regressed
with a linear interpolation for the whole extent of tests, showing how the neutral
axis is consistently located at mid-height. Figure 4.50.b shows the linear
regression38 at the center of both plateaus (highlighted in Figure 4.50.a). Results
show that the interpolated position of the neutral axis is exactly at mid-height, as
expected for a symmetrical cross-section made of an isotropic material.
38A linear regression is consistent with the Euler–Bernoulli hypothesis that plane sections
remain plane.
105/237
Chapter 4 – Experimental Tests
Figure 4.51 - comparison of long- and short-term measured strain on the lateral plies
after the failure of the central ply
Results in Figure 4.51 and Figure 4.52 show readings for those instruments which
are installed on plies that remain intact after the imposed failures. Strain gauges on
shattered glass plies may be affected by local effects which are arguably very
difficult to interpret. The maximum extent of induced strains appears to be
consistent among the specimen types B&C. Nonetheless, the double interlayer
thickness of type C specimens seems to play a role, allowing for a faster
stress/strain relaxation in the symmetrical PDLG configuration (Figure 4.51.b and
Figure 4.51.d).
106/237
Chapter 4 – Experimental Tests
Figure 4.52 - comparison of long- and short-term measured strain on the intact lateral ply
after the failure of the opposite lateral ply (the central one is already broken)
Figure 4.53 - short-term effects after central ply failure on all specimens type B and C
107/237
Chapter 4 – Experimental Tests
When a lateral ply fails, instruments on specimens of type C with double interlayer
thickness indicate that the curvature (hence the stress and strain) slightly increases
over time, at least for the first 24 hours. Figure 4.52, for type B specimens with
regular interlayer thickness, suggest that instantaneous effects of the lateral glass
failure are more consistent over time, except for instruments e11 and e15 installed
very close to the edges.
Figure 4.53 provides a graphical comparison of short-term effects of symmetric
glass failure for all specimens of type B and C, showing that results are consistent
among similar elements.
Dynamic tests have been analyzed similarly to those described in the previous
paragraph §4.3.3.1. Figure 4.54 shows spectral density results for accelerometers
(this time, on a linear scale, highlighting the spectrum peaks even more compared
to Figure 4.49). Results for different specimens do not display substantial
differences, a final synthesis of experimental natural frequencies is shown in Table
4.10 (p.113).
Bending tests II-4 and III-5 on PDLG, continue to showcase a substantial linear
elastic response, as did test I-2 on ULG specimens. Figure 4.55 shows curves for
applied load over vertical displacement at midspan: evidence suggests that no
permanent effects are induced, as the load and unload paths are essentially
undistinguishable (plots in Figure 4.55 are displaying both).
108/237
Chapter 4 – Experimental Tests
It is worth pointing out that the elastic response of PDLG elements vastly exceeds
the one expected only accounting for the undamaged glass plies contribution. As
this increased stiffness is ascribable to the collaboration of fractured elements and
interlayers, further analyses will be carried out in the following §5.2.
As observations of a rise in the neutral axis position in PDLG bending tests were
expected, strain-gauges data from the monitored sections were analyzed similarly
the procedure used to verify the isotropic behaviour of ULG specimen (cfr.
§4.3.3.1). Figure 4.56. shows how strain measurements gathered in tests II-4 and II-
5 (Table 4.7) were linearly regressed to interpolate the likely position of the neutral
axis.
Figure 4.56 - a) snapshot of interpolated strain distribution over the cross-section height at
and b) progression of the neutral axis position for increasing load
For stage II, when specimens are partially damaged with only the central ply
fractured, results show that the neutral axis continues to be steadily consistent with
the centroid, as it was previously observed for undamaged specimens in stage I.
Conversely, as damage progresses and one of the lateral plies has failed, evidence
from bending tests in stage III shows that the neutral axis rises even so slightly
towards the top of the beam, even for moderate loads. The evolution of the neutral
axis depth can be understood at best by looking at results of test IV-1, only
performed for specimen C2 (cfr. Table 4.7 p.101). Figure 4.57 shows that the neutral
axis settles rapidly upon a position slightly above of the centroid and does not
move significantly up to the ultimate failure39.
Table 4.9 lists interpolated positions of the neutral axes for all specimens during
bending tests I-2, II-4 and III-5, for the corresponding damage states 40. Despite
39
It is worth pointing out that, due to the lateral sag, torsional stress components are induced
alongside pure flexural ones. While this effect is not being recorded by the strain gauges installed
on section A-A (e5 to e8 in Figure 4.56), as they are designed to monitor strain in the direction of
the specimen’s axis, the sum of torsional actions to flexural ones may have an effect in anticipating
the ultimate failure.
40
As the bending tests procedure was designed to perform two consecutive load ramps, thus
including two loading “plateaus” where the load is held constant, individual values shown in Table
4.9 for specimens A1 to C3 have been computed as the averages between interpolations on both
plateaus.
109/237
Chapter 4 – Experimental Tests
some uniquely unexpected results (e.g. stage III for specimen B2), values of the
average rise over the cross-sectional height show remarkable consistency over
different specimens.
41 This choice was motivated by the attempt to test “undisturbed” FDLG specimens in the next
phase. Hence, a locally induced failure was preferred to a violent bending one, which was thought
to influence the FDLG specimens with different crack patterns and/or a pre-stress in the interlayers.
42
Operations included the installation of omega transducers, removal of strain-gauges cables and
full instruments check.
43
To prevent self-weight from inducing significant viscous deformations before the start of the
data acquisition, specimens were provisionally supported on three points at L/4, L/2 and 3/4L with
wooden wedges. Wedges were removed moments before the start of the tests.
110/237
Chapter 4 – Experimental Tests
Figure 4.58 - normalized moment over midspan sag for all specimens
where is the total mass of the specimen (glass plies + interlayers), is the
standard gravity, is the total span and the measured sag at midspan.
Ageing effects on SG-laminated specimens are shown in Figure 4.59 with a
comparison among tests IV-3 and IV-4 (cfr. Table 4.7 p.101). On specimen A1, the
IV-3 test has been repeated twice, to assess that the maximum load applied at
midspan during the first ramp (2.5 kN) was not yet significantly affecting the
mechanical properties of the fractured-glass-interlayer composite. Tests IV-4 were
performed on the same specimens, after a five-month storage interval at indoor
ambient temperature (20-25 °C) and RH.
Figure 4.59 - ageing effects for FDLG specimens for bending tests
Omega transducers data were used to compute average strain on top and bottom
of the specimens (cfr. Figure 4.44, Figure 4.45 and Figure 4.46 p.99-101). Figure 4.60
shows data for the specimen B2: measured strains on both sides of the mid-span
section allow for interpolating the likely position of the neutral axis, with the
average value yielding the final best estimate.
111/237
Chapter 4 – Experimental Tests
Figure 4.61 shows final results for all specimens. A clear majority of the cross-
section is consistently below the neutral axis, bearing tensile stresses. The
translation equilibrium along the direction of the axis is provided through a
smaller portion of the cross-section near the top of the beam, whose depth appears
to be approximately of the total height.
Figure 4.61 - evolution of the depth of the neutral axis in all FDLG specimens 44
44
Data for specimen B1 is not shown, as laboratory observations proved that omega transducers
did not work properly during that specific test.
112/237
Chapter 4 – Experimental Tests
After the conclusion of tests, Type A specimens were stored for five months inside
the laboratory, in a protected environment, and cut later to extract the smaller
material specimens described in §4.3.4. On the other hand, FDLG specimens Type
B & C have been stored outside the laboratory after the last bending tests, under
natural environmental conditions45. After six months, a remarkable difference in
the long-term effect was observed (Figure 4.62): all specimens were supported on
multiple points with bricks (4 or 5), while Type B laminates had completely lost the
original shape, the polycarbonate reinforcements in Type C specimens allowed
them to hold their original shape. Moreover, the limited deformation of the latter
also reduced penetration of moisture and rainfall inside cracks, allowing for a
significantly higher adhesion of fragments46.
45
About 2-4 hours of sun exposure per day, moderate exposure to wind and moisture, seasonal
precipitations and icing, with temperatures ranging approximately between -5 °C and 40 °C.
46
While fragments of Type B specimens could be removed with little effort, just rubbing a hand
over the beams, fragments of Type C beams were still firmly adherent to the interlayer and could
not be easily torn apart.
113/237
Chapter 4 – Experimental Tests
Figure 4.63 - specimens extracted for tensile and compression tests on the fractured-glass-
interlayer composite material
Figure 4.64 - setup for a) tensile and b) compression tests on material specimens
To preserve the integrity of the measuring instruments, tensile tests were initially
performed up to relatively small strain extents47; this choice allowed for repeated
tests on the same specimens, to gradually investigate the damage progression up
to slightly higher strain values ( ). Moreover, two different strain rates were
tested: (slow) and (fast), respectively.
Figure 4.65 shows experimental curves of tensile tests on material specimens.
Stresses were calculated assuming that the cross-section reduces to only the two
ionoplast interlayers. Similarly to observations in §4.1.2, experimental evidence
highlighted once again how the loading rate is an important parameter: the faster
the strain rate, the higher the stress peak. In addition, for strain cycling multiple
times between zero and appears not to give a major decay on the stress
peak and only a minor one for the elastic modulus, slightly lowering the initial
stiffness of this composite material.
47
about 4%, which is small compared to the full elongation potential of the interlayer material
itself, which is about 350% as previously shown in Figure 4.4 §4.1.2.
114/237
Chapter 4 – Experimental Tests
When the adhesion bond is compromised, the interface had been observed to
produce a slightly different gleam with reduced transparency; a contrast which can
been be further incremented with post-processing photography. Aiming to capture
the debonding progression of the glass-interlayer interface, a photographic follow-
up was performed. The setup was put in place during one of the “slow” tensile
tests on the composite FDLG material, similarly to the one described in §4.1.2 for
SG interlayers only.
Figure 4.66.a shows that a small initial debonding is already present in the
specimen, even before the start of the test. On the specific fragment selected for the
analysis, the area interested by the loss of adhesion bond started to expand rapidly
even at low strains, about as shown in Figure 4.66.b, and affected about 40%
of the total interface around . From there on, debonding on the selected
fragment was stable up to the end of tests.
Compression tests on FDLG material samples were also performed by imposing a
controlled displacement generating a strain rate of . Tests were carried
out to the ultimate failure of the samples, resulting in a violent destructive
procedure48 which didn’t allow for repeated tests, as it was instead for tensile ones.
Figure 4.67 shows the results of experimental compression tests: the mean stress
shown in the vertical axis is calculated on the total cross-sectional area, composed
by the two layers of SG and adherent glass fragments 49.
48
Final collapse happened with as a sudden explosion of the FDLG material specimens. The
explosion instantly and violently scattered a large number of glass fragments out of the specimens
body, reducing the load bearing capacity to a fraction of the one measured before the collapse.
49
Glass fragments alone are unable to bear any tensile actions. Therefore, assuming that tensile
stresses are concentrated in the interlayer seems a rather logical choice. On the other hand, glass
115/237
Chapter 4 – Experimental Tests
fragments are able to bear very significant compressive actions by exchanging contact forces one
another. Anyhow, the mechanic response measured during compressive tests on FDLG is greatly
influenced by the presence of interlayers to which the fragments adhere to. Without the interlayers
contribution, fragments would be unable to provide the same response on their own. With that in
mind and for sake of simplicity, an equivalent stress could be calculated on the interlayer surface
alone, similarly to tensile rests. Such a radical choice can be motivated by the ultimate goal to
provide a simple, macroscopic, overview of this composite material performance. This approach is
used in the following to model the stress distribution along the cross section of the beams in the
bending tests described in the previous paragraphs.
50
Experimental values are given for a three-layered LG arrangement, assembled with two 1.52
mm thick SG interlayers and 10 mm thick glass plies. It is likely that different arrangements may
yield different experimental values for similar tests, although the order of magnitude is not
expected to vary dramatically.
116/237
Chapter 5 – Critical Analysis and Discussion
5 CRITICAL ANALYSIS
AND DISCUSSION
The scientific community has yet to agree on reliable predicting models for
LG. Due to many factors involved, experimental tests and empirical models
are to this day the main toolkit at a designing engineer’s belt. Results
gathered within previous investigations are now examined in depth, aiming
to share more light on the subject.
Experimental tests in Chapter 4 have been divided into three parts. From a
thematic point of view, long-term viscoelastic properties and short-term tensile
response tests deal about interlayer properties alone, with limited interest on the
mechanics of LG composites. On the other hand, tests on structural LG beams
bypass the investigation on the intrinsic behaviour of interlayer, aiming instead to
a better understanding on the full-scale, post-failure performance.
In the following, topics will be conveniently united based on a thematic point of
view: a first part is dedicated to interlayer mechanics, merging the results
described in paragraphs §4.1 and §4.2, while the last part dedicated to analyses on
full-scale LG introduced in §4.3 and post-failure effects.
5.1 INTERLAYER PROPERTIES
Short-term tensile results are more suited to model effects such as impact loading,
blasts and generally any effect which takes place over a short amount of time (i.e.
seconds, minutes); those interlayer properties are discussed in the following
paragraphs. On the other hand, results and models provided for long-term
properties of interlayers are the core of the research and essential for the design of
LG elements; these are dealt with in §5.1.2.
5.1.1 SHORT-TERM PROPERTIES OF INTERLAYERS
Tests carried in out in §4.1 yielded a vast amount of data on tensile properties of
three among the most commonly used interlayer materials: PVB, SG and DG41.
Existing standard guidelines have been followed when possible and some early
results have already been discussed accordingly (cfr. §4.1 and sub-paragraphs).
In the following, an in-depth analysis of those results is carried out, a discussion
articulated around the experimental outcomes, with an eye on the state-of-the-art
research and the regulations in force. True stress-strain diagrams are produced by
interpreting the numerical results and the photographic follow-ups presented (cfr.
§4.1.2). Hyperelastic models are then calibrated and discussed for the supposedly
reversible network entanglement effects. Together with a novel interpretation of
the intermolecular cross-links response, those hyperelastic models are used to yield
a final decomposition according to the theory proposed by Haward and Thackeray
[76]. Finally, a discussion on the yielding (or “pseudo-yielding”) point for different
thermoplastics is undertaken, concluding with the proposal of a more universal
approach to the problem.
117/237
Chapter 5 – Critical Analysis and Discussion
Values of the Poisson ratio for climatic conditions corresponding to those for
tensile tests (cfr. §4.1) on the three investigated materials are shown in Table 5.1.
True strain is defined as [164]. While for materials with
homogeneous deformations such as PVB this equation can be applied in a
2
where and are the length of the segment during tests and its initial length
respectively. Figure 5.1 compares the strain readings from the testing device,
calculated as the average among the measure points in Figure 4.9, to
photogrammetric strain measurements resulting from the time-lapse. Results show
that the relative error in the early part is (i.e. for small average strain) in very
significant, while it seems to be consistently under 5% after the engineering stress
curve crosses once again the pseudo-yielding threshold (cfr. §). Complete results are
in Appendix E.
The maximum strain obtained by photogrammetric measurements is the best
estimate of the maximum strain within the body region. Although those
measurements still represent a local average, one which could be improved with
more tightly-packed marks (cfr. Figure 4.9 p.71), if this local average is a good
estimation of the effective punctual elongation, true strain has been computed
accordingly.
1As a general rule of thumb, the value of the Poisson’s ratio for polymers is known to be is
positively correlated with time and temperature [240] (i.e. the lower the temperature and/or elapsed
time, the lower the ratio). Greatest variations are observed across the glass transition temperature.
2
where it can be assumed that in the narrow area average and maximum strain are consistent.
118/237
Chapter 5 – Critical Analysis and Discussion
If the Poisson ratio for the three materials is consistent over large deformations 3
with the values given in Table 5.1, knowing the applied load as a function of
the true strain , true stress can now be calculated over the reduced cross-
sectional area 4.:
Results in Figure 5.2 show how formerly observed softening regions in SG and
DG41 were ascribable to the engineering representation of the curves, while true
stress-strain diagrams exhibit a consistent hardening response.
3
cfr. footnote 8 page 86.
4
If considering interlayers as incompressible materials (i.e. ) is an acceptable
approximation, the formula simplifies into .
119/237
Chapter 5 – Critical Analysis and Discussion
5
For instance, a neo-Hookean model has been used in [156] for studying true-stress strain curves
of a compressible thermoplastic material, while engineering stress-strain relations on
incompressible polymers have been studied by several other authors using various models
[222,226,230,234,282].
6
for approximately in SG and in DG41 tests (Figure 5.3).
120/237
Chapter 5 – Critical Analysis and Discussion
Experimental tensile tests results7 (cfr. §4.1.2) have been compared with several
known hyperelastic models8 (cfr. §3.1.2); while most among those could reliably
describe the shape of the engineering stress-strain curve for PVB, fitting the results
for both DG41 and SG proved to be a more challenging task.
Table 5.2 displays the different accuracy of four among the most commonly used
HE models to fit experimental results. The Yeoh model proved to be the simpler
among the two most reliable ones. Figure 5.3 shows final results for Yeoh
incompressible model with coefficients in Table 5.3.
Table 5.2 – accuracy of hyperelastic models for first and final parts of the curves
experimental data fit efficacy
Number of
Model SG DG41 PVB
coefficients
first final first final first final
Neo-Hookean 1 C C C C B C
Mooney-Rivlin 2 C B C A B A
Yeoh 3 C A C A A A
Ogden (N=3) 6 C A C A A A
A: good (average relative error lower than10%) - B: average (average relative error between 10
and 90%) - C: limited or very limited (average relative error greater than 100%)
7
Even though the engineering application and modelling are the ultimate focuses of this
manuscript, for the following analyses the true stress-strain curves have been preferred to
engineering ones (cfr. §5.1.1.1). True response diagrams are believed to be more suited for
understanding the microscopical origins of the experimentally measured macroscopic effects.
Furthermore, one can always go back easily to the engineering curves knowing the true
correlations, while the opposite is generally harder as it requires additional data.
8
The following models have been used during the analysis: Neo-Hookean (1943) [348], Mooney-
Rivlin (1951) [231,235], Ogden (1972) [236] and Yeoh (1990) [349].
121/237
Chapter 5 – Critical Analysis and Discussion
with strain and strain rate and , temperature and thermo-mechanical history
parameter . While the results discussed in this paragraph are obtained in
controlled and consistent climatic conditions, for the known dependency of
thermoplastics on the temperature , it seem safe to assume this may be a relevant
parameter for both and . The dependency on thermomechanically-reversible
aging effects is believed to be related to variations in density of intermolecular
actions and has not been proven to influence poorly cross-linked polymers [156].
The proposed decomposition aims to separate reversible effects from permanent
changes. Consistently with the above-mentioned studies, the following hypotheses
are made:
1. large strain effects are governed mostly by the network response and
short-range deformations primarily refer to intermolecular effects,
2. conformational changes of the network are fully-reversible,
3. the break of intermolecular bonds is responsible for non-reversible effects.
9
This choice is consistent with the fundamental interpretation of the hyperelastic effect, expecting
a greater influence the network response for high strains [10,14].
122/237
Chapter 5 – Critical Analysis and Discussion
Looking at the experimental results of true stress-strain diagrams in Figure 4.3, one
can observe that the two significantly cross-linked polymers, SG and DG41, exhibit
an initial response step which appears to be dependent on the strain rate.
Following both the HT and KD approaches, one can assume that, on a true stress
over strain diagram, intermolecular interactions will drive most of the initial
material response. Moreover, conclusions and models proposed by KD suggest
that intermolecular interactions may produce a steady response over time (the so
called rejuvenated yield stress cfr. Figure 2.5). Having excluded the presence of
softening branch in the true constitutive response, and considered the shape of the
experimental curves, it makes sense trying to model the initial materials reaction
with a non-decreasing function which tends to a finite limit. To this end,
exponentially asymptotic functions as simple as proved to be a
powerful and versatile tool: by means of just two coefficients and needed to
calibrate respectively the amplitude and rate of stress over strain functions.
Nevertheless, by further investigating the first part of the response curves,
experimental results showed that SG displays a very clear and abrupt slope
change, in the first part of the stress-strain diagram, as illustrated in Figure 5.4.c10
(this phenomenon, on a smaller scale, is also visible for DG41). On the other hand,
PVB only displays a single initial “load-step” before transitioning to a stress
plateau with minor rising slope, as shown in Figure 5.4.b (i.e. there is no point of
sudden slope change, but rather a gradual diminishing).
While all of this is – knowingly - speculative and still needs deeper investigations,
especially on a chemical level, such phenomenon of abrupt response changes in
early stages of deformation can be traced back to the presence of different chemical
interactions: primary interactions, weaker but very densely packed, may be
accountable for the first - very stiff - material response, limited to small stresses.
Besides, stronger secondary chemical interactions, able to withstand higher stresses
but concurrently scarcely distributed within the polymer matrix, are activated by
larger configurational changes. Following these assumptions, the intermolecular
response can be represented by a sum of terms with coefficients , 11:
10
It is not excluded that different materials, not included in this campaign, may show not just one
but multiple abrupt slope changes along the experimental curves.
11Ideally, each element of the series describes the contribution of a specific homogeneous
group of intermolecular bonds (cross-links). Nonetheless, the theoretical explanations
regarding chemical bonds forming within the three investigated materials within the
manuscript are speculative and lack the support of any chemical analysis. Future researches
123/237
Chapter 5 – Critical Analysis and Discussion
This phenomenological model for small strain will be referred to in the following
as the intermolecular series (IS) model.
Like in the example of Figure 5.4.b, PVB only showed one initial – and very
limited – load step, followed by the progressive rise of the response ascribable to
the network entanglement effect. While the network effect has been proved to be
well represented with an hyperelastic model, the initial response can be modelled
by means of one term of the IS model (Figure 5.5.c). On the other hand, the initial
phases DG41 and even more so SG showed some abrupt slope change mentioned
above of the curves for very small strains, as shown in Figure 5.4.c. With the
proposed model, such a slope change can be represented with a larger number of
terms in the series. For both materials, with only one slope change, the number of
terms needed turned up to be limited to two. Table 5.4 shows coefficients ,
manually calibrated to model the experimental curves. Overall results and close-
ups are shown in Figure 5.5.
on this matter, either to support or refuse the proposed theories, are deeply encouraged.
124/237
Chapter 5 – Critical Analysis and Discussion
With this convenient way for describing the initial response, which is assumed to
be ascribable to intermolecular interactions, we are able extract by subtraction the
effects of the entanglement of molecules network (dashed lines in Figure 5.6):
With IS model takes care of the description of the first part of experimental results,
the shape of the remainders in Figure 5.6 for the three materials suggest that
all can be well represented by numerous hyperelastic models. Among all the well-
known and widely used models that can have a fair shot in reproducing the shape
of these curves (cfr. Table 5.2), it is worth noting that the aforementioned simple
cubic equation , used by Lenci and Clementi in 2014 [284] could be a
simple option worth investigating, ultimately yielding to a complete description of
the materials response in the form of a generalized response (GR) model:
125/237
Chapter 5 – Critical Analysis and Discussion
Inverting the equations in §5.1.1.1 allows to transform the true curves back to the
original engineering representation:
Figure 5.8 and Figure 5.9 show the results using coefficients in Table 5.4. A
schematic recapitulation of the process is provided in Figure 5.7.
Figure 5.7 - scheme of the process used to produce the final models
12
Comparing results in Figure 5.36.a/b/c also suggest that the hyperelastic network response may
not always be significantly dependent on strain rate . For instance, DG41 results in Figure 5.36.b
show that the remainders for different strain rates are very consistent one another. Although
this might be a coincidental combination for multiple factors, this fact is worth investigating with
future experiments. The only logical deduction at this stage would be not to give for granted the
dependency of the model on all supposedly influential parameters.
126/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.8 - model true stress-strain results for a) SG b) DG41 and c) PVB
Figure 5.9 - model engineering stress-strain results for a) SG b) DG41 and c) PVB
127/237
Chapter 5 – Critical Analysis and Discussion
13
The same phenomenological observation suggested in §5.2.2 the use of a convenient
exponentially asymptotic mathematical formulation for the IS model.
14
For very small strains, engineering and true stress-strain diagrams are almost identical; the
engineering representation is chosen to match the approach used in the standards [66,264].
15
The coil conformation in undisturbed conditions is the state of lower entropy [263].
128/237
Chapter 5 – Critical Analysis and Discussion
A “yield point” for plastics is mentioned in the standards [66,264], but its very
definition seems only applicable to a small subset of materials, those for whom a
local maximum [is] observed during a tensile test16. This definition can be generally
applied to rigid plastics (Figure 5.10.a). Nonetheless, the same standards are very
often used to test semi-rigid and thermoplastic materials (cfr. §1.1.2) which to not
show a peak in the curve, but only a curvature change (Figure 5.10.b) or a flat and
very long plateau (Figure 5.10.c).
For engineering practice and a LG design perspective, the definition of a “yield
point” for interlayers must cope with the fact that, for such thermoplastic
materials, a distinction between elastic domain and permanent deformation is not
only a function of stress (like it is, for instance, for metals), but several other factors
come into play such as temperature, time and strain rate; playing a non-negligible
role in defining the threshold from which permanent deformations are significant.
For instance, permanent deformations can be produced in thermoplastics by minor
stresses, in long times [65]. Therefore, the traditional definition of yield seems not
straightforwardly applicable17.
By these remarks, a pseudo-yielding stress limit [16,289,290] has been
mentioned by several authors in the past, when talking specifically of
thermoplastics and elastomers. For instance, the pseudo-yielding limit for PVB has
been recently defined by X. Zhang [16] as:
“[the point] where material modulus changes abruptly”
All things considered, a different concept of yielding looks like a necessary
introduction in this field. Sadly, a clear definition seems still to be missing.
Many researchers showed that thermoplastic materials can exhibit a wide array of
mechanical responses to traditional tensile tests [16,114,153,154,291]. Looking at
various experimental outcomes, the case system seems to be ascribable to the three
typologies presented in Figure 4.35. Looking for a general tool to define this
“abrupt” change in the material response, an arbitrary but simple and generally
applicable criterion to identify the pseudo-yielding limit is proposed: on an
engineering stress-strain curve, the pseudo-yielding limit may be defined as the
point Y where the slope of the tangent line is 10% the secant modulus 18 ,
calculated for 19 .
129/237
Chapter 5 – Critical Analysis and Discussion
Results for the performed tests are given according to the proposed method in
Table 5.5 and Figure 5.11.
Table 5.5 - pseudo-yielding limit for the performed tests
claw displacement pseudo-strain
SG DG41 PVB
velocity point
0,060 0,00044 0.0012
5 mm/min
11,68 2,96 0,10
0,053 0,00038 0.0010
50 mm/min
16,20 7,02 0,16
0,040 0,00030 0.0007
500 mm/min
20,04 12,82 0,37
130/237
Chapter 5 – Critical Analysis and Discussion
5.1.2.1 TIME-SHIFT
By correlating properties of a given interlayer at a set temperature to those in a
different climatic condition, results can be effectively extended beyond the fixed
temperatures investigated in the campaign.
The shape of relaxation and – inverse - creep compliance curves on a log-log
modulus over time space shown in Figure 4.37 suggest that an horizontal shift of
the curves may be used to arrange a single function describing the viscoelastic
properties of the materials [114,115].
Choosing a reference temperature for each material (i.e. a segment of the curve
that will not be moved from its initial position), linear drifts in an -base
logarithmic time domain for the other -branches can be performed. The new time
scale for a branch shifted by a dimensionless factor is:
where is the original time of the branch and in our analysis an base has
been chosen (cfr. §3.1.3.5 and §4.2.4.4). An example is presented for .
Shift factors for creep and relaxation curves have been iterated to obtain the
smoothest functions, a process that will be further clarified in the following. While
this procedure may seem highly arbitrary, results from the final mastercurve - the
one that is produced from the time-shifted data - are tied by a biunivocal relation
to the shift factors used to produce the mastercurve itself (i.e. one cannot be used
without the other, and each version of one defines a unique parent of the second
among a family of possible curves). To clarify this reasoning, the scheme of the
implemented iterative process in this work is shown in Figure 5.13.
The interlocking correlation between WLF and mastercurves implies that the best
shift criterion is to produce mastercurves and shift factor distributions which are
prone to be well interpolated by the relative functions21: a Wiechert-model-based
Prony series and WLF equation respectively (cfr. §3.1.3.1 and §3.1.3.5).
131/237
Chapter 5 – Critical Analysis and Discussion
Shifted Linking
iterate branches procedure
no
Data polish
22
By doing so, the curves can be used as presented, without any shift, if the materials response is
needed at ambient temperature.
132/237
Chapter 5 – Critical Analysis and Discussion
The final choice of the shift factor distributions presented in Table 5.6 is made
according to the iterative procedure described in Figure 5.13. An algorithm was set
up to minimize the sum of square differences between the shift factors
distributions and interpolating WLF equations Γ 23 and the sum of square
differences between the experimental data of shifted branches and interpolating
mastercurves Γ 24. The input of the shift factors distribution was done manually
and so were the iterations by trial and error. For each attempt, the two sum of
square differences were stored recorded and normalized to Γ and Γ , with respect
of the corresponding number of data points 25 . The final choice of shift factors
(answering to the question “satisfying interpolation?” in the scheme of Figure 5.13)
was the one that minimized the combined error Γ δ Γ Γ .
Table 5.6 - shift factors 26
Temp. SG DG41 PVB
0 -4.524 -7.042
10 -2.434 -7.504 -2.076 -3.365
14 -1.367 -4.278
17 -0.688 -1.909
20 (ref.) 0 0 0 0 0 0
23 0.710 1.763
30 1.848 2.786 2.324 4.633 1.886 3.185
35 2.901 4.349
40 4.153 6.075 4.328 7.439 3.685 6.171
45 5.532 7.802
50 7.134 9.813 6.120 9.536 5.281 9.028
55 9.212 12.126
60 11.911 15.072
23
Using the proposed interpolating procedure described in §5.1.2.3.
24
Performed with the DOI method described in §5.1.2.4.
25
For the WLF, the number of data points is the number of different temperatures of the tests, for
the Prony regression the number of data points is the sum of points of the experimental branches.
26
The choice of the displayed shift factors will be clearer with the following discussions and
analyses presented in §5.1.2.3.
133/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.15 - shifted branches for a) exp. relaxation, b) exp. creep and c) averages
134/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.15 shows simultaneously average values of shifted branches for both
relaxation and creep moduli. Several noteworthy elements can be pointed out:
pink branches, both for relaxation and creep, are tests performed at the
reference temperature ( ); therefore, all curves remained on the
vertical of the same time frame (about ) as only the other branches
were shifted back and forward,
overall, the creep compliance modulus appears to be flatter than the
relaxation one (i.e. it is more stretched along the time axis), while values
for both moduli at the same temperature do not vary dramatically,
because of the two previous points, relaxation and creep curves always
cross one another at the reference temperature,
the transition temperature for the three materials, represented by the
maximum slope of the function, appears to be close to for PVB,
for DG41 and for SG,
to a first approximation, in a log-log modulus over time space, the six
curves can be roughly interpolated with two or three straight lines, thus in
linear scales with three different exponential curves.
135/237
Chapter 5 – Critical Analysis and Discussion
The linking function and boundary conditions can hence be expressed in a linear
time over moduli transformed environment, where the linking function is a third-
degree polynomial curve:
136/237
Chapter 5 – Critical Analysis and Discussion
From the linear regression, ordinates and of the target points and slopes
of the curves and were worked out, to be then used as boundary conditions
for the polynomial interpolation. Given:
Coefficients , , and were found by solving the linear system granting the
observance of the established boundary conditions. Numerical results in this
example are given to be later implemented in the general equation:
27
Note that, for explanatory reasons, the representation in Figure 5.18 has been intentionally
simplified by scattering the and expanding the neighbourhoods around the target points for which
the sets are defined. Actual sets generally consisted in thousands of elements, increasingly dense in
the ending parts of the curves, an effect due to the applied transformations.
137/237
Chapter 5 – Critical Analysis and Discussion
The plain black curve in Figure 5.18 is the polynomial solution for the
gap between target points and in the transformed variables , , where the
subscript is used to indicate a linking function. The solution in the original
over environment can be found by applying the following inverse
transformations:
for .
The same procedure has been applied in an analogous manner for all the gaps in
all moduli for the three materials in different intervals.
Dealing with experimental data, inherently prone to a variety of possible
approximations and errors, a critical approach in the selection of the target points
was preferred to a blind implementation of the procedure. Consequently, target
points yielding an overall smoother interpolation of the general trend of the
mastercurve were chosen, rather than slavishly connecting the terminal parts of all
the connected experimental branches. This approach also took into consideration
that the terminal parts of the experimental curves generally have a higher
regularity than the initial ones. This phenomenon is also discussed in §4.2.2 and
further highlighted in Figure 5.15 by the juxtaposition of branches which do not
strictly require a linking function (e.g. and , where the initial part of the
curves prove to be less regular than the following trend). With this approach,
target points and linking functions are presented in Figure 5.19.
138/237
Chapter 5 – Critical Analysis and Discussion
Portions of the linking branches have later been cut, to avoid unnecessary
superposition with the experimental data, in the - small - portions where the
curves were overlapping (experimental data were always preferred to points of the
linking functions).
The final operation on the data at this stage, to enhance both reliability and
computational load of Prony series interpolations, was to re-sample the data sets
among a fixed number of equally distributed points in logarithmic time scales.
Provided with a large distribution of points (where stands for either
or ), dense in some areas and sparser in others, a new time array is defined
to uniformly distribute the points along the curve:
139/237
Chapter 5 – Critical Analysis and Discussion
For a given set of shift factors and temperatures , the regression base is
described as:
For each element in (i.e. for each pair of constants , ), the value of the
right part of the WLF function is calculated for all temperatures , yielding
. The sum of square differences for all elements was then arranged
in .
For each iteration , the pair yielding minimum value of is found and
associated to a specific couple of constants and in .
minimum
28
The shape of the matrix arbitrary: any different space can have been adopted, without
jeopardizing the discussed procedure but possibly affecting its accuracy. In the following,
will be assumed, a choice motivated by the fact that there is no a-priori insight on the best fitting
values of any of the two constants.
140/237
Chapter 5 – Critical Analysis and Discussion
Curves in Figure 5.22 for the final shift-factor distributions show that the functions
are not all convergent to the same value. This ought to be expected, for the inherent
differences among the six regression bases : distributions which are more prone
to be fitted by a curve belonging to the WLF family will yield lower final values of
. Prior to the acquisition of final values for and two checks were
performed: all error curves must be asymptotic to a final value and the WLF
equations must be graphically fitting the trend of each regression base.
Due to the nature of the algorithm, notably due to the progressively increasing
density of the inquired mesh of points in . If and are the best possible
values of and respectively, the accuracy of the results at the end of the
iteration ( ) is:
141/237
Chapter 5 – Critical Analysis and Discussion
142/237
Chapter 5 – Critical Analysis and Discussion
143/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.24 –final mastercurve and elements of the series with a) 4-terms or b) 8-terms
144/237
Chapter 5 – Critical Analysis and Discussion
Gaps are meaningful to the DOI method, to the extent that the higher the number
of terms in the series, the smaller are the domains among which the constants
and are computed. If the domains are small enough and encompass one or
several wide gaps, the number of points in that specific domain may drop
dramatically, therefore severely affecting the accuracy of the regression in that
portion of the curve. For each interpolation, a limit to the minimum width of each
domain has been imposed as follow:
145/237
Chapter 5 – Critical Analysis and Discussion
For internal consistency among final results produced in this work, it was decided
to provide Prony series with equal number of terms for all moduli of all materials.
The number was then chosen to successfully fit all the curves with appropriate
accuracy and also avoid inappropriate “jiggling” (a technical term used in [243]
and visible in Figure 5.24) around the mean value.
For increasing number of terms, the accuracy was first estimated as the percental
error between experimental data29 and the generated Prony series, along the entire
time span of each mastercurve. For increasing number of terms, curves are shown
in Figure 5.25 and percental errors in Figure 5.26.
Error curves show that the accuracy initially increases if the number of terms also
increases: the maximum error between 5- and 15-terms series diminishes from
about to mere (cfr. Figure 5.26.a, b & c); nonetheless, data also
suggest that the error does not monotonically decrease with increasing number of
terms (cfr. Figure 5.26.d e).
To evaluate the accuracy of the Prony series fit to experimental data, a slightly
refined criterion may be used to knowingly choose the final number of terms for
the series. The normalized average of square root of square differences Γ , can be
defined as:
Γ Γ
being the number of elements in each data set and the value of the
modulus at time , calculated using the appropriate Prony series.
Figure 5.27 - error evaluation for Prony series with increasing number of terms
Figure 5.27 shows that a very high number of terms is not suited, if the purpose is
to maximise the accuracy of the fit.
All things considered, the choice of the conclusive number of elements for the
series has been chosen mediating among the three different criteria:
29The data as they are collected prior to the polish procedure, which has been used
to the sole purpose of running the fitting DOI algorithm and generate the series.
146/237
Chapter 5 – Critical Analysis and Discussion
This criterion was implemented in the general algorithm shown in Figure 5.13.
With the final mastercurves coefficients, the normalized sum of square differences
between the experimental data of shifted branches and interpolating Prony
mastercurves Γ was found, allowing to find and compare the overall measure Γ of
shift-factors accuracy among various distributions.
The shape of the 20-terms Prony series functions is displayed in Figure 5.28,
compared to experimental data and coefficients are shown in Table 5.8.
147/237
SG DG41 PVB
i
0 0,845 1,415 0,245 0,325 0,083 0,098
, time in seconds, moduli in MPa
1 21,90 3,50E+03 24,46 4,08E+03 2,570 5,06E+00 4,027 9,35E-05 6,609 1,47E-01 3,553 4,81E-04
2 11,52 2,22E+04 13,11 2,83E+04 2,719 2,12E+01 2,287 8,32E-04 2,929 6,01E-01 2,505 5,44E-03
3 13,56 1,15E+05 11,85 1,76E+05 5,437 6,60E+01 2,112 8,29E-03 2,439 2,31E+00 2,741 3,94E-02
4 10,98 5,13E+05 9,119 1,18E+06 5,282 2,11E+02 2,703 7,01E-02 2,119 8,21E+00 2,301 2,78E-01
Chapter 5 – Critical Analysis and Discussion
5 11,13 2,94E+06 8,660 4,73E+06 5,686 7,57E+02 3,148 6,91E-01 1,385 3,10E+01 1,379 1,85E+00
6 9,531 1,27E+07 8,004 5,90E+07 7,320 2,50E+03 3,282 8,30E+00 1,101 1,22E+02 0,720 1,87E+01
7 6,650 5,24E+07 8,698 3,53E+08 8,562 8,66E+03 3,533 7,88E+01 0,807 4,29E+02 0,565 1,55E+02
8 2,488 2,89E+08 4,398 2,78E+09 5,780 2,56E+04 5,220 6,50E+02 0,279 1,43E+03 0,452 1,15E+03
148/237
9 0,881 1,58E+09 2,791 1,47E+10 3,701 8,88E+04 5,224 5,04E+03 0,118 6,20E+03 0,154 7,56E+03
10 0,687 7,53E+09 1,378 1,10E+11 2,104 2,79E+05 5,506 5,77E+04 0,104 2,44E+04 0,097 8,09E+04
Table 5.8 - mastercurves coefficients,
11 0,523 4,40E+10 0,811 7,44E+11 0,445 9,25E+05 5,544 5,15E+05 0,093 9,11E+04 0,089 6,71E+05
12 0,495 1,66E+11 0,570 7,10E+12 0,154 3,38E+06 3,521 4,04E+06 0,082 3,40E+05 0,083 5,07E+06
13 0,234 8,99E+11 0,657 3,72E+13 0,088 1,17E+07 2,040 4,00E+07 0,057 1,29E+06 0,054 4,15E+07
14 0,356 4,66E+12 0,454 2,32E+14 0,056 3,91E+07 0,754 2,46E+08 0,053 4,96E+06 0,054 3,43E+08
15 0,315 2,20E+13 0,454 1,67E+15 0,083 1,31E+08 0,117 3,02E+09 0,055 1,79E+07 0,048 2,63E+09
16 0,247 1,12E+14 0,397 1,10E+16 0,045 4,54E+08 0,095 3,68E+10 0,035 6,90E+07 0,038 2,07E+10
17 0,219 5,55E+14 0,338 8,52E+16 0,044 1,69E+09 0,074 3,40E+11 0,024 2,88E+08 0,030 1,57E+11
18 0,203 2,75E+15 0,338 5,56E+17 0,043 5,39E+09 0,093 4,70E+12 0,033 1,22E+09 0,023 1,44E+12
19 0,157 1,30E+16 0,224 4,00E+18 0,038 1,74E+10 0,041 3,77E+13 0,021 4,28E+09 0,025 1,08E+13
20 0,091 5,92E+16 0,189 2,14E+19 0,029 5,94E+10 0,053 2,63E+14 0,003 1,45E+10 0,011 5,02E+13
Chapter 5 – Critical Analysis and Discussion
30
Tabulated times have been chosen according to the design instructions for load durations
described in prEN16612-2013 and used in CNR-DT 210/2013.
149/237
Chapter 5 – Critical Analysis and Discussion
§3.1.3.5) [114,291], application to direct shear tests as those performed in this work
must account for the actual testing conditions.
In chapter 4, it was discussed how experimental curves for relaxation and creep
(cfr. Figure 4.35 to 4.37) start to be reliable after a finite time interval. While short-
term effects for high and intermediate temperatures can be drawn through the
mastercurve gathering from experimental results at lower temperatures, the same
cannot be said for the lowest temperatures (i.e. for SG, for DG41 and
for PVB). Lowest temperatures must account for their own reliability, so short-term
properties cannot not be extrapolated from the series prior the reliability limit (this
justifies blank values in the top-left corner of Tables in Appendix G).
where and are the maximum and minimum reliability times for the
mastercurve at the reference temperature (Figure 5.30 and Table 5.9) and ,
are the WLF coefficients.
31
Coefficients and in Table 5.8 stand for this infinite modulus value (cfr. §3.1.3.1), for
relaxation and creep functions respectively.
150/237
Chapter 5 – Critical Analysis and Discussion
SG DG41 PVB
relax. creep relax. creep relax. creep
Note how the shape of the reliability limits is closely correlated to the shape of the
WLF equation; in fact, the same technique used to define the reliability limits can
be applied to plot the results on simple graphs, correlating time, temperature and
materials moduli. These graphical results are also shown in Appendix G.
An important aspect to evaluate concerns short-term properties of tested materials,
namely their known load/strain rate dependent properties (cfr. §3.2 and §4.2).
Analyses performed with LG specimens are dedicated to study viscous effects,
which are increasingly important with passing of time. Hence, the goal of
proposed viscoelastic models is to reliably replicate the effects of long-term
constant actions, while the rate at which these actions are applied has a limited
influence on the test outcome. Conversely, analyses described in §5.1.1 for the
tensile response of thermoplastic interlayers showed that the stiffness of such
materials for short-term effects is highly non-linear and influenced by the rate at
which actions are applied.
Instantaneous elastic properties are always superimposed to the viscoelastic effects
[114,282,291,292], therefore outcomes of viscoelastic analyses should always be
looked at with a critical eye. Results given in terms of Prony series mastercurves
and WLF coefficients are not meant to describe short-term properties of the
materials (indicatively, not much shorter than 5-10 minutes). Therefore, short-term
relaxation and creep moduli found with a linear viscoelastic approach should be
considered as conservative lower bound values of both creep and relaxation, as if
actions were applied at an infinitesimal stress or strain rate. In reality, the stiffness
of the materials for short-term effects is influenced by the strain-rate hardening
response described with the GR model in 5.1.1.2. In fact, the actual short-term
response of LG elements can possibly be better anticipated using the results of
short-term tensile tests on interlayer specimens, neglecting all viscous effects.
5.1.3 TEMPERATURE-DEPENTENT ADHESIVE PROPERTIES
Adhesive properties of interlayers in glass lamination are usually investigated with
TCT tests 32 : the ends of two-layered LG specimens, whose glass has been
artificially cracked on both sides, are pulled with a constant displacement rate. The
peak load measured during these tests can be correlated to the adhesion among the
materials. Nonetheless, this interpretation is intrinsically correlated to the
geometry of specimens and is not prone to be easily generalized [125,293]. While
adhesion is represented in many ways, what is measured is really the “practical
adhesion” which is defined as the force or the work required to detach a film or
coating from its substrate [294]. K. L. Mittal wrote [295]:
“A frequently asked question is: What is the best method for measuring adhesion?
The answer is: the best, rather most appropriate, method is one that simulates
usage stress conditions.”
32
Through Crack Tensile test
151/237
Chapter 5 – Critical Analysis and Discussion
“Using the TCT test, we have shown that steady state delamination can be
obtained only in a limited range of temperature and applied velocity.”
Alas, for specimens where delamination was measured at the end of tests, there is
no reliable data indicating the exact time frame when delamination process has
started, and how it has evolved over time. Nonetheless, in some tests the climatic
chamber was open for very brief periods after the end of the loading phase (i.e.
when presses are loaded by fastening of the loading bolt Figure 4.11 p.73) to
perform small operations and to check the setup. Observations during those
operations showed that some delamination had already occurred shortly after the
end of the loading phase (about 10 minutes to one hour). While the extent of the
delamination could not be measured accurately in those phases, neither compared
to the final delamination described in Appendix C, the simple observation that
some had occurred allows to ascribe to the loading phase itself a role in setting the
conditions for the surpass of the fracture toughness threshold. With this
consideration, an estimate of adhesion thoughtless can be given in terms of
supplied energy throughout the loading phase, similarly to more traditional tensile
tests [125]. An analysis in that direction has been performed for all specimens with
numerical integration of recorded force and displacement data. Figure 5.31
presents a graphical interpretation of this analysis for one of the specimens 33.
33
To help correlating these results to other parts of this manuscript, data is shown for specimen
S12 (one of the SG-laminates), one among those that had been previously chosen as a case study
for the data analysis in §4.2.2 and respective subparagraphs.
152/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.31 - a) displacement and b) load over time for specimen S12 (55°C), c)
instantaneous power and d) supplied energy by numerical integration
Figure 5.31.a shows the progression of the loading phase with the four imposed
displacement steps performed by design, to progressively load the specimens;
Figure 5.31.c highlights that the energy is supplied in the same number of
“chunks” corresponding to the manual bolts tightening.
Figure 5.32 - average delamination correlated with average loading energy density
Results show that most of the energy is supplied in the loading phase. The energy
density of the glass-interlayer interfaces can be computed dividing the supplied
energy for the area of the interlayers34. Details are in Appendix H.
34
approximately 2000 mm2, varying from one specimen to another for manufacturing tolerances.
153/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.32 condenses average data from the delamination analysis (Appendix C)
and outcomes of load energy density calculations (Appendix H); the resulting
correlation helps to draw some general conclusion on the fracture toughness of the
glass-interlayer interface. The maximum imposed shear stress in the loading
phase to the average extent of delamination for each material.
Table 5.10 – glass-interlayer fracture toughness estimate for varying temperatures
Temp. fracture toughness
PVB DG41 SG
0 >165 n/d n/d
10 >264 >100 n/d
14 n/d >100 n/d
17 n/d >97 n/d
20 >639 >106 >55
23 n/d >118 n/d
30 170 196 >788 >126
35 n/d n/d >148
40 <274 <518 >185
45 n/d n/d <115
50 <199 <333 160 190
55 n/d n/d <215
60 n/d n/d <161
This calculation requires once again the hypothesis that, glass being order of magnitudes stiffer
than interlayers, it can be considered relatively undeformable (cfr. §4.2.1). This hypothesis implies
that interlayer shear stress and strain can be approximated with uniform distributions.
154/237
Chapter 5 – Critical Analysis and Discussion
In PDLG specimens, a fictitious equivalent elastic modulus for fractured glass plies
can be estimated under two hypotheses: first, that fractured plies confined by
adherent interlayers can still be looked at as an isotropic material, able to bear
tensile and compressive stresses; the second, that the moment of inertia along the
strong axis of damaged and undamaged plies is the same. The validity the first
hypothesis is confirmed by experimental outcomes on the position of the neutral
axis in PDLG, at least when quasi-static imposed loads are modest (cfr. §4.3.3.2,
Table 4.9 p.110). For a simply supported LG beam having span , assembled with
three glass plies of thickness and inertia about the strongest axis,
with a concentrated load at mid-span, the equivalent modulus for one ( )
or two ( ) shattered plies is:
Table 5.11 - fictitious equivalent Young modulus [GPa] of damaged glass plies
unit A1 A2 B1 B2 B3 C1 C2 C3
stage I
† 66.3 67.8 66.4 65.8 68.5 66.7 67.1 69.1
stage II
29.9 28.8 18.2 16.3 22.3 16.3 22.9 14.3
stiffness increase % 22.5 21.2 13.7 12.4 16.3 12.2 17.1 10.3
stage III
16.3 16.9 16.7 15.3 13.6 9.09 10.7 11.0
stiffness increase % 49.2 50.0 50.3 46.5 39.7 27.3 31.9 31.8
† cfr. Table 4.8
where is a correction factor accounting for torsional effects. With two broken
plies, PDLG beams are bent laterally because of the asymmetric glass failure. As a
result, applying a vertical force at midspan generates non-negligible torsional
effects along with flexural ones. Throughout tests III-5 (Table 4.7), torsional
restraints35 (cfr. Figure 4.42 p98) were used to keep specimens in vertical position.
Accounting for the exact position of the restraints and the extent of the lateral sag,
the correction factor has been calibrated with a 3D FE linear elastic model.
5.2.2 DYNAMIC TESTS
With the same simplified approach used for the analysis of three-point bending
tests, the lower bound of expected frequencies for specimens in PDLG conditions
can be anticipated, by neglecting the mechanical contribution of broken glass plies
and interlayer while accounting for both of their masses 36 . An analytical
35
which have been installed as close as possible to midspan for those tests.
36
The exact weight of glass and interlayers has been determined experimentally with density tests
performed after of all the analyses described in §4.3. Density results are: glass
, SG and DG41 . The density of PC could
not be assessed experimentally, the literature value of was used [350] (although, for
this specific case study the volume of PC is so small that it could have been assimilated to DG41
interlayer with very limited consequences).
155/237
Chapter 5 – Critical Analysis and Discussion
formulation is known for the frequency of the first modal form of a uniform beam
of span and modulus , simply supported at both ends with uniform weight per
unit length (or mass per unit length ) with 37 [107,296]:
The expected natural frequency of specimens can be calculated in Table 5.12, using
experimental values of the glass elastic modulus from quasi-static three-point
bending tests (Table 4.8), along with material properties36 and geometries
described in §4.3.2. Such values have been compared to experimental data from
bending tests. An analysis of the influence of boundary conditions is performed in
Appendix F.
Table 5.12 - expected frequencies based on ULG three-point bending tests stiffness
†
L freq
A1 66,3 74,3
1,71 6,48 0,55 2,44 9,68 0,164 2,35
A2 67,8 75,1
B1 66,4 86,9
B2 65,8 5,35 8,40 0,85 2,44 1,08 0,212 2,7 86,5
B3 68,5 88,3
C1 66,7 85,1
C2 67,1 5,35 8,40 1,83* 2,44 1,08 0,225 2,7 85,3
C3 69,1 86,6
†
cfr. Table 4.8.
* an equivalent area of homogeneous DG41 interlayer is used, to account for the different
density of PC embedded within the interlayer matrix.
The approach used to compute the expected natural frequency in the ULG state
can be tweaked, using additional hypotheses, and applied for both PDLG states
(Stages II&III).
The minimum frequency expected, shown in Figure 5.33, can be computed by
neglecting the mechanical contribution of damaged plies and interlayers (i.e.
neglecting the TS effect), while accounting for all masses.
A better estimate of the experimental PDLG natural frequencies can be produced,
based on data from three-point bending tests elaborated in §5.2.1, which quantify
the TS effect. With an equivalent Young modulus for damaged plies (Table 5.11),
an equivalent modulus can be computed for the entire cross section by
homogenising material properties according to proportions shown below. Using
instead of in the previous formulation for yields results of the Equivalent
Stiffness Model (ESM) for stages II & III shown in Figure 5.33.
37
Validation of this analytical formulation, with calibration of the correction factor for this case
study is provided in Appendix F.
156/237
Chapter 5 – Critical Analysis and Discussion
157/237
Chapter 5 – Critical Analysis and Discussion
158/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.36 - flexibility method model for an expanding central fractured glass ply
38
The FE model was built using an equivalent non-isotropic temperature gradient: a numerical
gradient only applied in the two main dimensions of the specimen, but not across the thickness.
39
Also known as method of consistent deformations.
40
cfr. strain gauges for specimens Type B & C shown in Figure 4.45 p.104.
159/237
Chapter 5 – Critical Analysis and Discussion
41
Conventionally, the strain measured by instrument after thirty seconds from the failure of
the central ply was used. Figure 4.53 p.115 shows the progression of .
42
With the aforementioned hypotheses, such procedure is only rigorous for Stage II, a symmetrical
damage state. The same procedure is an approximation for the asymmetric Stage III, the curvature-
induced effects as with two broken plies are neglected.
160/237
Chapter 5 – Critical Analysis and Discussion
failure of the class ply. This should come to no surprise, as the ETEM with whom
the numerical analysis was carried was calibrated with experimental readings at
after failure of the glass ply.
Figure 5.37 - Stage II linear elastic finite element results using ETEM
161/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.39 – Stage III lateral deflection comparison (view from above)
Results from the numerical model highlight that a decrease of the interlayer’s
stiffness 43 , does not translate in a significant variation of the lateral deflection,
which is mainly affected by geometrical properties (i.e. the numerical model
showed that the thickness of the plies and the length of the specimens are the main
factors affecting the lateral deflection). Nonetheless, while the numerical model
shows a progressively diminishing sag for diminishing interlayer stiffness,
experimental measurements show instead a slight sag increase. This inconsistency
may be explained by pointing out that the ETEM model doesn’t account for the
delayed increase in the number of cracks, which in fact contributes to expanding
the lateral and central broken plies, thus inducing a further sag increase.
5.2.4 FULLY DAMAGED LG MODELS
Uniaxial tests performed on SG-laminated samples (cfr. §4.3.4 p.113), cut out of
Type A specimens after all tests had been performed, allow for a better
understanding of both TS effect and compressive response of glass fragments.
Experimental results described in Figure 4.65 can be compared with simple models
known in literature [38,55,63] and experimental outcomes on interlayer specimens
analysed in §4.1, whose results can be seen in Figure 4.4 p.66. The latter
comparison is shown in Figure 5.40. While interlayer results are shown for a strain
rate which is almost double the strain rate for FDLG44 (i.e. results of dotted curves
43
Either for conspicuous viscous relaxation or a choice of a different, softer, interlayer material.
44
The strain rate of has been calculated knowing the speed rate (
) and the span of dog-bone specimens between the clamps ( , cfr. §4.1.1).
162/237
Chapter 5 – Critical Analysis and Discussion
are overestimating the response one would record for the same specimens at
), a very significant TS effect is measured, due to the presence of
adherent glass fragments. While tests on FDLG could not be carried out up to very
high strain to preserve the instruments (cfr. §4.3.4), plausible forecasts can be
hypothesized, if glass fragments adhesion is eventually lost and FDLG ultimately
deteriorates to the interlayer only.
163/237
Chapter 5 – Critical Analysis and Discussion
164/237
Chapter 5 – Critical Analysis and Discussion
after the failure of the last glass ply. The decrease in stiffness and strength can be
explained by considering that glass cracks allowed for a widespread exposition of
the interlayer to moisture during the storage interval, triggering a decay in the
chemical bonds between glass and interlayers around each glass fragment [45].
A simple FDLG analytical beam model was defined using stress-strain curves
drawn from tensile and compression tests. For tensile stresses, “slow tests” were
used (cfr. Figure 4.65.a p.115), as the strain rate of three-point tests on FDLG beams
was the same imposed in those tests. The FDLG model curve in Figure 5.42.a was
obtained by computing the moment-curvature diagram with
experimental constitutive equations for tensile and compressive actions (Figure
5.42.b), then integrating the curvature twice over the beam axis. Boundary
conditions of a simply supported three-point bending test were considered.
Figure 5.42 shows how the analytical model can accurately predict the response of
five-month old FLDG beams (which ought to be expected, as the material
specimens used to perform tensile and compressive tests and draw the
experimental stress-strain correlations were extracted from those LG beams and
tested a few days after). On the other hand, the analytical model seems to
overestimate the ultimate load-bearing capacity and stiffness after the initial part of
the curves.
Differences between the model and 5-month experimental results can be explained
by comparing the experimentally recorded strain at the bottom and at the top of
beams (through “omega” transducers in Figure 4.44 p.99) with those of “material
specimens” in tension and compression tests (Figure 5.43).
Figure 5.43 – experimental compressive and tensile strain on top and bottom of the beam at
midspan, compared to FDLG analytical model results
In the tensile region (to the right in Figure 5.43), the model quite accurately
predicts the strain pattern in 5-month old FDLG beams. On the contrary, the model
is only accurate in the initial part of the compressed region (to the left in Figure
5.43), for loads approximately on third the ultimate capacity, while it vastly
underestimates the strain at top fibres for higher loads. As an example, at a load of
(85% of the ultimate capacity) the error is around in the tensile region,
but is about one order of magnitude higher in the compressed region.
Very high differences in the compressed region can be explained considering the
greater confinement of glass fragments in uniaxial compression tests, rather than in
the compressed region on top of the FDLG beams in three-point tests. In the latter
case, a progressive expulsion of shards and fragments was observed from top of
the beams, which can explain both the lower stiffness and the reduced value of the
ultimate load capacity with respect to the model.
165/237
Chapter 5 – Critical Analysis and Discussion
Figure 5.45 shows results of FE analyses, where it can be observed that mean TS
effect may be approximated using an exponential interpolation.
Figure 5.45 - increment in SG stiffness vs. delamination ratio from FE elastic analyses
For each delamination ratio, the figure shows the upper and lower values of the
stiffness increase of the interlayer due to TS, corresponding to the two cracking
patterns of Figure 5.44. For instance, the stiffness of the interlayer is about twice as
big for debonded surface ratios of ~40% and four times as big for ratios of ~15%. It
can be observed that at increasing debonded surface ratio the two limit values of
the stiffness increment approach one another and for high delamination ratios
(more than 80%) they are almost zero.
166/237
Chapter 5 – Critical Analysis and Discussion
167/237
Chapter 5 – Critical Analysis and Discussion
168/237
Chapter 6 – Applications
6 APPLICATIONS
To highlight and clarify some of the results discussed within this manuscript,
some examples are provided, showing how experimental findings can be used
to improve LG design. The following studies will try to focus on the main
aspects of novelty, whereas standard and well-established procedures may be
referenced and summarized but not dealt with extensively.
In the following paragraphs, two case-studies are shown: the first, dealing with a
cold-bending application of large glazing panes, the second with the post-failure
safety assessment of a compressed LG strut.
6.1 COLD-BENT GLAZING
A recent concept design for a museum in Parma, Italy is taken as a case study. A
ramp which contains part of an historical automotive exhibition grants access from
the ground floor to the 1st floor of the building. This part of the structure is
designed around two semi-circles: one for the external profile and one for an
internal court which serves as a light well (Figure 6.1).
The original design proposed the use of hot-bent LG panes (cfr. §1.1.1) to cover
both the internal and external glazing surfaces. The hot-bending-process being
very expensive and detrimental to the optical properties of glass, the aim of the
following analysis is to propose solutions for alternative cold-bent LG solution.
The geometry of both glazing elements is presented in Figure 6.2, Figure 6.3 and
Table 6.1.
169/237
Chapter 6 – Applications
Figure 6.2 - layout of the case study with details on LG plates sizes and curvatures
Figure 6.3 - geometry of LG panes for the internal and external glazing
Table 6.1 - median axis geometry for the cold-bent laminated elements
external glazing internal glazing
height mm mm
arc length mm mm
vertical supports distance mm mm
maximum sag mm mm
radius mm mm
arc rad rad
1
The additional boundary effects on vertical edges and the bi-dimensional structure of the LG
plates are neglected for these analyses, both assumptions going towards safety.
170/237
Chapter 6 – Applications
installation must be performed before the end of June and the opening to the
public must not be allowed earlier than September or it must be checked that
over hours above elapse between installation and use 3,4,
the building is a category C3 according to Eurocode 1, BSI 6399-1 and CNR-
DT210/2013, so that the horizontal design loads are: a uniformly distributed
and a linear load at a height ,
distributed loads of wind, including safety factors and dynamic effects, account
for a total of according to Italian standards [298,299],
the self-weight only induces a minor in-plane compression, whose effects will
be neglected for a safer design (slightly diminishing the tensile stress in glass).
The design maximum tensile stress in the glass plies is [38,300]:
with fixed , , , , ,
, and ; other coefficients may vary for different
load durations and are shown in Table 6.2 along with the final result for .
Table 6.2 – coefficients and values for the design tensile stress of tempered glass plies
service (50 yrs) bending (30’)* crowd (30”)* wind gust (5”)
0.26† 0.65† 0.78† 0.88†
0.80 0.80 1.00 1.00
0.89‡ 0.89‡ 0.89‡ 0.89‡
52 MPa 58 MPa 79 MPa 82 MPa
* a nominal duration of 30” is used for the crowd actions on balustrades in [38]
†
results from Linear Elastic Fracture Model (LEFM) are preferred to prEN16612
‡
values are different for the internal and external glazing; with a precautionary approach, the
lowest value is considered
2
This prescription as a design choice, to work with known values of the relaxation modulus.
3
This prescription in terms of months is intended to guarantee a minimum amount of time above
30°C between installation and opening to public. To this end, a comparison of weather station
LIMP (Parma Airport, N44°49.33' / E10°17.72', elevation 49.1 m MSL) data over the years 2015,
2016 and 2017 was performed. Daily logs for 2017 are shown:
4
For external glazing elements that are exposed to direct sunlight, this hypothesis is highly
conservative, considering that surface temperatures are likely to rise above ambient temperature.
171/237
Chapter 6 – Applications
The design values of the interlayer shear modulus have been calculated as:
Figure 6.4 - a) scheme of cold-bent LG and b) an extracted vertical stripe “beam”, used to
compute the effects of horizontal loads accounting for c) cold-bending pre-existing stresses
172/237
Chapter 6 – Applications
For the Poisson effect, cold-bending generates compressive and tensile stresses in
perpendicular directions ( and in Figure 6.4c). Such stresses diminish over
time as an effect of the decaying shear modulus and coupling abilities of all
interlayers.
For load cases 2, 3 and 4, cold bending effects are accounted for as pre-stress:
tensile effects induced by cold-bending add up with those produced by the live
loads. On the other hand, to further increase safety, compressive stresses induced
by cold-bending will not be accounted for.
Figure 6.5 - schemes for studying the effects of live-loads in Cases 2, 3 & 4
Cases 2, 3 and 4 contain both actions of cold-bending and live loads on glass, two
effects which refer to different times. In accordance with the CNR standard [38], a
procedure analogous to the Palmer-Miner linear damage hypothesis for fatigue is
adopted [301,302]. Therefore, verifications are satisfied when:
where and are the -th imposed and resistant stresses respectively.
The Enhanced Effective Thickness (EET) method is used to calculate the equivalent
thickness of the LG sandwich structure [113] (cfr. §2.2.2).
An example of the procedure for cases 1 and 2 of the external LG plate is described,
all other cases and the internal glazing can be dealt with in an analogous manner.
Case 1
The LG plate is made with two glass plies of the same thickness ,
bonded with a thick SG interlayer. For cold-bending, a soft
interlayer goes towards safety, hence and .
Considering a stripe of LG of unitary amplitude and the glass elastic modulus
, the coupling coefficient can be calculated as:
where:
Such value of very close to 1, hints that the behaviour of the plate is very close to
the monolithic limit. The equivalent thickness which can be used to model the
bending stiffness of the plate is therefore:
173/237
Chapter 6 – Applications
where:
The equivalent thickness for calculating the stresses in the -th glass ply is:
The horizontal stress and vertical stress for the -th ply in Case 1 are:
Yielding the new –slightly lower- values for cold-bending induced stresses :
The procedure for calculating the stresses induced by the horizontal live loads of
the crowd are also computed similarly. The coupling coefficient for the extracted
vertical stripe in Figure 6.4b is calculated with the appropriate value of :
and are calculated with the same formulae, while the maximum bending
moment, which is reached at a height (cfr. Figure 6.4b) is:
174/237
Chapter 6 – Applications
The verification must now account for the pre-existing stresses due to the cold-
bending process, hence the Palmer-Miner criterion must be satisfied for the four
stresses , , and :
Multiplications per zero are the result of the procedure of neglecting the beneficial
effects of compressive stresses in the glass plies5. Altogether, verifications for Case
2 are successfully checked with a maximum 62% use of the design resistance of
glass (for the outer ply in the vertical direction).
Other cases
Checks for the other two cases and for all other interlayer materials are not
reported in full, as the procedures are alike those shown above. Nonetheless, while
tensile effects of cold-bending and live loads were not addictive in case 2 for the
external glazing, it must be pointed out that those can in fact induce addictive
effects, depending on the direction of the bending and live loads.
6.1.3 RESULTS, DESIGN OPTIONS AND INTERLAYERS COMPARISON
Table 6.4 shows the results of verifications were carried out for several
combinations of glass plies and interlayer thicknesses. Numbers under the
verification columns represent the percentage of use of the tensile glass
mechanical resources according to the Palmer-Miner superposition criterion, which
was generally observed either for case 2 (0 is no stress, 100 is failure). Results are
found to be dependent on the mechanical properties and thickness of the interlayer
materials. The following statements can be made:
5
With the exception for those induced by the tempering process.
175/237
Chapter 6 – Applications
when dealing with combinations of cold-bending effects and live loads, stiffer
interlayers are not always beneficial (e.g. SG-laminated glazing with 15+15 mm
plies cannot be used, while softer interlayers allow for multiple solutions)
varying the interlayer thickness may or may not be beneficial (cfr. Table 6.4,
10+8 mm combinations using DG41 and PVB for external glazing: for DG41
stiffer is better, for PVB thinner is better).
Table 6.4 – analysis of glass/interlayer thickness (verifications are satisfied for )
V [%] V [%]
[mm] [mm] SG DG41 PVB SG DG41 PVB
3,04 97 97 fail 3,04 fail fail
10+10 1,52 fail fail fail 6+6 1,52 fail fail fail
0,76 fail fail fail 0,76 fail fail fail
3,04 88 86 fail 3,04 68 77 fail
10+12 1,52 93 94 fail 6+8 1,52 80 fail
0,76 98 98 fail 0,76 88 91 fail
3,04 88 84 fail 3,04 68 83 fail
12+10 1,52 93 fail 8+6 1,52 81 90 fail
0,76 98 98 fail 0,76 88 94 fail
3,04 92 76 99 3,04 54 66 fail
12+12 1,52 90 82 8+8 1,52 63 70 fail
0,76 89 86 94 0,76 68 72 96
3,04 98 71 89 3,04 49 52 80
12+15 1,52 95 78 84 8+10 1,52 51 56
0,76 94 84 81 0,76 55 57 72
3,04 98 76 78 3,04 49 56 96
15+12 1,52 96 82 77 10+8 1,52 51 58 90
0,76 95 86 80 0,76 55 59 82
3,04 fail 77 74 3,04 53 46 76
15+15 1,52 fail 84 74 10+10 1,52 50 47 71
0,76 fail 90 81 0,76 48 48 65
With such a vast array of design options, the final design can rely on different
criteria such as the economic edge or availability of interlayers; the most valuable
options have been circled in Table 6.4 6. Final options are then:
for internal glazing
12+10/1.52DG41 minimizes overall thickness and glass cost,
12+12/1.52 PVB minimizes interlayer cost,
6
Even if not always beneficial from a mechanical or economical standpoint, 1.52 mm interlayers
are most reliable, matching the thicknesses tested within the experimental campaigns described in
this work, and because scaling effects regarding interlayer thickness were not investigated.
Besides, this choice is made for most structural LG applications and has proven to be reliable.
176/237
Chapter 6 – Applications
7
This structure of the bridge was studied in 2013 with Prof. S. G. Morano at University of
Florence. The first analysis of the LG struts, which serves as reference for the following analyses,
was performed in my master thesis, supervised by Prof. M. Orlando and Ph.D. E. Cagnacci [110].
177/237
Chapter 6 – Applications
LG struts
arches
(hollow steel)
tie-rods
(steel cables)
transversal frames
(steel plates)
Figure 6.7 - portion of the deck showing the wind bracing made with LG struts and tie-rods
The object of this second case-study is the design of LG struts forming part of the
wind bracing trusses in the deck of the bridge. The design will account for the
event of accidental failure of any of the glass pies included in those elements.
The bridge is designed to be built in five phases, schematically illustrated in Figure
6.8. LG struts first, then LG plates for the causeway, are both designed to be
installed in the final part of the assembly. Therefore, the only actions on those
elements needing verification are referred to SLS and ULS of the finished structure.
Figure 6.8 –phase 1) groundwork and banks 2) assembly of main steel elements 3)
positioning by rotation 4) completion of the steel structure 5) LG elements installation
178/237
Chapter 6 – Applications
The design maximum tensile stress in the glass plies is calculated as previously
shown in §6.1, yielding the results shown in Table 6.6.
with , , , , ,
, and . Remaining coefficients may vary and are
reported in Table 6.2, along with the result for .
Table 6.6 –design tensile resistance of the glass plies in MPa used for the LG struts
glass 30 seconds 12 hours 50 years
tempered 81,7 74,5 68,4
float annealed 20,0 12,8 6,68
Figure 6.9 - view general scheme of the strut a) from above and b) in a lateral view
Figure 6.10 - a) old 3-plies design with PVB and b) new 5-plies design using DG41
179/237
Chapter 6 – Applications
Figure 6.11 - a) scheme of the strut and conventional axes b) partition used for validation
The model is validated by studying the case of failure of the two external plies
(Figure 6.10b). For this validation, analytical calculations are carried out under the
hypothesis of perfect coupling among the glass plies ( ). This assumption
implies that , which has been shown to be here a reasonable hypothesis: a
parametric analysis for a multi-layered LG of similar total thickness (composed by
five 7-mm thick glass plies and four 1.52-mm thick interlayers) was performed
with the EET method. Results show that a relative error between layered and
monolithic models of 5% for (Figure 6.12).
The action in Figure 6.13a, which annuls the expansion of the broken glass plies,
can be calculated using the equivalent thermal expansion model (cfr. §5.2.3). With
the average values of for type B specimens and corresponding average
for stage III (cfr. Table 5.13), one finds:
Figure 6.13 - external action equilibrating the expansion b) effects on the undamaged plies
180/237
Chapter 6 – Applications
The numerical model is assembled using the same properties for the glass plies
and an elastic modulus of , corresponding to with
. Results are shown in Figure 6.14, yielding:
Figure 6.15 – stresses after 50 yrs in half beam model in the event of failure of ply 5
181/237
Chapter 6 – Applications
Table 6.7 and Figure 6.16 show the results at midspan 8 for stresses in the
undamaged plies, because of the failure of one or both external plies.
Table 6.7 - maximum and minimum stresses at midspan in [MPa] and lateral sag
case ply 5 fails ply 4 fails plies 4 & 5 fail
time 30” 12h 50y 30” 12h 50y 30” 12h 50y
min -3.69 -1.88 -0.75 -1.16 -0.75 -1.07 -8.05 -4.63 -2.39
1
max -2.22 -0.26 +0.34 -0.41 +0.14 +0.47 -4.32 -0.12 +1.02
min -1.81 -1.44 -0.94 -0.26 +0.47 +0.58 -3.44 -1.75 -0.52
2
max -0.34 +0.17 +0.15 +0.49 +1.38 +2.13 +0.27 +2.74 +2.90
min +0.04 -1.79 -0.54 +0.60 -0.41 -1.77 +1.06 -3.84 -4.19
3
max +4.42 +3.03 +1.72 +2.86 +2.27 +2.83 +12.1 +9.56 +6.00
min +4.81 +5.86 +3.69
4
max +6.28 +7.49 +4.78
min +4.18 +4.67 +2.58
5
max +4.94 +5.57 +4.13
sag [mm] 0.817 0.824 0.533 0.420 0.650 1.27 2.18 2.55 2.35
Figure 6.16 – graphical representation of the evolution of stresses over time in the
intact glass plies (dotted lines are placed to signpost the position of damaged plies)
8
It has been observed that different geometries and mechanical properties of the anchorages can
produce high concentrated stresses (e.g. Figure 6.15, left part of the model). While these may be
significant, the analysis of local effects is deeply connected to the specific anchoring technology,
the study of which goes beyond the purpose of both this work and this specific example.
182/237
Chapter 6 – Applications
Figure 6.17 - scheme for maximum stresses verification in a loaded strut after partial failure
The resulting bending moment at midspan has been calculated with a step
by step analysis, whose results are shown in the figure below. Increments of the
compressive action are applied on a deformed configuration, which transitions
from a configuration to an equilibrium configuration . For the final FE
analysis the final midspan moment is then applied at both anchorages
as shown in Figure 6.17, which is a choice on the side of safety.
Figure 6.18 - result of the step-by-step non-linear analysis on the compressed strut
Despite the -small- bending moment, compressive actions are beneficial in this case
from the point of view of tensile stresses in the intact plies. On the other side, the
tensile force , induces stresses so small that the final effects are
almost indistinguishable from those found after the fractured glass expansion
alone. In light of these considerations, all verifications are met for loaded elements
in post-failure conditions.
For a bent strut such as the damaged LG element, a compressive action on doesn’t
lead to the risk of loss of stability, as the element is already in a deformed
configuration . Therefore, the compressive action will progressively
magnify the existing sag with no sudden changes of conformation. Nevertheless,
as the sag produced by failure of the external plies is indeed very small (cfr. Table
6.7), a verification with respect of the Euler’s critical load was performed
accounting for the three intact plies in the worst post-failure case scenario. The
equivalent stiffness was calculated using both the EET method and the FE model
using a safety coefficient of two for the interlayer modulus at 30”, with similar
results for both models: and . The
critical load is therefore , much higher than any design load
envisaged for these elements.
183/237
Chapter 6 – Applications
This analysis concluded that effects of failure of any of the lateral tempered glass
plies are not a threat to the safety for both ULS and SLS. Provided that failure only
occurs on one side of the LG strut, replacement of damaged elements can be
scheduled with no urgency from the point of view of structural safety.
184/237
Chapter 7 – Thesis and Developments
7 THESIS
AND DEVELOPMENTS
185/237
Chapter 7 – Thesis and Developments
especially for fully-damaged elements; in fact, the results of tensile tests were
used later in the experimental analyses1. While this part of the research may
seem less connected with practical LG applications, the study of tensile
properties and models is relevant as key draft standards propose a subdivision
of interlayers in “stiffness families”, based uniquely on the results of similar
tests (prEN 16613:2013 [303] and FprEN 16613:2016 [303]).
To some extent, results based on short-term tensile properties could be used for
LG design if loads have a duration comparable to the duration of the tests.
However, studies on viscous properties of interlayers allow to compare short-
term to long-term effects, suggesting that it would be unadvisable to extend
tensile test results to extrapolate the long-time response. In fact, such choice
would involve overlooking some of the most relevant mechanics governing the
interlayers response, possibly overestimating their properties.
2. A simple method for defining the pseudo-yielding limit of thermoplastics is
discussed, a definition providing an alternative to the one proposed in
standards such as “ISO 527:2012 Plastics - Determination of tensile properties” [264]
and “ASTM D638-14 Standard Test Method for Tensile Properties of Plastics” [66].
The proposed criterion allows to define a point on the engineering stress-strain
curves where the properties change abruptly, and can be applied even to those
materials which do not show a peak in the measured stress with a subsequent
softening branch (§5.1.1.3 p.128).
Converging on a definition for the pseudo-yielding criterion would allow to
identify limits in terms of stress and strain under which tensile results on
interlayers can be used in LG design using LE models.
Viscoelastic, long-term and adhesive properties of interlayers
3. Most applications and designs of LG involve actions and effects which take
place over a large amount of time. Thermoplastic materials are known to build
up conspicuous viscous effects which are in turn greatly dependent on the
operating temperatures. Tensile tests on interlayers are unfit to study the long-
term response of such materials within the layered structure of LG. Therefore,
the long-term viscoelastic properties have been studied with small specimens of
LG using three of the most widespread interlayers within the LG industry: PVB,
SG and DG41. The analyses ultimately provided analytical results in the form of
Prony series mastercurves and time-temperature superposition functions,
allowing to predict the effects of viscoelastic flow within well-defined reliability
limits. Results were provided in terms of both stress relaxation and creep
compliance moduli as a function of time and temperature. This double
formulation is needed for the design of structural LG elements: while these
quantities are often confused one another, they represent in fact two different
physical effects. As rigorous procedures for converting relaxation to creep
functions (and vice-versa) involve the introduction of Laplace transforms,
which are arguably not suited for practical engineering, separate formulations
based on empirical measurements have been deemed to be the most reliable
and effective way of proceeding.
The direct the method used within this work to derive viscoelastic properties of
interlayers 2 can be compared with a few other researches which used dynamic
1
To study the tension stiffening effect in fully-damaged LG specimens (cfr. §5.2.4 p.181).
2
A direct approach to this problem was called for by D. Delincé in his work [7], and by a handful
186/237
Chapter 7 – Thesis and Developments
techniques instead. While new results roughly agree with the –limited-
previously available data, the proposed approach is believed to be more reliable
and consistent with operating conditions of interlayers in LG elements (cfr. §0
p.148 and §5.1.2.6 p.149).
Experimental results gathered within the analyses of viscoelastic properties of
interlayers have been used in a case study for a cold-bending LG design
application, allowing for a comparison of the advantages and disadvantages of
different materials (cfr. §6.1 p.169).
Some researchers suggested that a high glass transition temperature for
interlayers is be desirable, to take advantage of a full range of applications in
the architectural industry [69]: “a higher glass transition temperature is suitable for
more demanding architectural applications”. While this belief has been observed to
be quite common throughout the pertinent literature, this research has shown
that one can hardly rely on this parameter alone to anticipate the effectiveness
of the coupling abilities of interlayers. For instance, the glass transition
temperature of DG41 is much closer to SG rather than it is to traditional PVB;
nevertheless, for medium- to long-time applied actions, the response of DG41 is
much closer to PVB than it is to SG for a wide range of temperatures.
4. The loss of adhesion between glass and interlayers, also known as delamination,
is a factor of major concern for LG, affecting durability, mechanical properties
and aesthetics. Within this work, a correlation between the strain energy and
the glass-interlayer interface fracture toughness is provided, based on
phenomenological observations. The correlation allows to assess the risk of
delamination between glass and different interlayer materials at varying
temperatures (cfr. §5.1.3 p.151).
Additional analyses on adhesive properties have been carried out for fully
damaged LG specimens, to provide a correlation between the loss of adhesion
which comes with the tension stiffening effect using numerical models.
Post-failure properties of LG elements
5. The response of three-ply, progressively damaged LG beams was investigated
with static and dynamic in-plane experimental analyses. Results showed that
the choice of interlayers does not have a significant influence for the in-plane
response in the pre-failure phase. On the other hand, the post-failure response
depends on the coupling ability of interlayers and glass fragments. In partial
damage conditions, the beams stiffness did not reduce to the stiffness of intact
glass plies alone; likewise, for fully damaged elements, the glass stiffness did
not vanish completely, but the small tempered glass fragments of were able to
effectively contribute to the stiffens and residual load-bearing capacity. This
behaviour was ascribed to the coupling of glass fragments and interlayers, a
phenomenon called tension stiffening effect, which has been observed to suffer a
reduction with the progression of the debonding. Experimental results
highlighted that SG and DG41 provide a similarly remarkable contribution to
the overall stiffness in partial damage conditions, while for fully damaged
conditions the ionoplast SG interlayer is vastly superior. Static and dynamic
tests have been compared in terms of in-plane bending stiffness, allowing to
validate an equivalent stiffness model for partially damaged LG beams.
of other researchers which also have tackled this problem in the past; such suggestion has been
firmly embraced and followed through in this work.
187/237
Chapter 7 – Thesis and Developments
188/237
Chapter 7 – Thesis and Developments
While this vision still seems to be out of reach, new researches should always go in
the direction of bridging the gap between different fields of inquiry. With a unified
model, simplified equations could also be brought up to quantify one specific
aspect of the polymer response that is of interest for one specific application, while
knowingly disregarding others.
In this work, a few wide-ranging conclusions can be given in the direction of a
unified model, with some advices for future developments.
Decomposition of intermolecular and network effects
Aside from the different hyperelastic models, two approaches have been follower
in material science for plastics and rubbers, to replicate the short-term correlations
between stress and strain: one is a statistical approach 3, in which elastic constraints
are linked to the distribution of interactions between polymer chains (i.e. crosslinks
density and chains length), the second is an empirical one, where experimental
data is used as a regression base for fitting functions that are exclusively based on
experience and intuition [304,305].
Arguments can be made to support the adoption of a mixed model, which
empirically accounts for the statistical distribution of cross-links within polymer
molecules. In that direction, the GR model provides an option for dividing in two
parts the response of thermoplastic materials, which are only supposedly ascribed
to intermolecular forces and network entanglement of molecules. While our results
fit perfectly within a the framework of such theory, there is currently no clear
evidence on whether or not such speculations are an accurate representation of
these micromechanical interactions. As such, it is useful to remind that our
calibrations are not yet meant to divide the materials response in that sense, but
only aim to a faithful reconstruction of the intrinsic response.
If it is proven that the reversible part of deformation can be ascribed to the network
effects and the plastic ones to the rupture of intermolecular cross-links, a
calibration of the decomposition could be achieved using cyclic load/hold/unload
tests at various speed rates and temperatures, using a pattern similar the one
shown in Figure 7.1.
3
Also referred to as the “Gaussian” or “molecular” approach.
4
Applicable Standards are ASTM D638 [66] and ISO 527 [264].
189/237
Chapter 7 – Thesis and Developments
Figure 7.2 - uniaxial tests on PVB for a) at various displacement rates (0.037 to 5.000
mm/s) and b) 0.33 mm/s for temperatures between -15 and +22°C (Kott and Vogel 2003)
These observations hints that the coefficients of the newly introduced GR model
may be in turn conveniently described as functions of strain rate, temperature and
possibly other parameters, to consistently replicate the outcome of tensile tests
with the same mathematical formulation.
5
At the time, the commercial acronym was SGP, rather than today’s SG (-P for “Plus”).
190/237
Chapter 7 – Thesis and Developments
very long amount of time. Nevertheless, predictions should always keep in mind
that properties of the materials are likely to change over the course of many years.
Therefore, predictions for years should be looked at with a critical eye,
acknowledging that changes in the chemical compositions of the materials could
lead to inaccurate predictions of the mechanical response.
191/237
Chapter 7 – Thesis and Developments
Step 2 find mastercurve coefficients using the DOI method, then compute the
sum of square differences for each for each of the distributions.
Step 3 use the WLF regression procedure on the random distributions of shift
factors and store final sum of least squares for each one.
Step 5 Produce small random variations of the saved distributions to fill the
population of the next generation with new shift factors distributions and
repeat from step 2 until the last generation is complete.
statistical ambition, this single result hints that the contribution of glass fragments
to the ultimate strength can be almost as significant as the contribution to stiffness.
6
Calculated using the characteristic strength for glass .
192/237
Chapter 7 – Thesis and Developments
Figure 7.4 - different scenarios for out-of plane bending of partially damaged LG:
a) top ply fractured and b) bottom ply fractured
193/237
Chapter 7 – Thesis and Developments
and solutions have been proposed. Yet, in many architectural applications, glazing
used for the envelope of buildings are made with insulated glass 7 rather than
single blocks of LG. In turn, the inner and outer panes comprising the insulated
element can be made with LG.
To my knowledge, cold-bent insulated glazing has not yet been studied nor tested.
The extension of results brought in this work to such applications call for an
analysis of the mechanical properties of edge seal spacers used for insulated glass
units. A first analysis in that direction has been performed by A. Hagl [308], but the
focus of that research was put to relative rotations of the edges, rather than the
shear displacements that can be induced by cold-bending (Figure 7.5).
Provided that the mechanical properties of the sealing unit can be correctly
understood, the mutual effects among LG panes can be modelled and added to the
cold-bending effects on the single panes.
7
Typicallt two glass panes separated by a vacuum or gas filled space, to reduce heat transfer.
194/237
Appendices
APPENDIX A
A theoretical approach is proposed for modelling the readings of the load
recording device during the opening of the climatic chamber. The problem is
related to thermal deformations of parts composing the device: dynamometric
ring, aluminium pads and LVDT transducer, which is itself composed by several
elements (external case, shaft, inner core, etc.). Parts are presented in Figure A.1.
where is the heat loss function, and are the mass of the object and the
specific heat capacity of the material and is the varying temperature over time .
In the case Newton's law for convective heat transfer applies, one can assume:
The loss heat function is defined as the rate of transferred heat . On the right
side of the equation, is the heat transfer coefficient, the external surface of the
object generating the convective flow and the ambient temperature. With the
applied convention the heat flow is positive when directed from the object to the
surrounding fluid (in our situation, an infinite volume of freely convective air).
Solving the differential equation with assigned boundary condition
and rearranging for yields:
195/237
Appendices
By subtracting this effect to the previous, one finds the cumulative effect for an
object subject to and temperature cycle between and between and :
To calibrate a model for the readings of the load recording devices in the effective
test situation, we can consider a conditioned test which has reached a steady
temperature of (i.e. the temperature of all the parts of the device, of
the chamber internal wall surface and the air contained inside the chamber itself
have all reached the set-point temperature for tests, therefore all parts are in
thermal equilibrium). The temperature of the air surrounding the metal parts at the
instant of the chamber opening is considered to be immediately dropping from
to 1 . The temperature at the closing of the chamber is
1
This is not an unreasonable hypothesis, given that the hot air flow is free to leave the volume of
the climatic chamber and spread into the surrounding room. The surrounding ambient being
extremely large, it will be considered uninfluenced by the mixture of the hot air coming from the
open chamber, in such a way that its temperature can be considered constant at all times.
2 This hypothesis may seem less rigorous than the previous, as the chamber needs time to build up
the climatic condition defined by the set-point. Nonetheless, readings of the in-built sensors of the
chamber suggest that the time needed to bring back the temperature to the set-point for a
previously conditioned chamber is very limited, in the range of to seconds, depending
on the temperature difference between the set-point and the ambient temperature.
196/237
Appendices
Results show that the thermal inertia of the ring is higher compared to the
transducer. Hence, the temperature variation of the transducer is greater, but this
elemental is also faster in coming back to the temperature and the effects of
the chamber opening are negligible approximately 10 minutes after it is closed.
Based on the separate temperature variations of the parts, we investigated whether
load readings were affected with positive or negative correlations between
temperature variations and load measurements. Starting from the dynamometric
ring itself, the correlation is clearly negative: if the temperature of the ring
diminishes, the structure contracts simulating a compression load, yielding a
positive value. However, this point of view alone may be misleading, as the
dynamometric ring and the transducer share one fixed point (on the upper side of
the ring). From the reference system originated from the shared fixed point: as
temperature lowers, the ring contracts, but because the transducer cools faster it
contracts even more. From the point of view of the measuring device this
phenomenon is perceived as an actual expansion, rather than a contraction (the
same logic applies in reverse if the elements are heating rather than cooling).
Therefore, it is safe to say that the load readings of the system ring-transducer
are affected by a negative correlation. This may be presented in the form of:
197/237
Appendices
The last step of the calibration procedure is finding the correlation among
coefficients and using the experimental data. Readings of the “virtual” load
registered upon chamber opening may be written as:
Figure A.3 - validation of the thermodynamic model for the load-recording devices
198/237
Appendices
APPENDIX B
A few criteria have been evaluated to model the response of testing devices used
for the campaign on long-term viscoelastic properties of interlayer materials (cfr.
§4.2). The problem is defined in terms of finding suitable mathematical models for
reproducing the response of the base contrast elements shown in Figure 4.22 in
terms of load over displacement curves. The final result of the analyses described
in the following in shown in Figure 4.23.
Many attempts for modelling the response of the setup did not immediately yield
satisfactory results, either for being poorly able to represent the shape of
experimental results shown in or for creating numerical glitches. Among others,
the following were evaluated:
δ δ
Where and δ are the loads and displacements, while , , , and are five
coefficients to be found to minimize the least-squares problem for assigned :
δ δ
This can be also be expressed as the root finding problem of the partial derivatives
of the function with respects of the coefficients to :
δ δ
Presented in this form, the problem is partially non-linear and cannot be solved
straightforwardly. To work around this issue, a slightly different approach has
been undertaken.
Data in the first part of the experimental curves was more suited for a polynomial
fit, while an exponential model is more appropriate in the final part. A polynomial
interpolation was performed on the initial part alone from the entire data set. After
several attempts, the boundaries for this - arbitrarily fixed - reduced data set have
199/237
Appendices
The polynomial fit is intended as the solution of the linear problem issued by:
Yielding the values of the variables , and and thus the function:
δ δ
With this part of the interpolation, values calculated with the function
have been subtracted from the original dataset, producing a new dataset
which was used to interpolate the exponential part of the curves:
The is presented in the following figures by the gray dot clouds. The
remaining part of the procedure was to interpolate the new set of data with the
exponential function using a - now possible - least square interpolation:
The problem can be made linear by recognizing that one can take the natural
logarithm of both sides of the equation, so that various methods of polynomial fit
approaches can be used:
Taking and :
Solving this system yields the values of the constants and . Values of the
constants have been calculated for the three devices used.
The final curves for modeling the response of the setup are shown in Figure B.2 for
the three devices.
200/237
Appendices
Figure B.1 - calibration of the nonlinear elastic correlation for device 048
Figure B.2 - calibration of the nonlinear elastic correlation all testing devices
201/237
Appendices
APPENDIX C
While specimens delamination was visible to the naked eye, photographing it
properly proved to be a very challenging task. Several techniques were put to the
test varying lighting, backgrounds and angles. Taking pictures from a square angle
proved to be an impractical solution, as there seems to be no way to avoid having a
disturbing camera reflection on the final image. A solution was found with the
setup in Figure C.1, using a dark blue cloth background, which helps the
delamination reflections to stand out, and a white paper sheet to intercept the
reflection from the glass surface.
b)
202/237
Appendices
203/237
Appendices
40°C 50°C
D05 – 1.28%
204/237
Appendices
S22 – 0,22%
205/237
Appendices
APPENDIX D
206/237
Appendices
APPENDIX E
207/237
Appendices
APPENDIX F
Hypothesis for the analytical formulation provided by W. Young and R. Budynas
[296], whose results were used for the modal analysis in §5.2.1, are not an exact
match of experimental conditions. Namely, supports at the ends of the specimen
are not located along the axis, but are supporting the bottom edge instead.
Moreover, a thin aluminium pad is attached do the glass plies, to serve as an
interface with the steel rollers and protect the specimen from the risk of local stress
concentrations. Moreover, the proposed analytical formulation neglects the
presence of – small – lateral protrusions of the specimens beyond the supports at
both ends.
Two LE numerical models were built: one to replicate the conditions under the
theoretical hypotheses and one that corresponds to the setup configurations. The
models were identical in every single aspect with the exception of the insertion of
the aluminium pads (the central mass was calibrated to replicate the mass of the
accelerometers).
A comparison of the results concluded that the differences between the two models
response in terms of natural frequencies are very limited.
Type A B C
calibrated using specimen data A1 B1 C1
Model 1 frequency [Hz] 75.03 88.01 86.15
Model 2 frequency [Hz] 74.97 87.97 86.08
correction factor 0.9992 0.9995 0.9991
208/237
Appendices
APPENDIX G
Part1 : tabulated results
Table G.1 - SG relaxation modulus G
Temp°C 5s 30s 10m 11h 1w 3mo. 6mo. 50yrs.
20 83,1* 56,5 35,9 19,2 15,3 4,56
25 75,6 41,8 23,6 9,90 7,82 3,71
†
30 83,0 57,6 27,3 11,5 5,51 4,89 2,99
35 77,3† 61,5† 38,6 12,6 5,62 3,98 3,71 2,42
40 53,5† 39,6† 20,0 5,48 3,86 2,98 2,78 1,86
45 29,3† 17,8† 6,53 3,63 2,76 2,25 2,11 1,34
†
50 8,29 5,35† 3,71 2,49 1,97 1,53 1,43 0,92
† †
55 3,71 3,11 2,39 1,60 1,23 0,96 0,91 0,86*
† †
60 2,18 1,82 1,36 0,92
____________________________
†
Numbers shown for and should always be looked at as lower conservative values, as
actual moduli are possibly significantly higher (cfr. §5.1.2.6)
* Higher and lower moduli that should be extracted from the mastercurve (cfr. §5.1.2.6).
209/237
Appendices
____________________________
†
Numbers shown for and should always be looked at as lower conservative values, as
actual moduli are possibly significantly higher (cfr. §5.1.2.6)
* Higher and lower moduli that should be extracted from the mastercurve (cfr. §5.1.2.6).
210/237
Appendices
211/237
Appendices
212/237
Appendices
213/237
Appendices
APPENDIX H
Average energy density for PVB, DG41 and SG specimens (continues to next page)
mat. temp. spec. delaminated area* [%] energy density average
PVB P03 - 165
0°C P04 - 154 154
P05 - 141
P14 - 264
10°C 251
P15 - 239
P16 - 639
20°C 634
P17 - 630
P11 - 170
30°C P12 1,18 196 197
P13 2,97 226
P08 1,35 175
40°C P09 1,32 243 231
P10 0,10 274
P06 - 251
50°C P07 0,54 199 226
P02 0,28 229
DG41 D13 - 66
10°C D14 - 84 83
D15 - 100
D19 - 100
14°C D20 - 78 81
D22 - 66
D23 - 97
17°C D25 - 88 85
D26 - 70
D16 - 76
20°C D17 - 95 92
D18 - 106
D10 - 109
23°C D11 - 118 110
D12 - 102
D06 - 696
30°C D07 - 760 748
D08 - 788
D01 3,75 612
40°C 565
D02 7,29 518
D03 1,52 510
50°C D04 1,73 333 408
D05 1,28 382
SG S27 - 48
20°C S28 - 55 52
S29 - 53
S24 0,02 125
30°C 126
S25 - 126
S17 - 118
35°C S18 - 130 132
S19 - 148
S23 - 185
40°C 181
S26 - 178
S13 0,97 115
45°C 146
S15 1,12 178
S20 2,50 227
50°C S21 - 160 192
S22 0,22 190
S09 7,65 260
55°C 237
S12 0,19 215
S14 6,06 161
60°C 196
S16 8,64 232
*cfr. Appendix C
214/237
Appendices
APPENDIX I
Step-by-step DOI method for building a Prony series, given experimental data set
, monotonous, decreasing, with an horizontal asymptote for .
Step 1 Step 8
Define the interval within which the Find the characteristic time of the -th
single exponential terms operate in group with linear interpolation among
terms of (the greater this value, the two consecutive log10-time values
the fewer elements are in the series). and which encompass the value of
Step 2 :
Step 7
For each group, normalize the variation
as:
215/237
Appendices
216/237
INDEX OF FIGURES
Chapter I – Structural glass introduction
Figure 1.1 - a) North Rose of Notre Dame (Paris), 1270 […] 1
Figure 1.2 - a) & b) Cité des Sciences et de l'Industrie greenhouses […] 2
Figure 1.3 - “The Snow House” concept by Santambrogio & Arosio […] 3
Figure 1.4 - worldwide interest in time for the query “laminated glass” 4
Figure 1.5 - stress distribution in chemically or thermal tempered glass 8
Figure 1.6 - typical failure patterns for a) annealed, b) heat-strenghtened […] 9
Figure 1.7 - tempered glass fragments expelled […] 9
Figure 1.8 - scheme of a) entanglement and b) intermolecular cross-links 12
Figure 1.9 - thermo-mechanical curves for a) elastic modulus [67] […] 13
Figure 1.10 - stress decomposition by Haward and Thackray [76] 14
Figure 1.11 - PVB-glass bonding mechanism (Keller and Mortelmans 1999 [79]) 17
Chapter II - State-of-the-art and aim of the research
Figure 2.1 - specimens used in the tensile tests (Santarsiero, 2017 [128]) 25
Figure 2.2 - scheme of the setup for shear tests by C Vallabhan et al. [135] 26
Figure 2.3 - schematic representation of the interlayer studies […] 27
Figure 2.4 - intermolecular and network interaction scenarios, Vegt 2002 [59] 28
Figure 2.5 - stress decomposition, Klompen 2005 [156] 29
Figure 2.6 - scheme of four-point bending tests […] 31
Figure 2.7 - 3D FE model by C. Bedon, 2017 [184] […] 33
Figure 2.8 - a) & b) embedded hollow/plain steel c) steel or GFRP […] 34
Figure 2.9 - the “ZAPPI1” beam specimen after a 4-points bend testing [201] 34
Figure 2.10 - five-points bending tests for SG-laminated beams […] 35
Figure 2.11 - cross section of the GFRP-reinforced FL […] 36
Figure 2.12 - results for four-point bending tests on P-FG (left) […] 36
Figure 2.13 - a) cross-section and lateral view of LG reinforced […] 37
Figure 2.14 - method for bonding the pre-tensioned tendons […] 37
Figure 2.15 - stress distribution in pre-stressed LG beams (Louter, 2014 [186]) 38
Figure 2.16 - desired load-carrying behaviour for LG beams […] 38
Figure 2.17 - qualitative stress distribution in cold-bent LG at short (left) […] 38
Figure 2.18 - Roof over the IJsei bus station in Amsterdam […] 39
Figure 2.19 - a) homogenized model problem by Galuppi-Royer […] 41
Figure 2.20 - a) interlayer specimens for tensile tests […] 42
Chapter III - Mechanics
Figure 3.1 - typical linear elastic and hyperelastic response 45
Figure 3.2 - hyperelastic models used within this manuscript […] 46
Figure 3.3 - a) function for varying , and b) hyperelastic NH models 48
Figure 3.4 - hyperelastic Mooney-Rivlin models for varying over ratio 48
Figure 3.5 - hyperelastic Yeoh model dependency from coefficients and 49
Figure 3.6 - conceptual visualization of a) creep and b) stress relaxation 51
Figure 3.7 - spring and dashpot elements 52
Figure 3.8 - rheological representation of a) the Maxwell element […] 52
Figure 3.9 – LG specimen and testing in a torsional rheometer […] 57
Figure 3.10 - relaxation & storage moduli plotted against and axes […] 57
Figure 3.11 - one-dimensional constitutive models for viscoplasticity [242,258] 59
Figure 3.12 - geometry of common adhesive joints for LG […] 61
Figure 3.13 - a) shear stress distribution in SG for varying temperatures […] 62
Figure 3.14 - tri-linear adhesion diagram, 2010 Cottone et al. [262] 62
217/237
Chapter IV - Experimental tests
Figure 4.1 - a) cutting of the specimens and b) manual finishing 64
Figure 4.2 - Type IV dog-bone specimen according to ASTM D638-14 64
Figure 4.3 - progression of tensile tests on PVB 66
Figure 4.4 - stress-strain curves: a) SG, b) DG41, c) PVB and d) averages 66
Figure 4.5 - construction for a) chord- and b) regression-slope tensile modulus 67
Figure 4.6 - mean engineering stress-strain curves […] 68
Figure 4.7 - experimental observations on non-homogeneous deformations 70
Figure 4.8 - setup for time-lapse photogrammetry 71
Figure 4.9 - measurement lines used for photogrammetric analysis 71
Figure 4.10 - SG tensile test at 500 mm/min time lapse, all shots 72
Figure 4.11 - testing frames 73
Figure 4.12 - testing devices inside the climatic chamber 74
Figure 4.13 - design two-phases test output for stress and strain […] 75
Figure 4.14 - scheme of the elements in contact with LG specimens 78
Figure 4.15 - four phases (clockwise) of installation of angular shims 78
Figure 4.16 - guiding dowels a) on a tested specimen and b) before a test 79
Figure 4.17 - thermal effects recorded by load-monitoring devices 80
Figure 4.18 - typical position and evolution of cracks […] 80
Figure 4.19 - displacement monitoring for the adjustment cycles 81
Figure 4.20 - a) setup for adjustment of deformations in lead foils […] 82
Figure 4.21 - setup for calibration of rings and base contrast elements 82
Figure 4.22 - a) load and b) displacement curves for base elements calibration 83
Figure 4.23 - non-linear deformations of the base contrast elements 83
Figure 4.24 - Gantt diagram of performed tests 85
Figure 4.25 - recorded loads trough dynamometric rings 86
Figure 4.26 - recorded displacements at the bottom of the central ply 86
Figure 4.27 - load measured by the reference dynamometric ring […] 87
Figure 4.28 - linear drift superimposed to load readings 87
Figure 4.29 - load readings corrected for linear drift 88
Figure 4.30 - displacement readings corrected for initial measurements 88
Figure 4.31 - displacement readings corrected for base deformation 89
Figure 4.32 - correction for thermal effects on dynamometric rings 89
Figure 4.33 - a) stress and b) strain in loading and viscoelastic recovery phases 90
Figure 4.34 - relaxation and creep compliance moduli for the case-study 91
Figure 4.35 - a) relaxation and b) inverse compliance moduli for SG interlayer† 92
Figure 4.36 - a) relaxation and b) inverse compliance moduli for DG41 92
Figure 4.37 - a) relaxation and b) inverse compliance moduli for PVB 92
Figure 4.38 - example of delamination analysis for specimen S09 93
Figure 4.39 - full-scale LG beams specimens type A, B and C [mm] 95
Figure 4.40 - damage progression imposed by design 96
Figure 4.41 - lateral setup scheme for a) type A and b) type B or C specimens 97
Figure 4.42 - picture of the setup arranged for specimen B1 98
Figure 4.43 –pads used for the setup interfaces (on stage II PDLG specimen B2) 98
Figure 4.44 - measure instruments installed on type A beams 99
Figure 4.45 - measure instruments installed on type B and C beams 100
Figure 4.46 - omega transducers installed for FDLG tests 101
Figure 4.47 - glass modulus evaluation over the 2nd load ramp […] 103
Figure 4.48 - normalized FFT response spectrum for undamaged specimen A2 104
Figure 4.49 - normalized FFT response spectrum for undamaged specimen B2 105
Figure 4.50 - a) measured strain for undamaged specimen B1 […] 105
218/237
Figure 4.51 - measured strain at failure of the central ply 106
Figure 4.52 - measured strain at failure of the first lateral ply 107
Figure 4.53 - short-term effects after central ply failure on specimens type B&C 107
Figure 4.54 – normalized spectral density for specimen B2 PDLG stages 108
Figure 4.55 - load-displacement curves for PDLG specimen B1 108
Figure 4.56 - a) snapshot of interpolated strain distribution […] 109
Figure 4.57 - neutral axis evolution for specimen C2 in test IV-1 110
Figure 4.58 - normalized moment over midspan sag for all specimens 111
Figure 4.59 - ageing effects for FDLG specimens for bending tests 111
Figure 4.60 - neutral axis interpolation over FDLG bending tests 112
Figure 4.61 - evolution of the depth of the neutral axis in all FDLG specimens 112
Figure 4.62 - long-term weathering effects on FDLG 113
Figure 4.63 - specimens extracted for tensile and compression tests […] 114
Figure 4.64 - setup for a) tensile and b) compression tests […] 114
Figure 4.65 - experimental curves of tensile tests on material specimens 115
Figure 4.66 - delamination observed from tensile tests photographic follow-up 115
Figure 4.67 - experimental curves of compressive tests on material specimens 116
Chapter V - Critical analysis and discussion
Figure 5.1 - photogrammetric strain analysis (SG at ) 119
Figure 5.2 - true stress-strain diagrams compared to engineering 119
Figure 5.3 - hyperelastic models for a) SG, b) DG41 and c) PVB 121
Figure 5.4 - initial response detail for very small strain 123
Figure 5.5 - intermolecular series model results 124
Figure 5.6 - obtained by subtraction for a) SG, b) DG41 and c) PVB 125
Figure 5.7 - scheme of the process used to produce the final models 126
Figure 5.8 - model true stress-strain results for a) SG b) DG41 and c) PVB 127
Figure 5.9 - model engineering stress-strain results for a) SG b) DG41 c) PVB 127
Figure 5.10 - a proposed criterion for defining the pseudo-yielding limit 130
Figure 5.11 - pseudo-yielding limit for the performed tests 130
Figure 5.12 - example of shifted experimental branches for SG creep modulus 132
Figure 5.13 - process used to define WLF and mastercurve coefficients 132
Figure 5.14 - shift factors for a) relaxation and b) creep 26 133
Figure 5.15 - shifted branches for a) relaxation, b) creep and c) averages 134
Figure 5.16 - arrangements for consecutive experimental branches 135
Figure 5.17 - procedure for overlapping areas 136
Figure 5.18 - scheme of the linking procedure27 137
Figure 5.19 - target points for linking branches 138
Figure 5.20 - example of a thirty-points re-sampling 139
Figure 5.21 - WLF functions for a) relaxation and b) creep mastercurves 142
Figure 5.22 - convergence of the WLF regression algorithm 142
Figure 5.23 - detail of data polish procedure outcome 143
Figure 5.24 - final mastercurve and elements of the series […] 144
Figure 5.25 - Prony fit comparison with different term series 144
Figure 5.26 - error for a) 5- , b) 10- , c) 15 -, d) 20- and e) 40-term series 145
Figure 5.27 - error evaluation for Prony series with increasing number of terms 146
Figure 5.28 - final 20-terms Prony series mastercurves 147
Figure 5.29 - DG41 relaxation modulus over time for 149
Figure 5.30 - reliability limits for mastercurves Prony series 150
Figure 5.31 - a) displacement and b) load over time for specimen S12 […] 153
Figure 5.32 - delamination correlated with average loading energy density 153
Figure 5.33 - experimental and model frequencies comparison 157
219/237
Figure 5.34 - fragmentation of broken tempered glass plies 158
Figure 5.35 - expansion of a confined tempered glass ply upon failure 159
Figure 5.36 - flexibility method for an expanding central fractured glass ply 159
Figure 5.37 - Stage II linear elastic finite element results using ETEM 161
Figure 5.38 - numerical model compared to experimental strain measures 161
Figure 5.39 - Stage III lateral deflection comparison (view from above) 162
Figure 5.40 - Experimental results of tensile tests on FDLG specimens 163
Figure 5.41 - FDLG compression tests and proposed constitutive law 163
Figure 5.42 - FDLG bending tests comparison to analytical model 164
Figure 5.43 - experimental compressive and tensile strain […] 165
Figure 5.44 - Schematic crack patterns: a) symmetric layout, b) maximum offset 166
Figure 5.45 - increment in SG stiffness vs. delamination ratio […] 166
Chapter VI - Applications
Figure 6.1 - external rendering and sketch of the semi-circular ramp 169
Figure 6.2 - layout of the case study with details on LG plates […] 170
Figure 6.3 - geometry of LG panes for the internal and external glazing 170
Figure 6.4 - a) scheme of cold-bent LG and b) an extracted vertical stripe […] 172
Figure 6.5 - schemes for studying the effects of live-loads in Cases 2, 3 & 4 173
Figure 6.6 - rendering of the pedestrian bridge 177
Figure 6.7 - portion of the deck showing the wind bracing […] 178
Figure 6.8 - phase 1) groundwork and banks 2) assembly […] 178
Figure 6.9 - view general scheme of the strut a) from above and b) lateral view 179
Figure 6.10 - a) old 3-plies design with PVB and b) new 5-plies design […] 179
Figure 6.11 - a) scheme of the strut and conventional axes […] 180
Figure 6.12 - external action equilibrating the expansion […] 180
Figure 6.13 - numerical model used for validation 181
Figure 6.14 – stresses after 50 yrs in half beam model […] 181
Figure 6.15 – graphical representation of the evolution of stresses […] 182
Figure 6.16 - scheme for maximum stresses verification in a loaded strut […] 183
Chapter VII - Thesis and developments
Figure 7.1 - cyclic tests for determination of plastic deformations 189
Figure 7.2 - uniaxial tests on PVB for a) at various displacement rates […] 190
Figure 7.3 - tensile tests on SG at room temperature (engineering) [307] 191
Figure 7.4 - scenarios for out-of plane bending of partially damaged LG […] 193
Figure 7.5 - a) insulated glass components, b) analyses […] 194
220/237
INDEX OF TABLES
Chapter I - Structural glass introduction
Table 1.1 - standard glass plies sizes and thicknesses in construction indstry 7
Table 1.2 - characteristic tensile strength for glass (min. required [38]) 8
Table 1.3 – list of the main properties of glass useful for LG design 9
Table 1.4 – glass transition and (melt) temperature for some interlayer materials 13
Chapter III - Mechanics
Table 3.1 - stress in uniaxial elongation for a Neo-Hookean material 47
Table 3.2 - uniaxial elongation stress for incompressible Mooney-Rivlin material 48
Table 3.3 - uniaxial elongation stress for incompressible Yeoh material 49
Table 3.4 - uniaxial elongation stress for incompressible Ogden material 50
Chapter IV- Experimental tests
Table 4.1 - number of specimens per material per test velocity 65
Table 4.2 – experimental minimum ultimate strain and engineering stress 67
Table 4.3 – mean chord- and (regression)-slope tensile modulus in MPa 69
Table 4.4 – SG tensile tests photographic time-lapse shooting details 71
Table 4.5 - list of performed tests and tested specimens 84
Table 4.6 - specimens type and nomenclature 95
Table 4.7 - tests performed on LG beam specimens 101
Table 4.8 – experimental elastic modulus of glass 104
Table 4.9 - neutral axis rise for bending tests 110
Table 4.10 - experimental natural frequencies [Hz] 113
Table 4.11 – tensile and compressive properties of FDLG 116
Chapter V - Critical analysis and discussion
Table 5.1 – 1h Poisson ratio for interlayers at 23°C/50%RH 118
Table 5.2 – accuracy of hyperelastic models for first and final parts of the curves 121
Table 5.3 - coefficients for hyperelastic Yeoh model 122
Table 5.4 - GR model coefficients for true stress-strain 126
Table 5.5 - pseudo-yielding limit for the performed tests 130
Table 5.6 - shift factors 133
Table 5.7 - WLF coefficients 141
Table 5.8 – mastercurves at , time in seconds, moduli in MPa 148
Table 5.9 – reliability times in seconds 151
Table 5.10 - glass-interlayer fracture toughness for varying temperatures 154
Table 5.11 - fictitious equivalent Young modulus [GPa] of damaged glass plies 155
Table 5.12 - expected frequencies based on ULG 3-point bending tests stiffness 156
Table 5.13 - equivalent modulus and gradient for PDLG models 160
Table 5.14 – effective modulus in FDLG with SG interlayer [MPa] 163
Table 5.15 – mechanical parameters for SG-FDLG made with tempered plies 164
Table 5.16 - noteworthy effects of progressive damage on LG beams 167
Chapter VI - Applications
Table 6.1 - median axis geometry for the cold-bent laminated elements 170
Table 6.2 - design tensile stress of tempered glass plies 171
Table 6.3 - interlayer relaxation modulus and tabulated results in [MPa] 172
Table 6.4 – analysis of glass/interlayer thickness […] 176
Table 6.5 - design loads for LG struts verification 179
Table 6.6 –design tensile resistance of the glass plies in MPa […] 179
Table 6.7 - maximum and minimum stresses at midspan and lateral sag 182
221/237
REFERENCES
[1] Barnes S. The design of caring environments and the quality of life of older people.
Ageing Soc 2002;22:775–89. doi:10.1017/S0144686X02008899.
[2] Sorensen S, Brunnstrom G. Quality of light and quality of life: An intervention study
among older people. Light Res Technol 1995;27:113–8.
doi:10.1177/14771535950270020501.
[3] Boyce PR, Veitch JA, Newsham GR, Jones CC, Heerwagen J, Myer M, et al. Lighting
quality and office work: two field simulation experiments. Light Res Technol
2006;38:191–223. doi:10.1191/1365782806lrt161oa.
[4] Krewinkel HW. Glass Buildings: Material, Structure, and Detail. Princeton
Architectural Press; 1998.
[5] Cagnacci E, Orlando M, Spinelli P. Il vetro come materiale strutturale. Polistampa
[Italian]; 2010.
[6] Ledbetter SR, Walker AR, Keiller AP. Structural use of glass. J Archit Eng
2006;12:137–49. doi:10.1061/(ASCE)1076-0431(2006)12:3(137).
[7] Delincé D. Experimental Approaches for Assessing Time and Temperature
Dependent Performances of Fractured Laminated Safety Glass. Gent University,
2014.
[8] Transparency Market Research. Construction Glass Market - Global Industry
Analysis, Size, Share, Trends and Forecast. 2016.
[9] Blascone M. Facciate in vetro strutturale: analisi tipologica e progettazione di una
facciata con pinne verticali in vetro. Univertità Degli Studi di Firenze, 2017.
[10] Carraher C. Polymer Chemistry. Polym Chem 2003:36–44. doi:10.1039/c4py01415d.
[11] Duser A Van, Jagota A, Bennison SJ. Analysis of Glass/Polyvinyl Butyral Laminates
Subjected to Uniform Pressure. J Eng Mech 1999;125:435–42.
doi:10.1061/(ASCE)0733-9399(1999)125:4(435).
[12] Preston DL, Tonks DL, Wallace DC. Model of plastic deformation for extreme
loading conditions. J Appl Phys 2003;93:211–20. doi:10.1063/1.1524706.
[13] Pelfrene J, Kuntsche J, Van Dam S, Van Paepegem W, Schneider J. Critical
assessment of the post-breakage performance of blast loaded laminated glazing:
Experiments and simulations. Int J Impact Eng 2016;88:61–71.
doi:10.1016/j.ijimpeng.2015.09.008.
[14] Larcher M, Arrigoni M, Bedon C, van Doormaal JCAM, Haberacker C, Hüsken G, et
al. Design of Blast-Loaded Glazing Windows and Facades: A Review of Essential
Requirements towards Standardization. Adv Civ Eng 2016;2016:1–14.
doi:10.1155/2016/2604232.
[15] Zhang X, Hao H, Wang Z. Experimental study of laminated glass window responses
under impulsive and blast loading. Int J Impact Eng 2015;78:1–19.
doi:10.1016/j.ijimpeng.2014.11.020.
[16] Zhang X, Hao H, Shi Y, Cui J. The mechanical properties of Polyvinyl Butyral (PVB)
at high strain rates. Constr Build Mater 2015;93:404–15.
doi:10.1016/j.conbuildmat.2015.04.057.
[17] Tilley RJD. Understanding Solids: The Science of Materials. John Wiley & Sons;
2013.
[18] Mihir J, Prabodh P, Shaan K. Flat Glass Industry At A Glance. 2015.
[19] Bansal NP, Doremus RH. Handbook of Glass Properties. Elsevier Science; 2013.
[20] Larsen AW, Merrild H, Christensen TH. Recycling of glass: accounting of
greenhouse gases and global warming contributions. Waste Manag Res 2009;27:754–
62. doi:10.1177/0734242X09342148.
[21] Ozer C, Boysan F, Cırık S. Technical Note, Recycling of Glass Waste - An
Application for the Esentepe Campus of Sakarya University. Procedia Earth Planet
Sci 2015;15:602–6. doi:10.1016/j.proeps.2015.08.111.
222/237
[22] Ward PD, Brownlee D. Rare Earth – why complex life is uncommon in the universe.
Int J Astrobiol 2003. doi:10.1017/S1473550404261924.
[23] Editor-in-Chief Haimei Zhang. Building Materials in Civil Engineering, 1st Edition.
Woodhead Publishing; 2011.
[24] Fanderlick I. Optical properties of glass. United States: Elsevier Science Pub.,New
York, NY; 1983.
[25] Vogel W, editor. Glass Chemistry. Berlin, Heidelberg: Springer Berlin Heidelberg;
1994. doi:10.1007/978-3-642-78723-2.
[26] Brodsley L, Frank C, Steeds JW. Prince Rupert’s Drops. Notes Rec R Soc 1986;41:1–
26. doi:10.1098/rsnr.1986.0001.
[27] Saint Gobain Glass. Manuale del Vetro 2006.
[28] BS EN 572-1:2012 Glass in building. Basic soda lime silicate glass products. Float
glass 2012.
[29] Pilkington NSG Group Flat Glass Bussines. Pilkington Glass Handbook 2010
2010:200.
[30] Varshneya AK. Fundamentals of Inorganic Glasses 1st Edition. New York:
Academic Press, Inc. An Imprint of Elsevier; 2013.
[31] Creyke WE, Sainsbury IE, Morrell R. Design with Non-Ductile Materials. Springer
Netherlands; 1982.
[32] Barsom JM. Fracture of Tempered Glass. J Am Ceram Soc 1968;51:75–8.
doi:10.1111/j.1151-2916.1968.tb11840.x.
[33] Varshneya AK. Chemical Strengthening of Glass: Lessons Learned and Yet To Be
Learned. Int J Appl Glas Sci 2010;1:131–42. doi:10.1111/j.2041-1294.2010.00010.x.
[34] Anderson TL. Fracture mechanics: fundamentals and applications. CRC press; 2017.
[35] Michalske TA, Bunker BC. The fracturing of glass. Sci Am 1987;257:122–9.
[36] Broek D. Elementary engineering fracture mechanics. Springer Science & Business
Media; 2012.
[37] Soules TF, Busbey RF, Rekhson SM, Markovsky A, Burke MA. Finite-Element
Calculation of Stresses in Glass Parts Undergoing Viscous Relaxation. J Am Ceram
Soc 1987;70:90–5. doi:10.1111/j.1151-2916.1987.tb04935.x.
[38] CNR-DT 210/2013. n.d.
[39] Overend M, De Gaetano S, Haldimann M. Diagnostic Interpretation of Glass Failure.
Struct Eng Int 2007;17:151–8. doi:10.2749/101686607780680790.
[40] Bennison SJ, Hayes RA, Pesek SC, Shaffer WA, Smith CA. Laminates comprising
ionomer interlayers with low haze and high moisture resistance 2013.
[41] Behr RA, Minor JE, Linden MP, Vallabhan CVG. Laminated Glass Units Under
Uniform Lateral Pressure. J Struct Eng 1985;111:1037–50. doi:10.1061/(ASCE)0733-
9445(1985)111:5(1037).
[42] Foraboschi P. Hybrid laminated-glass plate: Design and assessment. Compos Struct
2013;106:250–63. doi:10.1016/j.compstruct.2013.05.041.
[43] Hooper JA. On the bending of architectural laminated glass. Int J Mech Sci
1973;15:309–23. doi:10.1016/0020-7403(73)90012-X.
[44] Xu J, Li Y, Liu B, Zhu M, Ge D. Experimental study on mechanical behavior of PVB
laminated glass under quasi-static and dynamic loadings. Compos Part B Eng
2011;42:302–8. doi:10.1016/j.compositesb.2010.10.009.
[45] Louter C, Belis J, Veer F, Lebet J-P. Durability of SG-laminated reinforced glass
beams: Effects of temperature, thermal cycling, humidity and load-duration. Constr.
Build. Mater., vol. 27, 2012, p. 280–92. doi:10.1016/j.conbuildmat.2011.07.046.
[46] Stevanovic D, Kalyanasundaramb S, Lowe A, Jar P-Y. Mode I and mode II
delamination properties of glass/vinyl-ester composite toughened by particulate
modified interlayers. Compos Sci Technol 2003;63:1949–64. doi:10.1016/S0266-
3538(03)00162-3.
[47] Serafinavius T, Lebet JP, Louter C, Lenkimas T, Kuranovas A. Long-term laminated
223/237
glass four point bending test with PVB, EVA and SG interlayers at different
temperatures. Procedia Eng 2013;57:996–1004. doi:10.1016/j.proeng.2013.04.126.
[48] Teotia M, Soni RK. Polymer Interlayers for Glass Lamination-A Review. Int J Sci Res
2014;3:1264–70.
[49] Walley SM, Field JE, Blair PW, Milford AJ. The effect of temperature on the impact
behaviour of glass/polycarbonate laminates. Int J Impact Eng 2004;30:31–53.
doi:10.1016/S0734-743X(03)00046-0.
[50] Govaert LE, Timmermans PHM, Brekelmans WAM. The Influence of Intrinsic Strain
Softening on Strain Localization in Polycarbonate: Modeling and Experimental
Validation. J Eng Mater Technol 2000;122:177. doi:10.1115/1.482784.
[51] Tervoort TA, Smit RJM, Brekelmans WAM, Govaert LE. A Constitutive Equation for
the Elasto-Viscoplastic Deformation of Glassy Polymers. Mech Time-Dependent
Mater 1997;1:269–91. doi:10.1023/A:1009720708029.
[52] Royer-Carfagni G, Silvestri M. Fail-safe point fixing of structural glass. New
advances. Eng Struct 2009;31:1661–76. doi:10.1016/j.engstruct.2009.02.050.
[53] Biolzi L, Orlando M, Piscitelli LR, Spinelli P. Static and dynamic response of
progressively damaged ionoplast laminated glass beams. Compos Struct 2016;157.
doi:10.1016/j.compstruct.2016.09.004.
[54] Bennison SJ, Qin MHX, Davies PS. High-performance laminated glass for
structurally efficient glazing. Innov. Light. Struct. Sustain. Facades, Hong Kong:
2008, p. 1–12.
[55] Bennison SJ, Stelzer I, Davies PS, Sloan JG, Xiaokun H, Liu G. Calculation Methods
for the Structural Behavior of Laminated Glass. Glas Perform Days 2009 2009:433–4.
[56] Kennar GA. Interlayer for laminated safety glass 1977.
[57] Sharpe WN, editor. Springer Handbook of Experimental Solid Mechanics. Boston,
MA: Springer US; 2008. doi:10.1007/978-0-387-30877-7.
[58] Staudinger H. Über Polymerisation. Berichte Der Dtsch Chem Gesellschaft (A B Ser
1920;53:1073–85. doi:10.1002/cber.19200530627.
[59] Vegt a K Van Der, van der Vegt AK. From Polymers to Plastics. VSSD 2002:268.
doi:10.1017/CBO9781107415324.004.
[60] Fried JR. Polymer Science and Technology, 3rd Edition. Prentice Hall; 2014.
[61] Gilbert GK. The Moon’s Face: A Study of the Origin of Its Features. Philosophical
Society of Washington; 1893.
[62] Ivanov I V. Analysis, modelling, and optimization of laminated glasses as plane
beam. Int J Solids Struct 2006;43:6887–907. doi:10.1016/j.ijsolstr.2006.02.014.
[63] Bennison SJ. Structural Properties of Laminated Glass, Short Course. Glas. Process.
Days, Tampere (Fi): 2009.
[64] Galuppi L, Royer-Carfagni G. Laminated beams with viscoelastic interlayer. Int J
Solids Struct 2012;49:2637–45. doi:10.1016/j.ijsolstr.2012.05.028.
[65] Biolzi L, Cagnacci E, Orlando M, Piscitelli L, Rosati G. Long term response of glass-
PVB double-lap joints. Compos Part B Eng 2014;63:41–9.
doi:10.1016/j.compositesb.2014.03.016.
[66] International ASTM. ASTM D638-14, Standard Test Method for Tensile Properties of
Plastics. West Conshohocken, PA: 2014. doi:10.1520/D0638-14.
[67] Harper CA. Handbook of Plastics, Elastomers, and Composites, Fourth Edition.
New York, NY: McGraw-Hill; 2002.
[68] Askadskiĭ AA. Computational materials science of polymers. Cambridge Int Science
Publishing; 2003.
[69] Karagiannis A, Szlosek ML, Farhoudi Y. Multilayer interlayers having high Tg and
high modulus 2016.
[70] Meitzner R, Schulze SH. Method for determination of parameters for moisture
simulations in photovoltaic modules and laminated glass. Sol Energy Mater Sol
Cells 2016;144:23–8. doi:10.1016/j.solmat.2015.08.014.
224/237
[71] Dhaliwal AK, Hay JN. The characterization of polyvinyl butyral by thermal
analysis. Thermochim Acta 2002;391:245–55.
[72] Miller DC, Kempe MD, Glick SH, Kurtz SR. Creep in photovoltaic modules:
Examining the stability of polymeric materials and components. Conf Rec IEEE
Photovolt Spec Conf 2010:262–8. doi:10.1109/PVSC.2010.5615832.
[73] Cuddihy EF, Coulbert CD, Liang RH, Gupta a., Willis P, Baum B. Applications of
ethylene vinyl acetate as an encapsulation material for terrestrial photovoltaic
modules. NASA STI/Recon Tech Rep N 1983;83:29809.
[74] Agroui K, Collins G, Farenc J. Measurement of glass transition temperature of
crosslinked EVA encapsulant by thermal analysis for photovoltaic application.
Renew Energy 2012;43:218–23. doi:10.1016/j.renene.2011.11.015.
[75] Chen Z, Fang P, Wang H, Zhang S, Wang S. Property of ethylene vinyl acetate
copolymer in melting processing. J Appl Polym Sci 2006;101:2022–6.
doi:10.1002/app.23767.
[76] Haward RN, Thackray G. The Use of a Mathematical Model to Describe Isothermal
Stress-Strain Curves in Glassy Thermoplastics. Proc R Soc A Math Phys Eng Sci
1968;302:453–72. doi:10.1098/rspa.1968.0029.
[77] Delincé D, Callewaert D, Belis J, Van Impe R. Post-breakage behaviour of laminated
glass in structural applications. Proc Challenging Glas Pg 2008:459–67.
[78] Meijer HEH, Govaert LE. Mechanical performance of polymer systems: The relation
between structure and properties. Prog Polym Sci 2005;30:915–38.
doi:10.1016/j.progpolymsci.2005.06.009.
[79] Keller U, Mortelmans H. Adhesion in Laminated Safety Glass – What makes it
work? Glas. Process. Days, 1999, p. 353–6.
[80] Briccoli Bati S, Fagone M, Ranocchiai G. The elastic behavior of a rubber-like
material for composite glass. XIX Congr Dell’Associazione Ital Di Mecc Teor Ed
Appl 2009:1–10.
[81] Meissner M, Sackmann V. On the effect of artificial weathering on the shear bond
and the tear strength of two different interlayers of laminated glass.pdf. ISAAG -
Int. Symp. Appl. Archit. Glas., 2006, p. 223–33.
[82] Kott A. Zum Trag- und Resttragverhalten von Verbundsicherheitsglas 2006:189.
[83] Weller B, Kothe C, Kothe M, Wunsch J. Thermo Mechanical Behaviour of Polymeric
Interlayer Materials. Glas Perform Days 2009:734–7.
[84] Bolton NP, Smith WN. Laminated safety glass 1991.
[85] Weller B, Kothe C, Kothe M. Thermal Stability of Polymeric Interlayer Materials.
Challenging Glas 2 - Conf Archit Struct Appl Glas 2010:507–16.
[86] MacKnight WJ, Earnest TR. The structure and properties of ionomers. J Polym Sci
Macromol Rev 1981;16:41–122. doi:10.1002/pol.1981.230160102.
[87] Hausmann K, Hayes RA, Pesek SC, Shaffer WA, Smith CA. High-clarity ionomer
compositions and articles comprising the same 2010.
[88] Biolzi L, Cattaneo S, Rosati G. Progressive damage and fracture of laminated glass
beams. Constr Build Mater 2010;24:577–84. doi:10.1016/j.conbuildmat.2009.09.007.
[89] Calderone I, Davies PS, Bennison SJ, Xiaokun H, Gang L. Effective laminate
thickness for the design of laminated glass. Glas Perform Days 2009:1–5.
[90] Monserez S. Impact behavior of laminated glass for use in large glass facades.
Universiteit Gent, 2013.
[91] Karagiannis A, Szlosek ML, Farhoudi Y. High tg monolithic poly(vinyl acetal) sheet
2016.
[92] Snedeker RH, Garty KT, Skiermont FJ. Glass-polycarbonate resin laminates 1972.
[93] Savineau GF, Sauvineau GF, Savineau GF, Sauvineau GF. Fundamentals of
Laminating Process and Quality Requirements. Glas Process Days 1997:13–5.
[94] Adams RD, Wake WC. Adhesive Joints in Engineering. London and New York:
Elsevier Applied Science; 1984.
225/237
[95] ASTM-D3983-98. Standard test method for measuring strength and shear modulus
of nonrigid adhesives by the thick-adherend tensile-lap specimen. 2011.
[96] Louter C. Metal-to-glass bonding properties of an acrylate adhesive (DELO GB368)
and an ionoplast interlayer (SentryGlas) at 23 , -20 and 60 ° C. Glas. Perform. Days,
2009, p. 139–43.
[97] Vert M, Doi Y, Hellwich K-H, Hess M, Hodge P, Kubisa P, et al. Terminology for
biorelated polymers and applications (IUPAC Recommendations 2012). Pure Appl
Chem 2012;84. doi:10.1351/PAC-REC-10-12-04.
[98] Froli M, Lani L. Adhesion , creep and relaxation properties of PVB in laminated
safety glass. Glas. Perfomance Days, 2011, p. 1–4.
[99] Chapuis V, Pélisset S, Raeis-Barnéoud M, Li H-Y, Ballif C, Perret-Aebi L-E.
Compressive-shear adhesion characterization of polyvinyl-butyral and ethylene-
vinyl acetate at different curing times before and after exposure to damp-heat
conditions. Prog Photovoltaics Res Appl 2014;22:405–14. doi:10.1002/pip.2270.
[100] Suchy TA. Adhesion quality control of laminated safety glass using ultrasonic
velocity and attenuation measurements. University of Illinois at Urbana-
Champaign, 2011.
[101] Willham JEC, Brewster MH. Architectural panels with organic photovoltaic
interlayers and methods of forming the same 2015.
[102] Bennison SJ, Smith CA, Duser A Van, Jagota A. Structural Performance of
Laminated Glass Made with a “Stiff” Interlayer. Use Glas. Build., 100 Barr Harbor
Drive, PO Box C700, West Conshohocken, PA 19428-2959: ASTM International;
2002, p. 57-57–9. doi:10.1520/STP11056S.
[103] CEN. EN 1990 Eurocode : Basis of structural design. 2002.
[104] Galuppi L, Royer-Carfagni G. A homogenized model for the post-breakage tensile
behavior of laminated glass. Compos Struct 2016;154:600–15.
doi:10.1016/j.compstruct.2016.07.052.
[105] Delincé D, Callewaert D, Vanlaere W, Belis J, Depauw J. Plastic deformation of
polymer interlayers during post- breakage behavior of laminated glass. Int J Mod
Phys B 2008;22:5447–5452/.
[106] Castori G, Speranzini E. Structural analysis of failure behavior of laminated glass.
Compos Part B Eng 2017;125:89–99. doi:10.1016/j.compositesb.2017.05.062.
[107] Aenelle ML, Pelayo F, Ismael G. Damping of laminated glass beams with polyvinil
butiral interlayer. 7th Int. Oper. Modal Anal. Conf., 2017.
[108] Sable L, Skukis E, Japins G, Kalnins K. Correlation between Numerical and
Experimental Tests of Laminated Glass Panels with Visco-elastic Interlayer.
Procedia Eng 2017;172:945–52. doi:10.1016/j.proeng.2017.02.107.
[109] Delincé D, Belis J. Experimental assessment of polymers in glass constructions.
COST Action TU0905 Mid-term Conf. Struct. Glas., CRC Press/Balkema; 2013, p.
323–30.
[110] Piscitelli LR. Campagna sperimentale sul vetro stratificato: studio del
comportamento viscoelastico ed applicazione a un caso studio. Università degli
Studi di Firenze, 2014.
[111] Wu P, Zhou D, Liu W, Wan L, Liu D. Elasticity solution of two-layer beam with a
viscoelastic interlayer considering memory effect. Int J Solids Struct 2016;94–95:76–
86. doi:10.1016/j.ijsolstr.2016.05.007.
[112] Galuppi L, Royer-Carfagni G. The design of laminated glass under time-dependent
loading. Int J Mech Sci 2013;68:67–75. doi:10.1016/j.ijmecsci.2012.12.019.
[113] Galuppi L, Royer-Carfagni GF. Effective thickness of laminated glass beams: New
expression via a variational approach. Eng Struct 2012;38:53–67.
doi:10.1016/j.engstruct.2011.12.039.
[114] Ferry JD. Viscoelastic Properties of Polymers. New York: 1980.
doi:10.1080/00914037208075296.
[115] Aklonis JJ, Shaw MT, William MJ. Introduction to polymer viscoelasticity. vol. 59.
226/237
1972. doi:10.1016/j.matchar.2007.05.008.
[116] AS 1288-2006 Australian Australian Standard : Glass in buildings - Selection and
installation 2006.
[117] ASTM E1300-16 : Standard Practice for Determining Load Resistance of Glass in
Buildings. West Conshohocken, PA: ASTM International; 2016. doi:10.1520/E1300-
16.
[118] Giovanna R, Zulli F, Andreozzi L, Fagone M. Test Methods for the Determination of
Interlayer Properties in Laminated Glass. J Mater Civ Eng 2017;29:04016268.
doi:10.1061/(ASCE)MT.1943-5533.0001802.
[119] International Organization for Standardization. ISO 6721:2011 Plastics --
Determination of dynamic mechanical properties. 2011.
[120] Widman JC. Recent Developments in Penetration Resistance of Windshield Glass
1965. doi:10.4271/650474.
[121] Rieser RG, Chabal J. Safety Performance of Laminated Glass Configurations 1967.
doi:10.4271/670912.
[122] Yeomans D. The pre-history of the curtain wall. Constr Hist 1998;14:59–82.
doi:10.2307/41601861.
[123] Barakat MA, Hetherington JG. New architectural forms to reduce the effects of blast
waves and fragments on structures. Struct Under Shock Impact V 1998;32:53–62.
[124] Katsuki K, Shibaoka K, Onishi S, Miwa T. Laminated sheet composed of thin glass
and plastic sheets. J Non Cryst Solids 1980;42:619–36.
[125] Sha Y, Hui CY, Kramer EJ, Garrett PD, Knapczyk JW. Analysis of adhesion and
interface debonding in laminated safety glass. J Adhes Sci Technol 2012;11:49–63.
doi:10.1163/156856197X01010.
[126] Pelfrene J, Van Dam S, Van Paepegem W. Numerical and experimental study of the
peel test for assessment of the glass-PVB interface properties in laminated glass.
Challenging Glas. 4 COST Action TU0905 Final Conf., CRC Press; 2014, p. 513–9.
doi:10.1201/b16499-74.
[127] Dural E. Experimental and numerical treatment of delamination in laminated glass
plate structures. J Reinf Plast Compos 2016;35:56–70. doi:10.1177/0731684415613163.
[128] Santarsiero M, Louter C, Nussbaumer A. Laminated connections under tensile load
at different temperatures and strain rates. Int J Adhes Adhes 2017;79:23–49.
doi:10.1016/j.ijadhadh.2017.09.002.
[129] Serafinavicius T, Lebet J-P, Louter C, Kuranovas A, Lenkimas T. The effects of
environmental impacts on durability of laminated glass plates with interlayers (SG,
EVA, PVB). Challenging Glas. 4 COST Action TU0905 Final Conf., CRC Press; 2014,
p. 455–62. doi:10.1201/b16499-66.
[130] Santarsiero M, Louter C, Nussbaumer A. Laminated connections for structural glass
applications under shear loading at different temperatures and strain rates. Constr
Build Mater 2016;128:214–37. doi:10.1016/j.conbuildmat.2016.10.045.
[131] Andreozzi L, Briccoli Bati S, Fagone M, Ranocchiai G, Zulli F. Analysis of
environmental damage consequences on viscous thermo elastic properties of a
polymer interlayer. XXI Congr. AIMETA Assoc. Ital. di Mecc. Teor. e Appl., Torino,
Italy: 2013.
[132] Andreozzi L, Zulli F, Briccoli Bati S, Fagone M, Ranocchiai G. Weathering of
laminated glass – Degradation of interlayer mechanical properties. Challenging
Glas. 4 COST Action TU0905 Final Conf., CRC Press; 2014, p. 419–24.
doi:10.1201/b16499-61.
[133] Andreozzi L, Briccoli Bati S, Fagone M, Ranocchiai G, Zulli F. Weathering action on
thermo-viscoelastic properties of polymer interlayers for laminated glass. Constr
Build Mater 2015;98:757–66. doi:10.1016/j.conbuildmat.2015.08.010.
[134] Huntsberger JR. Adhesion of Plasticized Polyvinyl Butyral to Glass. J Adhes
1981;13:107–29.
[135] Vallabhan C, Das YC, Ramasamudra M. Properties of PVB Interlayer Used in
227/237
Laminated Glass. J Mater Civ Eng 1992;4:71–6. doi:10.1061/(ASCE)0899-
1561(1992)4:1(71).
[136] El-Din NMS, Sabaa MW. Thermal degradation of poly(vinylbutyral) laminated
safety glass. Polym Degrad Stab 1995;41:283–8. doi:10.1016/0141-3910(94)00118-R.
[137] Sobek W, Kutterer M, Messmer R. Shear stiffness of the interlayer in laminated
glass. Glas. Process. Days, Tampere, Finland: 1999, p. 360–5.
[138] D’Haene P. Flexural Stiffness of Laminated Glass Panels Under Dynamic
Conditions. GPD - Glas. Perform. Days, Tampere, Finland: 2003.
[139] Graf H, Schuler C, Albrecht G, Bucak O. The Influence of Various Support
Conditions on the Structural Behaviour of Laminated Glass. Glas Process Days
2003:408–11.
[140] Weller B, Wunsch J, Harth K. Experimental Study on Different Interlayer Materials
for Laminated Glass. Glas Process Days 2005:120–3.
[141] Delinc D, Belis J, Zarmati G, Parmentier B. Structural behaviour of laminated glass
elements – a step towards standardization. Glas Perform Days 2007:658–63.
[142] Callewaert D, Belis J, Van Impe R, Lagae G, Vanlaere W. Refined test set-up for pure
torsion of laminated glass. GPD - Glas. Perform. Days, Tampere, Finland: 2007, p.
118–21.
[143] D’Haene P, Savineau G. Mechanical properties of laminated safety glass - FEM
Study. Glas. Perform. Days, 2007, p. 594–8.
[144] Decourcelle R, Nugue JC, Levasseur F. Mechanical participation of interlayer on
laminated glass for building applications. Glas Perform Days 2009:718–21.
[145] Briccoli Bati S, Ranocchiai G, Reale C, Rovero L. Time-Dependent Behavior of
Laminated Glass. J Mater Civ Eng 2009;22:389–96. doi:10.1061/(ASCE)MT.1943-
5533.0000032.
[146] Biolzi L, Cagnacci E, Muciaccia G, Orlando M, Piscitelli LR. Experimental evaluation
on long-term behaviour of cold-bent laminated glass. XXVII A.T.I.V. Conf., Parma,
Italy: 2014.
[147] Hooper PA, Blackman BRK, Dear JP. The mechanical behaviour of poly (vinyl
butyral) at different strain magnitudes and strain rates. J Mater Sci 2012;47:3564–76.
[148] Andreozzi L, Briccoli Bati S, Fagone M, Ranocchiai G, Zulli F. Dynamic torsion tests
to characterize the thermo-viscoelastic properties of polymeric interlayers for
laminated glass. Constr Build Mater 2014;65:1–13.
doi:10.1016/j.conbuildmat.2014.04.003.
[149] Xu X, Liu B, Li Y. Experimental Studies on Viscoelasticity of Film Materials in
Laminated Glass Sheets. SAE Int J Mater Manuf 2015;8:2015-01–0709.
doi:10.4271/2015-01-0709.
[150] Nawar M, Salim H, Lafta M, Elshihy A. Quasi-Static Performance of Interlayer
Systems for Laminated Glass. J Mater Civ Eng 2016;28:04015084.
doi:10.1061/(ASCE)MT.1943-5533.0001348.
[151] Brinson HF, Brinson LC. Polymer engineering science and viscoelasticity: An
introduction, Second edition. 2015. doi:10.1007/978-1-4899-7485-3.
[152] Qi HJ, Boyce MC. Stress–strain behavior of thermoplastic polyurethanes. Mech
Mater 2005;37:817–39. doi:10.1016/j.mechmat.2004.08.001.
[153] Şerban DA, Weber G, Marşavina L, Silberschmidt V V., Hufenbach W. Tensile
properties of semi-crystalline thermoplastic polymers: Effects of temperature and
strain rates. Polym Test 2013;32:413–25. doi:10.1016/j.polymertesting.2012.12.002.
[154] Denault J, Vu-Khanh T, Foster B. Tensile properties of injection molded long fiber
thermoplastic composites. Polym Compos 1989;10:313–21. doi:10.1002/pc.750100507.
[155] van der Vegt AK, Govaert LE. Polymeren : van keten tot kunststof. 2003.
[156] Klompen ETJ. Mechanical properties of solid polymers. Eindhoven University, The
Netherlands, 2005.
[157] Hutchinson JM. Physical aging of polymers. Prog Polym Sci 1995;20:703–60.
doi:10.1016/0079-6700(94)00001-I.
228/237
[158] Behr RA, Minor JE, Norville SH. Structural Behavior of Architectural Laminated
Glass. J Struct Eng 1993;119:202–22. doi:10.1061/(ASCE)0733-9445(1993)119:1(202).
[159] Minor JE, Reznik PL, Fellow, ASCE. Failure Strengths of Laminated Glass. J Struct
Eng 1990;116:1030–9. doi:10.1061/(ASCE)0733-9445(1990)116:4(1030).
[160] Newmark NM, Siess CP, Viest IM. Tests and analysis of composite beams with
incomplete interaction. Proc Soc Exp Stress Anal 1951:75–92.
[161] Reissner E. Finite defections of sandwich plates. J Aeronaut Sci 1948;15:435–40.
[162] Wang C-T. Principle and application of complementary energy method for thin
homogeneous and sandwich plates and shells with finite defections. NACA Natl.
Advis. Comm. Aeronaut. TN 2620, Washington, D.C.: 1952.
[163] Pister KS, Dong SB. Elastic bending of layered plates. J Engrg Mech Div 1959;85:1–
10.
[164] Timoshenko S. History of Strength of Materials: With a Brief Account of the History
of Theory of Elasticity and Theory of Structures. Dover Publications; 1953.
[165] Minor JE. Developments in the Design of Architectural Glazing Systems. First Natl.
Struct. Eng. Conf. 1987 Prepr. Pap., Institution of Engineers, Australia; 1987, p. 77.
[166] Behr RA, Minor JE, Linden MP. Load Duration and Interlayer Thickness Effects on
Laminated Glass. J Struct Eng 1986;112:1441–53. doi:10.1061/(ASCE)0733-
9445(1986)112:6(1441).
[167] Vallabhan CVG, Minor JE, Nagalla SR. Stresses in Layered Glass Units and
Monolithic Glass Plates. J Struct Eng 1987;113:36–43. doi:10.1061/(ASCE)0733-
9445(1987)113:1(36).
[168] Vallabhan CVG, Das YC, Magdi M, Aşık MZ, Bailey JR. Analysis of Laminated
Glass Units. J Struct Eng 1993;119:1572–85. doi:10.1061/(ASCE)0733-
9445(1993)119:5(1572).
[169] Shutov AI, Novikov IA, Frank AN. Strength of Laminated Safety Glass for
Construction Purposes. Glas Ceram 2001;58:368–71. doi:10.1023/a:1012338931087.
[170] Pang S, Ma J, Ma Z, Liu X, Ren S. Correlation Between Temperature and Bending
Behavior of Laminated Glass. J Wuhan Univ Technol 2010;22:4.
[171] Louter C, Belis J, Bos F, Callewaert D, Veer F. Experimental investigation of the
temperature effect on the structural response of SG-laminated reinforced glass
beams. Eng Struct 2010;32:1590–9. doi:10.1016/j.engstruct.2010.02.007.
[172] Pankhardt K, Balázs GL. Temperature dependent load bearing capacity of
laminated glass panes. Civ Eng 2010;1:11–22. doi:10.3311/pp.ci.2010-1.02.
[173] Foraboschi P. Analytical model for laminated-glass plate. Compos Part B Eng
2012;43:2094–106. doi:10.1016/j.compositesb.2012.03.010.
[174] Galuppi L, Manara G, Royer Carfagni G. Practical expressions for the design of
laminated glass. Compos Part B Eng 2013;45:1677–88.
doi:10.1016/j.compositesb.2012.09.073.
[175] Aşık MZ. Laminated glass plates: Revealing of nonlinear behavior. Comput Struct
2003;81:2659–71. doi:10.1016/S0045-7949(03)00325-0.
[176] Aşık MZ, Tezcan S. A mathematical model for the behavior of laminated glass
beams. Comput Struct 2005;83:1742–53. doi:10.1016/j.compstruc.2005.02.020.
[177] Aşık MZ, Tezcan S. Laminated glass beams: Strength factor and temperature effect.
Comput Struct 2006;84:364–73. doi:10.1016/j.compstruc.2005.09.025.
[178] Aşık MZ, Bayar G. Free Vibration Analysis of the Laminated Curved Glass Beams.
2nd South East Eur. Conf. Comput. Mech., 2009, p. 256–62.
[179] Bedon C, Amadio C, Sinico A. Numerical and analytical investigation on the
dynamic buckling behavior of glass columns under blast. Eng Struct 2014;79:322–40.
doi:10.1016/j.engstruct.2014.08.024.
[180] Bedon C, Belis J, Luible A. Assessment of existing analytical models for the lateral
torsional buckling analysis of PVB and SG laminated glass beams via viscoelastic
simulations and experiments. Eng Struct 2014;60:52–67.
doi:10.1016/j.engstruct.2013.12.012.
229/237
[181] Bedon C, Louter C. Exploratory numerical analysis of SG-laminated reinforced glass
beam experiments. Eng Struct 2014;75:457–68. doi:10.1016/j.engstruct.2014.06.022.
[182] Bedon C, Amadio C. Buckling of flat laminated glass panels under in-plane
compression or shear. Eng Struct 2012;36:185–97. doi:10.1016/j.engstruct.2011.12.010.
[183] Amadio C, Bedon C. A buckling verification approach for monolithic and laminated
glass elements under combined in-plane compression and bending. Eng Struct
2013;52:220–9. doi:10.1016/j.engstruct.2013.02.022.
[184] Bedon C, Louter C. Finite Element analysis of post-tensioned SG-laminated glass
beams with adhesively bonded steel tendons. Compos Struct 2017;167:238–50.
doi:10.1016/j.compstruct.2017.01.086.
[185] Amadio C, Bedon C. Buckling of laminated glass elements in out-of-plane bending.
Eng Struct 2010;32:3780–8. doi:10.1016/j.engstruct.2010.08.022.
[186] Louter C, Cupac J, Lebet J-P. Exploratory experimental investigations on post-
tensioned structural glass beams. J Facade Des Eng 2014;2:3–18. doi:10.3233/FDE-
130012.
[187] Veer FA, Rijgersberg H, Ruytenbeek D, Louter PC, Zuidema J. Composite Glass
Beams, the Third Chapter. Glas Process Days 2003 2003:307–10.
[188] Feirabend S, Sobek W. Reinforced laminated safety glass (original title: Bewehrtes
Verbundsicherheitsglas). Stahlbau 2008;77:16–22. doi:10.1002/stab.200810027.
[189] Wellershoff F, Sedlacek G. Structural use of glass in hybrid elements, steel-glass-
beams, glass-GFRP-plates. Glas. Process. days, 2003, p. 268–70.
[190] Flinterhoff A. Load carrying behaviour of hybrid steel-glass beams in bending.
University of Dortmund, 2003.
[191] Freytag B. Glass-Concrete Composite Technology. Struct Eng Int 2004;14:111–7.
doi:10.2749/101686604777963991.
[192] Louter PC. Adhesively bonded reinforced glass beams. Heron 2007;52:31–57.
[193] Feirabend S, Sobek W. Reinforced laminated glass. Challenging Glas. Conf. Archit.
Struct. Appl., 2008, p. 469–77.
[194] Palumbo D, Palumbo M, Mazzuchelli M, Schiavon A, Rebeschini A. A New Roof for
the XIII th Century “ Loggia de Vicari ” ( Arquà Petrarca -PD – Italy ) Based on
Structural Glass Trusses : a Case Study. Proc Glas Process Days 2005 2005:1–3.
[195] Speranzini E, Agnetti S. Post-cracking behaviour of reinforced glass beams. Proc
COST Action TU 2013;905:285–92.
[196] Correia JR, Valarinho L, Branco FA. Post-cracking strength and ductility of glass-
GFRP composite beams. Compos Struct 2011;93:2299–309.
doi:10.1016/j.compstruct.2011.03.018.
[197] Speranzini E, Agnetti S. Flexural performance of hybrid beams made of glass and
pultruded GFRP. Constr Build Mater 2015;94:249–62.
doi:10.1016/j.conbuildmat.2015.06.008.
[198] Kozłowski M, Kadela M, Hulimka J. Numerical Investigation of Structural
Behaviour of Timber-Glass Composite Beams. Procedia Eng 2016;161:990–1000.
doi:10.1016/j.proeng.2016.08.838.
[199] Koz M, Hulimka J. Load-Bearing Capacity of Hybrid Timber-Glass Beams. Archit
Civ Eng Environ 2014;2:61–71.
[200] Kozłowski M. Experimental and numerical analysis of hybrid timber-glass beams.
Silesian University of Technology, 2014.
[201] Veer FA, van Liebergen MAC, de Vries SM. Designing and engineering transparent
building components with high residual strength. Proc 5th Glas Process Days,
Tampere, Finl 1997.
[202] Veer FA, Riemslag AC, Ting CN. Structurally Efficient Glass Laminated Composite
Beams. Glas Process Days 2001 2001:363–7.
[203] Louter C, Belis J, Veer F, Lebet JP. Structural response of SG-laminated reinforced
glass beams; experimental investigations on the effects of glass type, reinforcement
percentage and beam size. Eng Struct 2012;36:292–301.
230/237
doi:10.1016/j.engstruct.2011.12.016.
[204] Louter C. Fragile yet ductile, structural aspects of reinforced glass beams. Delft
University of Technology; 2011.
[205] Martens K, Caspeele R, Belis J. Performance of statically indeterminate reinforced
glass beams – Experimental comparison with determinate systems and effect of a
discontinuous glass section. Constr Build Mater 2017;146:251–9.
doi:10.1016/j.conbuildmat.2017.03.186.
[206] Martens K, Caspeele R, Belis J. Development of Reinforced & Post-Tensioned Glass
Beams - Review. J Struct Eng 2015;142:1–23. doi:10.1061/(ASCE)ST.1943-
541X.0001453.
[207] Louter C, Cupac J, Debonnaire M. Structural glass beams prestressed by externally
bonded tendons. Glas. Glob. Conf. Proc., 2014, p. 450–9.
[208] Veer FA. 10 years of Zappi research. Proc. Glas. Process. Days, 2005.
[209] Antonelli A, Cagnacci E, Giordano S, Orlando M, Spinelli P. Experimental and
Theoretical Analysis of C-FRP Reinforced Glass Beams. ISAAG 2008 - Int. Symp.
Appl. Archit. Glas., Munich: 2008.
[210] Cagnacci E, Orlando M, Spinelli P. Experimental investigation on Reinforced Glass
Beams. ICCST/7 - Seventh Int. Conf. Compos. Sci. Technol., Sharjah, United Arab
Emirates: 2008.
[211] Cagnacci E, Orlando M, Spinelli P. Experimental campaign and numerical
simulation of the behaviour of reinforced glass beams. Glas. Perform. Days,
Tampere (Fi): 2009, p. 484–7.
[212] Galuppi L, Royer-Carfagni G. Rheology of cold-lamination-bending for curved
glazing. Eng Struct 2014;61:140–52. doi:10.1016/j.engstruct.2014.01.003.
[213] Vakar LI, Gaal M. Cold bendable, laminated glass - New possibilities in design.
Struct Eng Int J Int Assoc Bridg Struct Eng 2004;14:95–7.
doi:10.2749/101686604777963946.
[214] Belis J, Inghelbrecht B, Van Impe R, Callewaert D. Cold bending of laminated glass
panels. Heron 2007;52:123–46.
[215] Baraldi D, Cecchi A, Foraboschi P. Broken tempered laminated glass: Non-linear
discrete element modeling. Compos Struct 2016;140:278–95.
doi:10.1016/j.compstruct.2015.12.050.
[216] Louter C. Bewehrte Glasträger. Stahlbau 2008;77:23–7. doi:10.1002/stab.200810028.
[217] Foraboschi P. Behavior and Failure Strength of Laminated Glass Beams. J Eng Mech
2007;133:1290–301. doi:10.1061/(ASCE)0733-9399(2007)133:12(1290).
[218] Oda J, Toru T, Yamamoto H. A Study on Impact Strength of Laminated Glass with
Initial Crack. Trans Japan Soc Mech Eng Ser A 2000;66:1308–13.
doi:10.1248/cpb.37.3229.
[219] Xu J, Sun Y, Liu B, Zhu M, Yao X, Yan Y, et al. Experimental and macroscopic
investigation of dynamic crack patterns in PVB laminated glass sheets subject to
light-weight impact. Eng Fail Anal 2011;18:1605–12.
doi:10.1016/j.engfailanal.2011.05.004.
[220] Xu X, Xu J, Chen J, Li P, Liu B, Li Y. Investigation of dynamic multi-cracking
behavior in PVB laminated glass plates. Int J Impact Eng 2017;100:62–74.
doi:10.1016/j.ijimpeng.2016.10.013.
[221] Huang X, Liu Q, Liu G, Zhou Z, Li G. Experimental investigation of multi-layered
laminated glass beams under in-plane bending. Struct Eng Mech 2016;60:781–94.
doi:10.12989/sem.2016.60.5.781.
[222] Ali A, Hosseini M, Sahari BB. A review and comparison on some rubber elasticity
models. J Sci Ind Res (India) 2010;69:495–500.
[223] Viola E. Scienza delle costruzioni. Pitagora; 1990.
[224] Treloar LRG. The physics of rubber elasticity. Oxford University Press, USA; 1975.
[225] Steinmann P, Hossain M, Possart G. Hyperelastic models for rubber-like materials:
Consistent tangent operators and suitability for Treloar’s data. Arch Appl Mech
231/237
2012;82:1183–217. doi:10.1007/s00419-012-0610-z.
[226] Nunes LCS. Mechanical characterization of hyperelastic polydimethylsiloxane by
simple shear test. Mater Sci Eng A 2011;528:1799–804. doi:10.1016/j.msea.2010.11.025.
[227] Bower AF. Applied Mechanics of Solids. CRC Press; 2009.
[228] Yeoh OH. Some Forms of the Strain Energy Function for Rubber. Rubber Chem
Technol 1993;66:754–71. doi:10.5254/1.3538343.
[229] Moreira DC, Nunes LCS. Comparison of simple and pure shear for an
incompressible isotropic hyperelastic material under large deformation. Polym Test
2013;32:240–8. doi:10.1016/j.polymertesting.2012.11.005.
[230] Tobajas R, Ibartz E, Gracia L. A comparative study of hyperelastic constitutive
models to characterize the behavior of a polymer used in automotive engines. Proc
2nd Int Electron Conf Mater 2016:A002. doi:10.3390/ecm-2-A002.
[231] Rivlin RS, Saunders DW. Large Elastic Deformations of Isotropic Materials. Collect.
Pap. R.S. Rivlin, New York, NY: Springer New York; 1997, p. 157–94.
doi:10.1007/978-1-4612-2416-7_12.
[232] Sadd MH. Elasticity: theory, applications, and numerics. Academic Press; 2009.
[233] Attard MM. Finite strain - isotropic hyperelasticity. Int J Solids Struct 2003;40:4353–
78. doi:10.1016/S0020-7683(03)00217-8.
[234] Hackett RM. Hyperelasticity Primer. Cham: Springer International Publishing; 2016.
doi:10.1007/978-3-319-23273-7.
[235] Mooney M. A Theory of Large Elastic Deformation. J Appl Phys 1940;11:582–92.
doi:10.1063/1.1712836.
[236] Ogden RW. Large Deformation Isotropic Elasticity - On the Correlation of Theory
and Experiment for Incompressible Rubberlike Solids. Proc R Soc A Math Phys Eng
Sci 1972;326:565–84. doi:10.1098/rspa.1972.0026.
[237] Bergstrom JS. Mechanics of solid polymers: theory and computational modeling.
William Andrew; 2015.
[238] Christensen R. Theory of viscoelasticity: an introduction. Elsevier; 2012.
[239] Perzyna P. Fundamental Problems in Viscoplasticity, 1966, p. 243–377.
doi:10.1016/S0065-2156(08)70009-7.
[240] Landel RF, Nielsen LE. Mechanical properties of polymers and composites. Crc
Press; 1993.
[241] Roylance D. Mechanics of materials. New York: Wiley; 1996.
[242] Simo JC, Hughes TJR. Computational inelasticity. vol. 7. Springer Science &
Business Media; 2006.
[243] Gant FS, Bower MV. Domain of Influence Method: A New Method for
Approximating Prony Series Coefficients and Exponents for Viscoelastic Materials. J
Polym Eng 1997;17:1–22. doi:10.1515/POLYENG.1997.17.1.1.
[244] Soussou JE, Moavenzadeh F, Gradowczyk MH. Application of Prony Series to
Linear Viscoelasticity. Trans Soc Rheol 1970;14:573–84. doi:10.1122/1.549179.
[245] Cost TL, Becker EB. A multidata method of approximate Laplace transform
inversion. Int J Numer Methods Eng 1970;2:207–19.
[246] Bradshaw RD, Brinson LC. A Sign Control Method For Fitting And Interconverting
Material Functions For Linearly Viscoelastic Solids "A Sign Control Method for
Fitting and Interconverting Material Functions for Linear Viscoelastic Solids
Introduction and Background. Mech Time Dep Matls 1997;1:85–108.
doi:10.1023/A:1009772018066.
[247] Tzikang C. Determining Viscoelastic Strain Data a Prony Material Series for a From
Time Varying. US Army Res Lab Veh Technol Dir Langley Res Center, Hampton,
Virginia 2000.
[248] Park SW, Kim YR. Fitting Prony-Series Viscoelastic Models with Power-Law
Presmoothing. J Mater Civ Eng 2001;13:26–32. doi:10.1061/(ASCE)0899-
1561(2001)13:1(26).
232/237
[249] Tarefder RA, Rahman ASMA, University of New Mexico A, Center SUT,
Technology O of the AS for R and, Transportation ND of, et al. Interconversion of
Dynamic Modulus to Creep Compliance and Relaxation Modulus: Numerical
Modeling and Laboratory Validation - Final Report 2016:61p.
[250] Schapery RA, Park SW. Methods of interconversion between linear viscoelastic
material functions. Part II—an approximate analytical method. Int J Solids Struct
1999;36:1677–99.
[251] Park SW, Schapery RA. Methods of interconversion between linear viscoelastic
material functions. Part I—A numerical method based on Prony series. Int J Solids
Struct 1999;36:1653–75.
[252] Mun S, Chehab GR, Kim YR. Determination of Time-domain Viscoelastic Functions
using Optimized Interconversion Techniques. Road Mater Pavement Des
2007;8:351–65. doi:10.1080/14680629.2007.9690078.
[253] Rouse PE. A Theory of the Linear Viscoelastic Properties of Dilute Solutions of
Coiling Polymers. J Chem Phys 1953;21:1272–80. doi:10.1063/1.1699180.
[254] Williams ML, Landel RF, Ferry JD. The Temperature Dependence of Relaxation
Mechanisms in Amorphous Polymers and Other Glass-forming Liquids. J Am Chem
Soc 1955;77:3701–3707. doi:10.1021/ja01619a008.
[255] Save MA, Massonnet CE. Plastic analysis and design of plates, shells and disks 1972.
[256] Służalec A. Theory of Metal Forming Plasticity. Berlin, Heidelberg: Springer Berlin
Heidelberg; 2004. doi:10.1007/978-3-662-10449-1.
[257] Bailey D. Ductile Fracture 1997;73:2003.
[258] Mirkhalaf SM, Andrade Pires FM, Simoes R. An elasto-viscoplastic constitutive
model for polymers at finite strains: Formulation and computational aspects.
Comput Struct 2016;166:60–74. doi:10.1016/j.compstruc.2016.01.002.
[259] Irgens F. Continuum mechanics. Springer Science & Business Media; 2008.
[260] CEN. EN ISO 12543-4, Glass in building. Laminated glass and laminated safety
glass. Part 4: Test methods for durability. 2011.
[261] Biolzi L, Cagnacci E, Orlando M, Piscitelli LR. Laminated glass specimen under
long-term push-out loads. XXVII A.T.I.V. Conf., 2012, p. 95–100.
[262] Cottone A, Giambanco G, Spada A. Valutazione della resistenza di adesione e della
lunghezza efficace di incollaggio nei giunti adesivi tramite le leggi dell’effetto scala.
XVIII GIMC Conf 2010:1–4.
[263] Moore DR, Turner S. Mechanical Evaluation Strategies for Plastics. Elsevier Science;
2001.
[264] International Organization for Standardization. ISO 527:2012 Plastics -
Determination of tensile properties. Genève, CH: 2012.
[265] Weber F. Curved glass structures. GPD - Glas. Perform. Days, 2009, p. 375–80.
[266] Ensslen F. Influences of weathering on the durability of laminated safety glass. Int.
Symp. Appl. Archit. Glas., Munich: 2006, p. 183–94.
[267] Wurm J. Glass Structures. Birkhauser; 2007.
[268] O’Regan C. Structural use of glass in buildings (Second edition) - The Institution of
Structural Engineers. 2014.
[269] Hai ND, Mutsuyoshi H. Structural behavior of double-lap joints of steel splice plates
bolted/bonded to pultruded hybrid CFRP/GFRP laminates. Constr Build Mater
2012;30:347–59. doi:10.1016/j.conbuildmat.2011.12.001.
[270] Visser R. Residual lifetime assessment of uPVC gas pipes. Enschede, Netherland:
2009. doi:10.3990/1.9789036529587.
[271] Morris JW. Laminating techniques 1966.
[272] Benedictus E. Verre feuilleté de sécurité Triplex (registered November 25, 1909).
405,881, 1906.
[273] Overend M, Zammit K. A computer algorithm for determining the tensile strength
of float glass. Eng Struct 2012;45:68–77. doi:10.1016/j.engstruct.2012.05.039.
233/237
[274] Brown WG, Research NRCCD of B. A Practicable Formulation for the Strength of
Glass and Its Special Application to Large Plates. National Research Council
Canada; 1974.
[275] Beason WL, Research TTUI for D, (U.S.) NSF. A Failure Prediction Model for
Window Glass. NTIS; 1980.
[276] Technology S of G, Gesellschaft DG. Glass Technology. The Society; 1985.
[277] Pagano NJ. Stress fields in composite laminates. Int J Solids Struct 1978;14:385–400.
doi:10.1016/0020-7683(78)90020-3.
[278] BS EN 572-9. 2004.
[279] DIN 53457. 1987.
[280] Charlton DJ, Yang J, Teh KK. A Review of Methods to Characterize Rubber Elastic
Behavior for Use in Finite Element Analysis. Rubber Chem Technol 1994;67:481–503.
doi:10.5254/1.3538686.
[281] Luible A. Stabilität von Tragelementen aus Glas 2004;3014:264. doi:10.5075/epfl-
thesis-3014.
[282] Yang LM, Shim VPW, Lim CT. A visco-hyperelastic approach to modelling the
constitutive behaviour of rubber. Int J Impact Eng 2000;24:545–60. doi:10.1016/S0734-
743X(99)00044-5.
[283] Kott A, Vogel T. Safety of laminated glass structures after initial failure. IABSE
Symp Rep 2004;88:67–72. doi:10.2749/222137804796302563.
[284] Lenci S, Clementi F, Warminski J. Nonlinear free dynamics of a two-layer composite
beam with different boundary conditions. Meccanica 2015;50:675–88.
doi:10.1007/s11012-014-9945-6.
[285] Ivanov I V., Velchev DS, Sadowski T, Kneć M. Computational Models of Laminated
Glass Plate under Transverse Static Loading (in: Altenbach H, Eremeyev V (eds)
Shell-Like Structures Non-classical Theories and Applications). Adv Struct Mater
Springer Berlin 2011:469–90. doi:10.1007/978-3-642-21855-2.
[286] Jia Z, Hui D, Yuan G, Lair J, Lau K, Xu F. Mechanical properties of an epoxy-based
adhesive under high strain rate loadings at low temperature environment. Compos
Part B Eng 2016;105:132–7. doi:10.1016/j.compositesb.2016.08.034.
[287] Chiu WK, Chalkley PD, Jones R. Effects of temperature on the shear stress-strain
behaviour of structural adhesives (FM73). Comput Struct 1994;53:483–9.
doi:10.1016/0045-7949(94)90095-7.
[288] Jia Z, Yuan G, Ma H, Hui D, Lau K. Tensile properties of a polymer-based adhesive
at low temperature with different strain rates. Compos Part B Eng 2016;87:227–32.
doi:10.1016/j.compositesb.2015.10.013.
[289] Giboreau A, Cuvelier G, Launay B. Rheological behaviour of three
biopolymer/water systems,with emphasis on yield stress and viscoelastic properties.
J Texture Stud 1994;25:119–38. doi:10.1111/j.1745-4603.1994.tb01321.x.
[290] Kolednik O. The yield stress gradient effect in inhomogeneous materials. Int J Solids
Struct 2000;37:781–808. doi:10.1016/S0020-7683(99)00060-8.
[291] Premamoy G. Polymer Science and Techology of Plastics and Rubbers. New Delhi:
McGraw-Hill; 1990.
[292] Graessley WW. The entanglement concept in polymer rheology. Entanglement
Concept Polym. Rheol., Berlin/Heidelberg: Springer-Verlag; n.d., p. 1–179.
doi:10.1007/BFb0031037.
[293] Elzière P, Dalle-Ferrier C, Creton C, Barthel É, Ciccotti M. Large strain viscoelastic
dissipation during interfacial rupture in laminated glass. Soft Matter 2017;13:1624–
33. doi:10.1039/C6SM02785G.
[294] Mittal KL, Kern W. Selected bibliography on adhesion measurement of films and
coatings. J Adhes Sci Technol 1987;1:247–62. doi:10.1163/156856187X00265.
[295] Mittal KL. Adhesion Measurement of Thin Films, Thick Films, and Bulk Coatings.
100 Barr Harbor Drive, PO Box C700, West Conshohocken, PA 19428-2959: ASTM
International; 1978. doi:10.1520/STP640-EB.
234/237
[296] Young WC, Budynas RG. Roark’s formulas for stress and strain. New York: 2002.
[297] UNI EN 12150: Vetro di silicato-calcico di sicurezza temperato termicamente. 2015.
[298] Ministero delle Infrastrutture e dei Trasporti. NTC - Norme Tecniche per le
Costruzioni (GU n. 29 04/02/2008 - Suppl. Ordinario n. 30), 2008.
[299] Ministero delle Infrastrutture e dei Trasporti. circolare applicativa del 2 febbraio
2009 , n. 617 per le “NTC - Norme Tecniche per le Costruzioni” (GU n.47 26/02/2009
- Suppl. Ordinario n.27). 2009.
[300] PREN 16612:2013 - Glass In Building - Determination Of The Load Resistance Of
Glass Panes By Calculation And Testing. Comite Europeen de Normalisation; 2013.
[301] Palmgren AZ. Die Lebensdauer von Kugellagcrn. Z Ver Deutsch Ing 1924;68:339.
[302] Miner MA. Cumulative Damage in Fatigue. J Appl Mech 1945;12:159.
[303] CEN. FprEN 16613:2016 Glass in building - Laminated glass and laminated safetu
glass - Determination of interlayer mechanical properties 2016.
[304] James AG, Green A, Simpson GM. Strain Energy Function of Rubber. I.
Characterization of Gum Vulcanizates. J Appl Polim Sci 1975;19:2033–58.
doi:10.1002/app.1975.070190723.
[305] Ericksen JL. Theoretical elasticity. Bull Am Math Soc 1955;61:362–4.
[306] Kott A, Vogel T. Remaining structural capacity of broken laminated safety glass.
Proc. 8th Int. Conf. Archit. Automot. Glas. (Glass Process. Days), Tampere, Finl.,
2003, p. 403–7.
[307] Belis J, Depauw J, Callewaert D, Delincé D, Van Impe R. Failure mechanisms and
residual capacity of annealed glass/SGP laminated beams at room temperature. Eng
Fail Anal 2009;16:1866–75. doi:10.1016/j.engfailanal.2008.09.023.
[308] Hagl A. Experimental and Numerical Analysis of Edge Seal Spacers of Insulated
Glass Units for Structural Sealant Glazing Applications 2012.
[309] Baldridge DB. Laminated glass having electrically operated instrument indicator
means embedded in the interlayer. US003317906A, 1964.
[310] Kenichi A, Naoki U, Hirofumi O. Interlayers for use in sound-insulating laminated
glasses. US005190826A, 1993.
[311] Klafka G, Wieme R. Interlayer structure for laminated glass. US006451435B1, 2002.
[312] Aben H, Anton J, Õis M, Viswanathan K, Chandrasekar S, Chaudhri MM. On the
extraordinary strength of Prince Rupert’s drops. Appl Phys Lett 2016;109:231903.
doi:10.1063/1.4971339.
[313] Hugh Chisholm, editor. The Encyclopaedia Britannica: A Dictionary of Arts,
Sciences, Literature and General Information (11th edition). Philadelphia:
Encyclopaedia Britannica Company; 1910.
[314] Cathers WP, Ferretti JS. Vacuum channel method of laminating glass sheets 1983.
[315] Moran JR. No autoclave process for forming a safety glass laminate 1996.
[316] Meneguzzo F, Ciriminna R, Albanese L, Pagliaro M. The great solar boom: a global
perspective into the far reaching impact of an unexpected energy revolution. Energy
Sci Eng 2015;3:499–509. doi:10.1002/ese3.98.
[317] Liao S-H, Jhuo H-J, Yeh P-N, Cheng Y-S, Li Y-L, Lee Y-H, et al. Single Junction
Inverted Polymer Solar Cell Reaching Power Conversion Efficiency 10.31% by
Employing Dual-Doped Zinc Oxide Nano-Film as Cathode Interlayer, 2015.
doi:10.1038/srep06813.
[318] Thompson BB, D’Errico JJ. Polymer interlayers comprising uv absorbers 2015.
[319] Kurtz S. Progress and Strategies for Testing of Materials for Solar Panels, Shangai,
China: 2017.
[320] Tripathy AR, Chen W, Kukureka SN, MacKnight WJ. Novel poly(butylene
terephthalate)/poly(vinyl butyral) blends prepared by in situ polymerization of
cyclic poly(butylene terephthalate) oligomers. Polymer (Guildf) 2003;44:1835–42.
doi:10.1016/S0032-3861(03)00029-6.
[321] Alvarez G, Flores JJ, Aguilar JO, Gómez-Daza O, Estrada CA, Nair MTS, et al.
235/237
Spectrally selective laminated glazing consisting of solar control and heat mirror
coated glass: preparation, characterization and modelling of heat transfer. Sol
Energy 2005;78:113–24. doi:10.1016/j.solener.2004.06.021.
[322] Kuhn TS. The Structure of Scientific Revolutions. vol. II. 1970. doi:10.1119/1.1969660.
[323] Wölfel E. Nachgiebiger Verbund Eine Näherungslösung und deren
Anwendungsmöglichkeiten. Stahlbau 1987;56:173–80.
[324] Nonlinear Viscoelasticity. Polym. Eng. Sci. Viscoelasticity, Boston, MA: Springer US;
2008, p. 327–64. doi:10.1007/978-0-387-73861-1_10.
[325] Koutsawa Y, Daya EM. Static and free vibration analysis of laminated glass beam on
viscoelastic supports. Int J Solids Struct 2007;44:8735–50.
doi:10.1016/j.ijsolstr.2007.07.009.
[326] Galuppi L, Royer-Carfagni G. The effective thickness of laminated glass:
Inconsistency of the formulation in a proposal of EN-standards. Compos Part B Eng
2013;55:109–18. doi:10.1016/j.compositesb.2013.05.025.
[327] El-Shami MM, Norville SH, Ibrahim YE. Stress analysis of laminated glass with
different interlayer materials. Alexandria Eng J 2012;51:61–7.
doi:10.1016/j.aej.2012.07.008.
[328] Molnár G, Vigh LG, Stocker G, Dunai L. Finite element analysis of laminated
structural glass plates with polyvinyl butyral (PVB) interlayer. Period Polytech Civ
Eng 2012;56:35–42. doi:10.3311/pp.ci.2012-1.04.
[329] Pelayo F, Lopez-Aenlle M, Ismael G, Fernandez-Canteli A. Buckling of multilayered
laminated glass beams: Validation of the effective thickness concept. Compos Struct
2017;169:2–9. doi:10.1016/j.compstruct.2017.01.040.
[330] Gräf H, Albrecht G, Sackmann V, Schuler C, Bucak Ö. Structural behaviour of point-
supported and clamped laminated glass. Struct Eng Int 2004;14:129–33.
[331] Freddi F, Royer-Carfagni G, Silvestri M. Full-Scale Experiments for Point-Fixing
Frameless Laminated Glass. Int J Appl Glas Sci 2010;1:257–72. doi:10.1111/j.2041-
1294.2010.00028.x.
[332] Pottmann H, Schiftner A, Bo P, Schmiedhofer H, Wang W, Baldassini N, et al.
Freeform surfaces from single curved panels. ACM Trans Graph 2008;27:76.
doi:10.1145/1360612.1360675.
[333] Schiftner A, Baldassini N, Bo P, Pottmann H. Architectural freeform structures from
single curved panels. Adv Archit Geom 2008;3:45–8.
[334] Eigensatz M, Kilian M, Schiftner A, Mitra NJ, Pottmann H, Pauly M. Paneling
architectural freeform surfaces. ACM Trans Graph 2010;29:45.
doi:10.1145/1778765.1778782.
[335] Press WH, Flannery BP, Teukolsky SA, Vetterling WT. Numerical recipes. vol. 3.
cambridge University Press, cambridge; 1989.
[336] Boyce MC, Parks DM, Argon AS. Large inelastic deformation of glassy polymers.
Part I: rate dependent constitutive model. Mech Mater 1988;7:15–33.
[337] Arruda EM, Boyce MC. A three-dimensional constitutive model for the large stretch
behavior of rubber elastic materials. J Mech Phys Solids 1993;41:389–412.
[338] Wu PD, Van Der Giessen E. On improved network models for rubber elasticity and
their applications to orientation hardening in glassy polymers. J Mech Phys Solids
1993;41:427–56.
[339] Ghorbel E. A viscoplastic constitutive model for polymeric materials. Int J Plast
2008;24:2032–58. doi:10.1016/j.ijplas.2008.01.003.
[340] Van Dommelen JAW, Parks DM, Boyce MC, Brekelmans WAM, Baaijens FPT.
Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline
polymers. J Mech Phys Solids 2003;51:519–41. doi:10.1016/S0022-5096(02)00063-7.
[341] Parks DM, Ahzi S. Polycrystalline plastic deformation and texture evolution for
crystals lacking five independent slip systems. J Mech Phys Solids 1990;38:701–24.
[342] Nikolov S, Doghri I. A micro/macro constitutive model for the small-deformation
behavior of polyethylene. Polymer (Guildf) 2000;41:1883–91.
236/237
[343] Drozdov AD, Gupta RK. Non-linear viscoelasticity and viscoplasticity of isotactic
polypropylene. Int J Eng Sci 2003;41:2335–61.
[344] Drabousky DP, Dwight D. Prony series representation and interconversion of
viscoelastic material functions of equine cortical bone. Case Western Reserve
University, 2009.
[345] Sarva SS, Deschanel S, Boyce MC, Chen W. Stress–strain behavior of a polyurea and
a polyurethane from low to high strain rates. Polymer (Guildf) 2007;48:2208–13.
doi:10.1016/j.polymer.2007.02.058.
[346] Blatz PJ, Ko WL. Application of Finite Elastic Theory to the Deformation of Rubbery
Materials. Trans Soc Rheol 1962;6:223–52. doi:10.1122/1.548937.
[347] Bourcier DP, Koran FA, Crommen JHL. Glass Lamination Process, 2003.
[348] Treloar LRG. The elasticity of a network of long-chain molecules. I. Trans Faraday
Soc 1943:36–41.
[349] Yeoh OH. Characterization of Elastic Properties of Carbon-Black-Filled Rubber
Vulcanizates. Rubber Chem Technol 1990;63:792–805. doi:10.5254/1.3538289.
[350] Clay SB. Characterization of crazing properties of polycarbonate. Blacksburg,
Virginia, 2000.
237/237