The state vector methods for space axisymmetric problems in multilayered piezoelectric media

The state vector equations for space axisymmetric problems of transversely isotropic piezoelectric media are established from the basic equations. Using the Hankel transform, the state vector equations are reduced to a system of ordinary differential equations. An analytical solution of the problems in the Hankel transform space is presented in the form of the product of initial state vector and transfer matrix. The transfer matrices are given for the three distinct eigenvalues. Applications of the solutions are discussed. An analytical solution for the transversely isotropic semi-infinite piezoelectric media subjected to concerted point loads on the surface z=0 is presented in the Hankel transform space. Using transfer matrix and the continuity conditions at the layer interfaces, the general solution formulation of N-layered transversely isotropic piezoelectric media is given. A selected set of numerical solutions is presented for a layered semi-infinite piezoelectric solid.     

Love solutions in the linear inhomogeneous transversely isotropic theory of elasticity

A general Love solution for the inhomogeneous transversely isotropic theory of elasticity with the elastic constants dependent on the coordinate z is proposed. This result may be considered as a generalization of the Love solutions we recently derived for the inhomogeneous isotropic theory of elasticity. The key steps of deriving the Love solution for the classical linear homogeneous transversely isotropic theory of elasticity are described for further use of the derivation procedure, which is then generalized to the inhomogeneous transversely isotropic case. Some particular cases of inhomogeneity traditionally used in the theory of elasticity are also examined. The significance of the derived solutions and their importance for the modeling of functionally graded materials are briefly discussed     

Spherically isotropic,elastic spheres subject to diametral point load strength test

This paper presents an analytic solution for the stress concentrations within a spherically isotropic, elastic sphere of radius R subject to diametral point load strength test. The method of solution uses the displacement potential approach together with the Fourier–Legendre expansion for the boundary loads. For the case of isotropic sphere, our solution reduces to the solution by Hiramatsu and Oka, 1966 and agrees well with the published experimental observations by Frocht and Guernsey (1953) . A zone of higher tensile stress concentration is developed near the point loads, and the difference between this maximum tensile stress and the uniform tensile stress in the central part of the sphere increases with E/E′ (where E and E′ are the Youngs moduli governing axial deformations along directions parallel and normal to the planes of isotropy, respectively) , G′/G (where G and G′ are the moduli governing shear deformations in the planes of isotropy and the planes parallel to the radial direction) , and ν̄/ν′ (where ν̄ and ν′ are the Poissons ratios characterizing transverse reduction in the planes of isotropy under tension in the same plane and under radial tension, respectively) . This stress difference, in general, decreases with the size of loading area and the Poissons ratio.     

Elasticity solutions for a transversely isotropic functionally graded circular plate subject to an axisymmetric transverse load qrk

This paper considers the bending of transversely isotropic circular plates with elastic compliance coefficients being arbitrary functions of the thickness coordinate, subject to a transverse load in the form of qr^k (k is zero or a finite even number). The differential equations satisfied by stress functions for the particular problem are derived. An elaborate analysis procedure is then presented to derive these stress functions, from which the analytical expressions for the axial force, bending moment and displacements are obtained through integration. The method is then applied to the problem of transversely isotropic functionally graded circular plate subject to a uniform load, illustrating the procedure to determine the integral constants from the boundary conditions. Analytical elasticity solutions are presented for simply-supported and clamped plates, and, when degenerated, they coincide with the available solutions for an isotropic homogenous plate. Two numerical examples are finally presented to show the effect of material inhomogeneity on the elastic field in FGM plates.     

Existence and Stability of Solutions for Steady Flows of Fibre Suspension Flows

We establish existence, uniqueness, convergence and stability of solutions to the equations of steady flows of fibre suspension flows. The existence of a unique steady solution is proven by using an iterative scheme. One of the restrictions imposed on the data confirms a well known fact proven in Galdi and Reddy (J Non-Newtonian Fluid Mech 83:205–230, 1999), Munganga and Reddy (Math Models Methods Appl Sci 12:1177–1203, 2002) and Munganga et al. (J Non-Newtonian fluid Mech 92:135–150, 2000) that the particle number N _p must be less than 35/2. Exact solutions are calculated for Couette and Poiseuille flows. Solutions of Poiseuille flows are shown to be more accurate than those of Couette flow when a time perturbation is considered.     

Scaling relations for elongational flow of polystyrene melts and concentrated solutions of polystyrene in oligomeric styrene

A consistent model of the rheology of polymer melts and concentrated solutions is presented, based on the idea that the pressures exerted by a polymer chain on the walls of an anisotropic confinement are anisotropic (Doi and Edwards. The Theory of Polymer Dynamics, Oxford University Press, 1986). In a tube model with variable tube diameter, chain stretch and tube diameter reduction are related, and at deformation rates larger than the inverse Rouse time τ _R, the chain is stretched and its confining tube becomes increasingly anisotropic. Tube diameter reduction leads to an interchain pressure in the lateral direction of the tube (Marrucci and Ianniruberto. Macromolecules 37:3934-3942, 2004). Chain stretch is balanced by interchain tube pressure in the lateral direction, which is proportional to the third power of stretch, and by a spring force in the longitudinal direction of the tube, which is linear in stretch. Analyzing elongational viscosity data of Huang et al. (Macromolecules 46:5026-5035, 2013a; ACS Macro Letters 2:741-744, 2013b) shows that dilution of polystyrene by oligomeric styrene does not change the relative interchain tube pressure. Based on this extended interchain pressure concept, scaling relations for linear viscoelasticity and elongational viscosity of polystyrene melts and concentrated solutions of polystyrene in oligomeric styrene are presented based exclusively on the relaxation modulus of a reference polymer melt, the volume fraction of polymer in the solution, and the time-molar-mass shift as well as the time-temperature shift caused by the reduction of the glass transition temperature T _g of the polymer in a solution relative to T _g of the melt.     

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).     

On the Energy Dissipation Rate of Solutions to the Compressible Isentropic Euler System

In this paper we extend and complement the results in Chiodaroli et al. (Global ill-posedness of the isentropic system of gas dynamics, 2014) on the well-posedness issue for weak solutions of the compressible isentropic Euler system in 2 space dimensions with pressure law p(ρ) = ρ ^γ , γ ≥ 1. First we show that every Riemann problem whose one-dimensional self-similar solution consists of two shocks admits also infinitely many two-dimensional admissible bounded weak solutions (not containing vacuum) generated by the method of De Lellis and Székelyhidi (Ann Math 170:1417–1436, 2009), (Arch Ration Mech Anal 195:225–260, 2010). Moreover we prove that for some of these Riemann problems and for 1 ≤ γ < 3 such solutions have a greater energy dissipation rate than the self-similar solution emanating from the same Riemann data. We therefore show that the maximal dissipation criterion proposed by Dafermos in (J Diff Equ 14:202–212, 1973) does not favour the classical self-similar solutions.     

Limit equilibrium of a transversely isotropic spherical shell with two surface cracks

An analog of the δ _c -model is used to reduce the limit-equilibrium problem for a transversely isotropic spherical shell with surface cracks to a system of integral equations. An algorithm for numerical solution of this system is proposed     

Singular Limiting Induced from Continuum Solutions and the Problem of Dynamic Cavitation

In the works of Pericak-Spector and Spector (Arch Rational Mech Anal. 101:293–317, 1988, Proc. Royal Soc. Edinburgh Sect A 127:837–857, 1997) a class of self-similar solutions are constructed for the equations of radial isotropic elastodynamics that describe cavitating solutions. Cavitating solutions decrease the total mechanical energy and provide a striking example of non-uniqueness of entropy weak solutions (for polyconvex energies) due to point-singularities at the cavity. To resolve this paradox, we introduce the concept of singular limiting induced from continuum solution (or slic-solution), according to which a discontinuous motion is a slic-solution if its averages form a family of smooth approximate solutions to the problem. It turns out that there is an energetic cost for creating the cavity, which is captured by the notion of slic-solution but neglected by the usual entropic weak solutions. Once this cost is accounted for, the total mechanical energy of the cavitating solution is in fact larger than that of the homogeneously deformed state. We also apply the notion of slic-solutions to a one-dimensional example describing the onset of fracture, and to gas dynamics in Langrangean coordinates with Riemann data inducing vacuum in the wave fan.     

Determination of boundary layers in the plane problem for three-layer strips. Part 2

In the first part of this paper, we considered the exact statement of the plane elasticity problem in displacements for strips made of various materials (problem A, an isotropic material; problem B, an orthotropic material with 2G ₁₂ < √E ₁ E ₂; problem C, an orthotropic material with 2G ₁₂ > √E ₁ E ₂). Further, we stated and solved the boundary layer problem (the problem on a solution decaying away from the boundary) for a sandwich strip of regular structure consisting of isotropic layers (problem AA). In the present paper, we use the solution of the plane problem to consider the problem for sandwich strips of regular structure with isotropic face layers and orthotropic filler (problem AB).     

Several three-dimensional solutions for transversely isotropic functionally graded material plate welded with circular inclusion

The problem of a transversely isotropic functionally graded material (FGM) plate welded with a circular inclusion is considered. The analysis starts with the generalized England-Spencer plate theory for transversely isotropic FGM plates, which expresses a three-dimensional (3D) general solution in terms of four analytic functions. Several analytical solutions are then obtained for an infinite FGM plate welded with a circular inclusion and subjected to the loads at infinity. Three different cases are considered, i.e., a rigid circular inclusion fixed in the space, a rigid circular inclusion rotating about the x-, y-, and z-axes, and an elastic circular inclusion with different material constants from the plate itself. The static responses of the plate and/or the inclusion are investigated through numerical examples.     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

Flexural-gravity wave resistances due to a surface-moving line source on a fluid covered by a thin elastic plate

Analytical solutions for the flexural-gravity wave resistances due to a line source steadily moving on the surface of an infinitely deep fluid are investigated within the framework of the linear potential theory. The homogenous fluid, covered by a thin elastic plate, is assumed to be incompressible and inviscid, and the motion to be irrotational. The solution in integral form for the wave resistance is obtained by means of the Fourier transform and the explicitly analytical solutions are derived with the aid of the residue theorem. The dispersion relation shows that there is a minimal phase speed c_min, a threshold for the existence of the wave resistance. No wave is generated when the moving speed of the source V is less than c_min while the wave resistances firstly increase to their peak values and then decrease when V ? c_min. The effects of the flexural rigidity and the inertia of the plate are studied.     

Alternative state space formulations for magnetoelectric thermoelasticity with transverse isotropy and the application to bending analysis of nonhomogeneous plates

By introducing two displacement functions and two stress functions, the governing equations of the linear theory of magneto-electro-thermo-elasticity with transverse isotropy are simplified. On selecting certain physical quantities as the basic unknowns, two new state equations are established. Each of them is order reduced when compared with the one reported recently in literature, leading to a higher numerical efficiency. The material inhomogeneity along the axis of symmetry (z-direction) can be taken into account and an approximate laminate model is employed to facilitate deriving analytical solutions. The validity of new formulations is examined by considering a laminated magneto-electro-elastic rectangular plate and good agreement is obtained with existent results. A plate with a functionally graded property is then analyzed. The effect of magnetoelectric coupling in a BaTiO₃–CoFe₂O₄ composite predicted from the micromechanics simulation is studied quantitatively.     

Normal Mode Analysis of the High Speed Channel Flow in a Bidisperse Porous Medium

The linear Darcy–Brinkman model of the high speed flow in a bidisperse porous medium proposed by Nield and Kuznetsov (Transport Phenomena in Porous Media, 2005) is revisited in this paper. For the steady unidirectional flow in a parallel plane channel the exact analytical solutions for the fluid velocities are worked out by the normal-mode reduction of the governing equations. The limiting cases of the weak and strong momentum transfer between the flows in the fracture and porous phases are discussed in detail. A comparison to the nonlinear Forchheimer extension of the model proposed recently by Nield and Kuznetsov (Transport Porous Media, 2013) shows that, in the considered parameter range, the nonlinear effect of the Forchheimer drag is negligibly small. Even the simplest zero-momentum transfer solution yields an acceptable approximation.     

Extrema of Young’s modulus for cubic and transversely isotropic solids

For a homogeneous anisotropic and linearly elastic solid, the general expression of Young’s modulus E(n), embracing all classes that characterize the anisotropy, is given. A constrained extremum problem is then formulated for the evaluation of those directions n at which E(n) attains stationary values. Cubic and transversely isotropic symmetry classes are dealt with, and explicit solutions for such directions n are provided. For each case, relevant properties of these directions and corresponding values of the modulus are discussed as well. Results are shown in terms of suitable combinations of elements of the elastic tensor that embody the discrepancy from isotropy. On the basis of such material parameters, for cubic symmetry two classes of behavior can be distinguished and, in the case of transversely isotropic solids, the classes are found to be four. For both symmetries and for each class of behavior, some examples for real materials are shown and graphical representations of the dependence of Young’s modulus on direction n are given as well.     

The penny-shaped crack at a bonded plane with localized elastic non-homogeneity

This paper examines the axisymmetric problem pertaining to a penny-shaped crack which is located at the bonded plane of two similar elastic halfspace regions which exhibit localized axial variations in the linear elastic shear modulus, which has the form G(z)=G₁+G₂e^±ζz. The equations of elasticity governing this type of non-homogeneity are solved by employing a Hankel transform technique. The resulting mixed boundary value problem associated with the penny-shaped crack is reduced to a Fredholm integral equation of the second kind which is solved in a numerical fashion to generate the crack opening mode stress intensity factor at the tip.