• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1997년도 가을 학술발표회 논문집
• /
• pp.42-49
• /
• 1997

In this paper, the method for soil-structure interaction analyses in the time domain is proposed. The far field soil region which is the outside of the artificial boundary is modeled by using explicit frequency-dependent two dimensional infinite elements which can include multiple wave components propagating into the unbounded medium. Since the dynamic stiffness matrix of the far field soil region using the proposed infinite elements is obtained explicitly in terms of exciting frequencies and constants in the frequency domain, the matrix can be easily transformed into the displacement unit-impulse response matrix, which corresponds to a convolution integral of it in the time domain. To verify the proposed method for soil-structure interaction analyses in the time domain, the displacement responses due to an impulse load on the surface of a soil layer with the rigid bed rock are compared with those obtained by the method in the frequency domain and those by models with extend finite element meshes. Good agreements have been found between them.

• Wind and Structures
• /
• 제10권1호
• /
• pp.61-82
• /
• 2007

This paper presents a novel finite element (FE) model for analyzing coupled flutter of long-span bridges using the commercial FE package ANSYS. This model utilizes a specific user-defined element Matrix27 in ANSYS to model the aeroelastic forces acting on the bridge, wherein the stiffness and damping matrices are expressed in terms of the reduced wind velocity and flutter derivatives. Making use of this FE model, damped complex eigenvalue analysis is carried out to determine the complex eigenvalues, of which the real part is the logarithm decay rate and the imaginary part is the damped vibration frequency. The condition for onset of flutter instability becomes that, at a certain wind velocity, the structural system incorporating fictitious Matrix27 elements has a complex eigenvalue with zero or near-zero real part, with the imaginary part of this eigenvalue being the flutter frequency. Case studies are provided to validate the developed procedure as well as to demonstrate the flutter analysis of cable-supported bridges using ANSYS. The proposed method enables the bridge designers and engineering practitioners to analyze flutter instability by using the commercial FE package ANSYS.

• Advances in Computational Design
• /
• 제5권3호
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• Journal of Advanced Marine Engineering and Technology
• /
• 제23권5호
• /
• pp.702-710
• /
• 1999

U-shaped bellows are usually used to piping system pressure sensor and controller for refriger-ator. Bellows subjected to internal pressure are designed for the purpose of absorbing deformation. Internal pressure on the convolution sidewall and end collar will be applied to an axial load tend-ing to push the collar away from the convolutions. To find out deformation behavior of bellow sub-jected to internal pressure the axisymmetric shell theory using the finite element method is adopted in this paper. U-shaped bellows can be idealized by series of conical frustum-shaped ele-ments because it is axisymmetric shell structure. The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displace-ments are added to r-z cylindrical coordinates of nodal points. The new stiffness matrix of the sys-tem using the new coordinates of nodal points is adopted to calculate the another increments of nodal displacement that is the step by step method is used in this paper. The force required to deflect bellows axially is a function of the dimensions of the bellows and the materials from which they are made. Spring constant is analyzed according to the changing geometric factors of U-shaped bellows. The FEM results were agreed with experiment. Using developed FORTRAN PROGRAM the internal pressure vs. deflection characteristics of a particu-lar bellows can be predicted by input of a few factors.

• 지구물리와물리탐사
• /
• 제5권3호
• /
• pp.169-177
• /
• 2002

일반 반사법 탐사 모델링에서 효율적인 주파수영역 수치모델링의 실용화를 위해 무엇보다 해결해야할 과제는 파장당 격자수를 줄이는 것이다. 본 연구에서는 이에 착안하여 수치분산 및 수치이방성을 최소화시키면서 한 파장당 필요한 격자수를 줄일 수 있는 가중평균 유한요소법을 개발하였다. 강성행렬과 질량행렬은 네 개의 사각형 유한요소로 확장하였으며, 모든 격자점이 포함되도록 유한요소를 배열하여 조합하였다. 확장된 네 개의 강성행렬과 질량행렬은 가중평균계수를 주어 선형결합하는 방법으로 가중평균하였다. 가중평균계수는 확장된 25점 평균차분법을 사용하여 가중평균계수를 결정하였다. 또한, 정확도 향상은 2차원 균질모델 과 수평층 모델에서 해석해와 한 파장당 4개의 격자점을 준 가중평균 유한요소법 수치해 비교를 통하여 검증하였다. 또한, 석유탐사에서 활용성이 높은 향사구조 모델을 선정하여 이의 반응을 관찰한 결과 지층경계면외에 네 개의 사각형 유한요소들의 구성으로 인한 인위적인 파의 도달이 인식되지 않았다. 따라서, 본 연구에서 고안된 가중평균 유한요소법은 주파수영역에서 폭 넓은 수치모델링연구을 가능하게 할 것이다.

• 한국농공학회지
• /
• 제22권2호
• /
• pp.83-100
• /
• 1980

This study was intended to computerize the structural analysis of two-dimensional continuum by finite element method, and to provide a preparatory basis for more sophisticated and more generalized computer programs of this kind. A computer program, applicable to any shape of two-dimensional continuum, was formulated on the basis of 16-degree-of- freedom rectangular element. Various computational aspects pertaining to the implementation of finite element method were reviewed and settled in the course of programming. The validity of the program was checked through several case studies. To assess the accuracy and the convergence characteristics of the method, the results computed by the program were compared with solutions by other methods, namely the analytical Navier's method and the framework method. Through actual programming and analysis of the computed results, the following facts were recognized; 1) The stiffness matrix should necessarily be assembled in a condensed form in order to make it possible to discretize the continuum into practically adequate number of elements without using back-up storage. 2) For minimization of solution time, in-core solution of the equilibrium equation is essential. LDLT decomposition is recommended for stiffness matrices condensed by the compacted column storage scheme. 3) As for rectangular plates, the finite element method shows better performances both in the accuracy and in the rate of convergence than the framework method. As the number of elements increases, the error of the finite element method approaches around 1%. 4) Regardless of the structural shape, there is a uniform tendency in convergence characteristics dependent on the shape of element. Square elements show the best performance. 5) The accuracy of computation is independent of the interpolation function selected.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 2002년도 ICCAS
• /
• pp.62.5-62
• /
• 2002

The main purpose of this paper is to identify the important parameters and to examine their relations to one another while gantry crane structure was modeled. The important elements of the structural analysis are included, such as the stiffness matrix and its relations to the degrees of freedom, the displacement, and frequency responses. To investigate these relations, the parametric modeling of a dynamic system is solved by using the finite element method (ANSYS-Program). Furthermore, EXPRESS schema and C-FAR (change favorable representation) are described how to change the frame length of gantry crane which influences other elements. Since this relationship is established, the results may...

• Architectural research
• /
• 제13권1호
• /
• pp.17-24
• /
• 2011

Isogeometric solution of Poisson equation is provided. NURBS (NonUniform B-spline Surface) is introduced to express both geometry of structure and unknown field of governing equation. The terms of stiffness matrix and load vector are consistently derived with very accurate geometric definition. The validity of the isogeometric formulation is demonstrated by using two numerical examples such as square plate and L-shape plate. From numerical results, the present solutions have a good agreement with analytical and finite element (FE) solutions with the use of a few cells in isogeometric analysis.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1997년도 추계학술대회논문집
• /
• pp.28-31
• /
• 1997

Buckling analysis of thin compressed beam has been carried out. Pre-buckling and post-buckling are simulated by finite element method incorporating with the incremental nonlinear theory and the Newton-Raphson solution technique. In order to find the bifurcation point, the determinent of the stiffness matrix is calculated at every iteration procedure. For post-buckling analysis, a small perturbed initial guess is given along the eigenvector direction at the bifurcation point. Nonlinear elastic buckling and elastic-plastic buckling of cantilever beam are analyzed. The buckling load and buckled shape of the two models are compared.

• 소음진동
• /
• 제5권4호
• /
• pp.565-575
• /
• 1995

Gearboxes are often used in the petrochemical and electrical power plants to transmit mechanical power between two branches of a machinery train rotating at different speeds. When the gearboxes are connected with rotors supported by journal bearings, bearing loads vary in magnitude and direction with rotor speed and torque transmitted by the gearboxes. In this study, dynamic characteristics of the system which consists of gearboxes and a rotor supported by journal bearings are investigated analytically and experimentally by employing the polynomial transfer matrix method and modal analysis under different speeds and torque levels. Journal bearing loads due to the transmitted torque are claculated analytically and the stiffness and damping coefficient of the journal bearings are obtained using finite element method. Comparison of the analytical and experimental results shows that the cross coupled stiffness coefficients increase with increasing rotor speed, while the cross coupled damping coefficients decrease. This generates the oil whirl instability in the journal bearings. As the transmitted torque level goes up, the stiffness coefficients of the journal bearing and the first horizontal natural frequency increase. High levels of the transmitted torque produce high bearing stiffness since the contact loads of the mating gear teeth increase. The logarithmic decrement, which is a stability indicator, is shown to decrease with increasing speed and decreasing torque. Thus, at the low torque level, the system become unstable even at the low shaft speed.