Some particular solutions for penny-shaped crack problem by using hypersingular integral equation or differential-integral equation

Summary A hypersingular integral equation or a differential-integral equation is used to solve the penny-shaped crack problem. It is found that, if a displacement jump (crack opening displacement COD) takes the form of (a ²−x ²−y ²)^1/2 x ^m y ⁿ , where a denotes the radius of the circular region, the relevant traction applied on the crack face can be evaluated in a closed form, and the stress intensity factor can be derived immediately. Finally, some particular solutions of the penny-shaped crack problem are presented in this paper. Received 1 July 1997; accepted for publication 13 October 1997     

Strain fields caused by doughnut-like and tubular inclusions with uniform eigenstrains

Inclusion problems for elongated toroidal inclusions are solved to discuss the strain fields caused by eigenstrained doughnut-like and tubular inclusions. For an infinitely long tubular inclusion, uniform eigenstrains do not generate elastic strains in the matrix region inside the tube. Effects of tube length on the elastic strains are shown for the purely dilatational eigenstrain ^*. For the tubular inclusion with the length 2b, the elastic strain at the center of the matrix region inside the tube becomes the maximum as much as ≈(a/b)^* when the tubular inclusion has the inner diameter ≈2(b−a) and the outer diameter ≈2(b+a).     

Applying the generalized step-function to solve the bending problems of nonuniform beams with variable cross section

The bending problems of nonuniform beams with variable cross-section can be approximated by that of a step beam under sectionally uniform load (including both concentrated forces and couples). In this paper, the concept of Heaviside function {x-a}⁰ will be generalized, and a new function {x-a}⁰, n=0,1,2…,will be defined, which may be named as a generalized step function. The rules of operation will also be given to {x-a}ⁿ{x-b}⁰. The reciprocal of the flexural of rigidity 1/EJ and the bending moment M(x) can all be expressed in terms of {x-a}ⁿ,and substituted into the differential equation of the elastic curve of the beam respectively. Thus we may establish a set of unified method to solve various types of bending problems of straight beams. The general solution of the deflection equation will be given.     

On Asymptotics of Solutions of Nonlinear Differential Equations of the Second Order

We study the problem of asymptotics of unbounded solutions of differential equations of the form y″ = α₀ p(t)ϕ(y), where α₀ ∈ {−1, 1}, p: [a, ω[→]0, +∞[, −∞ < a < ω ≤ +∞, is a continuous function, and ϕ: [y ₀, +∞[→]0, +∞[ is a twice continuously differentiable function close to a power function in a certain sense.__________Translated from Neliniini Kolyvannya, Vol. 8, No. 1, pp. 18–28, January–March, 2005.     

Plasticity of a high strength steel alloy

Experimentally determined plastic constitutive equations of the parabolic form (σ − σ_y = β(ε − ε_y)1/2) are presented for a high strength alloy steel. Two deformation moduli β were required to describe the quasistatic compression data, both of which, as well as the point of change, were predicted by a mode index and transition strain structure of a general theory of plasticity. Dynamic strain, duration of impact, and final strain distribution were measured on specimen rods subjected to axial, symmetric, and constant velocity impact. The dynamic yield stress was 16% higher than the quasistatic value. The dynamic response function, deduced from a simple wave propagation theory, was also parabolic and a single deformation modulus, equal to the initial quasistatic value, applied. Thus, it was shown that the form of the quasistatic response function was preserved in dynamic loading, and that the increase of the dynamic stress was due to the increase of the yield stress.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

Darboux problem for a differential equation with fractional derivative

We find conditions for the unique solvability of the problem u _xy (x, y) = f(x, y, u(x, y), (D ₀ ^r u)(x, y)), u(x, 0) = u(0, y) = 0, x ∈ [0, a], y ∈ [0, b], where (D ₀ ^r u)(x, y) is the mixed Riemann-Liouville derivative of order r = (r ₁, r ₂), 0 < r ₁, r ₂ < 1, in the class of functions that have the continuous derivatives u _xy (x, y) and (D ₀ ^r u)(x, y). We propose a numerical method for solving this problem and prove the convergence of the method. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 4, pp. 456–467, October–December, 2005.     

Effects of deformed volume, volume fraction and particle size on the deformation behaviour of W/Cu composites

Tungsten/copper (W/Cu) particle reinforced composites were used to investigate the scaling effects on the deformation and fracture behaviour. The effects of the volume fraction and the particle size of the reinforcement (tungsten particles) were studied. W/Cu-80/20, 70/30 and 60/40 wt.% each with tungsten particle size of 10 μm and 30 μm were tested under compression and shear loading. Cylindrical compression specimens with different volumes (D_S = H) were investigated with strain rates between 0.001 s⁻¹ and about 5750 s⁻¹ at temperatures from 20 °C to 800 °C. Axis-symmetric hat-shaped shear specimens with different shear zone widths were examined at different strain rates as well. A clear dependence of the flow stress on the deformed volume and the particle size was found under compression and shear loading. Metallographic investigation was carried out to show a relation between the deformation of the tungsten particles and the global deformation of the specimens. The size of the deformed zone under either compression or shear loading has shown a clear size effect on the fracture of the hat-shaped specimens.The quasi-static flow curves were described with the material law from Swift. The parameters of the material law were presented as a function of the temperature and the specimen size. The mechanical behaviour of the composite materials were numerically computed for an idealized axis-symmetric hat-shaped specimen to verify the determined material law.     

Self-similarity of nonstationary boundary-layer motions

In this paper we show that problems concerning the development of a boundary layer on a semi-infinite plate when the outer flow speed is of the form U = (1 + ct)^b a, and on a cylinder when the outer flow speed has the forms U = ctx^m and U = (1 + ct)^b ax^m, are self-similar. We present the results of numerical calculations for various values of , b, and m. We consider the problem of a stepwise nonstationary heating of a plate, impulsively set into motion in an incompressible fluid; we show that this problem is self-similar and obtain its solution numerically.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 122–125, July–August, 1975.     

Erratum to: Floquet’s Theorem and Stability of Periodic Solitary Waves

This paper is concerned with the spectrum the Hill operator L(y) = −y′′ + Q(x) y in L²_per[0, p]{L^2_{{\rm per}}[0, \pi]}. We show that the eigenvalues of L can be characterized by knowing one of its eigenfunctions. Applications are given to nonlinear stability of a class of periodic problems.     

SINGULAR PERTURBATIONS FOR A CLASS OF BOUNDARY VALUE PROBLEMS OF HIGHER ORDER NONLINEAR DIFFERENTIAL EQUATIONS

ＳＩＮＧＵＬＡＲ　ＰＥＲＴＵＲＢＡＴＩＯＮＳ　ＦＯＲ　Ａ　ＣＬＡＳＳ　ＯＦ　ＢＯＵＮＤＡＲＹ　ＶＡＬＵＥ　ＰＲＯＢＬＥＭＳ　ＯＦ　ＨＩＧＨＥＲ　ＯＲＤＥＲ　ＮＯＮＬＩＮＥＡＲ　ＤＩＦＦＥＲＥＮＴＩＡＬ　ＥＱＵＡＴＩＯＮＳＳｈｉＹｕａｎｉｎｇ（史玉明）...     

On nonlinear problems of mixed type: A qualitative theory using infinite-dimensional center manifolds

We consider the equation a(y)u_xx+div_y(b(y)_yu)+c(y)u=g(y, u) in the cylinder (–l,l)×, being elliptic where b(y)>0 and hyperbolic where b(y)<0. We construct self-adjoint realizations in L₂() of the operatorAu= (1/a) div_y(b_yu)+(c/a) in the case ofb changing sign. This leads to the abstract problem u_xx+Au=g(u), whereA has a spectrum extending to + as well as to –. For l= it is shown that all sufficiently small solutions lie on an infinite-dimensional center manifold and behave like those of a hyperbolic problem. Anx-independent cross-sectional integral E=E(u, u_x) is derived showing that all solutions on the center manifold remain bounded forx ±. For finitel, all small solutionsu are close to a solution on the center manifold such that u(x)-(x) Ce ^-(1-|x|) for allx, whereC and are independent ofu. Hence, the solutions are dominated by hyperbolic properties, except close to the terminal ends {±1}×, where boundary layers of elliptic type appear.     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

Floquet’s Theorem and Stability of Periodic Solitary Waves

This paper is concerned with the spectrum the Hill operator L(y) = −y′′ + Q(x) y in L²_per[0, p]{L^{2}_{\rm per}[0, \pi]} . We show that the eigenvalues of L can be characterized by knowing one of its eigenfunctions. Applications are given to nonlinear stability of a class of periodic problems.     

Non-linear oscillation of circular cylindrical shells

The method of multiple scales is used to analyze the non-linear forced response of circular cylindrical shells in the presence of a two-to-one internal (autoparametric) resonance to a harmonic excitation having the frequency Ω. If ω_r and a_r denote the frequency and amplitude of a flexural mode and ω_b and a_b denote the frequency and amplitude of the breathing mode, the steady-state response exhibits a saturation phenomenon when ω_b ≈ 2ω_r, if the excitation frequency Ω is near ω_b. As the amplitude ƒ of the excitation increases from zero, a_b increases linearly whereas a_r remains zero until a threshold is reached. This threshold is a function of the damping coefficients and ω_b−2ω_r. Beyond this threshold a_b remains constant (i.e. the breathing mode saturates) and the extra energy spills over into the flexural mode. In other words, although the breathing mode is directly excited by the load, it absorbs a small amount of the input energy (responds with a small amplitude) and passes the rest of the input energy into the flexural mode (responds with a large amplitude). For small damping coefficients and depending on the detunings of the internal resonance and the excitation, the response exhibits a Hopf bifurcation and consequently there are no steadystate periodic responses. Instead, the responses are amplitude- and phase-modulated motions. When Ω ≈ ω_r, there is no saturation phenomenon and at close to perfect resonance, the response exhibits a Hopf bifurcation, leading again to amplitude- and phase-modulated or chaotic motions.     

General high-order rational solutions and their dynamics in the (3+1)-dimensional Jimbo–Miwa equation

General high-order rational solutions are derived for the (3+1)-dimensional Jimbo–Miwa equation based on the Hirota bilinear form. The solutions are presented in terms of Gram determinants; the elements of determinants are connected to Schur polynomials and have simple algebraic expressions. Their dynamic behaviors are researched using three-dimensional imagery and contour plots. It is revealed that different kinds of solutions appear in (x, y) plane and (y, z) plane. When one of these internal parameters in the rational solutions is sufficiently large, in (x, y) plane Lump solutions appear with obvious geometric structures, which are deconstructed by a first-order Lump such as triangle, pentagon, and nonagon, among others; in (y, z) plane rational line soliton solutions with maximum background amplitude changing over time appear. These findings might help us comprehend the nonlinear wave propagation processes in the many nonlinear physical systems.

     

Thermal stresses in rectangular plates

Summary A method of determining the thermal stresses in a flat rectangular isotropic plate of constant thickness with arbitrary temperature distribution in the plane of the plate and with no variation in temperature through the thickness is presented. The thermal stress have been obtained in terms of Fourier series and integrals that satisfy the differential equation and the boundary conditions. Several examples have been presented to show the application of the method.Nomenclature x, y rectangular coordinates - _x, _y direct stresses - _xy shear stress - ø Airy's stress function - E Young's modulus of elasticity - coefficient of thermal expansion - T temperature - ² Laplace operator: - ⁴ biharmonic operator - 2a length of the plate - 2b width of the plate - a/b aspect ratio - a _mr, b_ms, c_nr, d_ns Fourier coefficients defined in equation (6) - _m=m/a m=1, 2, 3, ... _n=n/2a n=1, 3, 5, ... - _r=r/b r=1, 2, 3, ... _s=s/2b s=1, 3, 5, ... - A _m, B_m, C_n, D_n, E_r, F_r, G_s, H_s Fourier coefficients - K _rand L _s Fourier coefficients defined in equation (20) - direct stress at infinity - T ₁(x, y) temperature distribution symmetrical in x and y - T ₂(x, y) temperature distribution symmetrical in x and antisymmetrical in y - T ₃(x, y) temperature distribution antisymmetrical in x and symmetrical in y - T ₄(x, y) temperature distribution antisymmetrical in x and y     

Heat transfer and pressure drop for purely viscous non-Newtonian fluids in turbulent flow through rectangular passages

Experimental measurements of friction factor and heat transfer for the turbulent flow of purely viscous non-Newtonian fluids in a 21 rectangular channel are compared with results previously reported for the circular tube geometry. Comparisons are also made with available analytical and empirical predictions.It is found that the rectangular duct fully established friction factor measurements are within ± 5% of the Dodge-Metzner prediction if the Kozicki generalized Reynolds number is used. A modified form of the simpler explicit equation proposed by Yoo, [i.e.f=0.079n ^0.675(Re ^*)^–0.25], is found to yield predictions for both the rectangular duct and the circular tube geometries with approximately the same accuracy as the Dodge-Metzner equation.Fully developed Stanton numbers for the rectangular duct are in good agreement with the circular tube data over a range ofn from 0.37 to 0.88 for a given Prandtl number,Pr _a , when compared at a fixed value of the Reynolds number based on the apparent viscosity evaluated at the wall shear stress. In general, the experimental data are within ± 20% of Yoo's equation,St=0.0152Re _a ^–0.155 Pr _a ^–2/3 . A new equation is proposed to bring the prediction for circular pipes as well as rectangular channels into better agreement with generally accepted Newtonian heat transfer results.

---

Wärmeübergang und Druckverlust für viskose nicht-Newtonsche Fluide in turbulenter Strömung durch rechteckige Kanäle

Zusammenfassung Es werden Messungen des Reibungsfaktors und des Wärmeübergangs bei turbulenter Strömung viskoser nicht-Newtonscher Fluide in einem rechteckigen Kanal mit dem Seitenverhältnis 21 verglichen mit früheren Ergebnissen, die an runden Rohren gewonnen wurden. Weiterhin werden Vergleiche mit aus der Literatur verfügbaren analytischen und empirischen Beziehungen gemacht.Es zeigte sich, daß die Messungen des Reibungsfaktors im rechteckigen Kanal bei vollausgebildeter Strömung auf ± 5% mit der Vorhersage von Dodge-Metzner übereinstimmen, wenn die von Kozicki verallgemeinerte Reynolds-Zahl verwendet wird. Eine modifizierte Form der einfachen von Yoo vorgeschlagenen einfachen Gleichung in explizierter Form (f=0,079n ^0,675(Re ^*)^–0,25) bewies, daß sie sowohl für den rechteckigen Kanal als auch das runde Rohr die Werte mit fast der gleichen Genauigkeit wie die Methode von Dodge-Metzner vorhersagen kann.Die Stanton-Zahlen für den rechteckigen Kanal bei vollausgebildeter Strömung sind in guter Übereinstimmung mit den Werten für das runde Rohr in einem Bereich vonn= 0,37 – 0,88 für eine gegebene Prandtl-Zahl, wenn man den Vergleich bei einem vorgegebenen Wert der Reynolds-Zahl anstellt, die auf die scheinbare Viskosität — abgeleitet aus der Wandschubspannungbezogen ist. Generell läßt sich sagen, daß die Werte auf ± 20% mit der Gleichung von Yoo (St=0,0152Re _a ^–0,155 )Pr _a ^–2/3 ) übereinstimmen. Es wird eine neue Gleichung vorgeschlagen, welche sowohl die Werte für runde Rohre als auch die für rechteckige Kanäle in bessere Übereinstimmung bringt mit den in der Literatur üblichen Ergebnissen für den Wärmeübergang an Newtonsche Fluide.

---

Nomenclature a constant in Eq. (8) - A area of cross-section of channel [m²] - b constant in Eq. (8) - c _p specific heat of test fluid [J kg^–1 K^–1] - d capillary tube diameter [m] - D _h hydraulic diameter, 4A/P[m] - f Fanning friction factor, _w/(g9 V²/2) - h axially local (spanwise averaged) heat transfer coefficient,q _w /(T_wi-T_b) [Wm^–2 K^–1] - k _f thermal conductivity of test fluid [Wm^–1K^–1] - K consistency index of power law fluid - n power law index - Nu fully established, local (spanwise averaged) Nusselt numberh D _h /k _f - P perimeter of channel [m] - Pr _a Prandtl number based on apparent viscosjity, c _p /k _f - Pr ^* defined as (Re _a Pr _a )/Re ^* - q _w wall heat flux [Wm^–2] - Re _a Reynolds number based on apparent viscosity, VD _h/ - Re Metzner's generalized Reynolds number in Eq. (2) - Re ^* Reynolds number defined in Eq. (8) - St Stanton number,h/( V c_p) - T _b local bulk temperature of the fluid [K] - T _wi local inside wall temperature [K] - T _wo local outside wall temperature [K] - V bulk flow velocity [m s^–1] - x distance from the inlet of channel along flow direction [m] Greek symbols shear rate [s^–1] - apparent viscosity [Pa s] - density of test fluid [kg m^–3] - shear stress [Pa] - _w shear stress at the wall [Pa] Dedicated to Prof. Dr.-Ing. U. Grigull's 75th birthday