Series Editor Professor Brian Derby,Professor of Materials Science Manchester Materials Science Centre,Grosvenor Street,Manchester,MI 7HS,UK Other titles published in this series: Fusion Bonding of Polymer Composites C.Ageorges and L.Ye Composite Materials D.D.L.Chung Titanium G.Lutjering and J.C.Williams Corrosion of Metals H.Kaesche Corrosion and Protection E.Bardal Intelligent Macromolecules for Smart Devices L.Dai Microstructure of Steels and Cast Irons M.Durand-Charre Phase Diagrams and Heterogeneous Equilibria B.Predel,M.Hoch and M.Pool Failure in Fibre Polymer Laminates M.Knops Publication due January 2005 Materials for Information Technology E.Zschech,C.Whelan and T.Mikolajick Publication due March 2005 Gallium Nitride Processing for Electronics,Sensors and Spintronics S.J.Pearton,C.R.Abernathy and F.Ren Publication due March 2005 Thermoelectricity J.P.Heremans,G.Chen and M.S.Dresselhaus Publication due August 2005 Computer Modelling of Sintering at Different Length Scales I.Pan Publication due October 2005 Computational Quantum Mechanics for Materials Engineers L.Vitos Publication due January 2006 Fuel Cell Technology N.Sammes Publication due January 2006
Series Editor Professor Brian Derby, Professor of Materials Science Manchester Materials Science Centre, Grosvenor Street, Manchester, M1 7HS, UK Other titles published in this series: Fusion Bonding of Polymer Composites C. Ageorges and L. Ye Composite Materials D.D.L. Chung Titanium G. Lu¨tjering and J.C. Williams Corrosion of Metals H. Kaesche Corrosion and Protection E. Bardal Intelligent Macromolecules for Smart Devices L. Dai Microstructure of Steels and Cast Irons M. Durand-Charre Phase Diagrams and Heterogeneous Equilibria B. Predel, M. Hoch and M. Pool Failure in Fibre Polymer Laminates M. Knops Publication due January 2005 Materials for Information Technology E. Zschech, C. Whelan and T. Mikolajick Publication due March 2005 Gallium Nitride Processing for Electronics, Sensors and Spintronics S.J. Pearton, C.R. Abernathy and F. Ren Publication due March 2005 Thermoelectricity J.P. Heremans, G. Chen and M.S. Dresselhaus Publication due August 2005 Computer Modelling of Sintering at Different Length Scales J. Pan Publication due October 2005 Computational Quantum Mechanics for Materials Engineers L. Vitos Publication due January 2006 Fuel Cell Technology N. Sammes Publication due January 2006
M.M.Kaminski Computational Mechanics of Composite Materials Sensitivity,Randomness and Multiscale Behaviour ②Springer
M.M. Kamin´ski Computational Mechanics of Composite Materials Sensitivity, Randomness and Multiscale Behaviour
M.M.Kaminski,MSc,PhD Division of Mechanics of Materials,Technical University of Lodz, Al.Politechniki 6,93-590 Lodz,Poland British Library Cataloguing in Publication Data Kaminiski,M.M. Computational mechanics of composite materials Sensitivity,randomness and multiscale behaviour.- (Engineering materials and processes) 1.Composite materials-Mathematical models 2.Mechanics, Applied-Data processing I.Title 620.001518 ISBN1852334274 Library of Congress Cataloging-in-Publication Data Kaminiski,M.M.(Marcin M.),1969- Computational mechanics of composite materials:sensitivity,randomness,and Multiscale behaviour M.M.Kaminski p.cm.-(Engineering materials and processes,ISSN 1619-0181) Includes bibliographical references and index. ISBN 1-85233-427-4 (alk.Paper) 1.Composite materials-Mechanical properties-Mathematical models.I.Title.II. Series. TA418.9.C6K3452002 621.1'892-dc21 2002033327 Apart from any fair dealing for the purposes of research or private study,or criticism or review,as permitted under the Copyright,Designs and Patents Act 1988,this publication may only be repro- duced,stored or transmitted,in any form or by any means,with the prior permission in writing of the publishers,or in the case of reprographic reproduction in accordance with the terms of licences issued by the Copyright Licensing Agency.Enquiries concerning reproduction outside those terms should be sent to the publishers. Engineering Materials and Processes ISSN 1619-0181 ISBN1-85233-427-4 Springer Science+Business Media springeronline.com Springer-Verlag London Limited 2005 Printed in the United States of America Whilst we have made considerable efforts to contact all holders of copyright material contained in this book,we have failed to locate some of these.Should holders wish to contact the Publisher,we will be happy to come to some arrangement with them. The use of registered names,trademarks,etc.in this publication does not imply,even in the absence of a specific statement,that such names are exempt from the relevant laws and regulations and therefore free for general use. The publisher makes no representation,express or implied,with regard to the accuracy of the infor- mation contained in this book and cannot accept any legal responsibility or liability for any errors or omissions that may be made. Typesetting:Electronic text files prepared by author 69/3830-543210 Printed on acid-free paper SPIN 10791506
M.M. Kamin´ski, MSc, PhD Division of Mechanics of Materials, Technical University of Ło´dz, Al. Politechniki 6, 93 - 590 Ło´dz, Poland British Library Cataloguing in Publication Data Kamin´ski, M.M. Computational mechanics of composite materials : Sensitivity, randomness and multiscale behaviour. — (Engineering materials and processes) 1. Composite materials — Mathematical models 2. Mechanics, Applied — Data processing I. Title 620′.001518 ISBN 1852334274 Library of Congress Cataloging-in-Publication Data Kamin´ski, M. M. (Marcin M.), 1969– Computational mechanics of composite materials: sensitivity, randomness, and Multiscale behaviour / M.M. Kamin´ski p. cm.—(Engineering materials and processes, ISSN 1619-0181) Includes bibliographical references and index. ISBN 1-85233-427-4 (alk. Paper) 1. Composite materials—Mechanical properties—Mathematical models. I. Title. II. Series. TA418.9.C6 K345 2002 621.1′892—dc21 2002033327 Apart from any fair dealing for the purposes of research or private study, or criticism or review, as permitted under the Copyright, Designs and Patents Act 1988, this publication may only be reproduced, stored or transmitted, in any form or by any means, with the prior permission in writing of the publishers, or in the case of reprographic reproduction in accordance with the terms of licences issued by the Copyright Licensing Agency. Enquiries concerning reproduction outside those terms should be sent to the publishers. Engineering Materials and Processes ISSN 1619-0181 ISBN 1-85233-427-4 Springer Science+Business Media springeronline.com © Springer-Verlag London Limited 2005 Printed in the United States of America Whilst we have made considerable efforts to contact all holders of copyright material contained in this book, we have failed to locate some of these. Should holders wish to contact the Publisher, we will be happy to come to some arrangement with them. The use of registered names, trademarks, etc. in this publication does not imply, even in the absence of a specific statement, that such names are exempt from the relevant laws and regulations and therefore free for general use. The publisher makes no representation, express or implied, with regard to the accuracy of the information contained in this book and cannot accept any legal responsibility or liability for any errors or omissions that may be made. Typesetting: Electronic text files prepared by author 69/3830-543210 Printed on acid-free paper SPIN 10791506
Acknowledgements Chapter 2 includes material,Civil-Comp Press,2001 previously published in M. Kaminski,"Homogenization Method in Stochastic Finite Element Analysis of some 1D Composite Structures",in Proc.8h Int.Conf.on Civil Structural Engineering Computational Technology,B.H.V.Topping,ed.,(paper no 60),Civil-Comp Press,Stirling, United Kingdom,2001.This material is reprinted with permission from Civil-Comp Press, Stirling,United Kingdom. Chapter 5 includes material,Civil-Comp Press,2001 previously published in L.Figiel,M. Kaminski,"Mechanical and Thermal Fatigue of Curved Composite Beams",in Proc.8 Int. Conf.on Civil Struct.Engineering Computational Technology,B.H.V.Topping,ed., (paper no 61),Civil-Comp Press,Stirling,United Kingdom,2001.This material is reprinted with permission from Civil-Comp Press,Stirling,United Kingdom. Chapter 5 includes material,Civil-Comp Press,2002 previously published in L.Figiel,M. Kaminski,"Numerical Analysis of Fatigue Damage Evolution in Composite Pipe Joints",in Proc.6 Int.Conf.on Computational Structures Technology,B.H.V.Topping and Z. Bittnar,eds.,(paper no 134),Civil-Comp Press,Stirling,United Kingdom,2002.This material is reprinted with permission from Civil-Comp Press,Stirling,United Kingdom. Chapter 7 includes material,Civil-Comp Press,2002 previously published in M. Kaminski,"Multiresolutional Homogenization Technique in Transient Heat Transfer for Unidirectional Composites",in Proc.6th Int.Conf.on Computational Structures Technology,B.H.V.Topping and Z.Bittnar,eds.,(paper no 138),Civil-Comp Press, Stirling,United Kingdom,2002.This material is reprinted with permission from Civil- Comp Press,Stirling,United Kingdom. Figures 2.38-2.41 are reproduced from M.Kaminski,M.Kleiber,Stochastic finite element method in random non-homogeneous media,in Numerical Methods in Engineering '96,J.A. Desideri et al.,eds.pp.35-41,1996,John Wiley Sons Limited.Reproduced with permission. Figures 2.98-2.119 are reproduced from M.Kaminski,M.Kleiber,Numerical homogenization of n-component composites including stochastic interface defects,Int.J. Num.Meth.Engrg.,47:1001-1027,2000,John Wiley Sons Limited.Reproduced with permission. Figures 7.16-7.21 and 7.54 -7.60 are reproduced from M.Kaminski,Stochastic perturbation approach to wavelet-based multiresolutional analysis,Num.Linear Algebra with Applications,11(4):355-370,2004,John Wiley Sons Limited.Reproduced with permission. Figures 2.65 and 2.66 reprinted from International Journal of Engineering Science,Vol 38. Kaminski,M.,Homogenized properties of n-components composites,pp.405-427, Copyright (2000),with permission from Elsevier. Figures 2.4-2.13 reprinted from International Journal of Solids and Structures,Vol 33. Kaminski,M.and Kleiber,M.,Stochastic structural interface defects in composite materials, pp.3035-3056,Copyright(1996),with permission from Elsevier. Figures 2.1 -2.3 and 2.30-2.40 reprinted from Computers and Structures,Vol 78, Kaminski,M.and Kleiber,M.,Perturbation-based stochastic finite element method for homogenization of two-component elastic composites,pp.811-826,Copyright(2000),with permission from Elsevier
Acknowledgements Chapter 2 includes material, © Civil-Comp Press, 2001 previously published in M. KamiĔski, “Homogenization Method in Stochastic Finite Element Analysis of some 1D Composite Structures”, in Proc. 8th Int. Conf. on Civil & Structural Engineering Computational Technology, B.H.V. Topping, ed., (paper no 60), Civil-Comp Press, Stirling, United Kingdom, 2001. This material is reprinted with permission from Civil-Comp Press, Stirling, United Kingdom. Chapter 5 includes material, © Civil-Comp Press, 2001 previously published in à. Figiel, M. KamiĔski, “Mechanical and Thermal Fatigue of Curved Composite Beams”, in Proc. 8th Int. Conf. on Civil & Struct. Engineering Computational Technology, B.H.V. Topping, ed., (paper no 61), Civil-Comp Press, Stirling, United Kingdom, 2001. This material is reprinted with permission from Civil-Comp Press, Stirling, United Kingdom. Chapter 5 includes material, © Civil-Comp Press, 2002 previously published in à. Figiel, M. KamiĔski, “Numerical Analysis of Fatigue Damage Evolution in Composite Pipe Joints”, in Proc. 6th Int. Conf. on Computational Structures Technology, B.H.V. Topping and Z. Bittnar, eds., (paper no 134), Civil-Comp Press, Stirling, United Kingdom, 2002. This material is reprinted with permission from Civil-Comp Press, Stirling, United Kingdom. Chapter 7 includes material, © Civil-Comp Press, 2002 previously published in M. KamiĔski, “Multiresolutional Homogenization Technique in Transient Heat Transfer for Unidirectional Composites”, in Proc. 6th Int. Conf. on Computational Structures Technology, B.H.V. Topping and Z. Bittnar, eds., (paper no 138), Civil-Comp Press, Stirling, United Kingdom, 2002. This material is reprinted with permission from CivilComp Press, Stirling, United Kingdom. Figures 2.38 – 2.41 are reproduced from M. KamiĔski, M. Kleiber, Stochastic finite element method in random non-homogeneous media, in Numerical Methods in Engineering ’96, J.A. Desideri et al., eds. pp. 35 – 41, 1996, © John Wiley & Sons Limited. Reproduced with permission. Figures 2.98 – 2.119 are reproduced from M. KamiĔski, M. Kleiber, Numerical homogenization of n-component composites including stochastic interface defects, Int. J. Num. Meth. Engrg., 47: 1001-1027, 2000, © John Wiley & Sons Limited. Reproduced with permission. Figures 7.16 – 7.21 and 7.54 – 7.60 are reproduced from M. KamiĔski, Stochastic perturbation approach to wavelet-based multiresolutional analysis, Num. Linear Algebra with Applications, 11(4): 355-370, 2004, © John Wiley & Sons Limited. Reproduced with permission. Figures 2.65 and 2.66 reprinted from International Journal of Engineering Science, Vol 38, KamiĔski, M., Homogenized properties of n-components composites, pp. 405-427, Copyright (2000), with permission from Elsevier. Figures 2.4 – 2.13 reprinted from International Journal of Solids and Structures, Vol 33, KamiĔski, M. and Kleiber, M., Stochastic structural interface defects in composite materials, pp. 3035-3056, Copyright (1996), with permission from Elsevier. Figures 2.1 – 2.3 and 2.30 – 2.40 reprinted from Computers and Structures, Vol 78, KamiĔski, M. and Kleiber, M., Perturbation-based stochastic finite element method for homogenization of two-component elastic composites, pp. 811-826, Copyright (2000), with permission from Elsevier
Acknowledgements Figures 4.1-4.9 reprinted from International Journal of Engineering Science,Vol 41, Kaminski,M.,Homogenization of transient heat transfer problems for some composite materials,pp.1-29,Copyright(2003),with permission from Elsevier. Figures 2.67 -2.69,2.77,2.78,2.88,2.89 and 4.17-4.52 reprinted from Computer Methods in Applied Mechanics and Engineering,Vol 192,Kaminski,M.,Sensitivity analysis of homogenized characteristics of some elastic composites,pp.1973-2005 Copyright(2003),with permission from Elsevier. Figures 3.1 -3.12 reprinted from Computational Materials Science,Vol 22,Figiel,L. Kaminski,M.,Effective elastoplastic properties of the periodic composites,pp.221-239, Copyright(2001),with permission from Elsevier. Figures 2.129-2.140 reprinted from Computational Materials Science,Vol 11,Kaminski, M.,Probabilistic bounds on effective elastic moduli for the superconducting coils,pp.252- 260,Copyright(1998),with permission from Elsevier. Figures 5.1-5.4 and 5.66-5.73 reprinted from International Journal of Fatigue,Vol 24, Kaminski,M.,On probabilistic fatigue models for composite materials,pp.477-495, Copyright(2002),with permission from Elsevier. Figures 7.2-7.15 are reprinted from Computational Materials Science,Vol.27,Kaminski M.,Wavelet-based homogenization of unidirectional multiscale composites,pp.613-622, Copyright(2001),with permission from Elsevier. Figures 7.30-7.43 and 7.46-7.53 reprinted from Computer Methods in Applied,Mechanics and Engineering,Kaminski,M.,Homogenization-based finite element analysis of unidirectional composites by classical and multiresolutional techniques,in press,Copyright (2005),with permission from Elsevier. Figure 2.49 reprinted from Kaminski,M.,Stochastic computational mechanics of composite materials",in Advances in Composite Materials and Structures VIl,de Wilde,W.P.,Blain, W.R.and Brebbia,C.A.,eds.,pp.219-228,Copyright(2000)WIT Press,Ashurst Lodge, Ashurst,Southampton,UK.Used with permission. Figures 2.42 and 4.10-4.16 reprinted from Archives of Applied Mechanics,Material sensitivity analysis in homogenization of the linear elastic composites,Kaminski,M., 71(10):679-694,2001,copyright Springer-Verlag Heidelberg.Used with permission. Figure 2.143 is reproduced from Kaminski,M.,Stochastic finite element in homogenization of linear elastic composites.Arch.Civil Engrg.3(XLVID):291-325,2001.Copyright property of the Polish Academy of Science.Used with permission. Figures in Chapter 6 are reproduced from Kaminski,M.,Stochastic reliability in contact problems for spherical particle reinforced composites.Journal of Theoretical and Applied Mechanics,3(39):539-562,2001.Used with permission. Figure 7.1 appeared in Kaminski,M.,Multiresolutional wavelet-based homogenization of random composites in Proceedings of the European Conference on Computational Mechanics Cracow,26-29 June 2001
vi Acknowledgements Figures 4.1 – 4.9 reprinted from International Journal of Engineering Science, Vol 41, KamiĔski, M., Homogenization of transient heat transfer problems for some composite materials, pp. 1-29, Copyright (2003), with permission from Elsevier. Figures 2.67 – 2.69, 2.77, 2.78, 2.88, 2.89 and 4.17 – 4.52 reprinted from Computer Methods in Applied Mechanics and Engineering, Vol 192, KamiĔski, M., Sensitivity analysis of homogenized characteristics of some elastic composites, pp. 1973-2005, Copyright (2003), with permission from Elsevier. Figures 3.1 – 3.12 reprinted from Computational Materials Science, Vol 22, Figiel, à., KamiĔski, M., Effective elastoplastic properties of the periodic composites, pp. 221-239, Copyright (2001), with permission from Elsevier. Figures 2.129 – 2.140 reprinted from Computational Materials Science, Vol 11, KamiĔski, M., Probabilistic bounds on effective elastic moduli for the superconducting coils, pp. 252- 260, Copyright (1998), with permission from Elsevier. Figures 5.1 – 5.4 and 5.66 –5.73 reprinted from International Journal of Fatigue, Vol 24, KamiĔski, M., On probabilistic fatigue models for composite materials, pp. 477-495, Copyright (2002), with permission from Elsevier. Figures 7.2 – 7.15 are reprinted from Computational Materials Science, Vol. 27, KamiĔski M., Wavelet-based homogenization of unidirectional multiscale composites, pp. 613-622, Copyright (2001), with permission from Elsevier. Figures 7.30-7.43 and 7.46-7.53 reprinted from Computer Methods in Applied, Mechanics and Engineering, KamiĔski, M., Homogenization-based finite element analysis of unidirectional composites by classical and multiresolutional techniques, in press, Copyright (2005), with permission from Elsevier. Figure 2.49 reprinted from KamiĔski, M., Stochastic computational mechanics of composite materials”, in Advances in Composite Materials and Structures VII, de Wilde, W.P., Blain, W.R. and Brebbia, C.A., eds., pp. 219 – 228, Copyright (2000) WIT Press, Ashurst Lodge, Ashurst, Southampton, UK. Used with permission. Figures 2.42 and 4.10 – 4.16 reprinted from Archives of Applied Mechanics, Material sensitivity analysis in homogenization of the linear elastic composites, KamiĔski, M., 71(10): 679 – 694, 2001, copyright Springer-Verlag Heidelberg. Used with permission. Figure 2.143 is reproduced from KamiĔski, M., Stochastic finite element in homogenization of linear elastic composites. Arch. Civil Engrg. 3(XLVII): 291-325, 2001. Copyright property of the Polish Academy of Science. Used with permission. Figures in Chapter 6 are reproduced from KamiĔski, M., Stochastic reliability in contact problems for spherical particle reinforced composites. Journal of Theoretical and Applied Mechanics, 3(39): 539-562, 2001. Used with permission. Figure 7.1 appeared in KamiĔski, M., Multiresolutional wavelet-based homogenization of random composites in Proceedings of the European Conference on Computational Mechanics Cracow, 26 - 29 June 2001
Preface Composite materials accompanied the human activity from the beginning of the civilisation.Apart from natural composites,like the wood,applied in various structures people invented many multi-component materials even in ancient times. One of the most famous applications of the old-time composites is the Chinese Wall,whose durability and stability was ensured by contrastively different materials incorporated into a single structure.Next applications worked out and popularised in Central Europe in the Middle Ages was known as the Prussian wall combining the wooden skeleton filled with the bricks.One of the most significant milestones in the history of modern composites was the application of the concrete reinforced with the steel bars in France at the end of the nineteenth century. Nowadays composites play a very important role in engineering from aerospace technology and nuclear devices to microelectronics or structural engineering applications [37,128,203,286,298,351,367,389].Considering this fact and the growing role of numerical experiments in the designing of structures and industrial processes,one of the most important purposes of computational mechanics research and direction of progress appeared to be precise numerical modelling of these materials.On the other hand,experimental sciences prove that every structural parameter has a random,in fact stochastic,character.Thus,many probabilistic approaches and methodologies have emerged recently to simulate more accurately the real behaviour of mechanical systems and processes.These methods show that the random character of parameters discussed is very important for the systems simulated [14,121,357].This conclusion may lead us to the hypothesis,that the random character of the material and physical parameters should play an essential role in multi-component structures [32,34,151,154,275]. Modern computational mechanics of composite materials follows many various ways through different science domains from experimental materials science to advanced computational techniques and applied mathematics.They engage more and more complicated and precise testing methods and devices,stochastic and sensitivity analysis algorithms and multiscale domain theoretical solutions for partial and ordinary differential equations reflecting some practical engineering and physical problems.Commercial computer programs based on the Finite Element Method enable now visualisation of the multifield,multiphase and non- stationary physical and mechanical problems and even introducing uncertainty into computer simulation using random variables(ANSYS,for instance).The growth of computer power obtained from technological progress and advances in parallel numerical techniques practically eliminated the parameter of the cost of computational time in modelling,which resulted in the efficient implementation and use of Monte Carlo simulation. The basic idea behind this book was to collect relatively up-to-date approaches to the composite materials lying somewhere in between experimental measurements and their opportunities,theoretical advances in applied mathematics
Preface Composite materials accompanied the human activity from the beginning of the civilisation. Apart from natural composites, like the wood, applied in various structures people invented many multi-component materials even in ancient times. One of the most famous applications of the old-time composites is the Chinese Wall, whose durability and stability was ensured by contrastively different materials incorporated into a single structure. Next applications worked out and popularised in Central Europe in the Middle Ages was known as the Prussian wall combining the wooden skeleton filled with the bricks. One of the most significant milestones in the history of modern composites was the application of the concrete reinforced with the steel bars in France at the end of the nineteenth century. Nowadays composites play a very important role in engineering from aerospace technology and nuclear devices to microelectronics or structural engineering applications [37,128,203,286,298,351,367,389]. Considering this fact and the growing role of numerical experiments in the designing of structures and industrial processes, one of the most important purposes of computational mechanics research and direction of progress appeared to be precise numerical modelling of these materials. On the other hand, experimental sciences prove that every structural parameter has a random, in fact stochastic, character. Thus, many probabilistic approaches and methodologies have emerged recently to simulate more accurately the real behaviour of mechanical systems and processes. These methods show that the random character of parameters discussed is very important for the systems simulated [14,121,357]. This conclusion may lead us to the hypothesis, that the random character of the material and physical parameters should play an essential role in multi-component structures [32,34,151,154,275]. Modern computational mechanics of composite materials follows many various ways through different science domains from experimental materials science to advanced computational techniques and applied mathematics. They engage more and more complicated and precise testing methods and devices, stochastic and sensitivity analysis algorithms and multiscale domain theoretical solutions for partial and ordinary differential equations reflecting some practical engineering and physical problems. Commercial computer programs based on the Finite Element Method enable now visualisation of the multifield, multiphase and nonstationary physical and mechanical problems and even introducing uncertainty into computer simulation using random variables (ANSYS, for instance). The growth of computer power obtained from technological progress and advances in parallel numerical techniques practically eliminated the parameter of the cost of computational time in modelling, which resulted in the efficient implementation and use of Monte Carlo simulation. The basic idea behind this book was to collect relatively up-to-date approaches to the composite materials lying somewhere in between experimental measurements and their opportunities, theoretical advances in applied mathematics
viii Computational Mechanics of Composite Materials and mechanics,numerical algorithms and computers as well as the practical needs of the engineers.The methods are well-documented in the context of computer batch files,scripts and computer programs.It will enable the readers to start from this point and to continue and/or replace the ideas with newer,more accurate and efficient ones.The author believes that this book will appear to be useful for applied mathematicians,specialists in numerical methods and for engineers:civil, mechanical,aerospace and from related branches of industry.Some elements of probabilistic calculus and computation as well as general ideas can also be applied by students,who can incorporate these concepts into new research or into the existing well-documented knowledge dealing with composite materials. A primary version of the book was completed in Texas,during the author's postdoctoral research at Rice University in Houston in the academic year 1999/2000 under auspices of Prof.P.D.Spanos.The author would like to appreciate the help of many people,whose valuable comments and the time spent enabled finishing of the book.Special thanks are directed to Prof.Michal Kleiber from the Institute of Fundamental Technological Research,Polish Academy of Science in Warsaw,who expressed many precious ideas during a common research in random composites and who promoted this research.Prof.Tran Duong Hien from Technical University of Szczecin influenced the work in the area of stochastic finite elements.The cooperation with Prof.B.A.Schrefler from the University of Padua in Italy concerning numerical analysis of superconducting composites remarkably enhanced the relevant computational illustration included in the book. The help of Mr.Lukasz Figiel,M.Sc.and Mr.Marcin Pawlik,Dr.Eng.,two of my younger colleagues,was decisive for finishing of some computations devoted to heat transfer and fracture analysis.The author would like to express his respect to all the colleagues from Chair of Mechanics of Materials at the Technical University of Lodz for their advising voices,too.Last but not least,the role of the unknown reviewers,the editors and the people who commented and criticised this work is also appreciated. Layout of the Book Mathematical preliminaries open the book considerations and consist of basic definitions of random events,variables and probabilistic moments as well as description of the Monte Carlo simulation technique with the relevant statistical estimation theory elements.The stochastic perturbation approach (second order second central moment generalised to thenth order and higher moments technique) is explained using two examples:a transient heat transfer equation and the solution of the linear elastodynamic problem.The solution to these problems in terms of expected values and standard deviations as well as spatial and temporal cross- covariances is demonstrated and it illustrates the applicability of the method.An important part of this opening chapter is a probabilistic algebraic description of some transforms of random variables,which is necessary for further formulation and development of the stochastic interface defects model.Some of them are valid
viii Computational Mechanics of Composite Materials and mechanics, numerical algorithms and computers as well as the practical needs of the engineers. The methods are well-documented in the context of computer batch files, scripts and computer programs. It will enable the readers to start from this point and to continue and/or replace the ideas with newer, more accurate and efficient ones. The author believes that this book will appear to be useful for applied mathematicians, specialists in numerical methods and for engineers: civil, mechanical, aerospace and from related branches of industry. Some elements of probabilistic calculus and computation as well as general ideas can also be applied by students, who can incorporate these concepts into new research or into the existing well-documented knowledge dealing with composite materials. A primary version of the book was completed in Texas, during the author’s postdoctoral research at Rice University in Houston in the academic year 1999/2000 under auspices of Prof. P.D. Spanos. The author would like to appreciate the help of many people, whose valuable comments and the time spent from the Institute of Fundamental Technological Research, Polish Academy of Science in Warsaw, who expressed many precious ideas during a common research in random composites and who promoted this research. Prof. Tran Duong Hien from Technical University of Szczecin influenced the work in the area of stochastic finite elements. The cooperation with Prof. B.A. Schrefler from the University of Padua in Italy concerning numerical analysis of superconducting composites remarkably enhanced the relevant computational illustration included in the book. younger colleagues, was decisive for finishing of some computations devoted to heat transfer and fracture analysis. The author would like to express his respect to all the colleagues from Chair of Mechanics of Materials at the Technical unknown reviewers, the editors and the people who commented and criticised this work is also appreciated. Layout of the Book Mathematical preliminaries open the book considerations and consist of basic definitions of random events, variables and probabilistic moments as well as description of the Monte Carlo simulation technique with the relevant statistical estimation theory elements. The stochastic perturbation approach (second order second central moment generalised to the nth order and higher moments technique) is explained using two examples: a transient heat transfer equation and the solution of the linear elastodynamic problem. The solution to these problems in terms of expected values and standard deviations as well as spatial and temporal crosscovariances is demonstrated and it illustrates the applicability of the method. An important part of this opening chapter is a probabilistic algebraic description of some transforms of random variables, which is necessary for further formulation and development of the stochastic interface defects model. Some of them are valid enabled finishing of the book. Special thanks are directed to Prof. Michał Kleiber The help of Mr. Łukasz Figiel, M.Sc. and Mr. Marcin Pawlik, Dr. Eng., two of my University of Ło´ dz´ for their advising voices, too. Last but not least, the role of the
Preface ix for the Gaussian variates only,which essentially bounds the application.However, it leads to the specific formulae implemented further in the computer software attached.An important issue raised in this chapter is to show a difference between Gaussian and quasi-Gaussian random variables defined on some unempty and bounded real subsets. Elastic problems related to deterministic and probabilistic systems are collected in Chapter 2.They are divided into two essentially different parts-the first shows the linear elastic behaviour of some composite materials and structures in boundary value problems connected with their real microstructure.The other part contains description of the homogenisation technique together with the relevant numerical tests documenting the computational determination of so-called homogenisation functions,a posteriori error analysis related to homogenisation problems, probabilistic moments of effective material tensors and their variability with respect to some input parameters. The first part of this chapter starts from the mathematical model of composite, whose material characteristics are given arbitrarily as constant deterministic values or by using the first two probabilistic moments constant through the given component material region (or volume).Further,the stochastic interface defects concept is presented,which originated from some computational contact mechanics models.The interface defects are introduced as semicircles lying on the interface into a weaker material.The radii and total number of these defects are input cut-off Gaussian random variables defined using their expected values and the variances (or standard deviations)with elastic properties equal to 0.The modeling is performed through the following steps:(i)determination of the interphase-a thin film containing all the defects with thickness determined from defect probabilistic parameters,(ii)probabilistic spatial averaging of the defects over the interphase area,(iii)computational analysis of a new composite with the new extra component.Obviously,it is not possible to approximate the real composite with stochastic interface microdefects very accurately.However it can be and it is done intermediately-by comparison with the composites with the weakened interphase or interface,for instance.Computational experiments validating the model are performed using the system ABAQUS [1](in the deterministic approach)and the specially adapted academic package POLSAP(for the Stochastic Finite Element Method -SFEM needs)[183].All the results obtained for various composites and various combinations of interface defect parameters demonstrate a high level of structural uncertainty in the case of their presence as well as a significant increase of the structural state functions stresses and displacements around the interface region.The second part of the chapter concerns the homogenisation method both in deterministic and probabilistic context.Computational experiments dealing with a numerical solution of the homogenisation problem are done thanks to the FEM commercial system ANSYS [2],where most of the databases for these experiments are available from the author to be used in further extensions of mathematical and mechanical homogenisation model
Preface ix for the Gaussian variates only, which essentially bounds the application. However, it leads to the specific formulae implemented further in the computer software attached. An important issue raised in this chapter is to show a difference between Gaussian and quasi-Gaussian random variables defined on some unempty and bounded real subsets. Elastic problems related to deterministic and probabilistic systems are collected in Chapter 2. They are divided into two essentially different parts – the first shows the linear elastic behaviour of some composite materials and structures in boundary value problems connected with their real microstructure. The other part contains description of the homogenisation technique together with the relevant numerical tests documenting the computational determination of so-called homogenisation functions, a posteriori error analysis related to homogenisation problems, probabilistic moments of effective material tensors and their variability with respect to some input parameters. The first part of this chapter starts from the mathematical model of composite, whose material characteristics are given arbitrarily as constant deterministic values or by using the first two probabilistic moments constant through the given component material region (or volume). Further, the stochastic interface defects concept is presented, which originated from some computational contact mechanics models. The interface defects are introduced as semicircles lying on the interface into a weaker material. The radii and total number of these defects are input cut-off Gaussian random variables defined using their expected values and the variances (or standard deviations) with elastic properties equal to 0. The modeling is performed through the following steps: (i) determination of the interphase – a thin film containing all the defects with thickness determined from defect probabilistic parameters, (ii) probabilistic spatial averaging of the defects over the interphase area, (iii) computational analysis of a new composite with the new extra component. Obviously, it is not possible to approximate the real composite with stochastic interface microdefects very accurately. However it can be and it is done intermediately – by comparison with the composites with the weakened interphase or interface, for instance. Computational experiments validating the model are performed using the system ABAQUS [1] (in the deterministic approach) and the specially adapted academic package POLSAP (for the Stochastic Finite Element Method – SFEM needs) [183]. All the results obtained for various composites and various combinations of interface defect parameters demonstrate a high level of structural uncertainty in the case of their presence as well as a significant increase of the structural state functions stresses and displacements around the interface region. The second part of the chapter concerns the homogenisation method both in deterministic and probabilistic context. Computational experiments dealing with a numerical solution of the homogenisation problem are done thanks to the FEM commercial system ANSYS [2], where most of the databases for these experiments are available from the author to be used in further extensions of mathematical and mechanical homogenisation model
x Computational Mechanics of Composite Materials Interface defects model and probabilistic homogenisation using both Monte Carlo simulation techniques are analysed using the authors FEM implementation called MCCEFF.The results of simulation are compared in terms of expected values and variances with analogous results obtained through the stochastic second order perturbation methodology.The appendix to this chapter consists of necessary fundamental mathematical theorems and definitions for the asymptotic homogenisation theorem. Elastoplasticity of composites discussed in the next chapter is focused on the alternative homogenisation technique,where instead of periodicity conditions imposed on the external boundaries of the RVE,some combination of the symmetry conditions and strain fields are applied to this element.The application of this method to the homogenisation of a periodic superconducting coil cable is also shown-an effective elastoplastic constitutive law is determined numerically and shown as a function of the homogenising uniform strain applied at the RVE boundary.Analogously to the methods typical for elastostatic problems,the closed- form equations for effective yield stresses are formulated in various ways,which can next be extended on probabilistic analysis.This chapter is completed with the transformation matrices algebraic definition,which is the essence of the computational implementation of the method.Probabilistic moments of the effective elastoplastic constitutive law can be obtained as a conjunction of this method with the Monte Carlo simulation technique discussed in the previous chapter.The fundamental issue is however experimental determination of higher order probabilistic moments for the superconductor material characteristics; otherwise the analysis is useful in the context of the sensitivity of the homogenised characteristics with respect to the adopted level of input randomness only. Sensitivity analysis presented in Chapter 4 is entirely devoted to a relatively new research area-determination of the sensitivity gradients for homogenised material characteristics.For this purpose two essentially different homogenisation methods are used-algebraic approximation and asymptotic methodology.Starting from a traditional description of the effective parameters in both methods,the sensitivity gradients are determined by the symbolic calculus approach and,on the other hand,pure computational strategy based on the Finite Difference Method (FDM).The implementation and results obtained from these two methods demonstrate the basic limitations of the methods,i.e.necessity of closed-form equations for the symbolic approach and numerical instabilities in the FDM simulations.This knowledge is necessary for significant time savings in the extension of this study to the random composite sensitivity analysis where the heterogeneous periodic composites with probabilistically defined material properties are analysed.The probabilistic sensitivity of such structures is defined through the introduction of sensitivity gradients of probabilistic moments of the effective material parameters with respect to the appropriate moments of composite structure parameters-elastic properties of the constituents as well as interface defect data. Fracture and fatigue -the collection of various fatigue theories with special emphasis placed on the second order perturbation method application are discussed
x Computational Mechanics of Composite Materials Interface defects model and probabilistic homogenisation using both Monte Carlo simulation techniques are analysed using the authors FEM implementation called MCCEFF. The results of simulation are compared in terms of expected values and variances with analogous results obtained through the stochastic second order perturbation methodology. The appendix to this chapter consists of necessary fundamental mathematical theorems and definitions for the asymptotic homogenisation theorem. Elastoplasticity of composites discussed in the next chapter is focused on the alternative homogenisation technique, where instead of periodicity conditions imposed on the external boundaries of the RVE, some combination of the symmetry conditions and strain fields are applied to this element. The application of this method to the homogenisation of a periodic superconducting coil cable is also shown – an effective elastoplastic constitutive law is determined numerically and shown as a function of the homogenising uniform strain applied at the RVE boundary. Analogously to the methods typical for elastostatic problems, the closedform equations for effective yield stresses are formulated in various ways, which can next be extended on probabilistic analysis. This chapter is completed with the transformation matrices algebraic definition, which is the essence of the computational implementation of the method. Probabilistic moments of the effective elastoplastic constitutive law can be obtained as a conjunction of this method with the Monte Carlo simulation technique discussed in the previous chapter. The fundamental issue is however experimental determination of higher order probabilistic moments for the superconductor material characteristics; otherwise the analysis is useful in the context of the sensitivity of the homogenised characteristics with respect to the adopted level of input randomness only. Sensitivity analysis presented in Chapter 4 is entirely devoted to a relatively new research area – determination of the sensitivity gradients for homogenised material characteristics. For this purpose two essentially different homogenisation methods are used – algebraic approximation and asymptotic methodology. Starting from a traditional description of the effective parameters in both methods, the sensitivity gradients are determined by the symbolic calculus approach and, on the other hand, pure computational strategy based on the Finite Difference Method (FDM). The implementation and results obtained from these two methods demonstrate the basic limitations of the methods, i.e. necessity of closed-form equations for the symbolic approach and numerical instabilities in the FDM simulations. This knowledge is necessary for significant time savings in the extension of this study to the random composite sensitivity analysis where the heterogeneous periodic composites with probabilistically defined material properties are analysed. The probabilistic sensitivity of such structures is defined through the introduction of sensitivity gradients of probabilistic moments of the effective material parameters with respect to the appropriate moments of composite structure parameters – elastic properties of the constituents as well as interface defect data. Fracture and fatigue – the collection of various fatigue theories with special emphasis placed on the second order perturbation method application are discussed
Preface xi next.The crucial numerical illustration is presented in the case of the Paris Erdogan rule where some of the system input data are treated as random variables. Therefore,expected values are compared against the deterministic values and standard deviations are added,too.An analogous approach is used to reformulate the well-known fracture criteria applied for composite materials-Tsai-Wu and Tsai-Hill-and to use them in symbolic computations for probabilistic parameters of the composite material fracture parameters.The essential part of this chapter is devoted to the FEM modelling of fracture and fatigue of some composites where analytical solutions are not available.Computational illustrations consist of static fracture of curved composite under shear loading leading to the delamination, fatigue analysis of composite pipe joint as well as thermomechanical fatigue of the curved laminate under thermal and/or static quasistatic load varying in time with constant amplitude.Most of the frequently used theories and equations for fatigue analysis are collected in the appendix to this chapter. Reliability analysis is included in the Chapter 6 and it consists of a discussion of various order reliability computational approaches together with the Weibull Second Order and Third Moment model (W-SOTM).This methodology is used to compute the reliability index for the composite Hertz contact problem,where elastic spherical inclusion of the reinforcement is loaded by the force to remain in contact with the matrix.Further in this chapter a stochastic process description of the degradation phenomena is also given,which appears to be common for the homogeneous and heterogeneous structures and materials.It can find a broad field of applications together with efficient implementations of stochastic processes (with both spatial and temporal randomness)in the Finite Element Method (or BEM,FDM,meshless as well as hybrid method based)programs. An application of the wavelet-based multiresolutional approach to composite materials in terms of homogenisation of multiscale media is the extension of previous considerations and concludes the book.The traditional composite materials model consisting of two or three geometrical scales is now rewritten in view of practically infinite number of separate scales (resolutions)that can be linked using interscale wavelet projection (some mathematical transformation). The basic tool necessary for such an analysis development is the basic wavelet basis(a mathematical function varying rapidly in a given geometrical scale),which can be used now to transform between neighboring scales.The homogenised characteristics for the composite can be determined usually in the closed-form equation if and only if the limit of an infinite series of wavelet projections between all geometrical scales exists and is unique.As is illustrated by some wavelet function samples,such an analysis type can be some alternative for the random analysis,because the wavelet functions used in various scale makes,in the coarsest scale,the impression that the relevant material property demonstrates the great level of some kind of uncertainty.It is not underlined clearly that the main limitation of this methodology is that the wavelet projection between the neighbouring scales can be continued through the range of validity of the same physical laws.It is not possible to carry out the passage from the atomistic to the global scale of the composite using the same wavelet projection and,most
Preface xi next. The crucial numerical illustration is presented in the case of the Paris Erdogan rule where some of the system input data are treated as random va riables. Therefore, expected values are compared against the deterministic values and standard deviations are added, too. An analogous approach is used to reformulate the well-known fracture criteria applied for composite materials – Tsai-Wu and Tsai-Hill - and to use them in symbolic computations for probabilistic parameters of the composite material fracture parameters. The essential part of this chapter is devoted to the FEM modelling of fracture and fatigue of some composites where analytical solutions are not available. Computational illustrations consist of static fracture of curved composite under shear loading leading to the delamination, fatigue analysis of composite pipe joint as well as thermomechanical fatigue of the curved laminate under thermal and/or static quasistatic load varying in time with constant amplitude. Most of the frequently used theories and equations for fatigue analysis are collected in the appendix to this chapter. Reliability analysis is included in the Chapter 6 and it consists of a discussion of various order reliability computational approaches together with the Weibull Second Order and Third Moment model (W-SOTM). This methodology is used to compute the reliability index for the composite Hertz contact problem, where elastic spherical inclusion of the reinforcement is loaded by the force to remain in contact with the matrix. Further in this chapter a stochastic process description of the degradation phenomena is also given, which appears to be common for the homogeneous and heterogeneous structures and materials. It can find a broad field of applications together with efficient implementations of stochastic processes (with both spatial and temporal randomness) in the Finite Element Method (or BEM, FDM, meshless as well as hybrid method based) programs. An application of the wavelet-based multiresolutional approach to composite materials in terms of homogenisation of multiscale media is the extension of previous considerations and concludes the book. The traditional composite materials model consisting of two or three geometrical scales is now rewritten in view of practically infinite number of separate scales (resolutions) that can be linked using interscale wavelet projection (some mathematical transformation). The basic tool necessary for such an analysis development is the basic wavelet basis (a mathematical function varying rapidly in a given geometrical scale), which can be used now to transform between neighboring scales. The homogenised characteristics for the composite can be determined usually in the closed-form equation if and only if the limit of an infinite series of wavelet projections between all geometrical scales exists and is unique. As is illustrated by some wavelet function samples, such an analysis type can be some alternative for the random analysis, because the wavelet functions used in various scale makes, in the coarsest scale, the impression that the relevant material property demonstrates the great level of some kind of uncertainty. It is not underlined clearly that the main limitation of this methodology is that the wavelet projection between the neighbouring scales can be continued through the range of validity of the same physical laws. It is not possible to carry out the passage from the atomistic to the global scale of the composite using the same wavelet projection and, most