CN112464391A - 一种柴油发电机组哈密顿模型的建立方法 - Google Patents

一种柴油发电机组哈密顿模型的建立方法 Download PDF

Info

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
Application number
CN202011106209.8A
Other languages
English (en)
Other versions
CN112464391B (zh
Inventor
钱晶
曾云
吕顺利
王芳芳
梅宏
李丹
于凤荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202011106209.8A priority Critical patent/CN112464391B/zh
Publication of CN112464391A publication Critical patent/CN112464391A/zh
Application granted granted Critical
Publication of CN112464391B publication Critical patent/CN112464391B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling 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)将柴油机拉格朗日函数转换为哈密顿函数:
Figure BDA0002726584580000011
θm为广义位移,对应的广义速度为
Figure BDA0002726584580000012
ω,写为矩阵形式为
Figure BDA0002726584580000013
定义广义动量为p=[p1 p2 p3]T,广义动量可以通过拉格朗日函数
Figure BDA0002726584580000014
求取,如下:
Figure BDA0002726584580000015
Figure BDA0002726584580000016
Figure BDA0002726584580000021
根据动力学理论,定义如下的哈密顿函数H(1)
Figure BDA0002726584580000022
其中:拉格朗日函数为
Figure BDA0002726584580000023
Figure BDA0002726584580000024
是轴运动速度标么值;
Figure BDA0002726584580000025
是质量标么值;
Figure BDA0002726584580000026
是弹簧刚度标么值;x1是执行器位移,各个量取标么值;
对上式的哈密顿函数H(1)求偏导可获得:
Figure BDA0002726584580000027
利用这 组变换关系将柴油机微分方程式转化为广义哈密顿形式,如下:
Figure BDA0002726584580000028
其中:H(1)为柴油机哈密顿函数;
Figure BDA0002726584580000029
为执行器运动部分质量标么值;
Figure BDA00027265845800000210
为执行器位移标么值;
Figure BDA00027265845800000211
为执行器移动速度标么值
Figure BDA00027265845800000212
Figure BDA00027265845800000213
是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角 速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,
Figure BDA00027265845800000214
为阻尼系数标么值;
Figure BDA00027265845800000215
TD、Tb、;xnl为空载位移;
Figure BDA00027265845800000216
为执行器的输入;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.”中给出的模型如下。
发电机哈密顿能量函数:
Figure BDA00027265845800000217
哈密顿方程为:
Figure BDA0002726584580000031
其中,δ为发电机功角(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为励磁控制输出,以上参数中没有标明单位均为标么值;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得;
将柴油机与发电机集合为一个系统,重新选取状态变量:
Figure BDA0002726584580000032
v3=δ,v4=ω1,v5=Eq′, 其中,由于发电机与柴油机是同轴运行,在将二者的能量函数合成时,分项:
Figure BDA0002726584580000033
只能取 一次;包含两个子系统哈密顿函数为:
Figure BDA0002726584580000034
其中,由于发电机与柴油机是同轴运行,在将二者的能量函数合成时,分项:
Figure BDA0002726584580000035
只能取 一次;
(3)建立柴油机哈密顿方程:
将步骤(2)中获得的柴油机哈密顿方程与已知的发电机哈密顿方程联合得到柴油发电机 组哈密顿方程:
Figure BDA0002726584580000036
其中:
Figure BDA0002726584580000037
Figure BDA0002726584580000038
Figure BDA0002726584580000041
Figure BDA0002726584580000042
所述柴油机拉格朗日拉格朗日函数的建立方法,具体步骤如下:
(1)执行器的拉格朗日函数的建立:
设柴油机电磁执行器的输入信号为u,输出轴位移为x,衔铁运动在轴上产生的电磁力F 为二者的函数,表示为F(x,u),假设初始执行器工作在输入信号为u0,位移为x0的稳定状态, 则电磁力表示为F(x0,u0),当输入信号变化Δu,相应轴位移变化Δx,电磁力变化为:
F(x0,u0)+ΔF=F(x0+Δx,u0+Δu)
将上式右边采用泰勒级数展开得到:
Figure BDA0002726584580000043
Figure BDA0002726584580000044
则:F(x,u)=F(x0,u0)+kxΔx+kuΔu;其中
Figure BDA0002726584580000045
表示电磁力F随输入u的 变化梯度,
Figure BDA0002726584580000046
表示电磁力F随输出轴位移x的变化梯度;
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能, 即:
Figure BDA0002726584580000047
其中:v是轴运动速度,m/s;m1是质量,kg;k1是弹簧刚度,N/m;x是执行器位移,mm;
(2)柴油机旋转部分轴系统的拉格朗日函数的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转 机械角位移θm(rad)为广义坐标,轴系统的拉格朗日函数L2等于轴的旋转动能:
Figure BDA0002726584580000051
其中:J为转动部分的转动惯量,kg.m2;ωm为机械角速度为,rad/s;
(3)柴油机拉格朗日函数的建立:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将步骤1和步骤2中所获得的 执行器和本体的能量函数L1和L2相加,可得到柴油机的拉格朗日函数:
Figure BDA0002726584580000052
其中: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)
将上式右边采用泰勒级数展开得到:
Figure BDA0002726584580000061
Figure BDA0002726584580000062
则:F(x,u)=F(x0,u0)+kxΔx+kuΔu;其中
Figure BDA0002726584580000063
表示电磁力F随输入u的 变化梯度,
Figure BDA0002726584580000064
表示电磁力F随输出轴位移x的变化梯度;
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能, 即:
Figure BDA0002726584580000065
其中:v是轴运动速度(m/s),m1是质量,k1是弹簧刚度(N/m),x是执行器位移(mm);
执行器在运动过程中受到阻尼,设阻尼系数为c1(N.s/m),阻尼与速度成正比,即大小为c1v, 所以轴上的外力Q为电磁力和阻尼力之差:Q=F(x,u)-c1v
经典的拉格朗日方程一般形式为:
Figure BDA0002726584580000066
将L1和Q代入其中可得到轴的运动 方程:
Figure BDA0002726584580000067
在平衡状态下,输入控制信号和执行器轴位移恒定,即u=u0,x=x0,增量为零:Δu=0, Δx=0,代入上式方程有k1x0=F(x0,u0),即弹簧力与电磁力相等,方程可表示为:
Figure BDA0002726584580000071
其中,kx1=kx-k1,Δx=x-x0,Δu=u-u0
在平衡点处上式仍然要满足,所以将平衡点坐标(u=u0,x=x0)代入上式得:kx1x0=-kuu0, 于是上式可以改写为如下形式:
Figure BDA0002726584580000072
(i)式即是执行器的拉格朗日方程;
(2)柴油机旋转部分拉格朗日模型的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转 机械角位移θm(rad)为广义坐标,系统的拉格朗日函数L2等于轴的旋转动能:
Figure BDA0002726584580000073
其中:J为转动部分的转动惯量(kg.m2),ωm为机械角速度为(rad/s);
设柴油机轴力矩为M1,若柴油机输出端连接发电机,设发电机电磁力矩为M2,它代表的 是在柴油机轴上产生的制动性力矩,再计及发电机阻尼绕组产生的阻尼力矩Md,则作用在柴 油机轴上的非保守广义外力(力矩)为:
Q2=M1-M2-Md
其中,柴油机轴上的力矩M1与轴位移和角速度为线性关系,可表示为:
M1=kωmmB)-a1(x-xnl)
其中,kω、a1为系数,a1=M1r/(xr-xnl),M1r是柴油机额定转矩,xr是额定扭矩对应的轴位移,xnl是空载位移,即a1是柴油机力矩随轴位移变化的增益系数或斜率;ωmB是额定机械角速度;
阻尼转矩与角速度也成正比:Md=Dp(ωmmB),其中,D为柴油机转动阻尼系数;p为发电机磁极对数;所以Q2=kωmmB)-a1(x-xnl)-M2-Dp(ωmmB)
将L2、Q2代入拉格朗日方程式
Figure BDA0002726584580000074
得到:
Figure BDA0002726584580000075
其中r0=kω-Dp;r1=-r0ωmB+a1xn1;ωmB为额定机械角速度,M2为发电机电磁力矩;
(3)求取柴油机拉格朗日函数和拉格朗日方程:
(a)柴油机拉格朗日函数:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将所获得的执行器和本体的能 量函数L1和L2相加,可得到柴油机的拉格朗日函数:
Figure BDA0002726584580000081
(b)柴油机拉格朗日方程:
根据上述推导,将柴油机拉格朗日方程整理为标准的形式如下,设x1为执行器位移,x2为 执行器运动速度v,即x2=dx1/dt,代入(i)式得:
Figure BDA0002726584580000082
式(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:将柴油机拉格朗日函数转换为哈密顿函数:
Figure BDA0002726584580000083
θm为广义位移,对应的广义速度为
Figure BDA0002726584580000084
ω,写为矩阵形式为
Figure BDA0002726584580000085
定义广义动量为p=[p1 p2 p3]T,广义动量可以通过拉格朗日函数
Figure BDA0002726584580000086
求取,如下:
Figure BDA0002726584580000087
Figure BDA0002726584580000088
Figure BDA0002726584580000089
根据动力学理论,定义如下的哈密顿函数H(1)
Figure BDA0002726584580000091
其中:从步骤1可知拉格朗日函数为
Figure BDA0002726584580000092
Figure BDA0002726584580000093
是轴运动速度标么值;
Figure BDA0002726584580000094
是质 量标么值;
Figure BDA0002726584580000095
是弹簧刚度标么值;x1是执行器位移,各个量取标么值;
对上式哈密顿函数H(1)求偏导可获得:
Figure RE-GDA0002923202590000096
利用这组 变换关系将微分方程式转化为广义哈密顿形式,如下:
Figure BDA0002726584580000097
其中:H(1)为柴油机哈密顿函数;
Figure BDA0002726584580000098
为执行器运动部分质量标么值;
Figure BDA0002726584580000099
为执行器位移标么值;
Figure BDA00027265845800000910
为执行器移动速度标么值;
Figure BDA00027265845800000911
是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,
Figure BDA00027265845800000912
为阻尼系数 标么值;
Figure BDA00027265845800000913
TD、Tb、;xnl为空载位移;
Figure BDA00027265845800000916
为执行器的输入;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.”中给出的模型如下:
发电机哈密顿能量函数:
Figure BDA00027265845800000914
哈密顿方程为:
Figure BDA00027265845800000915
其中,δ为发电机功角(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中的,将柴油机与发电机集合为一个系统,重新选取状态变量:
Figure BDA0002726584580000101
v3=δ,v4=ω1,v5=Eq′,两个子系统哈密顿函数合并,其中,由于发电机与柴油机是同轴运行,在将二者的能量函数合成时,分项:
Figure BDA0002726584580000102
只能取一次,合并后的函数为:
Figure BDA0002726584580000103
将步骤(2)中获得的柴油机哈密顿方程(iv)与已知的发电机哈密顿方程(v)联立得到柴油发电 机组哈密顿方程:
Figure BDA0002726584580000104
其中:
Figure BDA0002726584580000105
Figure BDA0002726584580000106
Figure BDA0002726584580000111
Figure BDA0002726584580000112
其中,H为柴油机哈密顿函数;
Figure BDA0002726584580000113
为执行器运动部分质量标么值;
Figure BDA0002726584580000114
为执行器位移标么 值;
Figure BDA0002726584580000115
为执行器移动速度标么值;
Figure BDA0002726584580000116
是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角 速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,
Figure BDA0002726584580000117
为阻尼系数标么值;
Figure BDA0002726584580000118
TD、Tb、;xnl为空载位移;
Figure BDA0002726584580000119
为执行器的输入标么值;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)将柴油机拉格朗日函数转换为哈密顿函数:
Figure FDA0002726584570000011
θm为广义位移,对应的广义速度为
Figure FDA0002726584570000012
ω,写为矩阵形式为
Figure FDA0002726584570000013
定义广义动量为p=[p1 p2 p3]T,广义动量可以通过拉格朗日函数
Figure FDA00027265845700000122
求取,如下:
Figure FDA0002726584570000014
Figure FDA0002726584570000015
Figure FDA0002726584570000016
根据动力学理论,定义如下的哈密顿函数H(1)
Figure FDA0002726584570000017
其中:拉格朗日函数为
Figure FDA0002726584570000018
Figure FDA0002726584570000019
是质量标么值;
Figure FDA00027265845700000110
是执行器运动速度标么值;
Figure FDA00027265845700000111
是弹簧刚度标么值;x1是执行器位移标么值,各个量取标么值;
对上式的哈密顿函数H(1)求偏导可获得:
Figure FDA00027265845700000112
Figure FDA00027265845700000113
利用这组变换关系将柴油机微分方程式转化为广义哈密顿形式,如下:
Figure FDA00027265845700000114
其中:H(1)为柴油机哈密顿函数;
Figure FDA00027265845700000115
为执行器运动部分质量标么值;
Figure FDA00027265845700000116
为执行器位移标么值;
Figure FDA00027265845700000117
为执行器移动速度标么值;
Figure FDA00027265845700000118
是柴油机执行部分弹簧刚度标么值,ωB为发电机额定角速度;Tj为惯性时间常数;ω1为发电机角速度相对额定角速度的变化值,
Figure FDA00027265845700000119
为阻尼系数;
Figure FDA00027265845700000120
Figure FDA00027265845700000121
TD、Tb、;xnl为空载位移;
Figure FDA00027265845700000123
为执行器的输入标么值;m2为发电机电磁力矩;D为发电机等效阻尼系数;以上参数均为标么值;
(2)建立发电机的哈密顿能量函数:
发电机哈密顿能量函数为:
Figure FDA0002726584570000021
哈密顿方程为:
Figure FDA0002726584570000022
其中,δ为发电机功角(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)建立柴油发电机组哈密顿方程:
将柴油机与发电机合成为一个系统,重新选取状态变量:
Figure FDA0002726584570000023
包含两个子系统哈密顿函数为:
Figure FDA0002726584570000024
其中,由于发电机与柴油机是同轴运行,在将二者的能量函数合成时,分项:
Figure FDA0002726584570000025
只能取一次;
将步骤(2)中获得的柴油机哈密顿方程与已知的发电机哈密顿方程联合得到柴油发电机组哈密顿方程:
Figure FDA0002726584570000026
其中:
Figure FDA0002726584570000027
Figure FDA0002726584570000031
Figure FDA0002726584570000032
Figure FDA0002726584570000033
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)
将上式右边采用泰勒级数展开得到:
Figure FDA0002726584570000034
Figure FDA0002726584570000035
则:F(x,u)=F(x0,u0)+kxΔx+kuΔu;其中
Figure FDA0002726584570000036
表示电磁力F随输入u的变化梯度,
Figure FDA0002726584570000041
表示电磁力F随输出轴位移x的变化梯度;
执行器机械系统是典型的弹簧质量系统,则执行器的拉格朗日函数L1为动能减去势能,即:
Figure FDA0002726584570000042
其中:v是轴运动速度,m/s;m1是质量,kg;k1是弹簧刚度,N/m;x是执行器位移,mm;
(2)柴油机旋转部分轴系统的拉格朗日函数的建立:
柴油机本体旋转部分主要是指轴系旋转部分,功率从轴上输出,表现为轴力矩,取轴旋转机械角位移θm(rad)为广义坐标,轴系统的拉格朗日函数L2等于轴的旋转动能:
Figure FDA0002726584570000043
其中:J为转动部分的转动惯量,kg.m2;ωm为机械角速度为,rad/s;
(3)柴油机拉格朗日函数的建立:
柴油机拉格朗日函数由执行器和本体两部分能量函数构成,将步骤1和步骤2中所获得的执行器和本体的能量函数L1和L2相加,可得到柴油机的拉格朗日函数:
Figure FDA0002726584570000044
其中:L为拉格朗日函数;m1为执行器运动部分质量,kg;k1是执行器弹簧刚度,N/m;x1为执行器位移,mm;x2为执行器速度,m/s;J为转动部分的转动惯量,kg.m2;ωm为转动部分的机械角速度,rad/s;以上的各个参数通过设备额定参数、实验直接获得或经过已知参数计算获得。
CN202011106209.8A 2020-10-15 2020-10-15 一种柴油发电机组哈密顿模型的建立方法 Active CN112464391B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (10)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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