• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.155-163
• /
• 2012

This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.151-166
• /
• 2022

Peltier cells have low efficiency, but they are becoming attractive alternatives for affordable and environmentally clean cooling. In this line, the current article develops closed-form and semianalytical solutions to improve the temperature distribution of Bi₂Te₃ thermoelements. From the distribution, the main objective of the current work-the optimal electric intensity to maximize cooling-is inferred. The general one-dimensional differential coupled equation is integrated for linear and quadratic geometry of thermoelements, under temperature constant properties. For a general shape, a piece-wise solution based on heat flux continuity among virtual layers gives accurate analytical solutions. For variable properties, another piece-wise solution is developed but solved iteratively. Taking advantage of the formulae, the optimal intensity is directly derived with a minimal computational cost; its value will be of utility for more advanced designs. Finally, a parametric study including straight, two linear, barrel, hourglass and vase geometries is presented, drawing conclusions on how the shape of the thermoelement affects the coupled phenomena. A specially developed coupled and non-linear finite element research code is run taking into account all the materials of the cell and using symmetries and repetitions. These accurate results are used to validate the analytical ones.

• Journal of applied mathematics & informatics
• /
• v.9 no.2
• /
• pp.561-577
• /
• 2002

In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.147-168
• /
• 2019

This work aims to simulate three-dimensional heavy oil flow in a reservoir with heater-wells. Mass, momentum and energy balances, as well as correlations for rock and fluid properties, are used to obtain non-linear partial differential equations for the fluid pressure and temperature, and for the rock temperature. Heat transfer is simulated using a two-equation model that is more appropriate when fluid and rock have very different thermal properties, and we also perform comparisons between one- and two-equation models. The governing equations are discretized using the Finite Volume Method. For the numerical solution, we apply a linearization and an operator splitting. As a consequence, three algebraic subsystems of linearized equations are solved using the Conjugate Gradient Method. The results obtained show the suitability of the numerical method and the technical feasibility of heating the reservoir with static equipment.

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.54-63
• /
• 2002

In recent years, mechanical systems such as high speed vehicles and railway trains moving on elastic beam structures have become a very important issue to consider. In this paper, a general approach, which can predict the dynamic behavior of a constrained mechanical system moving on a flexible beam structure, is proposed. Various supporting conditions for the foundation support are considered for the elastic beam structure. The elastic structure is assumed to be a non-uniform and linear Bernoulli-Euler beam with a proportional damping effect. Combined differential-algebraic equation of motion is derived using the multi-body dynamics theory and the finite element method. The proposed equations of motion can be solved numerically using the generalized coordinate partitioning method and predictor-corrector algorithm, which is an implicit multi-step integration method.

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

• International Journal of Fuzzy Logic and Intelligent Systems
• /
• v.3 no.1
• /
• pp.13-17
• /
• 2003

A continuous solution of the Dirichlet boundary value problem for the heat equation $u_t$＝$a2u_{xx}$ using a fuzzy system is described. We first apply the Crank-Nicolson method to obtain a discrete solution at the grid points for the heat equation. Then we find a continuous function to represent approximately the discrete values at the grid points in the form of a bicubic spline function (equation omitted) that can in turn be represented exactly by a fuzzy system. We show that the computed values at non-grid points using the bicubic spline function is much smaller than the ones obtained by linear interpolations of the values at the grid points. We also show that the fuzzy rule table in the fuzzy system representation of the bicubic spline function can be viewed as a gray scale image. Hence, the fuzzy rules provide a visual representation of the functions of two variables where the contours of different levels for the function are shown in different gray scale levels

• Journal of Institute of Control, Robotics and Systems
• /
• v.9 no.1
• /
• pp.1-9
• /
• 2003

An active control of the lateral vibration of a translating tensioned Euler-Bemoulli beam is investigated. The dynamics of the translating tensioned beam is represented by a non-linear hyperbolic partial differential equation. A right boundary control law based upon the Lyapunov's second method is derived. The transverse motion of the translating tensioned beam is controlled by a time-varying external force besides a passive damping applied at the right boundary. Exponential stability of the closed loop system is proved. Simulation results demonstrate the effectiveness of the proposed controller.

• 전기의세계
• /
• v.22 no.1
• /
• pp.28-34
• /
• 1973

This paper treats with the sub-optimal control problem of noninear systems by approximation method. This method involves the approximation by linearization which provides the sub-optimal solution of non-linear control problems. The result of this work shows that, in the problem in which the controlled plant is characterized by an ordinary differential equation of first order, the solution obtained by this method coincides with the exact solution of problem. In of case of the second or higher order systems, it is proved analytically that this method of linearization produces the sub-optimal solution of the given problem. It is also shown that the sub-optimality of solution by the method can be evaluated by introducing the upper and lower bounded performance indices. Discussion is made on the procedure with some illustrative examples whose performance indices are given in the quadratic forms.

• Structural Engineering and Mechanics
• /
• v.28 no.4
• /
• pp.373-386
• /
• 2008

In this study stochastic analysis of non-linear dynamical systems under ${\alpha}$-stable, multiplicative white noise has been conducted. The analysis has dealt with a special class of ${\alpha}$-stable stochastic processes namely sub-Gaussian white noises. In this setting the governing equation either of the probability density function or of the characteristic function of the dynamical response may be obtained considering the dynamical system forced by a Gaussian white noise with an uncertain factor with ${\alpha}/2$- stable distribution. This consideration yields the probability density function or the characteristic function of the response by means of a simple integral involving the probability density function of the system under Gaussian white noise and the probability density function of the ${\alpha}/2$-stable random parameter. Some numerical applications have been reported assessing the reliability of the proposed formulation. Moreover a proper way to perform digital simulation of the sub-Gaussian ${\alpha}$-stable random process preventing dynamical systems from numerical overflows has been reported and discussed in detail.