A model for crack formation during active solid pyrolysis of a char-forming solid: crack patterns; surface area generation; volatile mass efflux

A model is developed for the formation and propagation of cracks in a material sample that is heated at its top surface, pyrolyses, and then thermally degrades to form char. In this work the sample is heated uniformly over its entire top surface by a hypothetical flame (a heat source). The pyrolysis mechanism is described by a one-step overall reaction that is dependent nonlinearly on the temperature (Arrhenius form). Stresses develop in response to the thermal degradation of the material by means of a shrinkage strain caused by local mass loss during pyrolysis. When the principal stress exceeds a prescribed threshold value, the material forms a local crack. Cracks are found to generally originate at the surface in response to heating, but occasionally they form in the bulk, away from ever-changing material boundaries. The resulting cracks evolve and form patterns whose characteristics are described. Quantities examined in detail are: the crack spacing in the pyrolysis zone; the crack length evolution; the formation and nature of crack loops which are defined as individual cracks that have joined to form loops that are disconnected from the remaining material; the formation of enhanced pyrolysis area; and the impact of all of the former quantities on mass flux. It is determined that the mass flux from the sample can be greatly enhanced over its nominal (non-cracking) counterpart. The mass efflux profile qualitatively resembles those observed in Cone Calorimeter tests.     

Closed-form solutions for a rectangular layer with central crack under anti-plane shear stress

We consider the problem of determining the stress distributionin a finite rectangular elastic layer containing a Griffithcrack which is opened by internal shear stress acting alongthe length of the crack. The mode III crack is assumed to belocated in the middle plane of the rectangular layer. The followingtwo problems are considered: (A) the central crack is perpendicularto the two fixed lateral surfaces and parallel to the othertwo stress-free surfaces; (B) all the lateral surfaces of therectangular layer are clamped and the central crack is parallelto the two lateral surfaces. By using Fourier transformations,we reduce the solution of each problem to the solution of dualintegral equations with sine kernels and a weight function whichare solved exactly. Finally, we derive closed-form expressionsfor the stress intensity factor at the tip of the crack andthe numerical values for the stress intensity factor at theedges of the cracks are presented in the form of tables.     

In-Situ Intrinsic Interface Strength Measurement at Elevated Temperatures and its Relationship to Interfacial Structure

A previously developed laser spallation technique has been modified to measure the tensile strength of thin film interfaces in-situ at temperatures up to 1100°C. Tensile strengths of Nb/A-plane sapphire, FeCrAl/A-plane sapphire and FeCrAlY/A-plane sapphire were measured up to 950°C. The measured strengths at high temperatures were substantially lower compared with their corresponding strengths at ambient temperature. For example, at 850°C, the interface tensile strength for the Nb/sapphire (151 ± 17 MPa), FeCrAl/sapphire (62 ± 8 MPa) and FeCrAlY/sapphire (82 ± 11 MPa) interface systems were lower by factors of approximately, 3, 5, and 8, respectively, over their corresponding ambient values. These results underscore the importance of using such in-situ measured values under operating conditions as the failure criterion in any life prediction or reliability models of such coated systems where local interface temperature excursions are expected. The results on alloy film interfaces also demonstrate that the presence of Y increases the strength of FeCrAl/Al₂O₃ interfaces.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of the Lazarus-Leblond constants in the asymptotic model of the interfacial wavy crack

The paper addresses the problem of a semi-infinite plane crack along the interface between two isotropic half-spaces. Two methods of solution have been considered in the past: Lazarus and Leblond [1998a. Three-dimensional crack-face weight functions for the semi-infinite interface crack-I: variation of the stress intensity factors due to some small perturbation of the crack front. J. Mech. Phys. Solids 46, 489-511, 1998b. Three-dimensional crack-face weight functions for the semi-infinite interface crack-II: integrodifferential equations on the weight functions and resolution J. Mech. Phys. Solids 46, 513-536] applied the “special” method by Bueckner [1987. Weight functions and fundamental fields for the penny-shaped and the half-plane crack in three space. Int. J. Solids Struct. 23, 57-93] and found the expression of the variation of the stress intensity factors for a wavy crack without solving the complete elasticity problem; their solution is expressed in terms of the physical variables, and it involves five constants whose analytical representation was unknown; on the other hand, the “general” solution to the problem has been recently addressed by Bercial-Velez et al. [2005. High-order asymptotics and perturbation problems for 3D interfacial cracks. J. Mech. Phys. Solids 53, 1128-1162], using a Wiener-Hopf analysis and singular asymptotics near the crack front.The main goal of the present paper is to complete the solution to the problem by providing the connection between the two methods. This is done by constructing an integral representation for Lazarus-Leblond's weight functions and by deriving the closed form representations of Lazarus-Leblond's constants.     

Eigenvalues of the antiplane-shear crack problem for a power-law material

This paper discusses the problem of finding the eigenvalue spectrum in determining the stress and strain fields at the tip of an antiplane-shear crack in a power-law material. It is shown that the perturbation method provides an analytical dependence of the eigenvalue on the material nonlinearity parameter and the eigenvalue of the linear problem. Thus, it is possible to find the entire spectrum of eigenvalues and not only the eigenvalue of the Hutchinson-Rice-Rosengren problem. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 173–180, January–February, 2008.     

Non-elliptic deformation field near the tip of a mixed-mode crack in a compressible hyperelastic material

This paper gives an asymptotic analysis of the deformation field near the tip of an arbitrary mixed-mode crack in a compressible hyperelastic harmonic material which loses ellipticity at sufficiently large deformations. It is found that the near-tip deformation field is characterized by a localized non-elliptic deformation band issuing from the crack-tip and bounded by two curves of discontinuous deformation gradient. Explicit expression for the near-tip deformation field is obtained both inside and outside the localized deformation band. In particular, a simple relation is derived that determines the orientation of the deformation band in terms of two complex governing parameters of the near-tip fields inside and outside the deformation band, respectively.     

Thermoelastoplastic deformation of a thick-walled cylinder with a radial crack

The thermoelastoplastic fracture mechanics problem of a thick-walled cylinder subjected to internal pressure and a nonuniform temperature field is solved by the method of elastic solutions combined with the finite-element method. The correctness of the solution is provided by using the Barenblatt crack model, in which the stress and strain fields are regular. The elastoplastic problem of a cracked cylinder subjected to internal pressure and a nonuniform temperature field are solved. The calculation results are compared with available data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 173–183, May–June, 2008.     

Development of a finite-element analysis methodology for the propagation of delaminations in composite structures

Analysing the collapse of skin-stiffened structures requires capturing the critical phenomenon of skin-stiffener separation, which can be considered analogous to interlaminar cracking. This paper presents the development of a numerical approach for simulating the propagation of interlaminar cracks in composite structures. A degradation methodology was introduced in MSC.Marc, which involved the modelling of a structure with shell layers connected by user-defined multiple-point constraints (MPCs). User subroutines were written that employ the virtual crack closure technique (VCCT) to determine the onset of crack growth and modify the properties of the user-defined MPCs to simulate crack propagation. Methodologies for the release of failing MPCs are presented and are discussed with reference to the VCCT assumption of self-similar crack growth. The numerical results obtained by using the release methodologies are then compared with experimental data for a double-cantilever beam specimen. Based on this comparison, recommendations for the future development of the degradation model are made, especially with reference to developing an approach for the collapse analysis of fuselage-representative structures. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 1, pp. 15–42, January–February, 2007.     

Measurement of unbound cocaine in blood, brain and bile of anesthetized rats using microdialysis coupled with liquid chromatography and verified by tandem mass spectrometry

To investigate the disposition of unbound cocaine in the rat blood, brain and bile, we demonstrate an in vivo multiple sampling microdialysis system coupled with liquid chromatography for cocaine assay and verified by tandem mass spectrometry. Three microdialysis probes were concurrently inserted into the jugular vein, bile duct and brain striatum of each anesthetized rat. After a period of 2 h post-surgical stabilization, cocaine (10 mg kg(-1)) was administered through the femoral vein. Separation of unbound cocaine from various biological fluids was applied to a reversed-phase C(18) column (250 x 4.6 mm I.D., 5 microm). The mobile phase consisted of acetonitrile--10 mm potassium dihydrogen phosphate buffer (25:75, v/v, pH 4.0) and 0.8% diethylamine at a flow rate of 1 mL min(-1). The UV detector wavelength was set at 235 nm. The results indicate that cocaine penetrates the blood--brain barrier with a rapid distribution. However, unbound cocaine in the bile dialysate was not detectable in the UV detection. We therefore use LC--tandem mass spectrometry to detect the bile fluid after cocaine administration (3 mg kg(-1), i.v.). The results indicate that cocaine goes through hepatobiliary excretion.