Investigation of the stress singularity of a magnetoelectroelastic bonded antiplane wedge

The order of the stress singularity of a magnetoelectroelastic bonded antiplane wedge is analyzed by complex potential function and eigenfunction expansion method. Contrary to the familiar problem of elastic anisotropic bonded wedges which always produce real values for the order of singularity, the results of the magnetoelectroelastic bonded wedges may be real or complex. Numerical results are presented for problems with different boundary conditions. In particular, special behaviors of the order of the stress singularity for some degenerate composite materials and for some special wedge angles are noted.     

Fully nonlinear curvature flow of axially symmetric hypersurfaces with boundary conditions

Inspired by earlier results on the quasilinear mean curvature flow, and recent investigations of fully nonlinear curvature flow of closed hypersurfaces which are not convex, we consider contraction of axially symmetric hypersurfaces by convex, degree-one homogeneous fully nonlinear functions of curvature. With a natural class of Neumann boundary conditions, we show that evolving hypersurfaces exist for a finite maximal time. The maximal time is characterised by a curvature singularity at either boundary. Some results continue to hold in the cases of mixed Neumann–Dirichlet boundary conditions and more general curvature-dependent speeds.     

Three-dimensional contact problems with friction for a composite elastic wedge

Contact problems for a composite elastic wedge in the form of two joined wedge-shaped layers with different aperture angles joined by a sliding clamp, where the layer under the punch is incompressible, are studied in a three-dimensional formulation. Conditions for a sliding or rigid clamp or the absence of stresses are set up on one face of the composite wedge. The integral equations of the problems are derived taking account of the friction forces perpendicular to the edge of the wedge. The method of non-linear boundary integral equations of the Hammerstein type is used when the contact area is unknown. A regular asymptotic solution is constructed for an elliptic contact area. By virtue of the incompressibility of the material of the layer in contact with the punch, this solution retains the well known root singularity in the boundary of the contact area when account is taken of friction.     

Dirichlet-Neumann-impedance boundary value problems arising in rectangular wedge diffraction problems

Boundary value problems originated by the diffraction of an electromagnetic (or acoustic) wave by a rectangular wedge with faces of possible different kinds are analyzed in a Sobolev space framework. The boundary value problems satisfy the Helmholtz equation in the interior (Lipschitz) wedge domain, and are also subject to different combinations of boundary conditions on the faces of the wedge. Namely, the following types of boundary conditions will be under study: Dirichlet-Dirichlet, Neumann-Neumann, Neumann-Dirichlet, Impedance-Dirichlet, and Impedance-Neumann. Potential theory (combined with an appropriate use of extension operators) leads to the reduction of the boundary value problems to integral equations of Fredholm type. Thus, the consideration of single and double layer potentials together with certain reflection operators originate pseudo-differential operators which allow the proof of existence and uniqueness results for the boundary value problems initially posed. Furthermore, explicit solutions are given for all the problems under consideration, and regularity results are obtained for these solutions.

     

On the determination of crack initiation directions under mixed mode loading

An important issue in the mechanics of multilayered structures is the knowledge of local stress fields. Light‐weight heterogeneous structures with crack‐ or notch‐like discontinuities possess great danger of stress concentration and possible crack initiation leading to failure. In the present study, computational procedure, namely the finite element method, has been used to analyze bimaterial medium with various notch opening angles under arbitrary loading conditions. The prediction of the direction of crack initiation is important for wedge type constructions to gather a first hand knowledge of a potential damage zones which may undermine the integrity of the structure as a whole. The application of computational procedure to predict potential directions of crack initiation gives the flexibility required for different structural configurations, lamina orientations and incorporation of various boundary conditions. Superposition of various external loads to simulate pure mode I, mode II and mixed mode cases can also be carried out with excellent results in computational procedure. All the analysis results are based on the hypothesis of Erdogan and Sih. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stress singularity at the vertex of homogeneous and composite cones for different boundary conditions

A method of constructing singular solutions for conical bodies is proposed and versions of its numerical realization based on the finite-elements methods are considered. The results of numerical experiments illustrating the effectiveness and reliability of the proposed method are presented, as well as the results of a calculation of the eigenvalues, which define the form of the stress singularity for homogeneous and composite circular and non-circular cones for different boundary conditions.     

A fast high-order finite difference algorithm for pricing American options

We describe an improvement of Han and Wu’s algorithm [H. Han, X.Wu, A fast numerical method for the Black–Scholes equation of American options, SIAM J. Numer. Anal. 41 (6) (2003) 2081–2095] for American options. A high-order optimal compact scheme is used to discretise the transformed Black–Scholes PDE under a singularity separating framework. A more accurate free boundary location based on the smooth pasting condition and the use of a non-uniform grid with a modified tridiagonal solver lead to an efficient implementation of the free boundary value problem. Extensive numerical experiments show that the new finite difference algorithm converges rapidly and numerical solutions with good accuracy are obtained. Comparisons with some recently proposed methods for the American options problem are carried out to show the advantage of our numerical method.     

On corner singularities in Reissner-Mindlin's theory of elastic plates

In the framework of linear elasticity, singularities occur in domains with non-smooth boundaries. Particularly in Fracture Mechanics, the local stress field near stress concentrations is of interest. In this work, singularities at re-entrant corners or sharp notches in Reissner-Mindlin plates are studied. Therefore, an asymptotic solution of the governing system of partial differential equations is obtained by using a complex potential approach which allows for an efficient calculation of the singularity exponent λ. The effect of the notch opening angle and the boundary conditions on the singularity exponent is discussed. The results show, that it can be distinguished between singularities for symmetric and antisymmetric loading and between singularities of the bending moments and the transverse shear forces. Also, stronger singularities than the classical crack tip singularity with free crack faces are observed. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Superelements for the finite element solution of two-dimensional elliptic problems with boundary singularities

A novel singular superelement (SSE) formulation has been developed to overcome the loss of accuracy encountered when applying the standard finite element schemes to two-dimensional elliptic problems possessing a singularity on the boundary arising from an abrupt change of boundary conditions or a reentrant corner. The SSE consists of an inner region over which the known analytic form of the solution in the vicinity of the singular point is utilized, and a transition region in which blending functions are used to provide a smooth transition to the usual linear or quadratic isoparametric elements used over the remainder of the domain. Solution of the finite element equations yield directly the coefficients of the asymptotic series, known as the flux/stress intensity factors in linear heat transfer or elasticity theories, respectively. Numerical examples using the SSE for the Laplace equation and for computing the stress intensity factors in the linear theory of elasticity are given, demonstrating that accurate results can be attained for a moderate computational effort.     

带裂纹十次对称二维准晶平面弹性的无摩擦接触问题

借助经典平面弹性复变函数方法，研究了单个刚性凸基底压头作用下，带任意形状裂纹十次对称二维准晶半平面弹性的无摩擦接触问题.利用十次对称二维准晶位移、应力的复变函数表达式, 带任意形状裂纹的准晶半平面弹性无摩擦接触问题被转换为可解的解析函数复合边值问题，进而简化成一类可解的Riemann边值问题.通过求解Riemann边值问题，得到了应力函数的封闭解, 并给出了裂纹端点处应力强度因子和压头下方准晶体表面任意点处接触应力的显式表达式.从压头下方接触应力的表达式可以看出, 接触应力在压头边缘和裂纹端点处具有奇异性.当忽略相位子场影响时, 该文所得结论与弹性材料对应结果一致.数值算例分别给出了单个平底刚性压头无摩擦压入带单个垂直裂纹和水平裂纹的十次对称二维准晶下半平面的结果.该文所得结论为准晶材料的应用提供了理论参考.     

Analytical solution of the quasi-static thermoelasticity problem in a pressurized thick-walled cylinder subjected to transient thermal loading

Thermoelasticity problem in a thick-walled cylinder is solved analytically using the finite Hankel transform. Time-dependent thermal boundary conditions are assumed to act on the inner surface of the cylinder. For the mechanical boundary conditions two different cases are assumed: Traction–displacement problem (traction is prescribed on the inner surface and the fixed displacement boundary condition on the outer one) and Traction–Traction problem (tractions are prescribed on both the inner and outer surfaces of the hollow cylinder). The quasi-static solution of the thermoelasticity problem is derived analytically, i.e., the transient thermal response of the cylinder is derived and then, quasi-static structural problem is solved and closed form relations are extracted for the thermal stresses in the two problems. The results show to be in accordance with that cited in the literature in the special cases.     

An exact analysis for mode III cracks between two dissimilar magnetoelectroelastic layers

This paper develops a closed-form solution for an interface crack in a layered magnetoelectroelastic strip of finite width. The strip is subjected to anti-plane mechanical and in-plane electric and magnetic fields. Explicit expressions for the stress, electric, and magnetic fields, together with their intensity factors, are obtained for two extreme cases of an impermeable and a permeable cracks. The stress intensity factor does not depend on the electromagnetic boundary conditions assumed for the crack. However, the electrically and magnetically permeable boundary conditions on the crack profile have a significant influence on the crack-tip electromagnetic field intensity factors. Solutions for some special cases, such as a central crack, an edge crack, two symmetric collinear cracks, and a row of collinear interface cracks, are also obtained in closed forms. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 44, No. 6, pp. 763–784, November–December, 2008.     

蠕变损伤材料中缺口尖端稳定扩展的应力场

研究了在拉伸载荷和反平面载荷作用下蠕变损伤材料缺口尖端稳定扩展的应力场．假设材料的应力及位移势函数,得到了缺口尖端场的各参数分量,进而在小范围蠕变条件下,建立了缺口尖端稳定扩展的蠕变损伤控制方程,并考虑缺口尖端蠕变钝化作用和问题的边界条件,对控制方程进行了数值分析,得到了缺口尖端的应力场,并讨论了缺口尖端应力场随各影响参数的变化规律．结果表明,缺口尖端的应力具有r1/(1-n)的奇异性,应力率具有rn/(1-n)的奇异性,n是蠕变指数．     

Scaled boundary finite element analysis of stress singularities in piezoelectric multi-material systems

The scaled boundary finite element method (SBFEM) in an extension for piezoelectric materials is used to analyze twoand three-dimensional stress singularities in piezoelectric multi-material systems. It is found to be an efficient tool for the analysis of singularity orders of such situations, that turn out to be rather complex. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Geometrical embeddings of distributions into algebras of generalized functions

We use spectral theory to produce embeddings of distributions into algebras of generalized functions on a closed (compact without boundary) Riemannian manifold. These embeddings are invariant under isometries and preserve the singularity structure of the distributions (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

GLOBAL SOLUTIONS TO THE NAVIER-STOKES EQUATIONS FOR COMPRESSIBLE HEAT-CONDUCTING FLOW WITH SYMMETRY AND FREE BOUNDARY

ABSTRACT

Global solutions of the multidimensional Navier-Stokes equations for compressible heat-conducting flow are constructed, with spherically symmetric initial data of large oscillation between a static solid core and a free boundary connected to a surrounding vacuum state. The free boundary connects the compressible heat-conducting fluids to the vacuum state with free normal stress and zero normal heat flux. The fluids are initially assumed to fill with a finite volume and zero density at the free boundary, and with bounded positive density and temperature between the solid core and the initial position of the free boundary. One of the main features of this problem is the singularity of solutions near the free boundary. Our approach is to combine an effective difference scheme to construct approximate solutions with the energy methods and the pointwise estimate techniques to deal with the singularity of solutions near the free boundary and to obtain the bounded estimates of the solutions and the free boundary as time evolves. The convergence of the difference scheme is established. It is also proved that no vacuum develops between the solid core and the free boundary, and the free boundary expands with finite speed.     

Moving mesh method for problems with blow-up on unbounded domains

This paper studies the numerical solutions of semilinear parabolic partial differential equations (PDEs) on unbounded spatial domains whose solutions blow up in finite time. There are two major difficulties usually in numerical solutions: the singularity of blow-up and the unboundedness. We propose local absorbing boundary conditions (LABCs) on the selected artificial boundaries by using the idea of unified approach (Brunner et al., SIAM J Sci Comput 31:4478–4496, (2010). Since the uniform fixed spatial meshes may be inefficient, we adopt moving mesh partial differential equation (MMPDE) method to adapt the spatial mesh as the singularity develops. Combining LABCs and MMPDE, we can effectively capture the qualitative behavior of the blow-up singularities in the unbounded domain. Moreover, the implementation of the combination consists of two independent parts. Numerical examples also illustrate the efficiency and the accuracy of the new method.     

Exterior Wedge Diffraction Problems with Dirichlet,Neumann and Impedance Boundary Conditions

Classes of problems of wave diffraction by a plane angular screen occupying an infinite 270 degrees wedge sector are studied in a Bessel potential spaces framework. The problems are subjected to different possible combinations of boundary conditions on the faces of the wedge. Namely, under consideration there will be boundary conditions of Dirichlet-Dirichlet, Neumann-Neumann, Neumann-Dirichlet, impedance-Dirichlet, and impedance-Neumann types. Existence and uniqueness results are proved for all these cases in the weak formulation. In addition, the solutions are provided within the spaces in consideration, and higher regularity of solutions are also obtained in a scale of Bessel potential spaces.     

三维切口尖端应力应变场

本文利用双重幂级数展开法分析三维切口尖端应力应变奇异性，通过切口边界条件导出切口特征方程，进而求得不同切口内角下特征值序列解答，最后推得切口尖端应力应变场。     

A comparison theorem

The existence of an explosive singularity in the unsteady boundary layer on a rotating disc in a counter-rotating fluid has now been shown incontrovertibly following the work of K. Stewartson, C. J. Simpson, and R. J. Bodonyi (J. Fluid. Mech.121 (1982), 507–515). Here, we develop some asymptotic results for the governing differential equations to help gain an understanding of the mechanism behind this phenomenon. No definite conclusions are possible, but the presence of inertial oscillations at the edge of the boundary layer could well play a definite role.