Solutions of thermoelastic crack problems in bonded dissimilar media or half-plane medium

The two-dimensional thermoelastic crack problem in bonded dissimilar media or in a half-plane medium is considered. The proposed method for solving this problem consists of two parts. In the first part, complex potential functions are derived which are enforced to satisfy the continuity conditions across the interface, while the second part consists of the derivation of singular integral equations by introducing the dislocation functions along the crack border which are solved numerically. For both half-plane and two bonded half-plane problems associated with an insulated crack, the thermal stress intensity factors are computed numerically by using the appropriate interpolation formulae. The results compared with those of the homogeneous case given in the literature show that the method proposed here is effective, simple and general.     

General solutions of coupled thermoelastic problem

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎｔｈｅｈｅａｔｃｏｎｄｕｃｔｉｏｎｅｑｕａｔｉｏｎｉｎｖｏｌｖｅｓｔｈｅｔｅｒｍｏｆｄｅｆｏｒｍａｔｉｏｎａｎｄｔｈｅｒｍｏｅｌａｓｔｉｃｅｑｕａｔｉｏｎｓｃｏｎｔａｉｎｔｅｍｐｅｒａｔｕｒｅ,ｔｈｉｓｋｉｎｄｏｆｐｒｏｂｌｅｍｉｓｃａｌｌｅｄｃｏｕｐｌｅｄｔｈｅｒｍｏｅｌａｓｔｉｃｐｒｏｂｌｅｍ,ｗｈｉｃｈｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄａｎｄｅｌａｓｔｉｃｆｉｅｌｄｍｕｓｔｂｅｓｏｌｖｅｄｓｉｍｕｌｔａｎｅｏｕｓｌｙ[1,2].Ｉｎｇｅｎｅｒａｌ,ｔｈｅｅｑｕａｔｉｏ…     

Two-dimensional fractal-like finite element method for thermoelastic crack analysis

The fractal-like finite element method has been proved to be very efficient and accurate in two-dimensional static and dynamic crack problems. In this paper, we extend our previous study to include the thermal effect for two-dimensional isotropic thermal crack problems. Both the temperature intensity factor and thermal stress intensity factor can be calculated directly. The temperature distribution is first found, which is imposed thereafter as a thermal load in the elastic problem. The transformation function used in the study has been found analytically. The effects of different thermal loading on the thermal stress intensity factor are presented. The numerical examples are compared with the results from other methods and find to be in good agreement.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

Some consequences of the inequality conditions in contact and crack problems

The importance of the inequalities and related side conditions that must be incorporated in contact and crack problems is emphasized and the ensuing consequences explored. An asymptotic analysis of the transitions from slip to separation, stick to slip, and stick to separation is carried out. The inequalities in contact problems make the contact pressure continuous for all levels of friction. They also make a direct transition from stick to separation impossible, unless the combination of materials is special. The inequalities in crack problems are less stringent, but they preclude certain singularities that appear to have flourished in the literature previously.

---

Résumé On appuie sur l'importance d'inéquations et d'autres conditions auxiliaires qui doivent être comprises dans les problèmes de contact et de fissure, et on en explore les conséquences. On emploie une analyse asympotique sur les transitions entre les zones glissement-décollement, adhérence-glussement et adhérence-décollement. Il s'ensuit que la contrainte normale du contact doit être continue pour toutes les valeurs du frottement. De plus, la transition directe adhérence-décollement est impossible, à moin que la combinaison des matériaux ne soit exceptionelle. Les inéquations dans les problèmes de fissure sont moins fortes, mais elles sont cependant suffisantes pour empêcher l'existence de certaines singularités qui apparaissent souvent dans les études précedentes.

     

热弹性平面问题的规则化边界积分方程

对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。     

On coupled mean fields in heterogeneous linear thermoelastic media

The governing equation of the coupled mean fields induced in the heterogeneous linear ther-moelastic media was derived. Discussion was made on the nonlocality of the governing equation. The relation of the mean fields was investigated in the uncoupled case when the acceleration term was disregarded.     

热冲击下物性与温度相关性对热弹性响应的渐进分析

计及材料物性与温度的相关性,基于Green-Naghdi能量无耗散广义热弹性理论(G-N II理论),对热冲击下具有变物性特征材料的热弹性响应进行了求解分析。借助Laplace正、反变换技术以及Krichhoff变换,在热物性参数随真实温度呈线性规律的前提下,推导了半无限大体受热冲击作用时热弹性响应的解析表达式,通过求解分析,得到了热冲击下热波、热弹性波的传播规律,位移场、温度场以及应力场的分布情况,以及物性随温度相关性对热弹性响应的影响效果。结果表明：当考虑材料物性随温度的变化时,热波、热弹性波的传播以及各物理场的分布均受到不同程度的影响,且物性随温度相关性对热弹性响应的作用效果将受到材料热-力耦合特性的影响。     

A coupled FEM/BEM approach and its accuracy for solving crack problems in fracture mechanics

The finite element (FEM) and the boundary element methods (BEM) are well known powerful numerical techniques for solving a wide range of problems in applied science and engineering. Each method has its own advantages and disadvantages, so that it is desirable to develop a combined finite element/boundary element method approach, which makes use of their advantages and reduces their disadvantages. Several coupling techniques are proposed in the literature, but until now the incompatibility of the basic variables remains a problem to be solved. To overcome this problem, a special super-element using boundary elements based on the usual finite element technique of total potential energy minimization has been developed in this paper. The application of the most commonly used approaches in finite element method namely quarter-point elements and J-integrals techniques were examined using the proposed coupling FEM–BEM. The accuracy and efficiency of the proposed approach have been assessed for the evaluation of stress intensity factors (SIF). It was found that the FEM–BEM coupling technique gives more accurate values of the stress intensity factors with fewer degrees of freedom.     

Nonlinear acceleration of coupled fluid–structure transient thermal problems by Anderson mixing

Conjugate heat‐transfer problems are typically solved using partitioned methods where fluid and solid subdomains are evaluated separately by dedicated solvers coupled through a common boundary. Strongly coupled schemes for transient analysis require fluid and solid problems to be solved many times each time step until convergence to a steady state. In many practical situations, a fairly simple and frequently employed fixed‐point iteration process is rather ineffective; it leads to a large number of iterations per time step and consequently to long simulation times. In this article, Anderson mixing is proposed as a fixed‐point convergence acceleration technique to reduce computational cost of thermal coupled fluid–solid problems. A number of other recently published methods with applications to similar fluid–structure interaction problems are also reviewed and analyzed. Numerical experiments are presented to illustrate relative performance of these methods on a test problem of rotating pre‐swirl cavity air flow interacting with a turbine disk. It is observed that performance of Anderson mixing method is superior to that of other algorithms in terms of total iteration counts. Additional computational savings are demonstrated by reusing information from previously solved time steps. Copyright © All rights reserved 2012 Rolls‐Royce plc.