Being a widely used structure, composite shells have been studied for a long time. The features of small thickness, heterogeneity, and anisotropy of composite shells have created many challenges for analyzing them. A number of theories have been developed for modeling composite shells, while they are either not practical for engineering use, or rely on assumptions that do not always hold. Consequently, a better theory is needed, especially for the application on challenging problems such as shells involving thermoelasticity, viscoelasticity, or viscoplasticity.<br><br>In this dissertation, a shell theory based on mechanics of structure genome (MSG), a unified theory for multiscale constitutive modeling, is developed. This theory is capable of handling fully anisotropy and complex heterogeneity, and because the derivation follows principle of minimum information loss (PMIL) and using the variational asymptotic method (VAM), high accuracy can be achieved. Both a linear version and a nonlinear version using Euler method combined with Newton-Raphson method are presented. This MSG-based shell theory is used for analyzing the curing process of composites, deployable structures made with thin-ply high strain composite (TP-HSC), and material nonlinear shell behaviors.<br><br>When using the MSG-based shell theory to simulate the curing process of composites, the formulation is written in an analytical form, with the effect of temperature change and degree of cure (DOC) included. In addition to an equivalent classical shell theory, a higher order model with the correction from initial geometry and transverse shear deformation is presented in the form of the Reissner-Mindlin model. Examples show that MSG-based shell theory can accurately capture the deformation caused by temperature change and cure shrinkage, while errors exist when recovering three-dimensional (3D) strain field. Besides, the influence of varying transverse shear stiffness needs to be further studied.<br><br>In order to analyze TP-HSC deployable structures, linear viscoelasticity behavior of composite shells is modeled. Then, column bending test (CBT), an experiment for testing the bending stiffness of thin panels under large bending deformation, is simulated with both quasi-elastic (QE) and direct integration (DI) implementation of viscoelastic shell properties. Comparisons of the test and analysis results show that the model is capable of predicting most of the measured trends. Residual curvature measured in the tests, but not predicted by the present model, suggests that viscoplasticity should be considered. A demonstrative study also shows the potential of material model calibration using the virtual CBT developed in this work. A deployable boom structure is also analyzed. The complete process of flattening, coiling, stowage, deployment and recovery is simulated with the viscoelastic shell model. Results show that major residual deformation happens in the hoop direction.<br><br>A nonlinear version of the MSG-based general purpose constitutive modeling code SwiftComp is developed. The nonlinear solving algorithm based on the combined Euler-Newton method is implemented into SwiftComp. For the convenience of implementing a nonlinear material model, the capability of using user material is also added. A viscoelastic material model and a continuum damage model is tested and shows excellent match when compared with Abaqus results with solid elements and UMAT. Further validation of the nonlinear SwiftComp is done with a nonlinear viscoelastic-viscoplastic model. The high computational cost is emphasized with a preliminary study with surrogate model.

<div>Polymer matrices present in composite materials are prone to have time-dependent behavior very sensitive to changes in temperature. The modeling of thermoviscoelasticity is fundamental for capturing the performance of anisotropic viscoelastic materials subjected to both mechanical and thermal loads, or for manufacturing simulation of composites. In addition, improved plate/shell and beam models are required to efficiently design and simulate large anisotropic composite structures. Numerical models have been extensively used to capture the linear viscoelasticity in composites, which can be generalized in integral or differential forms. The hereditary integral constitutive form has been adopted by many researchers to be implemented into finite element codes, but its formulation is complex and time consuming as it is function of the time history. The differential formulation provides faster computation times, but its applicability has been limited to capture the behavior of three-dimensional thermoviscoelastic orthotropic materials.</div><div><br></div><div>This work extends mechanics of structure genome (MSG) to construct linear thermoviscoelastic solid, plate/shell and beam models for multiscale constitutive modeling of three-dimensional heterogeneous materials made of time and temperature dependent constituents. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The thermal expansion creep is treated as inherent material behavior. Exact three-dimensional thermoviscoelastic homogenization solutions are also formulated for laminates modeled as an equivalent, homogeneous, anisotropic solid. The new model is implemented in SwiftComp, a general-purpose multiscale constitutive modeling code based on MSG, to handle real heterogeneous materials with arbitrary microstructures, mesostructures or cross-sectional shapes.</div><div><br></div><div>Three-dimensional representative volume element (RVE) analyses and direct numerical simulations using a commercial finite element software are conducted to verify the accuracy of the MSG-based constitutive modeling. Additionally, MSG-based plate/shell results are validated against thin-ply high-strain composites experimental data showing good agreement. Numerical cases with uniform and nonuniform cross-sectional temperature distributions are studied. The results showed that unlike MSG, the RVE method exhibits limitations to properly capture the long-term behavior of effective coefficients of thermal expansion (CTEs) when time-dependent constituent CTEs are considered. The analyses of the homogenized properties also revealed that despite the heterogeneous nature of the composite material, from a multiscale analysis perspective, the temperature dependencies of the effective stiffness and thermal stress properties are governed by the same shift factor as the polymer matrix. This conclusion remains the same for MSG-based solid, plate/shell and beam models with uniform temperature distributions.</div>

<div>Textile composites have been widely used due to the excellent mechanical performance and lower manufacturing costs, but the accurate prediction of the mechanical behaviors of textile composites is still very challenging due to the complexity of the microstructures and boundary conditions. Moreover, there is an unprecedented amount of design options of different textile composites. Therefore, a highly efficient yet accurate approach, which can predict the macroscopic structural performance considering different geometries and materials at subscales, is urgently needed for the structural design using textile composites.</div><div><br></div><div>Mechanics of structure genome (MSG) is used to perform multiscale modeling to predict various performances of textile composite materials and structures. A two-step approach is proposed based on the MSG solid model to compute the elastic properties of different two-dimensional (2D) and three-dimensional (3D) woven composites. The first step computes the effective properties of yarns at the microscale based on the fiber and matric properties. The effective properties of yarns and matrix are then used at the mesoscale to compute the properties of woven composites in the second step. The MSG plate and beam models are applied to thin and slender textile composites, which predict both the structural responses and local stress field. In addition, the MSG theory is extended to consider the pointwise temperature loads by modifying the variational statement of the Helmholtz free energy. Instead of using coefficients of thermal expansions (CTEs), the plate and beam thermal stress resultants derived from the MSG plate and beam models are used to capture the thermal-induced behaviors in thin and slender textile composite structures. Moreover, the MSG theory is developed to consider the viscoelastic behaviors of textile composites based on the quasi-elastic approach. Furthermore, a meso-micro scale coupled model is proposed to study the initial failure of textile composites based on the MSG models which avoids assuming a specific failure criterion for yarns. The MSG plate model uses plate stress resultants to describe the initial failure strength that can capture the stress gradient along the thickness in the thin-ply textile composites. The above developments of MSG theory are validated using high-fidelity 3D finite element analysis (FEA) or experimental data. The results show that MSG achieves the same accuracy of 3D FEA with a significantly improved efficiency.</div><div> </div><div>Taking advantage of the advanced machine learning model, a new yarn failure criterion is constructed based on a deep neural network (DNN) model. A series of microscale failure analysis based on the MSG solid model is performed to provide the training data for the DNN model. The DNN-based failure criterion as well as other traditional failure criteria are used in the mesoscale initial failure analysis of a plain woven composite. The results show that the DNN yarn failure criterion gives a better accuracy than the traditional failure criteria. In addition, the trained model can be used to perform other computational expensive simulations such as predicting the failure envelopes and the progressive failure analysis.</div><div> </div><div>Multiple software packages (i.e., texgen4sc and MSC.Patran/Nastran-SwiftComp GUI) are developed to incorporate the above developments of the MSG models. These software tools can be freely access and download through cdmHUB.org, which provide practical tools to facilitate the design and analysis of textile composite materials and structures.</div>

Quick and accurate predictions of equivalent properties for thin-walled composite structures are required in the preliminary design process. Existing literature provides analytical solutions to some structures but is limited to particular cases. No unified approach exists to tackle homogenization of thin-walled structures such as beams, plates, or three-dimensional structures using the thin-walled approximation. In this work, a unified approach is proposed to obtain equivalent properties for beams, plates, and three-dimensional structures for thin-walled composite structures using mechanics of structure genome. The adopted homogenization technique interprets the unit cell associated with the composite structures as an assembly of plates, and the overall strain energy density of the unit cell as a summation of the plate strain energies of these individual plates. The variational asymptotic method is then applied to drop all higher-order terms and the remaining energy is minimized with respect to the unknown fluctuating functions. This has been done by discretizing the two-dimensional unit cell into one-dimensional frame elements in a finite element description. This allows the handling of structures with different levels of complexities and internal geometry within a general framework. Comparisons have been made with other works to show the advantages which the proposed model offers over other methods.

<p>The Mechanics of Structural Genome (MSG) is a unified homogenization theory used to find equivalent constitutive models for beam, plate, and solid structures. It has been proven accurate for periodic structures. However, for certain applications such as non-prismatic wind turbine blades and helicopter flexbeams featuring ply drop-off, where there is no repeating structure and the periodic boundary condition cannot be used, MSG's accuracy is limited. In this work, we aim to extend MSG to find element stiffness matrices directly for aperiodic structures, instead of beam properties or three-dimensional (3D) solid material properties. Two finite elements based on MSG have been developed: Heterogeneous Beam Element (HBE) and Heterogeneous Solid Element (HSE).</p> <p><br></p> <p>For beam modeling, the beam-like structure is homogenized into a series of 3-node Heterogeneous Beam Elements (HBE) with 18×18 effective beam element stiffness matrices. These matrices are used as input for one-dimensional (1D) beam analysis using the Abaqus User Element subroutine (UEL). Using the macroscopic beam analysis results as input, we can also perform dehomogenization to predict the stresses and strains in the original structure. We use three examples (a prismatic composite beam, an isotropic homogeneous tapered beam, and a composite tapered beam) to demonstrate the capability of HBE and show its advantages over the MSG cross-sectional analysis approach. HBE can capture macroscopic behavior and detailed stresses due to non-prismatic geometry.</p> <p><br></p> <p>The Heterogeneous Solid Element (HSE) is developed based on MSG to model a heterogeneous body as an equivalent solid element using an effective element stiffness matrix. HSE modeling includes homogenization, macroscopic global analysis, and dehomogenization to recover local strains/stresses. HSE avoids the local periodicity assumption for traditional multiscale modeling techniques for composite structures that compute effective material properties instead. Abaqus composite solid element and MSG-based traditional multiscale modeling are used to validate the accuracy of HSE. All example results show that HSE is more accurate in predicting global structural behavior and local strains/stresses.</p> <p><br></p> <p>HBE and HSE provide a new concept for modeling aperiodic composite structures by modeling structures into equivalent beam or solid elements instead of beam properties of the reference line in 1D beam analysis or material properties of material points in solid structural analysis.</p>

Among the many active materials in use today, piezoelectric composite patches have enabled notable advances in emerging technologies such as disturbance sensing, control of flexible structures, and energy harvesting. The macro fiber composite (MFC), in particular, is well known for its outstanding performance. Multiscale models are typically required for smart-structure design with MFCs. This is due to the need for predicting the macroscopic response (such as tip deflection under a transverse load or applied voltage) while accounting for the fact that the MFC has microscale details. Current multiscale models of the MFC exclusively focus on predicting the macroscopic response with homogenized material properties. There are a limited number of homogenized properties available from physical experiments and various aspects of existing homogenization techniques for the MFC are shown here to be inadequate. Thus, new homogenized models of the MFC are proposed to improve smart-structure predictions and therefore improve device design. It is notable that current multiscale modeling efforts for MFCs are incomplete since, after homogenization, the local fields such as stresses and electric fields have not been recovered. Existing methods for obtaining local fields are not applicable since the electrodes of the MFC are embedded among passive layers. Therefore, another objective of this work was to find the local fields of the MFC without having the computational burden of fully modeling the microscopic features of the MFC over a macroscale area. This should enable smart-structure designs with improved reliability because failure studies of MFCs will be enabled. Large-scale 3D finite element (FE) models that included microscale features were constructed throughout this work to verify the multiscale methodologies. Note that after creating a free account on cdmhub.org, many files used to create the results in this work can be downloaded from https://cdmhub.org/projects/ernestocamarena.<br><br>First, the Mechanics of Structure Genome (MSG) was extended to provide a rigorous analytical homogenization method. The MFC was idealized to consist of a stack of homogeneous layers where some of the layers were homogenized with existing rules of mixtures. For the analytical model, the electrical behavior caused by the interdigitated electrodes (IDEs) was approximated with uniform poling and uniform electrodes. All other assumptions on the field variables were avoided; thus an exact solution for a stack of homogeneous layers was found with MSG. In doing so, it was proved that in any such multi-layered composite, the in-plane strains and the transverse stresses are equal in each layer and the in-plane electric fields and transverse electric displacement are constant between the electrodes. Using this knowledge, a hybrid rule of mixtures was developed to homogenize the entire MFC layup so as to obtain the complete set of effective device properties. Since various assumptions were avoided and since the property set is now complete, it is expected that greater energy equivalence between reality and the homogenized model has been made possible. The derivation clarified what the electrical behavior of a homogenized solid with internal electrodes should be—an issue that has not been well understood. The behavior was verified by large-scale FE models of an isolated MFC patch.<br> <br>Increased geometrical fidelity for homogenization was achieved with an FE-based RVE analysis that accounted for finite-thickness effects. The presented theory also rectifies numerous issues in the literature with the use of the periodic boundary conditions. The procedure was first developed without regard to the internal electrodes (ie a homogenization of the active layer). At this level, the boundary conditions were shown to satisfy a piezoelectric macrohomogeneity condition. The methodology was then applied to the full MFC layup, and modifications were implemented so that both types of MFC electrodes would be accounted for. The IDE case considered nonuniform poling and electric fields, but fully poled material was assumed. The inherent challenges associated with these nonuniformities are explored, and a solution is proposed. Based on the homogenization boundary conditions, a dehomogenization procedure was proposed that enables the recovery of local fields. The RVE analysis results for the effective properties revealed that the homogenization procedure yields an unsymmetric constitutive relation; which suggests that the MFC cannot be homogenized as rigorously as expected. Nonetheless, the obtained properties were verified to yield favorable results when compared to a large-scale 3D FE model.<br> <br>As a final test of the obtained effective properties, large-scale 3D FE models of MFCs acting in a static unimorph configuration were considered. The most critical case to test was the smallest MFC available. Since none of the homogenized models account for the passive MFC regions that surround the piezoelectric fiber array, some of the test models were constructed with and without the passive regions. Studying the deflection of the host substrate revealed that ignoring the passive area in smaller MFCs can overpredict the response by up to 20%. Satisfactory agreement between the homogenized models and a direct numerical simulation were obtained with a larger MFC (about a 5% difference for the tip deflection). Furthermore, the uniform polarization assumption (in the analytical model) for the IDE case was found to be inadequate. Lastly, the recovery of the local fields was found to need improvement.<br><br><br>

First, the theory of MSG is extended to aperiodic heterogeneous solid structures. Integral constraints are introduced to decompose the displacements and strains of the heterogeneous material into a fluctuating part and a macroscopic part, of which the macroscopic part represents the responses of the homogenized material. One advantage of this theory is that boundary conditions are not required. Consequently, it is capable of handling micro-structures of arbitrary shapes. In addition, periodic constraints can be incorporated into this theory as needed to model periodic or partially periodic materials such as textile composites. In this study, the newly developed method is employed to investigate the finite thickness effect of textile composites.<div><br></div><div>Second, MSG is enabled to deal with Timoshenko beam-like structures with spanwise heterogeneity, which provide higher accuracy than the previous available Euler–Bernoulli beam model. Its reduced form, the MSG beam cross sectional analysis, is found to be able to analyze generalized free-edge problems with arbitrary layups and subjected to general loads. In this method, the only assumption applied is that the laminate is long enough so that the Saint-Venant principle can be adopted. There is no limitation on the cross section of the laminate since no ad hoc assumption is involved with the microstructure geometry. This method solve the free-edge problem from a multiscale simulation point of view.<br></div><div><br></div>

Mechanics of structure genome (MSG) is a unified homogenization theory that provides constitutive modeling of three-dimensional (3D) continua, beams and plates. In present work, the author extends the MSG to study the buckling of structures such as stiffened and sandwich panels. Such structures are usually slender or flat and easily buckle under compressive loads or bending moments which may result in catastrophic failure.<div><br><div>Buckling studies of stiffened and sandwich panels are found to be scattered. Most of the existed theories employ unnecessary assumptions or only apply to certain types of structures. There are few unified approaches that are capable of studying the buckling of different kinds of structures altogether. The main improvements of current approach compared with other methods in the literature are avoiding unnecessary assumptions, the capability of predicting all possible buckling modes including the global and local buckling modes, and the potential in studying the buckling of various types of structures.<br></div><div><br></div><div>For global buckling that features small local rotations, MSG mathematically decouples the 3D geometrical nonlinear problem into a linear constitutive modeling using structure genome (SG) and a geometrical nonlinear problem defined in a macroscopic structure. As a result, the original structures are simplified as macroscopic structures such as beams, plates or continua with effective properties, and the global buckling modes are predicted on macroscopic structures. For local buckling that features finite local rotations, Green strain is introduced into the MSG theory to achieve geometrically nonlinear constitutive modeling. Newton’s method is used to solve the nonlinear equilibrium equations for fluctuating functions. To find the bifurcated fluctuating functions, the fluctuating functions are then perturbed under the Bloch-periodic boundary conditions. The bifurcation is found when the tangent stiffness associated with the perturbed fluctuating functions becomes singular. Moreover, the arc-length method is introduced to solve the nonlinear equilibrium equations for post-local-buckling predictions because of its robustness. The imperfection is included in the form of geometrical imperfection by superimposing the scaled buckling modes in linear perturbation analysis on mesh.<br></div><div><br></div><div>Extensive validation case studies are carried out to assess the accuracy of the MSG theory in global buckling analysis and post-global-buckling analysis, and assess the accuracy of the extended MSG theory in local buckling and post-local-buckling analysis. Results using MSG theory and extended MSG theory in buckling analysis are compared with direct numerical solutions such as 3D FEA results and results in literature. Parametric studies are performed to reveal the relative influence of selective geometric parameters on buckling behaviors. The extended MSG theory is also compared with representative volume element (RVE) analysis with Bloch-periodic boundary conditions using commercial finite element packages such as Abaqus to assess the efficiency and accuracy of the present approach.<br></div></div>

Face aux défis environnementaux actuels, les industriels ont mis en œuvre de nouveaux matériaux recyclables et permettant une réduction significative de la masse. Le développement de la résine thermoplastique Elium par ARKEMA s’inscrit dans cette problématique. L’utilisation de cette résine pour la fabrication de pièces composites qui peuvent être sujettes à des dommages d’impact, nécessite au préalable des études, dans le but de comprendre leurs mécanismes de ruine sous ce type de sollicitation. Ainsi, la présente thèse propose une contribution à l’analyse multi-échelle de la tenue à l’impact des composites stratifiés à base de la résine Elium. Une étude expérimentale préliminaire a permis de confirmer la meilleure résistance à l’impact des composites à matrice Elium acrylique, comparativement à celles des composites thermodurcissables conventionnels. Ensuite, les performances à l’impact des composites stratifiés ont été améliorées par l’introduction de copolymères à blocs dans la matrice. Ces derniers sont capables de former des micelles de tailles nanométriques et ainsi d’améliorer la ténacité de la matrice acrylique. Les effets de l’énergie d’impact, de la température et de la composition en nanocharges sur la réponse du matériau composite ont été analysés. Afin de proposer un outil d’aide à la prédiction de la réponse à l’impact des matériaux fibres de verre/Acrylique, deux stratégies de modélisation ont été retenues. La première modélisation (macroscopique) considère le pli tissé du stratifié comme un matériau homogène tandis que la seconde (mésoscopique) utilise une description géométrique de l’ondulation et de l’entrecroisement des torons noyés dans la résine Elium. Ces deux modèles considèrent des zones cohésives à l’interface entre les plis adjacents pour simuler le délaminage interlaminaire. Des essais de délaminage (expérimentaux et numériques) ont permis d’alimenter le modèle d’endommagement de l’interface interplis. D’autre part, des essais de caractérisation du comportement mécanique et de l’endommagement du matériau couplés à l’homogénéisation multi-échelle des matériaux par la Mécanique du Génome de Structure ont permis d’identifier les paramètres du modèle macroscopique. A l’échelle mésoscopique, le modèle géométrique a été réalisé grâce au logiciel Texgen. Ce logiciel permet d’obtenir une description approchée mais réaliste de l’ondulation des torons de fibres. La même description a servi à l’homogénéisation numérique multi-échelle des stratifiés étudiés. La simulation numérique de l’impact basse vitesse a été effectuée au moyen du logiciel d’éléments finis ABAQUS/Explicit. Les modèles de comportement du matériau ont été implémentés via la routine utilisateur VUMAT. Les résultats obtenus offrent une bonne corrélation avec les données expérimentales / In the race for light materials able of meeting modern environmental challenges, an acrylic resin (Elium) has been developed. Elium is a thermoplastic resin able to replace thermosetting matrices, which are widespread nowadays in the industrial world. The present study aims to evaluate the impact resistance and to understand the failure mechanisms of composite laminates based on acrylic matrix under impact loading. We provide a contribution to the multiscale analysis of the impact resistance of laminated composite.First, the impact resistance and the damage tolerance of the acrylic resin based composites were compared with those of conventional composites. Then, the impact performance of the laminated composites has been enhanced by adding copolymer blocks to the liquid acrylic resin. These copolymers are able to form micelles of nanometer sizes, which lead to the improvement of both the acrylic matrix fracture toughness and the impact resistance. The effects of the impact energy, temperature, and composition in nano-copolymers have also been investigated.In order to provide a numerical tool for the prediction of the impact response of the glass fiber/Acrylic laminates, two strategies have been analyzed. The first one, performed at the macroscopic scale, considers the woven ply of the laminate as homogeneous material, and the second one (at the mesoscopic scale), deals with a realistic geometrical description of the yarns undulation. Both models use cohesive zones at the interface between the adjacent plies, to simulate the delamination. For this purpose, experimental and numerical delamination tests were performed to feed the inter-ply damage model. Mechanical tests for material characterization were also performed on specimens in order to identify the ply-damage model parameters. The Mechanics of Structure Genome (MSG) and a finite element based micromechanics approaches were then conducted to evaluate the effective thermomechanical properties of the yarns and the plain woven composite laminate. The realistic topological and morphological textures of the composite were accounted through Texgen software. These numerical impact simulations were performed using the finite element software ABAQUS/Explicit. Both models were implemented through a user material subroutine VUMAT. The obtained results appear in a good agreement with the experimental data and confirm the relevance of the proposed approach.

<p>This study is concerned with the predictive modeling of the machining and the mechanical behaviors of additive manufactured (AMed) Ti6AlV/TiC composites and Ti6Al4V, respectively, using microstructure-based hierarchical multiscale modeling. The predicted results could constitute as a basis for optimizing the parameters of machining and AM of the current materials.</p> <p>Through hierarchical flow of material behaviors from the atomistic, to the microscopic and the macroscopic scales, multiscale heterogeneous models (MHMs) coupled to the finite element method (FEM) are employed to simulate the conventional and the laser assisted machining (LAM) of Ti6AlV/TiC composites. In the atomistic level, molecular dynamics (MD) simulations are used to determine the traction-separation relationship for the cohesive zone model (CZM) describing the Ti6AlV/TiC interface. Bridging the microstructures across the scales in MHMs is achieved by representing the workpiece by macroscopic model with the microscopic heterogeneous structure including the Ti6Al4V matrix, the TiC particles, and their interfaces represented by the parameterized CZM. As a result, MHMs are capable of revealing the possible reasons of the peculiar high thrust forces behavior during conventional machining of Ti6Al4V/TiC composites, and how laser assisted machining can improve this behavior, which has not been conducted before.</p> <p>Extending MHMs to predict the mechanical behaviors of AMed Ti6Al4V would require including the heterogeneous microstructure at the grain level, which could be computational expensive. To solve this issue, the extended mechanics of structure genome (XMSG) is introduced as a novel multiscale homogenization approach to predict the mechanical behavior of AMed Ti6Al4V in a computationally efficient manner. This is realized by embedding the effects of microstructure heterogeneity, porosity growth, and crack propagation in the multiscale calculations of the mechanical behavior of the AMed Ti6Al4V using FEM. In addition, the XMSG can predict the asymmetry in the Young’s modulus of the AMed Ti6Al4V under tensile and compression loading as well as the anisotropy in the mechanical behaviors. The applicability of XMSG to fatigue life prediction with valid results is conducted by including the energy dissipations associated with cyclic loading/unloading in the calculations of the cyclic response of the material.</p>