共查询到20条相似文献,搜索用时 0 毫秒
1.
2.
The goal of this paper is to present a flexible multibody formulation for Euler-Bernoulli beams involving large displacements. This method is based on a discretisation of internal and kinetic energies. The beam is represented by its line of centroids and each section is oriented by a frame defined by three Euler angles. We apply a finite element formulation to describe the evolution of these angles along the neutral fibre and describe the internal energy. The kinetic energy is approximated as the one of two rigid bars tangent to the neutral fibre at the ends of the beam. We derive the equations of motion from a Lagrange formulation. These equations are solved using the Newmark method or/and the Newton-Raphson technique. We solve some very classic problems taken from the literature as the curved beam presented by Simo [Simo, J. C., ‘A three-dimensional finite-strain rod model. the three-dimensional dynamic problem. Part I’, Comput. Meths. Appl. Mech. Engrg.
49, 1985, 55–70; Simo, J. C. and Vu-Quoc, L., ‘A three-dimensional finite-strain rod model, Part II: Computationals aspects’, Comput. Meths. Appl. Mech. Engrg.
58, 1988, 79–116] and Lee [Lee, Kisu, ‘Analysis of large displacements and large rotations of three-dimensional beams by using small strains and unit vectors’, Commun. Numer. Meth. Engrg.
13, 1997, 987–997] or the rotational rod presented by Avello [Avello, A., Garcia de Jalon, J., and Bayo, E., ‘Dynamics of flexible multibody systems using cartesian co-ordinates and large displacement theory’, Int. J. Num. Methods in Engineering
32, 1991, 1543–1563] and Simo [Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part I’ Jour. of Appl. Mech.
53, 1986, 849–854; Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part II’, Jour. of Appl. Mech.
53, 1986, 855–863]. 相似文献
3.
A wide variety of mechanical and structural multibody systems consist ofvery flexible components subject to kinematic constraints. The widelyused floating frame of reference formulation that employs linear modelsto describe the local deformation leads to a highly nonlinear expressionfor the inertia forces and can be applied to only small deformationproblems. This paper is concerned with the formulation and computerimplementation of spatial joint constraints and forces using the largedeformation absolute nodal coordinate formulation. Unlike the floatingframe of reference formulation that employs a mixed set of absolutereference and local elastic coordinates, in the absolute nodalcoordinate formulation, global displacement and slope coordinates areused. The nonlinear kinematic constraint equations and generalized forceexpressions are expressed in terms of the absolute global displacementsand slopes. In particular, a new formulation for the sliding jointbetween two very flexible bodies is developed. A surface parameter isintroduced as an additional new variable in order to facilitate theformulation of this sliding joint. The constraint and force expressionsdeveloped in this paper are also expressed in terms of generalizedCholesky coordinates that lead to an identity inertia matrix. Severalexamples are presented in order to demonstrate the use of theformulations developed in the paper. 相似文献
4.
In the general theory of continuum mechanics, the state of rotation and deformation of material points can be uniquely defined from the displacement field by using the nine independent components of the displacement gradients. For this reason, the use of the absolute rotation parameters as nodal coordinates, without relating them to the displacement gradients, leads to coordinate redundancy that leads to numerical and fundamental problems in many existing large rotation finite element formulations. Because of this fundamental problem, special measures that require modifications of the numerical integration methods were proposed in the literature in order to satisfy the principle of work and energy. As demonstrated in this paper, no such measures need to be taken when the finite element absolute nodal coordinate formulation is used since the principle of work and energy are automatically satisfied. This formulation does not suffer from the problem of coordinate redundancy and ensures the continuity of stresses and strains at the nodal points. In this study, the use of the implicit integration methods with the consistent Lagrangian elasto-plastic tangent moduli is examined when the absolute nodal coordinate formulation is used. The performance of different numerical integration methods in the dynamic analysis of large elasto-plastic deformation problems is investigated. It is shown that all these methods, in the case of convergence, yield a solution that satisfies the principle of work and energy without the need of taking any special measures. Semi-implicit integration methods, however, can lead to numerical difficulties in the case of very stiff problems due to the linearization made in these methods in order to avoid the iterative Newton--Raphson procedure. It is also demonstrated that the use of the consistent Lagrangian-plastic tangent moduli derived in this investigation using the absolute nodal coordinate formulation leads to better convergence of the iterative Newton--Raphson procedure used in the implicit integration methods. 相似文献
5.
6.
7.
Dynamic analysis of a rotating rigid-flexible coupled smart structure with large deformations 总被引:2,自引:0,他引:2
Based on Hamilton's principle,a new kind of fully coupled nonlinear dynamic model for a rotating rigid-flexible smart structure with a tip mass is proposed.The geometrically nonlinear effects of the axial,transverse displacement and rotation angle are considered by means of the first-order approximation coupling(FOAC)model theory, in which large deformations and the centrifugal stiffening effects are considered.Three kinds of systems are established respectively,which are a structure without piezoelectric layer,with piezoelectric layer in open circuit and closed circuit.Several simulations based on simplified models are presented to show the differences in characteristics between structures with and without the tip mass,between smart beams in closed and open circuit, and between the centrifugal effects in high speed rotating state or not.The last simulation calculates the dynamic response of the structure subjected to external electrical loading. 相似文献
8.
9.
有限元模型比知识的获取和形式化是有限元模型化专家系统建立的核心。本文以信息论为工具,在总结大量工程经验和理论研究成果的基础上,建立了工程结构模型化知识系统,它分为两个层次,一是复杂结构模型化,二是简单结构模型化。后者是前者的基础,包括力学模型化、有限元模型化和网络生成三部分。这些知识均已溶人了作者开发的专家系统之中。 相似文献
10.
三维液体非线性晃动动力学特性的数值模拟 总被引:7,自引:0,他引:7
主要讨论圆筒形贮腔中三维液体非线性晃动问题,将任意的拉格朗日-欧拉(即ArbitraryLagrangian-Eulerian,简称ALE)运动学描述引入到Navies-Stokes方程中,在时间域上采用一种速度和压力的分步计算格式进行时间离散;在空间域上利用Galerkin加权余量法对系统方程进行数值离散;得到了数值计算粘性不可压液体非线性晃动的ALE分步有限元法的计算格式,推导了三维液体自由液面上结点法向矢量的数值计算方法,模拟了圆筒形贮腔(包括带圆环形隔板的圆筒形贮腔)中三维液体的非线性晃动;并得到了一些重要的非线必不知所云 性,通过数值模拟结果与实验结果的比较,证明实了本文方法的可靠性与有效性。 相似文献
11.
火灾试验研究表明,火灾下无粘结预应力混凝土受弯构件的变形非常大,其中转动也很大,尤其是单面受火的板,是典型的材料非线性、几何非线性问题。如何对其进行精确、高效地求解是值得关注的问题。而基于S-R分解原理的更新拖带坐标有限元法有诸多优点:有利于跟踪变形物体中各点的变形;保证单元的质量守恒;在有限元增量法求解时,还可以避免对坐标的修正;而且将转动作为一个独立的自由度,提高了求解效率。在升温过程中,虽然预应力钢筋的应力处处相等,但其各点的温度不同,导致各点的温度应变、蠕变、塑性应变均不同。预应力筋必然产生滑移。本文采用该方法对火灾下无粘结预应力钢筋混凝土板进行编程分析。通过实际算例验证该算法的可靠性,该方法求解效率高,精度也比较好。 相似文献
12.
Finite Element Formulation of Viscoelastic Constitutive Equations Using Fractional Time Derivatives 总被引:9,自引:0,他引:9
Fractional time derivatives are used to deduce a generalization ofviscoelastic constitutive equations of differential operator type. Theseso-called fractional constitutive equations result in improvedcurve-fitting properties, especially when experimental data from longtime intervals or spanning several frequency decades need to be fitted.Compared to integer-order time derivative concepts less parameters arerequired. In addition, fractional constitutive equations lead to causalbehavior and the concept of fractional derivatives can be physicallyjustified providing a foundation of fractional constitutive equations.First, three-dimensional fractional constitutive equations based onthe Grünwaldian formulation are derived and their implementationinto an elastic FE code is demonstrated. Then, parameter identificationsfor the fractional 3-parameter model in the time domain as well as inthe frequency domain are carried out and compared to integer-orderderivative constitutive equations. As a result the improved performanceof fractional constitutive equations becomes obvious. Finally, theidentified material model is used to perform an FE time steppinganalysis of a viscoelastic structure. 相似文献
13.
等截面梁有限变形的传递函数解法 总被引:2,自引:0,他引:2
本文应用传递函数方法对等截面梁的有限变形进行了分析,对于等截面梁的有限变形问题,该方法从变分方程出发把问题表述为状态空间的形式,然后利用Gauss积分对轴力进行加速迭代求解,不需要进行增量迭代即可取得具有良好计算精度的数值结果,对简单受力的等截面梁情况该解可以看作是所讨论问题的精确解。对于受力比较复杂或者阶梯变截面梁情况,为减少运算量,可以和有限元法类似,采用多个单元进行拼接,从而得到问题的解。数值算例表明,本方法具有半解析、精度高、收敛快等特点。 相似文献
14.
15.
16.
In multibody system dynamics, the absolute nodal coordinate formulation(ANCF)uses power functions as interpolating polynomials to describe the displacement field. It can get accurate results for flexible bodies that undergo large deformation and large rotation. However, the power functions are irrational representation which cannot describe the complex shapes precisely, especially for circular and conic sections. Different from the ANCF representation,the rational absolute nodal coordinate formulation(RANCF) utilizes rational basis functions to describe geometric shapes, which allows the accurate representation of complicated displacement and deformation in dynamics modeling. In this paper, the relationships between the rational surface and volume and the RANCF finite element are provided, and the generalized transformation matrices are established correspondingly. Using these transformation matrices, a new four-node three-dimensional RANCF plate element and a new eight-node three-dimensional RANCF solid element are proposed based on the RANCF. Numerical examples are given to demonstrate the applicability of the proposed elements. It is shown that the proposed elements can depict the geometric characteristics and structure configurations precisely, and lead to better convergence in comparison with the ANCF finite elements for the dynamic analysis of flexible bodies. 相似文献
17.
18.
19.
SHINICHIRO MIURA KAZUHIKO KAKUDA NOIJUYOSHI TOSAFCA 《International Journal of Computational Fluid Dynamics》2013,27(3-4):291-299
In previous papers, we proposed finite element schemes based on the Petrov-Galerkin weak formulation using exponential weighting functions for solving accurately, and in a stable manner, the flow field of an incompressible viscous fluid. In this paper, we present the Petrov-Galerkin finite element scheme for turbulent flow fields based on large eddy simulation using the standard Smagorinsky model with the Van Driest damping function. The filtered incompressible Navier-Stokes equations are numerically integrated in time by using a fractional step strategy with second-order accurate Adams-Bashforth explicit differencing for both convection an diffusion terms. In order to evaluate more accurately a mass matrix, the well-known multi-pass algorithm was also adopted in this study. Numerical results obtained are compared through flow around a rectangular cylinder at Re = 22,000 with the experimental data and other existing numerical data. 相似文献
20.
ABSTRACT A new procedure for analysis of tracked-vehicle suspension systems under a track length constraint, based on a frictional contact problem formulation, is presented. Impenetrability, inextensibility, unilateral behavior of track, and friction conditions are expressed in complementarity form. The resulting problem is then amenable to efficient numerical algorithms, such as the modified simplex method. This approach is demonstrated through numerical applications to a model of a tracked-vehicle suspension system. 相似文献