A precise singular finite element for crack analysis

The multi-variable finite element algorithm based on the generalized Galerkin’smethod is more flexible to establish a finite element model in the continuum mechanies.Byusing this algorithm and numerical tests a new singular finite element for elasto-plasticfracture analysis has been formulated.The results of numerical tests show that the newelement possesses high accuracy and good performance.Some rules for formulating amulti-variable singular finite element are also discussed in this paper.     

A novel singular finite element of mixed-mode crack problems with arbitrary crack tractions

A novel singular finite element is presented to study cracked plates with arbitrary traction acting on crack surfaces. Firstly, the analytical solution around crack tips is determined using the symplectic dual approach. Subsequently, the solution is used to develop a novel singular finite element, which depicts accurately the characteristic of singular stresses field near crack tips. And the novel element can be applied to solve cracked plates, and both Mode I and Mode II stress intensity factors can be determined directly and accurately. Lastly, two numerical examples are given to illustrate the present method.     

压电材料轴对称有限元分析

给出了横观各向同性压电材料轴对称问题有限元法列式,叙述了压电材料有限元法分析中的一些特殊处理步骤,包括无量纲化和特征值问题计算中对电势自由度的凝聚。对压电材料轴对称问题的一些经典静力问题进行有限元计算,同时对压电圆板的轴对称自由振动进行了计算,和解析结果作了比较,二者符合良好。     

A finite element method on the plane elastic material analysis

Based on the finite element displacement method, a finite element method on the analysis of mechanical behaviour of plane elastic materials is proposed in this paper. By using this method and the corresponding computational program, the material behaviour of any unknown plane elastic material can be determined and all the elastic constants can be calculated.     

各向异性两相材料尖劈奇性场的非协调元分析

提出了一个基于位移的、分析柱状各向异性两相材料尖劈端部邻域的奇性位移场和应力场问题的非协调元特征分析法. 该方法从柱状扇区的散度定理出发，将柱状扇区控制方程的弱式化为一个与虚功原理相同形式的方程，采用一种新的非协调元技术把所导出的``虚功原理'转化为标准一阶特征方程的求解问题. 非协调元法中，尖劈端部邻域的位移场假定没有采用奇异变换技术，有限元的单元形式是一维的. 将柱状各向异性两相材料尖劈视为``广义平面应变'问题，位移场与坐标z无关，只关注界面端的幂奇异性而不考虑对数奇异性. 运用该方法给出了柱状各向异性两相材料尖劈端部奇性应力指数、奇性位移角分布和应力角分布的算例. 所有的计算结果表明，该方法使用的单元少而且精度较高.     

A novel hybrid finite element with a hole for analysis of plane piezoelectric medium with defects

A method is presented to formulate a novel hybrid finite element to obtain accurate distributions of mechanical and electrical quantities around a hole in plane piezoelectric mediums. The complex variable method is used in conjunction with Reissner's variational principle to formulate a hybrid special element with an elliptical hole. Detailed derivations are given and numerical examples are performed to demonstrate the accuracy and efficiency of the novel special element. Accurate results around the hole boundary are obtained for infinite and finite piezoelectric medium by using proposed special element.     

弹塑性板壳结构非线性有限元分析

本文根据塑性流动理论的基本公式，由隐式积分导出了与路径无关的变量更新算法和一致切线模量。采用单元广义应力应变直接离散塑性流动定律，构造了杂交应力单元一致切线刚度矩阵的显式表达式，编制了结构有限元程序ＳＡＦＥ，数值算例表明：本文的计算方法和计算程序是正确可靠的，可用于弹塑性板壳结构的非线性分析，计算结果屈曲临界载荷和极限承载能力。     

On compatible finite element models for elastic plastic analysis

Summary Some implicit or explicit assumptions on which the formulation of compatible finite element models for elastic plastic analysis is based are considered. It is shown that it is possible to formulate global elastic plastic element laws by independently assuming interpolation functions governing the displacement and the plastic strain fields. Some consistency requirements are discussed, which must be enforced in order to avoid artificial stiffening or fictitious redundancy. When these requirements are met, a «consistent» element is produced, whose properties are analyzed. It is also emphasized that, in spite of the use of compatible models, upper bounds to the collapse load are not always obtained for the limit analysis problem.For simplicity, beam elements are considered. However, the formulation holds for continua as well, as discussed in the last section.

---

Sommario Vengono analizzate alcune ipotesi esplicitmente o implicitamente assunte per la formulazione di elementi finiti usati nell'analisi elasto-plastica. Si mostra come l'assunzione indipendente di funzioni di interpolazione che governano i campi di spostamento e di deformazioni plastiche permetta di formulare legami costitutivi elasto-plastici globali per l'elemento. Vengono anche discusse alcune condizioni che devono essere rispettate per non irrigidire artificialmente l'elemento o non introdurvi una iperstaticità fittizia. Queste condizioni portano a definire elementi «consistenti», le cui proprietà vengono analizzate. Si fa anche notare come, nonostante l'uso di modelli congruenti, non vi sia sempre la garanzia di ottenere delimitazioni superiori al carico di collasso per problemi perfettamente plastici.Per semplicità, la trattazione viene svolta con riferimento ad elementi di trave. La possibilità di estensione al continuo è messa in luce nell'ultimo paragrafo.

---

---

Research supported by the C.N.R. (National (Italian) Research Council).     

复合材料中矩形夹杂角端部力学行为分析

提出了一种分析矩形夹杂角端部奇异应力场的新型杂交有限元方法,该方法在分析矩形夹杂角端部奇异应力场时,需要在夹杂端部构造一个超级单元。超级单元的刚度矩阵可以通过夹杂端部特征问题数值解建立。我们用这种方法计算了单向载荷作用下无限大均质板中单个矩形夹杂角端部奇异应力场,并与现有的数值解进行了比较。比较结果表明:本文提出的方法是可行的、有效的,而且数值结果精度高。为说明本文方法适用范围更广,文章最后讨论了各向异性弹性材料和横观各向同性压电材料中矩形夹杂角端部电弹性场行为。     

A transient finite element analysis of natural convection around a horizontal hot cylinder

This paper presents an advanced method for a 2-dimensional analysis of transient natural convection by finite element method. The present method, based on stream function—vorticity formulation, could get rid of numerical errors and constraint of perpendicular mesh subdivision, since we excluded a finite difference approximation of vorticity on no-slip boundaries. A considerable effect of upwind weighting function was examined. The method was successfully applied to a problem of natural convection around a horizontal hot cylinder.     

Plotting isoclinics for hybrid photoelasticity and finite element analysis

Displacement-based finite element method formulations are coupled with stress-based photoelasticity analysis. As the stress field is discontinuous at the interelement boundaries, the introduced smoothing procedure enables the generation of high-quality digital images acceptable for hybird experimental-numerical techniques. The proposed methods are applicable for the analysis of static and dynamic results of experimental photoelasticity.     

APPLICATION OF PENALTY FUNCTION METHOD IN ISOPARANIETRIC HYBRID FINITE ELEMENT ANALYSIS

IntroductionSince T.H.H.Pain firstly puts forward hybrid element method[1]in1964, the researchand application of hybrid element have got great development. T.H.H.Painet al.havemade important pioneer works in the domain of the research on incompatible displacementstructure hybrid elementstress pattern[2,3]in recentdecades. Reference [4] had provided theoptimizing design concept of hybrid element further and established standardization methodof incompatible displacement structure hybrid eleme…     

Notes on the finite element analysis of the axisymmetric elastic solid

A simple transformation of displacements considerably eases the explicit derivation of the finite element stiffness matrix for the axisymmetric elastic solid without causing a decline in the rate of convergence. The worsening of the condition of the global stiffness matrix caused by this transformation can be cured by scaling. A balanced numerical integration scheme maintaining the full rate of convergence is the one that integrates each term of the work and energy expressions to the order 2p ? 2, p being the degree of the complete polynomial in the shape functions.     

On “finite element analysis of Axisymmetric elastic body problems”

Regarding the calculation of the rigidity matrix of the linear triangular elements, there is really the eristence of the nonconvergent terms.     

The finite element method of singular perturbation problem

In this paper we construct new finite element subspace using polynomials of different degrees and the new finite element scheme is established. The convergence of the scheme and the stability of the reduced difference equation are proved.     

结构有限元分析中阶跃型弹性约束的一种有效处理方法

本文通过引入弹性约束刚度矩阵和结构位移约束列阵，提出了结构有限元分析中处理阶跃型弹性约束的一种有效方法。该法通过改变弹性约束系数及位移非约束量大小，可方便有效地处理结点常弹性约束，刚性约束，阶跃型弹性的约束及阶跃型刚性约束等问题。     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

三晶体中三晶交线附近弹性应变能集中的有限单元分析

提出应变能集中这一概念,对三晶体中三晶交线附近的弹性应变能的集中特性进行了计算,发现当加载方向平行于三晶交线时,三晶交线附近有弹性应变能集中,但是集中程度不是很高,当加载方向垂直于三晶交线时,根据三晶体中组元晶体的相对位置的不同和取向的不同,三晶交线附近或者没有弹性应变能集中,或者有弹性应变能集中,而且在一些情况下与平行于三晶交线方向加载的情况相比,弹性应变能的集中程度要高很多。