Isogeometric analysis is aimed to mitigate the gap between Computer-Aided Design (CAD) and analysis by using a unified geometric representation. Thanks to the exact geometry representation and high smoothness of adopted basis functions, isogeometric analysis demonstrated excellent mathematical properties and successfully addressed a variety of problems. In particular, it allows to solve higher order Partial Differential Equations (PDEs) directly omitting the usage of mixed approaches. Unfortunately, complex CAD geometries are often constituted by multiple Non-Uniform Rational B-Splines (NURBS) patches and cannot be directly applied for finite element analysis.parIn this work, we presents a dual mortaring framework to couple adjacent patches for higher order PDEs. The development of this formulation is initiated over the simplest 4th order problem-biharmonic problem. In order to speed up the construction and preserve the sparsity of the coupled problem, we derive a dual mortar compatible C1 constraint and utilize the Bezier dual basis to discretize the Lagrange multipler spaces. We prove that this approach leads to a well-posed discrete problem and specify requirements to achieve optimal convergence. After identifying the cause of sub-optimality of Bezier dual basis, we develop an enrichment procedure to endow Bezier dual basis with adequate polynomial reproduction ability. The enrichment process is quadrature-free and independent of the mesh size. Hence, there is no need to take care of the conditioning. In addition, the built-in vertex modification yields compatible basis functions for multi-patch coupling.To extend the dual mortar approach to couple Kirchhoff-Love shell, we develop a dual mortar compatible constraint for Kirchhoff-Love shell based on the Rodrigues' rotation formula. This constraint provides a unified formulation for both smooth couplings and kinks. The enriched Bezier dual basis preserves the sparsity of the coupled Kirchhoff-Love shell formulation and yields accurate results for several benchmark problems.Like the dual mortaring formulation, locking problem can also be derived from the mixed formulation. Hence, we explore the potential of Bezier dual basis in alleviating transverse shear locking in Timoshenko beams and volumetric locking in nearly compressible linear elasticity. Interpreting the well-known B projection in two different ways we develop two formulations for locking problems in beams and nearly incompressible elastic solids. One formulation leads to a sparse symmetric symmetric system and the other leads to a sparse non-symmetric system.

Cette thèse présente des travaux effectués dans le cadre de l'optimisation de forme de pièces mécaniques, sous critère dynamique, par approche isogéométrique. Pour réaliser une telle optimisation nous mettons en place dans un premier temps les éléments coque au travers des formulations Kirchhoff-Love puis Reissner-Minlin. Nous présentons une méthode permettant d'atteindre les vecteurs normaux aux fibres dans ces formulations au travers de l'utilisation d'une grille mixte de fonctions de base interpolantes, traditionnellement utilisées en éléments finis, et de fonction non interpolantes issues de la description isogéométrique des coques. Par la suite, nous détaillons une méthode pour le couplage de patch puis nous mettons en place la méthode de synthèse modale classique dans le cadre de structures en dynamique décrites par des éléments isogéometriques. Ce travail établit une base pour l'optimisation de forme sous critères dynamique de telles structures. Enfin, nous développons une méthode d'optimisation de forme basée sur le calcul du gradient de la fonction objectif envisagée. La sensibilité de conception est extraite de l'analyse de sensibilité au niveau même du maillage du modèle, qui est obtenue par l'analyse discrète de sensibilité. Des exemples d'application permettent de montrer la pertinence et l'exactitude des approches proposées. / Isogeometric method is a promising method in bridging the gap between the computer aided design and computer aided analysis. No information is lost when transferring the design model to the analysis model. It is a great advantage over the traditional finite element method, where the analysis model is only an approximation of the design model. It is advantageous for structural optimization, the optimal structure obtained will be a design model. In this thesis, the research is focused on the fast three dimensional free shape optimization with isogeometric shell elements. The related research, the development of isogeometric shell elements, the patch coupling in isogeometric analysis, the modal synthesis with isogeometric elements are also studied. We proposed a series of mixed grid Reissner-Minlin shell formulations. It adopts both the interpolatory basis functions, which are from the traditional FEM, and the non-interpolatory basis functions, which are from IGA, to approximate the unknown elds. It gives a natural way to define the fiber vectors in IGA Reissner-Mindlin shell formulations, where the non-interpolatory nature of IGA basis functions causes complexity. It is also advantageous for applying the rotational boundary conditions. A modified reduce quadrature scheme was also proposed to improve the quadrature eficiency, at the same time, relieve the locking in the shell formulations. We gave a method for patch coupling in isogeometric analysis. It is used to connect the adjacent patches. The classical modal synthesis method, the fixed interface Craig-Bampton method, is also used as well as the isogeometric Kirchhoff-Love shell elements. The key problem is also the connection between adjacent patches. The modal synthesis method can largely reduce the time costs in analysis concerning structural dynamics. This part of work lays a foundation for the fast shape optimization of built-up structures, where the design variables are only relevant to certain substructures. We developed a fast shape optimization framework for three dimensional thin wall structure design. The thin wall structure is modelled with isogeometric Kirchhoff-Love shell elements. The analytical sensitivity analysis is the key focus, since the gradient base optimization is normally more fast. There are two models in most optimization problem, the design model and the analysis model. The design variables are defined in the design model, however the analytical sensitivity is normally obtained from the analysis model. Although it is possible to use the same model in analysis and design under isogeomeric framework, it might give either a highly distorted optimum structure or a unreliable structural response. We developed a sensitivity mapping scheme to resolve this problem. The design sensitivity is extracted from the analysis model mesh level sensitivity, which is obtained by the discrete analytical sensitivity analysis. It provides exibility for the design variable definition. The correctness of structure response is also ensured. The modal synthesis method is also used to further improve the optimization eficiency for the built-up structure optimization concerning structural dynamics criteria.

The prediction of fracture in thin-walled structures is decisive for a wide range of applications. Modeling methods such as the phase-field method usually consider cracks to be constant over the thickness which, especially in load cases involving bending, is an imperfect approximation. In this contribution, fracture phenomena along the thickness direction of structural elements (plates or shells) are addressed with a phase-field modeling approach. For this purpose, a new, so called “mixed-dimensional” model is introduced, which combines structural elements representing the displacement field in the two-dimensional shell midsurface with continuum elements describing a crack phase-field in the three-dimensional solid space. The proposed model uses two separate finite element discretizations, where the transfer of variables between the coupled twoand three-dimensional fields is performed at the integration points which in turn need to have corresponding geometric locations. The governing equations of the proposed mixed-dimensional model are deduced in a consistent manner from a total energy functional with them also being compared to existing standard models. The resulting model has the advantage of a reduced computational effort due to the structural elements while still being able to accurately model arbitrary through-thickness crack evolutions as well as partly along the thickness broken shells due to the continuum elements. Amongst others, the higher accuracy aswell as the numerical efficiency of the proposed model are tested and validated by comparing simulation results of the new model to those obtained by standard models using numerous representative examples.