A probabilistic model for elastic-plastic fracture analysis of plane and axisymmetric structures

This work presents a methodology for elastic-plastic fracture reliability analysis of plane and axisymmetric structures. The structural reliability analysis is accomplished by means of the FORM analytical method. The virtual crack extension technique based on a direct minimization of potencial energy is utililized for the calculation of the energy release rate. Results are presented to illustrate the performance of the adopted methodology.

The rolling ball problem on the plane revisited

By a sequence of rollings without slipping or twisting along segments of a straight line of the plane, a spherical ball of unit radius has to be transferred from an initial state to an arbitrary final state taking into account the orientation of the ball. We provide a new proof that with at most 3 moves, we can go from a given initial state to an arbitrary final state. The first proof of this result is due to Hammersley ( 1983). His proof is more algebraic than ours which is more geometric. We also showed that generically no one of the three moves, in any elimination of the spin discrepancy, may have length equal to an integral multiple of 2 pi.

Effect of sample pre-cracking method and notch geometry in plane strain fracture toughness tests as applied to a PMMA resin

The influence of test method factors (notch shape, square or angular, and pre-cracking method, by tapping onto or pressing a razor blade) on the results obtained in plane strain fracture toughness test according to standard ASTM D5045 using SENB specimens made of a commercial PMMA resin were investigated. Results were analyzed quantitatively by comparing the obtained K-IC values and qualitatively by observing their effect on the Moire fringes observed using photoelasticity, showing that, at 95% significance level, the K-IC values are affected by the pre-cracking method, with the most conservative value being obtained when natural pre-cracks were introduced by tapping onto a razor blade (K-IC = 1.15 +/- 0.11 MPa.m(0.5)). This correlates with a perturbation in the stress field close to the pre-crack tip observed in the photoelasticity test sample when it was introduced by pressing the razor blade. Surprisingly, notch geometry only slightly affects the results. (C) 2012 Elsevier Ltd. All rights reserved.

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

Tight frames of k-plane ridgelets and the problem of representing objects that are smooth away from d-dimensional singularities in Rn

For each pair (n, k) with 1 ≤ k < n, we construct a tight frame (ρλ : λ ∈ Λ) for L2 (Rn), which we call a frame of k-plane ridgelets. The intent is to efficiently represent functions that are smooth away from singularities along k-planes in Rn. We also develop tools to help decide whether k-plane ridgelets provide the desired efficient representation. We first construct a wavelet-like tight frame on the X-ray bundle χn,k—the fiber bundle having the Grassman manifold Gn,k of k-planes in Rn for base space, and for fibers the orthocomplements of those planes. This wavelet-like tight frame is the pushout to χn,k, via the smooth local coordinates of Gn,k, of an orthonormal basis of tensor Meyer wavelets on Euclidean space Rk(n−k) × Rn−k. We then use the X-ray isometry [Solmon, D. C. (1976) J. Math. Anal. Appl. 56, 61–83] to map this tight frame isometrically to a tight frame for L2(Rn)—the k-plane ridgelets. This construction makes analysis of a function f ∈ L2(Rn) by k-plane ridgelets identical to the analysis of the k-plane X-ray transform of f by an appropriate wavelet-like system for χn,k. As wavelets are typically effective at representing point singularities, it may be expected that these new systems will be effective at representing objects whose k-plane X-ray transform has a point singularity. Objects with discontinuities across hyperplanes are of this form, for k = n − 1.

On the numerical solution of a Cauchy problem in an elastostatic half-plane with a bounded inclusion

We propose an iterative procedure for the inverse problem of determining the displacement vector on the boundary of a bounded planar inclusion given the displacement and stress fields on an infinite (planar) line-segment. At each iteration step mixed boundary value problems in an elastostatic half-plane containing the bounded inclusion are solved. For efficient numerical implementation of the procedure these mixed problems are reduced to integral equations over the bounded inclusion. Well-posedness and numerical solution of these boundary integral equations are presented, and a proof of convergence of the procedure for the inverse problem to the original solution is given. Numerical investigations are presented both for the direct and inverse problems, and these results show in particular that the displacement vector on the boundary of the inclusion can be found in an accurate and stable way with small computational cost.

Magnetoconfined levels in a parabolic quantum dot: An analytical solution of a three-dimensional Fock-Darwin problem in a tilted magnetic field

The energy spectrum of an electron confined in a quantum dot (QD) with a three-dimensional anisotropic parabolic potential in a tilted magnetic field was found analytically. The theory describes exactly the mixing of in-plane and out-of-plane motions of an electron caused by a tilted magnetic field, which could be seen, for example, in the level anticrossing. For charged QDs in a tilted magnetic field we predict three strong resonant lines in the far-infrared-absorption spectra.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

A contribution to the exact modal solution of in-plane beam structures

The exact vibration modes and natural frequencies of planar structures and mechanisms, comprised Euler-Bernoulli beams, are obtained by solving a transcendental. nonlinear, eigenvalue problem stated by the dynamic stiffness matrix (DSM). To solve this kind of problem, the most employed technique is the Wittrick-Williams algorithm, developed in the early seventies. By formulating a new type of eigenvalue problem, which preserves the internal degrees-of-freedom for all members in the model, the present study offers an alternative to the use of this algorithm. The new proposed eigenvalue problem presents no poles, so the roots of the problem can be found by any suitable iterative numerical method. By avoiding a standard formulation for the DSM, the local mode shapes are directly calculated and any extension to the beam theory can be easily incorporated. It is shown that the method here adopted leads to exact solutions, as confirmed by various examples. Extensions of the formulation are also given, where rotary inertia, end release, skewed edges and rigid offsets are all included. (C) 2008 Elsevier Ltd. All rights reserved.

Effect of alpha prime due to 475 degrees C aging on fracture behavior and corrosion resistance of DIN 1.4575 and MA 956 high performance ferritic stainless steels

The 475 degrees C embrittlement in stainless steels is a well-known phenomenon associated to alpha prime (alpha`) formed by precipitation or spinodal decomposition. Many doubts still remain on the mechanism of alpha` formation and its consequence on deformation and fracture mechanisms and corrosion resistance. In this investigation, the fracture behavior and corrosion resistance of two high performance ferritic stainless steels were investigated: a superferritic DIN 1.4575 and MA 956 superalloy were evaluated. Samples of both stainless steels (SS) were aged at 475 degrees C for periods varying from 1 to 1,080 h. Their fracture surfaces were observed using scanning electron microscopy (SEM) and the cleavage planes were determined by electron backscattering diffraction (EBSD). Some samples were tested for corrosion resistance using electrochemical impedance spectroscopy (EIS) and potentiodynamic polarization. Brittle and ductile fractures were observed in both ferritic stainless steels after aging at 475 degrees C. For aging periods longer than 500 h, the ductile fracture regions completely disappeared. The cleavage plane in the DIN 1.4575 samples aged at 475 degrees C for 1,080 h was mainly {110}, however the {102}, {314}, and {131} families of planes were also detected. The pitting corrosion resistance decreased with aging at 475 degrees C. The effect of alpha prime on the corrosion resistance was more significant in the DIN 1.4575 SS comparatively to the Incoloy MA 956.

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.

Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

Factors associated with pelvic asymmetry in transverse plane during gait in patients with cerebral palsy

The purpose of this study was to describe the patterns of pelvic rotational asymmetry in the transverse plane and identify the possible factors related to this problem. One thousand and forty-five patients with cerebral palsy (CP) and complete documentation in the gait laboratory were reviewed in a retrospective study. Pelvic asymmetry in the transverse plane was observed in 52.7% of the patients; and to identify the possible causes of pelvic retraction, clinical (Thomas test, popliteal angle, and gastrocnemius tightness) and dynamic parameters (mean rotation of the hip in stance, minimum hip flexion, minimum knee flexion, and peak ankle dorsiflexion) were evaluated. The association between these parameters and pelvic retraction was assessed statistically. The results showed that 75.7% of patients with asymmetric pattern of the pelvis had clinical diagnosis of diplegic spastic CP. Among the patients with asymmetrical CP, the most common pattern was pelvic retraction on the affected side. The relationship between pelvic retraction and internal hip rotation was stronger in patients with asymmetrical diplegic CP than in those with hemiplegic (P<0.001) or symmetrical diplegic CP (P=0.014). All of the patients exhibited a significant association among clinical parameters (Thomas test, popliteal angle, and gastrocnemius tightness) and pelvic retraction. In conclusion, pelvic retraction seems to be a multifactorial problem, and the etiology can change according to topographic classification, which must be taken into account during the decision-making process in patients with CP. J Pediatr Orthop B 18:320-324 (C) 2009 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins.

IN VITRO FRACTURE RESISTANCE OF GLASS-FIBER AND CAST METAL POSTS WITH DIFFERENT LENGTHS

Statement of problem. Dental fractures can occur in endodontically treated teeth restored with posts. Purpose. The purpose of this study was to evaluate the in vitro fracture resistance of roots with glass-fiber and metal posts of different lengths. Material and methods. Sixty endodontically treated maxillary canines were embedded in acrylic resin, except for 4 mm of the cervical area, after removing the clinical crowns. The post spaces were opened with a cylindrical bur at low speed attached to a surveyor, resulting in preparations with lengths of 6 mm (group 6 mm), 8 mm (group 8 mm), or 10 mm (group 10 mm). Each group was divided into 2 subgroups according to the post material: cast post and core or glass-fiber post (n=30). The posts were luted with dual-polymerizing resin cement (Panavia F). Cast posts and cores of Co-Cr (Resilient Plus) crowns were made and cemented with zinc phosphate. Specimens were subjected to increasing compressive load (N) until fracture. Data were analyzed with 2-way ANOVA and the Tukey-Kramer test (alpha=.05). Results. The ANOVA analysis indicated significant differences (P<.05) among the groups, and the Tukey test revealed no significant difference among the metal posts of 6-mm length (26.5 N +/- 13.4), 8-mm length (25.2 N +/- 13.9), and 10-mm length (17.1 N +/- 5.2). Also, in the glass-fiber post group, there was no significant difference when posts of 8-mm length (13.4 N +/- 11.0) were compared with the 6-mm (6.9 N +/- 4.6) and 10-mm (31.7 N +/- 13.1) groups. The 10-mm-long post displayed superior fracture resistance, and the 6-mm-long post showed significantly lower mean values (P<.001). Conclusions. Within the limitations of this study, it was concluded that the glass-fiber post represents a viable alternative to the cast metal post, increasing the resistance to fracture of endodontically treated canines. (J Prosthet Dent 2009;101:183-188)

Rock fracture mean trace length estimation and confidence interval calculation using maximum likelihood methods

There has been a resurgence of interest in the mean trace length estimator of Pahl for window sampling of traces. The estimator has been dealt with by Mauldon and Zhang and Einstein in recent publications. The estimator is a very useful one in that it is non-parametric. However, despite some discussion regarding the statistical distribution of the estimator, none of the recent works or the original work by Pahl provide a rigorous basis for the determination a confidence interval for the estimator or a confidence region for the estimator and the corresponding estimator of trace spatial intensity in the sampling window. This paper shows, by consideration of a simplified version of the problem but without loss of generality, that the estimator is in fact the maximum likelihood estimator (MLE) and that it can be considered essentially unbiased. As the MLE, it possesses the least variance of all estimators and confidence intervals or regions should therefore be available through application of classical ML theory. It is shown that valid confidence intervals can in fact be determined. The results of the work and the calculations of the confidence intervals are illustrated by example. (C) 2003 Elsevier Science Ltd. All rights reserved.