The applications of membranes are increasing rapidly in various fields of engineering and science. The geometric, material, force and contact non-linearities complicate their analysis, which increases the demand for computationally efficient methods and interpretation of counter-intuitive behaviours. The first part of the present work studies the free and constrained inflation of circular and cylindrical membranes. The membranes are assumed to be in contact with a soft substrate, modelled as a linear spring distribution.Adhesive and frictionless contact conditions are considered during inflation,while only adhesive contact conditions are considered during deflation. For a circular membrane, peeling of the membrane during deflation is studied, and a numerical formulation of the energy release rate is proposed. The second part of the thesis discusses the instabilities observed for fluid containing cylindrical membranes. Limit points and bifurcation points are observed on primary equilibrium branches. The secondary branches emerge from bifurcation points, with their directions determined by eigenvectors corresponding to zero eigenvalues at the bifurcation point. Symmetry has major implications on stability analysis of the structures, and the relationship between eigenvalue analysis and symmetry is highlighted in this part of the thesis. In the third part, wrinkling in the pressurized membranes is investigated,and robustness of the modified membrane theory and tension field theory is examined. The effect of boundary conditions, thickness variations, and inflating media on the wrinkling is investigated. It is observed that, with a relaxed strain energy formulation, the obtained equilibrium solutions are unstable due to the occurrence of pressure induced instabilities. A detailed analysis of pressure induced instabilities in the wrinkled membranes is described in the thesis. / <p>QC 20160518</p>

In recent years, a number of high profile accidental fires have occurred in several road and rail tunnels throughout the world. Many of these fires grew rapidly to catastrophic size and claimed many lives. The processes involved in the rapid growth and extremely severe of these fires are not adequately understood as yet. The introduction to this thesis reviews a number of these accidental fires and describes much of the previous experimental research which has brought about the current understanding of tunnel fire behaviour. A detailed review of the relevant parts of elementary fire dynamics is also presented. This thesis addresses two main questions: 1. What is the influence of longitudinal ventilation on fire size in tunnels? and 2. What is the influence of tunnel geometry on fire size? The answers to both these questions are determined using a probabilistic method called Bayes Theorem. This provides a method of answering the above two questions using the handful of experimental data which are available. It is found that the heat release rate (HRR) of a heavy goods vehicle (HGV) fire may be greatly increased in magnitude by longitudinal ventilation, for example by about a factor of 5 with a longitudinal ventilation velocity of 3ms-1. It is also found that longitudinal ventilation may cause a significant increase in the HRR of large pool fires, but may cause a decrease in the HRR of small pool fires and car fires. An equation is derived to predict the influence of tunnel geometry on HRR. It is found that HRR varies principally with the width of the tunnel and the width of the fire object. The HRR of a fire in a tunnel my be increased up to four times due to the geometry of the tunnel.

L'analyse des phénomènes locaux comme les effets de bord libre et le délaminage dans les structures multicouches nécessite des théories fines qui donnent une bonne description de la réponse locale. Étant donné que les approches tridimensionnelles sont, en général, très coûteuses en temps de calcul et en mémoire, des approches bidimensionnelles de type layerwise sont souvent utilisées. Dans ce travail de doctorat, un modèle layerwise en contrainte, appelé LS1, est appliqué au problème du multi-délaminage dans les plaques stratifiées invariantes dans le sens longitudinal. L'invariance dans la direction de la longueur nous permet d'aborder le problème analytiquement. Dans un premier temps, nous proposons une méthode analytique pour l'analyse des plaques multicouches multi-délaminées soumises à la traction uniaxiale. La singularité des contraintes interlaminaires aux bords libres et l'initiation du délaminage en mode III sont étudiées. Un modèle raffiné, nommé LS1 raffiné, est proposé pour améliorer les approximations dans les zones de singularités telles que les bords libres et les pointes de fissure. Les résultats du modèle raffiné sont validés en les comparant avec ceux obtenus par éléments finis tridimensionnels. Dans un deuxième temps, l'approche analytique proposée est étendue à la flexion cylindrique des plaques multicouches. La propagation du délaminage en modes I et II est étudiée et les approximations du modèle LS1 sont validées. À la fin, nous généralisons la méthode analytique proposée afin de prendre en considération tous les chargements invariants dans le sens longitudinal. L'approche finale permet d'analyser les plaques multicouches rectangulaires soumises à des charges invariantes sur les faces supérieure et inférieure, les forces ou les déplacements imposés sur les bords latéraux ainsi que quatre types de chargement sur les extrémités longitudinales: traction uniaxiale, flexion hors plan, torsion et flexion dans le plan. La solution analytique du modèle LS1 est obtenue pour une plaque stratifiée soumise à tous les chargements mentionnés ci-dessus. L'approche est validée en comparant avec la méthode des éléments finis tridimensionnels pour plusieurs types de chargement / The analysis of local phenomena such as free-edge effects and delamination in multilayered structures requires the accurate theories which can provide a good description of the local response. Since the three-dimensional approaches are generally very expensive in computational time and memory, the layerwise two-dimensional approaches are widely used. In this Ph.D. thesis, a stress layerwise model, called LS1, is applied to the multi-delamination problem in longitudinally invariant multilayered plates. The invariance in the longitudinal direction allows us to solve the problem analytically. At first, we propose an analytical method for the analysis of multi-delaminated multilayered plates subjected to the uniaxial traction. The free-edge interlaminaire stress singularities and the mode III delamination onset are investigated. A refined model, called Refined LS1, is proposed in order to improve the approximations in singularity zones such as free edges and crack tips. The results of the refined model are validated by comparing them with those obtained by a three-dimensional finite element model. Afterwards, the proposed analytical approach is extended to the cylindrical bending of the multilayered plates. The propagation of delamination in modes I and II is studied and the approximations of the LS1 model are validated. At last, we generalize the proposed analytical method to take into account all invariant loads in the longitudinal direction. The final approach allows us to analyze the rectangular multilayered plates subjected to invariant loads on the top and bottom surfaces, imposed displacements or forces at the lateral edges, and also four types of loading at the longitudinal ends: uniaxial traction, out-of-plane bending, torsion and in-plane bending. The analytical solution of the LS1 model is obtained for a laminated plate subjected to all the loads mentioned above. The approach is validated by comparison with the three-dimensional finite element method for various types of loading

Les incendies de végétation sont caractérisés par de nombreuses échelles de temps et d’espace. Une approche multi-physique et multi-échelle est donc nécessaire pour aborder la complexité de ces phénomènes. Ce travail de thèse est une contribution expérimentale à l’étude des effets d’échelle et du système sur la dynamique de combustion des végétaux et les fumées émises. Notre objectif principal a été de déterminer quels protocoles expérimentaux et plus particulièrement quelles échelles permettent de caractériser au mieux la combustion des végétaux en laboratoire. Nous avons ainsi étudié le comportement au feu des aiguilles de deux espèces de pin représentatives de l’écosystème méditerranéen, Pinus pinaster et Pinus laricio. Les litières d’aiguilles de pin représentent un enjeu important car elles entretiennent la propagation des incendies et elles participent à la transition d’un feu de surface à un feu total. Trois configurations expérimentales ont été étudiées au cours de cette thèse, l’échelle du cône calorimètre, l’échelle du grand calorimètre ou LSHR (permettant la combustion statique de masses plus importantes que le cône) et enfin, la propagation dans le LSHR permettant d’étudier l’effet du système sur la dynamique de la combustion et sur la production des fumées. Pour comparer ces trois configurations les protocoles expérimentaux ont été adaptés tout en maintenant la charge de combustible. Différents paramètres ont été étudiés pour analyser la dynamique de combustion tels que l’efficacité de la combustion, l’énergie dégagée ou encore la vitesse de perte de masse. La production des fumées a été étudiée par la mesure du coefficient d’extinction qui caractérise leur opacité et permet d’obtenir le facteur d’émission des suies. Les facteurs d’émissions des principaux composés émis lors de la combustion de ces deux types d’aiguilles de pin ont été mesurés en continu à partir d’un analyseur Infrarouge à Transformée de Fourier et d’un analyseur Infrarouge Non Dispersif. Des analyses par chromatographie en phase gazeuse couplée à un détecteur à ionisation de flamme et un spectromètre de masse nous ont permis de compléter ces mesures. Un bilan massique de carbone a également été réalisé afin de quantifier le carbone total mesuré dans nos analyses. Les principales contributions de notre travail sont les suivantes : l’étude du comportement au feu des aiguilles de P. pinaster a révélé des différences importantes pour la puissance dégagée aux échelles du cône calorimètre et du LSHR. En revanche, le système de combustion (propagation) n’influence pas cette grandeur. L’efficacité de la combustion est apparue peu dépendante de l’échelle et du système de combustion. Nous avons observé une influence de l’échelle de combustion sur la production totale des fumées (RSR) et sur le facteur d’émission des suies. Nous avons également montré que le système de combustion (la propagation) influence la dynamique et la valeur des facteurs d’émission de dioxyde et de monoxyde de carbone, composés majoritairement émis par ces combustions. Une influence de l’échelle et du système est également à noter sur les facteurs d’émissions des composés azotés et des COV émis pour les trois configurations expérimentales. Nous avons attribué les différences observées aux valeurs de températures des fumées. Enfin, une influence de la géométrie des particules a été mise en évidence par comparaison des combustions réalisées avec le cône calorimètre et le LSHR pour les deux types d’aiguilles de pin. La dynamique de combustion des aiguilles de Pinus laricio est moins affectée par le changement d’échelle que celle des aiguilles de Pinus pinaster (plus faible variation de la puissance de combustion). Nous avons également observé que les aiguilles de Pinus laricio, thermiquement plus fines que les aiguilles de Pinus pinaster présentent une valeur plus faible pour le facteur d’émission des suies. / Wildfires are characterized by a lot of scales of time and space. A multi-physics and multi-scale approach is required to consider the complexity of these phenomena. This thesis is an experimental contribution to the study of the scale effects and the effects of the system on the combustion dynamics of forest fuels and smoke emission. The aim of this work was to determine which experimental protocols and specifically which scales can be used to characterize the combustion of vegetal fuels in the laboratory. The reaction to fire of pine needles species representative of the Mediterranean ecosystem, (Pinus pinaster and Pinus laricio) has been studied. Litters of pine needles are an important issue because they sustain fire spread and are involved in the transition from a surface fire to a total fire.Three experimental configurations were studied in this thesis: the cone calorimeter scale; the large scale calorimeter or LSHR (allowing static combustion with larger masses than used with the cone); a fire spread in the LSHR, allowing to study the effect of the system on the dynamics of combustion and release of smoke. To compare these configurations, experimental protocols have been elaborated while keeping the same fuel load. Different parameters were studied to analyze the combustion dynamics such as combustion efficiency, heat released rate and mass loss rate. Smoke production was studied by measuring the coefficient of extinction to characterize their opacity and an estimation of the soot emission factor was derived. The emission factors of the main compounds emitted during the combustion of these two pines needles were measured with a Fourier Transform Infrared analyzer and a Non-dispersive infrared analyzer. Analysis by gas chromatography coupled with a flame ionization detector and a mass spectrometer allowed us to complete these measurements. A mass balance of carbon was also performed to quantify the total carbon measured through our analyzes.The main contributions of our work can be summarized as follow: the study of the burning of Pinus pinaster needles reveals significant differences for heat release rate (HRR) at both cone calorimeter and LSHR scales. However, the combustion system (fire spread) does not influence the HRR value at the LSHR scale. The combustion efficiency appeared to be independent with regard to the scale and the system. We observed a wide influence on the rate of smoke release and the emission factor of soot. We also shown that the combustion system (fire spread) influences the dynamics and value of emissions factors of carbon dioxide and carbon monoxide (major compound emitted for these combustions). An influence of the three experimental configurations on the emission factors of nitrogen compounds and VOC emissions was also noted. This difference was attributed to the level of smoke temperature. Finally, an effect of particles geometry was also pointed out by the comparison between the burnings performed with the cone calorimeter and the LSHR for both types of pine needles. The combustion dynamics of Pinus laricio needles was slighlty affected by changing scale in comparison to needles of Pinus pinaster (weak variation of HRR). We also observed that Pinus laricio’s needles, which are thermally thiner than Pinus pinaster ones have the lowest soot emission factor.

In this thesis an experimental method to calculate cohesive characteristics for an adhesive layer in a End Notched Flexure (ENF) specimen is presented and evaluated. The method is based on the path independent J-integral where the energy release rate (ERR) for the adhesive is derived as a function of the applied forces and the rotational displacements at the loading point and at the supports of the specimen. The major advantage with the method in comparison with existing theory known by the job initiator is that it is still applicable with ENF specimens that are subjected to yielding of the adherends. The structure of this thesis is disposed so that the theory behind the J-integral method is shortly described and then an evaluation of the method is performed by aid of finite element simulations using beam and cohesive elements. The finite element simulations indicates that the ERR can be determined with good accuracy for an ENF specimen where a small scale yielding of the adherends has occurred. However when a fully cross sectional yielding of the adherends is reached the ERR starts diverging from the exact value and generates a too high ERR according to input data in the simulations, i.e. the exact values. The importance in length of the adhesive process zone is also shown to be irrelevant to the ERR measured according to the J-integral method. Simulation performed with continuum elements indicates that a more reality based FE- simulation implies a higher value of the applied load in order to create crack propagation. This is an effect of that the specimen is allowed to roll on the supports which makes the effective length between the supports shorter than the initial value when the specimen is deformed. This results in a stiffer specimen and thus a higher applied force is needed to create crack propagation in the adhesive layer. An experimental set up of an ENF specimen is created and the sample data from the experiments are evaluated with the J-integral method. For measuring the rotational displacements of the specimen which are needed for the J-integral equation an image system is developed by the author and validated by use of linear elastic beam theory. The system calculates the three rotational displacements of the specimen by aid of images taken by a high resolution SLR camera and the system for measuring the rotations may also be used in other applications than for a specific ENF geometry. The validation of the image system shows that the rotations calculated by the image system diverge from beam theory with less than 2.2 % which is a quite good accuracy in comparison with the accuracies for the rest of the used surveying equipment. The results from the experiment indicates that the used, about 0.36 mm thick SikaPower 498, adhesive has an maximum shear strength of 37.3 MPa and a critical shear deformation of 482 µm. The fracture energy is for this thickness of the adhesive is determined as 12.9 kJ/m2. This report ends with a conclusion- and a suggested future work- chapter.

This paper deals with two-dimensional Finite Element Analysis of the End Notched Flexure (ENF) specimen. The specimen is known to be unstable if the crack length is shorter than some critical crack length acr. A geometric linear two-dimensional Finite Element (FE) analysis of the ENF specimen is performed to evaluate acr for isotropic and orthotropic elastic materials, respectively. Moreover, the Mode II Energy Release Rate (ERR) JII and the compliance of the specimen are calculated. The influence of anisotropy is studied. Comparisons are made with the results from beam theory. This work is an extension of previous work.

The fire development and flame spread in the railway carriages are investigated by performing a set of simulations using a widely accepted simulation software called &amp / #8216 / Fire Dynamics Simulator&amp / #8217 / . Two different rolling stock models / representing a train made up of physically separated carriages, and a 4-car train with open wide gangways / have been built to examine the effects of train geometry on fire development and smoke spread within the trains. The simulations incorporate two different ignition sources / a small size arson fire, and a severe baggage fire incident. The simulations have been performed incorporating variations of parameters including tunnel geometry, ventilation and evacuation strategies, and combustible material properties. The predictions of flame spread within the rolling stock and values of the peak heat release rates are reported for the simulated incident cases. In addition, for a set of base cases the onboard conditions are discussed and compared against the tenability criteria given by the international standards. The predictions of heat release rate and the onboard conditions from the Fire Dynamics Simulator case studies have been checked against the empirical methods such as Duggan&amp / #8217 / s method and other simulation softwares such as CFAST program.

Fire is a complex phenomenon including many parameters. The nature of fire makes it a very dangerous and hazardous. For many reasons the number of tunnels are increasing on earth and fire safety is one of the major problem related to tunnels. This makes important to predict and understand the behavior of fire, i.e., heat release rate, smoke movement, ventilation effect etc. The literature includes many experimental and numerical analyses for different conditions for tunnel fires. This study investigates pool fire of three different fuel sources: ethanol, gasoline and their mixture for different ventilation conditions, different geometries and different amounts. Combustion gases and the burning rates of the fuel sources are measured and analyzed. The numerical simulation of the cases is done with Fire Dynamics Simulator (FDS), a CFD code developed by NIST.

<p>In this thesis an experimental method to calculate cohesive characteristics for an adhesive layer in a End Notched Flexure (ENF) specimen is presented and evaluated. The method is based on the path independent J-integral where the energy release rate (ERR) for the adhesive is derived as a function of the applied forces and the rotational displacements at the loading point and at the supports of the specimen. The major advantage with the method in comparison with existing theory known by the job initiator is that it is still applicable with ENF specimens that are subjected to yielding of the adherends.</p><p>The structure of this thesis is disposed so that the theory behind the J-integral method is shortly described and then an evaluation of the method is performed by aid of finite element simulations using beam and cohesive elements. The finite element simulations indicates that the ERR can be determined with good accuracy for an ENF specimen where a small scale yielding of the adherends has occurred. However when a fully cross sectional yielding of the adherends is reached the ERR starts diverging from the exact value and generates a too high ERR according to input data in the simulations, i.e. the exact values. The importance in length of the adhesive process zone is also shown to be irrelevant to the ERR measured according to the J-integral method.</p><p>Simulation performed with continuum elements indicates that a more reality based FE- simulation implies a higher value of the applied load in order to create crack propagation. This is an effect of that the specimen is allowed to roll on the supports which makes the effective length between the supports shorter than the initial value when the specimen is deformed. This results in a stiffer specimen and thus a higher applied force is needed to create crack propagation in the adhesive layer.</p><p>An experimental set up of an ENF specimen is created and the sample data from the experiments are evaluated with the J-integral method. For measuring the rotational displacements of the specimen which are needed for the J-integral equation an image system is developed by the author and validated by use of linear elastic beam theory. The system calculates the three rotational displacements of the specimen by aid of images taken by a high resolution SLR camera and the system for measuring the rotations may also be used in other applications than for a specific ENF geometry. The validation of the image system shows that the rotations calculated by the image system diverge from beam theory with less than 2.2 % which is a quite good accuracy in comparison with the accuracies for the rest of the used surveying equipment.</p><p>The results from the experiment indicates that the used, about 0.36 mm thick SikaPower 498, adhesive has an maximum shear strength of 37.3 MPa and a critical shear deformation of 482 µm. The fracture energy is for this thickness of the adhesive is determined as 12.9 kJ/m².</p><p>This report ends with a conclusion- and a suggested future work- chapter.</p>

Over the last few decades, Shape Memory Alloys (SMAs) have been increasingly explored in order to take advantage of their unique properties (i.e., pseudoelasticity and shape memory effect), in various actuation, sensing and absorption applications. In order to achieve an effective design of SMA-based devices a thorough investigation of their behavior in the presence of cracks is needed. In particular, it is important to understand the effect of phase transformation on their fracture response. The aim of the present work is to study the effect of stress-induced as well as thermo-mechanically-induced phase transformation on several characteristics of the fracture response of SMAs. The SMA thermomechanical response is modeled through an existing constitutive phenomenological model, developed within the framework of continuum thermodynamics, which has been implemented in a finite element frame-work. The effect of stress-induced phase transformation on the mechanical fields in the vicinity of a stationary crack and on the toughness enhancement associated with crack advance in an SMA subjected to in-plane mode I loading conditions is examined. The small scale transformation assumption is employed in the analysis according to which the size of the region occupied by the transformed material forming close to the crack tip is small compared to any characteristic length of the problem (i.e. the size of the transformation zone is thirty times smaller than the size of the cracked ligament). Given this assumption, displacement boundary conditions, corresponding to the Irwin’s solution for linear elastic fracture mechanics, are applied on a circular region in the austenitic phase that encloses the stress-induced phase transformation zone. The quasi-static stable crack growth is studied by assuming that the crackpropagates at a certain critical level of the crack-tip energy release rate. The Virtual Crack Closure Technique (VCCT) is employed to calculate the energy release rate. Fracture toughness enhancement associated with transformation dissipation is observed and its sensitivity on the variation of key characteristic non-dimensional parameters related to the constitutive response is investigated. Moreover, the effect of the dissipation due plastic deformation on the fracture resistance is analyzed by using a Cohesive Zone Model (CZM). The effect of thermo-mechanically-induced transformation on the driving force for crack growth is analyzed in an infinite center-cracked SMA plate subjected to thermal actuation under isobaric mode I loading. The crack-tip energy release rate is identified as the driving force for crack growth and is measured over the entire thermal cycle by means of the VCCT. A substantial increase of the crack-tip energy release rate – an order of magnitude for some material systems – is observed during actuation as a result of phase transformation, i.e., martensitic transformation occurring during actuation causes anti-shielding that might cause the energy release rate to reach the critical value for crack growth. A strong dependence of the crack-tip energy release rate on the variation of the thermomechanical parameters characterizing the material response is examined. Therefore, it is implied that the actual shape of the strain- temperature curve is important for the quantitative determination of the change of the crack-tip energy release rate during actuation.