Dissertations / Theses on the topic 'Variational asymptotic method'

Modern technologies require the materials with combinations of properties that can not be met by conventional single phase materials. This requirement leads to the development of composite materials or other materials with engineered microstructures, such as polymer composites and nanotube. Though the well-established finite element analysis (FEA) has the ability to analyze a small portion of such material, for the whole structure, the total degrees of freedom of a finite element model can easily exceed the bearable time in analysis or the capability of the best mainstream computers. To reduce the total degrees of freedom and save the computational efforts, an efficient way is to use a simpler and coarser mesh at the structure level with the micro level complexities captured by a homogenization method. Throughout the dissertation, the homogenization is carried on by variational asymptotic method which has been developed recently as the Variational Asymptotic Method for Unit Cell Homogenization (VAMUCH). This methodology is also expandable to the structure analysis as long as a representative structural element (RSE) can be obtained from structure. In the present research, the following problems are handled: (1) Maximizing the flexibility of choosing a RSE; (2) Bounding the effective properties of a random RSE; (3) Obtaining the equivalent plate stiffnesses for a corrugated plate from a RSE; (4) Extending the shell element of relative degree of freedom to analyze thin-walled RSE. These problems covered some important topics in homogenization theory. Firstly, the rules need to be followed when choosing a unit cell from a structure that can be homogenized. Secondly, for a randomly packed structure, the efficient way to predict effective material properties is to predict their bounds. Then, the composite material homogenization and the structural homogenization can be unied from a mathematical point of view, thus the repeating structure can be always simplified by the homogenization method. Lastly, the efficiency of analyzing thin-walled structures has been enhanced by the new type of shell element. In this research, the first two topics have been solved numerically through the finite element method under the framework of VAMUCH. The third one has been solved both analytically and numerically, and in the last, a new type of element has been implemented in VAMUCH to adapt the characteristics of a thin-walled problem. Numerous examples have demonstrated VAMUCH application and accuracy as a general-purpose analysis tool.

This dissertation developed a method that can accurately and efficiently capture the response of a structure by rigorous combination of a reduced-dimensional beam finite element model with a model based on full two-dimensional (2D) or three-dimensional (3D) finite elements. As a proof of concept, a joint 2D-beam approach is studied for planar-inplane deformation of strip-beams. This approach is developed for obtaining understanding needed to do the joint 3D-beam model. A Matlab code is developed to solve achieve this 2D-beam approach. For joint 2D-beam approach, the static response of a basic 2D-beam model is studied. The whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 2D model. The free end part is constructed as a beam model. To assemble the two different dimensional model, a transformation matrix is used to achieve deflection continuity or load continuity at the interface. After the transformation matrix from deflection continuity or from load continuity is obtained, the 2D part and the beam part can be assembled together and solved as one linear system. For a joint 3D-beam approach, the static and dynamic response of a basic 3D-beam model is studied. A Fortran program is developed to achieve this 3D-beam approach. For the uniform beam constrained at the root end, similar to the joint 2D-beam analysis, the whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 3D model. The free end part is constructed as a beam model. To assemble the two different dimensional models, the approach of load continuity at the interface is used to combine the 3D model with beam model. The load continuity at the interface is achieved by stress recovery using the variational-asymptotic method. The beam properties and warping functions required for stress recovery are obtained from VABS constitutive analysis. After the transformation matrix from load continuity is obtained, the 3D part and the beam part can be assembled together and solved as one linear system. For a non-uniform beam example, the whole structure is divided into several parts, where the root end and the non-uniform parts are constructed as 3D models and the uniform parts are constructed as beams. At all the interfaces, the load continuity is used to connect 3D model with beam model. Stress recovery using the variational-asymptotic method is used to achieve the load continuity at all interfaces. For each interface, there is a transformation matrix from load continuity. After we have all the transformation matrices, the 3D parts and the beam parts are assembled together and solved as one linear system.

Thesis (Ph. D.)--Aerospace Engineering, Georgia Institute of Technology, 2008.<br>Committee Chair: Bauchau, Olivier; Committee Member: Craig, James; Committee Member: Hodges, Dewey; Committee Member: Mahfuz, Hassan; Committee Member: Volovoi, Vitali.

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted, moderately-thick beam having rectangular cross sections and made of transversely isotropic material. The beam is modelled with-out assumptions from 3-D elasticity. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method (VAM) is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimisation of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Composites are materials which cater to the present and future needs of many demanding industries, such as aerospace, as they are weight-sensitive for a given requirement of strength and stiff ness, corrosion resistant, potentially multi-functional and can be tailored according to the application. However, they are in particular difficult to join as they cannot be easily machined, without introducing damages which can eventually grow. Any structure is as strong as its weakest joint. Most of the joints belong to the category of mechanically-fastened joints and they pose enormous challenges in modeling due to contact phenomena, nonlinearity and stress concentration factors. It is therefore a necessity to construct an efficient model that would include all the relevant contact phenomena in the joints, as it has been pointed out in literature that damage typically initiates near the joint holes. The focus of this work is to describe the construction of an asymptotically-correct model using the Variational Asymptotic Method (VAM). Amongst its many potential applications, VAM is a well-established analytical tool for obtaining the stress and strain fields for beams and shells. The methodology takes advantage of the small parameter that is inherent in the problem, such as the ratio of certain characteristic dimensions of the structure. In shells and beams, VAM takes advantage of the dimension-based small parameter(s), thereby splitting the problem into 2-D + 1-D (for beams) and 1-D + 2-D (for shells), in turn offering very high computational efficiency with very little loss of accuracy compared to dimensionally unreduced 3-D models. In this work, the applicability of VAM is extended to two-dimensional (2-D) and three-dimensional (3-D) frictionless contact problems. Since a generalised VAM model for contact has not been pursued before, the `phantom0 step is adopted for both 2-D and 3-D models. The development of the present work starts with the construction of a 2-D model involving a large rectangular plate being pressed against a rigid frictionless pin. The differential equations governing the problem and the associated boundary conditions are obtained by minimizing the reduced strain energy, augmented with the appropriate gap function, by using a penalty method. The model is developed for both isotropic and orthotropic cases. The boundary value problem is solved numerically and the displacement field obtained is compared with the one obtained using commercial software (ABAQUSr) for validation at critical regions such as the contact surfaces. Banking on the validation of the 2-D model, a 3-D model with a pin and a finite annular cylinder was constructed. The strain energy for the finite cylinder was derived using geometrically exact 3-D kinematics and VAM was applied leading to the reduction in the strain energy for isotropic and orthotropic materials in rectangular and cylindrical co-ordinates. As in the 2-D case, the reduced strain energy, subject to the inequality constraint of the gap function, is minimized with respect to the displacement field and the corresponding boundary value problem is solved numerically. The displacements of the contact surface and the top surface of the annular cylinder are compared with those from ABAQUS and thus validated. The displacement fields obtained using the current 2-D and 3-D models show very good agreement with those from commercial finite element software packages. The model could be re ned further by using the gap function derived in this work and applying it to a plate model based on VAM, which could be explored in the future.

In aerospace applications, use of laminates made of composite materials as face sheets in sandwich panels are on the rise. These composite laminates have low transverse shear and transverse normal moduli compared to the in-plane moduli. It is also seen that the corresponding transverse strength values are very low compared to the in-plane strength leading to delaminations. Further, in sandwich structures, the core is subjected to significant transverse shear stresses. Therefore the interlaminar stresses (i.e., transverse shear and normal) can govern the design of sandwich structures. As a consequence, the first step in achieving efficient designs is to develop the ability to reliably estimate interlaminar stresses. Stress analysis of the composite sandwich structures can be carried out using 3-D finite elements for each layer. Owing to the enormous computational time and resource requirements for such a model, this process of analysis is rendered inefficient. On the other hand, existing plate/shell finite elements, when appropriately chosen, can help quickly predict the 2-D displacements with reasonable accuracy. However, their ability to calculate the thickness-wise distributions of interlaminar shear and normal stresses and 3-D displacements remains as a research goal. Frequently, incremental refinements are offered over existing solutions. In this scenario, an asymptotically correct dimensional reduction from 3-D to 2-D, if possible, would serve to benchmark any ongoing research. The employment of a mathematical technique called the Variational Asymptotic Method (VAM) ensures the asymptotical correctness for this purpose. In plates and sandwich structures, it is typically possible to identify (purely from the defined material distributions and geometry) certain parameters as small compared to others. These characteristics are invoked by VAM to derive an asymptotically correct theory. Hence, the 3-D problem of plates is automatically decomposed into two separate problems (namely 1-D+2-D), which then exchange relevant information between each other in both ways. The through-the-thickness analysis of the plate, which is a 1-D analysis, provides asymptotic closed form solutions for the 2-D stiffness as well as the recovery relations (3-D warping field and displacements in terms of standard plate variables). This is followed by a 2-D plate analysis using the results of the 1-D analysis. Finally, the recovery relations regenerate all the required 3-D results. Thus, this method of developing reduced models involves neither ad hoc kinematic assumptions nor any need for shear correction factors as post-processing or curve-fitting measures. The results are most general and can be made as accurate as desired, while the procedure is computationally efficient. In the present work, an asymptotically correct plate theory is formulated for composite sandwich structures. In developing this theory, in addition to the small parameters (such as small strains, small thickness-to-wavelength ratios etc.,) pertaining to the general plate theory, additional small parameters characterizing (and specific to) sandwich structures (viz., smallness of the thickness of facial layers com-pared to that of the core and smallness of elastic material stiffness of the core in relation to that of the facesheets) are used in the formulation. The present approach also satisfies the interlaminar displacement continuity and transverse equilibrium requirements as demanded by the exact 3-D formulation. Based on the derived theory, numerical codes are developed in-house. The results are obtained for a typical sandwich panel subjected to mechanical loading. The 3-D displacements, inter-laminar normal and shear stress distributions are obtained. The results are compared with 3-D elasticity solutions as well as with the results obtained using 3-D finite elements in MSC NASTRAN®. The results show good agreement in spite of the major reduction in computational effort. The formulation is then extended for thermo-elastic deformations of a sandwich panel. This thesis is organized chronologically in terms of the objectives accomplished during the current research. The thesis is organized into six chapters. A brief organization of the thesis is presented below. Chapter-1 briefly reviews the motivation for the stress analysis of sandwich structures with composite facesheets. It provides a literature survey on the stress analysis of composite laminates and sandwich plate structures. The drawbacks of the existing anlaytical approaches as opposed to that of the VAM are brought out. Finally, it concludes by listing the main contributions of this research. Chapter-2 is dedicated to an overview of the 3-D elasticity formulation of composite sandwich structures. It starts with the 3-D description of a material point on a structural plate in the undeformed and deformed configurations. Further, the development of the associated 3-D strain field is also described. It ends with the formulation of the potential energy of the sandwich plate structure. Chapter-3 develops the asymptotically correct theory for composite sandwich plate structure. The mathematical description of VAM and the procedure involved in developing the dimensionally reduciable structural models from 3-D elasticity functional is first described. The 1-D through-the-thickness analysis procedure followed in developing the 2-D plate model of the composite sandwich structure is then presented. Finally, the recovery relations (which are one of the important results from 1-D through-the-thickness analysis) to extract 3-D responses of the structure are obtained. The developed formulation is applied to various problems listed in chapter 4. The first section of this chapter presents the validation study of the present formulation with available 3-D elasticity solutions. Here, composite sandwich plates for various length to depth ratios are correlated with available 3-D elasticity solutions as given in [23]. Lastly, the distributions of 3-D strains, stresses and displacements along the thickness for various loadings of a typical sandwich plate structure are correlated with corresponding solutions using well established 3-D finite elements of MSC NASTRAN® commerical FE software. The developed and validated formulation of composite sandwich structure for mechanical loading is extended for thermo-elastic deformations. The first sections of this chapter describes the seamless inclusion of thermo-elastic strains into the 3-D elasticity formulation. This is followed by the 1-D through-the-thickness analysis in developing the 2-D plate model. Finally, it concludes with the validation of the present formulation for a very general thermal loading (having variation in all the three co-ordinate axes) by correlating the results from the present theory with that of the corresponding solutions of 3-D finite elements of MSC NASTRAN® FE commercial software. Chapter-6 summarises the conclusions of this thesis and recommendations for future work.

Composite materials show heterogeneity at different length scales. hence concurrent multiscale analysis is the only reliable method to analyze them. But unfortunately there is no concurrent multi-scale strategy that is efficient, and accurate while addressing all kinds of problems. This lack of reliability is partly because there is no micro-mechanical model which inherently keeps all relevent global information with it. This thesis tries to fill this gap. The presented micro-mechanical model not only homogenizes the micro-structure but also keeps the global information with it. Most of the micro-mechanical models in the literature extract the Representative Volume Element (RVE) from the continuum for analysis which results in loss of information and accuracy. In the present approach also, the RVE has been extracted from the continuum but with the major difference that all the macro/meso-scopic parameters are accounted for. Five macro/meso-scopic one dimensional parameters have been defined which completely define the effect of continuum. 11 for one dimensional stretch, _1 for torsion, __ (_ = 2, 3) for bending and _33 for uniform pressurization due to the presence of the continuum. Further, the above macro/meso-scopic parameters are proven, by the asymptotic, theory to be constant at a cross section but vary, in general, over the length of the fiber. Hence, the analysis is valid for any location and is not restricted to any local domain. Three major problems have been addressed: • Homogenization and analysis of RVE without any defects • Homogenization and analysis of RVE with fiber-matrix de-bonding • Homogenization and analysis of RVE with radial matrix cracking. Variational Asymptotic Method (VAM) has been used to solve the above mentioned problems analytically. The results have been compared against standard results in the literature and against 3D FEA. At the end, results for “Radial deformation due to torsion” problem will be presented which was solved “accidentally.”

This work aims at dimensional reduction of nonlinear material models in an asymptotically accurate manner. The three-dimensional(3-D) nonlinear material models considered include isotropic, orthotropic and dielectric compressible hyperelastic material models. Hyperelastic materials have potential applications in space-based inflatable structures, pneumatic membranes, replacements for soft biological tissues, prosthetic devices, compliant robots, high-altitude airships and artificial blood pumps, to name a few. Such structures have special engineering properties like high strength-to-mass ratio, low deflated volume and low inflated density. The majority of these applications imply a thin shell form-factor, rendering the problem geometrically nonlinear as well. Despite their superior engineering properties and potential uses, there are no proper analysis tools available to analyze these structures accurately yet efficiently. The development of a unified analytical model for both material and geometric nonlinearities encounters mathematical difficulties in the theory but its results have considerable scope. Therefore, a novel tool is needed to dimensionally reduce these nonlinear material models. In this thesis, Prof. Berdichevsky’s Variational Asymptotic Method(VAM) has been applied rigorously to alleviate the difficulties faced in modeling thin shell structures(made of such nonlinear materials for the first time in the history of VAM) which inherently exhibit geometric small parameters(such as the ratio of thickness to shortest wavelength of the deformation along the shell reference surface) and physical small parameters(such as moderate strains in certain applications). Saint Venant-Kirchhoff and neo-Hookean 3-D strain energy functions are considered for isotropic hyperelastic material modeling. Further, these two material models are augmented with electromechanical coupling term through Maxwell stress tensor for dielectric hyperelastic material modeling. A polyconvex 3-D strain energy function is used for the orthotropic hyperelastic model. Upon the application of VAM, in each of the above cases, the original 3-D nonlinear electroelastic problem splits into a nonlinear one-dimensional (1-D) through-the-thickness analysis and a nonlinear two-dimensional(2-D) shell analysis. This greatly reduces the computational cost compared to a full 3-D analysis. Through-the-thickness analysis provides a 2-D nonlinear constitutive law for the shell equations and a set of recovery relations that expresses the 3-D field variables (displacements, strains and stresses) through thethicknessintermsof2-D shell variables calculated in the shell analysis (2-D). Analytical expressions (asymptotically accurate) are derived for stiffness, strains, stresses and 3-D warping field for all three material types. Consistent with the three types of 2-D nonlinear constitutive laws,2-D shell theories and corresponding finite element programs have been developed. Validation of present theory is carried out with a few standard test cases for isotropic hyperelastic material model. For two additional test cases, 3-Dfinite element analysis results for isotropic hyperelastic material model are provided as further proofs of the simultaneous accuracy and computational efficiency of the current asymptotically-correct dimensionally-reduced approach. Application of the dimensionally-reduced dielectric hyperelastic material model is demonstrated through the actuation of a clamped membrane subjected to an electric field. Finally, the through-the-thickness and shell analysis procedures are outlined for the orthotropic nonlinear material model.

Research into advanced wind turbine design has shown that load alleviation strategies like bend-twist coupled blades and coned rotors could reduce costs. However these strategies are based on nonlinear aero-structural dynamics providing additional benefits to components beyond the blades. These innovations will require Multi-disciplinary Design Optimization (MDO) to realize the full benefits. This research expands the MDO capabilities of Horizontal Axis Wind Turbines. The early research explored the numerical stability properties of Blade Element Momentum (BEM) models. Then developed a provincial scale wind farm siting models to help engineers determine the optimal design parameters. The main focus of this research was to incorporate advanced analysis tools into an aero-elastic optimization framework. To adequately explore advanced designs with optimization, a new set of medium fidelity analysis tools is required. These tools need to resolve more of the physics than conventional tools like (BEM) models and linear beams, while being faster than high fidelity techniques like grid based computational fluid dynamics and shell and brick based finite element models. Nonlinear beam models based on Geometrically Exact Beam Theory (GEBT) and Variational Asymptotic Beam Section Analysis (VABS) can resolve the effects of flexible structures with anisotropic material properties. Lagrangian Vortex Dynamics (LVD) can resolve the aerodynamic effects of novel blade curvature. Initially this research focused on the structural optimization capabilities. First, it developed adjoint-based gradients for the coupled GEBT and VABS analysis. Second, it developed a composite lay-up parameterization scheme based on manufacturing processes. The most significant challenge was obtaining aero-elastic optimization solutions in the presence of erroneous gradients. The errors are due to poor convergence properties of conventional LVD. This thesis presents a new LVD formulation based on the Finite Element Method (FEM) that defines an objective convergence metric and analytic gradients. By adopting the same formulation used in structural models, this aerodynamic model can be solved simultaneously in aero-structural simulations. The FEM-based LVD model is affected by singularities, but there are strategies to overcome these problems. This research successfully demonstrates the FEM-based LVD model in aero-elastic design optimization.<br>Graduate<br>0548<br>pilot.mm@gmail.com