The use of energy methods and variational principles is widespread in many fields of engineering of which structural mechanics and curve and surface design are two prominent examples. In principle many different types of function can be used as possible trial solutions to a given variational problem but where piecewise polynomial behaviour and user controlled cross segment continuity is either required or desirable, B-splines serve as a natural choice. Although there are many examples of the use of B-splines in such situations there is no common thread running through existing formulations that generalises from the one dimensional case through to two and three dimensions. We develop a unified approach to the representation of the minimisation equations for B-spline based functionals in tensor product form and apply these results to solving specific problems in geometric smoothing and finite element analysis using the Rayleigh-Ritz method. We focus on the development of algorithms for the exact computation of the minimisation matrices generated by finding stationary values of functionals involving integrals of squares and products of derivatives, and then use these to seek new variational based solutions to problems in the above fields. By using tensor notation we are able to generalise the methods and the algorithms from curves through to surfaces and volumes. The algorithms developed can be applied to other fields where a variational form of the problem exists and where such tensor product B-spline functions can be specified as potential solutions.

Στα πλαίσια της παρούσας διατριβής προσεγγίστηκε το πρόβλημα της κρούσης των σύνθετων υλικών μέσω μιας συνολικής μεθοδολογίας. Ο όρος συνολική μπορεί να χαρακτηρίσει την παρούσα εργασία της διατριβής γιατί η αναλυτική και αριθμητική υλοποίηση του προβλήματος που παρουσιάσθηκε επιλύει πολλά παράπλευρα προβλήματα. Η επίλυση του προβλήματος της κρούσης χαμηλής ταχύτητας σε σύνθετα υλικά κινήθηκε σε τέσσερις φάσεις. Σε πρώτη φάση αναπτύχθηκε το κύριο εργαλείο με το οποίο εξελίσσεται η διαδικασία μοντελοποίησης του προβλήματος. Η κατάστρωση του μοντέλου επίλυσης φυσικών συχνοτήτων και η αριθμητική επίλυση αυτού δίνει εκτιμήσεις για την ακρίβεια της μεθόδου Ritz με τη χρήση των αποδεκτών συναρτήσεων p-Ritz. Η μέθοδος Ritz επιτρέπει την πλήρη μοντελοποίηση του πολύστρωτου λαμβάνοντας υπόψη και τους καμπτικούς όρους σύζευξης κάτι που μέχρι στιγμής δεν έχει παρουσιαστεί στη βιβλιογραφία. Κατά τη δεύτερη φάση και μετά την τεκμηρίωση της πρώτης φάσης και τον έλεγχο της ακρίβειας των υπολογισμών των φυσικών συχνοτήτων, διατυπώθηκε το δυναμικό πρόβλημα. Η ακρίβεια των υπολογισμών ήταν εξαιρετική και σε ικανοποιητική σύγκλιση με τα αποτελέσματα του κώδικα πεπερασμένων στοιχείων LSDYNA3D. Παράλληλα αποδείχτηκε ότι το δυναμικό μοντέλο, επιλύει το δυναμικό φαινόμενο για όλα τα σχήματα και τους συνδυασμούς συνοριακών συνθηκών. Σε εξέλιξη των ανωτέρω έγινε η διατύπωση του συζευγμένου προβλήματος της κρούσης και επιλύθηκε για διάφορα υλικά και συνοριακές συνθήκες που αποτελούν μοντέλα υποδείγματα στη διεθνή βιβλιογραφία. Η ικανότητα πρόβλεψης των μοντέλων χαρακτηρίζεται ιδιαίτερα ικανοποιητική ενώ παρουσιάστηκε ο τρόπος εισαγωγής μη γραμμικών νόμων επαφής, στο συζευγμένο σύστημα διαφορικών εξισώσεων. Με σκοπό τον πλήρη χαρακτηρισμό του φαινομένου αλλά και τη μείωση των παραμέτρων που χαρακτηρίζουν το πρόβλημα σε μια προσπάθεια χαρακτηρισμού των κύριων παραμέτρων που διέπουν την κρούση παρουσιάζεται ο τρόπος υπολογισμού των αδιάστατων ομάδων που περιγράφουν την κρούση. Στα πλαίσια αυτής της ανάλυσης υπολογίστηκε το Νομογράφημα Κρουστικής Απόκρισης. Το Νομογράφημα Κρουστικής Απόκρισης μπορεί να χρησιμοποιηθεί στην οργάνωση πειραμάτων χαρακτηρισμού της κρουστικής αντοχής πολύστρωτων διατάξεων. Για την επίλυση του κύριου προβλήματος προτάθηκε μια διαδικασία ιεραρχικής αντιμετώπισης του προβλήματος. Για να καθοριστεί η κρουστική απόκριση του πολύστρωτου η διάταξη υποβάλλεται σε δυναμικό και στατικό χαρακτηρισμό με τη χρήση αριθμητικών μοντέλων που βασίζονται στη μέθοδο Ritz. Υλοποιήθηκε μια μέθοδος που δύναται να επιλύει διάφορες γεωμετρίες σε ποικιλία συνοριακών συνθηκών. Στη βιβλιογραφία μέχρι στιγμής έχουν παρουσιαστεί μοντέλα για παραλληλόγραμμες πλάκες, ειδικώς ορθότροπες και σε συνοριακές συνθήκες απλής έδρασης. Υλοποιήθηκαν προγράμματα σε Fortran που επιλύουν στατικά και δυναμικά προβλήματα σε καμπτικά φορτία και για ποικίλους συνδυασμούς συνοριακών συνθηκών. Υλοποιήθηκε το συζευγμένο πρόβλημα της κρούσης πολύστρωτης διάταξης από σφαίρα (διεισδυτή) σε χαμηλές ταχύτητες Έγινε σύγκριση ελαστοπλαστικών μη γραμμικών μοντέλων επαφής με γραμμικά μοντέλα επαφής. Με τη χρήση της θεωρίας αδιαστατοποίησης καθορίστηκαν οι αδιάστατες ομάδες που διέπουν το πρόβλημα της κρούσης και δημιουργήθηκε το Νομογράφημα Κρουστικής Απόκρισης. Αποδείχτηκε ότι το Νομογράφημα Κρουστικής Απόκρισης ισχύει για οιαδήποτε συμμετρική πολύστρωτη διάταξη σε όποια συνοριακή συνθήκη και αν είναι και οιοδήποτε σχήμα πλάκας. Αποδείχτηκε ότι το Νομογράφημα Κρουστικής Απόκρισης μπορεί να χρησιμοποιηθεί ως χάρτης αντοχής και ανοχής σε κρούση με τη χρήση της αδιάστατου Κατωφλιού Δύναμης για Διαστρωματική Αποκόλληση. / Low-velocity impact in composite laminates is a destructive loading condition since it leads to significant internal damage, no detectable by visual inspection. Numerous researchers have made significant efforts to model the impact response of composite plates and facilitate the analysis and design for impact resistance, using simple linear spring-mass models or a combination of continuum mechanics models and contact laws. Many investigators have used simple engineering structures beams, plates and shells to demonstrate the impact response of composite structures. The basic concept for modelling the impact of an object to a target was proposed by Timoshenko (1913). He introduced the procedure where the description of the interaction between the impactor and the structure (Euler beam) was implemented using the Hertzian contact law. This approach was further extended in isotropic plates and shells. Lee, 1940; Greszcuk, 1982; Lee et al., 1983 and Shivakumar et al, 1985, proposed simple models of the low-velocity impact problem. Although many important contributions exist in the bibliography for the characterization of the impact response of laminated plates, analytical solutions are considered to be very few. Most of the models proposed are for analyzing specially orthotropic plates subjected to a local dominated impact. In addition, these models do not take into account the shear coupling terms of the bending stiffness matrix, e.g. D16 and D26. Sun and Chattopadhyay (1975), Dobyns (1981), and Ramkumar and Chen (1983) employed the first order shear deformation theory developed by Whitney and Pagano (1970), and used it in conjunction with the Hertzian contact law or plasticity contact law for characterizing the impact of laminated composite plates. In their analysis they have studied the impact response of a simply supported orthotropic plate subjected to central impact using the lamination theory that includes transverse shear deformations. It is evident that the contact force history must be computed as part of the solution of the dynamic response problem solving the nonlinear integral equation. Christoforou and Swanson (1991) and Carvalho (1996), obtained an analytical solution of the impact problem using the Laplace transform technique. Qian and Swanson (1990) obtained analytical solutions by linearising the contact deformation law and compared this with a Rayleigh-Ritz approach with numerical integration in time. The three dimensional finite element method with an explicit formulation where also used in the dynamic analysis of laminated plates subjected to low velocity impact. Moreover, many researchers have proposed simple models for characterizing the impact phenomenon. Shivakumar et.al. (1985) developed a simple model to predict the impact force and the duration of the impact phenomenon on composite plates. The composite plates were modeled as three springs while their stiffness was calculated by the plate properties and the contact parameters. Other researchers instead of a classical Hertzian law used the statical indentation law presented by Yang and Sun (1982) and Tan and Sun (1985). The non Hertzian Contact impact was considered for combining the overall deformation of the structure to the local deformation in the contact area and was used to predict the impact response of transversely isotropic beams and plates. In general, it is impossible to obtain the exact solution for an impact problem except by ignoring the local contact deformation or pre-assuming the contact force. The first step for understanding the problem is to predict the force applied by the projectile on the structure during impact. In order to characterize the contact force history, the model should account for the motion of the target, the effective stiffness of the projectile The present work deals with the development of a numerical scheme for the calculation of the dynamic response of any type of laminated composite plates subjected to any type (and/or combination) of boundary conditions under low-velocity impact. Using the p-version Ritz polynomials any shape that can be represented in Cartesian coordinates can be formulated. The governing second order differential equations are derived and allow furnishing solutions with a variety of boundary conditions along the edges of the plate shape represented in Cartesian coordinates. Statically determined non-linear contact laws, for loading and unloading conditions, are coupled with the partial differential equations governing the dynamic response of the composite plate. These nonlinear governing equations of the contact-impact problem are decoupled according to the second order terms by the method of principal transformation and then solved numerically. The dynamic response of fully clamped cross-ply and a variation of simply support and fixed boundary conditions for angle ply and cross ply composite laminates was investigated. Rectangular, circular and elliptical plates where analysed using the developed numetrical scheme. A total analysis method is used in order to check the efficiency of the results. Starting from the homogeneous formulation of the structure the eigenfrequencies and eigenshapes were calculated for all material and boundary conditions cases. Using the “transformation of the principal” method all cases were subjected in force and pressure pulses in order to check the efficiency of the dynamic solver developed. The efficiency of the present method was investigated comparing the respective results with well-known benchmark problems, experimental results and FEM analysis using MSC-NASTRAN and LSDYNA 3D codes. Following this exhaustive analysis, the coupled contact-impact problem was formulated and numerous contact schemes were applied to the problem. For the impact problem the type of the response is intimately related to the contact stiffness, impact energy, ply thickness and orientation, geometrical configuration and boundary conditions, as well as the orthotropic material properties and the longitudinal, transverse and interlaminar strength both in compression and in tension. Three types of dynamic responses were found to exist, local, global and transitional. All types are categorised and analysed. Using Buckingham Theorem the numerous parameters of the dynamic contact-impact problem were grouped to non-dimensional groups in order to provide a flexible analysis scheme for theoretical and experimental evaluation of the low velocity impact problem. A general normalised graph for characterising the impact response of structures was constructed. The numerous advantages for analysing impact phenomena using the normalised curve developed were presented. Moreover its use for predicting composite material damage area on low velocity impact conditions was presented.

The use of energy methods and variational principles is widespread in many fields of engineering of which structural mechanics and curve and surface design are two prominent examples. In principle many different types of function can be used as possible trial solutions to a given variational problem but where piecewise polynomial behaviour and user controlled cross segment continuity is either required or desirable, B-splines serve as a natural choice. Although there are many examples of the use of B-splines in such situations there is no common thread running through existing formulations that generalises from the one dimensional case through to two and three dimensions. We develop a unified approach to the representation of the minimisation equations for B-spline based functionals in tensor product form and apply these results to solving specific problems in geometric smoothing and finite element analysis using the Rayleigh-Ritz method. We focus on the development of algorithms for the exact computation of the minimisation matrices generated by finding stationary values of functionals involving integrals of squares and products of derivatives, and then use these to seek new variational based solutions to problems in the above fields. By using tensor notation we are able to generalise the methods and the algorithms from curves through to surfaces and volumes. The algorithms developed can be applied to other fields where a variational form of the problem exists and where such tensor product B-spline functions can be specified as potential solutions.

Please read the abstract in the section 00front of this document / Dissertation (MEng)--University of Pretoria, 2010. / Mechanical and Aeronautical Engineering / unrestricted

In this study, Ritz method has been employed to analyze the following problems: free vibrations of plates with curvilinear stiffeners, the lowest 100 frequencies of thick isotropic plates, free vibrations of thick quadrilateral laminates and free vibrations and static deformations of rectangular laminates, and sandwich structures. Admissible functions in the Ritz method are chosen as a product of the classical Jacobi orthogonal polynomials and weight functions that exactly satisfy the prescribed essential boundary conditions while maintaining orthogonality of the admissible functions. For free vibrations of plates with curvilinear stiffeners, made possible by additive manufacturing, both plate and stiffeners are modeled using a first-order shear deformation theory. For the thick isotropic plates and laminates, a third-order shear and normal deformation theory is used. The accuracy and computational efficiency of formulations are shown through a range of numerical examples involving different boundary conditions and plate thicknesses. The above formulations assume the whole plate as an equivalent single layer. When the material properties of individual layers are close to each other or thickness of the plate is small compared to other dimensions, the equivalent single layer plate (ESL) theories provide accurate solutions for vibrations and static deformations of multilayered structures. If, however, sufficiently large differences in material properties of individual layers such as those in sandwich structure that consists of stiff outer face sheets (e.g., carbon fiber-reinforced epoxy composite) and soft core (e.g., foam) exist, multilayered structures may exhibit complex kinematic behaviors. Hence, in such case, C_z⁰ conditions, namely, piecewise continuity of displacements and the interlaminar continuity of transverse stresses must be taken into account. Here, Ritz formulations are extended for ESL and layerwise (LW) Nth-order shear and normal deformation theories to model sandwich structures with various face-to-core stiffness ratios. In the LW theory, the C⁰ continuity of displacements is satisfied. However, the continuity of transverse stresses is not satisfied in both ESL and LW theories leading to inaccurate transverse stresses. This shortcoming is remedied by using a one-step well-known stress recovery scheme (SRS). Furthermore, analytical solutions of three-dimensional linear elasticity theory for vibrations and static deformations of simply supported sandwich plates are developed and used to investigate the limitations and applicability of ESL and LW plate theories for various face-to-core stiffness ratios. In addition to natural frequency results obtained from ESL and LW theories, the solutions of the corresponding 3-dimensional linearly elastic problems obtained with the commercial finite element method (FEM) software, ABAQUS, are provided. It is found that LW and ESL (even though its higher-order) theories can produce accurate natural frequency results compared to FEM with a considerably lesser number of degrees of freedom. / Doctor of Philosophy / In everyday life, plate-like structures find applications such as boards displaying advertisements, signs on shops and panels on automobiles. These structures are typically nailed, welded, or glued to supports at one or more edges. When subjected to disturbances such as wind gusts, plate-like structures vibrate. The frequency (number of cycles per second) of a structure in the absence of an applied external load is called its natural frequency that depends upon plate's geometric dimensions, its material and how it is supported at the edges. If the frequency of an applied disturbance matches one of the natural frequencies of the plate, then it will vibrate violently. To avoid such situations in structural designs, it is important to know the natural frequencies of a plate under different support conditions. One would also expect the plate to be able to support the designed structural load without breaking; hence knowledge of plate's deformations and stresses developed in it is equally important. These require mathematical models that adequately characterize their static and dynamic behavior. Most mathematical models are based on plate theories. Although plates are three-dimensional (3D) objects, their thickness is small as compared to the in-plane dimensions. Thus, they are analyzed as 2D objects using assumptions on the displacement fields and using quantities averaged over the plate thickness. These provide many plate theories, each with its own computational efficiency and fidelity (the degree to which it reproduces behavior of the 3-D object). Hence, a plate theory can be developed to provide accurately a quantity of interest. Some issues are more challenging for low-fidelity plate theories than others. For example, the greater the plate thickness, the higher the fidelity of plate theories required for obtaining accurate natural frequencies and deformations. Another challenging issue arises when a sandwich structure consists of strong face-sheets (e.g., made of carbon fiber-reinforced epoxy composite) and a soft core (e.g., made of foam) embedded between them. Sandwich structures exhibit more complex behavior than monolithic plates. Thus, many widely used plate theories may not provide accurate results for them. Here, we have used different plate theories to solve problems including those for sandwich structures. The governing equations of the plate theories are solved numerically (i.e., they are approximately satisfied) using the Ritz method named after Walter Ritz and weighted Jacobi polynomials. It is shown that these provide accurate solutions and the corresponding numerical algorithms are computationally more economical than the commonly used finite element method. To evaluate the accuracy of a plate theory, we have analytically solved (i.e., the governing equations are satisfied at every point in the problem domain) equations of the 3D theory of linear elasticity. The results presented in this research should help structural designers.

This research discusses the mechanical characteristics of laminated piezoelectric actuators that are manufactured at an elevated temperature, to cure the adhesive bonding the layers together, or to cure the layers made of polymeric composite material, and then cooled to a service temperature. Mainly discussed are actuators that are composed of layers of passive materials and a layer of piezoelectric material. THUNDER (THin layer UNimorph ferroelectric DrivER and sensor) and LIPCA (LIghtweight Piezo-composite Curved Actuator) actuators, which consist of layers of metal, adhesive and piezoelectric material, and carbon-epoxy, glass-epoxy and piezoelectric material, respectively, are studied and investigated in detail to understand the thermal effects due to the elevated manufacturing temperature. Owing to the large out-of-plane deformations of the THUNDER actuators as a result of cooling to the service temperature, inclusion of geometric nonlinearities in the kinematic relations is taken into consideration for prediction of the thermally-induced deformations and residual stresses. The deformations and residual stresses are predicted by using a 23-term Rayleigh-Ritz approach and more rigorous, time-consuming, finite-element analyses performed with ABAQUS. The thermally-induced deformations of THUNDER actuators can result in multiple room-temperature manufactured shapes, whereas those of LIPCA actuators (LIPCA-C1 and LIPCA-C2) exhibit single room-temperature manufactured shape. Actuation responses of these actuators caused by a quasi-static electric field applied to the piezoelectric layer are also studied with the Rayleigh-Ritz approach. It is shown that geometrical nonlinearities play an important role in the actuation responses, and these nonlinearities can be controlled by the choice of actuator geometry and the materials in the passive layers. In addition, blocking forces representing load-carrying capability of THUNDER and LIPCA actuators are determined. Support conditions and again geometrical nonlinearities are vital factor in load-resisting performances. Amongst the actuators considered, the actuated deflection and blocking forces are compared. Finally, based on the outcome of this study, new criteria for designing a new type of laminated piezoelectric actuators with improvement of performance characteristics are proposed. / Ph. D.

In the past few years a new class of layered piezoceramic and piezoceramic-composite actuators, known as RAINBOW and GRAPHBOW, respectively, that are capable of achieving 100 times greater out-of-plane displacements than previously available has been developed. Prior to the development of RAINBOW and GRAPHBOW, large stacks of piezoelectric actuators, requiring complicated electronic drive circuits, were necessary to achieve the displacement now possible through the use of a single RAINBOW actuator. The major issues with both RAINBOW and GRAPHBOW are the prediction of their room-temperature shapes after processing, and their deformation response under application of electric field. In this research, a methodology for predicting the manufactured shapes of rectangular and disk-style RAINBOW and GRAPHBOW is developed. All of the predictive analyses developed are based on finding approximate displacement responses that minimize the total potential energy of the devices through the use of variational methods and the Rayleigh-Ritz technique. These analyses are based on classical layered plate theory and assumed the various layers exhibited linear elastic, temperature-independent behavior. Geometric nonlinearities are important and are included in the strain-displacement relations. Stability of the predicted shapes is determined by examining the second variation of the total potential energy. These models are easily modified to account for the deformations induced by actuation of the piezoceramic. The results indicate that for a given set of material properties, rectangular RAINBOW can have critical values of sidelength-to-thickness ratio (Lx/H or Ly/H) below which RAINBOW exhibits unique, or single-valued, spherical or domed shapes when cooled from the processing temperature to room temperature. For values of sidelength-to-thickness ratio greater than the critical value, RAINBOW exhibits multiple room-temperature shapes. Two of the shapes are stable and are, in general, near-cylindrical. The third shape is spherical and is unstable. Similarly, disk-style RAINBOW can have critical values of radius-to-thickness ratios (R/H) below which RAINBOW exhibits axisymmetric room-temperature shapes. For values of R/H greater than the critical value, disk-style RAINBOW exhibits two stable near-cylindrical shapes and one unstable axisymmetric shape. Moreover, it is found that for the set of material properties used in this study, the optimal reduced layer thickness would be at 55%, since the maximum change in curvature is achieved under the application of an electric field, while the relationship between the change in curvatures and the electric field is kept very close to being linear. In general, good agreement is found for comparisons between the predicted and manufactured shapes of RAINBOW. A multi-step thermoelastic analysis is developed to model the addition of the fiber-reinforced composite layer to RAINBOW to make GRAPHBOW. Results obtained for rectangular RAINBOW indicate that if the bifurcation temperature in the temperature-curvature relation is lower than the composite cure temperature, then a unique stable GRAPHBOW shape can be obtained. If the RAINBOW bifurcation temperature is above the composite cure temperature, multiple room-temperature GRAPHBOW shapes are obtained and saddle-node bifurcations can be encountered during the cooling to room temperature of [0°/RAINBOW], [RAINBOW/0o], and [0o2/RAINBOW]. Rectangular [RAINBOW/0o/90o] seems to be less likely to encounter saddle-node bifurcations. Furthermore, the unstable spherical RAINBOW configuration is converted to a stable near-cylindrical configuration. For the case considered of disk-style GRAPHBOW, three stable room-temperature shapes are obtained and the unstable axisymmetric RAINBOW configuration is also converted to a stable near-cylindrical configuration. For both rectangular and disk-style GRAPHBOW, the relationship between the major curvature and the electric field is shown to be very close to being linear. This characteristic would aid any dynamic analysis of RAINBOW or GRAPHBOW. / Ph. D.

A pressurized arch-shell structural component made of flexible material is considered. The component is inflated with high internal pressure. The behavior of similar types of structures, such as a pair of leaning pressurized arches and pressurized arch-supported membrane shelters, has been investigated in the past. More recently, several types of pressurized structures have been incorporated as part of the framework for a variety of structural systems. Particularly, the U.S. Army has been investigating the use of large lightweight and transportable pressurized arch-shell structures to be used as maintenance shelters for vehicles, helicopters, and airplanes. The formulated equations using thin shell theory are applied to a pressurized arch-shell component. A numerical investigation based on the Rayleigh-Ritz method is utilized to determine the behavior of arch-shells under various types of loading. The types of loading include a uniformly distributed vertical load representing snow, a wind load, and a horizontal side load distributed along the arc length. Deflections, stress resultants, and moments at various locations are computed for two types of shapes: circular and non-circular arch-shells. / Master of Science

L’emploi de matériaux composites à matrice organique (CMO) pour la réalisation de structures aéronautiques « tièdes », peut exposer ces matériaux à l’action d’environnements agressifs, qui peuvent entrainer des phénomènes de vieillissement et de dégradation sévères associés à la diffusion d’espèces au sein du réseau macromoléculaire des matrices polymères. La complexité de la microstructure des CMO utilisés pour ces applications peut complexifier la compréhension de phénomènes de dégradation. Le vieillissement humide des CMO préoccupe particulièrement les industriels du secteur aéronautique ; la diffusion de l’eau dans la matrice polymère du composite peut entrainer des phénomènes de gonflement hygroscopique, des modifications des propriétés mécaniques. Des méthodes expérimentales existent pour la caractérisation de ces phénomènes et pour l’identification des paramètres associés : ces méthodes préconisent souvent l’emploi d’éprouvettes saturées en humidité,ce qui nécessite de longs temps de conditionnement et un nombre élevé d’échantillons. Cette thèse a pour but d’établir des protocoles d’identification rapide des propriétés diffuso-mécaniques de matériaux polymères et CMO pour applications aéronautiques. La démarche mise en place dans cette thèse s’organise autour de quatre chapitres. Le premier chapitre présente une étude bibliographique sur les outils de modélisation des couplages diffuso-mécaniques et sur les méthodes de caractérisation des propriétés diffuso-mécaniques. L’étude bibliographique permet de préciser le cadre de travail, qui prévoit l’emploi d’un modèle diffuso-mécanique faiblement couplé où la diffusion d’eau suit la loi de Fick et le comportement mécanique est hygroélastique linéaire, dépendant de la concentration en eau. Le second chapitre présente la mise en place et le développement d’une méthode d’identification rapide des propriétés de diffusion anisotrope, adaptée à des CMO à architecture complexe.La méthode s’appuie sur des mesures de prise de masse d’échantillons de CMO, l’anisotropie de diffusion de ces matériaux est obtenue par rotation des axes principaux d’orthotropie. La méthode proposée représente une extension de la « méthode de la pente » introduite par Shen et Springer pour l’identification des propriétés de diffusion de matériaux orthotropes, et basée sur l’exploitation des courbes gravimétriques aux temps courts. A travers cette méthode, les coefficients principaux et les axes principaux d’orthotropie peuvent être identifiés. Une discussion sur les conditions d’équivalence de la diffusion 3D à la diffusion 1D en fonction de l’épaisseur de l’échantillon est également présentée à la fin de ce chapitre. Le troisième chapitre explore à travers une étude numérique la possibilité d’identifier rapidement les propriétés mécaniques affectées par la concentration en eau de matériaux polymères à travers des essais mécaniques sur plaques minces avec des gradients de concentration en eau. Des essais de traction et de flexion sont considérés. Pour des matériaux isotropes, dans un cadre hygroélastique, il est montré que cette méthode permet l’identification du module d’Young et du coefficient de Poisson dépendant de la concentration en eau avec un gain remarquable de temps d’essai par rapport à des essais sur des échantillons saturés en humidité. Enfin, le dernier chapitre propose à travers une étude numérique une méthode d’identification rapide des propriétés diffuso-mécaniques de matériaux isotropes basée sur l’emploi de plaques sollicitées par un champ asymétrique de concentration en eau. L’identification est ainsi effectuée à partir du suivi des déflexions engendrées par les champs de concentration. Le coefficient de dilatation hygroscopique et le module d’Young dépendant de la concentration en eau peuvent être identifiés durant l’essai, à l’état transitoire du conditionnement, avec un gain remarquable du temps d’essai par rapport à des essais sur échantillons saturés en humidité. / The use of organic matrix composite materials (OMC) for the realization of "warm" aeronautical structures, may expose these materials to aggressive environments: wet or gaseous environments,high temperatures, which may promote severe aging and degradation phenomena related to species diffusion within the macromolecular network of the polymer matrices. The complexity of the OMC microstructure used for these applications can complicate the understanding of degradation phenomena : for example, species diffusion can be isotropic, orthotropic or anisotropic, depending on the texture of the fibrous reinforcement. Humid aging of OMC is of particular concern for the aeronautical industry ; the diffusion of water in the polymer matrix of the composite may promote hygroscopic swelling, changes in mechanical properties (stiffness, strength). Experimental methods exist for the characterization of these phenomena and for the identification of the associated parameters : these methods often recommend the use of moisture saturated specimens, which require long conditioning times, sometimes often a relevant number of samples (1 sample for each saturated state), high costs. The aim of this work is to establish protocols for fast identification of the diffusomechanical properties of polymers and polymer based OMC materials for aeronautical applications. The approach implemented in this thesis is organized in four chapters. The first chapter presents a bibliographic study on coupled diffuso-mechanics modeling tools and on methods of characterization/identification of diffuso-mechanical properties, more particularly for OMC for aeronautical applications. The bibliographic study allows specifying the framework of the present research, which foresees the employment of a weakly coupled diffuso-mechanical model, where water diffusion follows the Fick’s law and the mechanical behavior is linear hygroelastic, depending on water content. The second chap-ter presents the setting up and the development of a method for fast identification of anisotropic diffusion properties, suitable for OMC with complex architecture, such as for instance, 2D or 3Dwoven OMC. The method relies on mass-gain measures of OMC samples, the diffusion anisotropy ofthese materials is obtained by rotating the axes of orthotropy. The proposed method represents an extension of the "slope method" introduced by Shen and Springer in the 1970s for the identification of the diffusion properties of orthotropic materials (such as laminated composites), and is based onthe exploitation of gravimetric curves at short times. Through this method, the principal coefficients and the principal axes of orthotropy can be identified. A discussion about the transition from 3Dto 1D diffusion as a function of the sample geometry is also presented at the end of this chapter. The third chapter explores through a numerical study the possibility of identifying in a fast way the mechanical properties affected by moisture of polymeric materials by the use of mechanical tests on thin plates with water concentration gradients. Traction and bending tests are taken into account.For isotropic materials, in a hygroelastic setting, it is showed that this method allows identifyng the water concentration dependent Young’s modulus and the Poisson’s ratio with a remarkable time gain compared to tests on moisture saturated samples. Finally, the last chapter proposes through a numerical study a method for fast identification of the diffuso-mechanical properties of isotropic materials based on the use of plates loaded by an asymmetric water concentration field. The identification is thus carried out from the monitoring deflections generated by the concentration fields. The moisturedependent hygroscopic expansion coefficient and Young’s modulus can be identified during the test,by exploiting the transient state of conditioning, with a remarkable time gain compared with moisture saturated samples.

The aerospace industry is increasingly coming to the conclusion that physics-based high-fidelity models need to be used as early as possible in the design of its products. At the preliminary design stage of wing structures, though highly desirable for its high accuracy, a detailed finite element analysis(FEA) is often not feasible due to the prohibitive preparation time for the FE model data and high computation cost caused by large degrees of freedom. In view of this situation, often equivalent beam models are used for the purpose of obtaining global solutions. However, for wings with low aspect ratio, the use of equivalent beam models is questionable, and using an equivalent plate model would be more promising. An efficient method, Equivalent Plate Analysis or simply EPA, using an equivalent plate model, is developed in the present work for studying the static and free-vibration problems of built-up wing structures composed of skins, spars, and ribs. The model includes the transverse shear effects by treating the built-up wing as a plate following the Reissner-Mindlin theory (FSDT). The Ritz method is used with the Legendre polynomials being employed as the trial functions. Formulations are such that there is no limitation on the wing thickness distribution. This method is evaluated, by comparing the results with those obtained using MSC/NASTRAN, for a set of examples including both static and dynamic problems. The Equivalent Plate Analysis (EPA) as explained above is also used as a basis for generating other efficient methods for the early design stage of wing structures, such that they can be incorporated with optimization tools into the process of searching for an optimal design. In the search for an optimal design, it is essential to assess the structural responses quickly at any design space point. For such purpose, the FEA or even the above EPA, which establishes the stiffness and mass matrices by integrating contributions spar by spar, rib by rib, are not efficient enough. One approach is to use the Artificial Neural Network (ANN), or simply called Neural Network (NN) as a means of simulating the structural responses of wings. Upon an investigation of applications of NN in structural engineering, methods of using NN for the present purpose are explored in two directions, i.e. the direct application and the indirect application. The direct method uses FEA or EPA generated results directly as the output. In the indirect method, the wing inner-structure is combined with the skins to form an "equivalent" material. The constitutive matrix, which relates the stress vector to the strain vector, and the density of the equivalent material are obtained by enforcing mass and stiffness matrix equities with regard to the EPA in a least-square sense. Neural networks for these material properties are trained in terms of the design variables of the wing structure. It is shown that this EPA with indirect application of Neural Networks, or simply called an Equivalent Skin Analysis (ESA) of the wing structure, is more efficient than the EPA and still fairly good results can be obtained. Another approach is to use the sensitivity techniques. Sensitivity techniques are frequently used in structural design practices for searching the optimal solutions near a baseline design. In the present work, the modal response of general trapezoidal wing structures is approximated using shape sensitivities up to the second order, and the use of second order sensitivities proved to be yielding much better results than the case where only first order sensitivities are used. Also different approaches of computing the derivatives are investigated. In a design space with a lot of design points, when sensitivities at each design point are obtained, it is shown that the global variation in the design space can be readily given based on these sensitivities. / Ph. D.