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.

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.

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.

Quantum path control in few-optical-cycle regime

We theoretically show that selection of a single quantum path in high-order harmonics generation can be realized in a few-optical-cycle regime with two-color schemes. We also demonstrate, in theory as well, the generation of spectrally smooth and ultrabroad extreme ultraviolet supercontinuum in argon gas which can produce single similar to 79 as pulses with currently available ultrafast laser sources. Our finding can be beneficial for generating isolated sub-100 as extreme ultraviolet pulses.