Dissertations / Theses on the topic '(finite deformation theory)'

Deformation theory studies the behavior of mathematical objects, such as representations or modules, under small perturbations. This theory is useful in both pure and applied mathematics and has been used in the proof of many long-standing problems. In particular, in number theory Wiles and Taylor used universal deformation rings of Galois representations in the proof of Fermat's Last Theorem. The main motivation for determining universal deformation rings of modules for finite dimensional algebras is that deep results from representation theory can be used to arrive at a better understanding of deformation rings. In this thesis, I study the universal deformation rings of certain modules for algebras of dihedral type of polynomial growth which have been completely classied by Erdmann and Skowronski using quivers and relations. More precisely, let κ be an algebraically closed field and let λ be a κ-algebra of dihedral type which is of polynomial growth. In this thesis, first classify all λ-modules whose stable endomorphism ring is isomorphic to κ and which are given combinatorially by strings, and then I determine the universal deformation ring of each of these modules.

The theory is based on an assumed displacement field, in which the surface displacements are expanded in powers of the thickness coordinate up to the third order. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory accounts for small strains but moderately large displacements (i.e., von Karman strains). Exact solutions for certain cross-ply shells and finite-element models of the theory are also developed. The finite-element model is based on independent approximations of the displacements and bending moments (i.e., mixed formulation), and therefore only C°-approximations are required. Further, the mixed variational formulations developed herein suggest that the bending moments can be interpolated using discontinuous approximations (across inter-element boundaries). The finite element is used to analyze cross-ply and angle-ply laminated shells for bending, vibration, and transient response. Numerical results are presented to show the effects of boundary conditions, lamination scheme (i.e., bending-stretching coupling and material anisotropy) shear deformation, and geometric nonlinearity on deflections and frequencies. Many of the numerical results presented here for laminated shells should serve as references for future investigations.
Ph. D.

The continuum theory of large deformation elasto-plasticity is summarized as far as it is necessary for the numerical treatment with the Finite-Element-Method. Using the calculus of modern differential geometry and functional analysis, the fundamental equations are derived and the proof of most of them is shortly outlined. It was not our aim to give a contribution to the development of the theory, rather to show the theoretical background and the assumptions to be made in state of the art elasto-plasticity.

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Mechanical Engineering, 2007.
This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.
Includes bibliographical references.
Amorphous thermoplastic polymers are important engineering materials; however, their nonlinear, strongly temperature- and rate-dependent elastic-visco-plastic behavior has, until now, not been very well understood. The behavior has previously been modeled with mixed success by existing constitutive theories. As a result, there is currently no generally agreed upon theory to model the large-deformation, thermo-mechanically coupled, elasto-visco-plastic response of amorphous polymeric materials spanning their glass transition temperatures. What is needed is a unified constitutive framework that is capable of capturing the transition from a visco-elastic-plastic solidlike response below the glass transition temperature, to a rubbery-viscoelastic response above the glass transition temperature, to a fluid-like response at yet higher temperatures. We have developed a continuum-mechanical constitutive theory aimed to fill this need. The theory has been specialized to represent the salient features of the mechanical response of poly(methyl methacrylate) in a temperature range spanning room temperature to 60C above the glass transition temperature #g 110C of the material, in a strain-rate range of 10-4/s to 10-1/s, and under compressive stress states in which this material does not exhibit crazing. We have implemented our theory in the finite element program ABAQUS/Explicit. The numerical simulation capability of the theory is demonstrated with simulations of the micron-scale hot-embossing process for manufacture of microfluidic devices.
by Nicoli Margaret Ames.
Ph.D.

Theoretical formulations, analytical solutions, and finite element solutions for laminated composite plate and shell structures with smart material laminae are presented in the study. A unified third-order shear deformation theory is formulated and used to study vibration/deflection suppression characteristics of plate and shell structures. The von K??rm??n type geometric nonlinearity is included in the formulation. Third-order shear deformation theory based on Donnell and Sanders nonlinear shell theories is chosen for the shell formulation. The smart material used in this study to achieve damping of transverse deflection is the Terfenol-D magnetostrictive material. A negative velocity feedback control is used to control the structural system with the constant control gain. The Navier solutions of laminated composite plates and shells of rectangular planeform are obtained for the simply supported boundary conditions using the linear theories. Displacement finite element models that account for the geometric nonlinearity and dynamic response are developed. The conforming element which has eight degrees of freedom per node is used to develop the finite element model. Newmark's time integration scheme is used to reduce the ordinary differential equations in time to algebraic equations. Newton-Raphson iteration scheme is used to solve the resulting nonlinear finite element equations. A number of parametric studies are carried out to understand the damping characteristics of laminated composites with embedded smart material layers.

Le nouvel élément de poutre est une évolution d'un élément de Timoshenko poutre avec un nœud supplémentaire situé à mi-longueur. Ce nœud supplémentaire permet l'introduction de trois composantes supplémentaires de contrainte afin que la loi constitutionnelle 3D complète puisse être utilisée directement. L'élément proposé a été introduit dans un code d'éléments finis dans Matlab et une série d'exemples de linéaires/petites contraintes ont été réalisées et les résultats sont systématiquement comparés avec les valeurs correspondantes des simulations ABAQUS/Standard 3D. Ensuite, la deuxième étape consiste à introduire le comportement orthotrope et à effectuer la validation de déplacements larges / petites contraintes basés sur la formulation Lagrangienne mise à jour. Une série d'analyses numériques est réalisée qui montre que l'élément 3D amélioré fournit une excellente performance numérique. En effet, l'objectif final est d'utiliser les nouveaux éléments de poutre 3D pour modéliser des fils dans une préforme composite textile. A cet effet, la troisième étape consiste à introduire un comportement de contact et à effectuer la validation pour un nouveau contact entre 3D poutres à section rectangulaire. La formulation de contact est dérivée sur la base de formulation de pénalité et de formulation Lagrangian mise à jour utilisant des fonctions de forme physique avec l'effet de cisaillement inclus. Un algorithme de recherche de contact efficace, qui est nécessaire pour déterminer un ensemble actif pour le traitement de contribution de contact, est élaboré. Et une linéarisation constante de la contribution de contact est dérivée et exprimée sous forme de matrice appropriée, qui est facile à utiliser dans l'approximation FEM. Enfin, on présente quelques exemples numériques qui ne sont que des analyses qualitatives du contact et de la vérification de l'exactitude et de l'efficacité de l'élément de 3D poutre proposé
The new beam element is an evolution of a two nodes Timoshenko beam element with an extra node located at mid-length. That extra node allows the introduction of three extra strain components so that full 3D stress/strain constitutive relations can be used directly. The second step is to introduce the orthotropic behavior and carry out validation for large displacements/small strains based on Updated Lagrangian Formulation. A series of numerical analyses are carried out which shows that the enhanced 3D element provides an excellent numerical performance. Indeed, the final goal is to use the new 3D beam elements to model yarns in a textile composite preform. For this purpose, the third step is introducing contact behavior and carrying out validation for new 3D beam to beam contact with rectangular cross section. The contact formulation is derived on the basis of Penalty Formulation and Updated Lagrangian formulation using physical shape functions with shear effect included. An effective contact search algorithm is elaborated. And a consistent linearization of contact contribution is derived and expressed in suitable matrix form, which is easy to use in FEM approximation. Finally, some numerical examples are presented which are only qualitative analysis of contact and checking the correctness and the effectiveness of the proposed 3D beam element

Thesis (Ph.D)--Aerospace Engineering, Georgia Institute of Technology, 2009.
Committee Chair: Hanagud,Sathya; Committee Member: Apetre, Nicoleta; Committee Member: Kardomateas, George; Committee Member: McDowell, David L.; Committee Member: Ruzzene, Massimo. Part of the SMARTech Electronic Thesis and Dissertation Collection.

Une etude numerique et theorique de la localisation de la deformation en bandes de cisaillement dans les geomateriaux est presentee. Cette etude comprend deux parties distinctes : ' la premiere partie de la detection du moment d'apparition de la localisation et plus generalement des phenomenes de bifurcation. L'etude est realisee dans le cadre des grandes transformations en utilisant le modele cloe. Dans ce but, un modele cloe von mises et un algorithme de detection de perte d'unicite globale, permettant la resolution d'un probleme aux limites formule en vitesses, ont ete developpes. L'influence de l'initialisation de l'algorithme est plus particulierement etudiee. L'etude numerique d'un essai biaxial, a permis de mettre en evidence plusieurs modes de bifurcation (modes de flambage, bandes de cisaillement). Pour les modes localises, une comparaison entre les resultats numeriques et les predictions theoriques fournies par le critere de localisation de cloe ont permis d'illustrer la fiabilite de ce dernier. ' la seconde partie concerne le suivi de la localisation, de son declenchement, jusqu'a la ruine complete de l'echantillon. Dans ce but, un modele d'interface a ete developpe. Il presente la particularite d'assurer une transition continue entre le regime de pre et de post localisation (concept de consistance) d'une part et de decrire l'evolution specifique de la densite dans la zone localisee d'autre part. Le concept d'indice des vides critique est introduit explicitement dans ce modele. Les resultats obtenus sur un essai biaxial refletent assez fidelement les resultats experimentaux. Ce modele a ensuite ete integre dans un code de calcul par elements finis. Un element d'interface, fonctionnant en grandes transformations a ete developpe et valide. Enfin, une pre-etude du probleme de propagation des bandes de cisaillement au sein d'une structure a ete effectuee.

This thesis involves the derivation of a constitutive model under large deformationtheory using Updated Lagrange method applied on an orthotropic material.Thefollowing aspects are included in this thesis work: introduction, theory, FEM implementation, derivation of constitutive model, calibration, result, discussion, conclusion and the future work. This thesis studies the deformation behavior of wood, which is widely used as aconstruction material, in an advanced and more detailed way by analyzing the mechanical properties of wood from both, the application in laboratory and theoreticalcalculation under large deformation theory. A non-linear elastic constitutive model is proposed, derived and calibrated using asimple inverse analysis procedure. The calibration process was performed to identify8 constitutive parameters A1 − A8 of the constitutive model by performing inverseanalysis against relevant experimental data acquired using the Aramis system. Theresults in the comparison were extracted from the specimen when it is both intangential orientation and radial orientation. The project work will be dedicated to the development of a Finite Element Method(FEM) code implemented in MATLAB scripts which was directly used to study themechanical properties of the orthotropic wood material when hyper-elastic behavioris assumed. The results will contain three parts: 1) study of the influence of pith location onthe load required to deform the specimen specimen, 2) reaction force comparisonof the model results against experimental results, and, 3) comparison of the GreenLagrangian strain pattern over the specimen between the experimental data and themodel’s results.

Une etude theorique et numerique de la localisation de la deformation en bandes de cisaillement dans les geomateriaux est presentee. L'etude est realisee en deux volets : ? le premier volet concerne le phenomene de la localisation au moment de son declenchement en utilisant le modele cloe. Un critere de controlabilite en deformation plane a ete etabli pour le modele cloe. Une etude de sensibilite parametrique du critere de bifurcation de cloe roche a ete menee sur la marne de beaucaire. Une etude numerique de modelisation en grandes deformations de l'excavation d'un puits fore dans une couche des marnes raides (marne a hydrobies) a ete menee avec la loi cloe avec l'identification des parametres de cloe pour cette marne sur des essais elementaires de laboratoire. Le critere de bifurcation a permis de montrer que la localisation se declenchait en paroi, symetriquement ou non suivant l'isotropie du champ lointain, et plus tot dans le cas anisotrope. ? le second volet concerne le suivi de la localisation. Une etude theorique des milieux continus de second gradient unidimensionnels est realisee. La resolution analytique des equations constitutives, associees a un modele unidimensionnel developpe, est presentee dans le cadre des petites deformations. Une analyse numerique des equations d'equilibre a permis la determination d'une matrice de rigidite consistante ainsi que l'expression des forces hors equilibre. Un code elements finis unidimensionnel avec des elements finis conformes a ete developpe et valide. L'introduction du second gradient parait non satisfaisante pour assurer l'unicite des solutions. Un resultat important est l'independance du resultat vis a vis du maillage. Une etude theorique et une analyse numerique semblable a celle du cas monodimensionnel est ensuite realisee pour les milieux continus de second gradient bidimensionnels dans le cadre des grandes deformations.

Au cœur de cette thèse est une théorie de continuum alternatif connue comme la théorie micropolaire, qui est développée pour décrire des phénomènes lesquels on ne peut pas décrire en utilisant la théorie classique. Dans cette théorie, en complément du champ de déplacement, il existe aussi un autre champ indépendant, celui de microrotation, et afin de pouvoir décrire complètement un tel matériau, six paramètres des matériaux sont nécessaires. Dans le cadre de la modélisation par éléments finis, nouveaux éléments fondés sur la théorie micropolaire dans les régimes linéaire et géométriquement non linéaire sont développés. Dans le cadre de l'analyse linéaire, les problèmes bi- et tri-dimensionnels sont analysés. En 2D, les nouvelles familles des éléments triangulaires et quadrilatères sont développés avec l'interpolation liée des champs cinématiques. Ensuite, la forme faible est étendue aux 3D, et un élément fini hexaédrique du premier ordre, avec le champ de déplacement enrichi avec des modes incompatibles est dérivé. Il est constaté que l'interpolation liée et les modes incompatibles améliore la précision par rapport à la précision des éléments finis micropolaires conventionnels. Dans le part non-linéaire, les éléments de premier et deuxième ordre avec l'interpolation conventionnelle sont développés. Pour tester la performance des éléments présentés, une solution analytique non-linéaire de la flexion pure est dérivée. Il est observé que les éléments convergent vers la solution dérivée. Les éléments sont testés sur les autres exemples où la dépendance du sentier et l'invariance de déformation sont détectés. Une procédure pour résoudre ces anomalies est présentée
At the core of this thesis is an alternative continuum theory called the micropolar (Cosserat) continuum theory, developed in order to describe the phenomena which the classical continuum theory is not able to describe. In this theory, in addition to the displacement field, there also exists an independent microrotation field and, in order to completely describe such a material, six material parameters are needed. In the framework of the finite-element method, new finite elements based on the micropolar continuum theory in both linear and geometrically non-linear analysis are developed using the displacement-based approach. In the linear analysis, both two- and three-dimensional set-ups are analysed. In 2D new families of triangular and quadrilateral finite elements with linked interpolation of the kinematic fields are derived. In order to assure convergence of the derived finite elements, they are modified using the Petrov-Galerkin approximation. Their performance is compared against existing conventional micropolar finite elements on a number of micropolar benchmark problems. It is observed that the linked interpolation shows enhanced accuracy in the bending test when compared against the conventional Lagrange micropolar finite element. Next, the weak formulation is extended to 3D and a first-order hexahedral finite element enhanced with the incompatible modes is derived. The element performance is assessed by comparing the numerical results against the available analytical solutions for various boundary value problems, which are shown to be significant for the experimental verification of the micropolar material parameters. It is concluded that the proposed element is highly suitable for the validation of the methodology to determine the micropolar material parameters. In the non-linear part, first- and second-order geometrically nonlinear hexahedral finite elements with Lagrange interpolation are derived. In order to test the performance of the presented finite elements, a pure-bending non-linear micropolar analytical solution is derived. It is observed that the elements converge to the derived solution. The elements are tested on three additional examples where the path-dependence and strain non-invariance phenomena are detected and assessed in the present context. A procedure to overcome the non-invariance anomaly is outlined

Der Vortrag beschreibt die zugrunde liegende Theorie und die Softwarefunktionalität des in PTC Creo Simulate 3.0 eingeführten Kontaktmodells mit endlicher Reibung und vergleicht dieses mit den bis Creo Simulate 2.0 exklusiv verwendeten Kontaktmodellen (ideal reibungsfrei und unendlich reibungsbehaftet). An zwei Modellbeispielen (ein von zwei Bremsbacken geklemmtes Bremsschwert und ein verschraubtes Schwungrad) wird versucht, die Funktionsweise des neuen Modells zu demonstrieren. Wegen aktueller Qualitätsprobleme der Software wird die Brauchbarkeit der Kontaktmodelle für den Anwender bewertet (Stand Creo 3.0 M080 / Creo 2.0 M200) und umfangreiches Feedback an den Softwarehersteller PTC gegeben
The presentation describes the underlying theory and software functionality of the finite friction contact model introduced with PTC Creo Simulate 3.0. It is being compared with the friction-free and infinite friction contact model used exclusively until Creo Simulate 2.0. It is being tried to demonstrate the mode of operation of the new model with help of two examples: A brake sword clamped by two brake pads and a bolted flywheel. Because of actual software quality problems, the usability of the contact model for the user is being rated (status Creo 3.0 M080 / Creo 2.0 M200). Furthermore, comprehensive feedback is given to the software developer PTC

Smart materials are very important because of their potential applications in the biomedical, petroleum and aerospace industries. They can be used to build systems and structures that self-monitor to function and adapt to new operating conditions. In this study, we are mainly interested in developing a computational framework for the analysis of plate structures comprised of composite or functionally graded materials (FGM) with embedded or surface mounted piezoelectric sensors/actuators. These systems are characterized by thermo-electro-mechanical coupling, and therefore their understanding through theoretical models, numerical simulations, and physical experiments is fundamental for the design of such systems. Thus, the objective of this study was to perform a numerical study of smart material plate structures using a refined plate theory that is both accurate and computationally economical. To achieve this objective, an improved version of the Reddy third-order shear deformation theory of plates was formulated and its finite element model was developed. The theory and finite element model was evaluated in the context of static and dynamic responses without and with actuators. In the static part, the performance of the developed finite element model is compared with that of the existing models in determining the displacement and stress fields for composite laminates and FGM plates under mechanical and/or thermal loads. In the dynamic case, coupled and uncoupled electro-thermo-mechanical analysis were performed to see the difference in the evolution of the mechanical, electrical and thermal fields with time. Finally, to test how well the developed theory and finite element model simulates the smart structural system, two different control strategies were employed: the negative velocity feedback control and the Least Quadratic Regulator (LQR) control. It is found that the refined plate theory provides results that are in good agreement with the those of the 3-D layerwise theory of Reddy. The present theory and finite element model enables one to obtain very accurate response of most composite and FGM plate structures with considerably less computational resources.

In the literature, there are various material models proposed so as to model the constitutive behavior of hyperelastic materials for example, St. Venant-Kirchho_ model, Mooney-Rivlin model etc. The stability of such material models under various states of deformation is of important concern, and generally stability analysis is conducted in homogeneous states of deformation. Within hyperelasticity, instabilities can be broadly classified as geometrical and material types. Geometrical instabilities such as buckling, symmetric bifurcation etc. are of physical origin, and lead to multiple solutions at critical stretch. Material instability is a aw in the material model and leads to unphysical solutions at the onset. It is required that the constitutive model should be materially stable i.e., should not give unphysical results, and be able to predict correctly the onset of geometrical instabilities. Certain constitutive restrictions proposed in the literature are inadequate to characterize such instabilities. In the work, we propose stability criteria which will characterize geometrical as well as material instabilities. A new elasticity tensor is defined, which is found to characterize material instability adequately. In order to investigate the validity of proposed stability criteria, three important constitutive models of hyperelasticity viz., St. Venant-Kirchho_, compressible Mooney-Rivlin and compressible Ogden models are investigated for stability.

In the literature, there are various material models proposed so as to model the constitutive behavior of hyperelastic materials for example, St. Venant-Kirchho_ model, Mooney-Rivlin model etc. The stability of such material models under various states of deformation is of important concern, and generally stability analysis is conducted in homogeneous states of deformation. Within hyperelasticity, instabilities can be broadly classified as geometrical and material types. Geometrical instabilities such as buckling, symmetric bifurcation etc. are of physical origin, and lead to multiple solutions at critical stretch. Material instability is a aw in the material model and leads to unphysical solutions at the onset. It is required that the constitutive model should be materially stable i.e., should not give unphysical results, and be able to predict correctly the onset of geometrical instabilities. Certain constitutive restrictions proposed in the literature are inadequate to characterize such instabilities. In the work, we propose stability criteria which will characterize geometrical as well as material instabilities. A new elasticity tensor is defined, which is found to characterize material instability adequately. In order to investigate the validity of proposed stability criteria, three important constitutive models of hyperelasticity viz., St. Venant-Kirchho_, compressible Mooney-Rivlin and compressible Ogden models are investigated for stability.

碩士
國立交通大學
機械工程系所
106
A finite element method formulated on the basis of the simple first order shear deformation theory is presented for the vibration analysis of composite plates. In the simple first order shear deformation theory, the displacement field of the composite plate comprises four displacement components, namely, two in-plane displacement components (uo, vo), vertical deflection due to bending (wb), and vertical deflection due to through thickness shear deformation (ws). The sum of wb and ws gives the total vertical deflection w. In the finite element formulation, the finite element consists of four nodes and, at each node, there are four nodal displacements (uoi, voi, wbi, wsi) and two rotations ( , ). The in-plane displacements and vertical deflection due to through thickness deformation within an element are obtained from the corresponding nodal displacements via the linear interpolation. On the hand, the vertical deflection due to bending within an element is obtained using the nonconforming shape functions for formulating Kirchhoff plates. The finite element method is used to study the free vibration of rectangular composite plates with different layer arrangements, regular boundary conditions, aspect ratios, and length-to-thickness ratios. It has been shown that the proposed finite element method can produce acceptable modal characteristics (natural frequencies and mode shapes) for the plates with length-to-thickness ratio larger than 20 when compared with the results available in the literature. As for thin composite plates with length-to-thickness ratio larger than 100, the proposed finite element method is also capable of producing good results without having the problem of shear locking. Finally, the proposed finite element method is used to predict the natural frequencies of several elastically restrained composite plates. The accuracy of the proposed finite element method has been verified by the experimental results. The effects of the properties of the elastic restraints on the modal characteristics of the composite plates are studied using the proposed finite element method.

Advances in the design and manufacturing technologies have greatly enhanced the utility of fiber reinforced composite materials in aircraft, helicopter and space- craft structural components. The special characteristics of composites such as high strength and stiffness, light-weight corrosion resistance make them suitable sub- stitute for metals/metallic alloys. However, composites are very sensitive to the anomalies induced during their fabrication and service life. Also, they are suscepti- ble to the impact and high frequency loading conditions because the epoxy matrix is at-least an order of magnitude weaker than the embedded reinforced carbon fibers. On the other hand, the carbon based matrix posses high electrical conductivity which is often undesirable. Subsequently, the metal matrix produces high brittleness. Var- ious forms of damage in composite laminates can be identified as indentation, fiber breakage, matrix cracking, fiber-matrix debonding and interply disbonding (delam- ination). Among all the damage modes mentioned above, delamination has been found to be serious for all cases of loading. They are caused by excessive interlaminar shear and normal stresses. The interlaminar stresses that arise in the case of composite materials due to the mismatch in the elastic constants across the plies. Delamination in composites reduce it’s tensile and compressive strengths by consid- erable margins. Hence the knowledge of these stresses is the most important aspect to be looked into. Basic theories like the Euler-Bernoulli’s theory and Timoshenko beam theory are based on many assumptions which poses limitation to determine these stresses accurately. Hence the determination of these interlaminar stresses accurately requires higher order theories to be considered. Most of the conventional methods of determination of the stresses are through the solutions, involving the trigonometric series, which are available only to small and simple problems. The most common method of solution is by Finite Element (FE) Method. There are only few elements existing in the literature and very few in the commercially available finite element software to determine the interlaminar stresses accurately in the composite laminates. Accuracy of finite element solution depends on the choice of functions to be used as interpolating polynomials for the field variable. In-appropriate choice will manifest in the form of delayed convergence. This delayed convergence and accuracy in predicting these stresses necessiates a formulation of elements with a completely new concept. The delayed convergence is sometimes attributed to the shear locking phenomena, which exist in most finite element formulation based on shear deformation theories. The present work aims in developing finite elements based on higher order theories, that alleviates the slow convergence and achieves the solutions at a faster rate without compromising on the accuracy. The accuracy primarily depends on the theory used to model the problem. Thus the basic theories (such as Elementary Beam theory and Timoshenko Beam theory) does not suffice the condition to accuratley determine the interlaminar stresses through the thickness, which is the primary cause for delamination in composites. Two different elements developed on the principle of super-convergence has been presented in this work. These elements are subjected to several numerical experiments and their performance is assessed by comparing the solutions with those available in literature. Spacecraft and aircraft structures are light in weight and are also lightly damped because of low internal damping of the material of construction. This increased exibility may allow large amplitude vibration, which might cause structural instability. In addition, they are susceptible to impact loads of very short duration, which excites many structural modes. Hence, structural dynamics and wave propagation study becomes a necessity. The wave based techniques have found appreciation in many real world problems such as in Structural Health Monitoring (SHM). Wave propagation problems are characterized by high frequency loads, that sets up stress waves to propagate through the medium. At high frequency, the wave lengths are small and from the finite element point of view, the element sizes should be of the same order as the wave lengths to prevent free edges of the element to act as a free boundary and start reflecting the stress waves. Also longer element size makes the mass distribution approximate. Hence for wave propagation problems, very large finite element mesh is an absolute necessity. However, the finite element problems size can be drastically reduced if we characterize the stiffness of the structure accurately. This can accelerate the convergence of the dynamic solution significantly. This can be acheived by the super-convergent formulation. Numerical results are presented to illustrate the efficiency of the new approach in both the cases of dynamic studies viz., the free vibration study and the wave propagation study. The thesis is organised into five chapters. A brief organization of the thesis is presented below, Chapter-1 gives the introduction on composite material and its constitutive law. The details of shear locking phenomena and the interlaminar stress distribution across the thickness is brought out and the present methods to avoid shear locking has been presented. Chapter-2 presents the different displacement based higher order shear deformation theories existing in the literature their advantages and limitations. Chapter-3 presents the formulation of a super-convergent finite element formulation, where the effect of lateral contraction is neglected. For this element static and free vibration studies are performed and the results are validated with the solution available in the open literature. Chapter-4 presents yet another super-convergent finite element formulation, wherein the higher order effects due to lateral contraction is included in the model. In addition to static and free vibration studies, wave propagation problems are solved to demonstrate its effectiveness. In all numerical examples, the super-convergent property is emphasized. Chapter-5 gives a brief summary of the total research work performed and presents further scope of research based on the current research.

Advances in the design and manufacturing technologies have greatly enhanced the utility of fiber reinforced composite materials in aircraft, helicopter and space- craft structural components. The special characteristics of composites such as high strength and stiffness, light-weight corrosion resistance make them suitable sub- stitute for metals/metallic alloys. However, composites are very sensitive to the anomalies induced during their fabrication and service life. Also, they are suscepti- ble to the impact and high frequency loading conditions because the epoxy matrix is at-least an order of magnitude weaker than the embedded reinforced carbon fibers. On the other hand, the carbon based matrix posses high electrical conductivity which is often undesirable. Subsequently, the metal matrix produces high brittleness. Var- ious forms of damage in composite laminates can be identified as indentation, fiber breakage, matrix cracking, fiber-matrix debonding and interply disbonding (delam- ination). Among all the damage modes mentioned above, delamination has been found to be serious for all cases of loading. They are caused by excessive interlaminar shear and normal stresses. The interlaminar stresses that arise in the case of composite materials due to the mismatch in the elastic constants across the plies. Delamination in composites reduce it’s tensile and compressive strengths by consid- erable margins. Hence the knowledge of these stresses is the most important aspect to be looked into. Basic theories like the Euler-Bernoulli’s theory and Timoshenko beam theory are based on many assumptions which poses limitation to determine these stresses accurately. Hence the determination of these interlaminar stresses accurately requires higher order theories to be considered. Most of the conventional methods of determination of the stresses are through the solutions, involving the trigonometric series, which are available only to small and simple problems. The most common method of solution is by Finite Element (FE) Method. There are only few elements existing in the literature and very few in the commercially available finite element software to determine the interlaminar stresses accurately in the composite laminates. Accuracy of finite element solution depends on the choice of functions to be used as interpolating polynomials for the field variable. In-appropriate choice will manifest in the form of delayed convergence. This delayed convergence and accuracy in predicting these stresses necessiates a formulation of elements with a completely new concept. The delayed convergence is sometimes attributed to the shear locking phenomena, which exist in most finite element formulation based on shear deformation theories. The present work aims in developing finite elements based on higher order theories, that alleviates the slow convergence and achieves the solutions at a faster rate without compromising on the accuracy. The accuracy primarily depends on the theory used to model the problem. Thus the basic theories (such as Elementary Beam theory and Timoshenko Beam theory) does not suffice the condition to accuratley determine the interlaminar stresses through the thickness, which is the primary cause for delamination in composites. Two different elements developed on the principle of super-convergence has been presented in this work. These elements are subjected to several numerical experiments and their performance is assessed by comparing the solutions with those available in literature. Spacecraft and aircraft structures are light in weight and are also lightly damped because of low internal damping of the material of construction. This increased exibility may allow large amplitude vibration, which might cause structural instability. In addition, they are susceptible to impact loads of very short duration, which excites many structural modes. Hence, structural dynamics and wave propagation study becomes a necessity. The wave based techniques have found appreciation in many real world problems such as in Structural Health Monitoring (SHM). Wave propagation problems are characterized by high frequency loads, that sets up stress waves to propagate through the medium. At high frequency, the wave lengths are small and from the finite element point of view, the element sizes should be of the same order as the wave lengths to prevent free edges of the element to act as a free boundary and start reflecting the stress waves. Also longer element size makes the mass distribution approximate. Hence for wave propagation problems, very large finite element mesh is an absolute necessity. However, the finite element problems size can be drastically reduced if we characterize the stiffness of the structure accurately. This can accelerate the convergence of the dynamic solution significantly. This can be acheived by the super-convergent formulation. Numerical results are presented to illustrate the efficiency of the new approach in both the cases of dynamic studies viz., the free vibration study and the wave propagation study. The thesis is organised into five chapters. A brief organization of the thesis is presented below, Chapter-1 gives the introduction on composite material and its constitutive law. The details of shear locking phenomena and the interlaminar stress distribution across the thickness is brought out and the present methods to avoid shear locking has been presented. Chapter-2 presents the different displacement based higher order shear deformation theories existing in the literature their advantages and limitations. Chapter-3 presents the formulation of a super-convergent finite element formulation, where the effect of lateral contraction is neglected. For this element static and free vibration studies are performed and the results are validated with the solution available in the open literature. Chapter-4 presents yet another super-convergent finite element formulation, wherein the higher order effects due to lateral contraction is included in the model. In addition to static and free vibration studies, wave propagation problems are solved to demonstrate its effectiveness. In all numerical examples, the super-convergent property is emphasized. Chapter-5 gives a brief summary of the total research work performed and presents further scope of research based on the current research.

Tesis (Doctor en Matemática)--Universidad Nacional de Córdoba. Facultad de Matemática, Astronomía y Física, 2009.
La tesis trata sobre las álgebras de Hopf semisimples y las categorías tensoriales asociadas. En la primera parte de la tesis mostramos dos familias de ejemplos de álgebras de Hopf semisimples que no son simples, pero que admiten deformaciones simples. La primera es una deformaciones por twist del grupo simétrico. La segunda familia es una deformación por twist del producto directo de dos grupos diédricos generalizados cuyo orden es el producto de dos números primos. Posteriormente, damos condiciones necesarias y suficientes para que una deformación por twist de un álgebra de grupo sea un álgebra de Hopf simple. En la segunda parte de la tesis, introducimos y estudiamos la noción de categoría tensorial fuertemente graduada sobre un grupo G. El resultado principal es la descripción de las categorías módulo sobre una categoría fuertemente graduada por un grupo G, como categorías módulo inducidas desde subcategorías tensoriales asociadas con los subgrupos de G.
César Neyit Galindo Martínez.

Wave propagation is a phenomenon resulting from high transient loadings where the duration of the load is in µ seconds range. In aerospace and space craft industries it is important to gain knowledge about the high frequency characteristics as it aids in structural health monitoring, wave transmission/attenuation for vibration and noise level reduction. The wave propagation problem can be approached by the conventional Finite Element Method(FEM); but at higher frequencies, the wavelengths being small, the size of the finite element is reduced to capture the response behavior accurately and thus increasing the number of equations to be solved, leading to high computational costs. On the other hand such problems are handled in the frequency domain using Fourier transforms and one such method is the Spectral Finite Element Method(SFEM). This method is introduced first by Doyle ,for isotropic case and later popularized in developing specific purpose elements for structural diagnostics for inhomogeneous materials, by Gopalakrishnan. The general approach in this method is that the partial differential wave equations are reduced to a set of ordinary differential equations(ODEs) by transforming these equations to another space(transformed domain, say Fourier domain). The reduced ODEs are usually solved exactly, the solution of which gives the dynamic shape functions. The interpolating functions used here are exact solution of the governing differential equations and hence, the exact elemental dynamic stiffness matrix is derived. Thus, in the absence of any discontinuities, one element is sufficient to model 1-D waveguide of any length. This elemental stiffness matrix can be assembled to obtain the global matrix as in FEM, but in the transformed space. Thus after obtaining the solution, the original domain responses are obtained using the inverse transform. Both the above mentioned manuscripts present the Fourier transform based spectral finite element (FSFE), which has the inherent aliasing problem that is persistent in the application of the Fourier series/Fourier transforms. This is alleviated by using an additional throw-off element and/or introducing slight damping in to the system. More recently wave let transform based spectral finite element(WSFE) has been formulated which alleviated the aliasing problem; but has a limitation in obtaining the frequency characteristics, like the group speeds are accurate only up-to certain fraction of the Nyquist(central frequency). Currently in this thesis Laplace transform based spectral finite elements(LSFE) are developed for sandwich members. The advantages and limitations of the use of different transforms in the spectral finite element framework is presented in detail in Chapter-1. Sandwich structures are used in the space craft industry due to higher stiffness to weight ratio. Many issues considered in the design and analysis of sandwich structures are discussed in the well known books(by Zenkert, Beitzer). Typically the main load bearing structures are modeled as beam sand plates. Plate structures with kh<1 is analysed based on the Kirch off plate theory/Classical Plate Theory(CPT) and when the bending wavelength is small compared to the plate thickness, the effect of shear deformation and rotary inertia needs to be included where, k is the wave number and h is the thickness of the plate. Many works regarding the wave propagation in sandwich structures has been published in the past literature for wave propagation in infinite sandwich structure and giving the complete description of dispersion relation with no restriction on frequency and wavelength. More recently exact analytical solution or simply supported sandwich plate has been derived. Also it is seen by comparison of dispersion curves obtained with exact (3D formulation of theory of elasticity) and simplified theories (2D formulation as generalization of Timoshenko theory) made on infinite domain and concluded that the simplified theory can be reliably used to assess the waveguide properties of sandwich plate in the frequency range of interest. In order to approach the problems with finite domain and their implementation in the use of general purpose code; finite degrees of freedom is enforced. The concept of displacement based theories provides the flexibility in assuming different kinematic deformations to approach these problems. Many of the displacement based theories incorporate the Equivalent Single Layer(ESL) approach and these can capture the global behavior with relative ease. Chapter-2 presents the Laplace spectral finite element for thick beams based on the First order Shear Deformation Theory (FSDT). Here the effect of different choices of the real part of the Laplace variable is demonstrated. It is shown that the real part of the Laplace variable acts as a numerical damping factor. The spectrum and dispersion relations are obtained and the use of these relations are demonstrated by an example. Here, for sandwich members based on FSDT, an appropriate choice of the correction factor ,which arises due to the inconsistency between the kinematic hypothesis and the desired accuracy is presented. Finally the response obtained by the use of the element is validated with experimental results. For high shock loading cases, the core flexibility induces local effects which are very predominant and this can lead to debonding of face sheets. The ESL theories mentioned above cannot capture these effects due to the computation of equivalent through the thickness section properties. Thus, higher order theories such as the layer-wise theories are required to capture the local behaviour. One such theory for sandwich panels is the Higher order Sandwich Plate theory (HSaPT). Here, the in-plane stress in the core has been neglected; but gives a good approximation for sandwich construction with soft cores. Including the axial inertial terms of the core will not yield constant shear stress distribution through the height of the core and hence more recently the Extended Higher order Sandwich Plate theory (EHSaPT) is proposed. The LSFE based on this theory has been formulated and is presented in Chapter-4. Detailed 3D orthotropic properties of typical sandwich construction is considered and the core compressibility effect of local behavior due to high shock loading is clearly brought out. As detailed local behavior is sought the degrees of freedom per element is high and the specific need for such theory as compared with the ESL theories is discussed. Chapter-4 presents the spectral finite element for plates based on FSDT. Here, multi-transform method is used to solve the partial differential equations of the plate. The effect of shear deformation is brought out in the spectrum and dispersion relations plots. Response results obtained by the formulated element is compared and validated with many different experimental results. Generally structures are built-up by connecting many different sub-structures. These connecting members, called joints play a very important role in the wave transmission/attenuation. Usually these joints are modeled as rigid joints; but in reality these are flexible and exhibits non-linear characteristics and offer high damping to the energy flow in the connected structures. Chapter-5 presents the attenuation and transmission of wave energy using the power flow approach for rigid joints for different configurations. Later, flexible spectral joint model is developed and the transmission/attenuation across the flexible joints is studied. The thesis ends with conclusion and highlighting futures cope based on the developments reported in this thesis.

Wave propagation is a phenomenon resulting from high transient loadings where the duration of the load is in µ seconds range. In aerospace and space craft industries it is important to gain knowledge about the high frequency characteristics as it aids in structural health monitoring, wave transmission/attenuation for vibration and noise level reduction. The wave propagation problem can be approached by the conventional Finite Element Method(FEM); but at higher frequencies, the wavelengths being small, the size of the finite element is reduced to capture the response behavior accurately and thus increasing the number of equations to be solved, leading to high computational costs. On the other hand such problems are handled in the frequency domain using Fourier transforms and one such method is the Spectral Finite Element Method(SFEM). This method is introduced first by Doyle ,for isotropic case and later popularized in developing specific purpose elements for structural diagnostics for inhomogeneous materials, by Gopalakrishnan. The general approach in this method is that the partial differential wave equations are reduced to a set of ordinary differential equations(ODEs) by transforming these equations to another space(transformed domain, say Fourier domain). The reduced ODEs are usually solved exactly, the solution of which gives the dynamic shape functions. The interpolating functions used here are exact solution of the governing differential equations and hence, the exact elemental dynamic stiffness matrix is derived. Thus, in the absence of any discontinuities, one element is sufficient to model 1-D waveguide of any length. This elemental stiffness matrix can be assembled to obtain the global matrix as in FEM, but in the transformed space. Thus after obtaining the solution, the original domain responses are obtained using the inverse transform. Both the above mentioned manuscripts present the Fourier transform based spectral finite element (FSFE), which has the inherent aliasing problem that is persistent in the application of the Fourier series/Fourier transforms. This is alleviated by using an additional throw-off element and/or introducing slight damping in to the system. More recently wave let transform based spectral finite element(WSFE) has been formulated which alleviated the aliasing problem; but has a limitation in obtaining the frequency characteristics, like the group speeds are accurate only up-to certain fraction of the Nyquist(central frequency). Currently in this thesis Laplace transform based spectral finite elements(LSFE) are developed for sandwich members. The advantages and limitations of the use of different transforms in the spectral finite element framework is presented in detail in Chapter-1. Sandwich structures are used in the space craft industry due to higher stiffness to weight ratio. Many issues considered in the design and analysis of sandwich structures are discussed in the well known books(by Zenkert, Beitzer). Typically the main load bearing structures are modeled as beam sand plates. Plate structures with kh<1 is analysed based on the Kirch off plate theory/Classical Plate Theory(CPT) and when the bending wavelength is small compared to the plate thickness, the effect of shear deformation and rotary inertia needs to be included where, k is the wave number and h is the thickness of the plate. Many works regarding the wave propagation in sandwich structures has been published in the past literature for wave propagation in infinite sandwich structure and giving the complete description of dispersion relation with no restriction on frequency and wavelength. More recently exact analytical solution or simply supported sandwich plate has been derived. Also it is seen by comparison of dispersion curves obtained with exact (3D formulation of theory of elasticity) and simplified theories (2D formulation as generalization of Timoshenko theory) made on infinite domain and concluded that the simplified theory can be reliably used to assess the waveguide properties of sandwich plate in the frequency range of interest. In order to approach the problems with finite domain and their implementation in the use of general purpose code; finite degrees of freedom is enforced. The concept of displacement based theories provides the flexibility in assuming different kinematic deformations to approach these problems. Many of the displacement based theories incorporate the Equivalent Single Layer(ESL) approach and these can capture the global behavior with relative ease. Chapter-2 presents the Laplace spectral finite element for thick beams based on the First order Shear Deformation Theory (FSDT). Here the effect of different choices of the real part of the Laplace variable is demonstrated. It is shown that the real part of the Laplace variable acts as a numerical damping factor. The spectrum and dispersion relations are obtained and the use of these relations are demonstrated by an example. Here, for sandwich members based on FSDT, an appropriate choice of the correction factor ,which arises due to the inconsistency between the kinematic hypothesis and the desired accuracy is presented. Finally the response obtained by the use of the element is validated with experimental results. For high shock loading cases, the core flexibility induces local effects which are very predominant and this can lead to debonding of face sheets. The ESL theories mentioned above cannot capture these effects due to the computation of equivalent through the thickness section properties. Thus, higher order theories such as the layer-wise theories are required to capture the local behaviour. One such theory for sandwich panels is the Higher order Sandwich Plate theory (HSaPT). Here, the in-plane stress in the core has been neglected; but gives a good approximation for sandwich construction with soft cores. Including the axial inertial terms of the core will not yield constant shear stress distribution through the height of the core and hence more recently the Extended Higher order Sandwich Plate theory (EHSaPT) is proposed. The LSFE based on this theory has been formulated and is presented in Chapter-4. Detailed 3D orthotropic properties of typical sandwich construction is considered and the core compressibility effect of local behavior due to high shock loading is clearly brought out. As detailed local behavior is sought the degrees of freedom per element is high and the specific need for such theory as compared with the ESL theories is discussed. Chapter-4 presents the spectral finite element for plates based on FSDT. Here, multi-transform method is used to solve the partial differential equations of the plate. The effect of shear deformation is brought out in the spectrum and dispersion relations plots. Response results obtained by the formulated element is compared and validated with many different experimental results. Generally structures are built-up by connecting many different sub-structures. These connecting members, called joints play a very important role in the wave transmission/attenuation. Usually these joints are modeled as rigid joints; but in reality these are flexible and exhibits non-linear characteristics and offer high damping to the energy flow in the connected structures. Chapter-5 presents the attenuation and transmission of wave energy using the power flow approach for rigid joints for different configurations. Later, flexible spectral joint model is developed and the transmission/attenuation across the flexible joints is studied. The thesis ends with conclusion and highlighting futures cope based on the developments reported in this thesis.

The classical Cauchy-Boltzmann theory of continuum mechanics requires that the dimension, over which macroscopic gradients occur, are much larger than characteristic length scales of the microstructure. For this reason, the classical continuum theory comes to its limits for very small specimens or if material degradation leads to a localisation of deformations into bands, whose width is determined by the microstructure itself. Deviations from the predictions of the classical theory of continuum mechanics are referred to as size effects. It is well-known, that generalised continuum theories can describe size effects in principle. Especially micromorphic theories gain increasing popularity due its favorable numerical implementation. However, the formulation of the additionally necessary constitutive equations is a problem. For linear-elastic behavior, the number of material parameters increases considerably compared to the classical theory. The experimental determination of these parameters is thus very difficult. For nonlinear and history-dependent processes, even the qualitative structure of the constitutive equations can hardly be assessed solely on base of phenomenological considerations. Homogenisation methods are a promising approach to solve this problem. The present thesis starts with a critical review on the classical theory of homogenisation and the approaches on micromorphic homogenisation which are available in literature. On this basis, a theory is developed for the homogenisation of a classical Cauchy-Boltzmann continuum at the microscale towards a micromorphic continuum at the macroscale. In particular, the micro-macro-relations are specified for all macroscopic kinetic and kinematic field quantities. On the microscale, the corresponding boundary-value problem is formulated, whereby kinematic, static or periodic boundary conditions can be used. No restrictions are imposed on the material behavior, i. e. it can be linear or nonlinear. The special cases of the micropolar theory (Cosserat theory), microstrain theory and microdilatational theorie are considered. The proposed homogenisation method is demonstrated for several examples. The simplest example is the uniaxial case, for which the exact solution can be specified. Furthermore, the micromorphic elastic properties of a porous, foam-like material are estimated in closed form by means of Ritz' method with a cubic ansatz. A comparison with partly available exact solutions and FEM solutions indicates a qualitative and quantitative agreement of sufficient accuracy. For the special cases of micropolar and microdilatational theory, the material parameters are specified in the established nomenclature from literature. By means of these material parameters the size effect of an elastic foam structure is investigated and compared with corresponding results from literature. Furthermore, micromorphic damage models for quasi-brittle and ductile failure are presented. Quasi-brittle damage is modelled by propagation of microcracks. For the ductile mechanism, Gurson's limit-load approach on the microscale is extended by microdilatational terms. A finite-element implementation shows, that the damage model exhibits h-convergence even in the softening regime and that it thus can describe localisation.:1 Introduction 2 Literature review: Micromorphic theory and strain-gradient theory 2.1 Variational approach 2.1.1 Cauchy-Boltzmann continuum 2.1.2 Second gradient theory / Strain gradient theory 2.1.3 Micromorphic theory 2.1.4 Method of virtual power 2.2 Homogenisation approaches 2.2.1 Classical theory of homogenisation 2.2.2 Strain-gradient theory by Gologanu, Kouznetsova et al. 2.2.3 Micromorphic theory by Eringen 2.2.4 Average field theory by Forest et al. 2.3 Scope of the present thesis 3 Homogenisation towards a micromorphic continuum 3.1 Thermodynamic considerations and generalized Hill-Mandel lemma 3.2 Surface operator and kinetic micro-macro relations 3.3 Kinematic micro-macro relations 3.4 Porous material 3.5 Kinematic and periodic boundary conditions 3.6 Special cases 3.6.1 Strain-gradient theory / Second gradient theory 3.6.2 Micropolar theory 3.6.3 Microstrain theory 3.6.4 Microdilatational theory 4 Elastic Behaviour 4.1 Uniaxial case 4.2 Upper bound estimates by Ritz' Method 4.3 Isotropic porous material 4.4 Micropolar theory 4.5 Microdilatational theory 4.6 Size effect in simple shear 5 Damage Models 5.1 Quasi-brittle damage 5.2 Microdilatational extension of Gurson’s model of ductile damage 5.2.1 Limit load analysis for rigid ideal-plastic material 5.2.2 Phenomenological extensions 5.2.3 FEM implementation 5.2.4 Example 6 Discussion
Die klassische Cauchy-Boltzmann-Kontinuumstheorie setzt voraus, dass die Abmessungen, über denen makroskopische Gradienten auftreten, sehr viele größer sind als charakteristische Längenskalen der Mikrostruktur. Aus diesem Grund stößt die klassische Kontinuumstheorie bei sehr kleinen Proben ebenso an ihre Grenzen wie bei Schädigungsvorgängen, bei denen die Deformationen in Bändern lokalisieren, deren Breite selbst von der Längenskalen der Mikrostruktur bestimmt wird. Abweichungen von Vorhersagen der klassischen Kontinuumstheorie werden als Größeneffekte bezeichnet. Es ist bekannt, dass generalisierte Kontinuumstheorien Größeneffekte prinzipiell beschreiben können. Insbesondere mikromorphe Theorien erfreuen sich auf Grund ihrer vergleichsweise einfachen numerischen Implementierung wachsender Beliebtheit. Ein großes Problem stellt dabei die Formulierung der zusätzlich notwendigen konstitutiven Gleichungen dar. Für linear-elastisches Verhalten steigt die Zahl der Materialparameter im Vergleich zur klassischen Theorie stark an, was deren experimentelle Bestimmung sehr schwierig macht. Bei nichtlinearen und lastgeschichtsabhängigen Prozessen lässt sich selbst die qualitative Struktur der konstitutiven Gleichungen ausschließlich auf Basis phänomenologischer Überlegungen kaum erschließen. Homogenisierungsverfahren stellen einen vielversprechenden Ansatz dar, um dieses Problem zu lösen. Die vorliegende Arbeit gibt zunächst einen kritischen Überblick über die klassische Theorie der Homogenisierung sowie die im Schrifttum verfügbaren Ansätze zur mikromorphen Homogenisierung. Auf dieser Basis wird eine Theorie zur Homogenisierung eines klassischen Cauchy-Boltzmann-Kontinuums auf Mikroebene zu einem mikromorphen Kontinuum auf der Makroebene entwickelt. Insbesondere werden Mikro-Makro-Relationen für alle makroskopischen kinetischen und kinematischen Feldgrößen angegebenen. Auf der Mikroebene wird das entsprechende Randwertproblem formuliert, wobei kinematische, statische oder periodische Randbedingungen verwendet werden können. Das Materialverhalten unterliegt keinen Einschränkungen, d. h., dass es sowohl linear als auch nichtlinear sein kann. Die Sonderfälle der mikropolaren Theorie (Cosserat-Theorie), Mikrodehnungstheorie und mikrodilatationalen Theorie werden erarbeitet. Das vorgeschlagene Homogenisierungsverfahren wird für eine Reihe von Beispielen demonstriert. Als einfachstes Beispiel dient der einachsige Fall, für den die exakte Lösung angegebenen werden kann. Weiterhin werden die mikromorphen, elastischen Eigenschaften eines porösen, schaumartigen Materials mittels des Ritz-Verfahrens mit einem kubischen Ansatz in geschlossener Form abgeschätzt. Ein Vergleich mit teilweise verfügbaren exakten Lösungen sowie FEM-Lösungen weist eine qualitative und quantitative Übereinstimmung hinreichender Genauigkeit aus. Für die Sonderfälle mikropolaren und mikrodilatationalen Theorien werden die Materialparameter in der im Schrifttum üblichen Nomenklatur angegebenen. Mittels dieser Materialparameter wird der Größeneffekt in einer elastischen Schaumstruktur untersucht und mit entsprechenden Ergebnissen aus dem Schrifttum verglichen. Desweiteren werden mikromorphe Schädigungsmodelle für quasi-sprödes und duktiles Versagen vorgestellt. Quasi-spröde Schädigung wird durch das Wachstum von Mikrorissen modelliert. Für den duktilen Mechanismus wird der Ansatz von Gurson einer Grenzlastanalyse auf Mikroebene um mikrodilatationale Terme erweitert. Eine Finite-Elemente-Implementierung zeigt, dass das Schädigungsmodell auch im Entfestigungsbereich h-Konvergenz aufweist und die Lokalisierung beschreiben kann.:1 Introduction 2 Literature review: Micromorphic theory and strain-gradient theory 2.1 Variational approach 2.1.1 Cauchy-Boltzmann continuum 2.1.2 Second gradient theory / Strain gradient theory 2.1.3 Micromorphic theory 2.1.4 Method of virtual power 2.2 Homogenisation approaches 2.2.1 Classical theory of homogenisation 2.2.2 Strain-gradient theory by Gologanu, Kouznetsova et al. 2.2.3 Micromorphic theory by Eringen 2.2.4 Average field theory by Forest et al. 2.3 Scope of the present thesis 3 Homogenisation towards a micromorphic continuum 3.1 Thermodynamic considerations and generalized Hill-Mandel lemma 3.2 Surface operator and kinetic micro-macro relations 3.3 Kinematic micro-macro relations 3.4 Porous material 3.5 Kinematic and periodic boundary conditions 3.6 Special cases 3.6.1 Strain-gradient theory / Second gradient theory 3.6.2 Micropolar theory 3.6.3 Microstrain theory 3.6.4 Microdilatational theory 4 Elastic Behaviour 4.1 Uniaxial case 4.2 Upper bound estimates by Ritz' Method 4.3 Isotropic porous material 4.4 Micropolar theory 4.5 Microdilatational theory 4.6 Size effect in simple shear 5 Damage Models 5.1 Quasi-brittle damage 5.2 Microdilatational extension of Gurson’s model of ductile damage 5.2.1 Limit load analysis for rigid ideal-plastic material 5.2.2 Phenomenological extensions 5.2.3 FEM implementation 5.2.4 Example 6 Discussion

This thesis is concerned with deformation theory of finite subschemes of smooth varieties. Of central interest are the smoothable subschemes (i.e., limits of smooth subschemes). We prove that all Gorenstein subschemes of degree up to 13 are smoothable. This result has immediate applications to finding equations of secant varieties.We also give a description of nonsmoothable Gorenstein subschemes of degree 14, together with an explicit condition for smoothability.We prove that being smoothable is a local property, that it does not depend on the embedding and it is invariant under a base field extension. The above results are equivalently stated in terms of the Hilbert scheme of points, which is the moduli space for this deformation problem.We extensively use the combinatorial framework of Macaulay's inverse systems. We enrich it with a pro-algebraic group action and use this to reprove and extend recent classification results by Elias and Rossi. We provide a relative version of this framework and use it to give a local description of the universal family over the Hilbert scheme of points.We shortly discuss history of Hilbert schemes of points and provide a list of open questions.
Tematem rozprawy są deformacje podschematów skończonych gładkich rozmaitości. Koncentrujemy się na schematach wygładzalnych (tj. będących granicami schematów gładkich). Dowodzimy, że wszystkie schematy Gorensteina stopnia co najwyżej 13 są wygładzalne. To twierdzenie ma bezpośrednie zastosowanie dla znajdowania równań rozmaitości siecznych.Podajemy również opis niewygładzalnych schematów Gorensteina stopnia 14 wraz z warunkiem na wygładzalność.Dowodzimy, że wygładzalność jest własnością lokalną oraz że nie zależy ona od zanurzenia i że jest niezmienna przy rozszerzeniu ciała bazowego. Powyższe wyniki można równoważnie sformułować w terminach schematu Hilberta punktów, który jest przestrzenią moduli dla tego problemu deformacyjnego.Naszym podstawowym narzędziem kombintorycznym są systemy odwrotne Macaulaya. Wzbogacamy tę teorię o działanie pro-algebraicznej grupy i stosujemy ją do uogólnienia wyników klasyfikacyjnych Eliasa i Rossi. Podajemy relatywną wersję systemów odwrotnych Macaulaya i, używając jej, lokalny opis rodziny uniwersalnej nad schematem Hilberta punktów.Krótko dyskutujemy historię badań nad schematami Hilberta punktów i podajemy listę otwartych problemów.

The aim of this work is to develop and validate a continuum model for the simulation of the thermomechanical response of a shape memory polymer (SMP). Rather than integral type viscoelastic model, the approach here is based on the idea of two inter-penetrating networks, one which is permanent and the other which is transient together with rate equations for the time evolution of the transient network. We find that the activation stress for network breakage and formation of the material controls the gross features of the response of the model, and exhibits a "thermal Bauschinger effect". The model developed here is similar to a thermoviscoelastic model, and is developed with an eye towards ease of numerical solutions to boundary value problems. The primary hypothesis of this model is that the hysteresis of temperature dependent activation-stress plays a lead role in controlling its main response features. Validation of this hypothesis is carried out for the uniaxial response from the experimental data available in the literature for two different SMP samples: shape memory polyurethane and Veriflex, to show the control of the evolution of the temperature sensitive activation stress on the response. We extend the validated 1D model to a three dimensional small strain continuum SMP model and carry out a systematic parameter optimization method for the identification of the activation stress coefficients, with different weights given to different features of the response to match the parameters with experimental data. A comprehensive parametric study is carried out, that varies each of the model material and loading parameters, and observes their effect on design-relevant response characteristics of the model undergoing a thermomechanical cycle. We develop "response charts" for the response characteristics: shape fixity, shape recovery and maximum stress rise during cooling, to give the designer an idea of how the simultaneous variation of two of the most influential material parameters changes a specific response parameter. To exemplify the efficacy of the model in practical applications, a thermoviscoelastic extension of a beam theory model will be developed. This SMP beam theory will account for activation stress governed inelastic response of a SMP beam. An example of a three point bend test is simulated using the beam theory model. The numerical solution is implemented by using an operator split technique that utilizes an elastic predictor and dissipative corrector. This algorithm is validated by using a three-point bending experiment for three different material cases: elastic, plastic and thermoplastic response. Time step convergence and mesh density convergence studies are carried out for the thermoviscoelastic FEM model. We implement and study this model for a SMP beam undergoing three-point bending strain recovery, stress recovery and cyclic thermomechanical loading. Finally we develop a thermodynamically consistent finite continuum model to simulate the thermomechanical response of SMPs. The SMP is modeled as an isotropic viscoplastic material where thermal changes govern the evolution of the activation stress of the material. The response of the SMP in a thermomechanical cycle is modeled as a combination of a rubbery and a glassy element in series. Using these assumptions, we propose a specific form for the Helmholtz potential and the rate of dissipation. We use the technique of upper triangular decomposition for developing the constitutive equations of the finite strain SMP model. The resulting model is implemented in an ODE solver in MATLAB, and solved for a simple shear problem. We study the response of the SMP model for shear deformation as well as cyclic shear deformation at different initial temperatures. Finally, we implement the thermomechanical cycle under shear deformations and study the behavior of the model.

Looking at severe plastic deformation experiments, it seems that crystalline materials at yield behave as a special kind of anisotropic, compressible, highly viscous fluid. In the presented approach the plastic behaviour of crystalline solids is treated as a highly viscous material flow through an adjustable crystal lattice.The main purpose of this dissertation is to investigate model of a plastic flow of a highly viscous fluid that describes the ultra-fine structure formation induced by severe plastic deformation. The idea behind is to apply and further develop methods known from fluid mechanics to modelling crystal solids.We first provide the thermodynamic description of a fluid like - Eulerian model of crystal plasticity. The model derivation is based on the application of the Gibbs potential to obtain a rate type stress strain constitutive relation. The result is compared to the approaches of traditional plasticity.The second part of this thesis is devoted to the mathematical analysis of a simplified problems originating from the visco-elastic model derived to describe flows of crystal plastic materials. Even after simplifications, neglecting the plastic effects, the problem seems to be difficult to establish local-in-time existence of a smooth solution and the existence of a weak solution. We propose a regularisation of the stress evolution equation and prove the global-in-time existence of a weak solutions, in the case of two dimensional bounded domain, for general initial data. The results, obtained by the Galerkin method, hold true for the periodic case. We then give a thorough description of the numerical methods used in the dissertation. Necessary tools in order to discretize the fluid model of crystal plasticity both in space and time, such that: triangulation of a domain, proper finite dimensional spaces and time discretization schemes, are recalled. The discrete system resulting from the weak formulation of the considered system is solved by means of the Newton’s method.Further the results of numerical simulations are reported. Performed simulations aim to justify that the presented approach is capable of capturing large strains in typical experimental settings. We give a detailed description of performed numerical simulations. Two different deformation settings were considered: uniaxial compression and channel extrusion. In the case of a simple compression we deal with a free boundary problem. For this reason we employ the Arbitrary Lagrangian Eulerian (ALE) approach in the sense that we use the Eulerian formulation on a moving mesh which captures the free boundary. As a test example we analyse a plane strain compression in a channel die and an uniaxial compression of a pillar-shaped sample. The second case focuses on 2-turn equal channel angular pressing (ECAP) in two dimensions. We take advantage of the Eulerian formulation and obtain high strains in a single 2-turn ECAP pass.