• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2543-2554
• /
• 1993

In this paper, the kinematic and dynamic analyses of spatial complex kinematic chain are studied. Through the new method both using the set of identification numbers and applying the DenavitHartenberg link representation method to the spatial complex kinematic chain, the kinematic configuration of the chain is represented. Some link in the part of closed chain being fictitiously cutted, the complex kinematic chain is transformed to the branched chain. The kinematic constraint equations are derived from the constraint conditions which the cutted sections of the link have to satisfy. And the joint variables being partitioned in the independent joint variables and the dependent joint variables, the dependent variables are calculated from the independent variables by using the Newton-Raphson iterative method and the pseudoinverse matrix. The equations of motion are derived under the independent joint variables by using the principle of virtual work. Algorithms for dynamic analysis are presented and simulations are done to verify accuracy and efficiency of the algorithms.

• The Korean Journal of Rheology
• /
• v.8 no.3_4
• /
• pp.215-225
• /
• 1996

평판형의 열가소성 수지를 유리전이온도 이상으로 가열한 다음 압력을 가함으로써 원하는 형상의 제품을 성형하는 열성형공정은 대상 수지가 큰 변형을 일으킬 뿐만 아니라 비선형적 거동을 보이게 된다. 따라서 수지의 변형거동 예측과 최적성형조건의 설정에 많은 어려움과 시행착오를 거치게 되는 바, 열성형 공정의 최적화를 위한 연구의 일환으로 원형 평판위 수지를 대상으로 수지의 부풀림 거동과 이에 따른 두께 분포를 예측할수 있는유한요 소법의 수치모사 알고리듬을 개발하고자 하였다. Piola-Kirchhoff 응력 텐서와 Green 변형 텐서 및 lagrangian 변형 텐서를 사용하여 평판상의 응력-변형에 대한 비선형의 에너지 수 지식을 수립하고 Newton-Raphson 반복수렴법을 이용하여 근사적으로 해석하였으며 수지의 유변학적 구성방정식으로는 neo-Hookean 모델, Mooney-Rivlin 모델 및 Ogden 모델등의 초탄성 모델을 사용하여 그결과를 비교하였다. 수치모사에는 두께가 매우 얇기 때문에 두께 방향의 응력변화를 무시할수 있는 membrane 가정을 도입한 2차원적 해석과 두께 방향의 응력 변화를 고려하는 3차원적 해석을 모두 수행하고 그 차이를 비교하였으며 3차원적 해석 의 경우에는 penalty법을 이용하여 비 압축성을 만족하였다. 일차적으로 내압을 받는 두꺼 운 원통계에 대한 수치모사 해석을 수행하고 완전해와 비교함으로써 개발된 수치모사 알고 리듬의 열성형 공정에의 적용 타당성을 검증하였으며 이를 이용하여 원형 평판의 자유부풀 림거동을 예측한 결과 Treloar 등의 실험결과와 잘 부합함을 확인하였다. 또 간단한 형상의 금형이 있는 경우와 반지름 방향으로의 온도변화에 따른 수지의 변형거동을 해석함으로써 실제 열성형 공정에서 요구되고 있는 성형품의 두께 분포를 균일하게 하기 위한 방안을 제 시하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2298-2305
• /
• 1992

To obtain the convenient solution of the multibody dynamic systems composed of flexible beams, linear finite element technique is adopted and the nodal coordinates are interpolated in the global inertia frame. Mass matrix becomes an extremely simple constant matrix and the force vector also becomes extremely simple because Coriolis acceleration and centrifugal force are not required. And the elastic force is also simply computed from the moving frame attached to the material. To solve the global differential algebraic euation. an ODE technique is adopted after Lagrange multiplier is computed by the accelerated iterative technique, and the time demanding procedures such as Newton-Raphson iterations and decomposition of the big matrix are not required. The accuracy of the present solution is checked by a well-known example problem.

• Journal of Korea Water Resources Association
• /
• v.41 no.9
• /
• pp.895-905
• /
• 2008

To analyze hydrologic processes in a watershed requires both various geographical data and hydrological time series data. Recently, not only geographical data such as DEM(Digital Elevation Model) and hydrologic thematic map but also hydrological time series from numerical weather prediction and rainfall radar have been provided as grid data, and there are studies on hydrologic analysis using these grid data. In this study, GRM(Grid based Rainfall-runoff Model) which is physically-based distributed rainfall-runoff model has been developed to simulate short term rainfall-runoff process effectively using these grid data. Kinematic wave equation is used to simulate overland flow and channel flow, and Green-Ampt model is used to simulate infiltration process. Governing equation is discretized by finite volume method. TDMA(TriDiagonal Matrix Algorithm) is applied to solve systems of linear equations, and Newton-Raphson iteration method is applied to solve non-linear term. Developed model was applied to simplified hypothetical watersheds to examine model reasonability with the results from $Vflo^{TM}$. It was applied to Wicheon watershed for verification, and the applicability to real site was examined, and simulation results showed good agreement with measured hydrographs.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.185-191
• /
• 2004

A mooring system can be applied to keep the position of a floating structures. In this study, the structural analysis is carried out to analyze the dynamic characteristics of a mooring line for a floating breakwater. A three-dimensional equations of motion for a submerged chain are derived. Bending stiffness is considered for the necessary restoring force in the regions of zero tension. A fortran program is to be developed by employing finite difference method. In the algorithm, an implicit time integration and Newton-Raphson iteration are adopted. The results of simulation show good agreement in tension response pattern with the experimental results of a reference. The results of this study can contribute for the design of mooring system for a floating breakwater.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.1 no.1
• /
• pp.93-101
• /
• 1989

A geometrically nonlinear cable finite element is presented to use in the static or dynamic modeling of offshore and onshore structures such as guyed tower, tension leg platform or mooring buoy, submarine cable, cable-stayed bridge, suspension bridge, cable roof and so on. The cable finite element is derived directly from the compatibility equations and flexibility matrix of elastic catenary cable theory for the arbitary plane loading and geome try. A general and virsatile computer program has been developed to perform the analyses of cable member itself or cable guyed or suspened structures, in which Newmark-$\beta$ method is used to obtain a time domain solution and Newton-Raphson iteration method is used to solve the nonlinear system of compatibility equations of cable and algebraic static or dynamic equations at each time step. The results from the static and dynamic analysis of a cable member by the computer program are summarized and presented.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1998.04a
• /
• pp.17-17
• /
• 1998

연소반응을 가지는 후류(wake)는 가스터빈 연소실의 flame holder 등에서 발생한다. 후류유동의 안정성 혹은 불안정성은 이러한 유동에 있어서 많은 영향을 끼치므로 상당히 중요하다. 본 연구는 위와 같이 연소 반응에 의한 밀도구배를 가지는 후류유동에 대하여 불안정성 해석을 수행하였다. 교란에 대한 지배방정식은 Navier-Stokes 방정식에서 점성항을 제외한 Euler 방정식을 고려하였다. 충류유동의 압력은 일정하다고 가정하였다. 교란 방정식을 유도하기 위하여 충류 유동이 평행하여 유동 방향에 수직한 방향의 구배만이 존재한다고 가정하였다. 모든 변수들은 충류 유동의 값과 움직이는 파장의 형태를 가지는 작은 교란의 합으로 생각하여 압력에 대한 교란방정식을 구하며, shooting법과 Newton-Raphson법에 근거를 두는 반복법을 사용하여 풀었다. 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 비압축성의 경우 우선 Gaussian Profile 가정함과 동시에 연소반응을 포함하는 유동장을 보다 정확히 구하기 위하여 Navier-Stokes 방정식으로부터 구한 결과를 사용하였다. 연소반응을 포함하는 유동장을 구할 때에는 해석상 편의를 위해 예혼합물질은 이상기체로, 화학반응은 1단계의 비가역반응으로 가정하였으며, 반응열로 인한 부력의 효과는 무시하였다. 위와 같은 유동장을 가지고 불안정성 해석을 수행한 결과 후류유동은 두 개의 변곡점을 가지며 sinuous 모드와 varicose모드의 두 개의 불안정성 모드가 있음을 보였다. 밀도 변화가 있는 경우나 밀도 변화가 없는 경우 모두 sinuous 모드의 가장 불안정한 모드가 varicose 모드의 가장 불안정한 모드보다 더 불안정함을 보여주어 후류 유동은 자유 유동에 가까운 위상 속도를 가지는 sinuous 모드에 의해 지배될 것임을 예측할 수 있다. 연소반응이 완전연소에 가까울수록 그리고 압축성 효과가 클수록 유동내부의 온도가 증가하고 점성 또한 증가하여 후류유동은 안정됨을 알 수 있었다 유동변수들의 contour로부터 유동의 특성을 예측한 결과 baroclinic 항이 dilatational 항보다 상대적으로 크며, 중심선 상하에 생기는 vortex를 더욱 성장시킬 것으로 생각된다.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.2 no.2
• /
• pp.87-94
• /
• 1998

Behavioral characteristics of shallow circular arches with dynamic loading and different end conditions are analysed. Geometric nonlinearity is modelled using Lagrangian description of the motion. The finite element analysis procedure is used to solve the dynamic equation of motion, and the Newmark method is adopted in the approximation of time integration. The behavior of arches is analysed using the buckling criterion and non-dimensional time, load and shape parameters which Humphreys suggested. But a new deflection-ratio formula including the effect of horizontal displacement plus vertical displacement is presented to apply for the non-symmetric buckling problems. Through the model analysis, it's confirmed that fix-ended arches have higher buckling stability than hinge-ended arches, and arches with the same shape parameter have the same deflection ratio at the same time parameter when loaded with the same parametric load.

• Nuclear Engineering and Technology
• /
• v.20 no.2
• /
• pp.88-104
• /
• 1988

The FUel Rod Analysis(FURA) code is developed using two-dimensional finite element methods for axisymmetric and plane stress analysis of fuel rod. It predicts the thermal and mechanical behavior of fuel rod during normal and load follow operations. To evaluate the exact temperature distribution and the inner gas pressure, the radial deformation of pellet and clad, the fission gas release are considered over the full-length of fuel rod. The thermal element equation is derived using Galerkin's techniques. The displacement element equation is derived using the principle of virtual works. The mechanical analysis can accommodate various components of strain: elastic, plastic, creep and thermal strain as well as strain due to swelling, relocation and densification. The 4-node quadratic isoparametric elements are adopted, and the geometric model is confined to a half-pellet-height region with the assumption that pellet-pellet interaction is symmetrical. The pellet cracking and crack healing, pellet-cladding interaction are modelled. The Newton-Raphson iteration with an implicit algorithm is applied to perform the analysis of non-linear material behavior accurately and stably. The pellet and cladding model has been compared with both analytical solutions and experimental results. The observed and predicted results are in good agreement. The general behavior of fuel rod is calculated by axisymmetric system and the cladding behavior against radial crack is used by plane stress system. The sensitivity of strain aging of PWR fuel cladding tube due to load following is evaluated in terms of linear power, load cycle frequency and amplitude.