No description available.

No description available.

Dans la production d'électricité, un des leviers centraux pour réduire les détériorations et les pertes causées par l'humidité dans les turbines à vapeur est l'étude des films liquides. Ces films minces, sont créés par la déposition de gouttes et sont fortement cisaillés. Des gouttes peuvent ensuite être arrachées du film. A l'heure actuelle, aucun modèle complet et valide n'existe pour décrire ce phénomène. Un modèle 2D à formulation intégrale associé à des lois de fermetures a été dérivé pour représenter ce film. Comparé aux équations classiques de Saint-Venant, le modèle prend en compte davantage d'effets : le transfert de masse, l'impact des gouttes, le cisaillement à la surface libre, la tension de surface, le gradient de pression et la rotation. Une analyse des propriétés du modèle (hyperbolicité, entropie, conservativité, analyse de stabilité linéaire, invariance par translation et par rotation) est réalisée pour juger de la pertinence du modèle. Un nouveau code 2D est implémenté dans un module de développement libre du code EDF Code Saturne et une méthode de volumes finis pour un maillage non-structure a été développée. La vérification du code est ensuite effectuée avec des solutions analytiques dont un problème de Riemann. Le modèle, qui dégénère en modèle classique de Saint-Venant pour le cas d'un film tombant sur un plan inclinée, est validé par l'expérience de Liu and Gollub, 1994, PoF et comparé à des modèles de références (Ruyer-Quil and Manneville, 2000, EPJ-B et Lavalle, 2014, PhD thesis). Un autre cas d'étude met en scène un film cisaillé en condition basse-pression de turbine à vapeur et, est validé par l'expérience de Hammitt et al., 1981, I. Enfin, le code film est couplé aux données 3D du champ de vapeur autour d'un stator d'une turbine basse-pression du parc EDF, issues de Blondel, 2014, PhD thesis. Cette application industrielle montre la faisabilité d'une simulation d'un film en condition réelle du turbine à vapeur. / In the electricity production, one central key to reduce damages and losses due to wetness in steam turbines is the study of liquid films. These thin films are created by the deposition of droplets and are highly sheared. This film may then be atomized into coarse water. At the moment, no comprehensive and validated model exists to describe this phenomenon. A 2D model based on a integral formulation associated with closure laws is developed to represent this film. Compared to classical Shallow-Water equation, the model takes into account additional effect : mass transfer, droplet impact, shearing at the free surface, surface tension, pressure gradient and the rotation. The model properties (hyperbolicity, entropy, conservativity, linear stability, Galilean invariance and rotational invariance) has been analyzed to judge the pertinence of the model. A new 2D code is implemented in a free module of the code EDF Code Saturne and a finite volume method for unstructured mesh has been developed. The verification of the code is then carried out with analytical solutions including a Riemann problem. The model, which degenerates into classical Shallow-Water equations for the case of a falling liquid film on a inclined plane, is validated by the experiment of Liu and Gollub, 1994, PoF and compared to reference models (Ruyer-Quil and Manneville, 2000, EPJ-B et Lavalle, 2014, PhD thesis). Another study depicts a sheared film under low-pressure steam turbine conditions and is validated by the experiment of Hammitt et al., 1981, FiI. Lastly, the code film is coupled to 3D steam data around a fixed blade of a BP100 turbine, from Blondel, 2014, PhD thesis. This industrial application shows the feasibility of liquid film's simulation in real steam turbine condition.

Cette thèse s’inscrit dans le cadre d’une collaboration entre le Laboratoire Génie Civil et Bâtiment (LGCB)/Laboratoire de Tribologie et Dynamique des Systèmes (LTDS) et l’entreprise Arkema, l’entreprise routière Malet, l’Agence de l’Environnement et de la Maîtrise de l’Énergie (ADEME). L’objectif de la thèse est l’étude du comportement thermomécanique des enrobés bitumineux, et spécifiquement les enrobés tièdes avec et sans agrégat d’enrobé (RAP, Reclaimed Asphalt Pavement), dans le but de prévoir leur durée de vie et d’optimiser leur dimensionnement dans un contexte de développement durable.Sept types d’enrobés bitumineux sont étudiés. Les enrobés diffèrent par leur procédé de fabrication, le pourcentage de RAP, les additifs et le pourcentage d’additif utilisé. Deux types de fabrications sont étudiées : à chaud (Hot Mix –HM) qui sert de référence et tièdes (Warm Mix –WM). Trois pourcentages de RAP sont considérés : 0%, 30% et 50%. Deux types d’additifs fournis par ARKEMA sont utilisés, ETIMA (E) et ETIMA-BIO (B), avec deux pourcentages, 0% et 0.4%. L’utilisation des enrobés tièdes contribue à réduire considérablement la consommation d’énergie, l’impact sanitaire et l’impact environnemental. D’autre part, l’ajout des RAP dans les enrobés permet d’augmenter le taux de recyclage des matières non renouvelables. Ces innovations participent à la prise en compte des enjeux du développement durable et de la préservation de l’environnement.Quatre types d’essais ont été réalisés. Ce sont : l’essai de module complexe, l’essai de fatigue, l’essai de retrait thermique empêché (TSRST) et l’essai de propagation d’ondes. Dans le domaine du comportement viscoélastique linéaire, les essais de module complexe sont réalisés sur une large gamme de températures (de -25°C à 45°C) et de fréquences (de 0.03Hz à 10Hz). Chaque test est dupliqué. Les résultats sont modélisés à l’aide du modèle analogique 2S2P1D qui a été développé au laboratoire LGCB de l’ENTPE. Ce modèle permet de modéliser correctement le comportement tridimensionnel viscoélastique linéaire des enrobés. Ensuite, le comportement en fatigue est étudié grâce à des essais de fatigue. Pour chaque matériau, quatre ou cinq éprouvettes sont testées. Quatre critères sont utilisés pour évaluer la durée de vie des matériaux. Puis, le comportement à basse température est caractérisé à l’aide de l’essai de retrait thermique empêché (TSRST). Pour chaque matériau, trois éprouvettes sont testées. Enfin, des essais de propagation d’ondes sont réalisés. La méthode de détermination du temps de vol des ondes « P » et des ondes « S » ainsi que la méthode Impact Résonance sont utilisées. Ces essais fournissent des mesures non destructives et faciles à réaliser. A partir de ces méthodes, on peut calculer les valeurs de modules et de coefficients de Poisson des matériaux.Soulignons que les essais réalisés permettent d’obtenir certain des paramètres utilisés pour le dimensionnement des chaussées. Un résultat essentiel de l’étude est la bonne performance des enrobés tièdes contenant des agrégats d’enrobés (RAP). / This PhD dissertation is part of the collaboration between the French Agency for Environment and Energy Management (ADEME), the companies Arkema and Malet and the Ecole Nationale de Travaux Publics de l’Etat (ENTPE). The aim of this dissertation is to study the thermo-mechanical behavior of asphalt mixes, specifically warm mix asphalts with and without Reclaimed Asphalt Pavement (RAP) in order to predict their fatigue life and to optimize their design by targeting sustainable development. Seven types of asphalt mixes are studied. They differ in their mixing method, RAP content, additives used and their content. Two types of mixing methods are studied: Hot Mix Asphalt (HMA), which serves as a reference, and Warm Mix Asphalt (WMA). Three RAP contents are considered: 0%, 30% and 50%. Two types of additives provided by Arkema are used : ETIMA (E) and ETIMA-BIO (B) with two contents for each of them (0% and 0.4%). The use of warm mix asphalt helps significantly reducing both energy consumption and the impact on health and environment. These innovations are significant contributions with respect to issues of sustainable development and environmental preservation. Four types of laboratory tests were performed: complex modulus tests, fatigue tests, Thermal Stress Restrained Specimen Tests (TSRST) and wave propagation tests. In the Linear ViscoElastic (LVE) behavior domain, complex modulus tests were performed on a wide range of temperatures (from -25°C to 45°C) and frequencies (from 0.03 Hz to 10 Hz). For each material, two specimens were tested. Values of complex modulus and complex Poisson’s ratio were obtained and modeled using the 2S2P1D analogies model, developed in the Laboratoire Génie Civil et Bâtiement (LGCB) of ENTPE. It was observed that this model can correctly simulate the LVE behavior of asphalt mixes. Fatigue behavior was studied by analyzing fatigue tests results. For each material, four or five specimens were tested. Four different failure criteria were used in order to evaluate fatigue life of materials. Low temperature behavior of materials was characterized using TSRST. For each material, three replicates were performed. Finally, wave propagation tests were carried out. The method to determine the time of flight of « P » and « S » waves and the impact resonance method are presented. These tests provide a non-destructive method to characterize materials, which easy to perform. Using those methods, values of complex modulus and Poisson's ratio of tested materials can be estimated. It is to be underlined that the tests performed in this PhD allow to obtaining some parameters which are required for asphalt pavement design. An important result obtained in this study is that warm mix asphalt combined with RAP and additives (ETIMA ou ETIMA-BIO) could perform as well as hot mix asphalt.

De nombreux travaux dans la littérature se sont concentrés sur la modélisation vibro-acoustique de coques cylindriques raidies immergées, du fait des nombreuses applications industrielles, en particulier dans le domaine aéronautique ou naval. Cependant, peu d'entre elles prennent en compte des structures internes non-axisymétriques telles que des supports moteurs, des planchers ou des carlingages, qui peuvent avoir une influence importante sur le comportement vibro-acoustique du système. C'est pourquoi une méthode de sous-structuration baptisée CTF est présentée dans cette thèse. Elle est développée dans le cas général de deux structures minces couplées le long d'une ligne. Un ensemble de fonctions orthonormées, baptisées fonctions de condensation, est défini afin d'approximer les forces et déplacements à la jonction entre les sous-systèmes. Des fonctions de transfert condensées sont définies pour chaque sous-système découplé. L'utilisation du principe de superposition, de l'équilibre des forces et de la continuité des déplacements permet de déduire le comportement des sous-systèmes couplés. La méthode est d'abord développée et validée dans le cas de plaques, puis ensuite appliquée au cas d'une coque cylindrique raidie immergée couplée à des structures internes non-axisymétriques. Le système est dans ce cas décomposé en 3 familles de sous-systèmes : la coque cylindrique immergée décrite par une méthode semi-analytique basée sur la résolution des équations de Flügge dans le domaine des nombres d’onde, les structures internes axisymétriques (raidisseurs, cloisons) décrites par éléments finis axisymétriques et les structures non-axisymétriques décrites pas des modèles éléments finis. La méthode CTF est appliquée à différents cas tests afin de montrer l'influence des structures internes non-axisymétriques sur le comportement vibro-acoustique d'une coque cylindrique pour différents types d'excitations pertinents dans le domaine naval : une force ponctuelle, une onde plane acoustique et un champ de pression aléatoire (tel qu'un champ acoustique diffus ou une couche limite turbulente). / Many works can be found in the literature concerning the vibroacoustic modelling of submerged stiffened cylindrical shells, because of high interest in the industrial domain, in particular for aeronautical or naval applications. However, only a few of them take into account non-axisymmetric internal frames, as for instance engine foundations or floor partitions, that can play a role on the vibroacoustic behavior of the system. That is why a substructuring approach called the Condensed Transfer Function (CTF) approach is proposed in the first part of this thesis. The aim is to take advantage of both analytical models and element-based models, in order to be able to deal with the geometrical complexity, and to calculate at higher frequencies than with element-based methods only. The substructuring method is developed in the general case of thin mechanical structures coupled along curves. A set of orthonormal functions called condensation functions, which depend on the curvilinear abscissa along the coupling line, is considered. This set is then used as a basis for approximating and decomposing the displacements and the applied forces at the line junctions. Thanks to the definition and calculation of condensed transfer functions for each uncoupled subsystem and by using the superposition principle for passive linear systems, the behavior of the coupled subsystems can be obtained. The method is first developed and validated for plates and convergence criteria are defined in relation with the size of the basis of condensation functions. The CTF method is then applied to the case of a submerged stiffened cylindrical shell with non-axisymmetric internal frames. The system is partitioned in 3 types of subsystems: the submerged shell, the axisymmetric frames (stiffeners, bulkheads) and the non-axisymmetric frames. The submerged shell is described by a semi-analytical method based on the Flügge equations in the spectral domain. The axisymmetric frames are described by axisymmetric Finite Element models and the non-axisymmetric frames by Finite Element models. The CTF method is applied to different test cases in order to highlight the influence of non-axisymmetric internal frames on the vibroacoustic behavior of a submerged stiffened cylindrical shell, for different excitations particularly relevant in the naval domain: a point force, an acoustic plane wave, and a random pressure field (such as a diffuse sound field or a turbulent boundary layer for instance).

Dans cette thèse, nous nous intéressons à la modélisation mathématique des hétérogénéités de longueurs caractéristiques beaucoup plus petites que la longueur d'ondes. La thèse consiste en deux parties. La partie théorique est dédiée à l'obtention d'un développement asymptotique raccordé: la solution est décrite à l'aide d'un développement de champ proche au voisinage de l'obstacle et par un développement de champ lointain hors de ce voisinage. Le développement de champ lointain met en jeu des solutions singulières de l'équation des ondes tandis que le champ proche lui est régi par un modèle quasi-statique. Ces deux développements sont alors raccordés dans une zone intermédiaire dite de raccord. Nous obtenons alors des estimations d'erreurs permettant de rendre rigoureux ce développement asymptotique formel. La deuxième partie est numérique. Elle décrit à la fois la méthode de Galerkine discontinue, une méthode de raffinement de maillage espace-temps et propose une discrétisation des modèles asymptotiques obtenues précédemment. Elle est illustrée par un certain nombre de tests numériques. / In this thesis, we focus our attention on the modeling of heterogeneities which are smaller than the wavelength. The document is decomposed into two parts : a theoretical one and a numerical one. In the first part, we derive a matched asymptotic expansion composed of a far-field expansion and a near-field expansion. The terms of the far-field expansion are singular solutions of the wave equation whereas the terms of the near-field expansion satisfy quasistatic problems. These expansions are matched in an intermediate region. We justify mathematically this theory by proving error estimates. In the second part, we describe the Discontinuous Galerkin method, a local time stepping method and the implementation of the matched asymptotic method. Numerical simulations illustrate these results.

L’Emission Acoustique (EA) est une technique de contrôle non-destructif consistant en la mesure et l’interprétation de la signature acoustique de mécanismes d’endommagement. Dans l’approche conventionnelle (approche phénoménologique), l’interprétation des données issues des mesures par EA s’appuie sur des corrélations empiriques entre des caractéristiques de la source (le mécanisme d’endommagement) et le signal mesuré. Les modifications dues à la chaine d’acquisition de l’EA sont donc ignorées. Or, la propagation dans le matériau, la mesure par le capteur et le traitement par le système d’acquisition modifient la forme du signal et l’information qu’il transporte. Cela rend difficile l’identification de la source, et la comparaison des résultats issus d’essais effectués dans des conditions différentes. Une partie de la réponse à ces problèmes réside dans la prise en compte des étapes de transformation du signal d’EA. C’est l’objectif de l’approche quantitative de l’EA. Cette approche repose sur l’utilisation de techniques de modélisation pour évaluer l’impact de chaque étape de transformation sur le signal. Le premier volet de cette étude porte sur la validation des techniques utilisées pour simuler les étapes de transformation du signal d’EA. La méthode des éléments finis (MEF) est utilisée pour simuler la propagation du signal au sein du matériau. L’effet du capteur est quant à lui simulé par sa fonction de sensibilité, mesurée par la méthode de réciprocité, et utilisée comme fonction de transfert. Le second volet porte sur l’utilisation de ces techniques pour évaluer l’impact, sur le signal d’EA, des paramètres (position, temps de montée, amplitude) d’une source simple modélisée par des dipôles de force. Trois géométries d’éprouvette sont étudiées : une première éprouvette assimilable à une plaque, une seconde assimilable à une poutre de section rectangulaire et une dernière dont les dimensions sont intermédiaires entre une plaque et une poutre. Les résultats obtenus montrent que les signaux se propagent au sein des éprouvettes suivant des modes bien définis. Ces modes de propagation sont dépendants de la géométrie de l’éprouvette. Chaque source sollicite les modes différemment. Ainsi leur étude permet de réunir des informations sur la source afin de l’identifier. Par ailleurs, cette étude a permis de mettre en évidence des descripteurs pertinents pour l’identification des sources d’EA. Les descripteurs, corrélés entre eux, permettent une nette séparation des signaux en classes en fonction de la source. Ces résultats, obtenus en surface matériau, ne prennent pas en compte l’effet du capteur. Lorsque celui-ci est pris en compte, la signature modale des sources est modifiée ainsi que la valeur des descripteurs calculés. Cela conduit à un recouvrement des classes de signaux rendant plus difficile l’identification des sources. / Acoustic emission (AE) is a non-destructive testing technique consisting in measuring and interpreting the acoustic signature of damage mechanisms. In a conventional treatment approach (phenomenological approach), the interpretation of data measured by AE is based on empirical correlations between the source (the damage mechanism) parameters and the measured signal. Therefore, the modifications due to the acquisition chain of acoustic emission are ignored. Yet, propagation of the waves in the material, measures made by the sensor and signal treatments made by the acquisition system modify the signal and the information it carries. As a consequence, identification of the source and comparison with results from other tests made in different conditions are difficult. To find a solution to these problems, one can take into account the different steps of the acquisition chain. This is the goal of Quantitative Acoustic Emission (QAE). This approach is based on modelling techniques to evaluate the impact of each step of the acquisition chain on the AE signal. The first part of this study concerns the experimental validation of the modelling techniques that were used in simulating the steps of the acquisition chain. The Finite Element Method (FEM) is used in simulating the signal propagation inside the material. The sensor effect on the signal is simulated by its sensitivity function, measured by the reciprocity method and used as a transfer function. The second part deals with using these techniques to evaluate the impact of simple AE sources on the AE signal. These simple sources are considered as a point source and modelled by dipole forces. Three tensile specimen geometries are studied: a first specimen that can be compared to a plate, a second specimen that can be compared to a beam and a third specimen of intermediate dimensions. The obtained results show the mechanical waves propagate inside the specimens as modes. These modes depend on the specimen geometry. Each source excites the wave propagation modes in a different manner. Consequently, studying the excited modes, one can gather useful information on the AE source to identify it. In addition, this study highlighted relevant signal parameters to identify AE sources. The correlation of these parameters allows segregating the signals as a function of the source. These results obtained at the material surface don’t take into account the sensor modifications on the signal. The sensor modifies the modal signature of the sources as well as the value of the calculated parameters. This leads to more difficulties in identifying the AE sources.

Das Radar mit synthetischer Apertur (Synthetic Aperture Radar - SAR) liefert eine quasi-fotographische Abbildung der beleuchteten Bodenoberfläche mit zusätzlichen Informationen, die von der gesendeten und empfangenen Polarisation der Wellen abhängig sind. Eine nützliche Anwendung der polarimetrischen SAR-Daten liegt bei der Klassifizierung der Bodenstruktur anhand der polarimetrischen Streuungseigenschaften. In diesem Zusammenhang beschäftigt sich die vorliegende Arbeit mit der Entwicklung und Untersuchung neuer polarimetrischen Fokussierungsfunktion für die SAR-Datenverarbeitung mit Hilfe der polarimetrischen Rückstreuungseigenschaft, die zu einer alternativen quantitativen Auswertung der polarimerischen SAR-Daten führen kann. Die physikalische Optik Approximation wird für die numerische Berechnung der rückgestreuten elektrischen Felder der kanonischen Ziele unter SAR-Geometrie unter Berücksichtigung der Polarisationslage verwendet. Aus den rückgestreuten elektrischen Felder werden die polarimetrischen Radarrückstreuquerschnitte berechnet. Ein SAR-Simulator wird zur Datenverarbeitung der E-SAR des DLR entwickelt. Der Ansatz des polarimetrischen Radarrückstreuquerschnittes ermöglicht die approximierte numerische Berechnung der Rückstreuungseigenschaften der kanonischen Ziele sowohl im kopolaren als auch im kreuzpolaren Polarisationsbetrieb. Bei der SAR-Datenverarbeitung werden die Rohdatensätze durch die Referenzfunktion eines Punktzieles in der Entfernungsrichtung verarbeitet. Bei der Azimutkompression werden die vier Referenzfunktionen, das heißt die Referenzfunktion eines Punktzieles, die polarimetrische Fokussierungsfunktion einer flachen Platte, die polarimetrische Fokussierungsfunktion eines Zweifach-Reflektors und die polarimetrische Fokussierungsfunktion eines Dreifach-Reflektors, eingesetzt. Die qunatitativen Auswertung der SAR-Daten werden anhand des Pauli-Zerlegungstheorems, der differentiellen Reflektivität und des linearen Depolarisationsverhältnises durchgeführt.

The need to protect marine structures from the high-intensity impulsive loads created by underwater explosions has stimulated renewed interest in the mechanical response of sandwich structures. The objective of this combined numerical and experimental study is to analyze the dynamic response of composite sandwich structures and develop material-structure-property relations and design criteria for improving the blast-resistance of marine structures. Configurations analyzed include polymer foam core structures with planar geometries. A novel experimental facility to generate high-intensity underwater impulsive loads and carry out in-situ measurements of dynamic deformations in marine structures is developed. Experiments are supported by fully dynamic finite-element simulations which account for the effects of fluid-structure interaction, and the constitutive and damage response of E-glass/polyester composites and PVC foams. Results indicate that the core-density has a significant influence on dynamic deformations and failure modes. Polymeric foams experience considerable rate-effects and exhibit extensive shear cracking and collapse under high-magnitude multi-axial underwater impulsive loads. In structures with identical masses, low-density foam cores consistently outperform high-density foam cores, undergoing lesser deflections and transmitting smaller impulses. Calculations reveal a significant difference between the response of air-backed and water-backed structures. Water-backed structures undergo much greater damage and consequently need to absorb a much larger amount of energy than air-backed structures. The impulses transmitted through water-backed structures have significant implications for structural design. The thickness of the facesheets is varied under the conditions of constant material properties and core dimensions. The results reveal an optimal thickness of the facesheets which maximizes energy absorption in the core and minimizes the overall deflection of the structure.

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.