CN109873186B - 质子交换膜燃料电池准二维模型建立方法 - Google Patents
质子交换膜燃料电池准二维模型建立方法 Download PDFInfo
- Publication number
- CN109873186B CN109873186B CN201910162441.4A CN201910162441A CN109873186B CN 109873186 B CN109873186 B CN 109873186B CN 201910162441 A CN201910162441 A CN 201910162441A CN 109873186 B CN109873186 B CN 109873186B
- Authority
- CN
- China
- Prior art keywords
- flow channel
- node
- current density
- concentration
- model
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000000446 fuel Substances 0.000 title claims abstract description 23
- 239000012528 membrane Substances 0.000 title claims abstract description 12
- 238000009826 distribution Methods 0.000 claims abstract description 37
- 230000002441 reversible effect Effects 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 239000012495 reaction gas Substances 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000012546 transfer Methods 0.000 claims abstract description 6
- 239000007789 gas Substances 0.000 claims description 38
- 238000009792 diffusion process Methods 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 10
- 239000001301 oxygen Substances 0.000 claims description 10
- 229910052760 oxygen Inorganic materials 0.000 claims description 10
- 239000000376 reactant Substances 0.000 claims description 8
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 7
- 239000001257 hydrogen Substances 0.000 claims description 7
- 229910052739 hydrogen Inorganic materials 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 230000004907 flux Effects 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 abstract description 3
- 239000003054 catalyst Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000012827 research and development Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000003197 catalytic effect Effects 0.000 description 1
- 238000010351 charge transfer process Methods 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003487 electrochemical reaction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
Images
Classifications
-
- 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
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Landscapes
- Fuel Cell (AREA)
Abstract
本发明公开了一种建立质子交换膜燃料电池的准二维模型建立方法,该模型的建立主要包括:求解流道内传质方程、求解沿流道的局部电流密度和电池输出电压、顺逆流进气模式下,准二维模型的迭代求解等。由于电池运行时由于反应气体的消耗,沿流道方向存在反应气浓度的下降,催化层反应速率的变化使得电池沿流道方向各节点电流密度不同。模型考虑了沿流道方向的实际反应气体浓度分布,能够更好地反映电池运行的局部特性,使模型物理过程更完整,提高了模型的准确性和应用价值,特别有助于电堆及系统级燃料电池模型的开发,推动燃料电池实际应用领域的发展。模型的求解方法具有更好的收敛性和计算效率,有助于模型快速稳定地对电池运行情况进行预测。
Description
技术领域
本发明属于电化学燃料电池领域,具体涉及一种对质子交换膜燃料电池准二维模型的建立方法。
背景技术
质子交换膜燃料电池具有效率高、零排放、低温启动等优点,因此被广泛认为有很好的发展潜力和应用前景。相比燃料电池的实验测试,计算机仿真能够有效减少研发周期和成本,并且避免了实验测试的诸多干扰因素,因此将其作为研发手段之一。
目前,燃料电池模型的建立方式主要有两种:一种是通过现有的计算流体力学(CFD)软件,分析燃料电池内部的质量、动量、能量和电荷传输过程;另一种是在一定的假设基础上对电池的控制方程进行解析求解。然而,这两种模型建立方式都存在各自的弊端:前者虽然考虑的物理化学过程完整,且仿真工具较成熟,但模型的计算时间长、效率低,因此很难适用于燃料电池的研发;后者的计算效率高,能满足企业仿真工作者对于快速预测电池性能和运行情况的需求,但因模型做了不同程度的简化,其仿真结果的准确性受到一定影响,这直接决定了模型的实际应用价值。
为了解决上述问题,一种思路是开发出电化学反应机理和传输过程完整的解析模型。目前的解析模型多为一维模型,仅考虑了垂直于电池极板方向上的变化,而忽略了流道内部的气体浓度分布。工程实践证明,流道的结构设计和流道内进气方式的改变对于电池性能的影响是显著的,因此开发出一种考虑了流道配气的电池模型显得尤为重要。
发明内容
本发明的目的是提供一种质子交换膜燃料电池的准二维模型建立方法,能够对流道内的反应气体传输方程以及沿流道方向的局部电流密度分布进行求解,且能够对顺逆流两种不同配气方式下的电池运行情况进行分析。
建立质子交换膜燃料电池的准二维模型建立方法的具体步骤如下:
(1)求解流道内传质方程
稳态情况下,流道内反应气体的完整传输过程可以由对流扩散方程表示:
式中,x和y分别表示垂直电池极板方向和沿流道方向;CGDL(mol m-3)为气体扩散层(GDL)内的反应气摩尔浓度;Cch为流道内的气体摩尔浓度,由于不考虑流道截面上的气体浓度差异,此处Cch仅在y方向上变化;uch(m s-1)为流道内的气体流速;dch(m)为流道深度;和Dch(m3 s-1)分别为GDL和流道(ch)内气体的有效扩散系数。
等号左侧为垂直极板方向上(x方向)通过流道进入气体扩散层的反应气体流通量;等号右侧两项依次为流道内气体的对流项和扩散项;
从流道上游到下游,流道被划分为N个相同尺寸的节点,各个节点对应的电流密度不同,但均满足:
上式中,Ik为各节点对应的电流密度;qk(mol m-2s-1)为单位面积反应速率;F为法拉第常数;阳极氢气反应的n=2,阴极氧气反应的n=4;
式中,Lch为流道特征长度,取值为电池的极板长度;流道内对流强度和扩散强度的相对大小可以由传质佩克莱数Pe表征:
当流道内对流占主导时,忽略流道内气体扩散的作用,而将方程简化为:
由上式求出的流道气体浓度分布为线性分布,结果适用于每个流道节点内;表示各个节点入口界面处的气体浓度。在节点1入口处有:代入微分方程(1-5)可求得节点1内气体浓度分布;后续节点入口边界条件(即入口浓度)可以由上一节点求得的出口浓度给定;各节点内特征浓度可由该节点进出口平均浓度表示;
(2)求解沿流道的局部电流密度和电池输出电压
流道各节点的局部电流密度Ik和输出电压V均为待求量,共有N+1个变量,其满足两个确定关系:电池总的平均电流密度为给定值I;电池各节点处阴阳极电势差相等。这两个关系由下面的公式(2-1)和(2-2)表示:
联立公式(2-1)和(2-2)可以求出各节点的局部电流密度和输出电压;
(3)顺逆流进气模式下准二维模型的迭代求解
前述步骤1和2中的流道内反应气体浓度和电流密度Ik相互耦合,采用迭代方法进行求解,求解过程如下:
(3.1)迭代变量初始化
迭代初值的选取对迭代收敛性影响很大,采用给定的总平均电流密度I先代入一维模型(即公式(2-2))计算出一维稳态解,用计算出的稳态解对准二维模型中各节点参数进行初始化,各节点局部电流密度Ik用I初始化;
(3.2)代入各节点电流密度Ik,根据步骤(1)更新流道内各节点的平均浓度。
(3.4)由求解出的物理量,根据步骤(2)求出下一个迭代步的电流密度Ik *。
模型在顺逆流进气模式下,均默认y轴正方向为阳极入口到出口,阴阳极流道节点均由y轴正方向进行统一划分和编号。顺流模式下,阴极流道内氧气浓度的边界条件为第1节点的入口浓度已知,而在逆流模式下,该边界条件应调整为第5节点的入口浓度已知;
通过采用该方法,能够对顺逆流进气模式下的质子交换膜燃料电池性能和沿流道方向的物理量分布进行求解,且该求解方法能保证计算的收敛性和高效性;
由于电池运行时由于反应气体的消耗,沿流道方向必然存在反应气浓度的下降,从而导致催化层反应速率的变化,表现为电池沿流道方向各节点电流密度不同;
在工程实践上,实际的单电池流道为长流道,流道长度及流道结构对电池的影响往往很显著;另外,工程上往往更关注燃料电池电堆及系统的特性,而在燃料电池电堆层面,电池的局部问题会更加突出,且需要考虑流道的配气特性;
本发明的特点以及产生的有益效果是:
(1)该准二维模型的建立增加了流道方向的考虑,使模型物理过程更完整,提高了模型的准确性和应用价值;
(2)该准二维模型的建立有助于电堆及系统级燃料电池模型的开发,推动燃料电池实际应用领域的发展;
(3)模型的求解方法具有更好的收敛性和计算效率,采用的初始化方法,相比根据经验随意给定电池各参数初始值,模型求解的收敛性更好,且收敛速度更快,有助于模型快速稳定地对电池运行情况进行预测;
附图说明
图1为质子交换膜燃料电池结构和进气示意图,流道划分为五个节点;
图3为逆流进气模式下,在总的平均电流密度I取不同值时的局部电流密度分布图;
图4为顺逆流进气模式下不同阳极化学计量比STa时的阳极流道氢气浓度分布图;
图5为顺逆流进气模式下不同阴极化学计量比STc时的阴极流道氧气浓度分布图;
图2~图5均为利用本发明计算得到的实施效果;
具体实施方式
以下结合说明书附图并通过具体实施例,对本发明的方法和步骤进行详细的说明;
电池结构和相关的基本参数如下:总的平均电流密度I=1.5A cm-2;氢气/氧气在流道内的有效扩散系数分别为:1.055×10-4/2.652×10-5m2 s-1;进气压力为一个大气压,进气充分加湿条件下,氢气/氧气进气浓度为36.07/7.58mol m-3;进气化学计量比STa/STc为:1.5/2.0,2.0/2.0,2.5/2.0,2.0/1.5,2.0/2.5;流道特征长度Lch=0.25m;流道深度dch=1mm;
(1)选取迭代初始值,将流道划分为5个节点,每个节点局部电流密度Ik均由I初始化为1.5A cm-2,各节点代入局部电流密度对垂直极板方向的物理量分布进行求解;
(2)依次代入节点电流密度Ik,对流道内传质方程进行求解,更新流道内的浓度。
上式中,阳极氢气反应的n=2,阴极氧气反应的n=4,法拉第常数F=96487C mol-1;表示各个节点入口界面处的气体浓度;uch(m s-1)为流道内的气体流速,由给定的化学计量比和进气浓度确定;顺流进气模式下,流道边界条件为:在节点1入口处有逆流模式下,该边界条件为第5节点的入口浓度已知:=1处有根据该公式和流道边界条件可以得到流道各节点内的气体浓度分布。
(3)由各节点电流密度和更新后的流道内浓度计算对应x方向上的物理量分布;
(4)由求解出的物理量分布,求解下一个迭代步的电流密度Ik *,求解方程如下;
根据以上步骤可以得到电池不同运行工况下的局部电流密度分布和反应气体分布规律;
本发明所提出的质子交换膜燃料电池准二维模型的建立方法,考虑了沿流道方向的实际反应气体浓度分布,能够更好地反映电池运行的局部特性,如沿流道方向的局部电流密度和流道内反应气浓度分布等特征;
实施例图2显示了在STa/STc=2.0/2.0的顺流进气模式下,随着平均电流密度I增大时电池各节点的局部电流密度分布图;图中纵坐标为无量纲化的局部电流密度大小为各节点电流密度和流道入口节点电流密度的比值可以看到,顺流模式下电池上游的电流密度均高于下游,且分布均近似线性分布;随着平均电流密度I的提高,电流密度分布曲线更平缓;
实施例图3给出了STa/STc=2.0/2.0的逆流进气模式下,随着平均电流密度I增大时的局部电流密度分布;随着平均电流密度I的增大,电流密度分布逐渐趋于均匀,各节点电流密度均趋近于阴极入口节点值(节点5);逆流模式下的电流密度分布不为线性分布,且当电流密度I较高时呈现出“中间高两边低”的特征;相比顺流模式的电流密度分布(值最大相差0.5),逆流模式下的电流密度分布更均匀(值最大差异不超过0.2);当I增大到0.3A cm-2后,各节点电流密度的差异变得很小,此时值最大仅有0.02的差异;
实施例图4给出了在不同的阳极进气化学计量比时,顺逆流进气模式下流道内氢气的浓度分布,分布均近线性分布;从流道上游到下游区域,由于氢气的反应消耗,流道内氢气的浓度逐渐下降,在出口处达到最小值。阳极化学计量比越高,供气速率越大,因此流道内氢气浓度也越高;
实施例图5给出了在不同的阴极进气化学计量比时,顺逆流进气模式下流道内氧气的浓度分布;顺流模式下,节点1为流道入口节点,因此从节点1入口到节点5出口氧气的浓度逐渐下降;逆流模式下,节点5为流道入口节点,因此节点5到节点1氧气浓度逐渐降低;随着阴极化学计量比的提高,阴极流道内供气速率上升,因此氧气浓度也更高;
通过以上实施过程,可以量化说明本发明所建立的质子交换膜燃料电池准二维模型,可以用于快速高效地分析电池的局部特性。
Claims (1)
1.一种建立质子交换膜燃料电池的准二维模型建立方法,其特征是:所述模型建立的具体步骤如下:
(1)求解流道内传质方程
稳态情况下,流道内反应气体的完整传输过程由对流扩散方程表示:
式中,x和y分别表示垂直电池极板方向和沿流道方向;CGDL为GDL内的反应气摩尔浓度;Cch为流道内的气体摩尔浓度,不考虑流道截面上的气体浓度差异,此处Cch仅在y方向上变化;uch为流道内的气体流速;dch为流道深度;为GDL内气体的有效扩散系数,式中等号左侧为垂直极板的x方向上通过流道进入气体扩散层的反应气体流通量;等号右侧两项依次为流道内气体的对流项和扩散项,
从流道上游到下游,流道被划分为N个相同尺寸的节点,各个节点对应的电流密度不同,但均满足:
式中,Ik为各节点对应的电流密度;qk为单位面积反应速率;F为法拉第常数;阳极氢气反应的n=2,阴极氧气反应的n=4,
式中,Lch为流道特征长度,取值为电池的极板长度;流道内对流强度和扩散强度的相对大小可以由传质佩克莱数Pe表征;
当流道内对流占主导时,忽略流道内气体扩散的作用,将方程简化为:
由(1-5)式求出的流道气体浓度分布为线性分布,结果适用于每个流道节点内;在节点1入口处有:代入微分方程(1-5)可求得节点1内气体浓度分布;后续节点入口边界条件可以由上一节点求得的出口浓度给定;各节点内特征浓度可由该节点进出口平均浓度表示;
(2)求解沿流道的局部电流密度和电池输出电压
流道各节点的局部电流密度Ik和输出电压V作为待求量,共有N+1个变量,变量满足:电池总的平均电流密度为给定值I;电池各节点处阴阳极电势差相等,这两个关系由下式(2-1)和(2-2)表示:
对公式(2-1)和(2-2)联立求解得出各节点的局部电流密度和输出电压;
(3)顺逆流进气模式下,准二维模型的迭代求解
前述步骤(1)和(2)中的流道内反应气体浓度和电流密度Ik相互耦合,采用迭代方法进行求解:
(3.1)迭代变量初始化
采用给定的总平均电流密度I先代入一维模型计算出一维稳态解,用计算出的稳态解对准二维模型中各节点参数进行初始化,各节点局部电流密度Ik用I初始化;
(3.2)代入各节点电流密度Ik,根据步骤(1)更新流道内浓度边界条件;
(3.3)由各节点电流密度和更新后的流道内浓度边界条件计算对应x方向上的物理量分布;
(4)由求解出的物理量,根据步骤(2)求出下一个迭代步的电流密度Ik *;
模型在顺逆流进气模式下,默认沿流道方向y是阳极入口到出口,阴阳极流道节点均由阳极入口到出口方向进行划分和编号;顺流模式下,阴极流道内氧气浓度的边界条件为已知的第1节点的入口浓度,而在逆流模式下,该边界条件应调整为已知的第5节点的入口浓度;由此实现顺逆流进气模式下质子交换膜燃料电池沿流道方向的局部性能参数和物理量的求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910162441.4A CN109873186B (zh) | 2019-03-05 | 2019-03-05 | 质子交换膜燃料电池准二维模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910162441.4A CN109873186B (zh) | 2019-03-05 | 2019-03-05 | 质子交换膜燃料电池准二维模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109873186A CN109873186A (zh) | 2019-06-11 |
CN109873186B true CN109873186B (zh) | 2022-05-27 |
Family
ID=66919732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910162441.4A Active CN109873186B (zh) | 2019-03-05 | 2019-03-05 | 质子交换膜燃料电池准二维模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109873186B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110993035B (zh) * | 2019-11-22 | 2022-07-15 | 华中科技大学鄂州工业技术研究院 | 用于质子交换膜电解池特性分析的二维稳态模型建立方法 |
CN111199110B (zh) * | 2020-01-13 | 2022-04-01 | 清华大学 | 燃料电池低温启动性能预测方法及系统 |
CN113921847B (zh) * | 2021-08-23 | 2023-03-31 | 东风汽车集团股份有限公司 | 一种燃料电池堆流体流动分析建模方法及装置 |
CN114357806B (zh) * | 2022-03-11 | 2022-06-17 | 中国汽车技术研究中心有限公司 | 基于物质流接口的燃料电池电堆的双模式仿真方法和设备 |
CN116314952B (zh) * | 2022-09-08 | 2024-04-09 | 清华大学 | 燃料电池设计方案评估方法、装置、设备、介质和产品 |
CN118409217A (zh) * | 2024-04-22 | 2024-07-30 | 天津大学 | 基于电化学机理构建燃料电池阻抗谱分析模型的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281970A (zh) * | 2007-04-04 | 2008-10-08 | 通用汽车环球科技运作公司 | 改善水管理的显微结构化的燃料电池元件 |
CN106848351A (zh) * | 2016-12-26 | 2017-06-13 | 天津大学 | 质子交换膜燃料电池性能预测模型建立的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7829230B2 (en) * | 2007-07-17 | 2010-11-09 | Gm Global Technology Operations, Inc. | Method for optimizing diffusion media with spatially varying mass transport resistance |
-
2019
- 2019-03-05 CN CN201910162441.4A patent/CN109873186B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281970A (zh) * | 2007-04-04 | 2008-10-08 | 通用汽车环球科技运作公司 | 改善水管理的显微结构化的燃料电池元件 |
CN106848351A (zh) * | 2016-12-26 | 2017-06-13 | 天津大学 | 质子交换膜燃料电池性能预测模型建立的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109873186A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109873186B (zh) | 质子交换膜燃料电池准二维模型建立方法 | |
CN111199110B (zh) | 燃料电池低温启动性能预测方法及系统 | |
CN111180767B (zh) | 具有流体不均匀分配效应的燃料电池堆建模方法 | |
Chen et al. | Multi-criteria assessment and optimization study on 5 kW PEMFC based residential CCHP system | |
Xiao et al. | Control strategy optimization of electrolyte flow rate for all vanadium redox flow battery with consideration of pump | |
CN110335646B (zh) | 基于深度学习-预测控制的车用燃料电池过氢比控制方法 | |
CN114447378B (zh) | 一种质子交换膜燃料电池的参数优化方法 | |
CN112909303A (zh) | 质子交换膜燃料电池瞬态实时模型的建立方法 | |
CN109902435A (zh) | 质子交换膜燃料电池建模方法、存储介质及计算机设备 | |
Lan et al. | Analysis and optimization of high temperature proton exchange membrane (HT-PEM) fuel cell based on surrogate model | |
Hailun et al. | Performance analysis and prediction of ejector based hydrogen recycle system under variable proton exchange membrane fuel cell working conditions | |
WO2023216150A1 (zh) | 一种燃料电池的热管理方法 | |
Li et al. | A data-driven framework for performance prediction and parameter optimization of a proton exchange membrane fuel cell | |
Mao et al. | Multi-objective optimization of gradient porosity of gas diffusion layer and operation parameters in PEMFC based on recombination optimization compromise strategy | |
CN115470581A (zh) | 燃料电池气体流道优化设计方法、系统、电子设备及介质 | |
Guo et al. | Numerical study of Tesla valve flow field on proton exchange membrane fuel cell performance | |
Chen et al. | Development of flow field plates based on asymmetric leaf structure for PEM fuel cells | |
Lei et al. | Effect of cathode channel structure on the performance of open-cathode proton exchange membrane fuel cell in dry environments | |
Pan et al. | In-plane gradient design of flow fields enables enhanced convections for redox flow batteries | |
Fu et al. | Optimization of sinusoidal wave-like channel design for HT-PEMFCs based on genetic algorithm | |
CN111769312B (zh) | 一种基于压力补偿的燃料电池供给路解耦控制方法 | |
CN117878358A (zh) | 一种燃料电池输入控制方法、装置、设备及燃料电池系统 | |
Zhang et al. | Validation of 3D multi-physics equivalent resistance network model with flow field for VRFB stack and battery scale-up analysis | |
CN117195510A (zh) | 一种基于解析解的液流电池流道选择优化方法 | |
CN111682244B (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 |