• Refine Query
  • Source
  • Publication year
  • to
  • Language
  • 48
  • 15
  • 4
  • Tagged with
  • 68
  • 68
  • 32
  • 31
  • 20
  • 18
  • 12
  • 12
  • 10
  • 8
  • 8
  • 8
  • 8
  • 8
  • 8
  • About
  • The Global ETD Search service is a free service for researchers to find electronic theses and dissertations. This service is provided by the Networked Digital Library of Theses and Dissertations.
    Our metadata is collected from universities around the world. If you manage a university/consortium/country archive and want to be added, details can be found on the NDLTD website.
31

Contribution aux assemblages hybrides (boulonnés/collés) - Application aux jonctions aéronautiques

Paroissien, Eric 21 November 2006 (has links) (PDF)
Un assemblage est dit hybride quand il associe deux techniques d'assemblage différentes. La présente étude s'intéresse à la combinaison du collage et du boulonnage. Les jonctions longitudinales de fuselage d'avions civils, montées en simple cisaillement, constituent un cas d'application des travaux. Des outils paramétriques sous la forme de programmes informatiques sont mis au point et donnent la répartition des transferts de charge dans les substrats et dans les fixations ainsi que la répartition des contraintes adhésives, en fonction des paramètres géométriques et mécaniques de l'assemblage. Ces outils reposent sur le développement de modélisations analytiques 1D et 2D, qui supposent un comportement élastique et isotrope des matériaux, et qui sont basées sur la création de nouveaux éléments finis formulés pour simuler deux substrats collés. Un modèle 3D EF est développé afin d'accéder aux comportements mécaniques locaux, et dont l'utilisation est intéressante dans la détermination des raideurs des fixations utilisées dans les modélisations analytiques. La mesure expérimentale du taux de transfert aux fixations sous faibles charges à l'aide de boulons instrumentés valide et calibre les approches analytiques et numériques, tandis que les essais en fatigue estiment les gains en durée de vie.
32

Intégration et optimisation de l'assistance a la mise en forme par vibrations

KHAN, Armaghan 09 July 2013 (has links) (PDF)
Le travail présenté intègre l'utilisation des vibrations mécaniques dans les procédés de mise en forme à froid des matériaux de comportement viscoplastique. Les avantages de procédé de mise en forme assisté par vibrations se mesurent par une réduction de l'intensité de la force de mise en forme, une diminution de la contrainte d'écoulement, la réduction du frottement entre l'outil et lopin, une distribution uniforme de la déformation viscoplastique et une amélioration de l'état de surface de la pièce finale. A cet effet, un modèle analytique avec une seule source de vibration a été développé pour le procédé de forgeage. Les vibrations sinusoïdale et triangulaire ont été intégrées dans le modèle et il se trouve que les vibrations triangulaire est plus avantageux comparer à les vibrations sinusoïdal a cause du contrôle sur le rapport cyclique. On observe également que les paramètres principaux qui contrôlent la réduction de la force de forgeage sont : le rapport de vitesse définie par l'amplitude de la vitesse de vibration et la vitesse moyenne de l'outillage et la sensibilité à la vitesse de déformation. Le modèle analytique a été validé au moyen des simulations éléments finis réalisés avec le code de calcul FORGE2011 ® et les essais expérimentaux en dispositif expérimentale de pot vibrant a été conçu et les essais ont été réalisé sur les machines Lloyd LR30K ou ZWICK / Roell 1200. L'étude a ensuite été enrichi avec les plusieurs sources de vibrations. L'étude cinématique montre que les vibrations multiples générés une onde progressive. Cette nouvelle cinématique a été intégrée dans un modèle numérique sur le code de simulation par éléments finis FORGE2011 ® pour étudier le forgeage assisté par plusieurs sources des vibrations. L'analyse de résultats montre une réduction considérable de la force de forgeage et de la contrainte de l'écoulement ainsi qu'une modification de la loi de frottement générer par le changement de la direction du vecteur de vitesse.
33

Comportement mécanique des poutres cellulaires à ouvertures sinusoïdales : développement d'un modèle anlytique adapté

Durif, Sébastien 08 November 2012 (has links) (PDF)
L'objectif de ce travail de thèse est de développer une approche analytique permettant de définir la charge ultime d'une poutre cellulaire à ouvertures sinusoïdales. En effet, l'évolution des techniques de production a permis le développement d'une nouvelle forme de poutre cellulaire munie d'ouvertures sinusoïdales : la poutre AngelinaTM. Cette nouvelle forme d'ouverture implique de nouveaux modes de ruine. De ce fait, en vue de développer un modèle de calcul analytique adapté à cette nouvelle forme d'ouverture, une campagne d'essais expérimentaux a été menée sur des poutres cellulaires à ouvertures sinusoïdales à grande échelle (≈10m). Au travers de ces études expérimentales sur trois configurations de poutre, nous avons montré que le principal mode de ruine est lié à la flexion Vierendeel. En effet, la flexion locale des membrures de l'ouverture la plus sollicitée engendre soit la formation de 4 rotules plastiques aux 4 coins de l'ouverture, soit l'instabilité locale des parois d'âme comprimées. Un modèle aux éléments finis a été développé sur le logiciel SAFIR afin d'analyser les différents modes de ruines observés. Ce modèle a été validé sur la base des résultats expérimentaux et nous a permis d'identifier deux points particuliers : d'une part l'existence d'un maintien rotationnel entre le montant intermédiaire et la paroi d'ouverture et d'autre part, la ruine de l'ouverture ne se produit qu'au travers d'un mécanisme combinant les ruines des différents quarts d'ouverture. Une seconde campagne d'études expérimentales et numériques a ensuite été menée sur des parties isolées, extraites des poutres préalablement testées, afin d'étudier de manière locale le comportement à la flexion des quarts d'ouverture. Ces études ont servi à valider un second modèle aux éléments finis, développé sur le logiciel Cast3m. Celui-ci nous a permis, au travers d'une étude paramétrique, de quantifier le maintien rotationnel apporté par le montant intermédiaire sur la paroi d'âme d'ouverture adjacente. Cette étude a confirmé l'importance de la rigidité apportée par le montant intermédiaire aux parois d'âme adjacente. Ainsi, cet apport de rigidité doit être pris en compte dans l'approche analytique pour définir de manière réaliste la résistance au voilement local des différentes parties d'une ouverture sinusoïdale. Finalement, cette thèse a abouti au développement d'un nouveau modèle analytique de calcul de la résistance ultime des parois d'une ouverture sinusoïdale. Du fait des éventuelles instabilités locales, le modèle analytique s'est appuyé sur des éléments théoriques de stabilité des plaques. De plus, une étude numérique détaillée du mécanisme de ruine d'une ouverture isolée nous a permis de justifier une approche cinématique de ruine de l'ouverture sinusoïdale. Cette approche combine les résistances ultimes des différents quarts d'ouverture. Le modèle analytique proposé permet de considérer à la fois la résistance ultime de chaque partie de l'ouverture et leurs modes de ruine. Une étude comparative avec des résultats numériques a montré que ce modèle est fiable et représentatif de la réalité pour caractériser l'état limite ultime des poutres cellulaires à ouvertures sinusoïdales.
34

Conception, modélisation et réalisation de composants inductifs intégrés pour alimentations de faible puissance et microsystèmes

Artillan, Philippe 27 November 2008 (has links) (PDF)
L'intégration complète des dispositifs de conversion d'énergie destinés à créer des circuits d'alimentation compacts se heurte encore aujourd'hui à des contraintes technologiques fortes sur l'intégration des composants passifs inductifs et capacitifs. La fréquence de fonctionnement de ces dispositifs, d'ores et déjà de l'ordre du MHz, en augmentant pourrait réduire la taille de ces composants passifs et donc permettre leur intégration. Cependant, les pertes fortement liées à la fréquence et à la technologie, freinent encore cette augmentation.<br />Cette thèse a pour objectif l'étude systématique de structures d'inductances intégrées à travers le développement d'une modélisation accompagnée de simulations précises pour le développement d'une méthodologie de simulation valable sur la plage de fréquence de 0Hz à 1GHz. Pour cela, un modèle analytique basé sur la méthode PEEC (Partial Element Equivalent Circuit) a d'abord été développé. Ce type de modèle a du être adapté aux topologies des composants étudiés : symétrie cylindrique, conducteurs de section non négligeable, milieux non homogènes... Ce modèle décompose l'étude du composant en deux étapes : calculs électromagnétiques menés sur des sous-parties du composant (éléments partiels) et calcul de la réponse fréquentielle globale à partir de modèles électriques des éléments partiels. Un procédé complet de réalisation technologique des composants inductifs a également été mis au point en parallèle afin de valider expérimentalement la modélisation. L'amélioration de la résolution des motifs et l'augmentation de l'épaisseur des dépôts lors des étapes de croissance électrolytique du cuivre, ainsi que la planarisation des niveaux de conducteurs avec de la résine SU8 ont fait l'objet d'études spécifiques. Ces avancées ont permis la réalisation de selfs multibrin et multiniveau qui présentent des caractéristiques électriques intéressantes en haute fréquence. Enfin, la caractérisation des prototypes en basse fréquence a été rendue possible par la mise en oeuvre d'un banc de mesure d'impédance pour la gamme de fréquence de 40 Hz à 110MHz et des caractérisations en haute fréquence effectuées à l'aide d'un analyseur de réseau. Ces travaux ont été complétés par la mise au point d'une méthode analytique de calcul du champ magnétique. Les calculs permettent de prédire la répartition du champ magnétique généré par des actionneurs magnétiques planaires intégrés excités en courant continu.
35

Déploiement efficace de services complexes dans l'infrastructure de cloud

Tran, Khanh Toan 10 January 2013 (has links) (PDF)
Le but du travail réalisé dans cette thèse est de fournir aux fournisseurs de services une solution qui est capable de non seulement déployer les services dans le cloud de façon économique, automatique, mais aussi à grande échelle. La première contribution traite le problème de la construction d'un nouveau service demandé par le client à partir de services disponibles dans le cloud de manière à satisfaire les exigences en termes de qualité de service ainsi qu'en termes de coût. Nous présentons ce problème par un modèle analytique et proposons un algorithme heuristique dont la performance est améliorée de 20-30% par rapport aux autres approches. La seconde contribution est une solution pour déployer les services dans le cloud en considérant les demandes des utilisateurs finaux. Pour assurer qualité de services dans une grande échelle, le service demandé est dupliqué et distribué dans le réseau; chacun de ses réplicas servira les utilisateurs à proximité. Le plan d'approvisionnement selon lequel le service est dupliqué dépend de sa demande, ce qui ne cesse pas de changer en quantité ainsi qu'en distribution, ce qui rend le problème plus compliqué. Nous proposons une solution qui est capable de s'adapter aux changements dans le réseau, y compris ceux des demandes de ses utilisateurs. Enfin, nous proposons un système basé sur OpenStack qui permet de déployer les services complexes dans un cloud qui couvre différente locations (cloud multi-site). A partir d'une demande du client, le système automatiquement calcule le plan d'approvisionnement optimal et le déploie en respectant les contraintes du client.
36

Influence de la forme de la section transversale sur l'écoulement dans un canal: modélisation, simulation et expérience

Wu, Bo 23 January 2014 (has links) (PDF)
La modélisation des phénomènes physiologiques induits par un écoulement, tels que l'écoulement sanguin au travers d'une sténose ou l'écoulement d'air lors de la production de parole, repose souvent sur des théories quasi-unidimensionnelles ou bi-dimensionnelles. Cependant, il est établi que le développement des couches limites dépend de la section transversale. Le but de cette thèse est de modéliser, simuler et caractériser l'importance potentielle de la section transversale sur les écoulements laminaires, contrôlés en pression, en l'absence ou en présence d'une constriction. Des coordonnées de translation sont utilisées pour obtenir des solutions pour des écoulement visqueux au travers d'une section de forme arbitraire. Cette paramétrisation est appliquée à la résolution des équations physiques pour des formes à deux et à trois dimensions. Un modèle d'écoulement simplifié quasi-tridimensionnel, qui prend en compte les pertes dissipatives par convection, la viscosité et la forme de la section est présenté et appliqué à la description de l'écoulement le long d'une sténose. Des données expérimentales et issues de simulations numériques sont collectées afin de caractériser l'influence de la forme de la section transversale dans le cas d'une constriction. simulation numérique sont comparées.
37

Modélisation analytique et outils pour l'optimisation des transformateurs de puissance haute fréquence planars

Besri, Abdelhadi 26 May 2011 (has links) (PDF)
Avec un encombrement minimal et rendement de 99% les transformateurs planars haute fréquence sont des solutions technologiques importantes pour les systèmes de conversion et de conditionnement de l'énergie électrique en électronique de puissance. L'objectif de l'étude consiste à rechercher un modèle du composant permettant à la fois de simuler son comportement dans un circuit d'électronique de puissance et d'optimiser ses performances lors de sa conception. Les modèles proposés prennent en compte les propriétés physiques et géométriques du composant. Les validations sont basées sur des confrontations avec les mesures. Une autre motivation forte de ce travail est que les modélisations numériques doivent prendre en compte simultanément les effets des courants induits et les capacités parasites, ce qui conduit, en haute fréquence, à des besoins prohibitifs en termes de capacité mémoire ou de temps de calcul. En se limitant aux transformateurs planar. La modélisation (dite LEEC) présentée ici s'appuie sur une discrétisation à échelle intermédiaire : spire par spire, c'est-à-dire couche par couche. Elle assemble deux approches analytiques déjà introduites par l'équipe : la première traite les aspects électrocinétiques (incluant les courants induits) et magnétiques et l'autre les aspects électrostatiques. Les circuits à constantes localisées basés sur la méthode LEEC montrent un très bon accord avec toutes les mesures jusqu'à 40 MHz. Des outils numériques sont aussi développés pour faciliter l'obtention de ces différents circuits en partant : soit de la description du composant, soit de mesures d'impédances. Pour compléter, la représentation des capacités d'un transformateur quelconque, développée en plusieurs étapes par l'équipe, est résumée et des précautions expérimentales originales appliquées à l'impédancemétrie sont exposées dans le détail. D'autres travaux visant à élargir le champ d'application de la méthode LEEC sont également présentés.
38

Etude de qualité et de stabilité des réseaux aéronautiques embarqués HVDC / Quality and stability studies of aircraft embedded networks HVDC

Charrada, Madiha 21 February 2013 (has links)
Dans le contexte du concept de l’avion «plus électrique », l’électricité prend une part croissante dans les systèmes embarqués aéronautiques. Ainsi, de nombreuses évolutions ont eu lieu, concernant les réseaux de bord, notamment l’architecture du réseau qui est devenue de plus en plus complexe, intégrant de nombreux équipements de différentes natures ; ceci entraîne de multiples questions relatives à la qualité et la stabilité du réseau. C’est dans ce contexte que s’inscrivent ces travaux de thèse, qui aboutissent d’une part à la mise en œuvre de méthodes d’analyse de la stabilité et de la qualité du réseau aéronautique embarqué et d’autre part à l’analyse de couplage entre plusieurs équipements connectés sur un seul réseau. Les études menées au cours de ces travaux reposent sur l’utilisation des expressions analytiques représentant le comportement fréquentiel de réseaux continus. Une méthode de calcul basée sur les lois de la physique a été développée et validée par des approches de simulation et d’expérimentation. Par ailleurs, en absence d’un modèle analytique, nous choisissons de déterminer l'expression littérale du modèle équivalent du réseau à partir de l'identification des données issues de l’approche expérimentale en se basant sur des algorithmes d'optimisation. Ces modèles analytiques sont associés aux critères mathématiques de stabilité (Routh-Hurwitz, MiddleBrook), et aux contraintes de qualité imposées par les standards aéronautiques, dans le but de dimensionner les éléments de filtrage à placer entre le réseau et l’équipement. Le calcul des composants assurant le filtrage repose sur la recherche d’une solution optimale : nous introduisons les algorithmes d’optimisation utilisés en première étape, et nous proposons par la suite une formulation de la procédure de calcul des paramètres du filtre pour une meilleure illustration des résultats intéressants pour un concepteur de réseau. Enfin, nous clôturons ce travail par une analyse des couplages entre plusieurs équipements connectés sur un réseau continu et nous concluons par une approche expérimentale afin d’appuyer les résultats théoriques. / In the concept of the "more electric" aircraft, power takes an increasing proportion in aeronautic embedded systems. So many changes have occurred on board networks, including the network architecture that has become more and more complex, integrating many different types of equipment that causes quality and stability issues on the network. In this context that this thesis is inscribed, firstly it leads to the implementation of methods for stability and quality analyzing of the embedded network and secondly to coupling analysis of several devices connected to a single network. Studies conducted in this work are based on analytical expressions representing the frequency behavior of networks. A compute based on fundamental physical laws has been developed and validated by simulation and experimental approaches. Moreover, in the absence of an analytical model, we choose to determine the literal expression of the equivalent model of the network by identifying data, come from the experimental approach, using optimization algorithms. These analytical models are associated with mathematical criteria: Routh-Hurwitz, Middlebrook, and quality constraints imposed by aeronautic standards, in order to size the elements of the filter cell placed between the voltage source and the equipment to ensure the stability and quality of the entire system. The sizing of cell parameters based on optimal search by means of optimization algorithms in the first step, and then we propose a formulation of the calculation procedure of the filter parameters to illustrate some interesting results for designer’s network. Finally, we close this work with a coupling analysis of several devices connected to a single DC bus and we conclude by an experimental approach to support the theoretical results.
39

Comportement mécanique des poutres cellulaires à ouvertures sinusoïdales : développement d'un modèle anlytique adapté / Mechanical behaviour of cellular beams with sinusoidal openings : development of an adapted anlytic model

Durif, Sébastien 08 November 2012 (has links)
L’objectif de ce travail de thèse est de développer une approche analytique permettant de définir la charge ultime d’une poutre cellulaire à ouvertures sinusoïdales. En effet, l’évolution des techniques de production a permis le développement d’une nouvelle forme de poutre cellulaire munie d’ouvertures sinusoïdales : la poutre AngelinaTM. Cette nouvelle forme d’ouverture implique de nouveaux modes de ruine. De ce fait, en vue de développer un modèle de calcul analytique adapté à cette nouvelle forme d’ouverture, une campagne d’essais expérimentaux a été menée sur des poutres cellulaires à ouvertures sinusoïdales à grande échelle (≈10m). Au travers de ces études expérimentales sur trois configurations de poutre, nous avons montré que le principal mode de ruine est lié à la flexion Vierendeel. En effet, la flexion locale des membrures de l’ouverture la plus sollicitée engendre soit la formation de 4 rotules plastiques aux 4 coins de l’ouverture, soit l’instabilité locale des parois d’âme comprimées. Un modèle aux éléments finis a été développé sur le logiciel SAFIR afin d’analyser les différents modes de ruines observés. Ce modèle a été validé sur la base des résultats expérimentaux et nous a permis d’identifier deux points particuliers : d’une part l’existence d’un maintien rotationnel entre le montant intermédiaire et la paroi d’ouverture et d’autre part, la ruine de l’ouverture ne se produit qu’au travers d’un mécanisme combinant les ruines des différents quarts d’ouverture. Une seconde campagne d’études expérimentales et numériques a ensuite été menée sur des parties isolées, extraites des poutres préalablement testées, afin d’étudier de manière locale le comportement à la flexion des quarts d’ouverture. Ces études ont servi à valider un second modèle aux éléments finis, développé sur le logiciel Cast3m. Celui-ci nous a permis, au travers d’une étude paramétrique, de quantifier le maintien rotationnel apporté par le montant intermédiaire sur la paroi d’âme d’ouverture adjacente. Cette étude a confirmé l’importance de la rigidité apportée par le montant intermédiaire aux parois d’âme adjacente. Ainsi, cet apport de rigidité doit être pris en compte dans l’approche analytique pour définir de manière réaliste la résistance au voilement local des différentes parties d’une ouverture sinusoïdale. Finalement, cette thèse a abouti au développement d’un nouveau modèle analytique de calcul de la résistance ultime des parois d’une ouverture sinusoïdale. Du fait des éventuelles instabilités locales, le modèle analytique s’est appuyé sur des éléments théoriques de stabilité des plaques. De plus, une étude numérique détaillée du mécanisme de ruine d’une ouverture isolée nous a permis de justifier une approche cinématique de ruine de l’ouverture sinusoïdale. Cette approche combine les résistances ultimes des différents quarts d’ouverture. Le modèle analytique proposé permet de considérer à la fois la résistance ultime de chaque partie de l’ouverture et leurs modes de ruine. Une étude comparative avec des résultats numériques a montré que ce modèle est fiable et représentatif de la réalité pour caractériser l’état limite ultime des poutres cellulaires à ouvertures sinusoïdales. / The aim of this thesis is to develop an analytical approach so as to define the ultimate failure load of cellular beams with sinusoidal openings. Indeed, the evolution of the conception techniques led to the development of a new opening shape for cellular beams: the sinusoidal openings (AngelinaTM beam). This new opening shape involves new failure modes at ultimate limit state and cannot be calculated with existing methods. Thus, experimental tests have been carried out on full scale cellular beams with sinusoidal openings (≈10m span). It has been shown, through those experimental studies led on three beam configurations, that the main failure mode is linked to the Vierendeel mechanism. At ultimate limit state, the local bending of the members of the most stressed opening leads to either the formation of four plastic hinges or the local buckling of the compressed opening web panels. A finite element model has been developed on the software SAFIR in order to better understand the different observed failure modes. This model has been validated on the basis of the experimental results and allowed identifying two main points: firstly, the existence of a rotational restraint between the intermediate web-post and the opening web panel and secondly, the failure of the sinusoidal opening corresponds to a mechanism combining the failure of each opening quarters. A second experimental and numerical test campaign has been conducted on isolated parts of the previous tested beams, in order to study locally the behavior of the opening quarters under local bending. Those studies allowed validating a finite element model developed on Cast3m software. This model has been used for a parametrical study which allowed quantifying the rotational restraint supplied by the intermediate web-post on the adjacent opening panel. This study confirmed the importance of the stiffness supplied by the intermediate web-post to the opening panels. Thus, it has to be taken into account to get a representative analytical failure load of each opening quarter. Finally, this thesis led to the development of a new analytical model for the calculation of the ultimate strength of the different opening parts. Due to the possible local buckling, the model is based on theoretical elements on the stability of plates. Furthermore, a detailed numerical study of the mechanism of failure of an isolated opening allowed justifying a cinematic approach for the failure of the opening. This approach combines the ultimate strength of each opening parts. The proposed analytical model permits to consider each opening quarter failure mode and strength. A comparative study with finite element results has showed that this model is reliable and representative of reality for defining the ultimate limit state of cellular beams with sinusoidal openings.
40

The influence of the cross section shape on channel flow : modeling, simulation and experiment / Influence de la forme de section transversale sur l'écoulement dans un canal : modélisation, simulation et expérimentation

Wu, Bo 23 January 2014 (has links)
La modélisation des phénomènes physiologiques induits par un écoulement, tels que l'écoulement sanguin au travers d'une sténose ou l'écoulement d'air lors de la production de parole, repose souvent sur des théories quasi-unidimensionnelles ou bi-dimensionnelles. Cependant, il est établi que le développement des couches limites dépend de la section transversale. Le but de cette thèse est de modéliser, simuler et caractériser l'importance potentielle de la section transversale sur les écoulements laminaires, contrôlés en pression, en l'absence ou en présence d'une constriction. Des coordonnées de translation sont utilisées pour obtenir des solutions pour des écoulement visqueux au travers d'une section de forme arbitraire. Cette paramétrisation est appliquée à la résolution des équations physiques pour des formes à deux et à trois dimensions. Un modèle d'écoulement simplifié quasi-tridimensionnel, qui prend en compte les pertes dissipatives par convection, la viscosité et la forme de la section est présenté et appliqué à la description de l'écoulement le long d'une sténose. Des données expérimentales et issues de simulations numériques sont collectées afin de caractériser l'influence de la forme de la section transversale dans le cas d'une constriction. simulation numérique sont comparées. / Physical models of physiological flow-induced phenomena, such as blood flow through a stenosis or air flow during human speech production, often rely on a quasi-one-dimensional or two-dimensional flow model, so that details of the cross section shape are neglected. Nevertheless, boundary layer development is known to depend on the cross section shape. The aim of this thesis is to model, simulate and characterize the potential impact of the cross section shape for pressure-driven laminar channel flow without and with constriction. Stretched coordinates are introduced to obtain viscous flow solutions for channels with an arbitrary cross section. The proposed cross section shape parametrization is applied to solve physical equations for two-dimensional and three-dimensional shapes. A simplified quasi-three-dimensional flow model, which accounts for kinetic losses, viscosity and the cross section shape, is presented and applied to describe the flow through a stenosis. Finally, flow data are gathered experimentally and numerically in order to characterize the influence of the cross section shape in the case of a constricted channel. Modeled, experimental and numerical data are compared.

Page generated in 0.4531 seconds