A difference scheme for solving the Timoshenko beam equations with tip body

In this article, a Timoshenko beam with tip body and boundary damping is considered. A linearized three-level difference scheme of the Timoshenko beam equations on uniform meshes is derived by the method of reduction of order. The unique solvability, unconditional stability and convergence of the difference scheme are proved. The convergence order in maximum norm is of order two in both space and time. A numerical example is presented to demonstrate the theoretical results.     

Numerical Approximation and Error Analysis for the Timoshenko Beam Equations with Boundary Feedback

In this paper,the numerical approximation of a Timoshenko beam with bound- ary feedback is considered.We derived a linearized three-level difference scheme on uniform meshes by the method of reduction of order for a Timoshenko beam with boundary feedback.It is proved that the scheme is uniquely solvable,unconditionally stable and second order convergent in L_∞norm by using the discrete energy method. A numerical example is presented to verify the theoretical results.     

SEMI-DISCRETE AND FULLY DISCRETE PARTIALPROJECTION FINITE ELEMENT METHODS FOR THEVIBRATING TIMOSHENKO BEAM

1.IntroductionTheTimoshenkobeammodelisgivenbywherethebeamisconsidereddamped,drepresentsthebeamthicknessandI~[0,1].000istherotationofverticalfibersinthebeamandw(x)istheverticaldisplacementofthebeam'scenterline(underaverticalloadgivenbyg(x)).Analogoust...     

A difference scheme for a parabolic system modelling the thermoelastic contacts of two rods

In this article, we study a sequence of finite difference approximate solutions to a parabolic system, which models two dissimilar rods that each rod is fixed at one end and is free to expand or contact at the other end. A finite difference scheme is derived by the method of reduction of order on nonuniform mesh. The unique solvability, unconditional stability, and convergence of the difference scheme are proved. The convergence order is of order two in both time and space. The convergence of iterative algorithm for the difference scheme are also discussed. A numerical example is presented to demonstrate the theoretical results. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

ON LINEARIZED FINITE DIFFERENCE SIMULATION FOR THE MODEL OF NUCLEAR REACTOR DYNAMICS

1 hoeductIOuThe dynamics models of one--dimensional continuous medium nuclear reactor are the foelowing initial--boundary value problem of the formsubject to the innal conditionsand the boundary conditionsIn (1. 1), x denotes position along the reactor, which is regarded as a rod of length L, t denotes the time, u(t) the logarithm of the loud reactor POwer, v(x,t) the deviation of the temperature from equilibrium, a(x) the ratio of the temperature coefficient of reactivity to theynean life of…     

A finite difference scheme for solving the Timoshenko beam equations with boundary feedback

In this study, we derive a finite difference for a Timoshenko beam with boundary feedback by the method of reduction of order on uniform meshes. It is proved by the discrete energy method that the scheme is uniquely solvable, unconditionally stable and second order convergent in _L∞$L_{\infty }$ norm. Numerical results demonstrate the theoretical results.     

A fourth-order accurate method for fourth-order two-point nonlinear boundary value problems

A fourth-order accurate finite difference method is developed for a class of fourth order nonlinear two-point boundary value problems. The method leads to a pentadiagonal scheme in the linear cases, which often arise in the beam deflection theory. The convergence of the method is tested numerically on examples from the literature.     

A second‐order accurate difference scheme for the two‐dimensional Burgers' system

In this article, a Crank‐Nicolson‐type finite difference scheme for the two‐dimensional Burgers' system is presented. The existence of the difference solution is shown by Brouwer fixed‐point theorem. The uniqueness of the difference solution and the stability and L₂ convergence of the difference scheme are proved by energy method. An iterative algorithm for the difference scheme is given in detail. Furthermore, a linear predictor–corrector method is presented. The numerical results show that the predictor–corrector method is also convergent with the convergence order of two in both time and space. At last, some comments are provided for the backward Euler scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model

The numerical simulation of the dynamics of the molecular beam epitaxy (MBE) growth is considered in this article. The governing equation is a nonlinear evolutionary equation that is of linear fourth order derivative term and nonlinear second order derivative term in space. The main purpose of this work is to construct and analyze two linearized finite difference schemes for solving the MBE model. The linearized backward Euler difference scheme and the linearized Crank‐Nicolson difference scheme are derived. The unique solvability, unconditional stability and convergence are proved. The linearized Euler scheme is convergent with the convergence order of O(τ + h²) and linearized Crank‐Nicolson scheme is convergent with the convergence order of O(τ² + h²) in discrete L²‐norm, respectively. Numerical stability with respect to the initial conditions is also obtained for both schemes. Numerical experiments are carried out to demonstrate the theoretical analysis. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

带非线性源项的Riesz回火分数阶扩散方程的预估校正方法

本文针对带非线性源项的Riesz回火分数阶扩散方程,利用预估校正方法离散时间偏导数,并用修正的二阶Lubich回火差分算子逼近Riesz空间回火的分数阶偏导数,构造出一类新的数值格式.给出了数值格式在一定条件下的稳定性与收敛性分析,且该格式的时间与空间收敛阶均为二阶.数值试验表明数值方法是有效的.     

A linearized difference scheme for semilinear parabolic equations with nonlinear absorbing boundary conditions

A novel three level linearized difference scheme is proposed for the semilinear parabolic equation with nonlinear absorbing boundary conditions. The solution of this problem will blow up in finite time. Hence this difference scheme is coupled with an adaptive time step size, i.e., when the solution tends to infinity, the time step size will be smaller and smaller. Furthermore, the solvability, stability and convergence of the difference scheme are proved by the energy method. Numerical experiments are also given to demonstrate the theoretical second order convergence both in time and in space in L_∞-norm.     

A locking-free scheme for the LQR control of a Timoshenko beam

In this paper we analyze a locking-free numerical scheme for the LQR control of a Timoshenko beam. We consider a non-conforming finite element discretization of the system dynamics and a control law constant in the spatial dimension. To solve the LQR problem we seek a feedback control which depends on the solution of an algebraic Riccati equation. An optimal error estimate for the feedback operator is proved in the framework of the approximation theory for control of infinite dimensional systems. This estimate is valid with constants that do not depend on the thickness of the beam, which leads to the conclusion that the method is locking-free. In order to assess the performance of the method, numerical tests are reported and discussed.     

Compact difference scheme for two-dimensional fourth-order nonlinear hyperbolic equation

High-order compact finite difference method for solving the two-dimensional fourth-order nonlinear hyperbolic equation is considered in this article. In order to design an implicit compact finite difference scheme, the fourth-order equation is written as a system of two second-order equations by introducing the second-order spatial derivative as a new variable. The second-order spatial derivatives are approximated by the compact finite difference operators to obtain a fourth-order convergence. As well as, the second-order time derivative is approximated by the central difference method. Then, existence and uniqueness of numerical solution is given. The stability and convergence of the compact finite difference scheme are proved by the energy method. Numerical results are provided to verify the accuracy and efficiency of this scheme.     

二维线性双曲型方程Neumann边值问题的紧交替方向隐格式

对二维Neumann边界条件的线性双曲型方程建立了紧交替方向的隐格式.利用方程和边界条件得到在空间上的三阶与五阶导数的边界值,进而在内点、边界内点和边界角点分别建立9点、6点和4点紧差分格式;通过引进新的范数和L₂范数估计L_∞范数;借助能量估计、Gronwall不等式和Schwarz不等式等技巧,详细分析了差分格式在无穷范数下关于时间和空间分别为二阶和四阶收敛性,并给出了稳定性结果;通过数值算例,验证了理论分析结果.     

二维半线性反应扩散方程的交替方向隐格式

本文研究一类二维半线性反应扩散方程的差分方法.构造了一个二层线性化交替方向隐格式.利用离散能量估计方法证明了差分格式解的存在唯一性、差分格式在离散H~1模下的二阶收敛性和稳定性.最后给出两个数值例子验证了理论分析结果.     

Asymptotic analysis of Poisson's equation in a thin domain and its application to thin-walled elastic beams and tubes

We study the limit behaviour of solution of Poisson's equation in a class of thin two-dimensional domains, both simply connected or single-hollowed, as its thickness becomes very small. The method is based on a transformation of the original problem into another posed on a fixed domain, obtention of a priori estimates and convergence results when thickness parameter tends to zero. As an important application of abstract results we obtain the limit expressions for functions appearing in elastic beam theories as torsion and warping functions. In this way, we provide a mathematical justification and a correct definition of torsion, warping and Timoshenko functions and constants that should be used in the open and closed thin-walled elastic beam theories. © 1998 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Imposing boundary conditions in Sinc method using highest derivative approximation

Sinc approximate methods are often used to solve complex boundary value problems such as problems on unbounded domains or problems with endpoint singularities. A recent implementation of the Sinc method [Li, C. and Wu, X., Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives, Applied Mathematical Modeling 31 (1) 2007 1–9] in which Sinc basis functions are used to approximate the highest derivative in the governing equation of the boundary value problem is evaluated for structural mechanics applications in which interlaminar stresses are desired. We suggest an alternative approach for specifying the boundary conditions, and we compare the numerical results for analysis of a laminated composite Timoshenko beam, implementing both Li and Wu’s approach and our alternative approach for applying the boundary conditions. For the Timoshenko beam problem, we obtain accurate results using both approaches, including transverse shear stress by integration of the 3D equilibrium equations of elasticity. The beam results indicate our approach is less dependent on the selection of the Sinc mesh size than Li and Wu’s SIHD. We also apply SIHD to analyze a classical laminated composite plate. For the plate example, we experience difficulty in obtaining a complete system of equations using Li and Wu’s approach. For our approach, we suggest that additional necessary information may be obtained by applying the derivatives of the boundary conditions on each edge. Using this technique, we obtain accurate results for deflection and stresses, including interlaminar stresses by integration of the 3D equilibrium equations of elasticity. Our results for both the beam and the plate problems indicate that this approach is easily implemented, has a high level of accuracy, and good convergence properties.     

Finite element methods for Timoshenko beam,circular arch and Reissner-Mindlin plate problems

In this paper some finite element methods for Timoshenko beam, circular arch and Reissner-Mindlin plate problems are discussed. To avoid locking phenomenon, the reduced integration technique is used and a bubble function space is added to increase the solution accuracy. The method for Timoshenko beam is aligned with the Petrov-Galerkin formulation derived in Loula et al. (1987) and can be naturally extended to solve the circular arch and the Reissner-Mindlin plate problems. Optimal order error estimates are proved, uniform with respect to the small parameters. Numerical examples for the circular arch problem shows that the proposed method compares favorably with the conventional reduced integration method.     

A high order difference method for fractional sub-diffusion equations with the spatially variable coefficients under periodic boundary conditions

In this paper, we propose a difference scheme with global convergence order $O(\tau^{2}+h^4)$ for a class of the Caputo fractional equation. The difficulty caused by the spatially variable coefficients is successfully handled. The unique solvability, stability and convergence of the finite difference scheme are proved by use of the Fourier method. The obtained theoretical results are supported by numerical experiments.     

一个 油藏模型的二阶差分格式

本文应用伸缩变换和阶降法在动网格上对石油科学研究中提出的一个漂移聚集油藏模型建立了一个三层线性化弱耦合差分格式，证明了所构造的差分格式是唯一可解的且在ｌ２范数下以Ｏ（ｒ＾２＋ｈ＾２）阶收敛。所建立的差分方程组均是三对角的，可用追赶法求解。