• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.5
• /
• pp.899-909
• /
• 1992

An iterative solution technique is presented to analyze the dynamic systems of rigid bodies subjected to kinematic constraints. Lagrange multipliers associated with the constraints are iteratively computed by monotonically reducing an appropriately defined constraint error vector, and the resulting equation of motion is solved by a well-established ODE technique. Constraints on the velocity and acceleration as well as the position are made to be satisfied at joints at each time step. Time integration is efficiently performed because decomposition or orthonormalization of the large matrix is not required at all. An acceleration technique is suggested for the faster convergence of the iterative scheme.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.57 no.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.

• Journal of the Korean Institute of Telematics and Electronics C
• /
• v.35C no.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
• /
• v.18 no.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
• /
• v.20 no.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
• /
• v.45 no.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.06a
• /
• 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
• /
• v.65 no.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.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.4
• /
• pp.583-594
• /
• 1989

The strong nonliner dynamic behavior of mechanical systems in the presence of clearances are studied. The nonlinearity is induced from the assumed symmetric piecewise-linear characteristics for stiffness and damping by the contact and uncontact. Based on Stoker's assertion concering the reasoning beyond the occurrence of subharmonics, the nonlinear differential equation is converted to four nonlinear algebraic equations form the boundary conditions at the contact points. For a single contact per half exciting period, under the assumption of symmetric response, the steady-state solutions obtained are in agreement with those of numerical integration. Also a nondimen-sionalized formulation is made for the purpose of parametric studies.