积分方法和计算步长的选择
- 格式:ppt
- 大小:153.51 KB
- 文档页数:16
控制系统数字仿真题库一、填空题1. 定义一个系统时,首先要确定系统的边界;边界确定了系统的范围,边界以外对系统的作用称为系统的输入,系统对边界以为环境的作用称为系统的输出。
2.系统的三大要素为:实体、属性和活动。
3.人们描述系统的常见术语为:实体、属性、事件和活动。
4.人们经常把系统分成四类,它们分别为:连续系统、离散系统、采样数据系统和离散-连续系统。
5、根据系统的属性可以将系统分成两大类:工程系统和非工程系统。
6.根据描述方法不同,离散系统可以分为:离散时间系统和离散事件系统。
7. 系统是指相互联系又相互作用的实体的有机组合。
8.根据模型的表达形式,模型可以分为物理模型和数学模型二大类,其中数学模型根据数学表达形式的不同可分为二种,分别为:静态模型和动态模型。
9、采用一定比例按照真实系统的样子制作的模型称为物理模型,用数学表达式来描述系统内在规律的模型称为数学模型。
10.静态模型的数学表达形式一般是代数方程和逻辑关系表达式等,而动态模型的数学表达形式一般是微分方程和差分方程。
11.系统模型根据描述变量的函数关系可以分类为线性模型和非线性模型。
12 仿真模型的校核是指检验数字仿真模型和数学模型是否一致。
13.仿真模型的验证是指检验数字仿真模型和实际系统是否一致。
14.计算机仿真的三个要素为:系统、模型与计算机。
15.系统仿真的三个基本活动是系统建模、仿真建模和仿真试验。
16.系统仿真根据模型种类的不同可分为:物理仿真、数学仿真和数学-物理混合仿真。
17.根据仿真应用目的的不同,人们经常把计算机仿真应用分为四类,分别为:系统分析、系统设计、理论验证和人员训练。
18.计算机仿真是指将模型在计算机上进行实验的过程。
19. 仿真依据的基本原则是:相似原理。
20. 连续系统仿真中常见的一对矛盾为计算速度和计算精度。
21.保持器是一种将离散时间信号恢复成连续信号的装置。
22.零阶保持器能较好地再现阶跃信号。
第三章 化学反应动力学的计算化学反应的速度各不相同,有的反应速度极快,只要几个毫微秒就达到平衡(接近扩散速度,如无机酸碱中和),有的反应速度极慢,几乎看不到变化(如自然界的某些变化)。
大部分有机化学反应可用常规方法测量,对某些快速反应则可用停留法、驰豫法等测量。
不论反应速度的快慢,动力学方程都是类似的。
一、化学反应动力学方程反应物浓度随时间的变化绝大部分不是线性关系,而是一条曲线,见图3-1。
反应速度公式可用微分方程来表示。
具有简单级数的化学反应的反应速度公式可用积分式表示:一级 如:0AA1Adc A C =a, -=k c dt 生成物:,㏑C A =㏑a –K 1t 二级 A+A →产物 C A 0=a 2A 2A 2A d c 11-k C , =+k t d t c a对于反应 1-1k k A B 这一可逆反应初始条件 t=0 a 0 时间t 时 t=t a-x x达到平衡时,B 的浓度为X e ,则可逆反应的速度积分式为: 级数:1-1 1-10k A A e e 1A -1B k 0e 0C =a dc x xA B=-k C +k C : =kt dt a x -xC =0ln 1-21-10Ak0A e e e B 1A -1B C k e e 0CC =a dc x ax +x(a-x )A B+C C =0=-k C +k C C : =kt dt 2a-x a(x -x)C =0ln 二、常微分方程的解化学反应动力学方程是用微分方程表示的,对于简单的反应,可直接求得微分方程的解。
微分方程:()(1)(,,,......)......(1)n n y f x y y y -'=在区间a<x<b 的解,是指()y x ϕ=,这样一个函数,在所述区间内存在导数()(),(),......()n x x x ϕϕϕ'''。
且对于区间a<x<b 内的每一个x ,等式(1)都成立。
数值积分方法数值积分,又称为数值分析,是一种应用科学和数学技术来求解数学分析中几何或者微分方程的数学方法。
在实际应用中,有一系列的数值积分方法可以应用于解决某些数学问题,其中包括这些方法的微元法、有限元法、线性多项式插值法、指数插值法、函数拟合法和通用积分等方法。
通过合理的数值技术及其应用,可以有效地解决众多实际问题。
数值积分是数值分析中最基本的方法,指将数学分析中的连续函数或曲线所表示的求和问题离散化,以使其被数值计算机计算出来,也被称为数值积分。
当需要用数值积分方法求某函数的定积分时,首先必须找出该函数的积分表达式,然后对该表达式进行离散化,得到计算机可以处理的函数,最后根据具体的算法,得到数值积分的解。
数值积分方法具有多种形式,分别适用于不同实际问题。
首先,常用的数值积分方法有积分公式,如梯形公式、抛物线公式、Simpson 公式等,以及牛顿-拉夫逊多项式插值公式等,这些积分公式可以以直接的方式计算定积分,但是这种方法只适用于简单的定积分计算,在复杂定积分的计算中效果不佳。
其次,还有多元积分法,如变步长梯形法、双积分法等,这些积分法可以帮助求解一些复杂的定积分,但是计算时间较长。
此外,还有有限元法、隐式Runge-Kutta法、快速积分法等,这些积分方法能够帮助求解非定积分问题,其计算效率也相对较高。
数值积分方法在实际应用中得到了广泛的应用,如仿真求解有限元方程,求解复杂的拟合问题,估计系统的运行参数,计算力学分析等等都与数值积分技术有关。
另外,今天在这一领域,全球多家著名计算数值分析软件公司也在不断改进技术,开发出更加高效的数值积分软件,从而更好地服务于实际问题的求解。
总之,数值积分方法是一门重要的数值分析学科,可用于解决多种实际问题,广泛应用于科学和技术领域,具有重要的现实意义。
python数值积分Python是一种高级编程语言,广泛用于科学计算和数据分析。
在数学和科学计算领域,数值积分是一个重要的问题。
数值积分是指用数值方法计算函数的定积分,即给定一个函数$f(x)$和积分区间$[a,b]$,求$int_a^bf(x)dx$的近似值。
本文将介绍Python中常用的数值积分方法和库。
一、数值积分方法1.矩形法矩形法是最简单的数值积分方法之一。
它的思想是将积分区间$[a,b]$分成$n$个小区间,每个小区间的宽度为$h=frac{b-a}{n}$,然后用函数在小区间中点的函数值$f(frac{a+i*h}{2})$来近似表示小区间的面积,从而得到整个积分区间的近似值。
具体公式为:$$int_a^bf(x)dxapproxhsum_{i=0}^{n-1}f(frac{a+i*h}{2})$$矩形法的优点是简单易懂,容易实现。
但是它的精度较低,误差较大。
2.梯形法梯形法是一种比矩形法更精确的数值积分方法。
它的思想是将积分区间$[a,b]$分成$n$个小区间,每个小区间的宽度为$h=frac{b-a}{n}$,然后用小区间两端点的函数值$f(a+i*h)$和$f(a+(i+1)*h)$来近似表示小区间的面积,从而得到整个积分区间的近似值。
具体公式为:$$int_a^bf(x)dxapproxfrac{h}{2}sum_{i=0}^{n-1}[f(a+i*h)+f(a+(i+1)*h)]$$梯形法的优点是比矩形法更精确,误差较小。
但是它的计算量较大,对于复杂函数和大量数据,可能需要较长的计算时间。
3.辛普森法辛普森法是一种更加精确的数值积分方法。
它的思想是将积分区间$[a,b]$分成$n$个小区间,每个小区间的宽度为$h=frac{b-a}{n}$,然后用小区间两端点和中点的函数值$f(a+i*h)$,$f(a+(i+1)*h)$和$f(frac{a+i*h+a+(i+1)*h}{2})$来近似表示小区间的面积,从而得到整个积分区间的近似值。
1,荷载工况(load case):是对各种荷载类型的定义(define),然后通过指定(assign)建立模型中空间分布的力、位移或其他作用(例如:温度)。
这仅仅是建立了作用,荷载工况本身不在结构上产生响应。
2,分析工况(analysis case):是定义荷载作用方式(静力或动力)、结构的响应方式(线性或非线性)、分析方法(模态分析法或直接积分法)。
分析工况中包含荷载工况,分析工况可以对应一个荷载工况,也,可以是荷载的组合(多点风荷载、多维地震动)。
运行分析工况才能得到结构关于荷载的响应。
3,定义组合(define combination ):是将分析工况的计算结果进行组合(计算机运行减少人工进行计算的工作量),常用的组合形式是线性(linear)叠加或者包络(envelope)。
1.时程分析时用EL波,原始记录的波一般是以重力加速度g为单位,它的峰值为0.341g,也就是0.341*9.8m/s2.而你sap的单位用的是N/mm/s,也就是你的单位与原始波的单位相差1000*9.8个单位,那么你的系数要输入9800。
如果你sap的单位为N/m/s,那么你的系数取9.8即可。
2.规程中的8度罕遇要求是400g,这个g是单位gal的缩写字母,它的单位是cm/s2。
实际上就是0.4个重力加速度。
即400gal =0.4g,考虑第1点,那么你的系数应该取1000*9.8*(0.4/0.341)=11495.6。
3.定义时程函数时,单位无所谓,只要你的系数对应好就可以。
注:sap输入的地震函数本身是没有单位的,它的单位随着你sap 的右下角的单位走的。
所以才需要将这个单位和原始波单位对应。
1,将索得抗弯刚度设为极小值。
2,需作索的非线性分析,在作索得非线性分析需要打开大变形得选项。
3,加载需要分步加载,先加载预应力,再加载其它荷载。
4,在v9版本里面,可以直接用应变来直接模拟预应力,不用降温也可以。