• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1987년도 전기.전자공학 학술대회 논문집(I)
• /
• pp.185-189
• /
• 1987

An analysis and design of large-scale linear multivariable systems often requires to be block triangularized form for good sensitivity of the systems when their poles and zeros are varied. But the decomposition algorithms presented up to now need a procedure of permutation, rescaling and a solution of nonlinear algebraic equations, which are usually burden. To avoid these problem, in this paper we develop a newly alternative block triangular decomposition algorithm which used the generalized matrix sign function on the Z-plane. Also, the decomposition algorithm demonstrated using the fifth order linear model of a distillation tower system.

• 호남수학학술지
• /
• 제43권3호
• /
• pp.373-383
• /
• 2021

In this study, a numerical method deals with the Chebyshev wavelet collocation and Adomian decomposition methods are proposed for solving Korteweg-de Vries equation. Integration of the Chebyshev wavelets operational matrices is derived. This problem is reduced to a system of non-linear algebraic equations by using their operational matrix. Thus, it becomes easier to solve KdV problem. The error estimation for the Chebyshev wavelet collocation method and ADM is investigated. The proposed method's validity and accuracy are demonstrated by numerical results. When the exact and approximate solutions are compared, for non-linear or linear partial differential equations, the Chebyshev wavelet collocation method is shown to be acceptable, efficient and accurate.

• 한국산학기술학회논문지
• /
• 제14권10호
• /
• pp.4733-4738
• /
• 2013

본 연구에서는 구속조건을 가진 기계계의 동적 평형위치를 다물체 동역학 해석방법을 이용하여 계산하였다. 다물체계에서 얻어지는 시간 구속조건을 가진 구속조건식과 동역학식으로부터 독립좌표계로 이루어진 동적평형식을 유도하였다. 동적 평형식은 구속조건식과 함께 비선형 대수방정식의 형태로서 Newton-Raphson 방법을 이용하여 수치해를 구하였다. 제안된 동적 평형 계산 방법을 조속기에 적용하여 동적 평형위치를 구하였다. 해석결과는 상용 프로그램의 동역학해석을 통한 평형위치의 결과와 비교하여 타당성을 검증하였다. 조속기의 회전 각속도에 대한 평형위치를 계산하고 설계 파라미터에 대한 평형위치의 영향을 분석하였다.

• 한국해안해양공학회지
• /
• 제1권1호
• /
• pp.93-101
• /
• 1989

본 논문에서는 가이드 타워, 텐션 레그 프랫폼, 무어링 부이, 해저 케이블, 사장교, 현수교, 케이블 루프 등과 같은 해상 및 육상 구조물의 유한요소 모델에 사용하기 위한 기하학적 비선형 케이블 요소를 연구 제시하였으며, 케이블 요소는 평면내에서 임의의 하중과 기하형상을 갖는 케이블에 대한 탄성현수 케이블 이론으로부터의 적합방정식과 연성행렬을 직접 이용하여 유도하였다. 또한, 유도된 케이블 유한요소에 근거하여, 케이블 부재를 사용하는 구조물들의 유한요소 해석을 위해 전산 프로그램을 개발하였으며, 시간영역 동적 해석을 위해 뉴마크-베타의 직접적분법을 사용하였고, 각 시간간격에서의 비선형 평형방정식 및 적합방정식을 풀기 위한 방법으로서 뉴톤-랩슨의 반복법을 사용하였다. 이상과 같이 개발된 전산 프로그램을 이용하여 케이블 부재에 대한 정적 및 동적 해석을 수행한 후 그 결과를 분석ㆍ고찰하여 보았다.

• 대한기계학회논문집A
• /
• 제27권9호
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제13권1호
• /
• pp.833-847
• /
• 2021

The horizontal porous baffle and its effect as an anti-slosh device have been investigated intensively in a swaying and rolling rectangular tank. To accurately assess the level at which porous baffles reduce liquid sloshing, the Matched Eigenfunction Expansion Method (MEEM) has been utilized as an analytical tool. The velocity potentials in the horizontal baffle-covered fluid region are expressed by the sum of the homogeneous and particular solutions to avoid solving the complex dispersion equation. By applying an equivalent linearized quadratic loss model, the nonlinear algebraic equation is derived and solved by implementing the Newton-Raphson iterative scheme. To prove the validity of the present theoretical model, a series of experiments have been conducted with different centered horizontal porous baffles with varying porosities and submerged depths in a swaying and rolling rectangular tank. Reasonably good agreements are obtained regarding the analytical solutions and the experiment's findings. The influence of porosity, submerged depth, and length of a centered horizontal porous baffle on anti-slosh performance have been analyzed, especially at resonance modes. The developed predictive tool can potentially provide guidelines for optimal design of the horizontal porous baffle.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집D
• /
• pp.319-324
• /
• 2001

A finite element model for the process of squeeze casting for metal matrix composites (MMCs) in cylindrical mold is developed. The fluid flow and the heat transfer are the fundamental phenomena in the squeeze casing process. To describe heat transfer with solidification of molten aluminum, the energy equation in terms of temperature and enthalpy are applied to two dimensional axisymmetric model which is similar to the experimental system. And one dimensional flow model is employed to simulate the transient metal flow. The direct iteration technique was used to solve the resulting nonlinear algebraic equations. A computer program is developed to calculate the enthalpy, temperature and fluid velocity. Cooling curves and temperature distribution during infiltration and solidification are calculated for pure aluminum. The temperature is measured and recorded experimentally. At two points of the perform inside and one point of the mold outside, thermocouple wire are installed. The time-temperature data are compared with the calculated cooling curves. The experimental results show that the finite element model can estimate the solidification time and predict the cooling process.

• Journal of Electrical Engineering and Technology
• /
• 제8권6호
• /
• pp.1530-1541
• /
• 2013

The coexistence and extinction of species are important concepts for biological systems and can be distinguished by an investigation of stability. When determining local stability of nonlinear systems, Lyapunov indirect method based on the Jacobian linearization has been widely employed due to its simplicity. Despite such popularity, it is not applicable to singular systems whose Jacobian has at least one eigenvalue that is equal to zero. In such singular cases, an appropriate Lyapunov function should be sought to determine the stability of systems, which is rather difficult and quite involved. In this paper, we seek for a simple criterion to determine stability of the equilibrium that is located at the boundary of the positive orthant, when one of eigenvalues of the Jacobian is zero. The goal of the paper is to present a generalized condition for the equilibrium to attract all trajectories that starting from initial condition in the positive orthant and near the equilibrium. Unlike the Lyapunov direct method, the proposed method requires just a simple algebraic computation for checking the stability of the critical point. Our approach is applied to various biological systems to show the effectiveness of the proposed method.

• Wind and Structures
• /
• 제26권1호
• /
• pp.25-34
• /
• 2018

This paper presents a beam finite element model of a vibrate wind blade in large elastic deformation subjected to the aerodynamic, centrifugal, gyroscopic and gravity loads. The gyroscopic loads applied to the blade are induced by her simultaneous vibration and rotation. The proposed beam finite element model is based on a simplex interpolation method and it is mainly intended to the numerical analysis of wind blades vibration in large elastic deformation. For this purpose, the theory of the sheared beams and the finite element method are combined to develop the algebraic equations system governing the three-dimensional motion of blade vibration. The applicability of the theoretical approach is elucidated through an original case study. Also, the static deformation of the used wind blade is assessed by appropriate software using a solid finite element model in order to show the effectiveness of the obtained results. To simulate the nonlinear dynamic response of wind blade, the predictor-corrector Newmark scheme is applied and the stability of numerical process is approved during a large time of blade functioning. Finally, the influence of the modified geometrical stiffness on the amplitudes and frequencies of the wind blade vibration induced by the sinusoidal excitation of gravity is analyzed.

• 대한기계학회논문집A
• /
• 제26권10호
• /
• pp.2104-2113
• /
• 2002

A finite element model is developed for the process of squeeze casting of metal matrix composites (MMCs) in cylindrical molds. The fluid flow and the heat transit. are fundamental phenomena in squeeze casting. To describe heat transfer in the solidification of molten aluminum, the energy equation is written in terms of temperature and enthalpy are applied in an axisymmetric model which is similar to the experimental system. A one dimensional flow model simulates the transient metal flow. A direct iteration technique was used to solve the resulting nonlinear algebraic equations, using a computer program to calculate the enthalpy, temperature and fluid velocity. The cooling curves and temperature distribution during infiltration and solidification were calculated fer pure aluminum. Experimentally, the temperature was measured and recorded using thermocouple wire. The measured time-temperature data were compared with the calculated cooling curves. The resulting agreement shows that the finite element model can accurately estimate the solidification time and predict the cooling process.