High-order asymptotic analysis for the crack in nonlinear material

Accurate high-order asymptotic analyses were carried out for Mode II plane strain crack in power hardening materials. The second-order crack tip fields have been obtained. It is found that the amplitude coefficientk ₂ of the second term of the asymptotic field is correlated to the first order field as the hardening exponentn<n ^* (n ^*≈5), but asn≥n ^*,k ₂ turns to become an independent parameter. Our results also indicated that, the second term of the asymptotic field has little influence on the near-crack-tip field and can be neglected whenn<n ^*. In fact,k ₂ directly reflects the effects of triaxiality near the crack tip, the crack geometry and the loading mode, so that besidesJ-integral it can be used as another characteristic parameter in the two-parameter criterion. The project supported by National Natural Science Foundation of China     

Singular fields of a mode III interface crack in a power-law hardening bimaterial

A high order of asymptotic solution of the singular fields near the tip of a mode III interface crack for pure power-law hardening bimaterials is obtained by using the hodograph transformation. It is found that the zero order of the asymptotic solution corresponds to the assumption of a rigid substrate at the interface, and the first order of it is deduced in order to satisfy completely two continuity conditions of the stress and displacement across the interface in the asymptotic sense. The singularities of stress and strain of the zeroth order asymptotic solutions are −1/(n ₁+1) and −n/(n ₁+1) respectively. (n=n ₁,n ₂ is the hardening exponent of the bimaterials.) The applicability conditions of the asymptotic solutions are determined for both zeroth and first orders. It is proved that the Guo-Keer solution^[10] is limited in some conditions. The angular functions of the singular fields for this interface crack problem are first expressed by closed form. The project supported by National Natural Science Foundation of China     

A class of three-level explicit schemes with higher stability properties for a dispersive equation ut=auxxx

In this paper, a class of three level explicit schemes for a dispersive equation u_t=au_xxx with stability condition |r|=|α|Δt/(Δx)³≤2.382484, are considered. The stability condition for this class of schemes is much better than |r|≤0.3849 in [1], [2] and |r|≤0.701659 in [3], and |r|≤1.1851 in [4].     

The Decay Rate and Higher Approximation of Mild Solutions to the Navier�CStokes Equations

In this paper, we consider v(t) = u(t) − e ^tΔ u ₀, where u(t) is the mild solution of the Navier–Stokes equations with the initial data u₀ ? L²(\mathbb Rⁿ)?Lⁿ(\mathbb Rⁿ){u_0\in L^2({\mathbb R}^n)\cap L^n({\mathbb R}^n)} . We shall show that the L ² norm of D ^β v(t) decays like t^{-\frac |b|-1 2-\frac n4}{t^{-\frac {|\beta|-1} {2}-\frac n4}} for |β| ≥ 0. Moreover, we will find the asymptotic profile u ₁(t) such that the L ² norm of D ^β (v(t) − u ₁(t)) decays faster for 3 ≤ n ≤ 5 and |β| ≥ 0. Besides, higher-order asymptotics of v(t) are deduced under some assumptions.     

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

Some bounded results of θ(t)-type singular integral operators

Let B be a Banach space in UMD with an unconditional basis. The boundedness of the θ(t)-type singular integral operators in L _B ^p (R ⁿ), (1≤p<+∞) and H _B ¹ (R ⁿ) spaces are discussed. Foundation item: the Education Commission of Shandong Province (J98P51) Biography: Zhao Kai (1960-)     

Pseudo plane stress mode II crack growth in plastic materials

The near tip field of mode II crack that grows in thin bodies with power hardening or perfectly plastic behavior is analyzed. It is shown that for power hardening behavior, the pseudo plane stress field possesses the logarithm singularity, i.e. σ (ln r)²/(n−1), (ln r)^{2n/(n − 1)}, where r is the distance from the crack tip, n the hardening exponent is σⁿ. When n → ∞ the solution reduced to that for the perfectly plastic case.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

On Saint-Venant's Principle for a Homogeneous Elastic Arch-Like Region

In this paper, we consider a two-dimensional homogeneous isotropic elastic material state in the arch-like region a ≤ r ≤ b, 0 ≤ θ ≤ α, where (r, θ) denote plane polar coordinates. We assume that three of the edges r = a, r = b, θ = α are traction-free, while the edge θ = 0 is subjected to an (in plane) self-equilibrated load. We define an appropriate measure for the Airy stress function φ and then we establish a clear relationship with the Saint-Venant's principle on such regions. We introduce a cross-sectional integral function I(θ) which is shown to be a convex function and satisfies a second-order differential inequality. Consequently, we establish a version of the Saint-Venant principle for such a curvilinear strip, without requiring of any condition upon the dimensions of the arch-like region.     

Partial Regularity for Weak Heat Flows into a General Compact Riemannian Manifold

In this paper the partial regularity of the weak heat flow of harmonic maps from a Riemannian manifold M into a general compact Riemannian manifold N without boundary is considered. Partial results have been obtained for target manifolds that are spheres [12, 4] or homogeneous spaces [6]. The proofs in these special cases relied heavily on the geometry of these manifolds, and cannot be applied to the general case. We prove in this article that the singular set Sing(u) of the stationary weak heat flow satisfies H ⁿ _ρ (Sing(u))=0, with n=dimension M, where H ⁿ _ρ is the Hausdorff measure with respect to parabolic metric . (Accepted October 8, 2002) Published online March 12, 2003 Communicated by L. Ambrosio     

Perturbation for fractional-order evolution equation

Fractional calculus has gained a lot of importance during the last decades, mainly because it has become a powerful tool in modeling several complex phenomena from various areas of science and engineering. This paper gives a new kind of perturbation of the order of the fractional derivative with a study of the existence and uniqueness of the perturbed fractional-order evolution equation for ^CD^a-e₀₊u(t)=A ^CD^d₀₊u(t)+f(t),^{C}D^{\alpha-\epsilon}_{0+}u(t)=A~^{C}D^{\delta}_{0+}u(t)+f(t), u(0)=u _o , α∈(0,1), and 0≤ε, δ<α under the assumption that A is the generator of a bounded C _o -semigroup. The continuation of our solution in some different cases for α, ε and δ is discussed, as well as the importance of the obtained results is specified.     

Effect of initial conditions on decaying grid turbulence at low R λ

The transport equations for the second-order velocity structure functions 〈(δu)²〉 and 〈(δq)²〉 are used as a scale-by-scale budget to quantify the effect of initial conditions at low Reynolds numbers, typical of grid turbulence. The validity of these equations is first investigated via hot-wire measurements of velocity and transverse vorticity fluctuations. The transport equation for 〈(δq)²〉 is shown to be balanced at all scales, while anisotropy of the large scales leads to a significant imbalance in the equation for 〈(δu)²〉. The effect of using similarity to evaluate the transport equation is rigorously tested. This approach has the desirable benefit of requiring less extensive measurements to calculate the inhomogeneous term of the transport equation. The similarity form of the 〈(δq)²〉 equation produces nearly identical results as those obtained without the similarity assumption. In the case of the 〈(δu)²〉 equation, the similarity method forces a balance at large separation, although the imbalance due to large scale anisotropy remains. The initial conditions of the turbulence at constant R _M ≃ 10,400 (28≤ R _λ≤ 55) are changed by using three grids of different geometries. Initial conditions affect the shape and magnitude of the second- and third-order structure functions, as well as the anisotropy of the large scales. The effect of initial conditions on the scale-by-scale budget is restricted to the inhomogeneous term of the transport equations, while the dissipation term remains unaffected despite the low R _λ. Scales as small as λ are affected by the changes in initial conditions.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Asymptotic Stability of the Wave Equation on Compact Manifolds and Locally Distributed Damping: A Sharp Result

Let (M, g) be a n-dimensional ( ${n\geqq 2}Let (M, g) be a n-dimensional ( n\geqq 2{n\geqq 2}) compact Riemannian manifold with boundary where g denotes a Riemannian metric of class C ^∞. This paper is concerned with the study of the wave equation on (M, g) with locally distributed damping, described by

l utt - Dgu+ a(x) g(ut)=0,   on M×] 0,￥[ ,u=0 on ?M ×] 0,￥[, \left. \begin{array}{l} u_{tt} - \Delta_{{\bf g}}u+ a(x)\,g(u_{t})=0,\quad\hbox{on\ \thinspace}{M}\times \left] 0,\infty\right[ ,u=0\,\hbox{on}\,\partial M \times \left] 0,\infty \right[, \end{array} \right.      16. Direct displacement method in crack theory (numerical resolution)      Joseph J. Golecki 《Meccanica》2007,42(6):555-566 Solutions in crack theory can be defined directly by opening displacement v ○=v ○(x) and u ○=u ○(x) for the first and the second mode, respectively. In this case, the boundary conditions are expressed by singular integrals of the second order. Aiming to solve numerically the problems, we apply the finite-part definition for the singular integrals and the discretization procedure. Contributed to the memory of my Mother.      17. Bifurcations of Heteroclinic Orbits      Kenneth R. Meyer  Patrick McSwiggen  Xiaojie Hou 《Journal of Dynamics and Differential Equations》2010,22(3):367-380 The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ 1 as t → −∞ and u(t) → γ 2 as t → ∞ where γ 1, γ 2 are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.      18. Finite element analysis of mode III steady crack growth in anisotropic hardening materials      Mu Sun  Shou-Wen Yu 《Acta Mechanica Solida Sinica》1988,1(3):337-344 In this paper, the steady crack growth of mode III under small scale yielding conditions is investigated for anisotropic hardening materials by the finite element method. The elastic-plastic stiffness matrix for anisotropic materials is given. The results show the significant influences of anisotropic hardening behaviour on the shape and size of plastic zone and deformation field near the crack tip. With a COD fracture criterion, the ratio of stress intensity factorsk ss/kc varies appreciably with the anisotropic hardening parameterM and the hardening exponentN.      19. The optimal dimensions of convective-radiating circular fins      Jameel-ur-Rehman Khan  S. M. Zubair 《Heat and Mass Transfer》1999,35(6):469-478 The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r o/r) n ] is studied, while variation of thermal conductivity is of the form k=k o[1+ɛ((T−T ∞)/ (T b−T ∞)) m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r o)/(r e−r o)]λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999      20. A Density Result in Two-Dimensional Linearized Elasticity,and Applications   总被引：1，自引：0，他引：1   Antonin Chambolle 《Archive for Rational Mechanics and Analysis》2003,167(3):211-233 We show that in a two-dimensional bounded open set whose complement has a finite number of connected components, the vector fields uH 1 (Ωℝ2) are dense in the space of fields whose symmetrized gradient e(u) is in L 2 (Ωℝ4). This allows us to show the continuity of some linearized elasticity problems with respect to variations of the set, with applications to shape optimization or the study of crack evolution. (Accepted September 18, 2002) Published online February 4, 2003 Commmunicated by V. Šverák

Direct displacement method in crack theory (numerical resolution)

Solutions in crack theory can be defined directly by opening displacement v _○=v _○(x) and u _○=u _○(x) for the first and the second mode, respectively. In this case, the boundary conditions are expressed by singular integrals of the second order. Aiming to solve numerically the problems, we apply the finite-part definition for the singular integrals and the discretization procedure. Contributed to the memory of my Mother.     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

Finite element analysis of mode III steady crack growth in anisotropic hardening materials

In this paper, the steady crack growth of mode III under small scale yielding conditions is investigated for anisotropic hardening materials by the finite element method. The elastic-plastic stiffness matrix for anisotropic materials is given. The results show the significant influences of anisotropic hardening behaviour on the shape and size of plastic zone and deformation field near the crack tip. With a COD fracture criterion, the ratio of stress intensity factorsk _ss/k_c varies appreciably with the anisotropic hardening parameterM and the hardening exponentN.     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

A Density Result in Two-Dimensional Linearized Elasticity,and Applications

We show that in a two-dimensional bounded open set whose complement has a finite number of connected components, the vector fields uH ¹ (Ωℝ²) are dense in the space of fields whose symmetrized gradient e(u) is in L ² (Ωℝ⁴). This allows us to show the continuity of some linearized elasticity problems with respect to variations of the set, with applications to shape optimization or the study of crack evolution. (Accepted September 18, 2002) Published online February 4, 2003 Commmunicated by V. Šverák