• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Computers and Concrete
• /
• v.4 no.2
• /
• pp.135-156
• /
• 2007

Three 3D nonlinear finite-element models are developed to study the behavior of concrete beams and plates with and without external reinforcement by fibre-reinforced plastic (FRP). All three models are formulated based upon the 3D theory of elasticity. The stress model is modified from the element developed by Ramtekkar, et al. (2002) to incorporate material nonlinearity in the formulation. Both transverse stress and displacement components are used as nodal degrees-of-freedom to ensure the continuity of both stress and displacement components between the elements. The displacement model uses only displacement components as nodal degrees-of-freedom. The transition model has both stress and displacement components as nodal degrees-of-freedom on one surface, and only displacement components as nodal degrees-of-freedom on the opposite surface. The transition model serves as a connector between the stress and the displacement models. The developed models are validated by comparing the results of the analyses with an existing experimental result. Parametric studies of the effects of the externally reinforced FRP on the load capacity of reinforced concrete (RC) beams and concrete plates are performed to demonstrate the practicality and the efficiency of the proposed models.

• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• Structural Engineering and Mechanics
• /
• v.35 no.3
• /
• pp.315-333
• /
• 2010

The structural behaviors of circular concrete filled steel tube (CFT) structures are investigated by nonlinear finite element method. An efficient three-dimensional (3D) degenerated beam element is adopted. Based on those previous studies, a modified stress-strain relationship for confined concrete which introduces the influence of eccentricity on confining stress is presented. Updated Lagrange formulation is used to consider the geometrical nonlinearity induced by large deformation effect. The nonlinear behaviors of CFT structures are investigated, and the accuracy of the proposed constitutive model for confined concrete is mainly concerned. The results demonstrate that the confining effect in CFT elements subjected to combining action of axial force and bending moment is far sophisticated than that in axial loaded columns, and an appropriate evaluation about this effect may be important for nonlinear numerical simulation of CFT structures.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.38-42
• /
• 2000

In this study, the low velocity impact behavior of the composite laminates has been described by using 3 dimensional nonlinear finite elements. To describe the geometric nonlinearity due to large deformation, the dynamic contact problem is formulated using the exterior penalty finite element method on the base of Total Lagrangian formulation. The incremental decomposition is introduced, and the converged solution is attained by Newton-Raphson Method. The Newmark's constant-acceleration time integration algorithm is used. To make verification of the finite element program developed in this study, the solution of the nonlinear static problem with occurrence of large deformation is compared with ABAQUS, and the solution of the static contact problem with indentation is compared with the Hertz solution. And, the solution of low velocity impact problem for isotropic material is verificated by comparison with that of LS-DYNA3D. Finally the contact force of impact response from the nonlinear analysis are compared with those from the linear analysis.

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.33-42
• /
• 2017

We present in this paper a finite element formulation for nonlinear torsional analysis of 3D beams with arbitrary composite cross-sections. Since the proposed formulation employs a continuum mechanics based beam element with kinematics enriched by the extended St. Venant solutions, it can precisely account higher order warping effect and its 3D couplings. We propose a numerical procedure to calculate the extended St. Venant equation and the twisting center of an arbitrary composite cross-section simultaneously. The accuracy and efficiency of the proposed formulation are thoroughly investigated through representative numerical examples.

• Proceedings of the KIEE Conference
• /
• 2002.11d
• /
• pp.124-127
• /
• 2002

This paper presents a 3D shape optimization algorithm for electromagnetic devices using the design sensitivity analysis with finite element method. The structural deformation analysis based on the deformation theory of the elastic body under stress is used for mesh renewing. The design sensitivity and adjoint variable formulae are derived for the 3D nonlinear finite element method with edge element. The proposed algorithm is applied to the shape optimization of 3D electromagnet to get a uniform flux density at the air gap.

• Bulletin of the Korean Mathematical Society
• /
• v.38 no.2
• /
• pp.237-260
• /
• 2001

The approximation properties for $L^2$-projection, Raviart-Thomas projection, and inverse inequality have been derived in 3 dimensional space. h-p-mixed finite element methods for strongly nonlinear second order elliptic problems are proposed and analyzed in 3D. Solvability and convergence of the linearized problem have been shown through duality argument and fixed point argument. The analysis is carried out in detail using Raviart-Thomas-Nedelec spaces as an example.

• Geomechanics and Engineering
• /
• v.5 no.3
• /
• pp.223-240
• /
• 2013

In the present study an attempt was made to predict the complex nonlinear parameters of the soil-pile system subjected to the vertical vibration of rotating machines. A three dimensional (3D) finite element (FE) model was developed to predict the nonlinear dynamic response of full-scale pile foundation in a layered soil medium using ABAQUS/CAE. The frequency amplitude responses for different eccentric moments obtained from the FE analysis were compared with the vertical vibration test results of the full-scale single pile. It was found that the predicted resonant frequency and amplitude of pile obtained from 3D FE analysis were within a reasonable range of the vertical vibration test results. The variation of the soil-pile separation lengths were determined using FE analysis for different eccentric moments. The Novak's continuum approach was also used to predict the nonlinear behaviour of soil-pile system. The continuum approach was found to be useful for the prediction of the nonlinear frequency-amplitude response of full-scale pile after introducing the proper boundary zone parameters and soil-pile separation lengths.