Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Cellular Cell Bifurcation Of Cylindrical Detonations

Cellular cell pattern evolution of cylindrically-diverging detonations is numerically simulated successfully by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. From the simulation, three cell bifurcation modes are observed during the evolution and referred to as concave front focusing, kinked and wrinkled wave front instability, and self-merging of cellular cells. Numerical research demonstrates that the wave front expansion resulted from detonation front diverging plays a major role in the cellular cell bifurcation, which can disturb the nonlinearly self-sustained mechanism of detonations and finally lead to cell bifurcations.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

A theoretical analysis on stochastic behaviour of short cracks and fatigue damage of metals

In this paper, the closed form of solution to the stochastic differential equation for a fatigue crack evolution system is derived. and the relationship between metal fatigue damage and crack stochastic behaviour is investigated. It is found that the damage extent of metals is independent of crack stochastic behaviour ii the stochastic deviation of the crack growth rate is directly proportional to its mean value. The evolution of stochastic deviation of metal fatigue damage in the stage close to the transition point between short and long crack regimes is also discussed.

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.

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations (PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub-characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic, respectively. The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories, the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time, the methods of removing the remained ellipticity are further obtained from the nonlinear PSE.

Optimal bounded control of first-passage failure of strongly non-linear oscillators under combined harmonic and white-noise excitations

A procedure for designing the optimal bounded control of strongly non-linear oscillators under combined harmonic and white-noise excitations for minimizing their first-passage failure is proposed. First, a stochastic averaging method for strongly non-linear oscillators under combined harmonic and white-noise excitations using generalized harmonic functions is introduced. Then, the dynamical programming equations and their boundary and final time conditions for the control problems of maximizing reliability and of maximizing mean first-passage time are formulated from the averaged Ito equations by using the dynamical programming principle. The optimal control law is derived from the dynamical programming equations and control constraint. Finally, the conditional reliability function, the conditional probability density and mean of the first-passage time of the optimally controlled system are obtained from solving the backward Kolmogorov equation and Pontryagin equation. An example is given to illustrate the proposed procedure and the results obtained are verified by using those from digital simulation. (C) 2003 Elsevier Ltd. All rights reserved.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.