• 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
• /
• 제52권5호
• /
• pp.997-1017
• /
• 2014

According to the structure characteristics of the non-uniform beam bridge, a practical model for calculating the vibration equation of the non-uniform beam bridge is given and the application scope of the model includes not only the beam bridge structure but also the non-uniform beam with added masses and elastic supports. Based on the Bernoulli-Euler beam theory, extending the application of the modal perturbation method and establishment of a semi-analytical method for solving the vibration equation of the non-uniform beam with added masses and elastic supports based is able to be made. In the modal subspace of the uniform beam with the elastic supports, the variable coefficient differential equation that describes the dynamic behavior of the non-uniform beam is converted to nonlinear algebraic equations. Extending the application of the modal perturbation method is suitable for solving the vibration equation of the simply supported and continuous non-uniform beam with its arbitrary added masses and elastic supports. The examples, that are analyzed, demonstrate the high precision and fast convergence speed of the method. Further study of the timesaving method for the dynamic characteristics of symmetrical beam and the symmetry of mode shape should be developed. Eventually, the effects of elastic supports and added masses on dynamic characteristics of the three-span non-uniform beam bridge are reported.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 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.

• 제어로봇시스템학회논문지
• /
• 제5권6호
• /
• pp.647-655
• /
• 1999

For linear time-varying systems described by the triple (A(t),B(t),C(t)) where A(t),B(t),C(t) are the system, the input, and the output matrices, respectively, we propose concepts for measures of modal and gross controllability /observability. We introduce a differential algebraic eigenbvalue theory for linear time-varying systems to calculate the PD-eigenvalues and left and right PD-eigenvectors of the system matrix A(t) which will be used to derive the concepts for the measures. The time-dependent angle between the left PD-eigenvectors of the system matrix A(t) and the columns of the input matrix B(t), and the magnitude of the each element of the input matrix B(t) are used to propose the modal controllability measure. Similarly, the time-dependent angle between the right PD-eigenvectors of the system matrix A(t) and the rows of the output matrix C(t) are used to propose the madal observability measure. Gross measure of controllability of a mode from all inputs and its gross measure of observability in all outputs for the linear time-varying systems are also proposed. Numerical examples are presented to illustrate the proposed concepts.

• Advances in Computational Design
• /
• 제5권1호
• /
• pp.55-89
• /
• 2020

For the first time, an accurate analytical solution, based on coupled three-dimensional (3D) piezoelasticity equations, is presented for free vibration analysis of the angle-ply elastic and piezoelectric flat laminated panels under arbitrary boundary conditions. The present analytical solution is applicable to composite, sandwich and hybrid panels having arbitrary angle-ply lay-up, material properties, and boundary conditions. The modified Hamiltons principle approach has been applied to derive the weak form of governing equations where stresses, displacements, electric potential, and electric displacement field variables are considered as primary variables. Thereafter, multi-term multi-field extended Kantorovich approach (MMEKM) is employed to transform the governing equation into two sets of algebraic-ordinary differential equations (ODEs), one along in-plane (x) and other along the thickness (z) direction, respectively. These ODEs are solved in closed-form manner, which ensures the same order of accuracy for all the variables (stresses, displacements, and electric variables) by satisfying the boundary and continuity equations in exact manners. A robust algorithm is developed for extracting the natural frequencies and mode shapes. The numerical results are reported for various configurations such as elastic panels, sandwich panels and piezoelectric panels under different sets of boundary conditions. The effect of ply-angle and thickness to span ratio (s) on the dynamic behavior of the panels are also investigated. The presented 3D analytical solution will be helpful in the assessment of various 1D theories and numerical methods.

• Structural Engineering and Mechanics
• /
• 제55권3호
• /
• pp.655-674
• /
• 2015

In this paper, a new and simplified method is presented in which the natural frequencies of the uniform and non-uniform beams are calculated through simple mathematical relationships. The various vibration problems such as: Rayleigh beam under variable axial force, axial vibration of a bar with and without end discrete spring, torsional vibration of a bar with an attached mass moment of inertia, flexural vibration of the beam with laterally distributed elastic springs and also flexural vibration of the beam with effects of viscose damping are investigated. The governing differential equations are first obtained and then; according to a harmonic vibration, are converted into single variable equations in terms of location. Through repetitive integrations, the governing equations are converted into weak form integral equations. The mode shape functions of the vibration are approximated using a power series. Substitution of the power series into the integral equations results in a system of linear algebraic equations. The natural frequencies are determined by calculation of a non-trivial solution for system of equations. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained with those obtained from other published references and results of available finite element software.

• Steel and Composite Structures
• /
• 제52권5호
• /
• pp.501-513
• /
• 2024

Two-directional functionally graded materials (2D-FGMs) show extraordinary physical properties which makes them ideal candidates for designing smart micro-switches. Pull-in instability is one of the most critical challenges in the design of electrostatically-actuated microswitches. The present research aims to bridge the gap in the static pull-in instability analysis of microswitches composed of 2D-FGM. Euler-Bernoulli beam theory with geometrical nonlinearity effect (i.e. von-Karman nonlinearity) in conjunction with the modified couple stress theory (MCST) are employed for mathematical formulation. The micro-switch is subjected to electrostatic actuation with fringing field effect and Casimir force. Hamilton's principle is utilized to derive the governing equations of the system and corresponding boundary conditions. Due to the extreme nonlinear coupling of the governing equations and boundary conditions as well as the existence of terms with variable coefficients, it was difficult to solve the obtained equations analytically. Therefore, differential quadrature method (DQM) is hired to discretize the obtained nonlinear coupled equations and non-classical boundary conditions. The result is a system of nonlinear coupled algebraic equations, which are solved via Newton-Raphson method. A parametric study is then implemented for clamped-clamped and cantilever switches to explore the static pull-in response of the system. The influences of the FG indexes in two directions, length scale parameter, and initial gap are discussed in detail.

• 전자공학회논문지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.

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

• Structural Engineering and Mechanics
• /
• 제43권5호
• /
• pp.561-582
• /
• 2012

Despite popularity of FEM in analysis of static and dynamic structural problems and the routine applicability of FE softwares, analytical methods based on simple mathematical relations is still largely sought by many researchers and practicing engineers around the world. Development of such analytical methods for analysis of free vibration of non-prismatic beams is also of primary concern. In this paper a new and simple method is proposed for determination of vibration frequencies of non-prismatic beams under variable axial forces. The governing differential equation is first obtained and, according to a harmonic vibration, is converted into a single variable equation in terms of location. Through repetitive integrations, integral equation for the weak form of governing equation is derived. The integration constants are determined using the boundary conditions applied to the problem. The mode shape functions are approximated by a power series. Substitution of the power series into the integral equation transforms it into a system of linear algebraic equations. Natural frequencies are determined using a non-trivial solution for system of equations. Presented method is formulated for beams having various end conditions and is extended for determination of the buckling load of non-prismatic beams. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained to those obtained using available finite element software.