• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.119-139
• /
• 2014

A new method is proposed for random vibration anaylsis of hysteretic systems subjected to non-stationary random excitations. With the Bouc-Wen model, motion equations of hysteretic systems are first transformed into quasi-linear equations by applying the concept of equivalent excitations and decoupling of the real and hysteretic displacements, and the derived equation system can be solved by either the precise time integration or the Newmark-${\beta}$ integration method. Combining the numerical solution of the auxiliary differential equation for hysteretic displacements, an explicit iteration algorithm is then developed for the dynamic response analysis of hysteretic systems. Because the computational cost for a large number of deterministic analyses of hysteretic systems can be significantly reduced, Monte-Carlo simulation using the explicit iteration algorithm is now viable, and statistical characteristics of the non-stationary random responses of a hysteretic system can be obtained. Numerical examples are presented to show the accuracy and efficiency of the present approach.

• Journal of the Korean Society for Precision Engineering
• /
• v.14 no.3
• /
• pp.66-74
• /
• 1997

In a high speed and high precision vertical machining center, chatter vibration is easily generated due to unbalanced masses in rotating parts and changtes of cutting forces. In this paper, modal test is performed to obtain modal parameters of the vertical machining center. In order to predit the cutting force of endmilling process for various cutting conditions, a mathematical model is given and this model is based on chip load, cutting geometry, and relationship between cutting forces and the chip load. Specific cutting constants of the model are obtained by averaging forces of cutting tests. The interactions between the dy- namic characteristics and cutting dynamics of the vertical machining center make the primary and the secondary feedback loops, and we make use of the equations of system to predict the chatter vibration. The chatter prediction is formulated as linear differential-differene equations, and simulated for several cases. Trends of vibration as radial and axial depths of cut are changed are shown and compared.

• Structural Engineering and Mechanics
• /
• v.35 no.2
• /
• pp.141-173
• /
• 2010

In this paper a boundary element method is developed for the general flexural-torsional buckling analysis of Timoshenko beams of arbitrarily shaped cross section. The beam is subjected to a compressive centrally applied concentrated axial load together with arbitrarily axial, transverse and torsional distributed loading, while its edges are restrained by the most general linear boundary conditions. The resulting boundary value problem, described by three coupled ordinary differential equations, is solved employing a boundary integral equation approach. All basic equations are formulated with respect to the principal shear axes coordinate system, which does not coincide with the principal bending one in a nonsymmetric cross section. To account for shear deformations, the concept of shear deformation coefficients is used. Six coupled boundary value problems are formulated with respect to the transverse displacements, to the angle of twist, to the primary warping function and to two stress functions and solved using the Analog Equation Method, a BEM based method. Several beams are analysed to illustrate the method and demonstrate its efficiency and wherever possible its accuracy. The range of applicability of the thin-walled theory and the significant influence of the boundary conditions and the shear deformation effect on the buckling load are investigated through examples with great practical interest.

• Steel and Composite Structures
• /
• v.28 no.3
• /
• pp.349-362
• /
• 2018

This paper presents a semi-analytical solution for the creep analysis and life assessment of 304L austenitic stainless steel thick truncated conical shells using multilayered method based on the first order shear deformation theory (FSDT). The cone is subjected to the non-uniform internal pressure and temperature gradient. Damages are obtained in thick truncated conical shell using Robinson's linear life fraction damage rule, and time to rupture and remaining life assessment is determined by Larson-Miller Parameter (LMP). The creep response of the material is described by Norton's law. In the multilayer method, the truncated cone is divided into n homogeneous disks, and n sets of differential equations with constant coefficients. This set of equations is solved analytically by applying boundary and continuity conditions between the layers. The results obtained analytically have been compared with the numerical results of the finite element method. The results show that the multilayered method based on FSDT has an acceptable amount of accuracy when one wants to obtain radial displacement, radial, circumferential and shear stresses. It is shown that non-uniform pressure has significant influences on the creep damages and remaining life of the truncated cone.

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

• Journal of the Korean Society of Food Science and Nutrition
• /
• v.28 no.1
• /
• pp.178-190
• /
• 1999

This study was conducted to develop a simulation model for the growth dynamics of pigs and to describe quantitatively protein deposition depending on the amino acid composition of feed protein. In the model it is assumed that the essential processes that determine the utilization of feed protein in the whole body are protein synthesis, breakdown of protein, and oxidation of amino acid. Besides, it is also assumed that occurrence of protein deposition depends on genetic potential and amino acid composition of feed protein. The genetic potential for the protein deposition is the maximum capacity of protein synthesis, being dependent on the protein mass of the whole body. To describe the effect of amino acid composition of feed on the protein deposition, a factor, which consist of ten amino acid functions and lie between 0 and 1, is introduced. Accordingly a model was developed, which is described with 15 flux equations and 11 differential equations and is composed of two compartments. The model describes non linear structure of the protein utilization system of an organism, which is in non steady state. The objective function for the simulation was protein deposition(g/day) cal culated according to the empirical model, PAF(product of amino acid functions) of Menke. The mean of relative difference between the simulated protein deposition and PAF calculated values, lied in a range of 11.8%. The simulated protein synthesis and breakdown rates(g/day) in the whole body showed a parallel behavior in the course of growth.

• Nuclear Engineering and Technology
• /
• v.27 no.6
• /
• pp.845-852
• /
• 1995

We investigate the concurrent solution of differential equations by the waveform relaxation (WR) method, an iterative method for analyzing linear and nonlinear dynamical systems in the time do-main. The method, at each iteration, decomposes the dynamical system into several subsystems, each of which is analyzed for the entire given time interval. The method, when efficiently implemented, results in algorithms with a highly parallelizable concurrent fraction. In this paper, the waveform relaxation method is introduced and applied to two types of reactor dynamics problems. It is concluded that the U method can be applied to reactor dynamics equations, but that its parallel performance on the KMRR dynamics is only modest.

• Structural Engineering and Mechanics
• /
• v.28 no.6
• /
• pp.749-765
• /
• 2008

Free axisymmetric vibrations of layered cylindrical shells of variable thickness are studied using spline function approximation techniques. Three different types of thickness variations are considered namely linear, exponential and sinusoidal. The equations of axisymmetric motion of layered cylindrical shells, on the longitudinal and transverse displacement components are obtained using Love's first approximation theory. A system of coupled differential equations on displacement functions are obtained by assuming the displacements in a separable form. Then the displacements are approximated using Bickley-spline approximation. The vibrations of two-layered cylindrical shells, made up of several types of layered materials and different boundary conditions are considered. Parametric studies have been made on the variation of frequency parameter with respect to the relative layer thickness, length ratio and type of thickness variation parameter.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.2
• /
• pp.63-69
• /
• 2001

In this study, we shall be concerned with computing in parallel a few of the smallest eigenvalues and their corresponding eigenvectors of the eigenvalue problem, $Ax={\lambda}Bx$, where A is symmetric, and B is symmetric positive definite. Both A and B are large and sparse. Recently iterative algorithms based on the optimization of the Rayleigh quotient have been developed, and CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for small extreme eigenvalues. As in the case of a system of linear equations, successful application of the CG scheme to eigenproblems depends also upon the preconditioning techniques. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. The idea underlying the present work is a parallel computation of the Multi-Color Block SSOR preconditioning for the CG optimization of the Rayleigh quotient together with deflation techniques. Multi-Coloring is a simple technique to obatin the parallelism of order n, where n is the dimension of the matrix. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI(Message Passing Interface) library was adopted for the interprocessor communications. The test problems were drawn from the discretizations of partial differential equations by finite difference methods.

• Kyungpook Mathematical Journal
• /
• v.61 no.4
• /
• pp.831-843
• /
• 2021

In this paper we look at the three dimensional stagnation point flow problem over a rough rotating disk. We study the theoretical behaviour of the stagnation point flow, or forced flow, in the presence of a slip factor in which convective instability stationary modes appear. We make a numerical investigation of the effects of slip on the behaviour of the flow components of the stagnation point flow where the disk is rough. We provide, for the first time in the literature, a complete convective instability analysis and an energy analysis. Suitable similarity transformations are used to reduce the Navier-Stokes equations and the continuity equation into a system of highly non-linear coupled ordinary differential equations, and these are solved numerically subject to suitable boundary conditions using the bvp4c function of MATLAB. The convective instability analysis and the energy analysis are performed using the Chebyshev spectral method in order to obtain the neutral curves and the energy bars. We observe that the roughness of the disk has a destabilising effect on both Type-I and Type-II instability modes. The results obtained will be prominently treated as benchmarks for our future studies on stagnation flow.