• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.113-140
• /
• 2004

A mixed eight-node hexahedral element formulated via the Hu-Washizu principle as well as the field extrapolation technique is presented. The mixed element with only three translational degrees of freedom at each node can provide extremely accurate and reliable performance for popular benchmark problems such as spacial beams, plates, shells as well as general three-dimensional elasticity problems. Numerical calculations also show that when extremely skewed and coarse meshes and nearly incompressible materials are used, the proposed mixed element can still possess excellent behaviour. The mixed formulation starts with introduction of a parallelepiped domain associated with the given general eight-node hexahedral element. Then, the assumed strain field at the nodal level is constructed via the Hu-Washizu variational principle for that associated parallelepiped domain. Finally, the assumed strain field at the nodal level of the given hexahedral element is established by using the field extrapolation technique, and then by using the trilinear shape functions the assumed strain field of the whole element domain is obtained. All matrices involved in establishing the element stiffness matrix can be evaluated analytically and expressed explicitly; however, a 24 by 24 matrix has to be inverted to construct the displacement extrapolation matrix. The proposed hexahedral element satisfies the patch test as long as the element with a shape of parallelepiped.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.2
• /
• pp.365-371
• /
• 1998

In the present study, we attempted to add the incompatible functions as additional variable terms to the conventional u-p formulation. It is derived from the four-field generalized variational principle that encompasses velocity, pressure, velocity strains and stress fields as independent interpolated variables. As a severe test of the present formulation, we have investigated the driven cavity with the corner velocity singularity like leaky lid. Through the test, the present element performs very well without unstable oscillation of pressure distribution.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.645-653
• /
• 1996

In order to avoid pathological mesh dependency in finite element modelling of strain localization, an isotropic elasto-plastic model with a yield function depending on the Laplacian of the equivalent plastic strain is implemented in a 4-node quadrilateral finite element with one integration point based on a mixed formulation derived from Hu-Washizu principle. The evaluation of the Laplacian is based on a least square polynomial approximation of the equivalent plastic strain around each integration point. This non local approach allows to satisfy exactly the consistency condition at each integration point. Some examples are treated to illustrate the effectiveness of the method.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.375-396
• /
• 2009

In the present work, the elasto-viscoplastic behavior, interactions between grains, and the texture evolution in polycrystalline materials subjected to finite deformations are modeled using a multiscale analysis procedure within a finite element framework. Computational homogenization is used to relate the grain (meso) scale to the macroscale. Specifically, a polycrystal is modeled by a material representative volume element (RVE) consisting of an aggregate of grains, and a periodic distribution of such unit cells is considered to describe material behavior locally on the macroscale. The elastic behavior is defined by a hyperelastic potential, and the viscoplastic response is modeled by a simple power law complemented by a work hardening equation. The finite element framework is based on a Lagrangian formulation, where a kinematic split of the deformation gradient into volume preserving and volumetric parts together with a three-field form of the Hu-Washizu variational principle is adopted to create a stable finite element method. Examples involving simple deformations of an aluminum alloy are modeled to predict inhomogeneous fields on the grain scale, and the macroscopic effective stress-strain curve and texture evolution are compared to those obtained using both upper and lower bound models.

• Structural Engineering and Mechanics
• /
• v.26 no.6
• /
• pp.725-739
• /
• 2007

In recent years there are many plate bending elements that emerged for solving both thin and thick plates. The main features of these elements are that they are based on mix formulation interpolation with discrete collocation constraints. These elements passed the patch test for mix formulation and performed well for linear analysis of thin and thick plates. In this paper a member of this family of elements, namely, the Discrete Reissner-Mindlin (DRM) is further extended and developed to analyze both thin and thick plates with geometric nonlinearity. The Von K$\acute{a}$rm$\acute{a}$n's large displacement plate theory based on Lagrangian coordinate system is used. The Hu-Washizu variational principle is employed to formulate the stiffness matrix of the geometrically Nonlinear Discrete Reissner-Mindlin (NDRM). An iterative-incremental procedure is implemented to solve the nonlinear equations. The element is then tested for plates with simply supported and clamped edges under uniformly distributed transverse loads. The results obtained using the geometrically NDRM element is then compared with the results of available analytical solutions. It has been observed that the NDRM results agreed well with the analytical solutions results. Therefore, it is concluded that the NDRM element is both reliable and efficient in analyzing thin and thick plates with geometric non-linearity.