MIKE21教程
不错的MIKE21中文教程，主要介绍MIKE21水动力模块方面的内容Www.Zlvo.Com42.7风场( Wind forcing)···;;;·36注意:42.8冰盖( ce coverage)4.2.9引潮势( Tidal potential)42.10降水-蒸发( Precipitation- Evaporation)….4142.1波浪辐射应力( Wave radiation)424212源( Sources4342.13水工结构物( Structures)454.2.14初始条件( nitial conditions)42.15边界条件( Boundary conditions)6142.16温度/盐度模块(Tcmpcraturc/Salinity Module)6742.17湍流模块( Turbulence module)42.8解耜( Decoupling)….…6742.9输出( Outputs)…特别说明:本手册部分内容来源于网络。Www.Zlvo.Com第一章模型介绍11简介MIKE21是一个专业的工程软件包,用于模拟河流、湖泊、河口、海湾、海岸及海洋的水流、波浪、泥沙及环境。MIKE21为工程应用、海岸管理及规划提供了完备、有效的设计环境。高级图形川户界面与高效的计算引擎的结合使得MIKE2I在世界范围内成为了一个水流模拟专业技术人员不可缺少的工具。丹麦水力研究所开发的平面二维数学模型MIKE21,曾经在丹麦、埃及、澳洲、泰国及中国香港、台湾等国家和地区得到成功应用,在丬面二维白由表面流数值模拟方面具有强大的功能。目前该软件在中国的应用发展很快,并在一些大型工程中广泛应用,如:长江口综合治理工程、杭州湾数值模拟、南水北调工程、重庆市城市排污评价、太湖富营养模型、香港新机场工程建设等。12MIKE21软件特点(1)用户界面友好,属于集成的 Windows图形界面;(2)具有强大的前、后攵理功能。在前处理方面,能根据地形瓷料进行计算网格的划分;在后处理方面具有强大的分析功能,如流场动态演示及动画制作、计算断面流量、实测与计算过程的验证、不同方案的比较等;(3)多种计算网格、模块及许可选择确俫用户根据自身需求来选择模型(4)可以进行热启动,当用户因各种原因需暂时中断MIKE21模型时,只要在上次计算时设置了热启动文件,再次开始计算时将热启动文件调入便可继续计算,极大地方便了计算时间有限制的用户;(5)能进行干、湿节点和干、湿单元的设置,能较方便地进行滩地水流的模拟:(6)具有功能强大的卡片设置功能,可以进行多种控制性结构的设置,如桥墩、堰、闸、涵洞等(7)可广泛地应用于二维水力学现象的研究,潮汐、水流,风暴潮,传热、盐流,水质,波浪紊动,湖震,防浪堤布置,船运,泥沙侵蚀、输移和沉积等,Www.Zlvo.Com被推荐为河流、湖泊、河∏和海岸水流的二维仿真模拟工具。1.3水动力模块原理131控制方程模型是基于三向不可压缩和 Reynolds值均布的 Navier-SLokes方程,并服从于 Boussinesq假定和静水压力的假定。二维非恒定浅水方程组为Ch Chu chvhSChu ahauvan h6x+=1-a=欧h-a(1-22pa ax po po ph)+-(h12)+hu,Schv chuy chvfuh-ghan h apay po aygh ap2 Po ay Po po po、ax11)+hS式中:t为时间:x,y为笛卡尔坐标系坐标;n为水位;d为静止水深;h=n+d为总水深;tn,v分别为x,y方向上的速度分量;f是哥氏力系数,f=2 osin p,(为地球白转角速度,为当地纬度;g为重力加速度;p为水的密度;Sx、SS分别为辐射应力分量;S为源项;(uy,ν)为源项水流流速。字母上带横杠的是平均值。例如,矿、ν为沿水深平均的流速,由以下公式定义hu= udz, hvdzWww.Zlvo.Com为水平粘滞应力项,包括粘性力、紊流应力和水平对流,这些量是根据沿水深平均的速度梯度用涡流粘性方程得出的:T=2A2A13,2数值解法)空间离散计算区域的空间离散是用有限体积法( Finite volume method),将该连续统体细分为不重叠的单元,单元可以是任意形状的多边形,但在这里只考虑三角形和四边形单元。在MKE软件2007版本只能是三角形网格。浅水方程组的通用形式一般可以写成上(U)=S(U)(1-6)式中:U为守恒型物理向量:F为通量向量;S为源项在笛卡尔坐标系中,二维浅水方程组可以写为OU O(F-F)O(FY-Fy)S(1-7)式中:上标/和分别为无粘性的和粘性的通量。各项分别如下:0hCu+g(FhuyOu Cha0Fk=lhAolhugh42Www.Zlvo.Comadh2thPu cy Pogn+fuhpe gn opythiPo oy Po对方程(46)第i个单元积分,并运用 Gauss原理重写可得出「a(Fa)-JA(1-9)式中:A1为单元g2的面积;I;为单元的边界;ds为沿着边界的积分变量这里使用单太求积法来计算面积的积分,该求积点位于单元的质点,同时使用中点求积法水计算边界积分,方程(49)可以写为∑FnAT=S(1-10)式中:U和S分别为第个单元的U和S的平均值,并位于单元中心;NS是单元的边界数;^厂,为第j个单元的长度阶解法和二阶解法都可以用于空间离散求解。对于二维的情况,近似的Riemann解法可以用来计算单元界面的对流流动。使用Roc方法时,界面左边的和右边的相关变量需要估计取值。二阶方法中,空间准確度可以通过使用线性梯度重构的技术来获得。而平均梯度可以用由 jawahar和 Kamath于2000年提出的方法来估计,为了避免数值振荡,模型使用了二阶TVD格式。(2)时间积分考虑方程的一般形式aU=G(U)1-11)对于二维模拟,浅水方程的求解有两种方法:一种是低阶方法,另一种是高阶方法。低价方法即低阶显式的Euer方法Un1=Un+△G(Un)(1-12)式中:为时间步长。高阶的方法为以如下形式的使用了二阶的 Runge kuttaWww.Zlvo.Com方法n12=Un+△G(U,)Un+1=Un+△G(Un+12)(1-13)(3)边界条件1)闭合边界沿着闭合边界(陆地边界),所有垂直于边界流动的变量必须为0。对于动量方程,可以得知沿着陆地边界是完全平稳的。2)开边界开边界条件可以指定为流量过程或者是水位过程3)千湿边界处理动边界问题(T湿边界)的方法是基于赵棣华(1994)和 Sleigh(1998)的处理方式。当深度较小时,该问题可以被重新表述,通过将动量通量设置为零以及只考恳质量通量来实现。只有当深度足够小时,计算不考虑该网格屮元。每个单元的水深会被监测,并且单元会被定义为干、半干湿和湿。单元面也会被监测,以确定淹没边界。满足下面两个条件单元边界被定义为淹没边界:首先单元的一边水深必须小于hn,且另一边水深必须大于h;第二,水深小于hn的单元的静水深加上另一单元表面高程水位必须大于零。满足下亩两个条件单元会被定义为干单元:首先单元中的水深必须小于干水深hn;另外,该单元的三个边界中没有一个是淹没边界。被定义为干的单元在计算中会被忽略不计。单儿破定义为半干:如果单元水深介于h和hm之间,或是当水深小于hy但有一个边界是淹没边界。此时动量通量被设定为0,只有质量通量会被计算。单元会被定义为湿:如果单元水深大于ha。此时动量通量和质量通量都会在计算中被考虑。如果模型中的区域是处在τ湿边交替区,为了避免模型计算岀现不稳定性,使用者可以启用 Flood and Dry选项。在这个情形下使用者必须设定一个干水深Www.Zlvo.Com( drying depth),淹没深度( flooding water depth)和湿水深( wetting depth)者必须满足hn
- 2020-12-09下载
- 积分：1
鲁棒控制——LMI处理方法
俞立老师的经典书籍——《鲁棒控制——基于线性矩阵不等式的处理方法》振动论坛鲁棒控制—线性矩阵不等式处理方法俞著清华大学出版社振动论坛(京)新登字158号内容简介木书结合作者的研究工作,详细介绍了基于线性矩阵不等式的不确定系统食棒控制的概念、理论及设计方法。主要为謇括日前应用广泛的线性矩阵不等式的概念、理论、算法及相关软件:基于线性矩阵不等式处理方法的线性时不变系统性能分析和缤合方法:重点介绍了不确定系统的模型、曹棒性能分妡、鲁棒砂以H控制、LM区域及相应的区域极点配置方法,绪合二次型性能指标的保性能控制、鲁榨方差控制对将系统的分析与鲁棒控制器设计、不确定系统的鲁滤波问题及鲁棒滤波器设计本书反映了近年来鲁棒控制领域中的最新研究成果,系统介纲了线性矩阵不等式这一有效工具,它在应用中的典型处理方法及 MATLAB软件中的LMI具箱,本书可作为从事自动控制工作的科研员、工技术人员以及高等院校自动化及其他相关专业教师高年级学生和研究生的参老用书版权所有,翻印必究本书封面貼有清华大学出版杜激光防伪标签,无标签者不得销售。书名:鲁棒控制—线性矩阵不等式处理方法作者:俞立著出版者:清华大学出版社(北京清华大学学研大厦,的编100084)http://www.lup,lsingt:ua.educn费任编辑:朱英彪印刷者:北京通州区大中印刷发行者:新华书店总店北京发行所开本;787×10921:16即张;17.75字数:403千字版次:202年12月第1版20年12月第1次印刷书号:ISEN730245854-7/O·269印数:0001~400定价:26.00元振动论坛前言在实际工业控制牛,各种工业牛产过程、生产设备以及其他众多的被控对象,其动态牿性一般鄱难以用精确的数学模型来描述。有时即使能获得被控对象的精确数学模型,但白于过于短杂,使得难以对其进行有效的掉制性能分析和综合,因此必须进行适当的简化。另一方面,难着生产过程中主作条件和环境的变化,控制系统中元器件的老化或损坏,被控对象本身的特性也会随之发生变化。所有这些因素使得描述被控对象的数学模型和实际对象之间不可避免地具有误差。因此,在工程实践中,釆用基于精确数学模型的现代控制理论方法所设计的控制系统往往难以具有所期望的性能,甚至连系统的稳定性都难以得到保证。鲁棒控制理论结合系统模型参数不确定性和外部找动不确定性的考虑,研究系统的鲁棒性能分析和综合问题,弥衤了现代控制理论需要对象精确数学模型的缺陷,使得系统的分析和综合方法加有效、实用棒挖制自提出以来,很快受到了人们的广泛重视和研究,取得了一系列的研究结果和方法,并在一些工程领域中获得了成功的应用。特别地,随着线性短阵不等式及求解凸优化问题的内点法的提出,为许多控制问题的分析和求解提供了有效工具。 MATLAB软件中线性矩阵不等式工具箱的推出使得各种线性矩阵不等式问题的求解旻加方便、直接从而,进一步推动了线性矩阵不等式处理方法在系统和控制领域中的应用。本书系统介绍了线性矩阵不等式的概念、性质、求解线性矩阵不等式相关问题的算法以及 MATLAB软件中的线性矩阵不等式工具箱。结合作者的研究工作,介绍了基丁线性矩阵不等式的鲁棒控制性能分析和综合方法,并指出了一些在系统和控制中有广泛应用价值的线性矩阵不等式典型处理方法。作者努力将售棒控制的最新研究成果和方法反映存本书中,但限于篇幅,书中所包含的内容仅仅是鲁棒控制研究成果的很少-部分。振于作者的水平,书中不妥和错误之处在所难免,恳请广大读者批评指正本书中介绍的作者研究工作及本书的撰写得到了国家然科学基金和教育部高校优秀青年教帅教学科研奖励计划的资助,在此表示衷心的感谢。作者02年5月杭州振动论坛目录第1章引言…第2章线性矩阵不等式……621线性矩阵不等式的表示式口·■■■甲↓郾■■.■看■■bL■2Ll1线性矩阵不等式的般表示2l2可转化成线性矩阵不等式表小的问题213复线性矩阵不等式的处理…10214非严格线性矩阵不等式22一堂标准的线性矩阵不等式间题23求解线性矩阵不等式问题的算法23.1椭球法…Isn15232内点法↓乐24关于矩阵不等式的一些结论1824.1矩阵变量的消去法1T·平P1t幽山I82.4.2 S-procedureP甲■20第3章系统性能分析.连续时间系纷…2331.1系统增益指标∴233,12H性能甲日唱·■唱日4自·日血即141■4b1293.13H性能3132离散时间系统1鲁日日司即■命b命即■■■■■着■■■着L=,341第4章控制系统综合.4.1F控制4.l.1状态反馈H控制…d?■■424.],2输出反馈Hx控制■p即42H2控制43H2H控制平1■十■昏■b■m■■■4■qp6144设计示例64第5章不确定系统的分析与综合685]不确定模型68振动论坛會棒控制——一线性矩阵不等式处理方法511不确定状态空闻模型68512不确定线性分式模型■p■b血d卩■b■52鲁棒稳定性分忻P‘■■日音牛114a日4521二次稳定性522仿射二次稳定性I"P}4Pbb4b4}·53鲁棒性能分析.…■■昏■昏1山■d■晶t■■p聊8354鲁棒F2H控制■■■山昏晉昏晋■+口■口■口6541问题的描述和准备…86542H2H整制器设计91第6章区域极点配置1亠山品品61LMI区域9761.】LM区域的描述…,-11111976,2D-稳定性分析…〓P■■■■■■■■■■10062具有闭坏k域极点约束状态反馈控制器设计0463鲁棒D-稳定性分析日日日语■山■,1763无结构不确定性1963.2结构不确定性111464输出反馈控制器设计,第7章保性能控制1227】连续系统的保性能控制…··■■晶■上h香≠p甲日目甲日唱黑L聊■12272离散系统的保性能掉制.12了73其有闭环极点约束的保性能挡制131731鲁棒性能分析…132732二次D保性能控制器设计看自;血自自血即血日昏■是■如4135第B章鲁棒方差控制1418I连续系统的鲁棒方差控制.■d昏【山■血晶b14181.Ⅰ系统性能分析1418L.2状态反馈控制器设计14813输出反愤控制器设计14682离散系统的鲁棒方差控制152第9章时滞系统的分析和综合15891时滞系统的稳定性…15891.1时滞独立的稳定性条件……,1599.1.2时滞依赖的稳定性条件,+--114141609.1.3 Lurie时滞系统的稳定性分析T宁宁·自血日甲自自要振动论坛自录92时滞系统的鲁棒稳定性分析16992.!时滞独立的鲁棒稳定性条件P·q甲,甲■15992.2时滞依赖的鲁棒稳定性条件17493不确定滞系统的保性能控制.■■■■793.1鲁棒性能分析I789.32状态反馈俣性能控制器设计,4.4.418393.3输出反馈俣性能控制器设计1869.34不确定离散时滞系统的保性能控制I9394时滞系统的玨控制941时滞系统的性能分析,…199942H控制器设计943不确定离散时滞系统的鲁棒H控制207第10章寯波器设计213101B滤波器设计2I310.2FH滤波器设计2l9第11章大系统射分散控制22311时滞系统的分散稳定化控制22311,2离散关联系统的分散保性能控制229l1,21保性能分析■■■山■↓山 -r.E= PPE+■▲画2291121分散保性能控制器设计234附录ALM工具箱介绍■L■■■an·■■甲口■看■■口■山241A.F线性矩阵不等式及相关术语……24LA.2线性矩阵不等式的确定■■■■昏昏斷■+d山夏旷aaa■_画d■242A.3信息提取2499A4绒性矩阵不等式求解器■■■■■■日■司罾■严早■■■吾-■山■山m矿画maaa250A.5结果验证…258A6修改一个线性矩阵不等式系统A.7一些进一步的功能261A.8系统棋型描述…,267参考文献270振动论坛第1章引言自20世纪5年代末现代控制理论誕生以米,控制理论得到了飞速的发展,并在20世纪60年代然航天领域中得到成功的应用。但是,现代控制理论在随后的工业应用中却遇到了很大的困难。我们知道,现代控制理论的许多结果都是基于对象的一个数学模型,根据系统的性能要求,通过刈被控对象的数学模型进行分析来设计系统的控制律,进而将所得到的控制律应用于被控对象来保证闭环系统貝有所期望的性能。显然,当对象模型不能精地描述被控对象或在系统运行过程中模型和实际对象产生侃离时,基于这栟的模型设计的控制系统很难保证具有所期望的性能要求实际上,对于复杂物理系统的模型,存在以下两个问题描述物理系统的解析模型很难,甚至不可能精确地刻画,因此为了便于处理,不得不筒化模刑:2.一个模型,无论多么详细,都不可能是物理系统的一个濤确表示。因此,模型存在本质的不精确性。建模中的以上两个方面称为模型的不确定性。对于一个复杂系统,为了得到一个较为简单的模型,一种处珄方法是将其分解线性部分和非线性部分的组合,进而用一个更空易处理和分析的对象来替代这个非线性邺分,达到简化原文复杂系统棋型的目的考虑由以下非线性微分方程描述的复杂动态系统x=f(x, m)H(x, u)初始条件是x(0),x()、y()和()是向量值函数,∫和h是光滑的向量值函数。在特殊的运行点附近,可以将系统(1.)分解成一个线性部分和一个线性部分的组合特别地,可以在原点x,)=(0,0处进行这样的分解。定义系统Ax t Bu+gtr, Wv=Cx+IH+rx, uj其中:A、B、C和D是系统(11)的一个线性化近似g(x,叫)=f(x,)-Ax-Br(, u)-hr, u)-Cx-Du显然,这样定义的系统(1.2)和系统(1:)是等价的。因此,它们之间存在-—对应的关系。得到这样的等价系统的一种方式是将函数∫和h在原点处线性化,可得Cl(x.)=(0.0)x,)=(0,Q)0.0Ca(x,)=()振动论坛鲁棒控制——线性矩阵不等式处理方浊进一步可以将方程(12)受成以下等价的形式:x= Ar I Bu I wL1.3y=cx+Da+脚2(14)(w1,w2}=(g(x,L),r(x,)(1.5设G是宙(1.3)~(1.4)式确定的映射;对给定的初始条作x(0),門2,)以(x环,y)。Ω是由(1.5)式确定的映射:(x,n)卜(1,2)。因此,〈1.3(1〕式描述的系统可以用图11来表示图11系统分解容易看到,G是系统的线性部分,Q是静杰的非线性映射。这样就将系统的非线性部分分离出来,归入到映射g中,非线性部分和线性部分通过反馈关联联系起来更一般地,我们用这样的方法不仅可以处理系统的非线性特性,而且也可以处理系统的某些动态特性。考虑由以下方程组描述的系练f2(x1,x2(1y=H(x1,x2,采用前面系统分解思想,对系统(16)中的方程x=f(x1,x2,4)卩=(x,x2进行分解,并得到:=团x+B+8(x2f2(x1,x2,a)(1.)卩=C1x+D4+r(x,x2,a)进步,系统(17)中的方程等价于以下的线性方程:x=A1x1+BM+即8y=C,+Du+M其中,w2}=(g1(x,x2,),P(x,x2,Bm)x2=f1(x,x2,x)(1,9)设G是由方程(18)描述的线性系统;(w3,2,(x,,y),是由(9)式描述的系统:(x1,4)"(m,m2)。对这样定义的G和Q,图1也同样描述了系统(1.6)。
- 2020-12-06下载
- 积分：1
