• 대한기계학회논문집
• /
• 제16권5호
• /
• pp.899-909
• /
• 1992

본 연구에서는 대수 미분 방정식을 풀기위한 새로운 방법을 소개한다. 본 작업에서는 Lagrange multiplier의 값이 사전에 주어졌다고 생각하여, 즉 대수 미분 방정식을 순수한 상미분 방정식으로 변환하여, 잘 알려진 시간 적분법을 적용한다. 또 정확한 Lagrange Multiplier값은 반복 계산법(iterative scheme)에 의하여 계산한 다. 시간 적분의 정확도와 제한 조건의 정확도는 모두 보장된다. 특히 제한 조건 의 경우, 위치, 속도 및 가속도의 제한 조건이 모두 만족된다. 또 정확한 Lagrange multiplier의 값을 계산 가속기법(acceleration technique)에 의하여 대단히 빨리 계 산한다. 독립 좌표를 구할 필요가 없으므로 거대한 행열을 decomposition하는 등의 복잡한 절차가 불필요하며 N-R 반복법 역시 불필요하다. 이러한 사항들 및 Jacobian 행열의 sparsity로 인하여 경제적인 계산이 가능하게 된다.

• 전기학회논문지
• /
• 제57권6호
• /
• pp.915-923
• /
• 2008

This paper focuses on optimal tuning of nonlinear parameters of a dual-input power system stabilizer(dual-input PSS), which can improve the system damping performance immediately following a large disturbance. Until recently, various PSS models have developed to bring stability and reliability to power systems, and some of these models are used in industry applications. However, due to non-smooth nonlinearities from the interaction between linear parameters(gains and time constants of linear controllers) and nonlinear parameters(saturation output limits), the output limit parameters cannot be determined by the conventional tuning methods based on linear analysis. Only ad hoc tuning procedures('trial and error' approach) have been used. Therefore, the steepest descent method is applied to implement the optimal tuning of the nonlinear parameters of the dual-input PSS. The gradient required in this optimization technique can be computed from trajectory sensitivities in hybrid system modeling with the differential-algebraic-impulsive-switched(DAIS) structure. The optimal output limits of the dual-input PSS are evaluated by time-domain simulation in both a single machine infinite bus(SMIB) system and a multi-machine power system in comparison with those of a single-input PSS.

• 전자공학회논문지C
• /
• 제35C권4호
• /
• pp.38-49
• /
• 1998

A fossil power plant can be modeled by a lot of algebraic equations and differential equations. When we simulate a large, complicated fossil power plant by a computer such as workstation or PC, it takes much time until overall equations are completely calculated. Therefore, new processing systems which have high computing speed is ultimately needed for real-time or high-speed(faster than real-time) simulators. This paper presents an enhanced strategy in which high computing power can be provided by parallel processing of DSP processors with communication links. DSP system is designed for general purpose. Parallel DSP system can be easily expanded by just connecting new DSP modules to the system. General urpose DSP modules and a VME interface module was developed. New model and techniques for the task allocation are also presented which take into account the special characteristics of parallel I/O and computation. As a realistic cost function of task allocation, we suggested 'simulation period' which represents the period of simulation output intervals. Based on the development of parallel DSP system and realistic task allocation techniques, we cound achieve good efficiency of parallel processing and faster simulation speed than real-time.

• Smart Structures and Systems
• /
• 제18권6호
• /
• pp.1125-1143
• /
• 2016

Forced vibration analysis of a simple supported viscoelastic nanobeam is studied based on modified couple stress theory (MCST). The nanobeam is excited by a transverse triangular force impulse modulated by a harmonic motion. The elastic medium is considered as Winkler-Pasternak elastic foundation.The damping effect is considered by using the Kelvin-Voigt viscoelastic model. The inclusion of an additional material parameter enables the new beam model to capture the size effect. The new non-classical beam model reduces to the classical beam model when the length scale parameter is set to zero. The considered problem is investigated within the Timoshenko beam theory by using finite element method. The effects of the transverse shear deformation and rotary inertia are included according to the Timoshenko beam theory. The obtained system of differential equations is reduced to a linear algebraic equation system and solved in the time domain by using Newmark average acceleration method. Numerical results are presented to investigate the influences the material length scale parameter, the parameter of the elastic medium and aspect ratio on the dynamic response of the nanobeam. Also, the difference between the classical beam theory (CBT) and modified couple stress theory is investigated for forced vibration responses of nanobeams.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Structural Engineering and Mechanics
• /
• 제45권1호
• /
• pp.95-110
• /
• 2013

In the recent decade, practical of wavelet technique is being utilized in various domain of science. Particularly, engineers are interested to the wavelet solution method in the time series analysis. Fundamentally, seismic responses of structures against time history loading such as an earthquake, illustrates optimum capability of systems. In this paper, a procedure using particularly discrete Haar wavelet basis functions is introduced, to solve dynamic equation of motion. In the proposed approach, a straightforward formulation in a fluent manner is derived from the approximation of the displacements. For this purpose, Haar operational matrix is derived and applied in the dynamic analysis. It's free-scaled matrix converts differential equation of motion to the algebraic equations. It is shown that accuracy of dynamic responses relies on, access of load in the first step, before piecewise analysis added to the technique of equation solver in the last step for large scale of wavelet. To demonstrate the effectiveness of this scheme, improved formulations are extended to the linear and nonlinear structural dynamic analysis. The validity and effectiveness of the developed method is verified with three examples. The results were compared with those from the numerical methods such as Duhamel integration, Runge-Kutta and Wilson-${\theta}$ method.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 2005년도 ICCAS
• /
• pp.2410-2413
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAB, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. In the pre-processor of 3D FEM package, a new 3D mesh generating algorithm can make information on 3D Delaunay tetrahedral mesh elements for analyses of 3D FEM problems. The solver of the 3D FEM package offers three methods for solving the linear algebraic matrix equation, i.e., Gauss-Jordan elimination solver, Band solver, and Skyline solver. The post-processor visualizes the results for 3D FEM problems such as the deformed position and the stress. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.

• Structural Engineering and Mechanics
• /
• 제65권1호
• /
• pp.111-120
• /
• 2018

This paper presents a new and simple solution for determining the natural frequencies of framed tube combined with shear-walls and tube-in-tube systems. The novelty of the presented approach is based on the bending moment function approximation instead of the mode shape function approximation. This novelty makes the presented solution very simpler and very shorter in the mathematical calculations process. The shear stiffness, flexural stiffness and mass per unit length of the structure are variable along the height. The effect of the structure weight on its natural frequencies is considered using a variable axial force. The effects of shear lag phenomena has been investigated on the natural frequencies of the structure. The whole structure is modeled by an equivalent non-prismatic shear-flexural cantilever beam under variable axial forces. The governing differential equation of motion is converted into a system of linear algebraic equations and the natural frequencies are calculated by determining a non-trivial solution for the system of equations. The accuracy of the proposed method is verified through several numerical examples and the results are compared with the literature.

• 대한기계학회논문집
• /
• 제16권12호
• /
• pp.2298-2305
• /
• 1992

본 연구에서는 보의 종방향 강성치를 실제보다 훨씬 작은 횡방향 강성치 수준 으로 낮추어 시간 적분을 수행하는 대신, 종방향에 작용하는 힘을 분산시켜서 계산시 종방향 상대 변형을 억제하는 기법을 사용하였다. 또 여러개의 물체가 조인트에서 서로 연결되어 이루어진 다물체 도역학의 해는 운동 방정식과 기구학적 제한조건을 모 두 만족시켜야 하며, 상미분 방정식 해법용 시간 적분을 바로 적용할 수는 없다. 본 논문에서는 참고문헌(13)과 (14)에 소개된 반복계산법 및 계산가속기법에 의하여 Lag- range multiplier를 운동 방정식에서 소거한 후 시간 적분을 수행하는 방법을 사용한 다.

• 대한기계학회논문집
• /
• 제13권4호
• /
• pp.583-594
• /
• 1989

본 연구는 기계적 간극의 존재로 양면 접촉에 의한 강성 및 감쇠가 대칭적 으로 두번 생기는 강제 대칭 편적 선형(forced symmetric piecewise-linear vibarti- on)에 대하여 비선형진동의 분수조화 정상해 과정을 설명한 Stoker 가설로부터 분수조화진동을 포함한 복수의 정상해를 구하고 피로, 마멸 연구를 위한 접촉시간 및 접촉력을 계산하였다. 또한 강제 대칭 편적 선형진동자의 운동방정식을 무차원화 하고 각 무차원 시스템 변수에 따른 주파수응답특성을 살펴보았다.