THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓ,ｔｈｅｄｉｓｐｌａｃｅｍｅｎｔｅｌｅｍｅｎｔｓａｒｅｖｅｒｙｐｏｐｕｌａｒｂｅｃａｕｓｅｏｆｔｈｅｉｒｅａｓｙｃｏｎｓｔｒｕｃｔｉｏｎ .Ｈｏｗｅｖｅｒ,ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｓｓｕｍｅｔｈｅｃｏｎｔｉｎｕｏｕｓｄｉｓｐｌａｃｅｍｅｎｔｆｉｅｌｄｆｏｒｓａｔｉｓｆａｃｔｏｒｙｓｈｅｌｌｅｌｅｍｅｎｔｓｔｈａｔｉｓｆｒｅｅｆｒｏｍｓｈｅａｒｌｏｃｋｉｎｇ .Ｉｎ 1 964 ,ｔｈｅｈｙｂｒｉｄｓｔｒｅｓｓｅｌｅｍｅｎｔｗ…     

构造杂交应力单元的柔度矩阵H对角化方法

证明了杂交元柔度矩阵 H非奇异的充分必要条件是假设应力模式线性无关 ;以及等价应力模式形成相同的杂交元。在此基础上建立了假设应力模式的 Hilbert子空间 ,从而可以利用 Schmidt方法简单地得到等价的正交应力模式 ,实现了柔度矩阵 H对角化 ,使得杂交元形成过程中完全避免了繁杂的矩阵求逆运算 ,提高了杂交元分析的计算效率 ,特别在柔度矩阵不容易显式求逆的材料非线性分析中更具有实际意义     

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The resulting orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. Moreover, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.     

杂交元本征应力模式和应力子空间的性质研究

详细讨论了有限元本征应力模式和应力子空间的性质，并着重讨论和进一步完善了与杂交应力有限元应力子空间有关的一些定理，为提出新方法提供了理论基础，主要包括：（1）证明了杂交元特征值不大于对应位移元的特征值；（2）证明了矩阵H非奇异的充分必要条件是假设应力模式线性无关；（3）证明了杂交元所对应位移元的本征应力模式形成的杂交元与该位移元相同；（4）证明了等价假设应力模式形成相同的杂交元；（5）证明了确定杂交元本征应力模式的充分必要条件是其范数平方等于所形成杂交元的变形模态特征值；（6）证明了杂交元假设应力模式与变形模态的能量一一对应的充分必要条件是假设应力模式彼此正交且与所对应位移元的本征应力模式除了一一对应者之外都正交。     

A new hybrid stress element method for the analysis of elastoplastic solids

Based on a modified Hellinger/Reissner variational principle which includes the equivalent stress, equivalent plastic strain and non-conforming displacement increments as independent variables, a quadrilateral isoparametric hybrid stress element for the analysis of elastoplastic problem is proposed. By this formulation, the yield criterion and flow rule are satisfied in an average sense and greater accuracy can be obtained by using non-conforming displacement. A numerical example is presented to show that the present model has high accuracy and computational effectiveness.This project is supported by the Natural Science Foundation of the State Education Commission.     

DEFORMATION RIGIDITY OF ASSUMED STRESS MODES IN HYBRID ELEMENTS

The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.     

Coarse-mesh-accuracy improvement of bilinear Q<Subscript>4</Subscript>-plane element by the combined hybrid finite element method

The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q₄-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i. e., the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q₄-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q₄-element. Foundation item: the National Tianyuan Yourth Funds of China (TY10126027) Biography: Xie Xiao-ping (1970≈)     

Computation of super-convergent nodal stresses of timoshenko beam elements by EEP method

The newly proposed element energy projection (EEP) method has been applied to the computation of super-convergent nodal stresses of Timoshenko beam elements. Generalformul as based on element projection theorem were derived and illustrative numerical examples using two typical elements were given. Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions. The EEP method gives super-convergent nodal stresses, which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude. And in addition, it can overcome the “ shear locking“ difficulty for stresses even when the displacements are badly affected. This research paves the way for application of the EEP method to general onedimensional systems of ordinary differential equations.     

杂交元零能模式抑制的改进基本变形模式方法

基本变形模式方法是构造杂交元假设应力场的一种有效方法,但只给出了识别零能模式的必要条件.论文指出基本变形模式方法由于没有考虑耦合影响而存在应用上的局限,在此基础上利用柔度矩阵的正定性给出了判别零能机动模式的充要条件,并且通过内积向量无关性进一步完善了假设杂交元应力场的基本变形模式方法,可以适用于基本变形模式之间耦合的任...     

Stress analysis for an infinite strip weakned by periodic cracks

Stress analysis for an infinite stripcracks were assumed in a horizontal position,weakened by periodic cracks is studied. The and the strip was applied by tension “p“ in y-direction. The boundary value problem can be reduced into a complex mixed one. It is found that the EEVM ( eigenfunction expansion variational method) is efficient to solve the problem. The stress intensity factor at the crack tip and the T-stress were evaluated. From the deformation response under tension the cracked strip can be equivalent to an orthotropic strip without cracks. The elastic properties in the equivalent orthotropic strip were also investigated. Finally, numerical examples and results were given.     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli’s gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CON-TINUUM THEORIES (Ⅱ)—MICROMORPHIC CONTINUUM THEORY AND COUPLE STRESS THEORY

The purpose is to reestablish the balance laws of momentum, angular momentum and energy and to derive the corresponding local and nonlocal balance equations for micromorphic continuum mechanics and couple stress theory. The desired results for micromorphic continuum mechanics and couple stress theory are naturally obtained via direct transitions and reductions from the coupled conservation law of energy for micropolar continuum theory, respectively. The basic balance laws and equations for micromorphic continuum mechanics and couple stress theory are constituted by combining these results derived here and the traditional conservation laws and equations of mass and microinertia and the entropy inequality. The incomplete degrees of the former related continuum theories are clarified. Finally, some special cases are conveniently derived. Foundation items: the National Natural Science Foundation of China (10072024); the Research Foundation of Liaoning Education Committee (990111001) Biography: DAI Tian-min (1931≈)     

构造满足特征值要求的杂交元应力子空间

通过引进适当参数给出了杂交元应力空间特征值与参数之间的关系,从而可以通过调节参数来构造满足特征值即模态刚度要求的杂交元应力子空间。本在位移元本征应力模式基础上引进调节参数,同时,利用矩阵H对角化方法计算杂交元应力子空间的本征应力模式,然后由此方便有效地计算特征值,从而大大提高了计算效率。本通过建立Q4杂交应力元特征值与参数之间关系说明了这一方法是确定可行的。     

Mixed energy method for solution of quadratic programming problems and elastic-plastic analysis of truss structures

A new algorithm for the solution of quadratic programming problems is put forward in terms of the mixed energy theory and is further used for the incremental solution of elastic-plastic truss structures. The method proposed is different from the traditional one, for which the unknown variables are selected just in one class such as displacements or stresses. The present method selects the variables in the mixed form with both displacement and stress. As the method is established in the hybrid space, the information found in the previous incremental step can be used for the solution of the present step, making the algorithm highly efficient in the numerical solution process of quadratic programming problems. The results obained in the examples of the elastic-plastic solution of the truss structures verify what has been predicted in the theoretical analysis. Project supported by the National Natural Science Foundation of China (No. 50178916, No. 19732020 and No. 19872016), the National Key Basic Research Special Foundation (No. G1999032805), the Special Funds for Major State Basic Research Projects and the Foundation for University Key Teachers by the Ministry of Education of China.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅲ) —NOETHER’S THEOREM

The existing various couple stress theories have been carefully restudied.The purpose is to propose a coupled Noether’s theorem and to reestablish rather complete conservation laws and balance equations for couple stress elastodynamics. The new concrete forms of various conservation laws of couple stress elasticity are derived. The precise nature of these conservation laws which result from the given invariance requirements are established. Various special cases are reduced and the results of micropolar continua may be naturally transited from the results presented in this paper.     

Coupling of assumed stress finite element and boundary element methods with stress-traction equilibrium

In this study, the stress based finite element method is coupled with the boundary element method in two different ways. In the first one, the ordinary distribution matrix is used for coupling. In the second one, the stress traction equilibrium is used at the interface line of both regions as a new coupling process. This new coupling procedure is presented without a distribution matrix. Several case studies are solved for the validation of the developed coupling procedure. The results of case studies are compared with the distribution matrix coupling, displacement based finite element method, assumed stress finite element method, boundary element method, ANSYS and analytical results whenever possible. It is shown that the coupling of the stress traction equilibrium with assumed stress finite elements gives as accurate results as those by the distribution matrix coupling.     

Dynamical Behavior of Viscoelastic Cylindrical Shells Under Axial Pressures

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｒｙｏｆｖｉｓｃｏｅｌａｓｔｉｃｉｔｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｓｈａｖｅｂｅｅｎｂｅｃｏｍｉｎｇａｖｅｒｙａｔｔｒａｃｔｉｖｅｒｅｓｅａｒｃｈｆｉｅｌｄｓｉｎｃｅｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｗｅｒｅｗｉｄｅｌｙａｐｐｌｉｅｄｉｎｔｏｖａｒｉｏｕｓｅｎｇｉｎｅｅｒｉｎｇ .Ｂｕｔｔｏｔｈｅａｕｔｈｏｒｓ’ｋｎｏｗｌｅｄｇｅ,ｔｈｅｒｅａｒｅｒａｒｅｆｏｒｐａｐｅｒｓｃｏｎｃｅｒｎｉｎｇｗｉｔｈｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒｏｆｖｉｓｃｏｅ…     

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, KⅠ and KⅡ, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.     

ROBUST CONTROL OF PERIODIC BIFURCATION SOLUTIONS

The topological bifurcation diagrams and the coefficients of bifurcation equation were obtained by C-L method. According to obtained bifurcation diagrams and combining control theory, the method of robust control of periodic bifurcation was presented, which differs from generic methods of bifurcation control. It can make the existing motion pattern into the goal motion pattern. Because the method does not make strict requirement about parametric values of the controller, it is convenient to design and make it. Numerical simulations verify validity of the method.     

SINGULAR ANALYSIS OF BIFURCATION OF NONLINEAR NORMAL MODES FOR A CLASS OF SYSTEMS WITH DUAL INTERNAL RESONANCES

ＩｎｔｒｏｄｕｃｔｉｏｎＡｎｉｎｔｅｒｅｓｔｉｎｇｆｅａｔｕｒｅｉｎｔｈｅｆｒｅｅｖｉｂｒａｔｉｏｎｏｆａｎｏｎｌｉｎｅａｒｓｙｓｔｅｍｉｓｔｈｅｆａｃｔｔｈａｔｔｈｅｎｕｍｂｅｒｏｆｅｘｉｓｔｉｎｇｎｏｒｍａｌｍｏｄｅｓｍａｙｅｘｃｅｅｄｔｈｅｎｕｍｂｅｒｏｆｄｅｇｒｅｅｓｏｆｆｒｅｅｄｏｍ ,ａｐｈｅｎｏｍｅｎｏｎｎｏｔｅｎｃｏｕｎｔｅｒｅｄｉｎａｌｉｎｅａｒｓｙｓｔｅｍａｎｄｃａｕｓｅｄｂｙｍｏｄｅｂｉｆｕｒｃａｔｉｏｎ .Ｔｈｅｒｅｆｏｒｅｍｕｃｈｗｏｒｋｈａｓｂｅｅｎｄｏｎｅｏｎｔｈｅｓｔ…