Spallation analysis with a closed trans-scale formulation damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbers De* and De, Mach number M and damage number S. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is deter-mined by Deborah number De+. Deborah number De* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbers De* and De, Mach number M and damage number S is discussed.

Improved method makes a soft landing of well path

Research and field experience have shown that well-path control is important in many cases, not only to reach the desired coordinates, but also to arrive at the well completion target from the preferred trajectory.

An elastoplastic constitutive relation of whisker-reinforced composite for meso damage: Part I - formulation

An elastoplastic constitutive relation is developed for meso damage of whisker-reinforced composites. A model is constructed that includes orientation distribution of whiskers and slip systems as well as interface and crystal sliding. Evolution of damage will be addressed. Given in Part I is the formulation while examples will be illustrated in Part II.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

A geometry model for tortuosity of flow path in porous media

A simple geometry model for tortuosity of flow path in porous media is proposed based on the assumption that some particles in a porous medium are unrestrictedly overlapped and the others are not. The proposed model is expressed as a function of porosity and there is no empirical constant in this model. The model predictions are compared with those from available correlations obtained numerically and experimentally, both of which are in agreement with each other. The present model can also give the tortuosity with a good approximation near the percolation threshold. The validity of the present tortuosity model is thus verified.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Statistical formulation and experimental determination of growth rate of micrometre cracks under impact loading

A new statistical formulation and a relevant experimental approach to determine the growth rate of microcracks were proposed. The method consists of experimental measurements and a statistical analysis' on the basis of the conservation law of number density of microcracks in phase space. As a practical example of the method, the growth rate of microcracks appearing in an aluminium alloy subjected to planar impact loading was determined to be ca. 10 mu m/mu s under a tensile stress of 1470 MPa and load duration between 0.26 mu s and 0.80 mu s.

On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

Path Planner for Intelligent Laser Surface Modification

An optimal algorithm of manufacturing path planner for intelligent laser surface modification is presented. Elements included in the optimal objective have been analyzed. A 6-D manufacture trace that satisfies the requirements of special craft and 5-axis laser processing robot system has been generated from the path planner by method of parallel section in which combinations of modification spots size with curvature of processing surfaces and modification craft parameters are considered. Related experiments have been successfully carried out with the computer integrated multifunctional laser manufacturing system.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Electron quantum path tuning and isolated attosecond pulse emission driven by a waveform-controlled multi-cycle laser field

We investigate high-order harmonic emission and isolated attosecond pulse (IAP) generation in atoms driven by a two-colour multi-cycle laser field consisting of an 800 nm pulse and an infrared laser pulse at an arbitrary wavelength. With moderate laser intensity, an IAP of similar to 220 as can be generated in helium atoms by using two-colour laser pulses of 35 fs/800 nm and 46 fs/1150 nm. The discussion based on the three-step semiclassical model, and time-frequency analysis shows a clear picture of the high-order harmonic generation in the waveform-controlled laser field which is of benefit to the generation of XUV IAP and attosecond electron pulses. When the propagation effect is included, the duration of the IAP can be shorter than 200 as, when the driving laser pulses are focused 1 mm before the gas medium with a length between 1.5 mm and 2 mm.

Tunable high-order harmonic generation and the role of the folded quantum path

We theoretically demonstrate the selective enhancement of high-order harmonic generation (HHG) in two-color laser fields consisting of a single-cycle fundamental wave (800 nm wavelength) and a multicycle subharmonic wave (2400 nm wavelength). By performing time-frequency analyses based on a single-active-electron model, we reveal that such an enhancement is a result of the modified electron trajectories in the two-color field. Furthermore, we show that selectively enhanced HHG gives rise to a bandwidth-controllable extreme ultraviolet supercontinuum in the plateau region, facilitating the generation of intense single isolated attosecond pulses.

Precipitation of nanocrystals in glasses by electron irradiation: An alternative path to form glass ceramics?

This letter demonstrates an alternative method to form gallium silicate glass ceramics using high-energy electron irradiation. Compared with glass ceramics obtained from the conventional thermal treatment method, the distribution and crystal sizes of the precipitated Ga2O3 nanoparticles are the same. An advantage of this method is that the spatial distribution of the precipitated nanoparticles can be easily controlled. However, optically active dopants Ni2+ ions do not participate in the precipitation during electron irradiation. (c) 2007 American Institute of Physics.