功能梯度矩形板的非线性自由振动

研究了功能梯度矩形薄板的非线性自由振动问题.采用幂律分布规律描述功能梯度材料沿厚度的梯度性质,基于von Kámán理论,建立了功能梯度薄板的非线性振动控制方程.应用Bubnov-Galerkin法得到了功能梯度矩形薄板的单模态非线性振动的时域常微分方程,借助其势能函数分析了系统的周期振动状态.采用Lindstedt-Poincaré法和Runge-Kutta法分别获得了功能梯度矩形薄板单模态非线性周期振动的摄动解和数值解.研究表明:功能梯度薄板非线性振动控制方程中包含表征拉弯耦合效应的控制项,这导致其常微分方程中出现二次项;系统振幅在板横向的正负两个方向上是不相等的,其振动存在关于板中面的不对称性.     

轴对称自由振动功能梯度材料微圆板中的热弹性阻尼

基于经典弹性薄板理论和单向耦合热传导理论,研究了材料性质沿厚度连续变化的功能梯度微圆板的热弹性阻尼特性．首先,考虑热力耦合效应,建立了功能梯度微圆板轴对称横向自由振动微分方程．然后,忽略温度梯度在面内的变化,建立了单向耦合变系数一维热传导方程．采用分层均匀化近似方法,将变系数热传导方程转化为一系列常系数的微分方程,利用上下表面的热边界条件和层间连续性条件获得了微圆板温度场解析解．将所得温度场代入微圆板的自由振动微分方程,得到了包含热弹性阻尼的复频率,从而获得了反映热弹性阻尼水平的逆品质因子．最后,针对材料性质沿板厚按幂函数变化的陶瓷-金属功能梯度微圆板,定量地分析材料梯度指数、几何尺寸、边界条件、温度环境等对微圆板热弹性阻尼的影响．     

正交异性非线性变厚度环形板的分析

讨论了正交异性非线性变厚度混凝土环形板问题.建立了荷载为q(r)及q(r)*sinnθ条件下的环形板平衡条件式,导出了挠度表达的Euler方程.据此获得一组常系数微分方程和它们的特征值,可以求解正交异性非线性变厚度环形板的挠度及内力.     

功能梯度材料与结构的若干力学问题研究进展

功能梯度材料的宏观材料特性在空间上是连续变化的,因此即使在线弹性理论范围内,由于控制偏微分方程是变系数的,相应的力学分析具有很大的挑战性.综述了功能梯度材料与结构若干力学问题的最新研究进展,包括功能梯度材料梁、板、壳结构的解析解与半解析解以及简化理论的研究、功能梯度材料结构的数值计算方法研究、功能梯度材料的断裂力学研究.最后对未来功能梯度材料与结构的力学研究进行了展望.     

功能梯度矩形板问题三维理论的半解析解

针对功能梯度材料矩形板问题,基于三维弹性理论,将位移和应力分量作为基本变量,通过双三角级数将其控制微分方程转化为常微分方程组的边值问题。采用插值矩阵法直接对常微分方程组边值问题进行求解,得到了功能梯度材料矩形板三维位移、应力场的半解析解。通过算例给出了材料参数按指数形式和幂函数形式变化情况下的功能梯度板的弯曲问题。对比有限元法和状态空间法,结果表明：本文提出的状态空间与插值矩阵法相结合的半解析法能有效地分析材料参数按任意形式连续变化的功能梯度矩形板问题,且具有良好的精度,精度可达10-4量级,能够满足工程需要;与其他方法相比,本文方法具有实施便捷、计算量小等优点,根据其力学场分析结果可设计出满足各种不同需求的功能梯度材料。     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板,考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况,应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件,利用伽辽金法得到了横向非线性固有振动方程,并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程,得出非线性固有频率表达式．通过算例,分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明,非线性固有频率随金属含量的增加而降低;随转速和圆板厚度的增大而升高;随功能梯度圆板表面温度的升高而降低．     

常子午线曲率薄壳轴对称几何非线性问题的复变量方程

1.几何非线性问题的基本方程在本世纪初,Reissner H.和Meissner E.利用在线性薄壳理论中存在的静力-几何比拟关系,将线弹性薄壳轴对称问题,归结为以应力函数和转角为未知量的两个常微分方程。以后,人们利用这两个方程的相似性,引入复未知函数,把一些典型壳体的方程简化为一个二阶变系数常微分方程,为这些问题的求解带来极大的便利。本文将这一方法推广到薄壳大位移问题,导出用复未知函数表示的常子午线曲率壳体轴对称变形的非线性微分方程。从这个一般方程可以直接得到关于柱壳,锥壳,圆球壳,环壳和圆板几何非线性问     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板,考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况,应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件,利用伽辽金法得到了横向非线性固有振动方程,并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程,得出非线性固有频率表达式．通过算例,分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明,非线性固有频率随金属含量的增加而降低;随转速和圆板厚度的增大而升高;随功能梯度圆板表面温度的升高而降低．     

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

功能梯度材料明德林矩形微板的热弹性阻尼

功能梯度材料微板谐振器热弹性阻尼的建模和预测是此类新型谐振器热?弹耦合振动响应的新课题. 本文采用数学分析方法研究了四边简支功能梯度材料中厚度矩形微板的热弹性阻尼. 基于明德林中厚板理论和单向耦合热传导理论建立了材料性质沿着厚度连续变化的功能梯度微板热弹性自由振动控制微分方程. 在上下表面绝热边界条件下采用分层均匀化方法求解变系数热传导方程, 获得了用变形几何量表示的变温场的解析解. 从而将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的偏微分方程. 然后,利用特征值问题在数学上的相似性,求得了四边简支条件下功能梯度材料明德林矩形微板的复频率解析解, 进而利用复频率法获得了反映谐振器热弹性阻尼水平的逆品质因子. 最后, 给出了材料性质沿板厚按幂函数变化的陶瓷?金属组分功能梯度矩形微板的热弹性阻尼数值结果. 定量地分析了横向剪切变形、材料梯度变化以及几何参数对热弹性阻尼的影响规律. 结果表明, 采用明德林板理论预测的热弹性阻尼值小于基尔霍夫板理论的预测结果, 而且两者的差别随着相对厚度的增大而变得显著.      

Asymptotic Solutions of a First-Order Differential Equation with Deviating Argument and Slowly Varying Coefficients

The object of this paper is to study the problem of constructing an approximate solution of a first-order weakly nonlinear ordinary differential equation with deviating argument and slowly varying coefficients. On the basis of asymptotic techniques in nonlinear mechanics, we construct an algorithm for the asymptotic integration of the differential equation under consideration.__________Published in Neliniini Kolyvannya, Vol. 7, No. 4, pp. 475–486, October–December, 2004.     

Van der Pol type self-excited micro-cantilever probe of atomic force microscopy

A technique for dimensional reduction of nonlinear delay differential equations (DDEs) with time-periodic coefficients is presented. The DDEs considered here have a canonical form with at most cubic nonlinearities and periodic coefficients. The nonlinear terms are multiplied by a perturbation parameter. Perturbation expansion converts the nonlinear response problem into solutions of a series of nonhomogeneous linear ordinary differential equations (ODEs) with time-periodic coefficients. One set of linear nonhomogeneous ODEs is solved for each power of the perturbation parameter. Each ODE is solved by a Chebyshev spectral collocation method. Thus we compute a finite approximation to the nonlinear infinite-dimensional map for the DDE. The linear part of the map is the monodromy operator whose eigenvalues characterize stability. Dimensional reduction on the map is then carried out. In the case of critical eigenvalues, this corresponds to center manifold reduction, while for the noncritical case resonance conditions are derived. The accuracy of the nonlinear Chebyshev collocation map is demonstrated by finding the solution of a nonlinear delayed Mathieu equation and then a milling model via the method of steps. Center manifold reduction is illustrated via a single inverted pendulum including both a periodic retarded follower force and a nonlinear restoring force. In this example, the amplitude of the limit cycle associated with a flip bifurcation is found analytically and compared to that obtained from direct numerical simulation. The method of this paper is shown by example to be applicable to systems with strong parametric excitations.     

Analysis of velocity equation of steady flow of a viscous incompressible fluid in channel with porous walls

Steady flow of a viscous incompressible fluid in a channel, driven by suction or injection of the fluid through the channel walls, is investigated. The velocity equation of this problem is reduced to nonlinear ordinary differential equation with two boundary conditions by appropriate transformation and convert the two‐point boundary‐value problem for the similarity function into an initial‐value problem in which the position of the upper channel. Then obtained differential equation is solved analytically using differential transformation method and compare with He's variational iteration method and numerical solution. These methods can be easily extended to other linear and nonlinear equations and so can be found widely applicable in engineering and sciences. Copyright © 2009 John Wiley & Sons, Ltd.     

非半单分叉的Normal Form计算

在文［１］基础上，提出一种仅知道派生线性系统零实部特征值时求解非线性系统非半单分叉ＮｏｒｍａｌＦｏｒｍ的方法．通过适当的分类，将要求解的线性代数方程组分为若干相互独立的方程组．将所求系数向量按字典序列排列后，各独立方程组的系数矩阵是上三角矩阵．在非共振情形，各系数向量可按反字典序列递推求出．在共振情形，根据文中的二个定理，巧妙地由一简单的常数矩阵的最大秩子矩阵，定位其系数矩阵的满秩子矩阵，解决了这类方程组的降维简化．通过消元法，把简化后的方程化成类似于半单分叉ＮｏｒｍａｌＦｏｒｍ求解过程中方程的形式，其解法也类似．该方法非常易于在计算机代数软件平台上程序化．     

On the existence of non-polarized azimuthal-symmetric electromagnetic waves in circular dielectric waveguide filled with nonlinear isotropic homogeneous medium

The propagation of monochromatic nonlinear symmetric hybrid waves in a cylindrical nonlinear dielectric waveguide is considered. The physical problem is reduced to solving a transmission eigenvalue problem for a system of ordinary differential equations. Spectral parameters of the problem are propagation constants of the waveguide. The problem is reduced to the new type of nonlinear eigenvalue problem. The analytical method of solving this problem is presented. New propagation regime is found.     

Equations for nonlinear waves on shallow water

Steady-state free plane surface waves are investigated in the nonlinear formulation by means of the narrow streak method. The second-order differential equation derived can readily be integrated numerically. The first approximation of the equation is integrated. The next approximation leads to a nonlinear first-order differential equation which is traditional to the problem considered and whose coefficients have a simple physical meaning.     

1/3 subharmonic solution of elliptical sandwich plates

IntroductionInrecentyears,withtheessentialadvantageoflightweightandhighrigidity ,sandwichplatesandshellshavebeenusedasanimportantpatternofstructuralelementsinaeronautical,astronauticalandnavalengineering .However,nonlinearproblemsforsandwichplatesandshellsareonlyinvestigatedbyafewbecauseofthedifficultiesofnonlinearmathematicalproblems.LiuRen_huaiandXuJia_chu[1,2 ]andothershavemadesomeinvestigationsinthisfield .Bifurcationofnonlinearvibrationforsandwichplateshasnotyetbeeninvestigated .Inthisp…     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

Buckling load optimization of beams

In this paper, shape optimization is used to optimize the buckling load of a Euler–Bernoulli beam having constant volume. This is achieved by varying appropriately the beam cross section so that the beam buckles with the maximum or a prescribed buckling load. The problem is reduced to a nonlinear optimization problem under equality and inequality constraints as well as specified lower and upper bounds. The evaluation of the objective function requires the solution of the buckling problem of a beam with variable stiffness subjected to an axial force. This problem is solved using the analog equation method for the fourth-order ordinary differential equation with variable coefficients. Besides its accuracy, this method overcomes the shortcomings of a possible FEM solution, which would require resizing of the elements and recomputation of their stiffness properties during the optimization process. Several example problems are presented that illustrate the method and demonstrate its efficiency.     

Geometric Optics and Instability for Semi-Classical Schrödinger Equations

We prove some instability phenomena for semi-classical (linear or) nonlinear Schrödinger equations. For some perturbations of the data, we show that for very small times, we can neglect the Laplacian, and the mechanism is the same as for the corresponding ordinary differential equation. Our approach allows smaller perturbations of the data, where the instability occurs for times such that the problem cannot be reduced to the study of an ordinary differential equation.