• Computers and Concrete
• /
• v.7 no.4
• /
• pp.317-329
• /
• 2010

The paper describes localization of deformation in a bar under tensile loading. The material of the bar is considered as non-linear viscous elastic and the bar consists of two symmetric halves. It is assumed that the model represents behavior of the quasi-brittle viscous material under uniaxial tension with different loading rates. Besides that, the bar could represent uniaxial stress-strain law on a single plane of a microplane material model. Non-linear material property is taken from the microplane material model and it is coupled with the viscous damper producing non-linear Maxwell material model. Mathematically, the problem is described with a system of two partial differential equations with a non-linear algebraic constraint. In order to obtain solution, the system of differential algebraic equations is transformed into a system of three partial differential equations. System is subjected to loadings of different rate and it is shown that localization occurs only for high loading rates. Mathematically, in such a case two solutions are possible: one without the localization (unstable) and one with the localization (stable one). Furthermore, mass is added to the bar and in that case the problem is described with a system of four differential equations. It is demonstrated that for high enough loading rates, it is the added mass that dominates the response, in contrast to the viscous and elastic material parameters that dominated in the case without mass. This is demonstrated by several numerical examples.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.122-128
• /
• 1989

In this paper we propose a new method for solving a set of nonlinear algebraic equations encountered in the DC and transient analyses of electronic circuits. This method will be called Quadratic Newton-Raphson Method (QNRM), since it is based on the Newton-Raphson Method (NRM) but effectively takes into accoujnt the second order derivative terms in the Taylor series expansion of the nonlinear algebraic equations. The second order terms are approximated by linear terms using a carefully estimated solution at each iteration. Preliminary simulation results show that the QNRM saves the overall computational time significantly in the DC and transient analysis, compared with the conventional NRM.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.12
• /
• pp.2647-2655
• /
• 2002

This research presents an analytical model to investigate the stability due to the ball bearing waviness i n a rotating system supported by two ball bearings. The stiffness of a ball bearing changes periodically due to the waviness in the rolling elements as the rotor rotates, and it can be calculated by differentiating the nonlinear contact forces. The linearized equations of motion can be represented as a parametrically excited system in the form of Mathieu's equation, because the stiffness coefficients have time -varying components due to the waviness. Their solution can be assumed as a Fourier series expansion so that the equations of motion can be rewritten as the simultaneous algebraic equations with respect to the Fourier coefficients. Then, stability can be determined by solving the Hill's infinite determinant of these algebraic equations. The validity of this research is proved by comparing the stability chart with the time responses of the vibration model suggested by prior researches. This research shows that the waviness in the rolling elements of a ball bearing generates the time-varying component of the stiffness coefficient, whose frequency is called the frequency of the parametric excitation. It also shows that the instability takes place from the positions in which the ratio of the natural frequency to the frequency of the parametric excitation corresponds to i/2 (i=1,2,3..).

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.166-174
• /
• 2002

This paper presents an analytical method to Investigate the stability of a hydrodynamic journal bearing with rotating herringbone grooves. The dynamic coefficients of the hydrodynamic journal bearing are calculated using the FEM and the perturbation method. The linear equations of motion can be represented as a parametrically excited system because the dynamic coefficients have time-varying components due to the rotating grooves, even in the steady state. Their solution can be assumed as a Fourier series expansion so that the equations of motion can be rewritten as simultaneous algebraic equations with respect to the Fourier coefficients. Then, stability can be determined by solving Hill's infinite determinant of these algebraic equations. The validity of this research is proved by the comparison of the stability chart with the time response of the whirl radius obtained from the equations of motion. This research shows that the instability of the hydrodynamic journal bearing with rotating herringbone grooves increases with increasing eccentricity and with decreasing groove number, which play the major roles in increasing the average and variation of stiffness coefficients, respectively. It also shows that a high rotational speed is another source of instability by increasing the stiffness coefficients without changing the damping coefficients.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.181-189
• /
• 2002

This research presents an analytical model to investigate the stability due to the ball bearing waviness in a rotating system supported by two ball bearings. The stiffness of a ball bearing changes periodically due to the waviness in the rolling elements as the rotor rotates, and it can be calculated by differentiating the nonlinear contact forces. The linearized equations of motion can be represented as a parametrically excited system in the form of Mathieu's equation, because the stiffness coefficients have time-varying components due to the waviness. Their solution can be assumed as a Fourier series expansion so that the equations of motion can be rewritten as the simultaneous algebraic equations with respect to the Fourier coefficients. Then, stability can be determined by solving the Hill's infinite determinant of these algebraic equations. The validity of this research is proved by comparing the stability chart with the time responses of the vibration model suggested by prior researches. This research shows that the waviness in the rolling elements of a ball bearing generates the time-varying component of the stiffness coefficient, whose frequency is called the frequency of the parametric excitation. It also shows that the instability takes place from the positions in which the ratio of the natural frequency to the frequency of the parametric excitation corresponds to i/2 (i= 1,2,3..).

• Journal of Institute of Control, Robotics and Systems
• /
• v.21 no.7
• /
• pp.648-653
• /
• 2015

The kinematics with the instantaneous motion and statics of a manipulator has generally been proven algebraically. The algebraic solutions give very simple and straightforward results but the solutions do not have any meaning in physics or geometry. Therefore it is not easy to extend the algebraic results to design or control a robotic manipulator efficiently. Recently, geometrical approach to define the instantaneous motion or static relation of a manipulator is popularly researched and the results have very strong advantages to have a physical insight in the solution. In this paper, the instantaneous motion and static relation of a planar manipulator are described by geometrical approach, specifically by an axis screw and a line screw. The mass center of a triangle with weight and a perpendicular distance between the two screws are useful geometric measures for geometric analysis. This study provides a geometric interpretation of the kinematics and statics of a planar manipulator, and the method can be applied to design or control procedure from the geometric information in the equations.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.11
• /
• pp.943-950
• /
• 2013

This paper discusses a harmonic response estimation method on the L$\acute{e}$vy plate with two opposite edges simply supported and the other two edges having free boundary conditions. Since the equation of motion of the plate is not self-adjoint, the modes are not orthogonal to each other on the domain. Noting that the L$\acute{e}$vy plate can be expressed using one term sinusoidal function that is orthogonal to other sinusoidal functions, this paper suggested the calculation method that is equivalent to finding a least square error minimization solution of the finite number of algebraic equations. Example problems subjected to a distributed area loading and a distributed line loading are defined and their solutions are provided. The solutions are compared to those of the commercial code, ANSYS. According to the verification results, it is expected that the suggested method will be useful to predict the forced response on the L$\acute{e}$vy plate with the distributed area or line loading conditions.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.11
• /
• pp.176-184
• /
• 2000

Recently, it becomes a very important issue to consider the mechanical systems such as high-speed vehicle and railway train moving on a flexible beam structure. Using general approach proposed in the first part of this paper, it tis possible to predict planar motion of constrained mechanical system and elastic structure with various kinds of foundation supporting condition. Combined differential-algebraic equations of motion derived from both multibody dynamics theory and Finite Element Method can be analyzed numerically using generalized coordinate partitioning algorithm. To verify the validity of this approach, results from simply supported elastic beam subjected to a moving load are compared with exact solution from a reference. Finally, parameter study is conducted for a moving vehicle model on a simply supported 3-span bridge.

• Communications for Statistical Applications and Methods
• /
• v.28 no.5
• /
• pp.477-491
• /
• 2021

The Markov modulated Brownian motion is a substantial generalization of the classical Brownian Motion. On the other hand, the Markovian arrival process (MAP) is a point process whose family is dense for any stochastic point process and is used to approximate complex stochastic counting processes. In this paper, we consider a superposition of the Markov modulated Brownian motion (MMBM) and the Markovian arrival process of jumps which are distributed as the bilateral ph-type distribution, the class of which is also dense in the space of distribution functions defined on the whole real line. In the model, we assume that the inter-arrival times of the MAP depend on the underlying Markov process of the MMBM. One of the subjects of this paper is introducing how to obtain the first passage probabilities of the superposed process using a stochastic doubling algorithm designed for getting the minimal solution of a nonsymmetric algebraic Riccatti equation. The other is to provide eigenvalue and eigenvector results on the superposed process to make it possible to apply the GTH-like algorithm, which improves the accuracy of the doubling algorithm.

• Structural Engineering and Mechanics
• /
• v.85 no.5
• /
• pp.621-633
• /
• 2023

The major objective of this paper is to study the receding contact problem between two functional graded layers under a flat indenter. The gravity is assumed negligible, and the shear moduli of both layers are assumed to vary exponentially along the thickness direction. In the absence of body forces, the problem is reduced to a system of Fredholm singular integral equations with the contact pressure and contact size as unknowns via Fourier integral transform, which is transformed into an algebraic one by the Gauss-Chebyshev quadratures and polynomials of both the first and second kinds. Then, an iterative speediest descending algorithm is proposed to numerically solve the system of algebraic equations. Both semi-analytical and finite element method, FEM solutions for the presented problem validate each other. To improve the accuracy of the numerical result of FEM, a graded FEM solution is performed to simulate the FGM mechanical characteristics. The results reveal the potential links between the contact stress/size and the indenter size, the thickness, as well as some other material properties of FGM.