• Coupled systems mechanics
• /
• 제1권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.

• 대한토목학회논문집
• /
• 제9권2호
• /
• pp.11-22
• /
• 1989

등매개(等媒介) 요소(要素)를 사용(使用)하는 유한요소법(有限要素法)을 적용(適用)하여 직선상자항교(直線箱子桁橋) 강성해석(强性解析)을 수행(遂行)하였다. 이 방법(方法)의 유용성(有用性)과 정확성(正確性)을 보이기 위하여 하나의 연속상자형을 예로 들어 해석하였다. 여기서 구한 처짐과 종방향응력들은 다른 방법으로 구한 결과와 잘 일치함을 알 수 있었다. 본 방법은 임의의 경제조건과 하중을 받는 상자형을 해석하는데 효율적인 방법이며, 설계 해석에 직접 사용될 수 있는 실용적인 방법이다.

• Steel and Composite Structures
• /
• 제15권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
• /
• 제5권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
• /
• 제7권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.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2001년도 봄 학술발표회 논문집
• /
• 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.

• 한국정밀공학회지
• /
• 제16권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
• /
• 제1권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.

• 대한기계학회논문집
• /
• 제11권2호
• /
• pp.322-330
• /
• 1987

본 논문에서는 전자의 해석중에서 유한요소와 경계요소의 결합으로 생기는 문 제점 그리고 증가된 자요도로 인한 경제성의 손실 등을 고려하여 2차의 경계요소만을 사용한 경계요소법을 시도하였으며 다른 형태의 점탄성기본해 및 점탄성 경계적분식을 유도 사용하였다. 그리고 시도된 방법의 타당성을 보이기 위하여 사용된 예제들(평 면변형문제)의 이론해를 상응원리를 사용하여 유도하였으며 이를 수치결과와 비교하여 근사해의 정확함을 알 수가 있었다.

• 대한기계학회논문집
• /
• 제14권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.