Heteroclinic cycles in lattice dynamical systems

A criterion of spatial chaos occurring in lattice dynamical systems-heteroclinic cycle-is discussed. It is proved that if the system has asymptotically stable heteroclinic cycle, then it has asymptotically stable homoclinic point which implies spatial chaos.

Modified Generalized Beam Lattice Model Associated With Fracture Of Reinforced Fiber/Particle Composites

Lattice-type model can simulate in a straightforward manner heterogeneous brittle media. Mohr-Coulomb failure criterion has recently been involved into the generalized beam (GB) lattice model, and as a result, numerical experiments on concrete under various loading conditions can be conducted. The GB lattice model is further used to investigate the reinforced fiber/particle composites instead of only particle composites as the model did before. Numerical examples are given to show the effectiveness of the modeling procedure, and influences of inclusions (particle, fiber and rebar) on the fracture processes are also discussed. (c) 2008 Elsevier Ltd. All rights reserved.

Designing Of Partial Similarity Models And Evaluation Method In Polymer Flooding Experiment

Based on the scaling criteria of polymer flooding reservoir obtained in our previous work in which the gravity and capillary forces, compressibility, non-Newtonian behavior, absorption, dispersion, and diffusion are considered, eight partial similarity models are designed. A new numerical approach of sensitivity analysis is suggested to quantify the dominance degree of relaxed dimensionless parameters for partial similarity model. The sensitivity factor quantifying the dominance degree of relaxed dimensionless parameter is defined. By solving the dimensionless governing equations including all dimensionless parameters, the sensitivity factor of each relaxed dimensionless parameter is calculated for each partial similarity model; thus, the dominance degree of the relaxed one is quantitatively determined. Based on the sensitivity analysis, the effect coefficient of partial similarity model is defined as the summation of product of sensitivity factor of relaxed dimensionless parameter and its relative relaxation quantity. The effect coefficient is used as a criterion to evaluate each partial similarity model. Then the partial similarity model with the smallest effect coefficient can be singled out to approximate to the prototype. Results show that the precision of partial similarity model is not only determined by the number of satisfied dimensionless parameters but also the relative relaxation quantity of the relaxed ones.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

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.

Peeling Experiments Of Ductile Thin Films Along Ceramic Substrates - Critical Assessment Of Analytical Models

Two types of peeling experiments are performed in the present research. One is for the Al film/Al2O3 substrate system with an adhesive layer between the film and the substrate. The other one is for the Cu film/Al2O3 substrate system without adhesive layer between the film and the substrate, and the Cu films are electroplated onto the Al2O3 substrates. For the case with adhesive layer, two kinds of adhesives are selected, which are all the mixtures of epoxy and polyimide with mass ratios 1:1.5 and 1:1, respectively. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling process. The effects of the adhesive layer on the energy release rate are analyzed. Using the experimental results, several analytical criteria for the steady-state peeling based on the bending model and on the two-dimensional finite element analysis model are critically assessed. Through assessment of analytical models, we find that the cohesive zone criterion based on the beam bend model is suitable for a weak interface strength case and it describes a macroscale fracture process zone case, while the two-dimensional finite element model is effective to both the strong interface and weak interface, and it describes a small-scale fracture process zone case. (C) 2007 Elsevier Ltd. All rights reserved.

On The Thermally Induced Interfacial Delamination Of A Segmented Coating

The thermally induced interfacial delamination problem of a segmented coating is investigated using finite element method (FEM). The coating-substrate system, modeled as a coated semi-infinite medium with periodic segmentation cracks within coating, is assumed to be exposed to convective cooling from surface. The failure criterion based on the interfacial fracture toughness is adopted, in which the energy release rate for an interface crack is considered to be the driving force for interfacial delamination extension. The results confirm that a segmented coating has higher delamination resistance than an intact one under the same thermal transients, as the segmentation crack spacing is smaller than a critical value. Based on dimensional analysis, sensitivity analyses of the crack driving force are also obtained as a function of various dimensionless parameters such as time, convection severity and material constants. These results may provide some helpful references for the integrity of coating-substrate systems under thermal loading. (C) 2007 Elsevier B.V. All rights reserved.

Energy Dissipation In Fracture Of Bulk Metallic Glasses Via Inherent Competition Between Local Softening And Quasi-Cleavage

Compression, tension and high-velocity plate impact experiments were performed on a typical tough Zr41.2Ti13.8Cu10Ni12.5Be22.5 (Vit 1) bulk metallic glass (BMG) over a wide range of strain rates from similar to 10(-4) to 10(6) s(-1). Surprisingly, fine dimples and periodic corrugations on a nanoscale were also observed on dynamic mode I fracture surfaces of this tough Vit 1. Taking a broad overview of the fracture patterning of specimens, we proposed a criterion to assess whether the fracture of BMGs is essentially brittle or plastic. If the curvature radius of the crack tip is greater than the critical wavelength of meniscus instability [F. Spaepen, Acta Metall. 23 615 (1975); A.S. Argon and M. Salama, Mater. Sci. Eng. 23 219 (1976)], microscale vein patterns and nanoscale dimples appear on crack surfaces. However, in the opposite case, the local quasi-cleavage/separation through local atomic clusters with local softening in the background ahead of the crack tip dominates, producing nanoscale periodic corrugations. At the atomic cluster level, energy dissipation in fracture of BMGs is, therefore, determined by two competing elementary processes, viz. conventional shear transformation zones (STZs) and envisioned tension transformation zones (TTZs) ahead of the crack tip. Finally, the mechanism for the formation of nanoscale periodic corrugation is quantitatively discussed by applying the present energy dissipation mechanism.

Crack propagation in the power-law nonlinear viscoelastic material

An analysis on crack creep propagation problem of power-law nonlinear viscoelastic materials is presented. The creep incompressilility assumption is used. To simulate fracture behavior of craze region, it is assumed that in the fracture process zone near the crack tip, the cohesive stress sigma(f) acts upon the crack surfaces and resists crack opening. Through a perturbation method, i. e., by superposing the Mode-I applied force onto a referential uniform stress state, which has a trivial solution and gives no effect on the solution of the original problem, the nonlinear viscoelastic problem is reduced to linear problem. For weak nonlinear materials, for which the power-law index n similar or equal to 1, the expressions of stress and crack surface displacement are derived. Then, the fracture process zone local energy criterion is proposed and based on which the formulas of cracking incubation time t

Statistical analysis of synthetic earthquake catalogs generated by models with various levels of fault zone disorder

The stress release model, a stochastic version of the elastic rebound theory, is applied to the large events from four synthetic earthquake catalogs generated by models with various levels of disorder in distribution of fault zone strength (Ben-Zion, 1996) They include models with uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-scale heterogeneities (M). The results show that the degree of regularity or predictability in the assumed fault properties, based on both the Akaike information criterion and simulations, follows the order U, F, A, and M, which is in good agreement with that obtained by pattern recognition techniques applied to the full set of synthetic data. Data simulated from the best fitting stress release models reproduce, both visually and in distributional terms, the main features of the original catalogs. The differences in character and the quality of prediction between the four cases are shown to be dependent on two main aspects: the parameter controlling the sensitivity to departures from the mean stress level and the frequency-magnitude distribution, which differs substantially between the four cases. In particular, it is shown that the predictability of the data is strongly affected by the form of frequency-magnitude distribution, being greatly reduced if a pure Gutenburg-Richter form is assumed to hold out to high magnitudes.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Fracture statistics of brittle materials: Weibull or normal distribution

The fit of fracture strength data of brittle materials (Si3N4, SiC, and ZnO) to the Weibull and normal distributions is compared in terms of the Akaike information criterion. For Si3N4, the Weibull distribution fits the data better than the normal distribution, but for ZnO the result is just the opposite. In the case of SiC, the difference is not large enough to make a clear distinction between the two distributions. There is not sufficient evidence to show that the Weibull distribution is always preferred to other distributions, and the uncritical use of the Weibull distribution for strength data is questioned.

Influence of threshold stress on the estimation of the Weibull statistics

The influence of threshold stress on the estimation of the Weibull statistics is discussed in terms of the Akaike information criterion. Numerical simulations show that, if sample data are limited in number and threshold stress is not too large, the two-parameter Weibull distribution is still a preferred choice. For example, the fit of strength data of glass and ceramics to the two- and three-parameter Weibull distributions is compared.

Electrothermal stress in conductive body with collinear cracks

The temperature and stress field in a thin plate with collinear cracks interrupting an electric current field are determined. This is accomplished by using a complex function method that allows a direct means of finding the distribution of the electric current, the temperature and stress field. Temperature dependency for the heat-transfer coefficient, coefficient of linear expansion and the elastic modulus are considered. As an example, temperature distribution is calculated for an alloy (No. GH2132) plate with two collinear cracks under high temperature. Relationships between the stress, temperature, electric density and crack length are obtained. Crack trajectories emanating from existing crack are predicted by application of the strain energy density criterion which can also be used for finding the load carrying capacity of the cracked plate. (C) 2003 Elsevier Ltd. All rights reserved.

The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data

The stress release model, a stochastic version of the elastic-rebound theory, is applied to the historical earthquake data from three strong earthquake-prone regions of China, including North China, Southwest China, and the Taiwan seismic regions. The results show that the seismicity along a plate boundary (Taiwan) is more active than in intraplate regions (North and Southwest China). The degree of predictability or regularity of seismic events in these seismic regions, based on both the Akaike information criterion (AIC) and fitted sensitivity parameters, follows the order Taiwan, Southwest China, and North China, which is further identified by numerical simulations. (c) 2004 Elsevier Ltd. All rights reserved.