Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

Evolution-induced catastrophe and its predictability

Both earthquake prediction and failure prediction of disordered brittle media are difficult and complicated problems and they might have something in common. In order to search for clues for earthquake prediction, the common features of failure in a simple nonlinear dynamical model resembling disordered brittle media are examined. It is found that the failure manifests evolution-induced catastrophe (EIC), i.e., the abrupt transition from globally stable (GS) accumulation of damage to catastrophic failure. A distinct feature is the significant uncertainty of catastrophe, called sample-specificity. Consequently, it is impossible to make a deterministic prediction macroscopically. This is similar to the question of predictability of earthquakes. However, our model shows that strong stress fluctuations may be an immediate precursor of catastrophic failure statistically. This might provide clues for earthquake forecasting.

Evolution of three-dimensional coherent structures in compressible axisymmetric jet

A high order difference scheme is used to simulate the spatially developing compressible axisymmetric jet. The results show that the Kelvin-Helmholtz instability appears first when the jet loses its stability, and then with development of jet the increase in nonlinear effects leads to the secondary instability and the formation of the streamwise vortices. The evolution of the three-dimensional coherent structure is presented. The computed results verify that in axisymmetric jet the secondary instability and formation of the streamwise vortices are the important physical mechanism of enhancing the flow mixing and transition occurring.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Flood Routing Models in Confluent and Dividing Channels

By introducing a water depth connecting formula, the hydraulic equations in the dividing channel system were coupled and the relation of discharge distribution between the branches of the dividing channels can be yielded. In this manner, a numerical model for the confluent channels was established to study the variation of backwater effects with the parameters in the channel junction. The meeting of flood peaks in the mainstream and tributary can be analyzed with this model. The flood peak meeting is found to be a major factor for the extremely high water level in the mainstream during the 1998 Yangtze River flood. Subsequently the variations of discharge distribution and water level with channel parameters between each branch in this system were studied as well. As a result, flood evolution caused by Jingjiang River shortcut and sediment deposition in the entrance of dividing channels of the Yangtze River may be qualitatively elucidated. It is suggested to be an effective measure for flood mitigation to enhance regulation capability of reservoirs available upstream of the tributaries and harness branch entrance channels.

The Evolution of Fatigue bocks in Aluminum Alloy LY-12 under Shock beading

对LY12铝合金在低周疲劳条件下的裂纹情况和随后进行的动态拉伸条件下裂纹的发展给予了观察和统计分析。发现裂纹的累积数密度分布在损伤演化过程中保持指数形式,用NAG模型对实验结果进行分析,得出该材料裂纹演化发展方程的各种参数。

A further discussion on basic equations for two-phase flow: On collision stress of solid phase

The following points are argued: (i) there are two independent kinds of interaction on interfaces, i.e. the interaction between phases and the collision interaction, and the jump relations on interfaces can accordingly be resolved; (ii) the stress in a particle can also be divided into background stress and collision stress corresponding to the two kinds of interaction on interfaces respectively; (iii) the collision stress, in fact, has no jump on interface, so the averaged value of its derivative is equal to the derivative of its averaged value; (iv) the stress of solid phase in the basic equations for two\|phase flow should include the collision stress, while the stress in the expression of the inter\|phase force contains the background one only. Based on the arguments, the strict method for deriving the equations for two\|phase flow developed by Drew, Ishii et al. is generalized to the dense two\|phase flow, which involves the effect of collision stress.

On the Evolution of Simple Shear in Saturated Soils

Development of shear bands in saturated soils is a multi-stage process based on the theoretical and numerical investigations in this paper. The soil is initially in homogenous shear strain state, and the instability can be characterized by a dimensionless number D. The inhomogenous distribution of shear strains appears when D>1, and the shear band will initiate and develop gradually. Numerical solutions show that only single shear band that is finally formed in the central region of the specimen even several disturbances (distributed along the specimen) appear in the beginning.

Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Phase Evolution and its Effect on Magnetic Properties of Nd60Al10Fe20Co10 Bulk Metallic Glass

The thermal stability of nanocrystalline clusters, the phase evolution, and their effects on magnetic Propertieswere studied for as-cast Nd60Al10Fe20Co10 alloy using differential scanning calorimetry curves, x-ray diffraction patterns, scanning electron microscopy, and high-resolution transition electron microscopy. Thermomagnetic curves and hysteresis loops of the bulk metallic glass were measured during the annealing process. The high thermostability of the hardmagnetic properties of the samples observed is attributed to the stability of the nanocrystalline clusters upon annealing, while the slight enhancement in the magnetization is due to the precipitation of some Nd-rich metastable phases. The mechanism of thermostability of the nanocrystalline clusters and the formation of the metastable phases are discussed.

Numerical simulation of temperature field evolution of metallic alloy surface due to pulsed laser treatment

This paper performed a numerical simulation on temperature field evolution for the surface layer of a metallic alloy subjected to pulsed Nd:YAG laser treatment. The enthalpy method was adopted to solve the moving boundary problem, I.e. Stefan problem. Computational results were obtained to show the temperature field evolution. Effects of latent heat and mushy zone width on the temperature field were investigated. The results also show very high values of temperature gradient and cooling rate, which are typical characteristics during the solidification process.