The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

Stress concentration in a uniformly rotating circular plate weakened by an elliptic hole and two rectilinear cuts

We consider a uniformly rotating isotropic plate whose cross-section is a two-connected domain bounded from the outside by a circle of radius R and from inside by an ellipse with semiaxes a and b. The rectilinear cuts coming to the ellipse are located symmetrically on the real axis OX. On the imaginary axis OY, two point masses m are located at the distances ±id from the center of the plate. Such a plate can be considered as a plate under the action of bulk forces (forces of inertia), X, Y, and under the action of two lumped tensile forces P = mdw ² applied at the points ±id.     

Thermodynamic driving force in ferroelectric crystals with a rank-2 laminated domain pattern, and a study of enhanced electrostriction

The thermodynamic driving force for domain growth in a rank-2 laminated ferroelectric crystal is derived in this article, and we used it, together with a homogenization theory, to study the issue of enhanced electrostrictive actuation recently reported by Burcsu et al. [2004. Large electrostrictive actuation of barium titanate single crystals. J. Mech. Phys. Solids 52, 823-846]. We derived this force from the reduction of Gibbs free energy with respect to the increase of domain concentration. It is shown that both the free energy and the thermodynamic force consist of three parts: the first arises from the difference in M₀ and M₁, the linear electromechanical compliances of the parent and product domains, respectively, at a given level of applied stress and electric field, the second stems from the electromechanical work associated with the change of spontaneous strain and spontaneous polarization during domain switch, and the third from the internal energy due to the distribution of polarizations strain and electric polarization inside the crystal. We prove that the first term is substantially lower than the second one, and the third one is identically zero with compatible domain pattern. The second one is, however, not exactly equal to the commonly written sum of the products of stress with strain, and electric field with polarization during switch, unless both domains have identical moduli in the common global axes. We also show that, with compatible domain patterns and when M₁=M₀, this driving force is identical to Eshelby's driving force acting on a flat interface due to the jump of energy-momentum tensor. Applications of the theory to a BaTiO₃ crystal subjected to a fixed axial compression and decreasing electric field from the [0 0 1] state reveal that the crystal undergoes a three-stage switching process: (i) the 0→90° switch to form a rank-1 laminate, (ii) the 0→180° switch inside the 0° domain to form a laminate I with a concurrent 90°→−90° switch inside the 90° domain to form laminate II, creating a rank-2-laminated domain pattern, and (iii) finally the 90→180° switch. It is the exchange of stability between the 0, 90°, and 180° domains under compression and electric field that is the origin of the enhanced actuation. We illustrate these intrinsic features by showing the evolution of these domains, and demonstrate how the reported large actuation strain can be attained with a rank-2 laminate.     

Plate Bending Problem for a Finite Doubly Connected Domain with a Partially Unknown Boundary

The problem of bending of an isotropic elastic plate is solved for a finite doubly connected domain whose external boundary is a convex polygon and internal boundary is a smooth closed contour. The bending problem is reduced to the analytical solution of the Riemann—Hilbert problem for a ring. The deflection of the median surface and the shape of the plate's internal boundary are found     

Theoretical and experimental studies in fracture mechanics of crack-like defects under biaxial loading

It is a common point of view in fracture mechanics that, for any geometry of the body with a crack and any boundary conditions for the loading acting in the body plane, the stress and displacement components near the crack tip can be approximated in the framework of the theory of elasticity by a one-parameter or one-term representation, i.e., strictly in terms of the stress intensity coefficients K _I and K _II for an arbitrary failure crack [1, 2]. The authors of [2] specified the Westergaard function of the singular solution for a central crack under the biaxial loading of a plate. This approximate two-component solution has satisfactory accuracy. It is clear from [2] that this method cannot be admitted as a general statement [1], although it has long been assumed to be correct. The cause is that one cannot reasonably justify neglecting the second term in the Williams representation of the stress components in the plane case in the form of eigenfunction series; the contribution of this term in the rectangular coordinate system x, y is independent of the distance from the crack tip. This method may result in a serious mistake, from both the qualitative and quantitative viewpoints, in the prediction of local stresses, displacements, and related variables that are of interest. Apparently, this can best be demonstrated by an example of biaxial loading of a plate with a crack [1]. The unfounded neglect of the second term (whose contribution is independent of the distance from the crack tip) in the series representing the stress components is the source of the above-mentioned difficulties. In this problem, the influence of the load applied in the direction parallel to the crack plane manifests itself only in the second term of the series [3]. Therefore, this term should be clearly determined and studied in detail in the case of technological welding defects (faulty fusions, incomplete fusions, undercuts, and slag inclusions) and crack-like defects (scratches and cuts) in the base metal. The influence of the stress σ _OX along the crack axis on the stress tensor σ _x , σ _y , τ _xy and on the displacements u _x and u _y is confirmed by experimental studies of cracks by the photoelasticity method [4].     

Shape optimization of two equal holes in an infinite elastic plate

Abstract

The optimal design of the stress state in elastic plate structures with openings is a problem of great significance in engineering practice. Achieving proper shape of hole can reduce stress concentration around the boundaries remarkably. The optimal shape of a single hole in an infinite plate under uniform stresses has been obtained by complex variable method based on different optimal criteria. The complex variable method is particularly suitable for the hole shape optimization in infinite plate, in which the continuous hole boundary can be represented by the mapping function. It can also be used to solve the shape optimization problems of two or more holes. However, because of the difficulty of finding the mapping function for multi connected domain, the holes are mapped onto slits or separately mapped onto a circle. In this article, the two symmetrical and identical holes are mapped onto an annulus simultaneously by the newly found mapping function, which has a general form. The maximum tangential stress around the boundaries is minimized to achieve the optimal hole shape. And the coefficients of mapping function which describe the boundary are calculated by differential-evolution algorithm.     

Planar oscillations and dissipative heating of viscoelastic plates with a piezoeffect

Planar oscillations of thin piezoplates are important within the context of using this type of piezoelements as resonator frequency filters, frequency stabilizers, elements of piezotransformers, and other technological devices. In the publications currently known one usually considers piezoplates with elastic material behavior and linear governing equations. By their mechanical nature, however, a number of piezoelements, particularly piezoceramics, are viscoelastic, which, depending on the loading conditions, can lead to substantial dissipative heating of the piezoelement and confine its operation [3]. The use of piezopolymers and their composites raises particularly important issues of dissipative heating. At the present time the behavior of a piezoelement including heating can be described by the theory of thermoelectroviscoelasticity (TEVE) [2, 3], including the interaction between electromechanical and thermal fields. The complexity of TEVE problems leads to the necessity of using numerical methods to solve them, with the finite element method (FEM) being widely used in recent years. The present study is devoted to stating and solving TEVE problems concerning thin piezoceramic plates by the FEM. We treat a thin piezoceramic plate, confined by an arbitrary contour L and polarized across its thickness. A harmonic potential difference e^it is supplied to electrodes located on the smooth boundaries of the plate. Convective heat exchange with the surrounding media of temperatures T _k ^s and T^s is implemented at the contour surfaces and boundaries free of electrodes. The heat transfer coefficients equal, respectively, _k ^T and ^T. The initial plate temperature is T₀. The smooth boundary are free of mechanical loading. The mechanical forces at the contour surfaces are distributed symmetrically with respect to the mean plane of the plate.S. P. Timoshenko Institute of Mechanics, Ukrainian Academy of Sciences, Kiev (Ukraine). Translated from Prikladnaya Mekhanika, Vol. 30, No. 2, pp. 69–76, February 1994.     

Thermal stresses induced by a point heat source in a circular plate by quasi-static approach

The present paper deals with the determination of quasi-static thermal stresses due to an instantaneous point heat source of strength g_pi situated at certain circle along the radial direction of the circular plate and releasing its heat spontaneously at time t = τ. A circular plate is considered having arbitrary initial temperature and subjected to time dependent heat flux at the fixed circular boundary of r = b. The governing heat conduction equation is solved by using the integral transform method, and results are obtained in series form in terms of Bessel functions. The mathematical model has been constructed for copper material and the thermal stresses are discussed graphically.     

Exact solution of problem of deflection of an anisotropic plate with an aperture

The exact solution of the problem of the deflection of an anisotropic plate weakened by an aperture is known only for the case in which the aperture has the shape of a circle or an ellipse [1, 2]. An exact solution has not been derived for any other aperture shapes. Approximate methods [3–6] which are widespread for the case of multiply connected anisotropic plates [7] are applied to the determination of the bending moments in an anisotropic plate near an aperture differing little from an elliptical or circular one.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 168–177, September–October, 1977.     

Robust stabilization the Korteweg–de Vries–Burgers equation by boundary control

The problem of robust global stabilization by nonlinear boundary feedback control for the Korteweg–de Vries–Burgers equation on the domain [0,1] is considered. The main purpose of this paper is to derive nonlinear robust boundary control laws which make the system robustly globally asymptotically stable in spite of uncertainty in the system parameters. Furthermore, we show that the proposed boundary controllers guarantee L ₂-robust exponential stability, L _∞-robust asymptotic stability and boundedness in terms of both L ₂ and L _∞.     

Extended meshfree analysis of transverse and inplane loading of a laminated anisotropic plate of general planform geometry

An extended meshfree method is presented for the analysis of a laminated anisotropic plate under elastostatic loading. The plate may be of any planform shape with its thickness profile composed of perfectly bonded uniform thickness layers of distinct anisotropic materials. Both transverse and inplane loads are considered using a first order shear deformation theory for flexural behavior and generalized plane stress for the membrane behavior. In this extended meshfree method, a rectangular domain is initially considered with the plate of arbitrary geometry inscribed within it. A particular solution in the form of an analytic generalized Navier solution (a compound double Fourier series) is used to capture the response due to the loading within the rectangular domain. Then, a homogeneous solution by meshfree analysis is added to treat the augmented boundary conditions on the actual contour of the plate. These augmented conditions are composed of the prescribed values and that of the particular solution evaluated around the plate’s contour.Concentrated transverse and inplane loads in the form of uniform loads over a very small patch are considered with this generalized Navier solution representation. When a meshfree portion is added to account for the boundary conditions, such solutions constitute the Green’s functions for the plate. The viability of these double Fourier series representations is shown by the convergence rates for the kinematic and force/moment fields. An additional example of a two layer ±30° angleply circular plate is given to illustrate the capability of this extended meshfree method.     

Viscous flow in pipes whose cross-sections are doubly connected regions

Slow and steady viscous fluid flow in a pipe whose cross-section is bounded by an ellipse and a circle is considered. The method used is in deriving general solutions independently on the two boundaries using one or two conformal mapping functions and an attempt is made to make the complex potentials continuous in the doubly connected region. Numerical results are found and compared with those for the case of concentric circles and the case of concentric ellipses.     

Combined forced and free convection MHD row past a semi-infinite vertical plate

Effects of a transversely applied magnetic field on the forced and free convective flow of an electrically conducting fluid past a vertical semi-infinite plate, on taking into account dissipative heat and stress work, have been presented. Without magnetic field, it has been discussed by the authors [1] in an earlier paper. The effects of Gr (Grashof number, Gr>0 cooling of the plate by free convection currents, Gr<0 heating of the plate by free convection currents), Pr (Prandtl number), F (Froude number) and M² (the magnetic field parameter) are discussed. It is observed that reverse type of flow of air exists near the plate when Gr<0.

---

Gemischte erzwungene und freie Konvektions-MHD-Strömung an einer halbunendlichen senkrechten Platte

Zusammenfassung Die Wirkung eines transversalen Magnetfeldes auf die erzwungene und freie Konvektion einer elektrisch leitenden Flüssigkeit an einer halbunendlichen senkrechten Platte wurde unter Berücksichtigung der Dissipationswärme und der Kompressionsarbeit mitgeteilt. Das Problem wurde ohne Magnetfeld schon früher [1] behandelt. Diskutiert wurde die Wirkung der Grashof-Zahl Gr (Gr>0 Kühlung der Platte durch freie Konvektion), der Prandtl-Zahl Pr, der Froude-Zahl F und des magnetischen Feldparameters M₂. Bei Gr<0 wird Umkehr strömung in Plattennähe betrachtet.

     

Flow in the neighborhood of the trailing edge of a plate in a transonic flow of viscous gas

An investigation is made into the influence of the Mach number and the viscosity on the flow in the neighborhood of the trailing edge of a plate. The Mach number is assumed to satisfy m_∞ ² ? 1 = 0(R^?l/5), which corresponds to the regime of transonic interaction. It is shown that if the Mach number is such that ¦M_∞ ² ? 1¦ > O(R^?1/5) the problem in the region of free interaction can be reduced by an appropriate transformation to the already known solutions for an incompressible fluid [5] and supersonic flow [7].     

Interaction between an Inclined Shock and Boundary and High-Entropy Layers on a Flat Plate

The flow structure and heat exchange in the zone of interference between an inclined shock and the surface of a flat plate are investigated experimentally and theoretically as functions of many parameters, the interference being studied in both the presence and the absence of bluntness of the leading edge. The experiments were carried out at Mach numbers M = 6, 8, and 10 and the Reynolds numbers Re _L , calculated using the plate length L = 120 mm and the free-stream parameters, varied over the range from 0.24 ? 10⁶ to 1.31 ? 10⁶. The bluntness radius of the leading edge of the plate, the intensity of the impinging shock, and its location with respect to the leading edge were varied. The numerical simulation was carried out by solving the complete two-dimensional Navier-Stokes equations and averaged Reynolds equations using the q-ω turbulence model. The laminar boundary layer became turbulent inside the separation zone induced by the shock. It is shown that the plate bluntness significantly reduces the heat exchange intensity in the interference zone, this effect intensifying with increase in the Mach number.     

Dynamics as a mechanism preventing the formation of finer and finer microstructure

We study the dynamics of pattern formation in the one-dimensional partial differential equation $$u_u - (W'(u_x ))_x - u_{xxt} + u = 0{\text{ (}}u = u(x,t),{\text{ }}x \in (0,1),{\text{ }}t > 0)$$ proposed recently by Ball, Holmes, James, Pego & Swart [BHJPS] as a mathematical “cartoon” for the dynamic formation of microstructures observed in various crystalline solids. Here W is a double-well potential like 1/4((u _x )² ?1)². What makes this equation interesting and unusual is that it possesses as a Lyapunov function a free energy (consisting of kinetic energy plus a nonconvex “elastic” energy, but no interfacial energy contribution) which does not attain a minimum but favours the formation of finer and finer phase mixtures: $$E[u,u_t ] = \int\limits_0^1 {(\frac{{u_t^2 }}{2} + W(u_x ) + \frac{{u^2 }}{2})dx.}$$ Our analysis of the dynamics confirms the following surprising and striking difference between statics and dynamics, conjectured in [BHJPS] on the basis of numerical simulations of Swart & Holmes [SH]:

?While minimizing the above energy predicts infinitely fine patterns (mathematically: weak but not strong convergence of all minimizing sequences (u _nv_n) of E[u,v] in the Sobolev space W ^{1 p}(0, 1)×L²(0,1)), solutions to the evolution equation of ball et al. typically develop patterns of small but finite length scale (mathematically: strong convergence in W ₁ ^p(0,1)×L²(0,1) of all solutions (u(t),u_t(t)) with low initial energy as time t → ∞).

Moreover, in order to understand the finer details of why the dynamics fails to mimic the behaviour of minimizing sequences and how solutions select their limiting pattern, we present a detailed analysis of the evolution of a restricted class of initial data — those where the strain field u _x has a transition layer structure; our analysis includes proofs that

?at low energy, the number of phases is in fact exactly preserved, that is, there is no nucleation or coarsening

?transition layers lock in and steepen exponentially fast, converging to discontinuous stationary sharp interfaces as time t → ∞

?the limiting patterns — while not minimizing energy globally — are ‘relative minimizers’ in the weak sense of the calculus of variations, that is, minimizers among all patterns which share the same strain interface positions.

     

A Remark on the Existence of 2-D Steady Navier–Stokes Flow in Bounded Symmetric Domain Under General Outflow Condition

Let Ω be a 2-dimensional bounded domain, symmetric with respect to the x₂-axis. The boundary has several connected components, intersecting the x₂-axis. The boundary value is symmetric with respect to the x₂-axis satisfying the general outflow condition. The existence of the symmetric solution to the steady Navier–Stokes equations was established by Amick [2] and Fujita [4]. Fujita [4] proved a key lemma concerning the solenoidal extension of the boundary value by virtual drain method. In this note, we give a different proof via elementary approach by means of the stream function.     

The Two-Dimensional Euler Equations on Singular Domains

We establish the existence of global weak solutions of the two-dimensional incompressible Euler equations for a large class of non-smooth open sets. Loosely, these open sets are the complements (in a simply connected domain) of a finite number of obstacles with positive Sobolev capacity. Existence of weak solutions with L ^p vorticity is deduced from a property of domain continuity for the Euler equations that relates to the so-called γ-convergence of open sets. Our results complete those obtained for convex domains in Taylor (Progress in Nonlinear Differential Equations and their Applications, Vol. 42, 2000), or for domains with asymptotically small holes (Iftimie et al. in Commun Partial Differ Equ 28(1–2), 349–379, 2003; Lopes Filho in SIAM J Math Anal 39(2), 422–436, 2007).     

Interplay between polarization rotation and crack propagation in PMN-PT relaxor single crystals

Investigations on the interconnection between the polarization rotation and crack propagation are performed for [110]-oriented 74Pb(Mg_1/3Nb_2/3)O₃-26PbTiO₃ relaxor ferroelectric single crystal under electric loadings along [001] direction. The crystal is of predominantly monoclinic M_A phase with scatter distributed rhombohedral (R) phase under a moderate poling field of 900 V/mm in [001] direction. With magnitude of 800 V/mm, a through thickness crack is initiated near the electrode by electric cycling. Static electric loadings is then imposed to the single crystal. As the applied static electric field increases, domain switching in the monoclinic M_A phase and phase transition from M_A to R phase occur near the crack. The results indicate that the crack features a conducting one. Whether domain switching or phase transition occurs depends on the intensity of the electric field component that is perpendicular to the applied electric field.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.