• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.363-380
• /
• 2020

In this paper, a cylindrical shell is immersed in a non-viscous fluid using first order shell theory of Sander. These equations are partial differential equations which are solved by approximate technique. Robust and efficient techniques are favored to get precise results. Employment of the Rayleigh-Ritz procedure gives birth to the shell frequency equation. Use of acoustic wave equation is done to incorporate the sound pressure produced in a fluid. Hankel's functions of second kind designate the fluid influence. Mathematically the integral form of the Lagrange energy functional is converted into a set of three partial differential equations. Throughout the computation, simply supported edge condition is used. Expressions for modal displacement functions, the three unknown functions are supposed in such way that the axial, circumferential and time variables are separated by the product method. Comparison is made for empty and fluid-filled cylindrical shell with circumferential wave number, length- and height-radius ratios, it is found that the fluid-filled frequencies are lower than that of without fluid. To generate the fundamental natural frequencies and for better accuracy and effectiveness, the computer software MATLAB is used.

• Journal of the Korean Society of Food Science and Nutrition
• /
• v.24 no.3
• /
• pp.384-388
• /
• 1995

High content of polyunsaturated fatty acid in fish oil makes it very susceptible to oxidation, which prevent fish oil from successful application to food processing or functional foods. To resolve this problem, oxidation stability model of fish oil was developed using the following differential equation : $dp/dt=k{\cdot}p(t){\cdot}[P_{max}\;-\;p(t)]$. This differential equation can be intergrated using analytical techniques to give : $p(t)=P_{max}/[1\;+\;[(P_{max}/P_{(0)})\;-\;-1]{\cdot}EXP(-K_p{\cdot}t)]$. At 50, 60, 70 and $80^{\circ}C,\;K_p$ were 0.00535, 0.01345, 0.02516 and 0.04675, respectively. The proposed model was well agreed with the measured data except for some minor deviations. In addition, $K_p$ was expressed as a function of temperature : $K_p=(1/P_{max})EXP\;[1\;-\;(8148/T)+20.1]$. Where T is absolute temperature($^{o}K$).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.2A
• /
• pp.155-162
• /
• 2009

This paper deals with the strongest beams with the solid regular polygon cross-section, whose volumes are always held constant. The differential equation of the elastic deflection curve of such beam subjected to the concentrated and trapezoidal distributed loads are derived and solved numerically. The Runge-Kutta method and shooting method are used to integrate the differential equation and to determine the unknown initial boundary condition of the given beam. In the numerical examples, the simple beams are considered as the end constraint and also, the linear, parabolic and sinusoidal tapers are considered as the shape function of cross sectional depth. As the numerical results, the configurations, i.e. section ratios, of the strongest beams are determined by reading the section ratios from the numerical data related with the static behaviors, under which static maximum behaviors become to be minimum.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3A
• /
• pp.251-258
• /
• 2009

This paper deals with the strongest beams with the solid regular polygon cross-section, whose volumes are always held constant. The differential equation of the elastic deflection curve of such beam subjected to the concentrated and trapezoidal distributed loads are derived and solved by using the double integration method. The Simpson's formula was used to numerically integrate the differential equation. In the numerical examples, the clamped-clamped and clamped-hinged ends are considered as the end constraints and the linear, parabolic and sinusoidal tapers are considered as the shape function of cross sectional depth. As the numerical results, the configurations, i.e. section ratios, of the strongest beams are determined by reading the section ratios from the numerical data obtained in this study, under which static maximum behaviors become to be minimum.

• The KIPS Transactions:PartB
• /
• v.13B no.1 s.104
• /
• pp.27-34
• /
• 2006

In this study, we propose this new algorithm that generates score function in ICA(Independent Component Analysis) using entropy theory. To generate score function, estimation of probability density function about original signals are certainly necessary and density function should be differentiated. Therefore, we used kernel density estimation method in order to derive differential equation of score function by original signal. After changing formula to convolution form to increase speed of density estimation, we used FFT algorithm that can calculate convolution faster. Proposed score function generation method reduces the errors, it is density difference of recovered signals and originals signals. In the result of computer simulation, we estimate density function more similar to original signals compared with Extended Infomax and Fixed Point ICA in blind source separation problem and get improved performance at the SNR(Signal to Noise Ratio) between recovered signals and original signal.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.1
• /
• pp.9-16
• /
• 2020

The objective of this study is to propose closed-form solutions to the free vibration response of single-degree-of-freedom (SDOF) systems, as part of fundamental research on dynamic systems with Coulomb friction. The motion of a dynamic system with Coulomb friction is described by a nonlinear differential equation, and, due to the variation in the sign of friction force term with the direction of motion, it is difficult to obtain the closed-form solution. To solve this problem, the nonlinear differential equation is directly computed by numerical integration, or an approximated solution is indirectly obtained using a linear differential equation wherein the damping effect due to Coulomb friction is replaced by an equivalent viscous damping term. However, these conventional methods do not provide a closed-form solution from a mathematical point of view. In this regard, closed-form solutions to the free vibration response of SDOF systems with Coulomb friction are derived herein by considering that the sign of the friction force term is reversed in each half-cycle of motion and by expanding it to the entire time history using the power series function. In addition, for a given initial condition, both the number of free vibration half-cycles and the response at the instant when free vibration motion stops are predicted under the condition that the motion of free vibration is stopped when the amplitude of the friction force is higher than that of the restoring force due to stiffness.

• Journal of Korea Water Resources Association
• /
• v.39 no.7 s.168
• /
• pp.593-603
• /
• 2006

The basin response to storm is regarded as nonlinearity inherently. In addition, the consistent nonlinearity of hydrologic system response to rainfall has been very tough and cumbersome to be treated analytically. The thing is that such nonlinear models have been avoided because of computational difficulties in identifying the model parameters from recorded data. The parameters of nonlinear system considered as dynamic effects in the conceptual model are optimized as the sum of errors between the observed and computed runoff is minimized. For obtaining the optimal parameters of functions, the historical data for the Bocheong watershed in the Geum river basin were tested by applying the numerical methods, such as quasi-linearization technique, Runge-Kutta procedure, and pattern-search method. The estimated runoff carried through from the storage function with dynamic effects was compared with the one of 1st-order differential equation model expressing just nonlinearity, and also done with Nash model. It was found that the 2nd-order model yields a better prediction of the hydrograph from each storm than the 1st-order model. However, the 2nd-order model was shown to be equivalent to Nash model when it comes to results. As a result, the parameters of nonlinear 2nd-order differential equation model performed from the present study provided not only a considerable physical meaning but also a applicability to Korean watersheds.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1072-1078
• /
• 2002

The constrained motion requires the determination of constraint force acting on unconstrained systems for satisfying given constraints. Most of the methods to decide the force depend on numerical approaches such that the Lagrange multiplier method, and the other methods need vector analysis or complicated intermediate process. In 1992, Udwadia and Kalaba presented the generalized inverse method to describe the constrained motion as well as to calculate the constraint force. The generalized inverse method has the advantages which do not require any linearization process for the control of nonlinear systems and can explicitly describe the motion of holonomically and/or nongolonomically constrained systems. In this paper, an explicit equation to describe the constrained motion is derived by minimizing the performance index, which is a function of constraint force vector, with respect to the constraint force. At this time, it is shown that the positive-definite weighting matrix in the performance index must be the inverse of mass matrix on the basis of the Gauss's principle and the derived differential equation coincides with the generalized inverse method. The effectiveness of this method is illustrated by means of two numerical applications.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.1A
• /
• pp.1-7
• /
• 2011

This paper deals with buckling loads of the column with the constant surface area. The shape function of variable column depth is chosen as the linear taper. The ordinary differential equation governing buckled shapes of the column is derived based on the dynamic equilibrium equation of such column subjected to an axial load. Three kinds of end constraint of hinged-hinged, hinged-clamped and clamped-clamped are considered in numerical examples. Effects of the column parameters on buckling loads are extensively discussed. Especially, section ratios of the strongest column are calculated, under which the maximum, i.e. strongest, buckling loads are achieved. Also the buckled shapes are obtained for searching the nodal points where the inner transverse supports are simply installed to increase the buckling loads.