La versatilité de LiNbO3 est bien connue dans le domaine de l’électronique acoustique (filtres SAW), de l’optique non linéaire et de l’optique électronique. Alors que les applications de LiNbO3 congruent sont répandues dans le domaine de l’électronique acoustique et de l’optique où les puissances ne sont pas très élevées, le seuil de dommage optique étant relativement bas, les applications sont limitées dans le domaine des fortes puissances lumineuses. Le seuil de dommage optique de LiNbO3 stœchiométrique est plus haut et rend nécessaire l’utilisation d’une méthode de croissance plus difficile (flux). Les difficultés de la méthode de croissance par flux justifie la recherche en vue de développer des méthodes alternatives ou d’améliorer les méthodes élaborées. Dans ma thèse j’ai présenté des solvants alternatifs á la croissance de monocristaux de LiNbO3 à partir du système ternaire K2O-Li2O-Nb2O5. J’ai prouvé qu’il est possible de faire croître des cristaux stœchiométriques par la méthode HTTSSG à partir des systèmes ternaires X2O-Li2O-Nb2O5 (X = Rb, Cs). Les diagrammes de phase tracés et précisés au cours de ce travail sont importants et donnent la possibilité de composer le flux à partir duquel il est possible d’obtenir le cristal de son choix. Comme pour les matériaux typiquement non stœchiométrique, les paramètres physiques de LiNbO3 sont très sensibles aux moindres variations de la composition (± 0.01 mol%). La détermination précise de la composition est donc très importante. Dans ma thèse, j’ai élaboré une nouvelle méthode de haute précision basée sur le rapport des deux principaux pics vibrationnels OH- dans l’infrarouge. Cette méthode a une précision qui satisfait les exigences de qualification des cristaux de LiNbO3 / The versatility of the application of LiNbO3 is well-known in the practical fields of acousto-electronics and non-linear optics and electro-optics. The optical damage threshold of stoichiometric LiNbO3 is higher so it is more suitable for electo-optical and non-linear optical applications than the congruent crystal which can be grown easier. The growth difficulties justify the development of alternative preparation methods or the improvement of the existing methods. This emphasizes the importance of the growth of stoichiometric LiNbO3 by the HTTSSG method not only from the K2O - Li2O - Nb2O5 ternary system but from the X2O - Li2O - Nb2O5 (X = Rb, Cs) systems as well. Using the presented phase diagrams the composition of the solution in equilibrium with a custom-tailored crystal composition can be determined. Most of the physical parameters of LiNbO3 are composition dependent. Even a small change of the composition (± 0.01 mol%) causes large differences in the properties so the determination of the Li2O content of the crystal with the highest possible accuracy is very important. The accuracy of the composition determination – described in detail in my thesis – based on the intensity ratio of the two characteristic bands in the OH- vibrational spectra of the near stoichiometric LiNbO3 satisfies the requirements imposed by the applications

La modélisation multi-échelle des matériaux hétérogènes est un challenge en mécanique numérique. Dans le contexte non linéaire, les propriétés effectives des matériaux hétérogènes ne peuvent pas être obtenues par les techniques utilisées pour les milieux linéaires car le principe de superposition n'est plus valable. Ainsi, dans le contexte des éléments finis, une alternative au maillage de l'ensemble de la structure avec la prise en compte de toutes les hétérogénéités, est l'utilisation de la méthode d'éléments finis multi-échelles (EF2). Les techniques de ce type offrent de nombreux avantages, tels que la prise en compte : des grandes déformations au niveau micro et macro sont souvent résolus par les procédures classiques de Newton-Raphson, qui sont généralement adaptées à la résolution des problèmes non linéaires mais qui présentent des difficultés en présence d'instabilités. Dans cette thèse, la combinaison de la méthode des éléments finis multi-échelles (EF2) et la méthode asymptotique numérique (MAN), surnommée MAN multi-échelle, permet de mettre en œuvre une technique numérique efficace pour traiter les problèmes d'instabilités dans le cadre des matériaux hétérogènes. Ces instabilités peuvent survenir à la fois au niveau micro et au niveau macro. Différentes classes de comportement des matériaux ont été implantées dans notre procédure. Pour améliorer le conditionnement du problème multi-échelle à résoudre, une technique d'homogénéisation du second ordre a été également adaptée dans le cadre de la technique MAN multi-échelle. Par ailleurs, afin de réduire le temps de calcul, quelques techniques ont été proposées dans ce travail / The multiscale modelling of the heterogeneous materials is a challenge in computational mechanics. In the nonlinear case, the effective properties of heterogeneous materials cannot be obtained by the techniques used for linear media because the superposition principle is no longer valid. Hence, in the context of the finite element method, an alternative to mesh the whole structure, including all heterogeneities, is the use of the multiscale finite element method (FE2). These techniques have many advantages, such as taking into account : large deformations at the micro and macro scales, the nonlinear constitutive behaviors of the material, and microstructure evolution. The nonlinear problems in micro and macro scales are often solved by the classical Newton-Raphson procedures, which are generally suitable for solving nonlinear problems but have difficulties in the presence of instabilities. In this thesis, the combination of the multiscale finite element method (FE2) and the asymptotic numerical method (ANM), called Multiscale-ANM, allows one to obtain a numerical effective technique for dealing with the instability problems in the context of heterogeneous materials. These instabilities can occur at both micro and macro levels. Different classes of material constitutive relation have been implemented within our procedure. To improve the multiscale problem conditioning, a second order homogenization technique was also adapted in the framework of Multiscale-ANM technique. Furthermore, to reduce the computational time, some techniques been proposed in this work

Le présent travail visait à développer et à étudier une nouvelle méthode de fabrication des guides d'onde, « High Index Soft Proton Exchange » (HISoPE) , qui permet de réaliser des guides d'onde très confinés (dne = 0,1). Des caractérisations en génération d’harmonique localisée ont montré que les propriétés non linéaires de ces guides ne sont pas détruites mais les guides HISoPE réalisés sur coupe Z peuvent présenter des modes ayant des pertes élevées. Ces pertes peuvent être éliminés en utilisant des bains plus acides pour l’échange, mais cela implique des déformations plus importante dans les guides d'onde canal et la nature hybride des modes propageant. Dans le cadre du projet PhoXcry, nous avons essayé de réaliser un modulateur électro-optique très efficace en combinant des cristaux photoniques et des guides HISoPE sur coupe X. Dans les meilleurs guides d'onde fabriqués sur coupe X, les pertes à la propagation, expliquées par la nature hybride des modes sont de 1.75dB/cm, mais dans les guides d'ondes nanostructurés il n’a pas été possible d’identifier clairement une bande interdite photonique. HISoPE en combinaison avec l’échange inverse (Reverse Proton Exchange, RPE) a montré un grand potentiel pour la fabrication de guides d'onde enterrés. Dans une expérience de SHG et malgré des pertes élevées de 2dB/cm, nous avons pu estimer une efficacité de conversion de 160%/W*cm2. Un comportement du coupleur directionnel a été observé dans les guides d'onde enterrés fabriqué en raison d’une cinétique RPE différente dans les différentes parties du guide d'onde. Un développement ultérieur de la méthode HISoPE + RPE devrait permettre d'améliorer la qualité des guides enterrés. / The present work aimed to develop and to study a new method of waveguide fabrication, High Index Soft Proton Exchange (HISoPE), which allows realizing highly confining waveguides in LiNbO3 (dne=0.1). Characterizations by localized SHG experiments, showed that the nonlinear properties of the HISoPE waveguides are not destroyed, but modes with high propagation loss were observed for planar HISoPE waveguides on Z-cut wafers. These losses can be eliminated by performing the exchange in more acidic bath, but this results in more important deformations in channel waveguides and in the hybrid nature of the propagating modes. In the frame of the project PhoXcry, we tried to realize a highly efficient electro-optical modulator by combining photonic crystals and HISoPE waveguides on X-cut wafers of LiNbO3. The losses of the waveguides fabricated on X-cut, attributed to the hybrid nature of the propagating modes, were estimated to be around 1.75dB/cm. The nanostructured waveguides exhibited high losses and it was not possible to identify clear optical band gap. HISoPE in combination with reverse proton exchange (RPE) showed a great potential for buried waveguide fabrication. We used them in a SHG experiment and despite elevated losses of 2dB/cm, the conversion efficiency was estimated as high as 160%/W*cm2. A directional coupler behavior was observed in the buried waveguides due to different RPE kinetics in different parts of the waveguide. A further development of the HISoPE+RPE process will improve the quality of the buried waveguides.

Les réseaux réflecteurs imprimés connaissent un fort développement depuis la fin des années 80. Ce type d’antenne offre la possibilité de former des diagrammes de rayonnement complexes avec une relative simplicité, un faible coût de réalisation, de faibles pertes et un volume réduit. Cependant, il souffre encore de quelques défauts : - La non régularité de la géométrie de la cellule sur la surface du réseau, dans le cas d’une antenne passive, peut engendrer des dégradations sur le diagramme de rayonnement, surtout à la transition entre deux géométries extrêmes, lorsqu’un nouveau cycle de phase commence. - Le nombre relativement élevé de composants utilisés pour contrôler la phase de l’onde réfléchie, dans le cas d’une antenne reconfigurable, augmente le coût de fabrication de l’antenne et complexifie le circuit de commande des éléments reconfigurables. - La limitation en bande passante, qui a longtemps cantonné ce type d’antenne à des applications bande étroite, est principalement liée au comportement de la cellule unitaire constitutive du réseau. Dans cette thèse, nous nous intéressons donc à la conception de nouvelles topologies de cellules déphaseuses, passives et surtout reconfigurables, qui permettent, tout en conservant une relative simplicité de réalisation, d’offrir une large bande passante. De plus, le contrôle de la phase, dans le cas des cellules reconfigurables, doit être réalisé avec un nombre réduit de composants afin de respecter la contrainte de faible coût de fabrication. / Microstrip printed reflectarrays experience strong development since the late 80s. This type of antenna has the potential to form complex radiation patterns with relative simplicity, low cost, low losses and low profile. However, it still has some shortcomings: - The non-regularity of the geometry of the cell on the surface of the array, in the case of passive antenna, may cause degradation on the radiation pattern, especially at the transition between two extreme geometries, when a new phase cycle begins. - The high number of components used to control the phase of the reflected wave, in the case of reconfigurable antenna, increases the manufacturing cost of the antenna and complicates the control circuit of the reconfigurable elements. - The limited bandwidth has long confined this type of antenna to narrowband applications and is mainly due to the intrinsic bandwidth of the unit cell. In this thesis, we therefore focus on the design of new, passive and reconfigurable, phase-shifting cells that can provide a wide bandwidth while maintaining a simple implementation. In addition, the control of the phase, in the case of reconfigurable cells, must be made with a reduced number of components in order to comply the constraint of low manufacturing cost.

Le développement des nouvelles technologies, de l'informatique et d'internet ces dernières décennies s'est accompagné d'une demande croissante de supports de stockage de l'information. En particulier, le stockage optique de données. Les supports conventionnels (Cd-Rom, Blue-ray…), basés sur un stockage en surface du disque, ont atteint aujourd'hui leurs limites. Une nouvelle technologie en cours de développement, basée sur un stockage de données en trois dimensions, est une alternative prometteuse pour supplanter les supports conventionnels. Les matériaux doivent intégrer des entités aux propriétés photochromes (commutateur moléculaire) et d'absorption multi-photonique démontrées. Au cours de cette thèse, nous avons envisagé la synthèse de molécules fonctionnelles qui présenteraient ces deux caractéristiques. Dans un premier temps, nous nous sommes intéressés à la synthèse et à l'étude comparative de complexes homo et hétéroleptiques de ruthénium présentant une certaine potentialité pour l'absorption à deux photons. Nous avons montré que ces systèmes étaient très actifs et qu'ils permettaient d'accueillir une entité photochrome sans perte des propriétés bi-photoniques. Nous avons par la suite étudié les propriétés en photochromisme de complexes de ruthénium tris-bipyridine originaux comprenant un motif azobenzène. La complexation au métal bouleverse profondément le photochromisme de l'azobenzène avec des caractéristiques cinétiques notablement différentes de celles des ligands. Enfin, nous avons étudié les propriétés de complexes hybrides composés de ligands pour l'absorption à deux photons et d'autres pour le photochromisme, ces composés sont actifs dans ces deux domaines. Nous discutons finalement du potentiel du comportement optique de films composés de ces complexes. Nous avons procédé à des essais préliminaires de modulation du signal SHG sur ces films. Nous souhaitons optimiser toutes les composantes du processus afin de déterminer le potentiel de ces composés en stockage optique de données. / The development of new technologies , computer and internet in recent decades has been accompanied by an increasing demand for media storage information. In particular , the optical data storage .Conventional media ( CD-ROM , Blu-ray ... ) based on a storage disk surface, have now reached their limits. A new technology being developed , based on a data storage in three dimensions, is a promising alternative to replace conventional materials. Materials must include entities for photochromic properties (molecular switch) and multi- photon absorption demonstrated. In this thesis , we considered the synthesis of functional molecules which present these two characteristics. As a first step , we are interested in the synthesis and comparative study of homo and heteroleptic ruthenium complex having a certain potential for two-photon absorption . We have shown that these systems were very active and they allowed to host a photochromic entity lossless bi- photonic properties. We then studied the properties of originals photochromic ruthenium tris- bipyridine containing an azobenzene motif. The metal complexing profoundly changes the photochromism of azobenzene with significantly different from those ligands kinetics . Finally, we studied the properties of hybrid complexes of ligands for the two-photon absorption and others to photochromism , these compounds are active in both areas. We finally discuss the potential of the optical behavior of compounds of these complex films. We conducted preliminary tests of the SHG signal modulation on these films. We want to maximize all components of the process to determine the potential of these compounds in optical data storage .

La détermination des caractéristiques vibratoires est une étape cruciale lors du processus de certification d'une structure ou lors de l'évaluation de sa détérioration. Quelle que soit la technique utilisée, l'identification des modes d'une structure nécessite d'effectuer des essais de vibration. Durant ces essais, on effectue des mesures qui permettent d'observer la réponse du système à une excitation connue, partiellement connue ou inconnue. Les méthodes d'identification modale sont des méthodes globales. Elles permettent de déterminer les modes de la structure à partir d'un post-traitement des données mesurées. Certaines techniques d'identification modale, telles que les méthodes basées sur la minimisation d'erreur de modélisation ARX ou les méthodes sous-espaces, se sont fortement répandues dans le monde industriel ces quinze dernières années. Cependant, malgré ce succès, ces méthodes ne fourniront jamais les paramètres modaux exacts de la structure étudiée. En effet, à cause de nombreuses sources d'erreurs (e.g. bruits sur les capteurs, bruits environnementaux), l'identification modale à partir de données réelles ne peut mener qu'à des estimations des paramètres modaux. Une estimation est par définition entachée d'une erreur statistique, la variance. Le but de cette thèse est l'étude approfondie de l'estimation de cette erreur statistique. Le manuscrit est découpé en trois partie. La première concerne la description détaillée de deux familles d'identification modale (minimisation d'erreur de modélisation ARX et méthodes sous-espace). La seconde partie décrit le principe de l'estimation des incertitudes et fournit les chaînes de propagation des variances. La troisième partie présente une validation poussée de ces méthodes dans le but d'une utilisation opérationnelle. / Determinination of vibration characteristics is a crucial step for structure certification or health monitoring. Whatever the method, modes identification of a structure requires vibration tests. During these tests, measurements are made in order to observe responses of the system to known, partially known or unknown excitations. Modal identification methods are global methods. They allow the determination of modes from data post-treatment. Some modal identification methods, such as those based on ARX modelisation error minimisation or subspace methods, have been largely spread into the industry for about fifteen years. Despite this succes, these methods will never provide exact results. Indeed, due to numerous error sources (e.g. sensor noises, environmental noises), modal identifications based on measured data can only lead to modal parameters estimations. By definition, an estimation is subject to statistical errors, also called variances. The purpose of this thesis is to give an in-depth study of these statistical errors estimation. This document is divided in three parts. Firstly, a description of two kind of modal identification methods (ARX modelisation error minimisation and subspace methods) is given. Secondly, the uncertainty estimation principle is explained and propagation schemes of variances are given. Finally, the third part provides an in-depth validation of these methods in order to proceed to industrial utilisations.

La précision des calculs théoriques d'électrodynamique quantique dans l'atome d'hydrogène est actuellement limitée par la constante de Rydberg et la distribution de charge du proton. La comparaison entre ces calculs et les mesures expérimentales de deux fréquences de transition dans l'hydrogène permet d'extraire ces deux constantes. La mesure de la transition 1S-2S est la plus précise à ce jour avec une incertitude relative de 10-15. L'objectif de mon travail de recherche est d'améliorer la précision de mesure de la fréquence de la transition 1S-3S, pouvant être utilisée comme la deuxième mesure nécessaire.La transition 1S-3S est sondée par une excitation à deux photons à 205 nm, permettant de s'affranchir de l'effet Doppler du 1er ordre. Ce faisceau UV est produit par somme de fréquence dans un cristal non linéaire. L'onde lumineuse délivrée par un laser Titane-saphir à 894 nm est sommée avec un faisceau à 266 nm produit par doublage d'un laser Nd-YO4. Cette somme de fréquence délivre un faisceau continu à 205 nm d'une puissance de 15 mWLa distribution de vitesse du jet atomique, dont la connaissance est indispensable pour évaluer l'effet Doppler du 2ème ordre, est déterminée grâce à l'effet Stark motionnel où l'action d'un champ magnétique produit un décalage en fréquence quadratique en vitesse.Les fréquences des deux lasers sources sont mesurées à l'aide un peigne de fréquence optique.La fréquence de la transition 1S-3S est finalement déterminée avec une incertitude relative de 10-12. Sa valeur conduit à une valeur préliminaire du rayon du proton qui serait en contradiction avec celle préconisée par le CODATA. / The uncertainty of the Quantum Electrodynamics calculations for hydrogen atom is currently limited by the knowledge of the Rydberg constant and the proton charge radius. Those two quantities can be extracted from the comparison between the theoretical predictions and two different frequency measurements on hydrogen.The 1S-2S transition frequency is one measured with the highest resolution with a relative uncertainty of 10-15. The aim of this thesis is to improve the determination of the 1S-3S transition, which can be used as the second precise measurement. The 1S-3S two-photon transition is excited at 205 nm. This UV light beam is generated by frequency mixing in a non-linear crystal. An 894 nm light delivered by a Ti:Sa laser is mixed with a 266 nm light beam generated by a quadrupled Nd:YVO4 laser. A reliable 15 mW continuous radiation at 205 nm is then produced. The frequencies of both lasers are measured simultaneously using an optical frequency comb referenced to a cesium clock. To evaluate the second-order Doppler effect, the velocity distribution of the atomic beam is determined thanks to a motional Stark effect. This effect is realized with a static magnetic field which induces a velocity-dependent quadratic frequency shift. Finally, the frequency of the 1S-3S transition is determined with a relative uncertainty of 10-12 which is accurate enough to contribute to the “proton size puzzle”. However, depending on the velocity distribution used in the analysis, the obtained value agrees or not with the present recommended CODATA value.

Turbulence d’ondes dans les plaques minces en vibration : étude expérimentale et numérique de l’effet de l’amortissement La théorie de turbulence d'ondes a pour but de décrire le comportement à long terme de systèmes faiblement non linéaires hors équilibre. Pour les plaques minces en vibration, ce formalisme permet de prédire un spectre de Kolmogorov-Zakharov (KZ) avec un flux d'énergie transféré des échelles d'injection à celles de dissipation le long d’une fenêtre de transparence. Les études expérimentales antérieures à cette thèse ont mis en lumière des différences notables entre les spectres mesurés et les specres théoriques. La présence, dans un solide, de l’amortissement à toutes les échelles, est ici étudié afin d’expliquer ce désaccord. En contrôlant expérimentalement l’amortissement, il a été montré que la dissipation déterminait la forme des spectres. Par la caractérisation de l’amortissement, il a été trouvé que ce dernier peut être décrit en fonction de la fréquence par une loi de puissance. La dissipation expérimentale a alors pu être introduite directement dans une simulation numérique des équations de Föppl-von Kàrmàn. Cette démarche conduit à passer de la solution théorique KZ obtenue en l’absence de dissipation à des spectres très proches des spectres expérimentaux. Ces observations ne remettent pas en cause l’emploi de la turbulence d’ondes pour décrire la vibration d’une plaque mince sollicitée par un forçage d’une grande amplitude mais encouragent à étendre les outils théoriques à notre disposition en l’absence d’une fenêtre de transparence. En réalisant cette généralisation de façon phénoménologique, une nouvelle solution stationnaire, différente de KZ et valable pour toute loi d’amortissement, a pu être calculée numériquement. / Wave turbulence theory aims at describing the long time behavior of weakly non-linear, out-of-equilibrium systems. For thin vibrating plates, this framework allows predicting a Kolmogorov-Zakharov Spectrum (KZ) with an energy flux transfered from the injection to the dissipative scales along a transparency window. Previous experimental studies have pointed out some discrepancies between mesured and theoretical spectra. The fact that, in solid, damping acts at all scales, is here studied in order to explain this disagreement. By an experimental control of the dissipation, it is observed that dissipation determines the shape of spectra. Experimental measurement of the dissipation shows that damping can here be described, as a function of the frequency, by a power law. This behavior allows us to introduce directly damping in a numerical simulation of the Föppl-von Kàrmàn equations. It leads to pass from the theoretical solution KZ obtained without dissipation to spectra which are very closed to the experimental ones. These observations do not mean that wave turbulence theory should not be applied to thin plates excited by a strong forcing but encourage to extend our theoretical tools when there is no transparency window. By doing this in a phenomenological way, a new stationary solution, different from KZ and valid for any dissipation law, has been derived.

L'objectif est d'étudier la transformation du tonnerre (amplitude, spectre) pendant sa propagation depuis le canal d'éclair jusqu'à un détecteur en se basant sur deux approches complémentaires. Dans un premier temps, l'analyse des enregistrements acoustiques (audibles et infrasons) obtenus par un réseau acoustique de petite taille (50 mètres), durant une campagne d'observations qui s'est déroulée dans le Sud de la France (Automne 2012), nous fournit des cartes 3D des sources infrasonores présentes dans un canal de foudre. Nous montrons qu'elles sont parfaitement corrélées aux résultats fournis par des outils électromagnétiques très précis. Des spectres de tonnerre ont pu être associés aux différentes phases de la décharge. Des caractéristiques acoustiques (fréquence, amplitude) ont pu être précisément reliées aux différentes sources acoustiques identifiées à l'intérieur des éclairs. Il apparaît que la propagation de l'onde acoustique perturbe la localisation pour les distances supérieures à 20km. Dans un second temps le code de calcul Flhoward a été développé pour simuler la propagation non linéaire des ondes de choc acoustiques à travers un modèle de météo réaliste. Il permet d'étudier les formes d'ondes au cours de leur propagation. La simulation numérique avec Flhoward aide à comprendre l'impact des profils météorologiques sur la propagation des signatures de tonnerre. / The goal is to study the transformation of thunder (amplitude, spectrum) during its propagation from to lightning channel towards a detector, based on two complementary approaches.In a first time, the analysis of the acoustic records (audible and infrasounds) obtained with an acoustic network of small size (50 meters), during an observation campaign which took place in South of France (Autumn 2012), produces 3D maps of infrasounds sources present in a lightning channel.We show they are perfectly correlated to results obtained from very precise electromagnetic tools.Frequency contents were successfully associated to the different parts of the discharge.Acoustic parameters (frequency, amplitude) were clearly linked to the different acoustic sources identified inside the lightning channels.It appears that the propagation of the acoustic wave impacts the localization for distances greater than 20km.In a second step, the computing code Flhoward was developed to simulate the nonlinear propagation of acoustic shock waves through a realistic meteorological profile.It allows to study the waveforms during their propagation.The numerical simulation with Flhoward helps to understand the impact of the meteorological profiles on the propagation of the thunder signatures.

De par leur bonne tenue à l’impact, les matériaux composites tissés sont de bons candidats pour la conception de pièces aéronautiques. Toutefois, le manque de confiance dans les modèles se traduit par de lourdes campagnes expérimentales. L’augmentation de la part de simulation numérique et donc la réduction des coûts de certification passent par le développement d’outils permettant de dimensionner au plus juste les structures composites tissées à matrice organique (CMO) ou céramique (CMC). C’est dans ce cadre que s’inscrit cette thèse. Les formulations proposées par l’Onera a?n de décrire l’e?et de l’endommagement matriciel ont été adaptées et étendues pour prendre en compte les spéci?cités des matériaux de notre étude, en particulier la viscosité pour les CMO. Des essais sur pièces génériques ont été simulés et confrontés à l’expérience dans le but d’évaluer la pertinence des modèles développés. Ces comparaisons en partie validantes ont mis en évidence les limites des modèles dans le cadre de l’analyse de la tenue d’une structure présentant un gradient de contrainte. A?n d’améliorer les prévisions des simulations, les effets de la rupture progressive ont été pris en compte dans les formulations. Pour s’affranchir de la localisation numérique de l’endommagement, une approche originale de régularisation couplant modélisation non locale de l’endommagement et taux d’endommagement limité a été développée. Les confrontations essai/calcul ont mis en évidence l’apport d’une modélisation plus ?ne des mécanismes d’endommagement et de rupture sur l’étude de la tenue de la structure. Compte tenu des nombreuses sources d’incertitudes, quelle confiance accorder à la simulation ? Pour répondre à cette question, l’effet des incertitudes sur des quantités d’intérêts (contrainte à rupture) a été évalué. Par ailleurs, une analyse de sensibilité (décomposition de variance) a été entreprise pour l’étude de faisabilité d’une démarche de capitalisation. Dans l’optique d’une démarche d’analyse de la tolérance aux défauts, l’efficacité de notre approche à traiter diverses singularités (taille, forme) est démontrée. En?n, les limites de la modélisation macroscopique sont discutées. / Due to their high specific properties, the use of woven composite materials has become an important commercial issue in the aeronautical applications. However, the lack of confidence in classic simulation tools used for design and conception of new structures leads to huge experimental campaigns. The increase of simulations and so the decrease of certification costs requires tools development allowing a more accurate design of woven composite structures with polymer (PMC) or ceramic matrix (CMC). This present thesis is focused on that point. The current formulations are extensions of the damage approach developed at Onera to describe matrix damage. Particularly, the nonlinearity due to viscosity for PMC is taken into account. Various structural tests on generic parts has been simulated and compared to experimental results in order to evaluate the efficiency of the models. The simulations are in good agreement with the experiments except for structure with high stress gradient. In order to improve the simulations accuracy, the effect of progressive failures are introduced. To avoid the numerical problems such as damage localization or mesh sensitivity associated with the softening behavior, an original method coupling delay effect and non local approach is proposed. A good agreement between experimental and finite element calcultion results was shown. Given many uncertainties, how much can the simulation be trusted ? In order to answer this question, the influence of uncertainties on given quantities of interest (stress failure for example) has been evaluated. Moreover, a sensivity analysis has been performed for the study of an experiment capitalisation approach. In the prospect of a reliable damage tolerance analysis, the efficiency of our approach to deal with various singularities (size, form) is shown. The limits of a macroscopic modelling are ?nally discussed.