CN110083911B - 一种电磁振动能量回收系统的建模优化方法 - Google Patents

一种电磁振动能量回收系统的建模优化方法 Download PDF

Info

Publication number
CN110083911B
CN110083911B CN201910319079.7A CN201910319079A CN110083911B CN 110083911 B CN110083911 B CN 110083911B CN 201910319079 A CN201910319079 A CN 201910319079A CN 110083911 B CN110083911 B CN 110083911B
Authority
CN
China
Prior art keywords
electromagnetic vibration
recovery system
energy recovery
vibration energy
damping
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.)
Active
Application number
CN201910319079.7A
Other languages
English (en)
Other versions
CN110083911A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201910319079.7A priority Critical patent/CN110083911B/zh
Publication of CN110083911A publication Critical patent/CN110083911A/zh
Application granted granted Critical
Publication of CN110083911B publication Critical patent/CN110083911B/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/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Vibration Prevention Devices (AREA)
  • Control Of Electric Motors In General (AREA)
  • Feedback Control In General (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

一种电磁振动能量回收系统的建模优化方法,先将电磁式振动能量回收系统抽象为单自由度弹簧质量阻尼的动力学系统模型;再测量电磁式振动能量回收系统的内部电阻、反电动势系数、弹簧刚度、负重质量、增速器速比和系统机械阻尼;然后以外接电路电阻,激励频率比为优化对象,建立以外接电路阻尼和激励频率为参数的电磁式振动能量回收系统的俘能模型;再对俘能模型进行最优解析求解,计算最优外接电路阻抗,最优激励频率比,最优弹簧刚度,以及增速器速比;最后根据最优解析结果,对实际俘能系统进行参数优化调节;本发明具有实用性和准确性的优点。

Description

一种电磁振动能量回收系统的建模优化方法
技术领域
本发明涉及能量收集系统优化技术领域,具体涉及一种电磁振动能量回收系统的建模优化方法。
背景技术
电磁式振动能量俘获系统的结构通常有四部分:安装基础、配重、弹簧和能量转换部件,其利用配重在安装基础激励下的受迫振动驱动能量转换部件中的磁铁和线圈产生相对运动来产生电能。
传统的电磁振动能量回收系统的建模优化是以激励频率比和全电路等效阻尼为主要优化对象,根据经验调节激励频率比为共振点附近某值,根据模型得到俘能功率图像寻求最优外接电阻。但是在实际操作中,经验带入的人为不确定因素较多,全电路等效阻尼中的内部电阻不便于优化调参,根据传统模型和图像寻优有很大的误差。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种电磁振动能量回收系统的建模优化方法,具有实用性和准确性的优点。
为了达到上述目的,本发明采取的技术方案为:
一种电磁振动能量回收系统的建模优化方法,包括以下步骤:
第一步,将电磁振动能量回收系统抽象为单自由度弹簧质量阻尼的动力学系统模型;
第二步,测量电磁振动能量回收系统的系统参数,系统参数包括内部电阻、反电动势系数、弹簧刚度、负重质量、增速器速比和系统机械阻尼;
第三步,基于第一步抽象的动力学系统模型,以外接电路电阻,激励频率比为优化对象,建立以外接电路阻尼和激励频率为参数的电磁振动能量回收系统的俘能模型;
第四步,对俘能模型进行最优解析求解,代入第二步实测的系统参数,计算最优外接电路阻抗,最优激励频率比,最优弹簧刚度,以及增速器速比;
第五步,根据第四步得到的最优解析结果,对电磁振动能量回收系统进行参数优化调节。
所述的第三步中建立以外接电路阻尼和激励频率为参数的电磁式振动能量回收系统的俘能模型,表示为:
Figure BDA0002034081320000021
式中:P——俘获能量;cR——外接电路等效阻尼;ω——激励频率;m——负重质量;X0——激励位移;k——悬挂弹簧刚度;cm——机械阻尼;cr——内阻等效阻尼;
其中电路各阻尼为:
Figure BDA0002034081320000022
式中:R——外接电路总电阻;K——电机反电动势系数;r——电路总内阻;rd——旋转半径;N——增速器速比;
对俘能模型做无量纲化处理有:
Figure BDA0002034081320000031
式中:ωn——系统固有频率;ξm——机械阻尼比;α——激励频率比。
所述的第四步的具体计算为:
Figure BDA0002034081320000032
Figure BDA0002034081320000033
式中:ξ——系统总阻尼比。
所述的电磁振动能量回收系统的外接电路阻抗只显示电阻特性,不显示感抗或者容抗特性,或者感抗和容抗影响忽略不计。
所述的电磁振动能量回收系统的建模优化方法不仅适用于简谐位移激励下单自由度电磁振动能量回收系统,还适用于随机激励下电磁振动能量回收系统,或者多激励多自由度系统。
本发明的有益效果为:
1)本发明将电磁振动能量回收系统的俘能模型直接表达为实际中便于调节的激励频率比和外接电路阻抗的函数,可以直接应用于电磁振动能量回收系统的调试优化。
2)本发明避免了人为不确定因素的影响,可以准确的得到不同实际条件下最优激励频率比和最优外接电路阻抗。
3)本发明简化了最优功率点外部电路阻抗和激励频率比的优化过程,使得解析更加准确可靠易懂。
附图说明
图1是本发明的流程图。
图2是本发明实施例电磁振动能量回收系统的结构示意图。
图3是本发明实施例抽象为单自由度弹簧质量阻尼的动力学系统模型。
图4是本发明实施例的动力学系统模型仿真曲线。
图5是本发明实施例的实验阻抗匹配曲线。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细描述。
参照图1,一种电磁振动能量回收系统的建模优化方法,包括以下步骤:
第一步,将电磁振动能量回收系统抽象为单自由度弹簧质量阻尼的动力学系统模型,参照图2,图2是本实施例电磁振动能量回收系统的结构示意图,图2中齿轮、齿条将往复直线运动转变为旋转运动,通过传动轴经过增速器输入发电机发电;参照图3,图3是本实施例的动力学系统模型,其中系统等效质量为m,弹簧刚度为k,系统总阻尼为c,其中机械阻尼为cm,电路内阻等效阻尼cr,电路外阻等效阻尼cR,位移激励为X0=Asinωt,内部电路阻抗为r,外部电路阻抗为R;
第二步,测量电磁振动能量回收系统的系统参数,系统参数包括内部电阻r=2Ω、反电动势系数K=0.005V/(r/min)、弹簧刚度k=5000N/m、负重质量m=8.2kg和系统机械阻尼cm=8N/(m/s),齿轮动力输入半径rd=0.018m、增速器速比N=25。
第三步,基于第一步抽象的动力学系统模型,以外接电路电阻,激励频率比为优化对象,建立以外接电路阻尼和激励频率为参数的电磁振动能量回收系统的俘能模型;
电磁振动能量回收系统利用阻尼作用从外界获取能量,其中由外接电路阻尼俘获的能量能被有效利用,所以根据基本物理原理能够建立以外接电路阻尼和激励频率为参数的电磁式振动能量回收系统的俘能模型,表示为:
Figure BDA0002034081320000051
式中:P——俘获能量;cR——外接电路等效阻尼;ω——激励频率;m——负重质量;X0——激励位移;k——悬挂弹簧刚度;cm——机械阻尼;cr——内阻等效阻尼;
其中电路各阻尼为:
Figure BDA0002034081320000052
式中:R——外接电路总电阻;K——电机反电动势系数;r——电路总内阻;rd——旋转半径;N——齿轮箱速比;
对俘能模型做无量纲化处理有:
Figure BDA0002034081320000061
式中:ωn——系统固有频率;ξm——机械阻尼比;α——激励频率比;
第四步,对俘能模型进行最优解析求解,代入第二步实测的系统参数,计算最优外接电路阻抗,最优激励频率比,最优弹簧刚度,以及增速器速比;
Figure BDA0002034081320000062
Figure BDA0002034081320000063
式中:ξ——系统总阻尼比;
代入第二步实测数据计算外阻R=47.73Ω,参照图3,图3是本实施例的动力学系统模型仿真曲线,可以看出理论结果,仿真结果与实验结果相近在误差范围内;最优激励频率比计算类似;
第五步,根据第四步得到的最优解析结果,对电磁振动能量回收系统进行参数优化调节,参照图4,图4是本实施例的实验阻抗匹配曲线,其中实线为插值曲线,可以看到在50Ω附近实验得到最大输出功率。
所述的电磁振动能量回收系统的建模优化方法不仅适用于简谐位移激励下单自由度电磁振动能量回收系统,还适用于随机激励下电磁振动能量回收系统,或者多激励多自由度系统。

Claims (4)

1.一种电磁振动能量回收系统的建模优化方法,其特征在于,包括以下步骤:
第一步,将电磁振动能量回收系统抽象为单自由度弹簧质量阻尼的动力学系统模型;
第二步,测量电磁式振动能量回收系统的系统参数,系统参数包括内部电阻、反电动势系数、弹簧刚度、负重质量、增速器速比和系统机械阻尼;
第三步,基于第一步抽象的动力学系统模型,以外接电路电阻,激励频率比为优化对象,建立以外接电路阻尼和激励频率为参数的电磁式振动能量回收系统的俘能模型;
第四步,对俘能模型进行最优解析求解,代入第二步实测的系统参数,计算最优外接电路阻抗,最优激励频率比,最优弹簧刚度,以及增速器速比;
第五步,根据第四步得到的最优解析结果,对电磁振动能量回收系统进行参数优化调节;
所述的第三步中建立以外接电路阻尼和激励频率为参数的电磁式振动能量回收系统的俘能模型,表示为:
Figure FDA0002531381650000011
式中:P——俘获能量;cR——外接电路等效阻尼;ω——激励频率;m——负重质量;X0——激励位移;k——悬挂弹簧刚度;cm——机械阻尼;cr——内阻等效阻尼;
其中电路各阻尼为:
Figure FDA0002531381650000021
式中:R——外接电路总电阻;K——电机反电动势系数;r——电路总内阻;rd——旋转半径;N——增速器速比;
对俘能模型做无量纲化处理有:
Figure FDA0002531381650000022
式中:ωn——系统固有频率;ξm——机械阻尼比;α——激励频率比。
2.根据权利要求1所述的一种电磁振动能量回收系统的建模优化方法,其特征在于:所述的第四步的具体计算为:
Figure FDA0002531381650000023
Figure FDA0002531381650000024
式中:ξ——系统总阻尼比。
3.根据权利要求1所述的一种电磁振动能量回收系统的建模优化方法,其特征在于:所述的电磁振动能量回收系统的外接电路阻抗只显示电阻特性,不显示感抗或者容抗特性,或者感抗和容抗影响忽略不计。
4.根据权利要求1所述的一种电磁振动能量回收系统的建模优化方法,其特征在于:所述的电磁振动能量回收系统的建模优化方法不仅适用于简谐位移激励下单自由度电磁振动能量回收系统,还适用于随机激励下电磁振动能量回收系统,或者多激励多自由度系统。
CN201910319079.7A 2019-04-19 2019-04-19 一种电磁振动能量回收系统的建模优化方法 Active CN110083911B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910319079.7A CN110083911B (zh) 2019-04-19 2019-04-19 一种电磁振动能量回收系统的建模优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910319079.7A CN110083911B (zh) 2019-04-19 2019-04-19 一种电磁振动能量回收系统的建模优化方法

Publications (2)

Publication Number Publication Date
CN110083911A CN110083911A (zh) 2019-08-02
CN110083911B true CN110083911B (zh) 2020-10-27

Family

ID=67415808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910319079.7A Active CN110083911B (zh) 2019-04-19 2019-04-19 一种电磁振动能量回收系统的建模优化方法

Country Status (1)

Country Link
CN (1) CN110083911B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114760553A (zh) * 2021-01-11 2022-07-15 深圳市韶音科技有限公司 一种优化骨传导耳机工作状态的方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8866317B2 (en) * 2012-01-17 2014-10-21 Schlumberger Technology Corporation Broadband vibrational energy harvesting
CN105975673B (zh) * 2016-05-01 2019-06-25 上海大学 振动型能量采集器的非线性动力学参数辨识方法
CN109130757B (zh) * 2018-07-16 2020-10-27 西安交通大学 一种馈能式半主动悬架变阻尼系统与控制方法

Also Published As

Publication number Publication date
CN110083911A (zh) 2019-08-02

Similar Documents

Publication Publication Date Title
CN101949753B (zh) 一种高速柔性转子动平衡方法
CN101750216B (zh) 次同步振荡引起的汽轮发电机轴系疲劳损伤在线分析方法
CN111106783B (zh) 一种信号制作方法、信号制作装置、振动马达及触屏设备
CN102183366A (zh) 滚动轴承振动测量和故障分析装置及方法
CN104267597B (zh) 超精密运动平台机械谐振的抑制方法
CN109033488B (zh) 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法
CN105426644A (zh) 模态阻尼识别方法、装置和系统
CN104133176A (zh) 抽油机电动机动态负荷模拟加载系统及模拟加载方法
CN110083911B (zh) 一种电磁振动能量回收系统的建模优化方法
CN105429545A (zh) 一种基于软测量的无传感器永磁同步电机控制方法
Sang et al. Resonance control strategy for a slider crank WEC power take-off system
CN103308706A (zh) 一种叶片泵转速检测方法及其装置
CN112327957B (zh) 转子低频振动多阶线谱的控制方法及装置
CN103117693B (zh) 一种无需对转速微分运算的风轮机模拟器及其控制方法
CN102545762A (zh) 直驱式波浪发电系统的控制方法
CN105043667B (zh) 转子不平衡振动信号幅值、相位实时计算方法
Bonanno et al. A multiphysics approach to the design of a seawave energy conversion system
CN104200027B (zh) 一种考虑低电压穿越特性的风电场等效建模方法与装置
CN110094194A (zh) 电功图计算油井产液量方法
CN107798205B (zh) 双馈感应风力发电机组轴系模型参数的单独辨识方法
CN113890442B (zh) 一种粒子群滑模观测器风力发电机转子电流故障检测方法
CN104022705A (zh) 通过磁阻力补偿来减小拼接式伺服电机力矩波动的方法
CN113392473B (zh) 一种基于傅里叶级数的功图转换优化方法
CN105004993A (zh) 一种直线步进电机力频特性测试装置和方法
CN103412139B (zh) 一种同步发电机转速测量方法

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