• Structural Engineering and Mechanics
• /
• v.52 no.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.

• Journal of applied mathematics & informatics
• /
• v.11 no.1_2
• /
• pp.215-241
• /
• 2003

In this paper, solute transport in heterogeneous aquifers using a modified Fokker-Planck equation (MFPE) is investigated. This newly developed mathematical model is characterised with a time-, scale-dependent dispersivity. A two-dimensional finite volume quadrilateral mesh method (FVQMM) based on a quadrilateral background interpolation mesh is developed for analysing the model. The FVQMM transforms the coupled non-linear partial differential equations into a system of differential equations, which is solved using backward differentiation formulae of order one through five in order to advance the solution in time. Three examples are presented to demonstrate the model verification and utility. Henry's classic benchmark problem is used to show that the MFPE captures significant features of transport phenomena in heterogeneous porous media including enhanced transport of salt in the upper layer due to its parameters that represent the dependence of transport processes on scale and time. The time and scale effects are investigated. Numerical results are compared with published results on the some problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.1 s.256
• /
• pp.29-39
• /
• 2007

A second-moment closure is applied to the prediction of a homogeneous turbulent shear flow laden with mono-size particles. The closure is curried out based on a 'two-fluid' methodology in which both carrier and dispersed phases are considered in the Eulerian frame. To reduce the number of coupled differential equations to be solved, Reynolds stress transport equations and algebraic stress models are judiciously combined to obtain the Reynolds stress of carrier and dispersed phases in the mean momentum equation. That is, the Reynolds stress components for carrier and dispersed phases are solved by modelled transport equations, but the fluid-particle velocity covariance tensors are treated by the algebraic models. The present predictions for all the components of Reynolds stresses are compared to the DNS data. Reasonable agreements are observed in all the components, and the effects of the coupling of carrier and dispersed phases are properly captured in every aspects.

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

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.263-285
• /
• 2002

This paper presents a numerical procedure for solving initial-value problems using the special functions which belong to a class of Rvachev's basis functions $R_{bf}$ based on algebraic and trigonometric polynomials. Because of infinite derivability of these functions, derivatives of all orders, required by differential equation of the problem and initial conditions, are used directly in the numerical procedure. The accuracy and stability of the proposed numerical procedure are proved on an example of a single degree of freedom system. Critical time step was also determined. An algorithm for solving multiple degree of freedom systems by the collocation method was developed. Numerical results obtained by $R_{bf}$ functions are compared with exact solutions and results obtained by the most commonly used numerical procedures for solving initial-value problems.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.1769-1775
• /
• 2004

In presence of collision between two rigid bodies, they exhibit impulsive behavior to generate physically feasible state. When the frictional impulse is involved, collision resolution can not be easily made based on a simple Newton's law or Poisson's law, mainly due to possible change of collision mode during collision, For example, sliding may change to sticking, and then sliding resumes. We first examine two conventional methods: the method of mode evolution by differential equation, and the other by linear complementarity programming. Then, we propose a new method for mode evolution by solving only algebraic equations defining mode changes. Further, our method attains the original nonlinear impulse cone constraint. The numerical simulation will elucidate the advantage of the proposed method as an alternative to conventional ones.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.1-8
• /
• 2004

For the linearized differential algebraic equation of the nonlinear constrained system, exact initial values of the acceleration are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental quantities since the mass matrix contains the zero element in the diagonal. This fact makes the mass matrix impossible to be positive definite. To overcome this singularity phenomenon the mass matrix needs to be modified to allow the feasible application of predictor and corrector in the iterative computation. In this paper the proposed numerical algorithm based on the modified mass matrix combines the conventional implicit algorithm, Newton-Raphson method and Newmark method. The numerical example presents reliabilities for the proposed algorithm via comparisons of the 4th order Runge-kutta method. The proposed algorithm seems to be satisfactory even though the acceleration, Lagrange multiplier, and energy show unstable behaviour. Correspondingly, it provides one important clue to another algorithm for the enhancement of the numerical results.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.6
• /
• pp.28-35
• /
• 1998

Computer simulation technique has been applied on the various engineering fields to reduce cost and development period. On this paper, we introduce a crane analysis program. Using this program, we can predict reaction force of each part or supporting force of truck crane on a personal computer system with out exclusive knowledge of multi-body dynamics. In order to consider the effect of boom flexibility according to each working condition, flexible dynamic theory is applied to the program. Actual crane model is analyzed on special work condition using this program and the results are compared with those of rigid boom model.

• Journal of Mechanical Science and Technology
• /
• v.17 no.4
• /
• pp.538-544
• /
• 2003

Although many methodologies exist for determining the constrained equations of motion, most of these methods depend on numerical approaches such as the Lagrange multiplier's method expressed in differential/algebraic systems. In 1992, Udwadia and Kalaba proposed explicit equations of motion for constrained systems based on Gauss's principle and elementary linear algebra without any multipliers or complicated intermediate processes. The generalized inverse method was the first work to present explicit equations of motion for constrained systems. However, numerical integration results of the equation of motion gradually veer away from the constraint equations with time. Thus, an objective of this study is to provide a numerical integration scheme, which modifies the generalized inverse method to reduce the errors. The modified equations of motion for constrained systems include the position constraints of index 3 systems and their first derivatives with respect to time in addition to their second derivatives with respect to time. The effectiveness of the proposed method is illustrated by numerical examples.