1问题一：什么是非等温试验？
通常有等温法(也称静态法)和非等温法(也称动态法)，
等温法是较早研究化学动力学时普遍采用的方法，该法的缺点在于比较费时，并且研究物质分解时，往往在升到一定的试验温度之前物质己发生初步分解，使得结果不很可靠。在非等温法中，试样温度随时间按线性变化，它在不同温度下的质量由热天平连续记录下来。非等温法是从反应开始到结束的整个温度范围内研究反应动力学，测得的一条热重曲线与不同温度下测得的多条等温失重曲线提供的数据等同，相比于等温法，非等温法只需一个微量的试验样品，消除了样品间的误差以及等温法将样品升至一定温度过程中出现的误差，并节省了试验时间。在目前的热重分析中常采用非等温法来进行动力学的研究。
问题二：文献中常用热解动力学表达式
d (a)/dt=kf(a) ——（1）
a为t时刻的分解率（材料的失重百分率）又称转化率。a=(m0-m)/(m0-m∞)
k=A exp(-E/RT)——（2）β=dT／dt ——（3）
采用coats-Readferm积分法推到
Ln[g(a)/T2]=ln(AR/βE)-E/RT f(a)=(1-a)2 f(a)为分饵的固体反应物与反应速率的函数关系。设Y= Ln[g(a)/T2] X=1/T 做X，Y直线曲线，求出斜率即可得到活化能E，同时得到结局求出指前因子A。
确定g(a)的值就能得到活化能E，常用g(a)的形式很多，有的是模型，有的是反应级数，总之尝试多种方法，找到最合适的，得到更精确的线性关系。
问题三：
1单条升温速率曲线的Coats-Redfern法，跟上述方程表达式一样，可得，
ln[－ln( 1 －a)/T 2] = ln[AR/βE( 1－2RT/ E) ]－E/RT( n = 1) ，（4）
ln[－( 1 －a)1 －n/T2( 1 －n ) ] = ln [AR/βE (1－2RT/ E) ]－E/RT( n≠1) . （5）
因为，一般活化能 E 的数值远大于温度T，所以(1?2RT/E)≈1，则式(4)和式(5)右端第1项几乎是常数。因此，可分别取n等于0.5, 0.6, 0.7, 0.8, 1.0, 1.2和1.5，结合热重实验的数据得到式(4)和式(5)的左端数值，并对1/T作图，得到这些直线的线性相关系数和标准误差数据，通过对比确定出线性较好的直线，由其斜率得到活化能E。
2，多条升温速率曲线的Flynn-Wall-Ozawa 法
Flynn-Wall-Ozawa(FWO)法通过多条升温速率曲线确定动力学参数，是等转化率法、积分法的一种。
根据式（1）（2）（3）进行移项积分得到，
Logβ=log[AE/RG(a)]-2.315-0.4567E/RT
由不同升温速率βi的TG 实验数据，在同一反应深度a下，找到相应的温度Ti，则lgβi 与Ti可以拟合得到一条直线，由其斜率可以得到活化能E，并且可以得到活化能随反应深度a的变化关系。（例如excel蒙古栎的四种升温速率）
计算流体动力学(简称CFD)是建立在经典流体动力学与数值计算方法基础之上的一门新型独立学科，通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。它兼有理论性和实践性的双重特点。 第一章节 流体流动现象大量存在于自然界及多种工程领域中，所有这些过程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。本章向读者介绍这些守恒定律的数学表达式，在此基础上提出数值求解这些基本方程的思想，阐述计算流体力学的任务及相关基础知识，最后简要介绍目前常用的计算流体动力学商用软件。 计算流体动力学((Computational Fluid Dynamics简称CFD)是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看做是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制卜对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 1.1.2计算流体动力学的工作步骤 采用CFD的方法对流体流动进行数值模拟，通常包括如下步骤: (1)建立反映工程问题或物理问题本质的数学模型。具体地说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数 学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 (2}}寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是c}}的核心。 (3})编制程序和进行计算。这部分工作包括计算网格划分、初始条件和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Na}ier-Stakes方程就是一个讨，分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲.数值模拟又叫数值试验。应该指出，这部分工作不是轻而易举就可以完成的。 4})显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。 以上这些步骤构成了CFD数值模拟的全过程。其中数学模型的建立是理论
机械系统动力学报告 题目：电梯机械系统的动态特性分析 姓名： 专业： 学号：
电梯机械系统的动态特性分析 一、课题背景介绍 随着社会的快速发展，城市人口密度越来越大，高层建筑不断涌现，因此，现在对电梯的提出了更高的要求，随着科技的进步，在满足客观需求的基础上，电梯向着舒适性，高速，高效的方向发展。在电梯的发展过程中，安全性和功能性一直是电梯公司首要考虑的因素，其中舒适性也要包含在电梯的设计中，避免出现速度或者加速度出现突变，或者电梯运行过程中的振动引起人们的不适。因此，在电梯的设计过程中，对电梯进行动态特性分析是十分必要的。 二、在MATLAB中编程、绘图。 通过同组小伙伴的努力，已经得到了该系统的简化模型与运动方程。因此进行编程： 该系统的微分方程：[][][]{}[]Q x k x c x M= + ? ? ? ? ? ? + ? ? ? ? ? ?? ? ? ，其中矩阵[M]、 [C]、[K]、[Q]都已知。 该系统的微分方程是一个二阶一元微分方程，在MATLAB中，提供有求解常微分方程数值解的函数，其中在MATLAB中常用的求微分方程数值解的有7个：ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb 。 ode是MATLAB专门用于解微分方程的功能函数。该求解器有变步长（variable-step）和定步长（fixed-step）两种类型。不同类型有着不同的求解器，其中ode45求解器属于变步长的一种，采用Runge-Kutta
算法；和他采用相同算法的变步长求解器还有ode23。 ode45表示采用四阶，五阶Runge-Kutta单步算法,截断误差为(Δx)^3。解决的是Nonstiff(非刚性)常微分方程。 ode45是解决数值解问题的首选方法，若长时间没结果，应该就是刚性的，可换用ode23试试。 Ode45函数调用形式如下：[T,Y]=ode45(odefun,tspan,y0) 相关参数介绍如下： 通过以上的了解，并对该微分方程进行变换与降阶，得出程序。MATLAB程序： （1）建立M函数文件来定义方程组如下： function dy=func(t,y) dy=zeros(10,1); dy(1)=y(2); dy(2)=1/1660*(-0.006*y(2)+0.003*y(4)-0.0006*y(10)-1.27*10^7*y(1)+1.27*10^7*y (3)+2.54*10^6*y(9)); dy(3)=y(4); dy(4)=1/1600*(+0.03*y(2)-0.007*y(4)+0.003*y(6)+1.27*10^7*y(1)-7.274*10^8*y(3 )+1.27*10^7*y(5)); dy(5)=y(6);
《计算流体动力学分析》学习报告 计算流体力学基础： 本章主要讲解流体动力学的核心思想以及流体动力学的控制方程。 1、计算流体动力学（Computational Fluid Dynamic ）基本思想：把原来在时间和空间上的连续的物理量，用一系列离散点上的变量值来代替，通过一定的原则和方式建立变量之间的代数方程式，求解之后获得变量的近似值。 2、CFD 控制方程： 质量守恒方程 0)·=?+??u t ρρ（ 动量守恒方程（Navier-Stokes 方程） Fz z y x z u w div t w F z y x y u v div t v F z y x x u u div t u zz zx zx y zy yy xy x zx yx xx +??+??+??+??-=+??+??+??+??+??-=+??+??+??+??+??-=+??τττρρρτττρρρτττρρρ)()()()()()( 能量守恒方程 T p S gradT c k div T u div t +=+??)()(T ( ρρ） S T 为粘性耗散项。 方程含有u ，v ，w ，p ，T 和ρ六个未知量，所以还需要一个方程组，才能使其封闭，而这个方程组就是联系P 和ρ的状态方程组：P=（ρ，T ）。 组分质量守恒方程（在一个系统中，可能存在质的交换，或者存在化学组分时使用。） ()s s s s S c grad D div c u div t +=+??)()(c (s ρρρ ） 为便于对控制方程进行计算和分析，对CFD 控制方程写成通用格式： ()S z z y y x x z w y v x u t S grad div u div t +??Γ??+??Γ??+??Γ??=??+??+??+??+Γ=+??)()()()()()())()(φφφφρφρφρρφφφρρφ 依次为瞬态项，对流项，扩散项和源项。 3、湍流控制方程 三维的N-S 方程无论对于层流还是湍流都是是使用的，但由于直接求解三维瞬态的控制方程，对计算机的内存和速度要求很高，因此在工程上广为采用的方法是对瞬态的N-S 方程进行实践平均处理，同时补充反应湍流特性的其他方程，例如湍动能方程以及湍流耗散率方程
第二章：动力学系统的微分方程模型 利用计算机进行仿真时，一般情况下要给出系统的数学模型，因此有必要掌握一定的建立数学模型的方法。在动力学系统中，大多数情况下可以使用微分方程来表示系统的动态特性，也可以通过微分方程可以将原来的系统简化为状态方程或者差分方程模型等。在这一章中，重点介绍建系统动态问题的微分方程的基本理论和方法。 在实际工程中，一般把系统分为两种类型，一是连续系统；其数学模型一般是高阶微分方程；另一种是离散系统，它的数学模型是差分方程。 §2.1 动力学系统统基本元件 任何机械系统都是由机械元件组成的，在机械系统中有3种类型的基本机械元件：惯性元件、弹性元件和阻尼元件。 1 惯性元件：惯性元件是指具有质量或转动惯量的元件，惯量可以定义为使加速度（或角加速度）产生单位变化所需要的力（或力矩）。 惯量（质量）= ） 加速度（力（2 /) s m N 惯量（转动惯量）= ） 角加速度（力矩（2/) s rad m N ? 2 弹性元件：它在外力或外力偶作用下可以产生变形的元件，这种元件可以通过外力做功来储存能量。按变形性质可以分为线性元件和非线性元件，通常等效成一弹簧来表示。 对于线性弹簧元件，弹簧中所受到的力与位移成正比，比例常数为弹簧刚度k 。 x k F ?= 这里k 称为弹簧刚度，x ?是弹簧相对于原长的变形量，弹性力的方向总是指向弹 簧的原长位移，出了弹簧和受力之间是线性关系以外，还有所谓硬弹簧和软弹簧，它们的受力和弹簧变形之间的关系是一非线性关系。 3 阻尼元件：这种元件是以吸收能量以其它形式消耗能量，而不储存能量，可以形象的表示为一个活塞在一个充满流体介质的油缸中运动。阻尼力通常表示为： α x c R = 阻尼力的方向总是速度方向相反。当1=α，为线性阻尼模型。否则为非线性阻 尼模型。应注意当α等于偶数情况时，要将阻尼力表示为： ||1--=αx x c R 这里的“-”表示与速度方向相反
CFD 流体动力学软件介绍 CFD—计算流体动力学，因历史原因，国一直称之为计算流体力学。其结构为： 提出问题—流动性质（流、外流；层流、湍流；单相流、多相流；可压、不可压等等），流体属性（牛顿流体：液体、单组分气体、多组分气体、化学反应气体；非牛顿流体） 分析问题—建模—N-S方程（连续性假设），Boltzmann方程（稀薄气体流动），各类本构方程与封闭模型。 解决问题—差分格式的构造/选择，程序的具体编写/软件的选用，后处理的完成。 成果说明—形成文字，提交报告，赚取应得的回报。 CFD实现过程： 1.建模——物理空间到计算空间的映射。 主要软件： 二维： AutoCAD： 大家不要小看它，非常有用。一般的网格生成软件建模都是它这个思路，很少有参数化建模的。相比之下AutoCAD的优点在于精度高，草图处理灵活。可以这样说，任何一个网格生成软件自带的建模工具都是非参数化的，而对于非参数化建模来说，AutoCAD应该说是最好的，毕竟它发展了很多很多年！ 三维： CATIA：航空航天界CAD的老大，法国人的东西，NB，实体建模厉害，曲面建模独步武林。本身可以生成有限元网格，前几天又发布了支持ICEM-CFD的插件ICEM-CFD CAA V5。有了它和ICEM-CFD，可以做任何建模与网格划分！ UG：总觉得EDS脑袋进水了，收了I-deas这么久了，也才发布个几百M的UG NX 2.0，还被大家争论来争论去说它如何的不好用！其实，软件本身不错，大公司用得也多，可是就这么打市场，早晚是走下坡路。按CAD建模的功能来说它排不上第一，也不能屈居第二，尤其是加上了I－DEAS更是如虎添翼。现
微分形式的连续性方程
连续方程是流体力学的基本方程之一，流体运动的连续方程，反映流体运动和流体质量分布的关系，它是在质量守恒定律在流体力学中的应用。 重点讨论不同表现形式的流体连续方程。
用一个微六面体元控制体建立微分形式的连续性方程。 设在流场中取一固定不动的微平行六面体（控制体），在直角坐标系oxyz 中，六面体的边长取为dx ，dy ，dz 。 先看x 轴方向的流动，流体从ABCD 面流入六面体，从EFGH 面流出。 在x 轴方向流出与流入质量之差 ()()[]x x x x u u u dx dydzdt u dydzdt dxdydzdt x x ρρρρ??+-=??
用同样的方法，可得在y 轴方向和z 轴方向的流出与流入 质量之差分别为 ()y u dxdydzdt y ρ??() z u dxdydzdt z ρ??这样，在dt 时间内通过六面体的全部六个面净流出的质量为： ()()()[]y x z u u u dxdydzdt x x x ρρρ???++???
在dt 的时间内，六面体内的质量减少了 ， 根据质量守恒定律，净流出六面体的质量必等于六面体内所减少的质量 ()dxdydzdt t ρ?-?()()()[]y x z u u u dxdydzdt dxdydzdt x y z t ρρρρ ????++=-????()()()0y x z u u u x y z t ρρρρ ????+++=????这就是直角坐标系中流体运动的微分形式的连续性方程。 这就是直角坐标系中流体运动的微分形式的连续性方程。 代表单位时间内，单位体积的质量变化 代表单位时间内，单位体积内质量的净流出
课后习题 第一章 1.计算流体动力学的基本任务是什么 计算流体动力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。 2.什么叫控制方程？常用的控制方程有哪几个？各用在什么场合？ 流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。如果流动包含有不同组分的混合或相互作用，系统还要遵守组分守恒定律。如果流动处于湍流状态，系统还要遵守附加的湍流输运方程。控制方程是这些守恒定律的数学描述。 常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组分质量守恒方程。质量守恒方程和动量守恒方程任何流动问题都必须满足，能量守恒定律是包含有热交换的流动系统必须满足的基本定律。组分质量守恒方程，在一个特定的系统中，可能存在质的交换，或者存在多种化学组分，每种组分都需要遵守组分质量守恒定律。 4.研究控制方程通用形式的意义何在？请分析控制方程通用形式中各项的意义。 建立控制方程通用形式是为了便于对各控制方程进行分析，并用同一程序对各控制方程进行求解。
各项依次为瞬态项、对流项、扩散项、源项。 6.CFD商用软件与用户自行设计的CFD程序相比，各有何优势？常用的商用CFD软件有哪些？特点如何？ 由于CFD的复杂性及计算机软硬件条件的多样性，用户各自的应用程序往往缺乏通用性。 CFD商用软件的特点是 功能比较全面、适用性强。 具有比较易用的前后处理系统和其他CAD及CFD软件的接口能力，便于用户快速完成造型、网格划分等工作。 具有比较完备的容错机制和操作界面，稳定性高。 可在多种计算机、多种操作系统，包括并行环境下运行。 常用的商用CFD软件有PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。PHOENICS除了通用CFD软件应该拥有的功能外，PHOENICS软件有自己独特的功能：开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。 CFX除了可以使用有限体积法外，还采用基于有限元的有限体积法。用于模拟流体流动、传热、多相流、化学反应、燃烧问题。其优势在于处理流动物理现象简单而几何形状复杂的问题。 SRAR-CD基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
非流体、热动专业CFD新手入门 首先掌握流体力学基本原理，丁祖荣主编的流体力学这本教材，仔细看两天，这样就会知道gambit中为什么会有边界层设置，边界层厚度如何设置；雷诺系数如何确定来判断层流与湍流；马赫数如何确定来判断流体是可压还是不可压，这样就能解决Fluent，是基于压力还是基于密度求解。能够对实际中一些看似简单的流体现象有深刻的认识，能够准确判断是定常流还是非定常流。 CFD网格划分 网格划分对于初学者所接触案例，其实非常简单。但实际工程中，大项目，特别涉及到整套工程，如环保，飞机，网格质量与数量都要求非常高，往往服务器类的PC才能解决问题，所谓的内存128G，CPU四核主频3.0以上。初学者，简单的管道，一般的机器还是没问题。有机械三维软件基础的，对于gambit建模就非常容易了。往往大项目，复杂的结构gambit 建模显得力不从心，所以对于流体工作者来说，学习三维软件对于建模有莫大的帮助，如Proe。 1.1Gambit介绍 网格的划分使用Gambit软件，首先要启动Gambit，在Dos下输入Gambit
第三章 流体动力学基础 习 题 一、单选题 1、在稳定流动中，在任一点处速度矢量是恒定不变的，那么流体质点是 （ ） A ．加速运动 B ．减速运动 C ．匀速运动 D ．不能确定 2、血管中血液流动的流量受血管内径影响很大。如果血管内径减少一半，其血液的流量将变为原来的（ ）倍。 A ．21 B ．41 C ．81 D ．161 3、人在静息状态时，整个心动周期内主动脉血流平均速度为0.2 m/s ，其内径d =2×10-2 m ，已知血液的粘度η =×10-3 Pa·S，密度ρ=×103 kg/m 3 ，则此时主动脉中血液的流动形态处于（ ）状态。 A ．层流 B ．湍流 C ．层流或湍流 D ．无法确定 4、正常情况下，人的小动脉半径约为3mm ，血液的平均速度为20cm/s ，若小动脉某部分被一硬斑阻塞使之变窄，半径变为2mm ，则此段的平均流速为（ ）m/s 。 A ．30 B ．40 C ．45 D ．60 5、有水在同一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为 S B =5cm 2，A 、B 两点压强差为1500Pa ，则A 处的流速为（ ）。 A ．1m/s B ．2m/s C ．3 m/s D ．4 m/s 6、有水在一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为S B =5cm 2 ，A 、B 两点压强之差为1500Pa ，则管道中的体积流量为（ ）。 A ．1×10-3 m 3 /s B ．2×10-3 m 3 /s C ．1×10-4 m 3 /s D ．2×10-4 m 3 /s 7、通常情况下，人的小动脉内径约为6mm ，血流的平均流速为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，测得此处血流的平均流速为80cm/s ，则小动脉此处的内径应为（ ）mm 。 A ．4 B ．3 C ．2 D ．1 8、正常情况下，人的血液密度为×103 kg/m 3 ，血液在内径为6mm 的小动脉中流动的平均速度为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，此处内径为4mm ，则小动脉宽处与窄处压强之差（ ）Pa 。 二、判断题
1问题一：什么是非等温试验？ 通常有等温法(也称静态法)和非等温法(也称动态法)， 等温法是较早研究化学动力学时普遍采用的方法，该法的缺点在于比较费时，并且研究物质分解时，往往在升到一定的试验温度之前物质己发生初步分解，使得结果不很可靠。在非等温法中，试样温度随时间按线性变化，它在不同温度下的质量由热天平连续记录下来。非等温法是从反应开始到结束的整个温度范围内研究反应动力学，测得的一条热重曲线与不同温度下测得的多条等温失重曲线提供的数据等同，相比于等温法，非等温法只需一个微量的试验样品，消除了样品间的误差以及等温法将样品升至一定温度过程中出现的误差，并节省了试验时间。在目前的热重分析中常采用非等温法来进行动力学的研究。 问题二：文献中常用热解动力学表达式 d (a)/dt=kf(a) ——（1） a为t时刻的分解率（材料的失重百分率）又称转化率。a=(m0-m)/(m0-m∞) k=A exp(-E/RT)——（2）β=dT／dt ——（3） 采用coats-Readferm积分法推到 Ln[g(a)/T2]=ln(AR/βE)-E/RT f(a)=(1-a)2 f(a)为分饵的固体反应物与反应速率的函数关系。设Y= Ln[g(a)/T2] X=1/T 做X，Y直线曲线，求出斜率即可得到活化能E，同时得到结局求出指前因子A。 确定g(a)的值就能得到活化能E，常用g(a)的形式很多，有的是模型，有的是反应级数，总之尝试多种方法，找到最合适的，得到更精确的线性关系。 问题三： 1单条升温速率曲线的Coats-Redfern法，跟上述方程表达式一样，可得， ln[－ln( 1 －a)/T 2] = ln[AR/βE( 1－2RT/ E) ]－E/RT( n = 1) ，（4） ln[－( 1 －a)1 －n/T2( 1 －n ) ] = ln [AR/βE (1－2RT/ E) ]－E/RT( n≠1) . （5） 因为，一般活化能 E 的数值远大于温度T，所以(1?2RT/E)≈1，则式(4)和式(5)右端第1项几乎是常数。因此，可分别取n等于0.5, 0.6, 0.7, 0.8, 1.0, 1.2和1.5，结合热重实验的数据得到式(4)和式(5)的左端数值，并对1/T作图，得到这些直线的线性相关系数和标准误差数据，通过对比确定出线性较好的直线，由其斜率得到活化能E。 2，多条升温速率曲线的Flynn-Wall-Ozawa 法 Flynn-Wall-Ozawa(FWO)法通过多条升温速率曲线确定动力学参数，是等转化率法、积分法的一种。 根据式（1）（2）（3）进行移项积分得到， Logβ=log[AE/RG(a)]-2.315-0.4567E/RT 由不同升温速率βi的TG 实验数据，在同一反应深度a下，找到相应的温度Ti，则lgβi 与Ti可以拟合得到一条直线，由其斜率可以得到活化能E，并且可以得到活化能随反应深度a的变化关系。（例如excel蒙古栎的四种升温速率）
第10章系统动力学模型 系统动力学模型(System Dynamic)是社会、经济、规划、军事等许多领域进行战略研究的重要工具，如同物理实验室、化学实验室一样，也被称之为战略研究实验室，自从问世以来，可以说是硕果累累。 1 系统动力学概述 2 系统动力学的基础知识 3 系统动力学模型 第1节系统动力学概述 概念 系统动力学是一门分析研究复杂反馈系统动态行为的系统科学方法，它是系统科学的一个分支，也是一门沟通自然科学和社会科学领域的横向学科，实质上就是分析研究复杂反馈大系统的计算仿真方法。 系统动力学模型是指以系统动力学的理论与方法为指导，建立用以研究复杂地理系统动态行为的计算机仿真模型体系，其主要含义如下： 1 系统动力学模型的理论基础是系统动力学的理论和方法； 2 系统动力学模型的研究对象是复杂反馈大系统； 3 系统动力学模型的研究内容是社会经济系统发展的战略与决策问题，故称之为计算机仿真法的“战略与策略实验室”； 4 系统动力学模型的研究方法是计算机仿真实验法，但要有计算机仿真语言DYNAMIC的支持，如：PD PLUS，VENSIM等的支持；
5 系统动力学模型的关键任务是建立系统动力学模型体系； 6 系统动力学模型的最终目的是社会经济系统中的战略与策略决策问题计算机仿真实验结果，即坐标图象和二维报表； 系统动力学模型建立的一般步骤是：明确问题，绘制因果关系图，绘制系统动力学模型流图，建立系统动力学模型，仿真实验，检验或修改模型或参数，战略分析与决策。 地理系统也是一个复杂的动态系统，因此，许多地理学者认为应用系统动力学进行地理研究将有极大潜力，并积极开展了区域发展，城市发展，环境规划等方面的推广应用工作，因此，各类地理系统动力学模型即应运而生。 发展概况 系统动力学是在20世纪50年代末由美国麻省理工学院史隆管理学院教授福雷斯特（）提出来的。目前，风靡全世界，成为社会科学重要实验手段，它已广泛应用于社会经济管理科技和生态灯各个领域。福雷斯特教授及其助手运用系统动力学方法对全球问题，城市发展，企业管理等领域进行了卓有成效的研究，接连发表了《工业动力学》，《城市动力学》，《世界动力学》，《增长的极限》等着作，引起了世界各国政府和科学家的普遍关注。 在我国关于系统动力学方面的研究始于1980年，后来，陆续做了大量的工作，主要表现如下： 1）人才培养 自从1980年以来，我国非常重视系统动力学人才的培养，主要
1738年瑞士数学家：伯努利在名著《流体动力学》中提出了伯努利方程。 1755年欧拉在名著《流体运动的一般原理》中提出理想流体概念，并建立了理想流体基本方程和连续方程，从而提出了流体运动的解析方法，同时提出了速度势的概念。 1781年拉格朗日首先引进了流函数的概念。 1826年法国工程师纳维，1845年英国数学家、物理学家斯托克思提出了著名的N-S方程。 1876年雷诺发现了流体流动的两种流态：层流和紊流。 1858年亥姆霍兹指出了理想流体中旋涡的许多基本性质及旋涡运动理论，并于1887年提出了脱体绕流理论。 19世纪末，相似理论提出，实验和理论分析相结合。 1904年普朗特提出了边界层理论。 20世纪60年代以后，计算流体力学得到了迅速的发展。流体力学内涵不断地得到了充实与提高。 理想势流伯努利方程 （3-14） 或（3-15） 物理意义:在同一恒定不可压缩流体重力势流中，理想流体各点的总比能相等即在整个势流场中，伯努利常数C 均相等。 （应用条件：“”所示） 符号说明 二、沿流线的积分
1.只有重力作用的不可压缩恒定流，有 2.恒定流中流线与迹线重合： 沿流线（或元流）的能量方程： （3-16） 注意：积分常数C，在非粘性、不可压缩恒定流流动中，沿同一流线保持不变。一般不同流线各不相同（有旋流）。 （应用条件：“”所示，可以是有旋流） 流速势函数（势函数）观看录像>> ?存在条件：不可压缩无旋流，即或 必要条件存在全微分d 直角坐标 （3-19） 式中：——无旋运动的流速势函数，简称势函数。 ?势函数的拉普拉斯方程形式 对于不可压缩的平面流体流动中，将（3-19）式代入连续性微分方程（3-18），有： 或（3-20） 适用条件：不可压缩流体的有势流动。 点击这里练习一下 极坐标 （3-21） 流函数
第9章 质点动力学基本方程 一、是非题(正确的在括号内打“√”、错误的打“×”) 1. 凡是适合于牛顿三定律的坐标系称为惯性参考系。 ( √ ) 2. 一质点仅受重力作用在空间运动时，一定是直线运动。 ( × ) 3. 两个质量相同的物体，若所受的力完全相同，则其运动规律也相同。 ( × ) 4. 质点的运动不仅与其所受的力有关，而且还和运动的初始条件有关。 ( √ ) 5. 凡运动的质点一定受力的作用。 ( × ) 6. 质点的运动方向与作用于质点上的合力方向相同。 ( × ) 二、填空题 1．质点是指大小可以忽略不计，但具有一定质量的物体。 — 2．质点动力学的基本方程是∑= i m F a ，写成自然坐标投影形式为∑=τF dt s d m 2 2 ∑= n F v m ρ 2 ∑ =b F 0。 3．质点保持其原有运动状态不变的属性称为惯性。 4．质量为m 的质点沿直线运动，其运动规律为0ln(1)v t x b b =+，其中0v 为初速度，b 为常数。则作用于质点上的力=F 20 2 0() mbv b v t - +。 5．飞机以匀速v 在铅直平面内沿半径为r 的大圆弧飞行。飞行员体重为P ，则飞行员对座椅的最大压力为2 (1)v P gr +。 三、选择题 1．如图所示，质量为m 的物块A 放在升降机上， 当升降机以加速度a 向上运动时，物块对地板的压力等于( B )。 (A) mg (B) )(a g m + (C) )(a g m - (D) 0 2．如图所示一质量弹簧系统，已知物块的质量为m ，弹簧的刚度系数为c ，静伸长量为s δ，原长为0l ，若以弹簧未伸长的下端为坐标原点，则物块的运动微分方程可写成( B )。 ， (A) 0=+x m c x (B) 0)(=-+s x m c x δ 、 、
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
习题13－1图 *第13章 动力学普遍方程 和第二类拉格朗日方程 13－1 图示均质细杆OA 长为l ，重力为P ，在重力作用下可在铅垂平面内摆动，滑块O 质量不计，斜面倾角θ，略去各处摩擦，若取x 及?为广义坐标，试求对应于x 和?的广义力。 解：应用几何法，令0δ=x ；0δ≠? 则：??????sin 2 1δδ2sin δδPl l P W Q -=-='= 令0δ≠x ；0δ=? 则：θθsin δδsin δδP x x P x W Q x -=-=''= 13－2 图示在水平面内运动的行星齿轮机构，已知固定齿轮半径为R ，均质行星齿轮半径为r ，质量为m ，均质杆OA 质量为m 1，杆受矩为M 的力偶作用而运动，若取?为广义坐标，试求相应的广义力。 解：应用几何法，设对应于?的虚位移0δ≠? 则： M M W Q === ? ???δδδδ 13－3 在图示系统中，已知：均质圆柱A 的质量为M 、半径为R ，物块B 的质量为m ，光滑斜面的倾角为β，滑轮质量忽略不计，并假设斜绳段平行斜面。若以θ 和y 为广义坐标，试分别用动力学普遍方程和第二类拉格朗日方程求： （1）系统运动微分方程； （2）圆柱A 的角加速度和物块B 的加速度。 解：（1）在系统上施加惯性力如图（a ）所示。 其中：)(I θ R y M F A -=；y m F B =I θ θ 2I 2 1MR J M A A == 应用动力学普遍方程， δ)sin (δ)sin (I I I I =+-+-- -θββR Mg M R F y Mg F F mg A A A B 可得系统运动微分方程： 0sin )(=----βθMg R y M y m mg 0sin 2 1)(2=+--R Mg MR R R y M βθθ 整理后有：0)sin ()(=-+-+g m M MR y M m βθ 0sin 2 3=--βθg y R 习题13-2图 习题13-3图 F I
计算流体动力学概述 作者：王福军 1 什么是计算流体动力学 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为：把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值CFD可以看做是在流动基本方程（质量守恒方程飞动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（如速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。CFD方法与传统的理论分析方法、实验测量方法组成了研究流体流动问题的完整体系，图1给出了表征三者之间关系的“三维”流体力学示意图理论分析方法的优点在于所得结果具有普遍性，各种影响因素清晰可见，是指导实验研究和验证新的数值计算方法的理论基础。但是，它往往要求对计算对象进行抽象和简化，才有可能得出理论解。对于非线性情况，只有少数流动才能给出解析结果。 “三维”流体力学示意图 实验测量方法所得到的实验结果真实可信，它是理论分析和数值方法的基础，其重要性不容低估。然而，实验往往受到模型尺寸、流场扰动、人身安全和测量精度的限制，有时可能很难通过试验力一法得到结果。此外，实验还会遇到经费投入、人力和物力的巨大耗费及周期长等许多困难。 而CFD方法恰好克服了前面两种方法的弱点，在计算机上实现一个特定的计算。就好像在
第十章 质点动力学基本方程 10-3 半径为R 的偏心轮绕O 轴以匀角速度ω转动，推动导板沿铅直轨道运动，如图所示。导板顶部放有一质量为m 的物块A ，设偏心距e OC =，开始时OC 沿水平线。求：（1）物块对导板的最大压力；（2）使物块不离开导板的ω最大值。 解：建立如图所示直角坐标系Oxy ，导板与物块均沿y 轴线作直 线运动，导板作平动，其运动规律为 t e R y ωsin += 对时间求二阶导数得 t e a y ωωsin 2-= 物块A 受重力m g 和导板的约束反力N F 作用如图)a (。 物块对导板的压力与N F 等值、反向、共线。由图(a)得物块A 的运动微分方程在y 轴的投影式为 ) sin (2N N t e g m F ma mg F y ωω-==- 1）物块对导板的最大压力 )(2N ωe g m F += 2）要使物块不离开导板，则应有 0)(2min N ≥-=ωe g m F 即 2ωe g ≥ 故 e g =max ω 10-7 销钉M 的质量为0.2 kg ，水平槽杆带动，使其在半径为mm 200=r 的固定半圆槽内运动。设水平槽杆以匀速mm/s 400=v 向上运动，不计摩擦。求在图示位置时圆槽对销钉M 的作用力。 解：以水平槽为动系，速度分析如图)a (，v v =e 3 24.02 330cos e a ?==?=v v v 受力与加速度分析如图(b)， 2222a n m/s 07.132.044.04 3=??=?==r v r v a M r t n a a a =+M M 向铅直方向投影，得 2t n 2 n t t n m/s 23.13079.09238.030sin 30cos m/s 616.03 30cos 30sin =+=?+?====?-?M M Mx M M M M a a a a a a a 设水平槽对M 的反力为F N ，圆槽对M 的反力为F ，则
2012～2013学年第1学期 12级研究生《计算流体力学》结课作业 适用专业：供热供燃气通风及空调工程 一、结合某一具体学科，阐述纯理论方法、实验方法及数值方法在科学研究中的各自优缺点，在此基础上论述数值模拟方法的发展前景。（不少于4千字）。 流体力学是力学的一个重要分支, 是研究流体(液体和气体)的力学运动规律及其应用的学科, 主要研究在各种力的作用下，流体本身的静止状态和运动状态特征，以及流体和相邻固体界面有相对运动时的相互作用和流动规律。在人们的生活和生产活动中随时随地都可遇到流体，流体力学与人类的日常生活和生产事业密切相关。按其研究内容的侧重点不同，分为理论流体力学和工程流体力学。其中理论流体力学主要采用严密的数学推理方法，力求准确性和严密性，工程流体力学侧重于解决工程实际中出现的问题，而不追求数学上的严密性。当然由于流体力学研究的复杂性，在一定程度上，两种方法都必须借助于实验研究，得出经验或半经验的公式。 在实际工程的诸多领域流体力学都起着十分重要的作用。如气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，都广泛地用到流体力学知识。许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了流体力学自身的不断发展。1950年后，计算机的发展给予流体力学以极大的推动作用。 目前，解决流体力学问题的方法主要有实验方法、理论分析方法和数值方法三种。 实验方法 同物理学、化学等学科一样，流体力学的研究离不开实验，尤其是对新的流体运动现象的研究。实验能显示运动特点及其主要趋势，有助于形成概念，检验理论的正确性。二百年来流体力学发展史中每一项重大进展都离不开实验。流体力学实验研究方法有实物实验、比拟研究和模型研究三类：实物实验是用仪器实测原型系统的流动参数，适用于较小的原型；比拟实验是利用电场和磁场来模拟流场，实施起来限制条件较多；模型研究是实验流体力学最常用的研究方法。 实验研究的一般过程是：在相似理论的指导下建立实验模型，用流体测量技术测量流动参数，处理和分析实验数据。建立实验模型要求模型与原型满足相似理论，即满足两个流场
机械系统动力学报告 学院：机械工程学院 专业：机械电子工程 姓名： 学号：
机械系统动力学 1 机械系统动力学简介 随着现代工业对机械设备及机械传动系统的要求越来越高，机械设备及机械传动系统向着大型化、高速化、轻量化、构件柔性化方向发展。人们对生产率的不断追求，使得机械的运转速度不断提高；与此同时，人们总是希望使用的机器轻巧一些，材质的改善使得构件的截面可以设计得更小一些，这样就减轻了重量、节省了材料；速度高了使得机器中的惯性力增大，截面小了使得构件的柔性加大，这样使得系统更容易产生振动，振动降低了机械的精度和寿命，恶化了劳动条件。 由于动力学研究的复杂性，人们常常引入一些假定，使问题的研究过程简化。随着生产实践的发展对动力学分析的准确度提出了新的要求；而科学技术的发展又为动力学分析提供了新的理论和分析手段。动力学的发展趋势是：逐步将这些假定抛弃，使得分析更接近客观实际。 对于低速机械，运动中产生的惯性力可以忽略不计。随着机械速度的提高，惯性力不能再被忽略，此时可根据达朗伯原理将惯性力加入静平衡方程进行求解，这种方法就称为动态静力方法。 为了求出惯性力，就必须知道构件的加速度。因此在动态静力分析之前首先要进行运动学分析，而在运动学分析时总是假定构件是按某一给定的理想运动规律运动，多数驱动构件均被假定做等速回转运动。由于采用了等速回转这一假定，在动态静力分析中便不涉及原动机的特性，因而，着本质上是一种理想运动状态下的力学分析。现在在许多速度较高的机械中，用动态静力分析代替了静力分析。 在力的作用下，机械很难维持“驱动件等速回转”这种假定。尽管这种假定在许多情形下是允许的，但在实际运动中常常需要知道系统的真实运动规律，因而进行动力分析就是求出在外力作用下系统的真实运动，用于解决动力学的正问题。由于分析的对象是整个机械系统，所以又称为机械系统动力学。 在高速情况下，动态精度与静态精度有很大的区别。精密机床的动态性能研究、高速间歇机构的动态定位精度研究就是这样发展起来的。 高速旋转机械可以采用静态设计，制造出来后通过动平衡减少振动，还要使运转速度避开临界转速。但是，随着转速的进一步提高和柔性转子的出现，就必须采取全方位的综合措施，不仅在设计时要进行认真的动力学分析，而且在运行过程中
目录 1 刚体系统 (1) 2 弹性系统动力学 (6) 3 高速旋转体动力学 (10)
1 刚体系统 一般力学研究的对象，是由两个或两个以上刚体通过铰链等约束联系在一起的力学系统，为一般力学研究对象。自行车、万向支架陀螺仪通常可看成多刚体系统。人体在某种意义上也可简化为一个多刚体系统。现代航天器、机器人、人体和仿生学中关于动物运动规律的研究都提出了多刚体系统的一系列理论模型作为研究对象。多刚体系统按其内部联系的拓扑结构，分为树型和非树型（包含有闭链）；按其同外界的联系情况，则有有根和无根之别。利用图论的工具可以一般地分析多刚体系统的构造，建立系统的数学模型和动力学方程组。也可从分析力学中的高斯原理出发，用求极值的优化算法直接求解系统的运动和铰链反力。依照多刚体系统动力学的理论和方法，广泛采用电子计算机对这些模型进行研究，对于精确地掌握这些对象的运动规律是很有价值的。 1.1 自由物体的变分运动方程 任意一个刚体构件i ，质量为i m ，对质心的极转动惯量为i J '，设作用于刚体的所有外力向质心简化后得到外力矢量i F 和力矩i n ，若定义刚体连体坐标系y o x '''的原点o '位于刚体质心，则可根据牛顿定理导出该刚体带质心坐标的变分运动方程： 0][][=-'+-i i i i i i i T i n J F r m r φδφδ&&&& (1-1) 其中，i r 为固定于刚体质心的连体坐标系原点o '的代数矢量，i φ为连体坐标系相对于全局坐标系的转角，i r δ与i δφ分别为i r 与i φ的变分。 定义广义坐标： T i T i i r q ],[φ= (1-2) 广义： T i T i i n F Q ],[= (1-3) 及质量矩阵： ),,(i i i i J m m diag M '= (1-4) 体坐标系原点固定于刚体质心时用广义力表示的刚体变分运动方程：