第二类拉格朗日方程在无初速释放动力学问题中的应用

第二类拉格朗日方程在求解复杂动力学问题中有着广泛的应用,其特点是：只要求出系统在一般位置时的动能及相应于各广义坐标的广义力,经过程序化的求导运算就可获得控制系统的动力学方程。正是因为需要对系统的动能进行求导运算,第二类拉格朗日方程不能在系统的特殊位置求写系统的动能,而求写系统在一般位置时的动能很多情况下是一件不容易的事情,这正是应用拉格朗日方程的最大障碍。我们经过理论分析发现对无初速释放动力学问题,第二类拉格朗日方程提供了简便的求解途径。     

考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

冷压状态方程计算的新方法和材料相图的研究

本文运用位力(virial)定理严格地给出了材料的冷压状态方程.本文的理论没有引入对系统势能的任何假定,只需计算系统的总动能,因而在状态方程的研究上得到了新的进展.本文采用了经典TFD统计模型,给出了计算材料状态方程的新方法.此方法可以运用于整个凝聚体系统,从零压开始计算材料各个压缩度下的压力,考虑到经典TFD模型的统计与真实量子力学统计之间的差异.此方法对每种元素的总动能计算作了修正.用此状态方程的新方法计算,可以得到在整个压缩度范围(包括零压附近)压力与实测值一致的冷压状态方程, 对Li,Na,Al,Fe,Ag和U等元素以及Fe-Ni,Al-Cu等二元合金进行了数值计算,结果与实测值符合得较好. 文中还探讨了材料相图的理论研究.冷压状态方程计算的新方法和相图的理论研究将为材料的设计提供依据.      

NONLINEAR INTERNAL RESONANCE OF FUNCTIONALLY GRADED CYLINDRICAL SHELLS USING THE HAMILTONIAN DYNAMICS

Internal resonance in nonlinear vibration of functionally graded （FG） circular cylin- drical shells in thermal environment is studied using the Hamiltonian dynamics formulation. The material properties are considered to be temperature-dependent. Based on the Karman-Donnell＇s nonlinear shell theory, the kinetic and potential energy of FG cylindrical thin shells are formu- lated. The primary target is to investigate the two-mode internal resonance, which is triggered by geometric and material parameters of shells. Following a secular perturbation procedure, the underlying dynamic characteristics of the two-mode interactions in both exact and near resonance cases are fully discussed. It is revealed that the system will undergo a bifurcation in near resonance case, which induces the dynamic response at high energy level being distinct from the motion at low energy level. The effects of temperature and volume fractions of composition on the exact resonance condition and bifurcation characteristics of FG cylindrical shells are also investigated.     

自由场空泡溃灭过程能量转化机制研究

综合应用实验与数值模拟方法, 深入讨论了自由场空泡溃灭过程中的能量转化机制. 在实验研究中, 应用纹影法记录了空泡溃灭的演变过程, 提取了空泡在溃灭过程中的半径, 溃灭速度等数据, 结合空泡势能和动能方程, 描述了空泡能量的转化过程. 在开展数值模拟分析时, 运用弱可压缩流体质量守恒方程和动量方程, 建立了三维数值模型用以模拟空泡在自由场中的溃灭过程, 并且由结果中获取了空泡溃灭过程中的压力及速度变化规律, 揭示了空泡在溃灭过程中能量转化机制. 研究结果表明: (1) 自由场空泡在溃灭过程中, 空泡势能与空泡半径具有相同的演化趋势, 空泡动能与势能变化趋势相反; 当空泡达到最大半径处时, 空泡势能最大, 流场动能为零. (2) 溃灭后期在空泡周围会形成高压区域, 该区域的压力梯度与速度梯度较高, 随着空泡收缩, 高压区域面积逐渐减小. (3) 空泡在自由场中发生溃灭时, 空泡势能不断转化为流场动能, 在溃灭时刻可以明显观察到冲击波现象, 空泡的大部分能量会在此时转化为冲击波的波能.      

超晶格的位错模型与系统的动力学稳定性

在经典力学框架内和Seeger方程基础上,讨论了超晶格界面附近的位错动力学行为,指出了由于系统的分叉或混沌将导致位错的运动与堆积,造成了超晶格的分层或断裂;同时,也指出了,将生长过程中的超晶格置于适当的声场中将应力减至最小,或者适当调节系统参数就可最大限度的保证系统的动力学稳定性.首先,引入阻尼项,把描述一般位错运动的Seeger方程化为了超晶格系统的广义摆方程.利用Jacobian椭圆函数和椭圆积分分析了无扰动系统的相平面特征,并解析地给出了系统的解和粒子振动周期.其次,利用Melnikov方法分析了系统相平面上三类轨道的分叉性质和进入Smale马蹄意义下的混沌行为,找到了系统的全局分叉与系统进入混沌的临界条件.结果表明,系统的临界条件与它的物理参数有关,只需适当调节这些参数就可以原则上避免、控制分叉或混沌的出现,进一步保证生长过程的稳定性和超晶格材料的完整性.     

软体连续机械臂动力学建模与仿真

由于软体机械臂的质量是沿臂的长度连续分布,因此采用拉格朗日方法建立软体机械臂的动力学模型时,涉及计算复杂的积分运算,采用离散化的集中质量模型降低了计算的复杂性,但准确性不足.为了提高软体机械臂动力学建模与仿真的准确度和计算效率,本文采用模态方法对软体机械臂进行运动学描述,并从能量的角度分析软体机械臂动力学特性,研究发现...     

A new dynamic subgrid-scale model using artificial neural network for compressible flow

The subgrid-scale (SGS) kinetic energy has been used to predict the SGS stress in compressible flow and it was resolved through the SGS kinetic energy transport equation in past studies. In this paper, a new SGS eddy-viscosity model is proposed using artificial neural network to obtain the SGS kinetic energy precisely, instead of using the SGS kinetic energy equation. Using the infinite series expansion and reserving the first term of the expanded term, we obtain an approximated SGS kinetic energy, which has a high correlation with the real SGS kinetic energy. Then, the coefficient of the modelled SGS kinetic energy is resolved by the artificial neural network and the modelled SGS kinetic energy is more accurate through this method compared to the SGS kinetic energy obtained from the SGS kinetic energy equation. The coefficients of the SGS stress and SGS heat flux terms are determined by the dynamic procedure. The new model is tested in the compressible turbulent channel flow. From the a posterior tests, we know that the new model can precisely predict the mean velocity, the Reynolds stress, the mean temperature and turbulence intensities, etc.     

多层复合壳体三维振动分析的谱——微分求积混合法

对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求. 建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值. 本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 数值验证结果表明该方法收敛性好,计算精度高.      

Kinetic and potential energy in the first law of thermodynamics

Kinetic and potential energy are included in the first law of thermodynamics in quite a contradictory way. Whereas in thermodynamics the total energy is understood as the sum of internal, kinetic and potential energy, the total energy in continuum mechanics incorporates only internal and kinetic energy, the potential energy being part of the work. The Gibbs ' fundamental equation is also occasionally extendend to contain a term for the potential energy. Some serious contradictions may result from this. As is first shown, kinetic and potential energy do not have any influence on the internal energy as long as relativistic effects are excluded. The Gibbs' fundamental equation therefore describes exchange processes between the “internal variables? of a system and its surroundings. Proceeding from this result one obtains a general definition of heat in open systems, including electromagnetic reactions, surfaceeffects and variable mole numbers. Exchange processes between the “external variables? of a system and its surroundings and hence also the influence of kinetic and potential energy are described by another independent equation, i.e. the energy equation of mechanics. Addition of both equations leads to the heat definition which is usually but under some further neglects given in textbooks. This definition has considerable disadvantages compared to the one derived before. In particular it is no longer possible to realize how mechanical and thermal energy are transformed into each other, which may give rise to errors.     

Kinetic solutions of the Boltzmann-Peierls equation and its moment systems

The evolution of heat in crystalline solids is described at low temperatures by the Boltzmann-Peierls equation, which is a kinetic equation for the phase density of phonons.In this study, we solve initial value problems for the Boltzmann-Peierls equation in relation to the following issues: In thermodynamics, a given kinetic equation is usually replaced by a truncated moment system, which in turn is supplemented by a closure principle so that a system of PDEs results for some moments as thermodynamic variables. A very popular closure principle is the maximum entropy principle, which yields a symmetric hyperbolic system. In recent times, this strategy has led to serious studies on two problems that might arise: 1. Do solutions of the maximum entropy principle exist? 2. Is the physics that is embodied by the kinetic equation more or less equivalently displayed by the truncated moment system? It was Junk who proved for the BOLTZMANN equation of gases that maximum entropy solutions do not exist. The same failure appears for the Fokker-Planck equation, which was proved by means of explicit solutions by Dreyer, Junk, and Kunik.This study has two main objectives:1. We give a positive existence result for the maximum entropy principle if the underlying kinetic equation is the Boltzmann-Peierls equation. In other words we show that the maximum entropy principle can be used here to establish a closed hyperbolic moment system of PDEs. However, the intent of the paper is by no means a general justification of the maximum entropy principle.2. We develop an approximative method that allows the solutions of the kinetic equations to be compared with the solutions of the hyperbolic moment systems. To this end we introduce kinetic schemes that consists of free flight periods and two classes of update rules. The first class of rules is the same for the kinetic equation as well as for the maximum entropy system, while the second class of update rules contains additional rules for the maximum entropy system. It is shown that if a sufficient number of moments are taken into account, the two solutions converge to each other. However, in terms of numerical effort, the presented solver for the kinetic equation clearly outperforms the one for the maximum entropy principle.Received: 15 August 2003, Accepted: 8 November 2003, Published online: 11 February 2004PACS: 02.30.Jr, 02.60.Cb, 05.30.Jp, 44.10. + i, 63.20.-e, 66.70. + f, 65.40.Gr Correspondence to: M. Herrmann     

新型自然循环自动清洗式蒸发器及设计原理

以高效蒸发器为目标，旨在以自然循环替代高能耗的强制循环泵、以自动清洗取代加热面盐垢的周期性停车清洗．主要研究沸腾室内的不平衡汽化过程机理，提出动力温度的新概念，采用显热与汽化热转换平衡的研究方法，得到自然循环推动力的工程设计计算式，并且由此引出推动力强化的基本原则，进而获得中试的成功．还提出加热室必须采用能够低流速下自转清洗又高效传热强化的螺旋齿带，才能保障循环总阻力低于推动力．     

Two Theories of Complex Bodies, one Lagrangean, the other Kinetic, Have a Common Ground

I have found in previous works that most special models proposed to represent bodies with some type of microstructure can be classified easily under the general umbrella of a theory where each element of the continuum is thought of as a Lagrangian system. To study phenomena in ‘kinetic’ continua I proposed an apparently different approach; the outcome is again a set of evolution equations. They mimic equations familiar in continua with affine microstructure: a Cauchy’s equation and an equation of balance of tensor moment of momentum, with the addition, however, of an equation of balance for a ‘Reynolds’ tensor’, an equation which, in a sense, shifts the boundary between kinetic and thermal properties of matter. I will show that there is no contrast between the two approaches. The latter one is based on an adequate and appropriately justified expression of the kinetic energy of the continuum, comprising the trace of the quoted Reynolds’ tensor and thus importing into the mechanical energy a term usually accounted by additional heat.     

基于比例边界有限元法的高阶透射边界应用

基于比例边界有限元法和连分式展开推导了无限域弹性动力分析的求解方程,实现了一种局部的高阶透射边界. 采用改进的连分式法求解无限域的动力刚度矩阵,克服了原连分式算法可能会造成矩阵运算病态的问题. 该局部高阶透射边界在时域里表示为一阶常微分方程组,其稳定性取决于其系数矩阵的广义特征值问题. 如果出现虚假模态,采用移谱法来校正系数矩阵以消除虚假模态. 通过两个算例验证了该高阶透射边界的精确性、鲁棒性.     

近似不可压软材料动力分析的完全拉格朗日物质点法

物质点法(MPM)在模拟非线性动力问题时具有很好的效果,其已被广泛应用于许多大变形动力问题的分析中.然而传统的MPM在模拟不可压或近似不可压材料的动力学行为时会产生体积自锁,极大地影响模拟精度和收敛性.本文针对近似不可压软材料的大变形动力学行为,提出一种混合格式的显式完全拉格朗日物质点法(TLMPM).首先基于近似不可压软材料的体积部分应变能密度,引入关于静水压力的方程;之后将该方程与动量方程基于显式物质点法框架进行离散,并采用完全拉格朗日格式消除物质点跨网格产生的误差,提升大变形问题的模拟精度;对位移和压强场采用不同阶次的B样条插值函数并通过引入针对体积变形的重映射技术改进了算法,提升算法的准确性.此外,算法通过实施一种交错求解格式在每个时间步对位移场和压强场依次进行求解.最后,给出几个典型数值算例来验证本文所提出的混合格式TLMPM的有效性和准确性,计算结果表明该方法可以有效处理体积自锁,准确地模拟近似不可压软材料的大变形动力学行为.     

DISCUSSION ON THE APPLICATION OF NUMERICAL EXPERIMENT IN EXPERIMENTAL TEACHING OF FLUID MECHANICS

为强化学生对流体力学课程理论知识及实验教学环节的掌握,本文针对传统实验教学过程中受限于硬件设施、模型等问题,将数值模拟技术与流体力学实验教学相结合。以伯努利实验为例,应用数值模拟软件进行物理建模、网格划分、变参数数值模拟等实验内容,通过压力云图、流线等模拟结果对能量转换特性进行了分析,增强了学生对流体力学理论知识及相关参数的物理意义、变化特征等的理解,使学生在实验中产生强烈的“参与感”,有效地提高了教学水平和教学质量。     

Matlab求解理论力学问题系列(四)乒乓球的滚动

本文讨论了乒乓球运动的理论解和数值解。在考虑摩擦、滚动摩阻、空气阻力等条件的一般情况下,只能获得部分解析解,而采用数值计算的方法可获得系统所有的运动规律。文章介绍了如何利用Matlab处理分段函数问题,包括求解微分代数方程、利用积分中断条件确定分段的时间点,以及利用部分解析解验证数值计算的准确性等,最后通过大量计算获得乒乓球向前运动及向后运动的分界线,以及乒乓球正好回到出发位置的条件。     

粗糙颗粒动理学及稠密气固两相流动的数值模拟

考虑颗粒碰撞过程中摩擦作用,给出了粗糙颗粒碰撞动力学.引入颗粒相拟总温来表征颗粒平动和转动脉动能量的特征.基于气体分子运动论,建立颗粒碰撞中平动和旋转共同作用的粗糙颗粒动理学,给出了颗粒相压力和黏度等输运参数计算模型.运用基于颗粒动理学的欧拉-欧拉气固两相流模型,数值模拟了流化床内气体颗粒两相流动特性,分析了颗粒旋转流动对颗粒碰撞能量交换和耗散的影响.模拟得到的流化床内径向颗粒浓度和提升管内颗粒轴向速度与他人实验结果相吻合.模拟结果表明随着颗粒浓度的增加,颗粒相压力和能量耗散逐渐增加,而颗粒拟总温先增加后下降.随着颗粒粗糙度系数的增加,床内平均颗粒相拟总温和能量耗散增加,表明颗粒旋转产生的摩擦将导致颗粒旋转脉动能量的改变,影响床内气体-颗粒两相宏观流动特性.      

Mathematical modeling for dynamic stability of sandwich beam with variable mechanical properties of core

The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is formulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton’s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.     

高应变率下延性多孔介质中孔洞的动态演化

本文提出了一个新的材料延性动态损伤模型，模型中不但包括了率效应，同时还考虑了惯性效应，孔洞表面能变化和材料硬化对孔洞演化的影响。此外，在模型中同时考虑了体应力和偏应力对孔洞演化的作用，从孔洞演化方程地接到了孔洞增长和压缩应力临界表达式，Ｃａｒｒｏｌｌ和Ｈｏｌｔ结果作为该表达式的一个特例而得出，模型的数值分析得出以下结论：①延性孔洞的动太增长对率效应十分敏感，应变率越高，孔洞增长越快；②惯性效应在主