• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.1
• /
• pp.79-88
• /
• 1997

Very stiff floor system in a residential-commercial building causes some problems in the numerical analysis procedure due to significant difference in stiffness with adjacent elements. Static analysis of structure with a stiff transfer-floor can be performed approximately in two steps for upper and lower parts for the structure. However, it is impossible to perform dynamic analysis in two steps with separate models. An efficient method for dynamic analysis of a structure with a right floor system is proposd in this study. The matrix condensation technique is employed to reduce the degree of freedom for upper and lower parts of the structure and a beam elements with rigid bodies at both ends are introduce to model the rigid floor system. Efficiency and accuracy of the proposed method are verified through analysis of several example structures.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.2
• /
• pp.250-255
• /
• 2008

Springs are widely utilized in machine element. To find out stiffness of frustum-shaped coil spring, the space beam theory using the finite element method is adopted in this paper. In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. To find out load vector of coil spring subjected to distributed compression. principle of virtual work is adapted. And this theory was programming using MATLAB software. To compare FEM using MATLAB software was applied MSC. Nastran software. The geometry model for MSC. Patran was produced by 3-D design modeling software. Finite element model was produced by MSC. Patran. Finite element was applied tetra (CTETRA) having 10 node. The analysis results of the MATLAB and MSC. Nastran are fairly well agreed with those of various experiments. Using MATLAB program proposed in this paper and MSC. Nastran, spring constants and stresses can be predicted by input of few factors.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.169-179
• /
• 2005

An improved method for evaluating effective buckling lengths of beam-column members in plane frames is newly proposed based on system inelastic buckling analysis. To this end, the tangent stiffness matrix of be am-column elements is first calculated using stability functions and then the inelastic buckling analysis method is presented. The scheme for determining effective length of individual members is also addressed. Design examples and numerical results ?uc presented to show the validity of the proposed method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1A
• /
• pp.31-43
• /
• 2009

In this study, the numerical method is presented, which can consider the various train types and can solve the equations of motion for a vehicle-bridge interaction analysis by non-iteration procedure through formulating the coupled equations of motion. The coupled equations of motion for the vehicle-bridge interaction are solved by the Newmark ${\beta}$ of a direct integration method, and by composing the effective stiffness matrix and the effective force vector according to a analysis step, those can be solved with the same manner of the solving procedure of equilibrium equations in static analysis. Also, the effective stiffness matrix is reconstructed by the Skyline method for increasing the analysis effectiveness. The Cholesky's matrix decomposition scheme is applied to the analysis procedure for minimizing the numerical errors that can be generated in directly calculating the inverse matrix. The equations of motion for the conventional trains are derived, and the numerical models of the conventional trains are idealized by a set of linear springs and dashpots with 16 degrees of freedom. The bridge models are simplified by the 3 dimensional space frame element which is based on the Euler-Bernoulli theory. The rail irregularities of vertical and lateral directions are generated by the PSD functions of the Federal Railroad Administration (FRA). The results of the vehicle-bridge interaction analysis are verified by the experimental results for the railway plate girder bridges of a span length with 12 m, 18 m, and the experimental and analytical data are applied to the low pass filtering scheme, and the basis frequency of the filtering is a 2 times of the 1st fundamental frequency of a bridge bending.

• Journal of the Korea Concrete Institute
• /
• v.17 no.6 s.90
• /
• pp.939-946
• /
• 2005

The stiffness properties of heterogeneous concrete materials and their degradation were investigated at different-levels of observations with aids of the opportunities and limitations of multi-resolution wavelet analysis. The successive Haw transformations lead to a recursive separation of the stiffness properties and the response into coarse-and fine-scale features. In the limit, this recursive process results in a homogenization parameter which is an average measure of stiffness and strain energy capacity at the coarse scale. The basic concept of multi-resolution analysis is illustrated with one and two-dimensional model problems of a two-phase particulate composite representative of the morphology of concrete materials. The computational studies include the meso-structural features of concrete in the form of a hi-material system of aggregate particles which are immersed in a hardened cement paste taking due to account of the mismatch of the two elastic constituents.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.10 no.6 s.52
• /
• pp.19-28
• /
• 2006

Most conventional model updating methods must use mathematical objective function with experimental modal matrices and analytical system matrices or must use information about the gradient or higher derivatives of modal properties with respect to each updating parameter. Therefore, most conventional methods are not appropriate for complex structural system such as bridge structures due to stability problem in inverse analysis with ill-conditions. Sometimes, moreover, the updated model may have no physical meaning. In this paper, a new FE model updating method based on a hybrid optimization technique using genetic algorithm (GA) and Holder-Mead simplex method (NMS) is proposed. The performance of hybrid optimization technique on the nonlinear problem is demonstrated by the Goldstein-Price function with three local minima and one global minimum. The influence of the objective function is evaluated by the case study of a simulated 10-dof spring-mass model. Through simulated case studies, finally, the objective function is proposed to update mass as well as stiffness at the same time. And so, the proposed hybrid optimization technique is proved to be an efficient method for FE model updating.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.3
• /
• pp.344-352
• /
• 2002

In the analysis of the 3 dimensional plate structures by the finite element method, the triangular elements are generally used for the global stiffness matrix of the analyzed system. But the triangular elements of the plates have some problems in the process of formulation and in the precision of analysis. The formulation of the finite element method to analyze 3 dimensional plate structures using quadrilateral elements is presented in this paper. The degree of freedom off nodal point is 6, that is, the displacements in the direction off-y-z is and the rotations about x-y-z axis and then the degree of freedom off element is 24. For the comparison of the analysis using triangular elements and quadrilateral elements, the rectangular plates subjected to the uniform load and a concentrated load on the centroid of the plate, for which the theoretical solutions have been obtained, are analyzed. The calculated deflections of the rectangular plates using the finite element method by the triangular elements and the quadrilateral elements are also compared with the deflections of the plates calculated by theoretical solutions. The defections of the rectangular plates calculated by the finite element method using the quadrilateral elements are closer to the theoretical solutions than the defections calculated by the finite element method using the triangular elements. The deflection of the centroid of plate, calculated by the finite element method, converges to that of theoretical solution as the number of elements is increased. This convergence is much more rapid for the case of using the quakrilateral elements than fir the case of using triangular elements.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.25-30
• /
• 1992

최근 전자계산기를 이용한 진동해석 방법이 눈부시게 발달하여, 일반 구조물 이나 기계 구조물 등의 동특성을 설계 단계에서 정도 높게 예측하는 것이 가능하게 되었다. 그러나 종래의 구조해석은 주어진 시스템의 동특성을 위한 것으로 얻어진 동특성으로부터 질량, 관성제원 및 스프링상수값 등의 설계상 수값을 규명하는 연구는 미미한 실정이다. 이것에 대한 해결방법으로 크게 해석적인 방법과 실험적인 방법으로의 접근이 있어 왔다. 해석적인 방법으로 유한요소해석에서 얻은 모드좌표를 물리좌표로 변환하는 방법으로 Guyan의 정축소와 같은 절점축소를 행하는 방법이 고찰되었다. 실험적인 방법으로 가 진실험에서 얻은 전달함수나 모드파라미터로부터 [M], [K] 행렬을 결정하는 연구가 있었지만 어떤것도 질량, 스프링상수 등의 설계상수를 완전히 규명하 지는 못하였다. 또한, 설계 단계에서 필요한 질량, 관성제원 또는 스프링상수 등의 최적한 값이나, 원하는 시스템특성을 얻을 수 있는 설계상수의 적정한 폭을 구하는 연구는 설계자의 경험과 반복된 시행착오에 의존하는 실정이다. 감도해석은 이러한 문제점을 개선하는 수단으로 설계변수에 대한 동특성의 변화율을 구하는 것이다. 감도해석을 수행하는 것은 어느 설계변수를 수정하 는 것이 주어진 동특성에 부합되는 지를 알려주고, 어느 것을 수정하는 것이 원하는 방향의 동특성변화에 가장 효과적인지를 알려주는 것이다. 따라서 감 도해석을 이용하여 설계의 최적화 프로그램을 만들수 있고, 이것은 설계자가 요구하는 동특성을 목적함수로 하여 주어진 구조물을 최적화하는 설계상수 값을 얻을 수 있게 한다. 본 논문에서는 강체모델의 동특성으로부터 모델의 설계 상수를 규명하고, 동특성의 개선을 위하여 설계변수의 변경량을 물리좌 표계에서 얻는것을 목적으로 한다. 강체 마운트계의 관성제원 및 마운트강성 의 규명을 위하여 임으로 주어진 설계상수를 모델데이타로 하여 관성제원과 스프링 강성을 구하였다. 관성제원의 규명은 주어진 모델의 관성값을 모르는 것으로 하여 임의의 초기 관성값으로 감도해석에 의해 주어진 계의 관성값 을 물리 좌표계에서 규명하였다. 마운트 강성의 규명도 관성제원의 규명과 같은 방법으로 임의의 강성값으로 감도해석을 하여 강성값을 규명하였다. 또 한 감도해석에 의한 동특성 변경은 특정한 고유진동 수의 변경이 필요할 때, 고유진동수의 이동을 위한 관성제원의 변경 및 마운트 강성변경값을 예측할 수 있다. 본 연구수행의 기본적인 흐름도는 Fig.1.1과 같다. 위와 같은 작업 으로 엔진 마운트와 같은 강체 모델의 시스템 규명을 행하는 경우에 유한요 소해석 및 가진 실험으로 얻은 고유진동수의 정보 또는 원하는 고유진동수 의 특성을 기본으로 실제 설계에서 사용이 가능하도록 물리 좌표계에서 관 성 제원 및 스프링상수를 구할 수 있을 것이다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.307-314
• /
• 2009

The present paper briefly describes the space frame element and the fundamental strategies in computational elastic bifurcation theory of geometrically nonlinear, single load parameter conservative elastic spatial structures. A method for large deformation(rotation) analysis of space frame is based on an eulerian formulation, which takes into consideration the effects of large joint translations and rotations with finite deformation(rotation). The local member force-deformation relationships are based on the beam-column approach, and the change in member chord lengths caused by axial strain and flexural bowing are taken into account. and the derived geometric stiffness matrix is unsymmetric because of the fact that finite rotations are not commutative under addition. To detect the singular point such as bifurcation point, an iterative pin-pointing algorithm is proposed. And the path switching mode for bifurcation path is based on the non-negative eigen-value and it's corresponding eigen-vector. Some numerical examples for bifurcation analysis are carried out for a plane frame, plane circular arch and space dome structures are described.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.301-308
• /
• 2013

In this paper, an accuracy analysis of parallel method based on non-overlapping domain decomposition method is carried out. In this approach, proposed by Tak et al.(2013), the decomposed subdomains do not overlap each other and the connection between adjacent subdomains is determined via simple connective finite element named interfacial element. This approach has two main advantages. The first is that a direct method such as gauss elimination is available even in a singular problem because the singular stiffness matrix from floating domain can be converted to invertible matrix by assembling the interfacial element. The second is that computational time and storage can be reduced in comparison with the traditional finite element tearing and interconnect(FETI) method. The accuracy of analysis using proposed method, on the other hand, is inclined to decrease at cross points on which more than three subdomains are interconnected. Thus, in this paper, an accuracy analysis for a novel non-overlapping domain decomposition method with a variety of subdomain numbers which are interconnected at cross point is carried out. The cause of accuracy degradation is also analyze and establishment of countermeasure is discussed.