CN108845594A - 基于传递函数的振动谐波迭代控制方法 - Google Patents
基于传递函数的振动谐波迭代控制方法 Download PDFInfo
- Publication number
- CN108845594A CN108845594A CN201810598169.XA CN201810598169A CN108845594A CN 108845594 A CN108845594 A CN 108845594A CN 201810598169 A CN201810598169 A CN 201810598169A CN 108845594 A CN108845594 A CN 108845594A
- Authority
- CN
- China
- Prior art keywords
- vibration
- signal
- acceleration
- frequency
- transmission function
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000001133 acceleration Effects 0.000 claims abstract description 93
- 239000004615 ingredient Substances 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims description 28
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 8
- 238000006073 displacement reaction Methods 0.000 description 7
- 230000003321 amplification Effects 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D19/00—Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase
- G05D19/02—Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase characterised by the use of electric means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Vibration Prevention Devices (AREA)
- Apparatuses For Generation Of Mechanical Vibrations (AREA)
Abstract
本发明公开了一种基于传递函数的振动谐波迭代控制方法,包括计算当前输入量;驱动振动激励器产生振动加速度信号;检测输出振动加速度信号;计算谐波成分及波形失真度值;判断波形失真度值是否满足设定要求;计算与各阶谐波成分对应的输入补偿量;计算修正输入量,并代替当前输入量共7个步骤,通过对振动激励器输出振动信号中的谐波成分进行开环控制,达到降低振动信号波形失真度的目的。本发明公开的基于传递函数的振动谐波迭代控制方法所需控制系统组成结构简单、操作流程简便、适用范围广泛,具有较好的振动激励器输出振动信号波形失真控制稳定性及控制效果。
Description
技术领域
本发明属于振动谐波控制技术领域,具体涉及一种基于传递函数的振动谐波迭代控制方法。
背景技术
振动激励器广泛用于振动环境模拟、测振传感器标定及其它需产生振动激励信号的技术领域。为提高信噪比,振动激励器输出振动信号的位移幅值应随频率的降低而逐渐增大。但是,振动激励器中弹性回复装置、磁路结构等系统的非线性也将随位移幅值的增大而逐渐增大,进而导致输出振动信号叠加高次谐波成分,产生波形失真。输出振动信号的波形失真会严重制约振动激励器在需产生高精度单频振动激励信号场合的应用。
为降低振动激励器输出低频振动信号的波形失真度(控制谐波成分),通常需构建闭环反馈控制系统。反馈控制对波形失真的控制效果与待控制系统的频率响应特性及反馈深度(反馈通道上的放大系数)相关。一般情况下,振动激励器在低频存在较大的非线性特性,为实现较好的波形失真控制效果,通常需施加深度负反馈控制。而且,振动激励器一般存在相位滞后响应特性。随着反馈深度及相位滞后的增大,控制系统将出现自激振荡或正反馈等不稳定现象,严重影响控制效果。为达到更好的反馈控制效果,通常需设计结构复杂的控制器(如PID控制器),并通过对相关参数的反复调整完成系统校正。然而,引入控制器不可避免地会增加系统的复杂性及参数调整难度。
发明内容
为有效解决通过闭环反馈控制实现振动激励器输出振动信号波形失真控制中存在的系统结构复杂、参数调整难度大及稳定性要求高等问题,本发明提出基于传递函数的振动谐波迭代控制方法,通过对振动激励器输出振动信号中的谐波成分进行开环控制,达到降低振动信号波形失真度的目的。
基于传递函数的振动谐波迭代控制方法,具体步骤为:
1)计算当前输入量
在需施加波形失真控制的频率范围内,基于振动激励单元的加速度传递函数及待控制振动激励器需产生的目标振动加速度信号频率、幅值及相位要求,以该目标振动加速度信号作为加速度传递函数的输出量,计算得到相应的输入量频率、幅值及相位,并以该输入量作为当前输入量;
2)驱动振动激励器产生振动加速度信号
由程控信号源产生当前输入量,并输入功率放大器驱动振动激励器产生振动加速度信号;
3)检测输出振动加速度信号
采用振动加速度传感器及其放大器或适配器检测得到振动激励器的输出振动加速度信号,由数据采集卡采集该输出振动加速度信号并送入计算机;
4)计算谐波成分及波形失真度值
采用具有对信号实施FFT分析及失真度计算能力的谐波分析软件,计算得到输出振动加速度信号中包含的基频及各阶谐波成分的频率、幅值、相位及对应的波形失真度值;
5)判断波形失真度值是否满足设定要求
判断输出振动加速度信号的波形失真度值是否满足小于或等于设定的波形失真度要求,波形失真度要求指待控制振动激励器产生的输出振动加速度信号波形失真度的最大允许值,若满足,则保持当前输入量驱动振动激励器产生振动加速度信号;若不满足,则顺序执行步骤6)、7),之后返回步骤2);
6)计算与各阶谐波成分对应的输入补偿量
基于振动激励单元的加速度传递函数,分别以输出振动加速度信号中各阶谐波成分作为加速度传递函数的输入量,计算得到与各阶谐波成分对应的输入补偿量的频率、幅值及相位;
7)计算修正输入量,并代替当前输入量
将输出振动加速度信号各阶谐波对应的输入补偿量反相后与当前输入量叠加,得到修正输入量,并用修正输入量代替当前输入量。
所述的程控信号源为任意波形信号发生器。
步骤1)、6)所述的振动激励单元包括振动激励器、功率放大器、振动加速度传感器及其放大器或适配器。
步骤1)、6)所述的振动激励单元的加速度传递函数辨识步骤为:
(a)在一定频率范围内,选取多个测试频率点,由程控信号源和功率放大器驱动振动激励器产生各频率点对应的振动输出信号;该频率范围的下限值远远小于振动激励器需施加波形失真控制的下限频率,上限值远远大于振动激励器需施加波形失真控制的上限频率;
(b)由振动加速度传感器及其放大器或适配器检测得到振动激励器在各频率点的输出加速度信号,并由数据采集卡同步采集该输出加速度信号及程控信号源产生的输入信号,然后,将该输出加速度信号及输入信号送入计算机;
(c)分别计算各频率点对应的输出加速度信号与输入信号间的幅值比及相位差,进一步得到振动激励单元各频率点对应的加速度幅频及相频特性值;
(d)基于上述各选定频率点的加速度幅频及相频特性值,通过MATLAB系统辨识工具辨识得到振动激励单元的加速度传递函数。
本发明所述方法可实现对振动激励器输出振动信号中谐波成分的开环控制,所需控制系统组成结构简单、操作流程简便、适用范围广泛,具有较好的振动激励器输出振动信号波形失真控制稳定性及控制效果。
附图说明
图1为本发明振动激励单元的加速度传递函数辨识流程图。
图2为本发明基于传递函数的振动谐波迭代控制方法流程图。
具体实施方式
以下结合附图对本发明进行详细说明:
以采用振动加速度传感器检测振动激励器的输出振动信号为例,用于驱动振动激励器的功率放大器放大倍数及用于检测振动加速度输出信号的振动加速度传感器(包括放大器或适配器)灵敏度在整个测试频段内可以设置为常数。所以,可以将振动激励器、功率放大器、振动加速度传感器及其放大器或适配器整体定义为振动激励单元。进一步,基于如图1所示的振动激励单元的加速度传递函数辨识流程,辨识得到振动激励单元的加速度传递函数,具体步骤为:
(a)在一定频率范围内,选取多个测试频率点,由程控信号源和功率放大器驱动振动激励器产生各频率点对应的振动输出信号;该频率范围的下限值远远小于振动激励器需施加波形失真控制的下限频率,上限值远远大于振动激励器需施加波形失真控制的上限频率;
(b)由振动加速度传感器及其放大器或适配器检测得到振动激励器在各频率点的输出加速度信号,并由数据采集卡同步采集该输出加速度信号及程控信号源产生的输入信号,然后,将该输出加速度信号及输入信号送入计算机;
(c)分别计算各频率点对应的输出加速度信号与输入信号间的幅值比及相位差,进一步得到振动激励单元各频率点对应的加速度幅频及相频特性值;
(d)基于上述各选定频率点的加速度幅频及相频特性值,通过MATLAB系统辨识工具辨识得到振动激励单元的加速度传递函数。
基于上述振动激励单元的加速度传递函数,实现振动激励器输出振动加速度信号中谐波成分的开环控制。如图2所示,具体步骤为:
1)计算当前输入量
在需施加波形失真控制的频率范围内,基于振动激励单元的加速度传递函数及待控制振动激励器需产生的目标振动加速度信号频率、幅值及相位要求,以该目标振动加速度信号作为加速度传递函数的输出量,计算得到相应的输入量频率、幅值及相位,并以该输入量作为当前输入量;
2)驱动振动激励器产生振动加速度信号
由程控信号源产生当前输入量,并输入功率放大器驱动振动激励器产生振动加速度信号;
3)检测输出振动加速度信号
采用振动加速度传感器及其放大器或适配器检测得到振动激励器的输出振动加速度信号,由数据采集卡采集该输出振动加速度信号并送入计算机;
4)计算谐波成分及波形失真度值
采用具有对信号实施FFT分析及失真度计算能力的谐波分析软件,计算得到输出振动加速度信号中包含的基频及各阶谐波成分的频率、幅值、相位及对应的波形失真度值;
5)判断波形失真度值是否满足设定要求
判断输出振动加速度信号的波形失真度值是否满足小于或等于设定的波形失真度要求,波形失真度要求指待控制振动激励器产生的输出振动加速度信号波形失真度的最大允许值,若满足,则保持当前输入量驱动振动激励器产生振动加速度信号;若不满足,顺序执行步骤6)、7),之后返回步骤2);
6)计算与各阶谐波成分对应的输入补偿量
基于振动激励单元的加速度传递函数,分别以输出振动加速度信号中各阶谐波成分作为加速度传递函数的输入量,计算得到与各阶谐波成分对应的输入补偿量的频率、幅值及相位;
7)计算修正输入量,并代替当前输入量
将输出振动加速度信号各阶谐波对应的输入补偿量反相后与当前输入量叠加得到修正输入量,并由修正输入量代替当前输入量。
所述的程控信号源为任意波形信号发生器,可受计算机软件控制输出具有设定频率、幅值及相位的单频信号或多个单频信号的叠加信号。
步骤1)、6)所述的振动激励单元包括振动激励器、功率放大器、振动加速度传感器及其放大器或适配器。
所述的振动加速度传感器可以由振动速度传感器或振动位移传感器代替,当采用振动速度传感器时,以上所述振动激励单元的传递函数为速度传递函数;当采用振动位移传感器时,以上所述振动激励单元的传递函数为位移传递函数。
与所述的振动加速度传感器、振动速度传感器或振动位移传感器对应,基于步骤1)-7)所述的迭代控制过程,可分别实现在整个需施加波形失真控制的频段内将振动激励器输出振动加速度信号的波形失真度、输出振动速度信号的波形失真度或输出振动位移信号的波形失真度控制到设定范围内。
总体而言,本实施例基于振动激励单元的传递函数,构建振动激励器输出振动信号中谐波成分的补偿量,并叠加到输入量中,通过迭代控制实现对谐波成分的有效控制,从而达到降低振动激励器输出振动信号波形失真度的目的。
本实施例所述方法可实现对振动激励器输出振动信号中谐波成分的开环控制,所需控制系统组成结构简单、操作流程简便、适用范围广泛,具有较好的振动激励器输出振动信号波形失真控制稳定性及控制效果。
说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (4)
1.基于传递函数的振动谐波迭代控制方法,其特征在于:具体步骤为:
1)计算当前输入量
在需施加波形失真控制的频率范围内,基于振动激励单元的加速度传递函数及待控制振动激励器需产生的目标振动加速度信号频率、幅值及相位要求,以该目标振动加速度信号作为加速度传递函数的输出量,计算得到相应的输入量频率、幅值及相位,并以该输入量作为当前输入量;
2)驱动振动激励器产生振动加速度信号
由程控信号源产生当前输入量,并输入功率放大器驱动振动激励器产生振动加速度信号;
3)检测输出振动加速度信号
采用振动加速度传感器及其放大器或适配器检测得到振动激励器的输出振动加速度信号,由数据采集卡采集该输出振动加速度信号并送入计算机;
4)计算谐波成分及波形失真度值
采用具有对信号实施FFT分析及失真度计算能力的谐波分析软件,计算得到输出振动加速度信号中包含的基频及各阶谐波成分的频率、幅值、相位及对应的波形失真度值;
5)判断波形失真度值是否满足设定要求
判断输出振动加速度信号的波形失真度值是否满足小于或等于设定的波形失真度要求,波形失真度要求指振动激励器产生的输出振动加速度信号波形失真度的最大允许值,若满足,则保持当前输入量驱动振动激励器产生振动加速度信号;若不满足,则顺序执行步骤6)、7),之后返回步骤2);
6)计算与各阶谐波成分对应的输入补偿量
基于振动激励单元的加速度传递函数,分别以输出振动加速度信号中各阶谐波成分作为加速度传递函数的输入量,计算得到与各阶谐波成分对应的输入补偿量的频率、幅值及相位;
7)计算修正输入量,并代替当前输入量
将输出振动加速度信号各阶谐波对应的输入补偿量反相后与当前输入量叠加,得到修正输入量,并用修正输入量代替当前输入量。
2.根据权利要求1所述的基于传递函数的振动谐波迭代控制方法,其特征在于:程控信号源为任意波形信号发生器。
3.根据权利要求1所述的基于传递函数的振动谐波迭代控制方法,其特征在于:步骤1)、6)所述的振动激励单元包括振动激励器、功率放大器、振动加速度传感器及其放大器或适配器。
4.根据权利要求1所述的基于传递函数的振动谐波迭代控制方法,其特征在于:步骤1)、6)所述的振动激励单元的加速度传递函数辨识步骤为:
(a)在一定频率范围内,选取多个测试频率点,由程控信号源和功率放大器驱动振动激励器产生各频率点对应的振动输出信号;该频率范围的下限值远远小于振动激励器需施加波形失真控制的下限频率,上限值远远大于振动激励器需施加波形失真控制的上限频率;
(b)由振动加速度传感器及其放大器或适配器检测得到振动激励器在各频率点的输出加速度信号,并由数据采集卡同步采集该输出加速度信号及程控信号源产生的输入信号,然后,将该输出加速度信号及输入信号送入计算机;
(c)分别计算各频率点对应的输出加速度信号与输入信号间的幅值比及相位差,进一步得到振动激励单元各频率点对应的加速度幅频及相频特性值;
(d)基于上述各选定频率点的加速度幅频及相频特性值,通过MATLAB系统辨识工具辨识得到振动激励单元的加速度传递函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810598169.XA CN108845594B (zh) | 2018-06-12 | 2018-06-12 | 基于传递函数的振动谐波迭代控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810598169.XA CN108845594B (zh) | 2018-06-12 | 2018-06-12 | 基于传递函数的振动谐波迭代控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108845594A true CN108845594A (zh) | 2018-11-20 |
CN108845594B CN108845594B (zh) | 2021-05-28 |
Family
ID=64211675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810598169.XA Active CN108845594B (zh) | 2018-06-12 | 2018-06-12 | 基于传递函数的振动谐波迭代控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108845594B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
CN113359892A (zh) * | 2021-06-21 | 2021-09-07 | 歌尔股份有限公司 | 振动电机的传递函数生成方法、装置及存储介质 |
CN114001193A (zh) * | 2021-09-18 | 2022-02-01 | 上海华兴数字科技有限公司 | 作业机械的启停减振控制方法、装置及电子设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101180789A (zh) * | 2005-05-31 | 2008-05-14 | 三菱电机株式会社 | 电动机控制装置 |
CN101702090A (zh) * | 2009-10-19 | 2010-05-05 | 大连海事大学 | 一种基于卡尔曼滤波器的电液伺服系统波形再现控制方法 |
CN101702093A (zh) * | 2009-10-09 | 2010-05-05 | 江苏大学 | 一种振动结构辐射噪声的有源控制系统及控制方法 |
CN101741318A (zh) * | 2009-12-31 | 2010-06-16 | 中国航空工业集团公司北京长城计量测试技术研究所 | 正弦振动台失真度控制方法 |
CN103425051A (zh) * | 2013-08-21 | 2013-12-04 | 北京航空航天大学 | 一种磁悬浮旋转机械的不平衡量辨识与振动抑制控制系统 |
US20150312962A1 (en) * | 2007-10-01 | 2015-10-29 | Timothy Russell Mullen | Wireless electric heat trace and vibration control and monitoring system |
CN105959884A (zh) * | 2016-05-24 | 2016-09-21 | 陈菁 | 平面振膜复合伺服式扬声器系统及其控制方法 |
-
2018
- 2018-06-12 CN CN201810598169.XA patent/CN108845594B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101180789A (zh) * | 2005-05-31 | 2008-05-14 | 三菱电机株式会社 | 电动机控制装置 |
US20150312962A1 (en) * | 2007-10-01 | 2015-10-29 | Timothy Russell Mullen | Wireless electric heat trace and vibration control and monitoring system |
CN101702093A (zh) * | 2009-10-09 | 2010-05-05 | 江苏大学 | 一种振动结构辐射噪声的有源控制系统及控制方法 |
CN101702090A (zh) * | 2009-10-19 | 2010-05-05 | 大连海事大学 | 一种基于卡尔曼滤波器的电液伺服系统波形再现控制方法 |
CN101741318A (zh) * | 2009-12-31 | 2010-06-16 | 中国航空工业集团公司北京长城计量测试技术研究所 | 正弦振动台失真度控制方法 |
CN103425051A (zh) * | 2013-08-21 | 2013-12-04 | 北京航空航天大学 | 一种磁悬浮旋转机械的不平衡量辨识与振动抑制控制系统 |
CN105959884A (zh) * | 2016-05-24 | 2016-09-21 | 陈菁 | 平面振膜复合伺服式扬声器系统及其控制方法 |
Non-Patent Citations (1)
Title |
---|
刘一江等: "《基于迭代控制的电液振动台控制系统》", 《控制工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112505747A (zh) * | 2020-12-22 | 2021-03-16 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
CN112505747B (zh) * | 2020-12-22 | 2021-10-01 | 吉林大学 | 基于多信号发生器协同可控震源振动畸变抑制系统及方法 |
CN113359892A (zh) * | 2021-06-21 | 2021-09-07 | 歌尔股份有限公司 | 振动电机的传递函数生成方法、装置及存储介质 |
CN113359892B (zh) * | 2021-06-21 | 2022-06-21 | 歌尔股份有限公司 | 振动电机的传递函数生成方法、装置及存储介质 |
CN114001193A (zh) * | 2021-09-18 | 2022-02-01 | 上海华兴数字科技有限公司 | 作业机械的启停减振控制方法、装置及电子设备 |
CN114001193B (zh) * | 2021-09-18 | 2024-03-15 | 上海华兴数字科技有限公司 | 作业机械的启停减振控制方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108845594B (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108845594A (zh) | 基于传递函数的振动谐波迭代控制方法 | |
CN108827457A (zh) | 基于传递函数的振动幅值迭代控制方法 | |
CN108716974A (zh) | 基于传递函数的振动幅值迭代控制系统 | |
CN102221401B (zh) | 一种振动试验台和一种振动传感器的测试方法 | |
Wu et al. | Linear and nonlinear aeroelastic analysis frameworks for cable-supported bridges | |
CN101741318B (zh) | 正弦振动台失真度控制方法 | |
Rowley et al. | Model-based control of cavity oscillations. ii-system identification and analysis | |
CN102749873B (zh) | 电磁驱动式相控震源相位补偿控制系统及补偿控制方法 | |
CN111220263A (zh) | 一种检测马达谐振频率的方法和系统 | |
CN103346692A (zh) | 振动主动控制中压电作动器迟滞非线性频域补偿方法 | |
dos Santos et al. | On limit cycle oscillations of typical aeroelastic section with different preset angles of incidence at low airspeeds | |
CN109062036A (zh) | 基于传递函数的振动谐波迭代控制系统 | |
CN107797387A (zh) | 工件台运动控制方法 | |
CN106354016A (zh) | 一种金属管夯实高度的控制方法 | |
CN106053000A (zh) | 一种导弹与运载火箭模态试验方法 | |
CN101937198A (zh) | 一种振动状态控制装置及其控制方法 | |
CN104655111B (zh) | 一种激光探测器模拟方法及装置 | |
CN204807089U (zh) | 一种陀螺仪检测电路、陀螺仪及电子设备 | |
CN106989880A (zh) | 模拟地震波作用下构件受力测试方法与装置 | |
Wang et al. | Parameter tuning method for dither compensation of a pneumatic proportional valve with friction | |
CN104615126A (zh) | 一种基于控制系统的频率响应特性测试系统及方法 | |
CN107479376A (zh) | 基于模糊和比例微分控制动态切换的风洞支杆抑振方法 | |
CN102507225B (zh) | 一种舵系统动态测试方法 | |
CN103423239B (zh) | 基于自适应颤振的电液伺服系统死区的测定装置及方法 | |
CN113624498B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240104 Address after: Room 103, 1st Floor, Comprehensive Office Building, No. 512 Yingbin West Street, Huitong Industrial Park, Jinzhong Development Zone, Shanxi Comprehensive Reform Demonstration Zone, Jinzhong City, Shanxi Province, 030600 (multiple locations per photo) Patentee after: Shanxi Hydrogen Power Technology Co.,Ltd. Address before: 030024 No. 79 West Main Street, Wan Berlin District, Shanxi, Taiyuan, Yingze Patentee before: Taiyuan University of Technology |