A new finite element formulation of three-dimensional beam theory based on interpolation of curvature

A new finite element formulation of the `kinematically exact finite-strain beam theory' is presented. The finite element formulation employs the generalized virtual work in which the main role is played by the pseudo-curvature vector. The solution of the governing equations is found by using a combined Galerkin-collocation algorithm

Locking-free two-layer Timoshenko beam element with interlayer slip

A new locking-free strain-based finite element formulation for the numerical treatment of linear static analysis of two-layer planar composite beams with interlayer slip is proposed. In this formulation, the modified principle of virtual work is introduced as a basis for the finite element discretization. The linear kinematic equations are included into the principle by the procedure, similar to that of Lagrangian multipliers. A strain field vector remains the only unknown function to be interpolated in the finite element implementation of the principle. In contrast with some of the displacement-based and mixed finite element formulations of the composite beams with interlayer slip, the present formulation is completely locking-free. Hence, there are no shear and slip locking, poor convergence and stress oscillations in these finite elements. The generalization of the composite beam theory with the consideration of the Timoshenko beam theory for the individual component of a composite beam represents a substantial contribution in the field of analysis of non-slender composite beams with an interlayer slip. An extension of the present formulation to the non-linear material problems is straightforward. As only a few finite elements are needed to describe a composite beam with great precision, the new finite element formulations is perfectly suited for practical calculations. (c) 2007 Elsevier B.V. All rights reserved

A finite membrane element formulation for surfactants

Surfactants play an important role in various physiological and biomechanical applications. An example is the respiratory system, where pulmonary surfactants facilitate the breathing and reduce the possibility of airway blocking by lowering the surface tension when the lung volume decreases during exhalation. This function is due to the dynamic surface tension of pulmonary surfactants, which depends on the concentration of surfactants spread on the liquid layer lining the interior surface of the airways and alveoli. Here, a finite membrane element formulation for liquids is introduced that allows for the dynamics of concentration-dependent surface tension, as is the particular case for pulmonary surfactants. A straightforward approach is suggested to model the contact line between liquid drops/menisci and planar solid substrates, which allows the presented framework to be easily used for drop shape analysis. It is further shown how line tension can be taken into account. Following an isogeometric approach, NURBS-based finite elements are used for the discretization of the membrane surface. The capabilities of the presented computational model is demonstrated by different numerical examples - such as the simulation of liquid films, constrained and unconstrained sessile drops, pendant drops and liquid bridges - and the results are compared with experimental data.Comment: Some typos are removed. Eqs. 13 and 105 are modified. Eqs. 64 and 73 are added; thus, the rest of equations is renumbered. All the numerical experiments are repeated. The example of Sec. 6.3 is slightly modifie

Finite element formulation for linear thermoviscoelastic materials

Report presents the finite difference equations in time and finite element matrix equations in space for general linear thermovisoelastic problems. The equations are derived for a general three-dimensional body but are applicable to one- and two-dimensional configurations with minor changes

The strain-based beam finite elements in multibody dynamics

We present a strain-based finite-element formulation for the dynamic analysis of flexible elastic planar multibody systems, composed of planar beams. We consider finite displacements, rotations and strains. The discrete dynamic equations of motion are obtained by the collocation method. The strains are the basic interpolated variables, which makes the formulation different from other formulations. The further speciality of the formulation is the strong satisfaction of the cross-sectional constitutive conditions at collocation points. In order to avoid the nested integrations, a special algorithm for the numerical integration over the length of the finite element is proposed. The midpoint scheme is used for the time integration. The performance of the formulation is illustrated via numerical examples, including a stiff multibody system. (c) 2007 Elsevier Ltd. All rights reserved

Non-negative mixed finite element formulations for a tensorial diffusion equation

We consider the tensorial diffusion equation, and address the discrete maximum-minimum principle of mixed finite element formulations. In particular, we address non-negative solutions (which is a special case of the maximum-minimum principle) of mixed finite element formulations. The discrete maximum-minimum principle is the discrete version of the maximum-minimum principle. In this paper we present two non-negative mixed finite element formulations for tensorial diffusion equations based on constrained optimization techniques (in particular, quadratic programming). These proposed mixed formulations produce non-negative numerical solutions on arbitrary meshes for low-order (i.e., linear, bilinear and trilinear) finite elements. The first formulation is based on the Raviart-Thomas spaces, and is obtained by adding a non-negative constraint to the variational statement of the Raviart-Thomas formulation. The second non-negative formulation based on the variational multiscale formulation. For the former formulation we comment on the affect of adding the non-negative constraint on the local mass balance property of the Raviart-Thomas formulation. We also study the performance of the active set strategy for solving the resulting constrained optimization problems. The overall performance of the proposed formulation is illustrated on three canonical test problems.Comment: 40 pages using amsart style file, and 15 figure

A micromechanics-enhanced finite element formulation for modelling heterogeneous materials

In the analysis of composite materials with heterogeneous microstructures, full resolution of the heterogeneities using classical numerical approaches can be computationally prohibitive. This paper presents a micromechanics-enhanced finite element formulation that accurately captures the mechanical behaviour of heterogeneous materials in a computationally efficient manner. The strategy exploits analytical solutions derived by Eshelby for ellipsoidal inclusions in order to determine the mechanical perturbation fields as a result of the underlying heterogeneities. Approximation functions for these perturbation fields are then incorporated into a finite element formulation to augment those of the macroscopic fields. A significant feature of this approach is that the finite element mesh does not explicitly resolve the heterogeneities and that no additional degrees of freedom are introduced. In this paper, hybrid-Trefftz stress finite elements are utilised and performance of the proposed formulation is demonstrated with numerical examples. The method is restricted here to elastic particulate composites with ellipsoidal inclusions but it has been designed to be extensible to a wider class of materials comprising arbitrary shaped inclusions.Comment: 28 pages, 12 figures, 2 table