非连续问题中单元分割的模板方法

扩展有限元法 (extended finite element method, XFEM) 因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一. 但该方法在工程实践中存在单元被裂纹分割的几何困难 —— 现有精确几何分割方法实现复杂、计算量大、鲁棒性差. 为克服这一困难, 本文提出一种基于单元水平集的模板分割方法, 用于非连续单元子剖分和数值积分. 首先, 遍历单元水平集值所有形态并建立标准单元分割模板库; 然后, 根据单元水平集值, 对非标准单元进行形态查询和模板插值; 最后, 套用标准单元分割模板实现单元高效分割和子剖分. 将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中. 算例分析表明, 本文提出的模板分割方法具有较高计算精度. 由于不引入复杂几何操作, 该模板分割方法同时具有较高计算效率和鲁棒性, 故可为XFEM类方法在实际工程应用中提供有效支撑.      

Numerical analysis of a moving boundary problem in coastal hydrodynamics

A finite element model to tackle the moving boundary problem of wave run-up on moderately steep slopes is developed. The special aspects considered in this study are (1) the modification of shallow water equations to accommodate the effect of vertical accelerations and (2) the use of Lagrangian acceleration coupled with an element that adapts itself to the moving boundary closely. The pressure term in the one-dimensional momentum equation is derived using the Eulerian equation in the vertical direction. This takes care of the vertical accelerations which are significant during the motion of a wave on moderately steep slopes. The element near the boundary is allowed to change its dimension so that the fluid boundary is closely followed. Such a flexible element precludes the need for approximation of the variables with regard to the indefinite position of the boundary. This element is split into two when its dimension becomes unduly large compared to the unchanging elements. The need for such a splitting is shown by an examination of the entries in the global matrix. Results of water profile as a wave runs up a structure are given. A brief history of the work on similar problems is outlined.     

Theoretical analysis of the dynamic behavior of hysteresis elements in mechanical systems

Many machine elements in common engineering use exhibit the characteristic of “hysteresis springs”. Plain and rolling element bearings that are widely used in motion guidance of machine tools are typical examples. The study of the non-linear dynamics caused by such elements becomes imperative if we wish to achieve accurate control of such machines.

This paper outlines the properties of rate-independent hysteresis and shows that the calculation of the free response of a single-degree-of-freedom (SDOF) mass-hysteresis-spring system is amenable to an exact solution. The more important issue of forced response is not so, requiring other methods of treatment. We consider the approximate describing function method and compare its results with exact numerical simulations. Agreement is good for small excitation amplitudes, where the system approximates to a linear mass-spring-damper system, and for very large amplitudes, where some sort of mass-line is approached. Intermediate values however, show high sensitivity to amplitude variations, and no regular solution is obtained by either approach. This appears thus to be an inherent property of the system pointing to the need for developing further analysis methods.     

Nanoscale finite element models for vibrations of single-walled carbon nanotubes: atomistic versus continuum

By the atomistic and continuum finite element models, the free vibration behavior of single-walled carbon nanotubes （SWCNTs） is studied. In the atomistic finite element model, the bonds and atoms are modeled by the beam and point mass elements, respectively. The molecular mechanics is linked to structural mechanics to determine the elastic properties of the mentioned beam elements. In the continuum finite element approach, by neglecting the discrete nature of the atomic structure of the nanotubes, they are modeled with shell elements. By both models, the natural frequencies of SWCNTs are computed, and the effects of the geometrical parameters, the atomic structure, and the boundary conditions are investigated. The accuracy of the utilized methods is verified in comparison with molecular dynamic simulations. The molecular structural model leads to more reliable results, especially for lower aspect ratios. The present analysis provides valuable information about application of continuum models in the investigation of the mechanical behaviors of nanotubes.     

THE DESIGN AND ANALYSIS OF VIBRATION STRUCTURE OF VERTICAL DYNAMIC BALANCING MACHINE

I. INTRODUCTION The shape of some ?ying objects (called workpieces) is similar to that of shafts. These ?ying objectsoften have a large size. Their dispersion e?ect a?ects the ?ying control and the precision in hitting atarget. The dispersion e?ect is the statistical deviation when ?ying objects depart from their preconcerted?ying orbits. As the unevenness of mass distribution is the major cause of the dispersion e?ect, thedouble-plane vertical dynamic balancing technique for ?ying ob…     

一种改进的水下三阵元辅助惯导系统的定位算法

针对水下惯导系统定位误差随时间发散的缺陷,设计了一种基于DVL、声学定位辅助惯导的组合导航系统。声学定位系统由安装在水下潜器上的换能器阵和布放在海底的三阵元定位基阵组成,传统长基线定位水下至少要布放四个阵元,三个阵元无法直接定位求解。根据三阵元测得的距离信息和惯导系统输出的位置信息建立距离耦合模型,给出了惯导、DVL和声学定位系统的状态方程和量测方程;同时针对水下声速垂直分布特性,提出了一种等效声速的修正方法。仿真结果表明:在三阵元定位基阵作用范围内,组合导航系统能有效修正惯导位置误差,定位误差小于4 m,具有良好的适用性。     

Experimental methodology for the tappet characterization of timing system in I.C.E.

The aim of present work is the containment of the inertia forces, the stiffness components optimization and the fit tolerances of valve train in internal combustion engines (I.C.E.) 4T. The proposed methodology allows, through the development of a test machine, the evaluation of axial stiffness of tappet depending on eccentricity of the cam tappet contact, performing a functional analysis that simulate the behaviour of the system in operational condition, even if, some adjustment of tolerances of the fit between tappet and his guide, occurred. The dynamic study of the valve train, through modern computer codes, is performed by connecting lumped masses, springs and dampers that characterize each element. In numerical models the tappet is represented as constituted by the tappet and by the hydraulic element. Each of these elements is characterized by stiffness and mass. The structural rigidity of the tappet has, in fact, important effects on the dynamic behaviour of the entire valve train. The test machine makes possible the choice of the dimensional and geometrical tolerances of the fit between tappet and his guide; allows furthermore the evaluation of errors occurred during construction and integration phase. In addition, the test machine is also suitable for reverse engineering applications, makes it possible to automatically draw the cam profile in polar coordinates.     

Computation of the pressure on the surface of a fuselage with a wing in an ideal fluid

The flows around complex three-dimensional bodies by an ideal fluid are computed by methods [1–3] using approximation of the surface by a set of plane elements. A layer of surface singularities, whose intensity is found by solving a system of linear algebraic equations of very high order, is distributed continuously over each element. Evaluation of the system coefficients and its solution require significant machine time expenditures on powerful electronic computers. If in the method of [2] the total system of equations is separated successfully into several subsystems by simplifications and an approximate solution of the problem is obtained more rapidly than by the method in [1], then the same author practically used the method in [1] to design specific fuselages in [3]. A method [4] developed for a fuselage is expanded in this paper to design a wing-fuselage combination. This method turns out to be less tedious, without being inferior in accuracy, by being different from the method in [1] in the means of solving the fundamental integral equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–115, May–June, 1977.     

新型立式动平衡机的设计与分析

在分析传统的双面立式动平衡机的振动结构的基础上，设计了一种新的立式动平衡机的振动结构(即摆架)，采用新摆架结构的动平衡机可以实现静不平衡和偶不平衡的有效分离。文中详细分析了新摆架结构的特点，并在ANSYS7．0平台上对摆架进行了有限元模态分析和谐响应分析。由于新的摆架结构采用多块簧板，在多个部位刚度匹配，且采用多个传感器安装于测量点处，这样新的动平衡机可以直接测量静不平衡和偶不平衡，且相互影响极小。实践证明：采用新摆架结构的动平衡机特别适合于飞行物体的惯性测量，且克服了传统的双面立式动平衡机长期以来平面分离不佳的缺点，具有很高的推广价值。     

考虑变形耦合的几何非线性空间梁单元

以"精确几何模型梁单元"为代表的很多几何非线性梁单元,在构造过程中分别对描述截面转动的转角和描述截面形心位置的位移进行了独立插值,由此引起了诸如运动学描述冗余和剪切闭锁等困难。其根本原因在于单元形函数没能体现细长梁中的变形耦合关系。本文对这类传统单元进行了改造,通过深入研究单元变形之间的内在联系,提出了一种变形场完全满足Bernoulli梁变形耦合关系的新单元,避免了构造过程中对转动矢量的插值,并通过数值算例检验了单元的有效性。     

Geometrical shapes of waterflooded zones in an inhomogeneous reservoir

The geometrical shapes of waterflooded zones are determined for an individual injection well in a powerful reservoir whose structure is vertically inhomogeneous. It is found that with decrease in the permeability from the top of the reservoir to its bottom the side spreading of water injected in the neighborhood of the bottom also decreases appreciably. In the cases of the power and exponential vertical variations in the permeability and voidage analytic expressions for the pressure distribution and the geometrical shapes of volumes are obtained.     

基于纤维单元的钢筋混凝土桥墩地震损伤评价

基于纤维模型的粱柱单元特别是柔度法单元有着较高的计算精度,已被较多地用于评价钢筋混凝土桥墩的位移（廷性）能力。随着基于性能／位移抗震设计理论的发展,相继提出了残余位移、极限曲率及曲率延性系数、纵筋和混凝土的最大应变、纵筋低周疲劳损伤等桥墩地震损伤量化指标。选用刚度法和柔度法纤维单元,考虑材料非线性、几何非线性和结点锚固...     

Oblique water entry of a cone by a fully three-dimensional nonlinear method

The hydrodynamic problem of a cone entering the water surface obliquely has been analyzed by the three-dimensional (3-D) incompressible velocity potential theory with the fully nonlinear boundary conditions on the moving free surface and body surface boundary. The time stepping method is used in the stretched coordinate system defined as the ratio of the physical system to the distance that the cone has travelled into water. The boundary element method is used to solve the potential at each time step. Both triangular element mesh and quadrilateral element mesh have been used. Discretisation of the body surface and the free surface is applied regularly during the simulation to account for their change and deformation, and data from the old mesh is transferred into the new one through interpolation. Both the dynamic and kinematic free surface boundary conditions are satisfied through the Eulerian form. In particular the free surface elevation and potential variation are traced at a given azimuth of the cylindrical coordinate system, in the direction parallel to the body surface or perpendicular to the free surface to avoid multi-valued function. Detailed convergence study with respect to time step and element size has been undertaken and high accuracy has been achieved. Results for the cone in vertical entry are compared with those obtained from the 2-D axisymmetric method and good agreement is found. Simulations are made for cones of various deadrise angles and different oblique entries and detailed results are provided.     

Dynamic behavior of long-span box girder bridges subjected to moving loads: Numerical analysis and experimental verification

In this study we conducted three-dimensional dynamic analyses of long-span box girder bridges subjected to moving loads, using four-node Lagrangian and Hermite finite elements. In finite element formulation, a 6 × 6 transformation matrix is derived to transform the system element matrices before assembly. The usual 5 degrees of freedom per node are appended with an additional drilling degree of freedom in order to fit the transformation. The numerical results show good agreement with the experimental data from an existing two-span prestressed concrete box girder bridge under travelling vehicles. Parametric studies are focused on the various effects of moving loads on the dynamic behavior for different locations on the cross-section of box girder bridges.     

Numerical modelling of structure vibrations under inertial moving load

Inertial loading of structures by mass travelling with near-critical velocity has been intensively debated. In the literature a moving mass is replaced by an equivalent force or an oscillator that is in permanent contact with the structure. A direct mass matrix modification method frequently implemented in the finite element approach gives reasonable results only in the range of relatively low velocities and for low mass value if compared with the mass of a structure. However, existing solutions are incorrect and are not implemented in commercial computer codes. In this paper we present the space–time finite element approach to the problem. The interaction of the moving mass/supporting structure is described in a local coordinate system of the space-time finite element domain. Resulting characteristic matrices include inertia, Coriolis and centrifugal forces. Simple modification of matrices in the discrete equations of motion allows us to gain accuracy in a wide range of velocity, up to the over-critical speed. Numerical examples of string and beam vibrations prove the simplicity and efficiency of the method.     

Adaptive non‐conformal mesh refinement and extended finite element method for viscous flow inside complex moving geometries

Two different techniques to analyze non‐Newtonian viscous flow in complex geometries with internal moving parts and narrow gaps are compared. The first technique is a non‐conforming mesh refinement approach based on the fictitious domain method (FDM), and the second one is the extended finite element method (XFEM). The refinement technique uses one fixed reference mesh, and to impose continuity across non‐conforming regions, constraints using Lagrangian multipliers are used. The size of elements locally in the high shear rate regions is reduced to increase accuracy. FDM is shown to have limitations; therefore, XFEM is applied to decouple the fluid from the internal moving rigid bodies. In XFEM, the discontinuous field variables are captured by using virtual degrees of freedom that serve as enrichment and by applying special integration over the intersected elements. The accuracy of the two methods is demonstrated by direct comparison with results of a boundary‐fitted mesh applied to a two‐dimensional cross section of a twin‐screw extruder. Compared with non‐conforming FDM, XFEM shows a considerable improvement in accuracy around the rigid body, especially in the narrow gap regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Optimization of the response times of the feed kinematical linkages of the numerical control machine tools in order to minimize the path error

The main topic refers to enhancing the machining accuracy by minimizing the path error while machining through linear, circular or other interpolations on CNC machines for milling, boring and milling, grinding or turning. Currently the path error minimizing process while machining on CNC machine tools is a little explored, because of the absence of the dependence relations between the path error and the parameters of the kinematical axes taking part to the interpolation of the work piece contour. This work establishes and analyzes the relations of dependence between the path error and the response times of the kinematical axes that take part to the linear and circular interpolation. The theoretical results are experimentally verified by performing several machining operations through circular interpolation, where the moving speed along the contour and the response times of the kinematical axes are modified. The results of this research are intended for the machine tool designers and those who use such machines.     

基于Bézier提取的三维等几何分析

等几何分析(IGA)将非均匀有理B样条(NURBS)函数作为有限元形函数,具有几何精确、高阶连续和精度高等优点。与常规有限元法C0连续的形函数不同,高阶IGA基函数不是定义在一个单元上,而是跨越由几个单元组成的参数空间,因而编程复杂且无法嵌入现有的有限元法计算框架及相应算法。本文建立了基于Bézier提取的三维IGA,将NURBS函数分解成伯恩斯坦多项式的线性组合,从而实现把NURBS单元分解为C0连续Bézier单元,这些单元与Lagrange单元相似,使IGA的实现和常规有限元一样,以便将IGA分析嵌入现有的有限元软件中。两个三维算例结果表明,基于Bézier提取的IGA和传统IGA的收敛性和精度相同。     

Investigation of stress-strain behavior of single walled carbon nanotube/rubber composites by a multi-scale finite element method

The excellent properties of carbon nanotubes have generated technological interests in the development of nanotube/rubber composites. This paper describes a finite element formulation that is appropriate for the numerical prediction of the mechanical behavior of rubber-like materials which are reinforced with single walled carbon nanotubes. The considered composite material consists of continuous aligned single walled carbon nanotubes which are uniformly distributed within the rubber material. It is assumed that the carbon nanotubes are imperfectly bonded with the matrix. Based on the micromechanical theory, the mechanical behavior of the composite may be predicted by utilizing a representative volume element. Within the representative volume element, the reinforcement is modeled according to its atomistic microstructure. Therefore, non-linear spring-based line elements are employed to simulate the discrete geometrical structure and behavior of the single-walled carbon nanotube. On the other hand, the matrix is modeled as a continuum medium by utilizing solid elements. In order to describe its behavior an appropriate constitutive material model is adopted. Finally, the interfacial region is simulated via the use of special joint elements of variable stiffness which interconnect the two materials in a discrete manner. Using the proposed multi-scale model, the stress-strain behavior for various values of reinforcement volume fraction and interfacial stiffness is extracted. The influence of the single walled carbon nanotube addition within the rubber is clearly illustrated and discussed.     

直接基于单元平衡的梁非线性分析方法

以单元力平衡条件为基础建立的梁单元在近年来受到广泛关注,与传统的以单元位移相容条件为基础建立的梁单元相比,基于平衡的梁单元可获得较高的精度。针对基于平衡的梁单元,给出了一种几何非线性分析方法,使用该方法进行几何非线性分析时不需进行单元内部位移场的构造,避免了使用现有位移场构造方法时可能引起的振荡问题,同时降低了有限元程...