• Structural Engineering and Mechanics
• /
• 제51권6호
• /
• pp.923-937
• /
• 2014

The rational finite element method is different from the standard finite element method, which is constructed using basic solutions of the governing differential equations as interpolation functions in the elements. Therefore, it is superior to the isoparametric approach because of its obvious physical meaning and accuracy; it has successfully been applied to the isotropic elasticity problem. In this paper, the formulation of rational finite elements for plane orthotropic elasticity problems is deduced. This method is formulated directly in the physical domain with full consideration of the requirements of the patch test. Based on the number of element nodes and the interpolation functions, different approaches are applied with complete polynomial interpolation functions. Then, two special stiffness matrixes of elements with four and five nodes are deduced as a representative application. In addition, some typical numerical examples are considered to evaluate the performance of the elements. The numerical results demonstrate that the present method has a high level of accuracy and is an effective technique for solving plane orthotropic elasticity problems.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2009년도 춘계학술대회 논문집
• /
• pp.573-578
• /
• 2009

Starting from the rotating beam finite element in which the interpolating shape functions satisfies the governing static homogeneous differential equation of Euler-Bernoulli rotating beams, we derived new shape functions that satisfies the governing differential equation which contains the terms of hub radius and setting angle. The shape functions are rational functions which depend on hub radius, setting angle, rotational speed and element position. Numerical results for uniform and tapered cantilever beams with and without hub radius and setting angle are compared with the available results. It is shown that the present element offers an accurate method for solving the free vibration problems of rotating beam.

• 한국소음진동공학회논문집
• /
• 제19권6호
• /
• pp.591-598
• /
• 2009

Starting from the rotating beam finite element in which the interpolating shape functions satisfy the governing static homogeneous differential equation of Euler-Bernoulli rotating beams, we derived new shape functions that satisfy the governing differential equation which contains the terms of hub radius and setting angle. The shape functions are rational functions which depend on hub radius, setting angle, rotational speed and element position. Numerical results for uniform and tapered cantilever beams with and without hub radius and setting angle are compared with the available results. It is shown that the present element offers an accurate method for solving the free vibration problems of rotating beams.

• Structural Engineering and Mechanics
• /
• 제28권2호
• /
• pp.239-257
• /
• 2008

Double-sided punched metal plate timber fasteners present projections on both sides, which offer improved joint fire resistance and better joint aesthetics. In this paper, 3-D nonlinear finite element models were developed to simulate double-sided nail plate fastener timber joints. The models, incorporating orthotropic elasticity, Hill's yield criterion and elasto-plasticity and contact algorithms, are capable of simulating complex contact between the tooth and the timber and between the base plate and the timber in a fastener. Using validated models, parametric studies of the double-sided nail plate joints was undertaken to cover the tooth length and the tooth width. Optimal configuration was assumed to have been attained when increase in nail plate tooth width did not result in a raise in joint capacity, in conjunction with the optimum tooth length. This paper presents the first attempt to model and optimise tooth profile of double-sided nail plate fastener timber joints, which offers rational designs of such fasteners.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• 한국정밀공학회지
• /
• 제23권6호
• /
• pp.128-135
• /
• 2006

Rational dynamic modeling and analysis method f3r complex structures are studied with special attention to slide way joints. For modeling of slide way joints, a general modeling technique is used by using the influence coefficients method which is applied to the conversion of detailed finite element model to the equivalent reduced joint model. The theoretical part of this method is illustrated and the method is applied to the structure with slide way joint. In this method, the non-linearity of the contact surfaces is considered within a proper range and the boundary effect of the joint model could be eliminated. The proposed method was applied to finite element modal analysis of a clamp jointed cantilever beam and slide way joints of the vertical type lathe. The method can also be used to other kinds of joint modeling. The results of these analysis were compared with those of Yoshimura models and rigid joint models, which demonstrated the practical applicability of the proposed method.

• Structural Engineering and Mechanics
• /
• 제25권1호
• /
• pp.91-106
• /
• 2007

A new class of finite elements is described for dealing with mesh gradation. The approach employs the moving least square (MLS) scheme to devise a class of elements with an arbitrary number of nodal points on the parental domain. This approach generally leads to elements with rational shape functions, which significantly extends the function space of the conventional finite element method. With a special choice of the nodal points and the base functions, the method results in useful elements with polynomial shape functions for which the $C^1$ continuity breaks down across the boundaries between the subdomains comprising one element. Among those, (4 + n)-noded MLS based finite elements possess the generality to be connected with an arbitrary number of linear elements at a side of a given element. It enables us to connect one finite element with a few finite elements without complex remeshing. The effectiveness of the new elements is demonstrated via appropriate numerical examples.

• Structural Engineering and Mechanics
• /
• 제30권1호
• /
• pp.119-129
• /
• 2008

A small strain and elastoplastic formulation of Polygonal Element Method (PEM) is developed for efficient analysis of elastoplastic solids. In this work, the polygonal elements are constructed based on traditional triangular finite meshes. The construction method of polygonal mesh can directly utilize the sophisticated triangularization algorithm and reduce the difficulty in generating polygonal elements. The Wachspress rational finite element basis function is used to construct the approximations of polygonal elements. The incremental variational form and a von Mises type model are used for non-linear elastoplastic analysis. Several small strain elastoplastic numerical examples are presented to verify the advantages and the accuracy of the numerical formulation.

• Advances in Computational Design
• /
• 제9권2호
• /
• pp.81-96
• /
• 2024

Modeling and analyzing the dynamic behavior of fluid-soil-structure interaction problems are crucial in structural engineering. The solution to such coupled engineering systems is often not achievable through analytical modeling alone, and a numerical solution is necessary. Generally, the Finite Element Method (FEM) is commonly used to address such problems. However, when dealing with coupled problems with complex geometry, the finite element method may not precisely represent the geometry, leading to errors that impact solution quality. Recently, Isogeometric Analysis (IGA) has emerged as a preferred method for modeling and analyzing complex systems. In this study, IGA based on Non-Uniform Rational B-Splines (NURBS) is employed to analyze the seismic behavior of concrete gravity dams, considering fluid-structure-foundation interaction. The performance of IGA is then compared with the classical finite element solution. The computational efficiency of IGA is demonstrated through case studies involving simulations of the reservoir-foundation-dam system under seismic loading.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집A
• /
• pp.324-329
• /
• 2001

A novel method for non-matching interfaces on the boundaries of the finite elements in partitioned domains is presented by introducing interface elements in this paper. The interface element method (IEM) satisfies the continuity conditions exactly through interfaces without recourse to the Lagrange multiplier technique. The moving least square (MLS) approximation in the present study is implemented to construct the shape functions of the interface elements. Alignment of the boundaries of sub-domains in the MLS approximation and integration domains provides a consistent numerical integration due to one form of rational functions in an integration domain. The compatibility of displacements on the boundaries of the finite elements and the interface elements is always preserved in this method, and the completeness of the shape functions of the interface elements guarantees the convergence of numerical solutions. The numerical examples show that the interface element method is a useful tool for the analysis of a partitioned system and for a global-local analysis.