CN112464391A - 一种柴油发电机组哈密顿模型的建立方法 - Google Patents
一种柴油发电机组哈密顿模型的建立方法 Download PDFInfo
- Publication number
- CN112464391A CN112464391A CN202011106209.8A CN202011106209A CN112464391A CN 112464391 A CN112464391 A CN 112464391A CN 202011106209 A CN202011106209 A CN 202011106209A CN 112464391 A CN112464391 A CN 112464391A
- Authority
- CN
- China
- Prior art keywords
- generator
- actuator
- diesel engine
- function
- hamiltonian
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims description 41
- 238000013016 damping Methods 0.000 claims description 20
- 230000001052 transient effect Effects 0.000 claims description 18
- 230000033001 locomotion Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000005381 potential energy Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008278 dynamic mechanism Effects 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开一种柴油发电机组哈密顿模型的建立方法,属于发电机组稳定分析与控制技术领域,柴油发电机组由柴油机和发电机组成,是一个复杂的非线性系统,本发明从动力学层面研究其建模方法,为进一步对其进行控制奠定基础,首先利用分析动力学原理将柴油机拉格朗日函数转换为哈密顿函数,并将其加入到发电机哈密顿函数中,获得柴油发电机组整体的哈密顿函数,基于所获得的哈密顿函数和柴油机本体及执行器的微分方程,推导柴油机哈密顿方程,并和发电机的哈密顿方程一起构成整个发电机组的哈密顿模型。
Description
技术领域
本发明属于柴油发电机技术领域,具体是柴油发电机组哈密顿建模技术,更为具体的是涉 及一种柴油机发电机组哈密顿模型的建立方法。
背景技术
柴油发电机组是一种灵活、快速的供电电源,被广泛应用于各种用电场所,尤其是作为一 种备用电源,被各行各业采用,当前随着大规模新能源的开发与利用,特别是微电网的快速发 展,利用柴油发电机组作为调节电源,平抑微电网功率波动是一种有效的方法,使得柴油发电 机组的应用得到快速发展。在此背景下,柴油发电机组与其它设备的协同控制技术成为研究的 热点,要求对柴油发电机组运行控制的内部机理进行更深入的探索。
柴油发电机组是由柴油机和发电机两个子系统组成的非线性系统,由于发电机建模问题比 较成熟,研究重点主要集中于柴油机系统。目前还没有关于柴油发电机组哈密顿建模的方法。
发明内容
本发明将柴油机拉格朗日系统演变得到广义哈密顿系统,进而连接发电机构成柴油发电机 组的广义哈密顿模型,提出柴油发电机组的哈密顿建模方法,为基于系统结构动力学的柴油发 电机组运行控制研究提供基础。
本发明提出一种柴油发电机组的哈密顿建模方法,柴油发电机组包含柴油机和发电机,本 发明在柴油机拉格朗日方程的基础上,采用动力学方法将拉格朗日函数转换为哈密顿函数,获 得柴油机的哈密顿函数,并将其加入到发电机哈密顿函数中,获得柴油发电机组整体的哈密顿 函数并结合微分方程,推导柴油发电机组哈密顿方程,具体建模方法如下。
一种柴油发电机组哈密顿模型的建立方法,柴油发电机组包括柴油机和发电机,包含以下 几个步骤:
(1)将柴油机拉格朗日函数转换为哈密顿函数:
根据动力学理论,定义如下的哈密顿函数H(1):
其中:H(1)为柴油机哈密顿函数;为执行器运动部分质量标么值;为执行器位移标么值; 为执行器移动速度标么值 是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角 速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,为阻尼系数标么值; TD、Tb、;xnl为空载位移;为执行器的输入;m2为发电机电磁力矩;D为发电机等 效阻尼系数;以上参数均为标么值;以上的各个参数通过设备额定参数、实验直接获得或经过 已知参数计算获得;
(2)建立柴油发电机组的哈密顿能量函数:
发电机单机无穷大系统哈密顿模型采用文献“Zeng Yun etc.,The generalizedHamiltonian model for the shafting transient analysis of the hydro turbinegenerating sets[J].Nonlinear Dynamics, 2014,76(4),1921-1933.”中给出的模型如下。
发电机哈密顿能量函数:
哈密顿方程为:
其中,δ为发电机功角(rad),Eq′为q轴瞬变电动势相对值,Us为无穷大系统电压相对值, Xd∑=Xd+XT+XL,Xd为发电机d轴电抗,XT为变压器等效电抗,XL为线路等效电抗, Xq∑=Xq+XT+XL,Xq为发电机q轴电抗,Xd∑′=Xd′+XT+XL,Xd′为发电机d轴次暂态电抗,Xf为励 磁绕组电抗,Xad为d轴电枢反应电抗,各电抗均为相对值,Td0′为d轴开路暂态时间常数(秒), mt为水轮机力矩,Ef为励磁控制输出,以上参数中没有标明单位均为标么值;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得;
(3)建立柴油机哈密顿方程:
将步骤(2)中获得的柴油机哈密顿方程与已知的发电机哈密顿方程联合得到柴油发电机 组哈密顿方程:
所述柴油机拉格朗日拉格朗日函数的建立方法,具体步骤如下:
(1)执行器的拉格朗日函数的建立:
设柴油机电磁执行器的输入信号为u,输出轴位移为x,衔铁运动在轴上产生的电磁力F 为二者的函数,表示为F(x,u),假设初始执行器工作在输入信号为u0,位移为x0的稳定状态, 则电磁力表示为F(x0,u0),当输入信号变化Δu,相应轴位移变化Δx,电磁力变化为:
F(x0,u0)+ΔF=F(x0+Δx,u0+Δu)
将上式右边采用泰勒级数展开得到:
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能, 即:
其中:v是轴运动速度,m/s;m1是质量,kg;k1是弹簧刚度,N/m;x是执行器位移,mm;
(2)柴油机旋转部分轴系统的拉格朗日函数的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转 机械角位移θm(rad)为广义坐标,轴系统的拉格朗日函数L2等于轴的旋转动能:
其中:J为转动部分的转动惯量,kg.m2;ωm为机械角速度为,rad/s;
(3)柴油机拉格朗日函数的建立:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将步骤1和步骤2中所获得的 执行器和本体的能量函数L1和L2相加,可得到柴油机的拉格朗日函数:
其中:L为拉格朗日函数;m1为执行器运动部分质量,kg;k1是执行器弹簧刚度,N/m;x1为执行器位移,mm;x2为执行器速度,m/s;J为转动部分的转动惯量,kg.m2;ωm为转动 部分的机械角速度,rad/s;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数 计算获得。
本发明的有益效果是:
1、柴油发电机组哈密顿模型与微分方程模型比,给出了结构和阻尼矩阵,揭示了柴油发电 机组运行控制的内部动力学机理。
2、基于该发明,可以实现复杂的非线性控制,设计出基于动力学内部耦合机理的非线性控 制器。
3、该模型可以用于涉及柴油发电机组稳定分析和仿真的商业软件系统中。
4、本发明提供了柴油发电机组哈密顿建模的方法和哈密顿模型,所获得的哈密顿模型提供 了系统的结构矩阵和阻尼矩阵,它揭示了系统内部关联耦合的机理,为进一步基于哈密顿模型 的控制提供了方便。本发明为柴油发电机组的建模提供了一种新方法。
附图说明
图1为柴油发电机组与电网的连接图;
图2为柴油发电机哈密顿函数随时间的变化曲线;
图3为柴油发电机哈密顿函数与有功功率的变化曲线;
图中,1-柴油机,2-发电机,3-变压器,4-输电线路,5-电网。
具体实施方案
为了进一步表述本发明的内容,结合具体例子进行进一步分析。
实施例1
一种柴油发电机组哈密顿模型的建立方法,柴油发电机组包括柴油机和发电机,包含以下 几个步骤:
步骤1:建立柴油机拉格朗日函数:
(1)执行器拉格朗日模型的建立过程:
执行器是柴油机对外部控制信号的执行机构,通过它实现对柴油机的控制,各种类型的柴 油机执行器抽象为数学模型后都是相同的。
设柴油机电磁执行器的输入信号为u,输出轴位移为x,衔铁运动在轴上产生的电磁力F 可表示为F(x,u),假设初始执行器工作在输入信号为u0,位移为x0的稳定状态,则电磁力可表 示为F(x0,u0),当输入信号变化Δu,相应轴位移变化Δx,电磁力变化为:
F(x0,u0)+ΔF=F(x0+Δx,u0+Δu)
将上式右边采用泰勒级数展开得到:
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能, 即:
其中:v是轴运动速度(m/s),m1是质量,k1是弹簧刚度(N/m),x是执行器位移(mm);
执行器在运动过程中受到阻尼,设阻尼系数为c1(N.s/m),阻尼与速度成正比,即大小为c1v, 所以轴上的外力Q为电磁力和阻尼力之差:Q=F(x,u)-c1v
在平衡状态下,输入控制信号和执行器轴位移恒定,即u=u0,x=x0,增量为零:Δu=0, Δx=0,代入上式方程有k1x0=F(x0,u0),即弹簧力与电磁力相等,方程可表示为:
其中,kx1=kx-k1,Δx=x-x0,Δu=u-u0;
在平衡点处上式仍然要满足,所以将平衡点坐标(u=u0,x=x0)代入上式得:kx1x0=-kuu0, 于是上式可以改写为如下形式:
(i)式即是执行器的拉格朗日方程;
(2)柴油机旋转部分拉格朗日模型的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转 机械角位移θm(rad)为广义坐标,系统的拉格朗日函数L2等于轴的旋转动能:
其中:J为转动部分的转动惯量(kg.m2),ωm为机械角速度为(rad/s);
设柴油机轴力矩为M1,若柴油机输出端连接发电机,设发电机电磁力矩为M2,它代表的 是在柴油机轴上产生的制动性力矩,再计及发电机阻尼绕组产生的阻尼力矩Md,则作用在柴 油机轴上的非保守广义外力(力矩)为:
Q2=M1-M2-Md
其中,柴油机轴上的力矩M1与轴位移和角速度为线性关系,可表示为:
M1=kω(ωm-ωmB)-a1(x-xnl)
其中,kω、a1为系数,a1=M1r/(xr-xnl),M1r是柴油机额定转矩,xr是额定扭矩对应的轴位移,xnl是空载位移,即a1是柴油机力矩随轴位移变化的增益系数或斜率;ωmB是额定机械角速度;
阻尼转矩与角速度也成正比:Md=Dp(ωm-ωmB),其中,D为柴油机转动阻尼系数;p为发电机磁极对数;所以Q2=kω(ωm-ωmB)-a1(x-xnl)-M2-Dp(ωm-ωmB)
其中r0=kω-Dp;r1=-r0ωmB+a1xn1;ωmB为额定机械角速度,M2为发电机电磁力矩;
(3)求取柴油机拉格朗日函数和拉格朗日方程:
(a)柴油机拉格朗日函数:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将所获得的执行器和本体的能 量函数L1和L2相加,可得到柴油机的拉格朗日函数:
(b)柴油机拉格朗日方程:
根据上述推导,将柴油机拉格朗日方程整理为标准的形式如下,设x1为执行器位移,x2为 执行器运动速度v,即x2=dx1/dt,代入(i)式得:
式(iii)即为柴油机的拉格朗日方程;
其中:L为拉格朗日函数;m1为执行器运动部分质量,kg;k1是执行器弹簧刚度,N/m;x1为执行器位移,mm;x2为执行器速度,m/s;J为转动部分的转动惯量,kg.m2;ωm为转动 部分的机械角速度,rad/s;u为执行器的输入信号;c1为执行器的阻尼系数;kx1、ku、kω、 a1为系数;r0=kω-Dp;r1=-r0ωmB+a1xn1;D为柴油机转动阻尼系数;p为发电机磁极对数; ωmB为柴油机额定机械角速度;xnl为柴油机额定转速下空载位移;M2为柴油机负荷力矩;以 上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得;
步骤2:将柴油机拉格朗日函数转换为哈密顿函数:
根据动力学理论,定义如下的哈密顿函数H(1):
其中:H(1)为柴油机哈密顿函数;为执行器运动部分质量标么值;为执行器位移标么值;为执行器移动速度标么值;是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,为阻尼系数 标么值;TD、Tb、;xnl为空载位移;为执行器的输入;m2为发电机电磁力矩;D为发电机等效阻尼系数;以上参数均为标么值;以上的各个参数通过设备额定参数、实验直接获 得或经过已知参数计算获得;
步骤3:建立柴油发电机组的哈密顿能量函数:
发电机单机无穷大系统哈密顿模型采用文献“Zeng Yun etc.,The generalizedHamiltonian model for the shafting transient analysis of the hydro turbinegenerating sets[J].Nonlinear Dynamics, 2014,76(4),1921-1933.”中给出的模型如下:
发电机哈密顿能量函数:
哈密顿方程为:
其中,δ为发电机功角(rad),Eq′为q轴瞬变电动势相对值,Us为无穷大系统电压相对值, Xd∑=Xd+XT+XL,Xd为发电机d轴电抗,XT为变压器等效电抗,XL为线路等效电抗, Xq∑=Xq+XT+XL,Xq为发电机q轴电抗,Xd∑′=Xd′+XT+XL,Xd′为发电机d轴次暂态电抗,Xf为 励磁绕组电抗,Xad为d轴电枢反应电抗,各电抗均为相对值,Td0′为d轴开路暂态时间常数(秒), mt为水轮机力矩,Ef为励磁控制输出,以上参数中没有标明单位均为标么值;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得;
步骤4:建立柴油发电机组哈密顿函数及方程:
柴油发电机组与电网的连接如图1所示,包括柴油机1、发电机2、变压器3、输电线路4、 电网5;柴油机1与发电机2连接,组成柴油发电机组,发电机2与变压器3连接,变压器3通过输电线路与电网5连接,将图1中的,将柴油机与发电机集合为一个系统,重新选取状态变量:v3=δ,v4=ω1,v5=Eq′,两个子系统哈密顿函数合并,其中,由于发电机与柴油机是同轴运行,在将二者的能量函数合成时,分项:只能取一次,合并后的函数为:
将步骤(2)中获得的柴油机哈密顿方程(iv)与已知的发电机哈密顿方程(v)联立得到柴油发电 机组哈密顿方程:
其中,H为柴油机哈密顿函数;为执行器运动部分质量标么值;为执行器位移标么 值;为执行器移动速度标么值;是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角 速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,为阻尼系数标么值; TD、Tb、;xnl为空载位移;为执行器的输入标么值;m2为发电机电磁力矩;D为发 电机等效阻尼系数;以上参数均为标么值;δ为发电机功角(rad),Eq′为q轴瞬变电动势相对值, Xd∑=Xd+XT+XL,Xd为发电机d轴电抗,XT为变压器等效电抗,XL为线路等效电抗, Xq∑=Xq+XT+XL,Xq为发电机q轴电抗,Xd∑′=Xd′+XT+XL,Xd′为发电机d轴次暂态电抗,Xf为励 磁绕组电抗,Xad为d轴电枢反应电抗,各电抗均为相对值,Td0′为d轴开路暂态时间常数(秒), mt为水轮机力矩,Ef为励磁控制输出,以上参数中没有标明单位均为标么值。
本实施例完成了柴油发电机组哈密顿函数的建立。
实施例2
为了验证实施例1哈密顿函数的正确性,从两个方面进行仿真,(1)哈密顿函数随时间的 变化;(2)哈密顿函数与发电机有功功率的变化。
仿真参数:
某柴油发电机组主要参数:柴油机额定功率1250kW,额定转速n=1500r/min,m1=0.8kg, 转动惯量J=71.822kg.m2,发电机磁极对数p=2;发电机阻尼系数D=0.1753;直轴电枢电抗Xad=1.96,发电机时间常数Td0=2.8;直轴同步电抗Xd=2.053;交流同步电抗Xq=1.0;直轴暂态电 抗X′d=0.213;交轴暂态电抗X′q=0.168,励磁绕组电抗Xf=2.09;线路电抗XL=0.1;变压器等 效电抗XT=0.15,无穷大母线电压US=1.0。
仿真工况一:在0秒时给定0.5(pu)的输入阶跃扰动,哈密顿能量函数的变化曲线,如 图2,图2表明暂态过程中哈密顿函数值从初始稳态值开始出现震荡,经历暂态过程后进入新 的平衡点;暂态过程中哈密顿函数保持为正值,并且随着输入的稳定而趋近于新的稳态值,其 变化趋势较好地反映了暂态过程中系统能量的变化。也说明所建立的柴油发电机哈密顿函数是 合理的。
仿真工况二:发电机输出功率从空载到额定值,哈密顿函数的变化曲线,如图3所示,图 3表明,稳态工况下,在柴油发电机组功率范围内哈密顿函数的变化与有功功率的变化一致, 且系统的哈密顿函数均大于零,H为正定函数,说明所建立的系统的哈密顿函数是合理的。
图2和图3采用两种工况,从暂态和稳态两个方面分析了哈密顿函数的变化,仿真结果说 明所建立的柴油发电机组的哈密顿函数是合理的。
Claims (2)
1.一种柴油发电机组哈密顿模型的建立方法,柴油发电机组包括柴油机和发电机,其特征在于,包含以下几个步骤:
(1)将柴油机拉格朗日函数转换为哈密顿函数:
根据动力学理论,定义如下的哈密顿函数H(1):
其中:H(1)为柴油机哈密顿函数;为执行器运动部分质量标么值;为执行器位移标么值;为执行器移动速度标么值;是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,为阻尼系数; TD、Tb、;xnl为空载位移;为执行器的输入标么值;m2为发电机电磁力矩;D为发电机等效阻尼系数;以上参数均为标么值;
(2)建立发电机的哈密顿能量函数:
发电机哈密顿能量函数为:
哈密顿方程为:
其中,δ为发电机功角(rad),Eq′为q轴瞬变电动势相对值,Us为无穷大系统电压相对值,Xd∑=Xd+XT+XL,Xd为发电机d轴电抗,XT为变压器等效电抗,XL为线路等效电抗,Xq∑=Xq+XT+XL,Xq为发电机q轴电抗,Xd∑′=Xd′+XT+XL,Xd′为发电机d轴次暂态电抗,Xf为励磁绕组电抗,Xad为d轴电枢反应电抗,各电抗均为相对值,Td0′为d轴开路暂态时间常数(秒),mt为水轮机力矩,Ef为励磁控制输出,以上参数中没有标明单位均为标么值;
(3)建立柴油发电机组哈密顿方程:
将步骤(2)中获得的柴油机哈密顿方程与已知的发电机哈密顿方程联合得到柴油发电机组哈密顿方程:
2.根据权利要求1所述柴油发电机组哈密顿模型的建立方法,其特征在于,所述柴油机拉格朗日拉格朗日函数的建立方法,具体步骤如下:
(1)执行器的拉格朗日函数的建立:
设柴油机电磁执行器的输入信号为u,输出轴位移为x,衔铁运动在轴上产生的电磁力F为二者的函数,表示为F(x,u),假设初始执行器工作在输入信号为u0,位移为x0的稳定状态,则电磁力表示为F(x0,u0),当输入信号变化Δu,相应轴位移变化Δx,电磁力变化为:
F(x0,u0)+ΔF=F(x0+Δx,u0+Δu)
将上式右边采用泰勒级数展开得到:
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能,即:
其中:v是轴运动速度,m/s;m1是质量,kg;k1是弹簧刚度,N/m;x是执行器位移,mm;
(2)柴油机旋转部分轴系统的拉格朗日函数的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转机械角位移θm(rad)为广义坐标,轴系统的拉格朗日函数L2等于轴的旋转动能:
其中:J为转动部分的转动惯量,kg.m2;ωm为机械角速度为,rad/s;
(3)柴油机拉格朗日函数的建立:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将步骤1和步骤2中所获得的执行器和本体的能量函数L1和L2相加,可得到柴油机的拉格朗日函数:
其中:L为拉格朗日函数;m1为执行器运动部分质量,kg;k1是执行器弹簧刚度,N/m;x1为执行器位移,mm;x2为执行器速度,m/s;J为转动部分的转动惯量,kg.m2;ωm为转动部分的机械角速度,rad/s;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011106209.8A CN112464391B (zh) | 2020-10-15 | 2020-10-15 | 一种柴油发电机组哈密顿模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011106209.8A CN112464391B (zh) | 2020-10-15 | 2020-10-15 | 一种柴油发电机组哈密顿模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464391A true CN112464391A (zh) | 2021-03-09 |
CN112464391B CN112464391B (zh) | 2022-09-06 |
Family
ID=74834039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011106209.8A Active CN112464391B (zh) | 2020-10-15 | 2020-10-15 | 一种柴油发电机组哈密顿模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464391B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113110023A (zh) * | 2021-04-02 | 2021-07-13 | 昆明理工大学 | 一种基于柴油机哈密顿模型的结构修正控制方法 |
CN113987988A (zh) * | 2021-10-11 | 2022-01-28 | 昆明理工大学 | 一种逆变器虚拟发电机哈密顿模型的建立方法 |
CN114004063A (zh) * | 2021-10-11 | 2022-02-01 | 昆明理工大学 | 一种异步电机哈密顿模型的建立方法 |
CN115296330A (zh) * | 2022-06-17 | 2022-11-04 | 昆明理工大学 | 一种孤网模式下的哈密顿阻尼协同控制方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813053A (zh) * | 2010-03-10 | 2010-08-25 | 昆明理工大学 | 改善水轮发电机组功角振荡的互联结构修正控制方法 |
CN101915203A (zh) * | 2010-08-06 | 2010-12-15 | 昆明理工大学 | 改善水轮发电机组功角振荡的阻尼注入控制方法 |
US8121708B1 (en) * | 2007-03-21 | 2012-02-21 | Sandia Corporation | Control system design method |
US20120221301A1 (en) * | 2011-02-28 | 2012-08-30 | Fujitsu Limited | Matrix generation technique and plant control technique |
CN105162377A (zh) * | 2015-08-19 | 2015-12-16 | 昆明理工大学 | 一种建立发电机励磁控制器能量函数的方法 |
CN106547957A (zh) * | 2016-10-19 | 2017-03-29 | 天津大学 | 一种旋转环状周期结构参激弹性振动分析方法 |
CN107147133A (zh) * | 2017-05-24 | 2017-09-08 | 国网辽宁省电力有限公司 | 一种提升新能源消纳的电网分时多频控制方法 |
CN109347141A (zh) * | 2018-10-19 | 2019-02-15 | 华北水利水电大学 | 一种双馈风力发电系统网侧终端滑模控制器设计方法 |
US10295973B1 (en) * | 2016-07-12 | 2019-05-21 | Lockheed Martin Corporation | System and method for control and guidance of an object |
US20190164079A1 (en) * | 2017-11-28 | 2019-05-30 | International Business Machines Corporation | Cost function deformation in quantum approximate optimization |
-
2020
- 2020-10-15 CN CN202011106209.8A patent/CN112464391B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8121708B1 (en) * | 2007-03-21 | 2012-02-21 | Sandia Corporation | Control system design method |
CN101813053A (zh) * | 2010-03-10 | 2010-08-25 | 昆明理工大学 | 改善水轮发电机组功角振荡的互联结构修正控制方法 |
CN101915203A (zh) * | 2010-08-06 | 2010-12-15 | 昆明理工大学 | 改善水轮发电机组功角振荡的阻尼注入控制方法 |
US20120221301A1 (en) * | 2011-02-28 | 2012-08-30 | Fujitsu Limited | Matrix generation technique and plant control technique |
CN105162377A (zh) * | 2015-08-19 | 2015-12-16 | 昆明理工大学 | 一种建立发电机励磁控制器能量函数的方法 |
US10295973B1 (en) * | 2016-07-12 | 2019-05-21 | Lockheed Martin Corporation | System and method for control and guidance of an object |
CN106547957A (zh) * | 2016-10-19 | 2017-03-29 | 天津大学 | 一种旋转环状周期结构参激弹性振动分析方法 |
CN107147133A (zh) * | 2017-05-24 | 2017-09-08 | 国网辽宁省电力有限公司 | 一种提升新能源消纳的电网分时多频控制方法 |
US20190164079A1 (en) * | 2017-11-28 | 2019-05-30 | International Business Machines Corporation | Cost function deformation in quantum approximate optimization |
CN109347141A (zh) * | 2018-10-19 | 2019-02-15 | 华北水利水电大学 | 一种双馈风力发电系统网侧终端滑模控制器设计方法 |
Non-Patent Citations (4)
Title |
---|
ZENG YUN: "Hamiltonian control strategy with modified interconnection structure for a hydroelectric generating set", 《CONTROL THEORY& APPLICATIONS》, 31 December 2009 (2009-12-31), pages 795 - 799 * |
张奕黄: "由哈密顿方程建立电机模型", 《北方交通大学学报》, 28 February 1997 (1997-02-28), pages 95 - 97 * |
曾云: "电站局部多机条件下五阶发电机哈密顿模型", 《中国电机工程学报》, 25 January 2014 (2014-01-25), pages 415 - 422 * |
曾云等: "发电机单机无穷大系统动力学模型的理论研究", 《中国电机工程学报》, no. 17, 15 June 2008 (2008-06-15), pages 140 - 145 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113110023A (zh) * | 2021-04-02 | 2021-07-13 | 昆明理工大学 | 一种基于柴油机哈密顿模型的结构修正控制方法 |
CN113110023B (zh) * | 2021-04-02 | 2024-02-27 | 昆明理工大学 | 一种基于柴油机哈密顿模型的结构修正控制方法 |
CN113987988A (zh) * | 2021-10-11 | 2022-01-28 | 昆明理工大学 | 一种逆变器虚拟发电机哈密顿模型的建立方法 |
CN114004063A (zh) * | 2021-10-11 | 2022-02-01 | 昆明理工大学 | 一种异步电机哈密顿模型的建立方法 |
CN115296330A (zh) * | 2022-06-17 | 2022-11-04 | 昆明理工大学 | 一种孤网模式下的哈密顿阻尼协同控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112464391B (zh) | 2022-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112464391B (zh) | 一种柴油发电机组哈密顿模型的建立方法 | |
CN109217362A (zh) | 一种双馈风机并网系统低频振荡扰动源定位系统及方法 | |
Rocha | A sensorless control for a variable speed wind turbine operating at partial load | |
CN112542855B (zh) | 一种双馈风力发电系统相量模型建模及仿真方法 | |
CN105488315B (zh) | 风力发电机并网产生次同步振荡的分析方法 | |
CN108011381A (zh) | 一种风储一体化系统调频控制方法 | |
CN105449699A (zh) | 双馈感应风电机组非线性分数阶自抗扰阻尼控制方法 | |
Rui et al. | Fractional‐order sliding mode control for hybrid drive wind power generation system with disturbances in the grid | |
CN111478365B (zh) | 一种直驱风电机组虚拟同步机控制参数的优化方法及系统 | |
CN115882762A (zh) | 构网型风电并网系统的频率优化控制方法 | |
Gao et al. | A fast high-precision model of the doubly-fed pumped storage unit | |
CN117175617A (zh) | 一种综合惯性控制下的双馈风机等效惯量评估方法及系统 | |
CN106294959B (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
Jia et al. | Research on ship power plant simulation system based on LabVIEW and MATLAB mixed programming | |
Shi et al. | Model predictive control of direct-driven surface-mounted permanent magnet synchronous generator based on active disturbance rejection controller | |
CN115800844A (zh) | 基于降阶pi观测器的永磁同步电机无模型滑模控制方法 | |
CN113783183B (zh) | 弱电网下双馈风机在故障穿越期间的暂态稳定性评估方法 | |
CN113612238A (zh) | 一种分析风电并网引发火电机组轴系小干扰振荡的改进阻尼转矩方法 | |
Zhang et al. | Research of coordination control system between nonlinear robust excitation control and governor power system stabilizer in multi-machine power system | |
CN106655261A (zh) | 针对电铁引起的负序对风电机组影响的抑制方法及装置 | |
Jizhen et al. | Dynamic modeling of wind turbine generation system based on grey-box identification with genetic algorithm | |
Huang et al. | Analysis and Control of Nonlinear Torsional Vibration of Direct-Drive Permanent Magnet Wind Turbine Shaft System | |
Yan et al. | Research on governor parameter optimization to suppress ultra-low frequency oscillation of power system caused by hydropower unit | |
CN113113908A (zh) | 适用于现代大电网频率响应的时域解析方法与系统 | |
Subudhi et al. | Sliding mode approach to torque and pitch control for a wind energy system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |