On a symplectic analytical singular element for cracks under thermal shock considering heat flux singularity

In a precise numerical modelling of cracks under thermal shock, the singularity issue resulted from heat flux should also be considered in addition to the one resulted from stress. The assumptions of constant temperature distribution usually adopted in the existing studies may lead to significant error. The concerned problem involves the discretization in both space and time domains. Numerical error resulted from the singularity issues in the space domain may be accumulated in the time domain. Hence, a unified framework which integrates reliable methods for both space and time domains are desired. In the present contribution, the classic thermal stress problem is restudied under the Hamiltonian system and the eigen functions are obtained analytically. A symplectic analytical singular element (SASE) for thermal stress analysis is reformulated based on the existing ones for thermal conduction and stress analyses. The singularity issues of both stress and heat flux are considered. A unified framework is formed with the precise time domain expanding algorithm (PTDEA) for the time domain and the formulated SASE for the space domain. A self-adaptive technique is used for the PTDEA to improve the numerical efficiency. The time dependent fracture parameters i.e., heat flux intensity factors (HFITs) and the mixed mode thermal stress intensity factors (TSIFs) can be solved accurately without any post-processing. Numerical examples are given for verification and validation of the proposed method.     

Asymptotics of the elastic solution of a plane problem near the crack tip terminating at a nonideal bimaterial interface (mode I and II)

The local stress-strain state in the vicinity of the crack tip in a composite is studied, taking into account the mechanical and geometric features of the nearest interface. The modeling of Mode I and II problems for a semi-infinite crack terminating normally at a nonideal interface in the bimaterial plane is considered. The constituents, of the composite are assumed to be elastic, homogeneous, and isotropic. The intermediate zone between the constituents is modeled by interfacial conditions in the form: [ _n ]=0, [u]=r _n , where [u] and [ _n ] are jumps of the vectors of displacements and tractions along the interface. The diagonal matrix with nonnegative components and the parameter, 0 are defined by the mechanical and geometric characteristics of the intermediate zone, respectively. Thus, the case =0 corresponds to the usual ideal contact conditions along the interface. Using the method of integral transformations, the corresponding problems are reduced to systems of functional equations, and later to systems of integral equations with fixed point singularities. The solvability of the systems of integral equations is proved and the asymptotics of their solutions is found. Based on these results, the local distributions of the displacements and stresses near the crack tip are obtained. It is shown that the interfacial parameters and greatly influence the stress not only qualitatively (the character of the stress singularity near the crack tip changes), but also quantitatively (number of singular terms in the asymptotics increases). The graphs illustrating these results are presented as the values of the interfacial parameters and , as well as the ratio of the shear moduli ₀/₁ of the constituents.Presented at the 10th International Conference on the Mechanics of Composite Materials (Riga, April 20–23, 1998).Polytechnical Institute Poland. Translated from Mekhanika Kompozitnykh Materialov, Vol. 34, No. 5, pp. 621–642, September–October, 1998.     

Problem of thermoelectromagnetoelasticity for a piezoelectric/piezomagnetic bimaterial with interface crack

We consider a plane strain problem for a piezoelectric/piezomagnetic bimaterial space with a crack in the region of the interface of the materials. At infinity, tensile and shear stresses and heat, electric, and magnetic flows are set. Using representations for all mechanical, thermal, and electromagnetic factors in terms of piecewise analytic functions, we formulate problems of linear conjugation that correspond to a model of an open crack and models taking into account the contact zone in the vicinity of a crack tip. Exact analytic solutions of the indicated problems are constructed. Expressions for stresses, the electric and magnetic inductions, jumps of derivatives of displacements, and electric and magnetic potentials on the interface are written. The coefficients of intensities of the indicated factors are presented. We derive a transcendental equation for the determination of the real length of the contact zone. The dependences of this length and the coefficients of intensity on the set external influences are investigated.     

A mathematical model and analytical solution for buckling of inclined beam-columns

The differential equation that governs the buckling behavior of an inclined beam-column is obtained using the energy method, and the use of a suitable change of variable reduces the various geometrical and physical parameters into a single dimensionless length parameter (X). An exact solution is presented via the use of some members of the family of generalized hypergeometric functions. One type of boundary condition (pinned-ends) is presented and analyzed, and others boundary conditions may be easily studied following the same process. The analysis shows the singular behavior of the inclined beam-column and demonstrates the procedure used to obtain the critical values of the axial force. An important application of this model is its usefulness in the analysis of buckling of drillstrings within curved boreholes.     

A finite element method for solving singular boundary-value problems

It is proved that under certain assumptions on the functions q(t) and f(t), there is one and only one function u₀(t) ∈ at which the functional

A finite element method for cohesive crack modelling

The present contribution is concerned with the computational modelling of cohesive cracks, whereby the discontinuity is not limited to interelement boundaries, but is allowed to propagate freely through the elements. Inelastic material behaviour is described by a discrete constitutive law, formulated in terms of tractions and displacements at the surface. Details on the implementation and numerical examples are given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A cohesive element approach for viscoelastic interface properties

This paper introduces a time dependent extension of the nonlinear traction separation law for an interface element. The constitutive equations of the load history dependent behaviour of the material are depicted and derived according to a generalized Maxwell‐model. This finite linear, viscoelastic approach allows the consideration of long term loading and time dependent material behaviour in thin layers. The implementation of the presented element formulation and the material approach are verified by numerical examples. The paper gives an outlook on further work and research topics. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A generalized cohesive element technique for arbitrary crack motion

A computational method for arbitrary crack motion through a finite element mesh, termed as the generalized cohesive element technique, is presented. In this method, an element with an internal discontinuity is replaced by two superimposed elements with a combination of original and imaginary nodes. Conventional cohesive zone modeling, limited to crack propagation along the edges of the elements, is extended to incorporate the intra-element mixed-mode crack propagation. Proposed numerical technique has been shown to be quite accurate, robust and mesh insensitive provided the cohesive zone ahead of the crack tip is resolved adequately. A series of numerical examples is presented to demonstrate the validity and applicability of the proposed method.     

A BE based finite element for crack propagation processes

A boundary element based finite macro element for the simulation of 3D crack propagation in the framework of linear elastic fracture mechanics is presented. While the major part of the numerical model is discretized with finite elements, a small domain containing the crack is meshed with boundary elements. By means of the Symmetric Galerkin BEM a stiffness formulation for the cracked BE domain is obtained which enables a direct FEM/BEM coupling. All necessary operations for the crack propagation are carried out within this boundary element based finite macro element and exploit the potential of the boundary integral formulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A partially penalty immersed interface finite element method for anisotropic elliptic interface problems

In this article, we develop a partially penalty immersed interface finite element (PIFE) method for a kind of anisotropy diffusion models governed by the elliptic interface problems with discontinuous tensor‐coefficients. This method is based on linear immersed interface finite elements (IIFE) and applies the discontinuous Galerkin formulation around the interface. We add two penalty terms to the general IIFE formulation along the sides intersected with the interface. The flux jump condition is weakly enforced on the smooth interface. By proving that the piecewise linear function on an interface element is uniquely determined by its values at the three vertices under some conditions, we construct the finite element spaces. Therefore, a PIFE procedure is proposed, which is based on the symmetric, nonsymmetric or incomplete interior penalty discontinuous Galerkin formulation. Then we prove the consistency and the solvability of the procedure. Theoretical analysis and numerical experiments show that the PIFE solution possesses optimal‐order error estimates in the energy norm and norm.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1984–2028, 2014     

A conductive crack and a remote electrode at the interface between two piezoelectric materials

An interaction of a tunnel conductive crack and a distant strip electrode situated at the interface between two piezoelectric semi-infinite spaces is studied. The bimaterial is subject by an in-plane electrical field parallel to the interface and by an anti-plane mechanical loading. Using the presentations of electromechanical quantities at the interface via sectionally-analytic functions the problem is reduced to a combined Dirichlet-Riemann boundary value problem. Solution of this problem is found in an analytical form excepting some one-dimensional integrals calculations. Closed form expressions for the stress, the electric field and their intensity factors, as well as for the crack faces displacement jump are derived. On the base of these presentations the energy release rate is also found. The obtained solution is compared with simple particular case of a single crack without electrode and the excellent agreement is found out. An auxiliary plane problem for open and closed cracks between two isotropic materials is also considered. The mathematical model of this problem is identical to the above one, therefore, the obtained solution is used for this model. It is compared with finite element solution of a similar problem and good agreement is found out.     

A cut finite element method for a Stokes interface problem

We present a finite element method for the Stokes equations involving two immiscible incompressible fluids with different viscosities and with surface tension. The interface separating the two fluids does not need to align with the mesh. We propose a Nitsche formulation which allows for discontinuities along the interface with optimal a priori error estimates. A stabilization procedure is included which ensures that the method produces a well conditioned stiffness matrix independent of the location of the interface.     

A double degenerated finite element for modeling the crack tip singularity

The objective of this paper is to propose a modified finite element called double quarter point finite element (DQPE) for modeling the singularity near the crack tip. Two techniques of evaluation (displacement correlation technique DCT and quarter point displacement technique QPDT) were used to estimate numerically the calibration factor for CN specimen. This study appears that the DQPE element is more effective than the QPE element. Not only that, but the length of the double quarter point finite element (DQPE) has little impact on the results. Numerical experiments are provided to demonstrate the effectiveness and robustness of the proposed element.     

A novel finite element model for helical springs

A general and accurate finite element model for helical springs subject to axial loads (extension or/and torsion) is developed in this paper. Due to the establishment of precise boundary conditions, only a slice of the wire cross-section needs to be modelled; hence, more accurate results can be achieved. An example application to a circular cross-sectional spring is analysed in detail.     

A second order isoparametric finite element method for elliptic interface problems

A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.     

A new robust C-type nonconforming triangular element for singular perturbation problems

In this paper, a new robust C⁰ triangular element is proposed for the fourth order elliptic singular perturbation problem with double set parameter method and bubble function technique, and a general convergence theorem for C⁰ nonconforming elements is presented. The convergence of the new element is proved in the energy norm uniformly with respect to the perturbation parameter. Numerical experiments are also carried out to demonstrate the efficiency of the new element.     

A symplectic analytical wave propagation model for damping and steady state forced vibration of orthotropic composite plate structure

An analytical wave propagation model is proposed in this paper for damping and steady state forced vibration of orthotropic composite plate structure by using the symplectic method. By solving an eigen-problem derived in the symplectic dual system of free bending vibration of orthotropic rectangular thin plates, the wave shape of plate is obtained in symplectic analytical form for any combination of simple boundary conditions along the plate edges. And then the specific damping capacity of wave mode is obtained symplectic analytically by using the strain energy theory. The steady state forced vibration of built-up plates structure is calculated by combining the wave propagation model and the finite element method. The vibration of the uniform plate domain of the built-up plates structure is described using symplectic analytical waves and the connector with discontinuous geometry or material is modeled using finite elements. In the numerical examples, the specific damping capacity of orthotropic rectangular thin plate with three different combinations of boundary condition is first calculated and analyzed. Comparisons of the present method results with respect to the results from the finite element method and from the Rayleigh–Ritz method validate the effectiveness of the present method. The relationship between the specific damping capacity of wave mode and that of modal mode is expounded. At last, the damped steady state forced vibration of a two plates system with a connector is calculated using the hybrid solution technique. The availability of the symplectic analytical wave propagation model is further validated by comparing the forced response from the present method with the results obtained using the finite element method.     

Notes on a novel finite element for anisotropy at large strains

A locking phenomenon can be observed in the case of anisotropic elasticity, due to high stiffness in preferred directions. In this contribution we propose a formulation with the goal to overcome this locking problem. We apply a Simplified Kinematic for the Anisotropic part of the free-energy by means of a constant approximation of the right Cauchy-Green tensor. For the tested boundary value problem the SKA-element performs excellent and behaves extremely robust. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Notes on a novel finite element for anisotropy at large strains

A locking phenomenon can be observed in the case of anisotropic elasticity, due to high stiffness in preferred directions. In this contribution we propose a formulation with the goal to overcome this locking problem. We apply a S implified K inematic for the A nisotropic part of the free-energy by means of a constant approximation of the right Cauchy-Green tensor. For the tested boundary value problem the SKA -element performs excellent and behaves extremely robust. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)