• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.115-131
• /
• 2012

Paper presents an improved solution algorithm based on Finite Element Method to analyse piled raft foundation. Piles are modelled as beam elements with soil springs. Finite element analysis of raft is based on the classical theory of thick plates resting on Winkler foundation that accounts for the transverse shear deformation of the plate. Four node, isoparametric rectangular elements with three degrees of freedom per node are considered in the development of finite element formulation. Independent bilinear shape functions are assumed for displacement and rotational degrees of freedom. Effect of raft thickness, soil modulus and load pattern on the response is considered. Significant improvement in the settlements and moments in the raft is observed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.2
• /
• pp.11-22
• /
• 1989

The finite element utilizing isoparametric plate element is applied for the elastic analysis of straight box girder bridges. A continuous box girder is analyzed as an example to verify the validity and accuracy of this method. It is indicated that the deflections and longitudinal stresses obtained by this method agree well with those from the orther methods. This theory may, therefore, be directly used as an efficient tool to analyze and design box girder bridges subjected to arbitrary loading and boundary conditions.

• Steel and Composite Structures
• /
• v.15 no.2
• /
• pp.175-202
• /
• 2013

The super convergent laminated composite beam element is newly derived for the lateral stability analysis. For this, a theoretical model of the laminated composite beams is developed based on the first-order shear deformation beam theory. The present laminated beam takes into account the transverse shear and the restrained warping induced shear deformation. The second-order coupling torque resulting from the geometric nonlinearity is rigorously derived. From the principle of minimum total potential energy, the stability equations and force-displacement relationships are derived and the explicit expressions for the displacement parameters are presented by applying the power series expansions of displacement components to simultaneous ordinary differential equations. Finally, the member stiffness matrix is determined using the force-displacement relationships. In order to show accuracy and superiority of the beam element developed by this study, the critical lateral buckling moments for bisymmetric and monosymmetric I-beams are presented and compared with other results available in the literature, the isoparametric beam elements, and shell elements from ABAQUS.

• Structural Engineering and Mechanics
• /
• v.5 no.4
• /
• pp.471-490
• /
• 1997

In this paper, the classical Irons' patch tests which have been generally accepted for the convergence proof of a finite element are performed for Mindlin plate bending elements with a special emphasis on the nonconforming elements. The elements considered are 4-node and 8-node quadrilateral isoparametric elements which have been dominantly used for the analyses of plate bending problems. It was recognized from the patch tests that some nonconforming Mindlin plate elements pass all the cases of patch tests even though nonconforming elements do not preserve conformity. Then, the clues for the Mindlin plate element to pass the Irons' patch tests are investigated. Also, the convergent characteristics of some nonconforming Mindlin plate elements that do not pass the Irons' patch tests are examined by weak patch tests. The convergence tests are performed on the benchmark numerical problems for both nonconforming elements which pass the patch tests and which do not. Some conclusions on the relationship between the patch test and convergence of nonconforming Mindlin plate elements are drawn.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.5-25
• /
• 2018

A fully-coupled thermodynamic-based transient finite element formulation is proposed in this article for electric, magnetic, thermal and mechanic fields interactions limited to the linear case. The governing equations are obtained from conservation principles for both electric and magnetic flux, momentum and energy. A full-interaction among different fields is defined through Helmholtz free-energy potential, which provides that the constitutive equations for corresponding dual variables can be derived consistently. Although the behavior of the material is linear, the coupled interactions with the other fields are not considered limited to the linear case. The implementation is carried out in a research version of the research computer code FEAP by using 8-node isoparametric 3D solid elements. A range of numerical examples are run with the proposed element, from the relatively simple cases of piezoelectric, piezomagnetic, thermoelastic to more complicated combined coupled cases such as piezo-pyro-electric, or piezo-electro-magnetic. In this paper, some of those interactions are illustrated and discussed for a simple geometry.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.03a
• /
• pp.577-584
• /
• 2001

For construction and design of tunnels, groundwater flow models are used to find the influence of groundwater to the stability of tunnels considering the geological condition around the tunnels and the materials used in tunnel linings. For the analysis of tunnel flow, some commercial programs, e.g. MODFLOW, SEEP/W etc., are used. These programs have limitations that MODFLOW could not define curved surface smoothly in three dimensional flow media and SEEP/W is the 2-dimensional flow model. In this paper, the ability of a finite element program developed for analyzing 3-dimensional groundwater flow is examined. Confined steady state groundwater flow solution in non-homogeneous media is obtained using isoparametric element with eight trilinear hexahedron nodes and is compared with the result of MODFLOW. It was found that the solution yielded a good result with the three dimensional flow studied.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.6
• /
• pp.148-159
• /
• 1999

A numerical method is presented for the dynamic analysis of spur gears rotating with very high angular speeds. For an efficient computation each gear is assumed to consist of a rotating rigid disk and an elastic tooth having mass, and finite element formulations are used for the equations of motion of the tooth. The geometric constraint is imposed between the rigid disk and the elastic tooth to fix them, and contact condition is imposed between the meshing teeth of the gears. At each iteration of each time step the Lagrange multiplier and contact force are revised by using the constraint error vector, and then the whole equations of motion are time integrated with the given Lagrange multiplier and contact force. For the accurate solution the velocity and acceleration constraints as well as the displacement constraint are satisfied by the monotone reductions of the constraint error vectors. Computing procedures associated with the iterative schemes are explained and numerical simulations are conducted with the spur gears.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.322-330
• /
• 1987

Isotropic linear viscoelasticity problems are analyzed numerically in time domain by Boundary Element Method with quadratic isoparametric boundary elements. Viscoelastic fundamental solutions are newly derived by using the elastic-viscoelastic correspondence principle and corresponding boundary integral equations are also presented. Numerical results of two examples are compared with the derived exact solutions to verify the accuracy and validity of the method. A detailed study on the accuracy of displacement and stress in terms of time integration step is given.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.310-323
• /
• 1990

A three-noded, with three degree-of-freedom at each node, in-plane curved beam element is formulated and employed in eigen-analysis of constant curvature beam. The conventional quadratic shape functions used in a three noded C .deg. type curved beam element produce such an undesirable large stiffness that a significant error is introduced in displacements and stresses. These phenomena are called 'Stiffness Locking Phenomena', which result from spurious strain energy due to inappropriate assumptions on independent isoparametric quadratic interpolation functions. Stiffness locking phenomena can be alleviated by using modified interpolation functions which get rid of spurious constraints of conventional interpolation functions. Eigenvalues and their modes as well as displacements and stresses may be locked because they are related to stiffness. Using modified curved beam element in eigenvalue problem of cantilever and arch, the property and performance of modified curved beam element are examined by numerical experimentations. In these eigen-analyses, mass matrices are calculated by using both modified and unmodified curved beam element, are compared with theoretical solutions. These comparisons show that the performance of the modified curved beam element is better than that of the unmodified curved beam element.