La modélisation du trafic routier dans un réseau dense et de grande étendue nécessite un grand nombre de données, si bien qu'une modélisation par arcs est impossible en pratique. Pour simplifier le problème, une idée est d'agréger les tronçons du réseau en un continuum sur lequel le trafic routier s'écoule comme un fluide surfacique. Cette modélisation est qualifiée de bidimensionnelle. Même si la structure géométrique détaillée du réseau est perdue, une telle modélisation évite la description très fine du trafic sur un réseau dans lequel les points de mesure ne sont pas en nombre suffisant pour permettre une évaluation exhaustive de l'état du trafic. Une série de travaux commencée dans les années 1980 a permis de dégager quelques concepts importants. Cependant, ces travaux n'ont pas résolu les problèmes de modélisation fondamentaux : comment déduire et modéliser des comportements globaux à partir de comportements locaux (flux sur un axe, interactions aux intersections) ? Deux modèles bidimensionnels de trafic sont développés. Le premier modèle est statique. Le trafic s'écoule dans des directions de propagation privilégiées (orthogonales). Le modèle prend en compte l'équilibre entre l'offre de déplacement du réseau et une demande élastique de déplacement des usagers. Les principales sorties sont constituées en chaque point et pour chaque destination par les débits directionnels et les coûts de déplacement. Le deuxième modèle est dynamique. L'écoulement du trafic est décrit au niveau de cellules élémentaires du réseau dans lesquelles on définit les notions d'offre et de demande. À partir d'une loi comportementale obtenue sur un réseau discret, on écrit l'équation dynamique de conservation du trafic routier en tout point d'un réseau anisotrope

La modélisation du trafic routier dans un réseau dense et de grande étendue nécessite un grand nombre de données, si bien qu'une modélisation par arcs est impossible en pratique. Pour simplifier le problème, une idée est d'agréger les tronçons du réseau en un continuum sur lequel le trafic routier s'écoule comme un fluide surfacique. Cette modélisation est qualifiée de bidimensionnelle. Même si la structure géométrique détaillée du réseau est perdue, une telle modélisation évite la description très fine du trafic sur un réseau dans lequel les points de mesure ne sont pas en nombre suffisant pour permettre une évaluation exhaustive de l'état du trafic. Une série de travaux commencée dans les années 1980 a permis de dégager quelques concepts importants. Cependant, ces travaux n'ont pas résolu les problèmes de modélisation fondamentaux : comment déduire et modéliser des comportements globaux à partir de comportements locaux (flux sur un axe, interactions aux intersections) ? Deux modèles bidimensionnels de trafic sont développés. Le premier modèle est statique. Le trafic s'écoule dans des directions de propagation privilégiées (orthogonales). Le modèle prend en compte l'équilibre entre l'offre de déplacement du réseau et une demande élastique de déplacement des usagers. Les principales sorties sont constituées en chaque point et pour chaque destination par les débits directionnels et les coûts de déplacement. Le deuxième modèle est dynamique. L'écoulement du trafic est décrit au niveau de cellules élémentaires du réseau dans lesquelles on définit les notions d'offre et de demande. À partir d'une loi comportementale obtenue sur un réseau discret, on écrit l'équation dynamique de conservation du trafic routier en tout point d'un réseau anisotrope / Traffic road modelling in the dense network of a wide area needs a large amount of data. It renders a keen modelling unmanageable in practice. To simplify the problem, an idea is to aggregate the network links as continuous medium where traffic road flows as a fluid on a surface. This modelling is called two-dimensional. Even if the detailed geometric structure of the network is lost, such a modelling avoids the traffic keen description on a network where measure points are not numerous enough to allow an exhaustive evaluation of traffic state. A series of articles started in the 80's have highlighted relevant concepts. Nevertheless, these previous works have not solved fundamental modelling issues: how to deduce and model global behaviours basing on local behaviours (flow on an axe, interactions at intersections). Two models are developed. The first model is static. Traffic road flows on privileged directions (orthogonal). The model takes into account the balance between the network trip supply and the users' elastic trip demand. The principal outputs are, for each point of the network and each destination, the directional flows and the trip costs. The second model is dynamical. The description of traffic flows stands at the level of elementary traffic cells, where concepts of supply and demand are defined. With a behaviour law deduced from a discrete network, we establish the conservation dynamic equation of road traffic at each point of an anisotropic network

Treten im Beton mehraxiale Spannungszustände auf, führen diese gegenüber einer einaxialen Beanspruchung zu einer signifikanten Änderung des Materialverhaltens. Neben einer festigkeitssteigernden bzw. -abmindernden Wirkung ergeben sich ebenfalls große Unterschiede im Spannungs-Dehnungs-Verhalten. Zur effizienten Konzipierung von Betonstrukturen unter komplexen Beanspruchungszuständen ist daher die Kenntnis des veränderten Materialverhaltens notwendig. Zur experimentellen Bestimmung des Spannungs-Dehnungs-Verhaltens eines Ultrahochleistungsbetons mit einer einaxialen Druckfestigkeit von über 170 N/mm² wurden mehraxiale Belastungsversuche an würfelförmigen Probekörpern durchgeführt. Die Untersuchung umfasste insgesamt 35 zwei- und dreiaxiale Spannungsverhältnisse unter proportionaler Laststeigerung mit vorrangiger Betrachtung von Zug-Druck-Druck-Beanspruchungen. Für die Einleitung der Zugbeanspruchungen in die Prüfkörper wurde eine neue Methode entwickelt, bei der mittels einbetonierter Schrauben die Belastung auf den Beton übertragen wird. Die Bestimmung des Verformungsverhaltens erfolgte im Inneren der Probekörper mit sechs tetraederförmig angeordneten Faser-Bragg-Gittern. Die somit direkt gemessenen Dehnungen ermöglichen die nachträgliche Berechnung der Komponenten des Dehnungstensors des Bezugskoordinatensystems. Für den untersuchten Ultrahochleitsungsbeton fallen die auf die einaxiale Druckfestigkeit bezogenen mehraxialen Festigkeitswerte mit zunehmendem hydrostatischen Druckspannungsanteil der Beanspruchung geringer aus als bei Normalbetonen. Weiterhin weist das Verformungsverhalten eine größere Sprödigkeit gegenüber Normalbetonen auf, so dass auch unter dreiaxialen Druckspannungszuständen die Probekörper schlagartig versagen. Aus den gemessenen Spannungs-Dehnungs-Linien werden neben den maximalen Festigkeiten die Festigkeitswerte an der Elastizitätsgrenze, der Affinitätsgrenze sowie beim Volumenminimum der Probekörper bestimmt. Zur Approximation dieser charakteristischen Werte wurde eine Grenzflächenbeschreibung entwickelt und an den Versuchsergebnissen kalibriert. Des Weiteren erfolgte die Zusammenstellung einer Datenbank mit in der Literatur verfügbaren mehraxialen maximalen Festigkeitswerten von Betonen mit einaxialen Druckfestigkeiten von 10 N/mm² bis 180 N/mm² und die Kalibrierung des entwickelten Modells zur Grenzflächenbeschreibung in Abhängigkeit der einaxialen Druckfestigkeit. Die bei der Kalibrierung der Grenzfläche für einzelne Betonfestigkeitsklassen bestimmten Freiwerte hängen dabei stark von den vorliegenden Versuchsdaten und speziell vom Wertebereich der hydrostatischen Spannungsanteile der maximalen Beanspruchungen ab. Die Approximation des Spannungs-Dehnungs-Verhaltens der mehraxial beanspruchten Probekörper erfolgt mittels eines schädigungsbasierten Materialgesetzes. Hierbei wird für den anfänglich isotropen Beton zum einen eine lastinduzierte isotrope Schädigung und zum anderen eine lastinduzierte orthotrope Schädigung angenommen, die von den auftretenden Hauptdehnungen abhängig ist. Mit dem entwickelten Materialgesetz werden sehr gute Übereinstimmungen mit den gemessenen Spannungs-Dehnungs-Linien erreicht, so dass sich ebenfalls eine gute Vorhersage der maximalen Festigkeitswerte ergibt. / Concrete under multiaxial stress states shows significant changes of the material behaviour compared to uniaxial loading. Besides strength increasing and decreasing effects, also great differences in the stress-strain behaviour occur. In order to design concrete structures efficiently concerning complex stress states, the knowledge about the modified material behaviour is necessary. To determine experimentally the stress-strain behaviour of an ultra high performance concrete with a uniaxial compressive strength of about 170 N/mm², multiaxial loading tests on cubic-shaped specimens were carried out. Altogether, the investigation contained 35 biaxial and triaxial stress ratios under proportionally increasing load with primarily tension-compression-compression loadings. Applying the tensile load on the specimen, a new method was developed, which uses screws embedded in the concrete to transfer the loading. The deformations were measured by using six tetrahedron-shaped arranged Fibre Bragg Gratings inside the concrete specimen. Subsequently, with the directly measured strains the components of the strain tensor of the reference coordinate system could be determined. For the investigated ultra high performance concrete the increase of the multiaxial strength, referring to the uniaxial compressive strength, decreases compared to normal strength concrete with the increasing hydrostatic stress component of the load. Moreover, the deformation behaviour shows an increased brittleness compared to normal strength concrete, so that even under triaxial compressive stress states the specimens fail abruptly. Besides the ultimate strength, from the measured stress-strain curves the strength at the proportional limit, at the limit of affinity as well as at the minimum volume of the specimen is determined. To approximate these characteristic values, a description of a hypersurface is developed and calibrated with the test results. Furthermore, a database with multiaxial ultimate strength values of concretes with uniaxial compressive strengths between 10 N/mm² to 180 N/mm² available from literature was compiled and a calibration of the developed hypersurface model depending on the uniaxial compressive strength was carried out. Thereby, the obtained values of arbitrary parameters of individual concrete strength classes depend severely on the available test results, especially on the range of values of the hydrostatic stress component of the ultimate strength. The approximation of the stress-strain behaviour of the multiaxial loaded specimens is carried out by means of a damage-based material law. For this purpose, concerning the initially isotropic concrete, a load-induced isotropic and orthotropic damage depending on the principle strains is assumed. With the developed material law, very good accordance with the measured stress-strain curves could be achieved, so that also results in a good approximation of the ultimate concrete strength.

La technique des goujons collés dans les structures bois allie performances mécaniques et esthétisme. Elle répond au besoin de conservation du bâti pour le domaine de la rénovation et au besoin d’assemblages de plus en plus rigides et résistants en construction neuve avec l’avènement des structures bois de grandes ampleurs et l’arrivée de nouveaux produits tels que le bois lamellé croisé (CLT). De nombreuses études de caractérisation mécanique ont été menées au cours de ces 35 dernières années. L’enquête menée auprès des professionnels et experts en vue de l’introduction de cette technique dans l’Eurocode 5 a néanmoins révélé des attentes et lacunes particulières vis-à-vis des connaissances sur le comportement des assemblages multi-tiges, le choix des adhésifs, le comportement à long terme du collage et les techniques de contrôle qualité. Les travaux de cette thèse, menés en collaboration entre l’université de Bordeaux et l’Université Laval (Québec, Canada), ont pour objectif d’apporter des réponses sur le comportement mécanique des adhésifs, et assemblages sous sollicitation thermique et d’initier les développements sur le comportement à long terme et les assemblages multi-tiges. Les principaux résultats sont issus de travaux expérimentaux et numériques valorisés sous la forme d’articles scientifiques. Les essais mécaniques dynamiques menés sur une colle polyuréthane (PUR) et une résine époxy (EPX) entre 30°C et 120°C mettent en avant des différences importantes de raideur et de dégradations entre les deux adhésifs. Exposés à la chaleur, les assemblages sont eux aussi affectés avec des pertes de raideur et de résistance dès 40°C, en amont de la température de transition vitreuse et de la dégradation des adhésifs. Amené à suivre le régime de température extérieur, l’assemblage doit par sa conception être isolé afin de prévenir les risques aux états limites de service (ELS) et aux états limites ultimes (ELU). Très peu de données sont disponibles dans la littérature sur le comportement à long terme des assemblages goujons collés. L’étude du fluage ajoute une dimension temporelle aux problématiques de raideur. La base de données expérimentales collectée dans ces travaux repose sur 12 campagnes d’essais de 1 à 2 mois en conditions climatiques régulées en température et humidité relative, ou variables (service classe 1). Chargement et climat influent sur le glissement, mais aussi sur la durée de vie. Si des disparités sont présentes entre EPX et PUR, les études menées à 20°C et 50°C révèlent à nouveau l’importance de prendre en considération la température, notamment aux ELS. L’étude en régime variable illustre de son côté la sensibilité du système aux variations de température et d’humidité. L’étude de la répartition des efforts en fonction de la raideur des goujons est choisie comme fil conducteur pour investiguer le passage vers des éléments multi-tiges. Un modèle 3D mono-tige aux éléments finis est tout d’abord développé pour prendre en compte l’orthotropie du bois et regarder l’influence du centre de moelle jusqu’à présent négligée. Par la suite ce modèle est étendu aux configurations multi-tiges afin de simuler différentes configurations d’essais. Un prototype d’essai de caractérisation mécanique des assemblages multi-tiges est finalement présenté. / Glued-in rods have successfully been used for connections or reinforcement of timber structures. With the development of tall and large timber buildings and new products such as Cross Laminated Timber (CLT), there is an increasing need for connections that provide high stiffness and strength. Timber connections using glued in rods have a general aesthetic appeal, and take advantage of the structural adhesives that provide a high stiffness and load capacity. During the last 30 years, multiple applications have been developed for renovation and new construction. In parallel, numerous investigations have characterized the mechanical performance of these connections, but harmonized design rules are not available. A lack of knowledge on the multiple-rod connections, on the choice of adhesives, on long-term effects and on the control of quality is identified by a recent study and must be investigated before a new submission to Eurocode 5. This thesis is performed in the framework of collaboration between Université de Bordeaux (France) and Université Laval (Canada). Objectives are focused on the mechanical performance of adhesives and glued-in rod connections under elevated temperature, on the development of creep tests and the investigations of multiple-rod connections. This study combines experimental and finite element modelling results, which are presented in the form of scientific articles. A Dynamic Mechanical Analysis (DMA) on one polyurethane (PUR) and one epoxy (EPX) adhesives followed by static tensile tests on the connections with small-diameter steel glued-in rods have been conducted at different temperatures. High differences in stiffness are observed between the two adhesives. Glued-in rod connections with the EPX and PUR adhesives demonstrate significant losses of stiffness and resistance beyond 40°C, before the glass transition temperature of the adhesive. Following the outside temperature, connections must be insulated to prevent excessive slip and risk of failure at the Serviceability Limits State (SLS) and at the Ultimate Limit States (ULS). Few data on the long-term effects are available in the literature. In this work, 12 campaigns of creep tests, adding time effects to previous work, have been performed. Glued-in rod connections were tested during one to two months in a conditioning room where temperature and relative humidity were controlled or in a room with a variable climate representative of service class 1. Disparities were observed between EXP and PUR. Creep tests conducted at 20°C and 50°C revealed an important role of the temperature, particularly for SLS. Creep tests in variable climate illustrated the sensibility of connections to the temperature and humidity variations. The stiffness and stress distribution in multiple-rod connections were investigated. First, a 3D model was developed to observe the influence of the wood orthotropy on the stiffness and stress distribution in a connection with a single rod. Then, the model was extended to multiple-rod connections to simulate different loads and boundary conditions. Finally, a prototype of a connection with multiple glued-in rods was tested.

Treten im Beton mehraxiale Spannungszustände auf, führen diese gegenüber einer einaxialen Beanspruchung zu einer signifikanten Änderung des Materialverhaltens. Neben einer festigkeitssteigernden bzw. -abmindernden Wirkung ergeben sich ebenfalls große Unterschiede im Spannungs-Dehnungs-Verhalten. Zur effizienten Konzipierung von Betonstrukturen unter komplexen Beanspruchungszuständen ist daher die Kenntnis des veränderten Materialverhaltens notwendig. Zur experimentellen Bestimmung des Spannungs-Dehnungs-Verhaltens eines Ultrahochleistungsbetons mit einer einaxialen Druckfestigkeit von über 170 N/mm² wurden mehraxiale Belastungsversuche an würfelförmigen Probekörpern durchgeführt. Die Untersuchung umfasste insgesamt 35 zwei- und dreiaxiale Spannungsverhältnisse unter proportionaler Laststeigerung mit vorrangiger Betrachtung von Zug-Druck-Druck-Beanspruchungen. Für die Einleitung der Zugbeanspruchungen in die Prüfkörper wurde eine neue Methode entwickelt, bei der mittels einbetonierter Schrauben die Belastung auf den Beton übertragen wird. Die Bestimmung des Verformungsverhaltens erfolgte im Inneren der Probekörper mit sechs tetraederförmig angeordneten Faser-Bragg-Gittern. Die somit direkt gemessenen Dehnungen ermöglichen die nachträgliche Berechnung der Komponenten des Dehnungstensors des Bezugskoordinatensystems. Für den untersuchten Ultrahochleitsungsbeton fallen die auf die einaxiale Druckfestigkeit bezogenen mehraxialen Festigkeitswerte mit zunehmendem hydrostatischen Druckspannungsanteil der Beanspruchung geringer aus als bei Normalbetonen. Weiterhin weist das Verformungsverhalten eine größere Sprödigkeit gegenüber Normalbetonen auf, so dass auch unter dreiaxialen Druckspannungszuständen die Probekörper schlagartig versagen. Aus den gemessenen Spannungs-Dehnungs-Linien werden neben den maximalen Festigkeiten die Festigkeitswerte an der Elastizitätsgrenze, der Affinitätsgrenze sowie beim Volumenminimum der Probekörper bestimmt. Zur Approximation dieser charakteristischen Werte wurde eine Grenzflächenbeschreibung entwickelt und an den Versuchsergebnissen kalibriert. Des Weiteren erfolgte die Zusammenstellung einer Datenbank mit in der Literatur verfügbaren mehraxialen maximalen Festigkeitswerten von Betonen mit einaxialen Druckfestigkeiten von 10 N/mm² bis 180 N/mm² und die Kalibrierung des entwickelten Modells zur Grenzflächenbeschreibung in Abhängigkeit der einaxialen Druckfestigkeit. Die bei der Kalibrierung der Grenzfläche für einzelne Betonfestigkeitsklassen bestimmten Freiwerte hängen dabei stark von den vorliegenden Versuchsdaten und speziell vom Wertebereich der hydrostatischen Spannungsanteile der maximalen Beanspruchungen ab. Die Approximation des Spannungs-Dehnungs-Verhaltens der mehraxial beanspruchten Probekörper erfolgt mittels eines schädigungsbasierten Materialgesetzes. Hierbei wird für den anfänglich isotropen Beton zum einen eine lastinduzierte isotrope Schädigung und zum anderen eine lastinduzierte orthotrope Schädigung angenommen, die von den auftretenden Hauptdehnungen abhängig ist. Mit dem entwickelten Materialgesetz werden sehr gute Übereinstimmungen mit den gemessenen Spannungs-Dehnungs-Linien erreicht, so dass sich ebenfalls eine gute Vorhersage der maximalen Festigkeitswerte ergibt. / Concrete under multiaxial stress states shows significant changes of the material behaviour compared to uniaxial loading. Besides strength increasing and decreasing effects, also great differences in the stress-strain behaviour occur. In order to design concrete structures efficiently concerning complex stress states, the knowledge about the modified material behaviour is necessary. To determine experimentally the stress-strain behaviour of an ultra high performance concrete with a uniaxial compressive strength of about 170 N/mm², multiaxial loading tests on cubic-shaped specimens were carried out. Altogether, the investigation contained 35 biaxial and triaxial stress ratios under proportionally increasing load with primarily tension-compression-compression loadings. Applying the tensile load on the specimen, a new method was developed, which uses screws embedded in the concrete to transfer the loading. The deformations were measured by using six tetrahedron-shaped arranged Fibre Bragg Gratings inside the concrete specimen. Subsequently, with the directly measured strains the components of the strain tensor of the reference coordinate system could be determined. For the investigated ultra high performance concrete the increase of the multiaxial strength, referring to the uniaxial compressive strength, decreases compared to normal strength concrete with the increasing hydrostatic stress component of the load. Moreover, the deformation behaviour shows an increased brittleness compared to normal strength concrete, so that even under triaxial compressive stress states the specimens fail abruptly. Besides the ultimate strength, from the measured stress-strain curves the strength at the proportional limit, at the limit of affinity as well as at the minimum volume of the specimen is determined. To approximate these characteristic values, a description of a hypersurface is developed and calibrated with the test results. Furthermore, a database with multiaxial ultimate strength values of concretes with uniaxial compressive strengths between 10 N/mm² to 180 N/mm² available from literature was compiled and a calibration of the developed hypersurface model depending on the uniaxial compressive strength was carried out. Thereby, the obtained values of arbitrary parameters of individual concrete strength classes depend severely on the available test results, especially on the range of values of the hydrostatic stress component of the ultimate strength. The approximation of the stress-strain behaviour of the multiaxial loaded specimens is carried out by means of a damage-based material law. For this purpose, concerning the initially isotropic concrete, a load-induced isotropic and orthotropic damage depending on the principle strains is assumed. With the developed material law, very good accordance with the measured stress-strain curves could be achieved, so that also results in a good approximation of the ultimate concrete strength.

L’objectif de ce travail est de développer des outils numériques utilisables dans la détermination de manière exacte des propriétés modales des structures sandwichs viscoélastiques composites au vue de la conception des structures sandwichs viscoélastiques légères mais à haut pouvoir amortissant. Pour cela nous avons tout d’abord développé un outil générique implémenté en Matlab pour la détermination des propriétés modales en vibration libre des plaques sandwichs viscoélastiques dont les faces sont en stratifié de plusieurs couches orientées dans diverses directions. L’intérêt de cet outil, basé sur une formulation éléments finis, réside dans sa capacité à prendre en compte l’anisotropie des couches composites, la non linéarité matérielle de la couche viscoélastique traduit par diverses lois viscoélastiques dépendant de la fréquence ainsi que diverses conditions aux limites. La résolution du problème aux valeurs propres non linéaires complexes se fait par le couplage entre la technique d’homotopie, la méthode asymptotique numérique et la différentiation automatique. Ensuite pour permettre une étude continue des effets d’un paramètre de modélisation sur les propriétés modales des sandwichs viscoélastiques, nous avons proposé une méthode générique de résolution de problème résiduel non linéaire aux valeurs propres complexes possédant en plus de la dépendance en fréquence introduite par la couche viscoélastique du coeur, la dépendance du paramètre de modélisation qui décrit un intervalle d’étude bien spécifique. Cette résolution est basée sur la méthode asymptotique numérique, la différentiation automatique, la technique d’homotopie et la continuation et prend en compte diverses lois viscoélastiques. Nous proposons après cela, deux formulations distinctes pour étudier les effets, sur les propriétés amortissantes, de deux paramètres de modélisation qui sont importants dans la conception de sandwichs viscoélastiques à haut pouvoir amortissement. Le premier est l’orientation des fibres des composites dans la référence du sandwich et le second est l’épaisseur des couches qui lorsqu’elles sont bien définies permettent d’obtenir non seulement des structures sandwichs à haut pouvoir amortissant mais très légères. Les équations fortement non linéaires aux valeurs propres complexes obtenues dans ces formulations sont résolues par la nouvelle méthode de résolution d’équation résiduelle développée. Des comparaisons avec des résultats discrets sont faites ainsi que les temps de calcul pour montrer non seulement l’utilité de ces deux formulations mais également celle de la méthode de résolution d’équations résiduelles non linéaires aux valeurs propres complexes à double dépendance / Modeling and design of composite viscoelastic structures with high damping powerThe aim of this thesis is to develop numerical tools to determine accurately damping properties of composite sandwich structures for the design of lightweight viscoelastic sandwichs structures with high damping power. In a first step, we developed a generic tool implemented in Matlab for determining damping properties in free vibration of viscoelastic sandwich plates with laminate faces composed of multilayers. The advantage of this tool, which is based on a finite element formulation, is its ability to take into account the anisotropy of composite layers, the material non-linearity of the viscoelastic core induiced by the frequency-dependent viscoelastic laws and various boundary conditions . The nonlinear complex eigenvalues problem is solved by coupling homotopy technic, asymptotic numerical method and automatic differentiation. Then for the continuous study of a modeling parameter on damping properties of viscoelastic sandwichs, we proposed a generic method to solve the nonlinear residual complex eigenvalues problem which has in addition to the frequency dependence introduced by the viscoelastic core, a modeling parameter dependence that describes a very specific study interval. This resolution is based on asymptotic numerical method, automatic differentiation, homotopy technique and continuation technic and takes into account various viscoelastic laws. We propose after that, two separate formulations to study effects on the damping properties according to two modeling parameters which are important in the design of high viscoelastic sandwichs with high damping power. The first is laminate fibers orientation in the sandwich reference and the second is layers thickness which when they are well defined allow to obtain not only sandwich structures with high damping power but also very light. The highly nonlinear complex eigenvalues problems obtained in these formulations are solved by the new method of resolution of eigenvalue residual problem with two nonlinearity developed before. Comparisons with discrete results and computation time are made to show the usefulness of these two formulations and of the new method of solving nonlinear complex eigenvalues residual problem of two dependances