5
5
5
in Laminated Composites
A Meta-model Based Approach
Sudip Dey
Mechanical Engineering Department
National Institute of Technology, Silchar, India
Tanmoy Mukhopadhyay
Department of Engineering Science
University of Oxford
Oxford, UK
Sondipon Adhikari
Chair of Aerospace Engineering
College of Engineering
Swansea University, Wales, UK
p,
p,
A SCIENCE PUBLISHERS BOOK
A SCIENCE PUBLISHERS BOOK
MATLAB® and Simulink® are trademarks of The MathWorks, Inc. and are used with permission. The
MathWorks does not warrant the accuracy of the text or exercises in this book. This book’s use or
discussion of MATLAB® and Simulink® software or related products does not constitute endorsement
or sponsorship by The MathWorks of a particular pedagogical approach or particular use of the
MATLAB® and Simulink® software.
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts have been
made to publish reliable data and information, but the author and publisher cannot assume responsibility for the
validity of all materials or the consequences of their use. The authors and publishers have attempted to trace the
copyright holders of all material reproduced in this publication and apologize to copyright holders if permission to
publish in this form has not been obtained. If any copyright material has not been acknowledged please write and let
us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmitted,
or utilized in any form by any electronic, mechanical, or other means, now known or hereafter invented, includ-
ing photocopying, microfilming, and recording, or in any information storage or retrieval system, without written
permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.com
(http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Danvers,
MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and registration for a variety
of users. For organizations that have been granted a photocopy license by the CCC, a separate system of payment
has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used only for
identification and explanation without intent to infringe.
Library of Congress Cataloging-in-Publication Data
Names: Dey, Sudip, author. | Mukhopadhyay, Tanmoy, author. | Adhikari,
Sondipon, author.
Title: Uncertainty quantification in laminated composites : a meta-model
approach / Sudip Dey, Asst. Professor, Mechanical Engineering Department,
National Institute of Technology, Silchar, India, Tanmoy Mukhopadhyay,
Department of Engineering Science, University of Oxford, Oxford, UK,
Sondipon Adhikari, Chair of Aerospace Engineering, College of Engineering,
Swansea University, Wales, UK.
Description: Boca Raton, FL : CRC Press, Taylor & Francis Group, [2018] | "A
science publishers book." | Includes bibliographical references and index.
Identifiers: LCCN 2018022527 | ISBN 9781498784450 (hardback)
Subjects: LCSH: Laminated materials--Mathematical models.
Classification: LCC TA418.9.L3 D49 2018 | DDC 620.1/18--dc23
LC record available at https://lccn.loc.gov/2018022527
Who should read it and why? Natural and engineered materials have acquired
an unprecedented role in the history of human civilisation. Entire eras have been
named after predominant materials, such, as the stone age, bronze age, and iron age
to mention a few. Composite materials are the defining advanced materials in the
current era with the exceptional promise of applicability in a number of high-end
areas such as aerospace, marine, automotive, construction and defence sectors.
Composite materials convincingly demonstrate that their mechanical properties
can be tailored to specific engineering demands and therefore perform superiorly
compared to conventional metallic materials. From a design perspective, this
advantage arises through a fundamental mathematical fact that the mechanical
properties of composite materials are functions of significantly more parameters
compared to their metallic counterparts. However, the advantages of composite
materials come at a cost that the designers and analysts have to deal with many
parameters. This gives rise to two major problems. Firstly, considering many
parameters simultaneously makes the design process computationally more
expensive. Secondly, perhaps more importantly, the increase in the number of
parameters leads to an unavoidable escalation of uncertainty associated with these
parameters. One way to address both of these problems simultaneously is to use
metamodels which effectively ‘replace’ the original physics-based computationally
expensive model with a data-based computationally inexpensive model. The aim
of this book is to introduce predominant techniques in this direction.
This book is the first comprehensive text on the treatment of uncertainties in
the modelling and analysis of composite materials. The authors have drawn on their
considerable research experience to produce this book. The text is written from
an engineering standpoint, comprising fundamental and complex theories that are
relevant across a wide range of metamodeling techniques. The book introduces
faculties, researchers and students about the confluence of composite structure
theory, uncertainty modelling and propagation and metamodeling approaches.
The pedagogical objective of this book is to systematically present the latest
developments in the metamodeling techniques and explain how they can be used
in conjunction with composite structures. The focus has been on the mathematical
iv Uncertainty Quantification in Laminated Composites: A Meta-Model Approach
contributions and discussions at different times. The authors are grateful to Professor
J. E. Cooper (University of Bristol, UK), Dr R. Chowdhury (IIT Roorkee, India),
Dr A. Chakrabari (IIT Roorkee, India), Dr S. Chakraborty (University of Notre
Dame, USA), Professor S. K. Sahu (NIT Rourkela, India), Dr. G. Li (Princeton
University, USA), Professor H. Rabitz (Princeton University, USA), Professor
E. Carrera and Dr A. Pagani (Politecnico di Torino, Italy), Dr C. Scarth (University
of Bath, UK) and Professor R. Banerjee (City University London, UK) for
many stimulating discussions contributing towards the intellectual content of
this book. The authors would like to gratefully acknowledge the contribution of
Dr S. Naskar (University of Aberdeen, UK) in preparing the initial two chapters of
this book. Beside the names taken here, we are thankful to many colleagues, fellow
researchers and students working in this field of research around the world, whose
name cannot be listed here for page limitations. The lack of explicit mentions by
no means implies that their contributions are any less. The opinions presented in
the book are entirely ours, and none of our colleagues, students, collaborators and
mentors has any responsibility for any shortcomings.
We have been fortunate to receive grants from various companies, charities
and government organisations including an Advanced Research Fellowship from
UK Engineering and Physical Sciences Research Council (EPSRC), the Wolfson
research merit award from The Royal Society, the Philip Leverhulme Prize from The
Leverhulme Trust, research grants from Embraer, NRN Wales and a Zienkiewicz
Scholarship from Swansea University. Without these funding, it would have been
impossible to conduct the works leading to this book.
Sudip Dey
Tanmoy Mukhopadhyay
Sondipon Adhikari
February 2018
.
Contents
Preface iii
1. Introduction 1
2. Uncertainty Quantification in Composite Structures: An Overview 17
3. Stochastic Dynamic Analysis of Laminated Composite Plates 37
4. Stochastic Free Vibration Analysis of Pre-twisted Singly Curved 67
Composite Shells
5. Stochastic Dynamic Analysis of Doubly Curved Composite Shells 90
6. Effect of Operational Uncertainties on the Stochastic Dynamics 122
of Composite Laminates
7. Effect of Environmental Uncertainties on the Free Vibration Analysis 145
of Composite Laminates
8. Effect of Cutout on the Stochastic Free Vibration Analysis of Twisted 165
Composite Panels
9. Stochastic Dynamic Stability Analysis of Composite Laminates 191
10. Uncertainty Quantification for Skewed Laminated Soft-core Sandwich 220
Panels
11. Fuzzy Uncertainty Quantification in Composite Laminates 250
12. Comparative Assessment of Kriging Variants for Stochastic Analysis 272
of Composites
13. A Comparative Study on Metamodel Based Stochastic Analysis of 302
Composite Structures
Bibliography 331
Index 363
CHAPTER 1
Introduction
“Uncertainty is the only certainty there is, and knowing how to live
with insecurity is the only security”.
― John Allen Paulos
If we think deeply, sun-rise or sun-set never ever occurs. In reality, there is no new
year or end of year and thus there is no century. Actually, all of us have framed
everything as per our convenience. We always try to formulate our surroundings
in the form of a theory-law-hypothesis or an equation or algorithm. The perpetuity
of these deterministic frameworks based on the model of certainty does not exist in
the universal-continuum of time scale. The Sun and earth rotate to their own tune
of rotation. In many of the cases, we cannot even identify the variability/change
in pattern due to the unimaginable vastness of such systems both in terms of space
and time. Since the dawn of civilization all models introduced are meant for some
specific purposes, such as the solar model initially proposed by Pythagoras and
Aristotle (500 BC) was meant for the purpose of assessment of weather needed
for crop production. Later on, Ptolemy (300 BC) proposed a new solar model
which claimed to have corrected the previous solar model. Even though the earlier
models served the intended purpose to some extent, both the models are found to be
wrong as compared to the present solar model (refer to Fig. 1.1). Thus there exists
a deficiency in the accuracy of our understanding of the physical systems in many
cases. We human beings have formulated all the uncertain facts and transformed
them into the cage of certainty for our convenience. Hence, many of the theories,
laws, hypothesis, and formulae (from Aristotle to Hawking) become questionable
if we change its domain, assumptions and boundary conditions.
The convenient simplified models have their limitations in accurate predictions
and often differ from the true values. For centuries, researchers have been abiding
2 Uncertainty Quantification in Laminated Composites: A Meta-Model Approach
a b
by and spending their valuable time either to demonstrate or to validate those certain
(/deterministic) models, which are actually uncertain in the domain of space and
time. However, uncertainty being an inevitable characteristics of this universe, it is
almost impossible to model a system accurately in a deterministic manner. Rather,
it is a rational idea to try to quantify the effect of such uncertainty by considering it
as an integral part of the model and deal with the consequences. This approach of
modelling physical systems equips us to be prepared for the possible combination of
outcomes by providing a detailed account of the variation from deterministic values
of responses (i.e., the output quantity of interest). In case of engineering systems,
if the design is carried out by considering the effect of source-uncertainties and
prospective service-life conditions (such as environmental effects and damages)
instead of blatantly avoiding them, the possibility of failure can be minimized (/
controlled) based on a strong scientific foundation.
Let us take a simple example to understand the influence of source-uncertainty.
When we flip a fair coin, the probability of getting head or tail is 50%; but what
happens if there is a soft muddy floor. It may lead to some possibilities wherein
the coin will fall neither head nor tail, i.e., the coin may fall vertical. In such
cases, it will not be 50% probability of getting either head or tail. Moreover, due
to manufacturing uncertainties the coin may not be perfectly unbiased. This may
also deviate the outcome from the general expectation. Likewise, the probability
of getting one in dice having six faces is 1/6. But if the similar conditions (such
as soft muddy surface, manufacturing uncertainty of the dice) are imposed here,
there can be some possibilities wherein the dice will not show any specific face, or
the probability of different outcomes would vary. Hence, the exact probability of
getting one in those cases can’t be predicted. In a similar fashion, all the physical
or mathematical models are formulated based on certain boundary conditions and
assumptions, which are often not strictly valid.
Introduction 3
Tim e
...
what could happen, but how likely it is to happen, i.e., the probability of occurrence.
x1(ro)
x2(ro) y(ro)
x3(ro)
Stocbastic
system
... J ( \ ~utput probability
\_stribution
Fig. 1.6 Schematic representation for Monte Carlo simulation based analysis of a stochastic system
with three input parameters (xi, i = 1, 2, 3) and one output parameter (y). Here w– represents the
stochastic character of a parameter.
Introduction 9
– ) and y(x(w
Fig. 1.7 Metamodel based analysis of stochastic systems. (Here x(w – )) are the symbolic
representation of stochastic input parameters and output responses respectively. w – denotes the
stochasticity of parameters.)
preceding sections, this book is written to address both the aspects of modeling
different forms of uncertainties in composites as well as efficient computational
approaches for uncertainty propagation. After providing a general overview
of uncertainty quantification and stochastic analysis of composite structures in
Chapters 1–2, rest of the chapters in this book are concentrated on various specific
aspects related to uncertainty analysis of laminated composites and metamodel
based algorithms for uncertainty quantification. The Chapters 3–10 of this book
are focused on probabilistic approach of uncertainty quantification in composite
laminates and sandwich structures, while the Chapter 11 deals with a non-
probabilistic approach of uncertainty quantification in composites. The aspects of
uncertainty and sensitivity of the material, geometric, environmental and operational
factors for different responses of composite structures are analyzed in the Chapters
3–11. Various application-specific requirements (such as cutouts, twist and skewed
geometry) in modern high-performance structural systems are analyzed. The effect
of noise is analyzed in surrogate based uncertainty quantification algorithms for
composites. Chapter 12 provides a critical assessment of different kriging model
variants for the uncertainty quantification of composite structures. Different meta-
models are used for formulating the uncertainty quantification schemes in various
chapters of this book. Finally, comparative assessment of these meta-models is
presented in Chapter 13. A concise summary explaining the contributions of each
of the chapters in this book is given below.
The Chapter 2 of this book gives an overview of uncertainty quantification and
a general review of the literature related to uncertainty quantification in composite
structures. Different probabilistic and non-probabilistic methods of stochastic
structural analysis are briefly presented in this chapter. After providing a concise
review of the deterministic models for analyzing composite structures, the recent
works on the aspect of uncertainty quantification in composite structures are
discussed. The present standing of research in this area is assessed critically and the
contribution of this book is justified in that context. Even though a comprehensive
literature review is presented in Chapter 2, Chapters 3–13 also cite relevant studies
specific to the topic of a particular chapter.
The effect of material and geometric uncertainty on the dynamic responses
of composite plates is investigated in Chapter 3. A bottom up surrogate based
approach is employed to quantify the variability in free vibration responses of
composite cantilever plates due to uncertainty in ply orientation angle, elastic
modulus and mass density. The finite element method is employed incorporating
effects of transverse shear deformation based on Mindlin’s theory in conjunction
with a random variable approach. Parametric studies are carried out to determine
the stochastic frequency response functions (SFRF) along with stochastic natural
frequencies and modeshapes. In this study, a surrogate based approach using General
High Dimensional Model Representations (GHDMR) is employed for achieving
computational efficiency in quantifying uncertainty. This chapter also presents
an uncertainty quantification scheme using commercial finite element software
Introduction 13
analyses are carried out to illustrate the results of the first three stochastic natural
frequencies and buckling load.
A non-probabilistic uncertainty propagation approach (fuzzy) for composites is
presented in Chapter 11. Probabilistic descriptions of uncertain model parameters
are not always available due to lack of data. This chapter investigates on the
uncertainty propagation in dynamic characteristics (such as natural frequencies,
frequency response function and mode shapes) of laminated composite plates by
using fuzzy approach. A non-intrusive Gram–Schmidt polynomial chaos expansion
(GPCE) method is adopted in the uncertainty propagation, wherein the parameter
uncertainties are represented by fuzzy membership functions. A domain in the
space of input data at zero-level of membership functions is mapped to a zone of
output data with the parameters determined by D-optimal design. The obtained
meta-model (GPCE) can also be used for higher α-levels of fuzzy membership
function. The most significant input parameters such as ply orientation angle, elastic
modulus, mass density and shear modulus are identified and then fuzzified. Fuzzy
analysis of the first three natural frequencies is presented to illustrate the results
and its performance. The proposed fuzzy approach is applied to the problem of
fuzzy modal analysis for frequency response function of a simplified composite of
cantilever plates. The fuzzy mode shapes are also depicted for a typical laminate
configuration. The GPCE based approach is found more efficient compared to the
conventional global optimization approach in terms of computational time and cost.
Chapter 12 presents a critical comparative assessment of Kriging model
variants for surrogate based uncertainty propagation considering stochastic natural
frequencies of laminated composite shells. The five Kriging model variants
studied here are: Ordinary Kriging, Universal Kriging based on pseudo-likelihood
estimator, Blind Kriging, Co-Kriging and Universal Kriging based on marginal
likelihood estimator. First three stochastic natural frequencies of the composite
shell are analysed by using a finite element model that includes the effects of
transverse shear deformation based on Mindlin’s theory in conjunction with a
layer-wise random variable approach. The comparative assessment is carried out
to address the accuracy and computational efficiency of the five Kriging model
variants. Comparative performance of different covariance functions is also
studied. Subsequently the effect of noise in uncertainty propagation is addressed
by using the Stochastic Kriging. Representative results are presented for both
individual and combined stochasticity in layer-wise input parameters to address
performance of various Kriging variants for low dimensional and relatively higher
dimensional input parameter spaces. The error estimation and convergence studies
are conducted with respect to original Monte Carlo Simulation to justify merit
of the present investigation. The study reveals that Universal Kriging coupled
with marginal likelihood estimate yields the most accurate results, followed by
Co-Kriging and Blind Kriging. As far as computational efficiency of the Kriging
models is concerned, it is observed that for high-dimensional problems, CPU time
required for building the Co-Kriging model is significantly less as compared to
other Kriging variants.
16 Uncertainty Quantification in Laminated Composites: A Meta-Model Approach