• Title/Summary/Keyword: linear differential equations

Search Result 539, Processing Time 0.022 seconds

Intelligent computer modelling and simulation for the large amplitude of nano systems

  • Yi, Wenjuan
    • Advances in nano research
    • /
    • v.13 no.1
    • /
    • pp.63-75
    • /
    • 2022
  • The nonlinear dynamic behavior of a nonuniform small-scale nonlocal beam is investigated in this work. The nanobeam is theoretically modeled using the nonlocal Eringen theory, as well as a few of Von-nonlinear Kármán's theories and the classical beam theory. The Hamilton principle extracts partial differential equations (PDE) of an axially functionally graded (AFG) nano-scale beam consisting of SUS304 and Si3N4 throughout its length, and an elastic Winkler-Pasternak substrate supports the tapered AFG nanobeam. The beam thickness is a function of beam length, and it constantly varies throughout the length of the beam. The numerical solution strategy employs an iteration methodology connected with the generalized differential quadratic method (GDQM) to calculate the nonlinear outcomes. The nonlinear numerical results are presented in detail to examine the impact of various parameters such as nonlinear amplitude, nonlocal parameter, the component of the elastic foundation, rate of cross-section change, and volume fraction parameter on the linear and nonlinear free vibration characteristics of AFG nanobeam.

A REVIEW ON THE MATHEMATICAL ASPECTS OF FLUID FLOW PROBLEMS IN AN INFINITE CHANNEL WITH ARBITRARY BOTTOM TOPOGRAPHY

  • Chakrabarti, A.;Martha, S.C.
    • Journal of applied mathematics & informatics
    • /
    • v.29 no.5_6
    • /
    • pp.1583-1602
    • /
    • 2011
  • A special system of partial differential equations (PDEs) occur in a natural way while studying a class of irrotational inviscid fluid flow problems involving infinite channels. Certain aspects of solutions of such PDEs are analyzed in the context of flow problems involving multiple layers of fluids of different constant densities in a channel associated with arbitrary bottom topography. The whole analysis is divided into two parts-part A and part B. In part A the linearized theory is employed along with the standard Fourier analysis to understand such flow problems and physical quantities of interest are derived analytically. In part B, the same set of problems handled in part A are examined in the light of a weakly non-linear theory involving perturbation in terms of a small parameter and it is shown that the original problems can be cast into KdV type of nonlinear PDEs involving the bottom topography occurring in one of the coefficients of these equations. Special cases of bottom topography are worked out in detail and expressions for quantities of physical importance are derived.

Dynamic load concentration caused by a break in a Lamina with viscoelastic matrix

  • Reza, Arash;Sedighi, Hamid M.;Soleimani, Mahdi
    • Steel and Composite Structures
    • /
    • v.18 no.6
    • /
    • pp.1465-1478
    • /
    • 2015
  • The effect of cutting off fibers on transient load in a polymeric matrix composite lamina was studied in this paper. The behavior of fibers was considered to be linear elastic and the matrix behavior was considered to be linear viscoelastic. To model the viscoelastic behavior of matrix, a three parameter solid model was employed. To conduct this research, finite difference method was used. The governing equations were obtained using Shear-lag theory and were solved using boundary and initial conditions before and after the development of break. Using finite difference method, the governing integro-differential equations were developed and normal stress in the fibers is obtained. Particular attention is paid the dynamic overshoot resulting when the fibers are suddenly broken. Results show that considering viscoelastic properties of matrix causes a decrease in dynamic load concentration factor and an increase in static load concentration factor. Also with increases the number of broken fibers, trend of increasing load concentration factor decreases gradually. Furthermore, the overshoot of load in fibers adjacent to the break in a polymeric matrix with high transient time is lower than a matrix with lower transient time, but the load concentration factor in the matrix with high transient time is lower.

Free Vibrations of Clamped Circular Arches with Linear Variable Cross-Section (1차원 변화단면을 갖는 양단고정 원호아치의 자유진동 해석)

  • Lee, Byoung Koo
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.9 no.4
    • /
    • pp.1-8
    • /
    • 1989
  • The main purpose of the present paper is to present both the fundamental frequency and some higher free vibration frequencies for circular arches with variable section, in which rotatory inertia is included. The differential equations are derived for the in-plan free vibration of elastic circular arches with variable section. These equations were solved numerically for the linear variable circular cross-section with clamped-clamped end constraint. As the numerical results, the four lowest nondimensional natural frequencies presented as functions of the nondimensional system parameters : the end moment of inertia to crown moment of inertia ratio, the slenderness ratio, and the opening angle. The effect of rotatory inertia on the nondimensional natural frequency is also reported.

  • PDF

Dynamic Behavior of Rotating Shaft System Corresponding to Operating Modes (운전모드에 따른 회전축계의 동적거동)

  • Kim, Sang-Hwan
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.20 no.9
    • /
    • pp.2744-2751
    • /
    • 1996
  • In case of limited power supply, a rotating shaft system may not reach its operating speed that is greater than its critical speed, but the speed oscillates with small ampllitude near critical speed. As a result, it is considered that the operating mode plays an important role in the smooth start of machines. In order to investigate the dynamic behaviors of the rotating shaft system at the beginning stage, one has derived the equations of motion whose degrees of freedom is three, two translations and one rotation. The simultaneous differential equations are numerically solved by using runge-Kutta method, and thus the small time step length could be required corresponding to the stability of solution. Three types of operating modes dependent upon the driving torque rate have been numerically investigated according to the maximum displacement of shaft center. The first type of relation is linear, the second type is composed of two linear curves recommended by machine manufacturer, and the last one is the proposed torque curve reflecting the frequency response curve of one degree of freedom system. For the second type of modes, it is found that the optimal range of intermediate speed to the critical speed lies between 0.8 and 0.9. In addition to that, the maximum displacement can be reduced more if the third type of mode is utilized.

A Study on the Analysis and Design of Nonlinear Control Systems using Personal Computer (개인용 컴퓨터를 이용한 비선형 제어 시스템의 해석 및 설계에 관한 연구)

  • Nam, Moon-Hyun;Jeong, Cheol
    • Proceedings of the KIEE Conference
    • /
    • 1987.11a
    • /
    • pp.82-85
    • /
    • 1987
  • The objective of this paper is to develop computer programs to aid in the design and analysis of control systems in which nonlinear characteristics exist. Control systems are dynamic systems, which can be described using various mathematical models. A convenient model for digital computer simulation is the state model in which described using a set of linear and non linear first order differential equations. The digital simulation was performed on a IBM PC/XT personal computer, and the computer language was BASIC. There are four possible configurations from which a user may choose. When running a program, the user is asked to enter the system parameters according to a specified control system configurations are; 1. A control system with a nonlinear element followed by a plant in a feedback configurations(NLSVF1). 2. A control system with a nonlinear device situated between two plants in a feedback configurations(NLSVF2). 3. A control system with a nonlinear element followed by a plant, followed by a the dealy in feedback configurations(TLAG). 4. A motor and load with a backlash nonlinearity between dynamic portions of the motor/load configurations (BACKLASH). The matrix from state equations are integrated using combination the trapezoidal method and fixed point iteration. Several cases which have nonlinearity were implemented on the computer and the results were discussed.

  • PDF

A novel method for solving structural problems: Elastoplastic analysis of a pressurized thick heterogeneous sphere

  • Abbas Heydari
    • Advances in Computational Design
    • /
    • v.9 no.1
    • /
    • pp.39-52
    • /
    • 2024
  • If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

INFLUENCE OF THERMAL CONDUCTIVITY AND VARIABLE VISCOSITY ON THE FLOW OF A MICROPOLAR FLUID PAST A CONTINUOUSLY MOVING PLATE WITH SUCTION OR INJECTION

  • Salem, A.M.;Odda, S.N.
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.9 no.2
    • /
    • pp.45-53
    • /
    • 2005
  • This paper investigates the influence of thermal conductivity and variable viscosity on the problem of micropolar fluid in the presence of suction or injection. The fluid viscosity is assumed to vary as an exponential function of temperature and the thermal conductivity is assumed to vary as a linear function of temperature. The governing fundamental equations are approximated by a system of nonlinear ordinary differential equations and are solved numerically by using shooting method. Numerical results are presented for the distribution of velocity, microrotation and temperature profiles within the boundary layer. Results for the details of the velocity, angular velocity and temperature fields as well as the friction coefficient, couple stress and heat transfer rate have been presented.

  • PDF

A Sudy on the Undamped Forced Vibration of Nonlinear Two-Degree-of-Freedom Systems (비선형 2자유도계의 비감쇠 강제진동 연구)

  • 박철희;박선재;윤영석
    • Transactions of the Korean Society of Mechanical Engineers
    • /
    • v.12 no.2
    • /
    • pp.193-199
    • /
    • 1988
  • The forced vibrations of nondissipative nonlinear two-degree-of-freedom system, subjected to periodic forcing functions, are investigated by use of the method of slowly changing phase and amplitude. The first order differential equations are derived for nonrationally solutions and the coupled nonlinear algebraic equations for stationary solutions. Through investigating the response curves of the system, which are obtained numerically by using Newton-Raphson method, it is found that the resonances can occur at more than the number of degree-of-freedom of the system depending on the relation between the nonlinear spring parameters, which has no counterpart in linear systems.

CALCULATION METHODS OF SOLAR ATMOSPHERIC MODEL (태양대기모델 계산법)

  • KIM KAP-SUNG
    • Publications of The Korean Astronomical Society
    • /
    • v.15 no.spc2
    • /
    • pp.65-71
    • /
    • 2000
  • We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

  • PDF