Mixed-mode stress intensity factors of 3D interface crack in fully coupled electromagnetothermoelastic multiphase composites

This contribution presents an extended hypersingular intergro-differential equation (E-HIDE) method for modeling the 3D interface crack problem in fully coupled electromagnetothermoelastic anisotropic multiphase composites under extended electro-magneto-thermo-elastic coupled loads through theoretical analysis and numerical simulations. First, based on the extended boundary element method, the 3D interface crack problem is reduced to solving a set of E-HIDEs coupled with extended boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended singular stress indices around the interface crack front terminating at the interface is analyzed by the extended main-part analysis. The extended stress intensity factors near the crack front are defined. In addition, a numerical method for a 3D interface crack problem subjected to extended loads is proposed, in which the extended displacement discontinuities are approximated by the product of basic density functions and polynomials. Finally, the radiation distribution of extended stress intensity factors at the interface crack surface are calculated, and the results are presented toward demonstrating the applicability of the proposed method.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

Analysis of an interfacial crack in a piezoelectric bi-material via the extended Green’s functions and displacement discontinuity method

Based on the extended Stroh formalism, we first derive the extended Green’s functions for an extended dislocation and displacement discontinuity located at the interface of a piezoelectric bi-material. These include Green’s functions of the extended dislocation, displacement discontinuities within a finite interval and the concentrated displacement discontinuities, all on the interface. The Green’s functions are then applied to obtain the integro-differential equation governing the interfacial crack. To eliminate the oscillating singularities associated with the delta function in the Green’s functions, we represent the delta function in terms of the Gaussian distribution function. In so doing, the integro-differential equation is reduced to a standard integral equation for the interfacial crack problem in piezoelectric bi-material with the extended displacement discontinuities being the unknowns. A simple numerical approach is also proposed to solve the integral equation for the displacement discontinuities, along with the asymptotic expressions of the extended intensity factors and J-integral in terms of the discontinuities near the crack tip. In numerical examples, the effect of the Gaussian parameter on the numerical results is discussed, and the influence of different extended loadings on the interfacial crack behaviors is further investigated.     

Some compression flows in nonaxisymmetric ring nozzles

Busemann's problem concerning fully developed conical flow in an axisymmetric nozzle of special type is extended to include certain nonaxisymmetric ring nozzles. The constructed flows contain strong discontinuities in the form of developable surfaces (in Busemann's solution, strong discontinuities have the form of a circular-cone surface).     

An h-adaptive modified element-free Galerkin method

In this work an h-adaptive Modified Element-Free Galerkin (MEFG) method is investigated. The proposed error estimator is based on a recovery by equilibrium of nodal patches where a recovered stress field is obtained by a moving least square approximation. The procedure generates a smooth recovered stress field that is not only more accurate then the approximate solution but also free of spurious oscillations, normally seen in EFG methods at regions with high gradient stresses or discontinuities.The MEFG method combines conventional EFG with extended partition of unity finite element (EPUFE) methods in order to create global shape functions that allow a direct imposition of the essential boundary conditions.The re-meshing of the integration mesh is based on the homogeneous error distribution criterion and upon a given prescribed admissible error. Some examples are presented, considering a plane stress assumption, which shows the performance of the proposed methodology.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.     

Accounting for local capillary effects in two-phase flows with relaxed surface tension formulation in enriched finite elements

This paper introduces a numerical method able to deal with a general bi-fluid model integrating capillary actions. The method relies first on the precise computation of the surface tension force. Considering a mathematical transformation of the surface tension virtual work, the regularity required for the solution on the evolving curved interface is weakened, and the mechanical equilibrium of the triple line can be enforced as a natural condition. Consequently, contact angles of the liquid over the solid phase result naturally from this equilibrium. Second, for an exhaustive representation of capillary actions, pressure jumps across the interface must be accounted for. A pressure enrichment strategy is used to properly compute the discontinuities in both pressure and gradient fields. The resulting method is shown to predict nicely static contact angles for some test cases, and is evaluated on complex 3D cases.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

Numerical investigations of the XFEM for solving two-phase incompressible flows

This paper discusses the application of the extended finite element method (XFEM) to solve two-phase incompressible flows. The Navier–Stokes equations are discretised using the Taylor–Hood finite element. To capture the different discontinuities across the interface, kink or jump enrichments are used for the velocity and/or pressure fields. However, these enrichments may lead to an inappropriate combination of interpolations. Different polynomial enrichment orders and different enrichment functions are investigated; only the stable combination will be used afterward.

In cases with a surface tension force, the accuracy mainly relies on the precise computation of the normal and curvature. A novel method for computing normal vectors to the interface is proposed. This method employs successive mesh refinements inside the cut elements. Comparisons with analytical and numerical solutions demonstrate that the method is effective. Moreover, the mesh refinement improves the sub-integration in the XFEM and allows for a precise re-initialisation procedure.     

On a Sharp Sobolev‐Type Inequality on Two-Dimensional Compact Manifolds

Motivated by the asymptotic analysis of double vortex condensates in the Chern‐Simons‐Higgs theory, we construct a suitable minimizing sequence for a sharp Sobolev inequality “à la Moser” for two‐dimensional compact manifolds. As a consequence, we first obtain a direct proof of the sharp character of such an inequality. Secondly, and more interestingly, we use such minimizing sequence to show that for the flat torus the corresponding extremal problem attains its infimum. (Accepted April 6, 1998)     

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.     

Relative Entropy and the Stability of Shocks and Contact Discontinuities for Systems of Conservation Laws with non-BV Perturbations

We develop a theory based on relative entropy to show the uniqueness and L ² stability (up to a translation) of extremal entropic Rankine?CHugoniot discontinuities for systems of conservation laws (typically 1-shocks, n-shocks, 1-contact discontinuities and n-contact discontinuities of large amplitude) among bounded entropic weak solutions having an additional trace property. The existence of a convex entropy is needed. No BV estimate is needed on the weak solutions considered. The theory holds without smallness conditions. The assumptions are quite general. For instance, strict hyperbolicity is not needed globally. For fluid mechanics, the theory handles solutions with vacuums.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

Extended three-dimensional digital image correlation (X3D-DIC)

A correlation algorithm is proposed to measure full three-dimensional displacement fields in a three-dimensional domain. The chosen kinematic basis for this measurement is based on continuous finite-element shape functions. It is furthermore proposed to account for the presence of strong discontinuities, similarly to extended finite element schemes, with a suited enrichment of the kinematics with discontinuities supported by a (crack) surface. An optimization of the surface geometry is proposed based on correlation residuals. The procedure is applied to analyze one loading cycle of a fatigue-cracked nodular graphite cast iron sample by using computed tomography pictures. Subvoxel crack openings are revealed and measured. To cite this article: J. Réthoré et al., C. R. Mecanique 336 (2008).     

Some extremum properties of finite-step solutions inelastoplasticity

For the nonholonomic elastic–plastic problem under a given external action history overa time interval, an extremal formulation is given in terms of the complete solution over the wholeinterval. The assumed elastic–plastic behaviour is of the associated type with piecewiselinearized yield surface and linear hardening.When the loading history is reduced to an infinitesimal increment of the external actions(incremental problem) or when the material behaviour is assumed to be of the holonomic type (finite holonomic step) problem, the functional of the extremal formulation may be split into thesum of two other simpler functionals (previously introduced) whose minimum, for both of them,gives the problem solution under less constraints than in the original problem.For general non-holonomic loading histories the above splitting is shown to be still possiblewhen a particular change of the complementarity condition of the constitutive law is considered,which leads to a new class of holonomic problems.It is shown that some problems of this new class, together with a suitable time discretization,represent the schematization of the original problem corresponding to well known numericalintegration schemes.     

Closed form solutions of Euler–Bernoulli beams with singularities

The problem of the integration of the static governing equations of the uniform Euler–Bernoulli beam with discontinuities is studied. In particular, two types of discontinuities have been considered: flexural stiffness and slope discontinuities. Both the above mentioned discontinuities have been modeled as singularities of the flexural stiffness by means of superimposition of suitable distributions (generalized functions) to a uniform one dimensional field. Closed form solutions of governing differential equation, requiring the knowledge of the boundary conditions only, are proposed, and no continuity conditions are enforced at intermediate cross-sections where discontinuities are shown. The continuity conditions are in fact embedded in the flexural stiffness model and are automatically accounted for by the proposed integration procedure. Finally, the proposed closed form solution for the cases of slope discontinuity is compared with the solution of a beam having an internal hinge with rotational spring reproducing the slope discontinuity.     

Exact Solution of the Problem on a Six‐Constant Jeffreys Model of Fluid in a Plane Channel

An exact analytic solution of the problem of a generalized viscoelastic Jeffreys fluid flow in a plane channel under the action of a pressure gradient is found. The velocity profiles are obtained in a parametric form with a velocity gradient taken as a parameter. The critical values of the pressure gradient are determined, which, when exceeded, lead to weak tangential discontinuities in the longitudinal velocity profile. When the pressure gradient changes smoothly over some range of parameters, a hysteresis loop emerges on the graph of the flow rate versus the pressure gradient.     

计及表面能的应变梯度理论本构参数识别

本文提出利用静态位移信息对一种计及表面能的应变梯度理论本构参数进行识别的求解策略.基于Vardoulakis和Sulem的计及表面能的简单线性应变梯度理论,文献[13]给出了伯努力-欧拉梁弯曲问题的正演解析模型,本文将其反演归结为两个带有不等式约束的非线性规划问题.在此基础上,采用黄金分割一维搜索方法进行求解,给出了数值验证,讨论了信息误差对反演结果的影响.结果显示,这种方法可以用来对应变梯度理论本构参数进行识别,即使在体积和表面能常数非常小的情况下,仍然能够得到满意的结果.     

Numerical investigation of the slope discontinuities in large deformation finite element formulations

In many multibody system applications, the system components are made of structural elements that can have different orientations, leading to slope discontinuities. In this paper, a numerical investigation of a new procedure that can be used to model structures with slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is presented. This procedure can be applied to model slope discontinuities in the case of commutative rotations of gradient deficient elements that are used for modeling thin beam and plate structures. An important special case to which the proposed procedure can be applied is the case of all planar gradient deficient ANCF finite elements. The use of the proposed method leads to a constant orthogonal element transformation that describes an arbitrary initial configuration. As a consequence, one obtains, in the case of large commutative rotations and large deformations, a constant mass matrix for structures which have complex geometry. The procedure used in this investigation to model slope discontinuities requires the use of the concept of the intermediate finite element coordinate system. For each finite element, a new set of gradient coordinates that define, at the discontinuity node, the element deformation with respect to the intermediate element coordinate system is introduced. These new gradient coordinates are assumed to be equal for the two finite elements at the point of intersection. That is, the change of the gradients of two elements at the intersection point from their respective intermediate initial reference configuration is assumed to be the same. This procedure leads to a set of linear algebraic equations that define the orthogonal transformation matrix for the finite element. Numerical examples are presented in order to demonstrate the use of the proposed procedure for modeling slope discontinuities.     

Some problems of optimizing shell shape and thickness distribution on the basis of a genetic algorithm

We consider problems related to designing axisymmetric shells of minimal weight (mass) and the development of efficient nonlocal optimization methods. The optimization problems under study consist in simultaneous search for the optimal geometry and the shell thickness optimal distribution from the minimal weight condition under strength constraints and additional geometric constraints imposed on the thickness function, the transverse cross-section radii distribution, and the volume enclosed by the shell. Using the method of penalty functions, we reduce the above optimal design problem to a nonconvex minimization problem for the extended Lagrange functional. To find the global optimum, we apply an efficient genetic algorithm. We present the results of numerical solution of the optimal design problem for dome-like shells of revolution under the action of gravity forces. We present some data characterizing the convergence of the method developed here.