CN103366048B - 建立车辆与地面耦合的一体式动力学模型的方法 - Google Patents
建立车辆与地面耦合的一体式动力学模型的方法 Download PDFInfo
- Publication number
- CN103366048B CN103366048B CN201310257524.4A CN201310257524A CN103366048B CN 103366048 B CN103366048 B CN 103366048B CN 201310257524 A CN201310257524 A CN 201310257524A CN 103366048 B CN103366048 B CN 103366048B
- Authority
- CN
- China
- Prior art keywords
- ground
- coordinate
- matrix
- formula
- vehicle
- 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.)
- Expired - Fee Related
Links
Landscapes
- Vehicle Body Suspensions (AREA)
Abstract
本发明涉及一种建立车辆和地面耦合的一体式动力学模型的方法,该方法将车辆的多体动力学模型与地面的有限单元模型都集成到一个模型内,它一方面从数学上涵盖了车辆与地面的力学关系;另一方面,从求解过程中也反映出它们之间的耦合关系,再以此计算地面的坐标和总应变,然后再根据地面的力学特性、屈服方程、流动准则、地面的变形状态去计算地面的弹性力矢量Qs等,最后再将这些弹性力矢量带回式六中进行新一轮的计算,通过这样的循环迭代去完成全时域的车辆和地面的所有变形、运动和受力状态等的计算过程。实现了与连续体地面力学模型进行一体化耦合,真正解决了车辆与地面的动力学耦合问题,适用于解决地面的大变形情形,如塑形变形问题。
Description
技术领域
本发明涉及一种建立车辆和地面耦合的一体式动力学模型的方法,属于建立车辆与地面耦合动力学模型的方法技术领域。
背景技术
车辆的研发过程中,需要充分考虑车辆与不同地面之间的耦合特性,从而使得车辆能够适应各种地面的特性,减少颠簸,打滑等情况。
车辆—地面耦合动力学模型是分析和研究车辆动力特性、舒适性等的重要基础。
目前所有商业多体动力学软件都不能真正地与连续体地面力学模型进行一体化耦合,只能简单地用一些力学中的经验公式或模型(比如BEKKER公式或模型)进行力和位移的计算,而且这些商业多体动力学软件都不适合大变形的地面情形。具体来说:
1)目前所有商业多体动力学软件都不能与连续体地面力学模型进行一体化耦合,它们通常需要两个独立的计算程序来解决车辆与地面的耦合问题,因此这些商业多体动力学软件只能在两个软件间交换状态变量和力信息,却不能统一处理必须满足多体动力学算法中有关位置、速度、加速度要求的代数约束方程,这样就不可能真正解决车辆与地面的动力学耦合问题。
2)而且这些商业多体动力学软件都不适合解决地面的大变形情形,如塑形变形问题。
发明内容
本发明需要解决的技术问题是:目前所有商业多体动力学软件都不能与连续体地面力学模型进行一体化耦合,它们通常需要两个独立的计算程序来解决车辆与地面的耦合问题,因此这些商业多体动力学软件只能在两个软件间交换状态变量和力信息,却不能统一处理必须满足多体动力学算法中有关位置、速度、加速度要求的代数约束方程,这样就不可能真正解决车辆与地面的动力学耦合问题。
本发明采取以下技术方案:
一种建立车辆与地面耦合的一体式动力学模型的方法,包括以下步骤:
a)建立绝对节点坐标系:
单元j上任意点位置矢量rj表示成全局坐标系XYZ成rj=Sj(xj,yj,zj)ej(t),在该方程中xj,yj,和zj是单元的空间坐标,Sj是形状函数矩阵,ej是时刻t时单元节点坐标矢量,节点坐标矢量ejk在节点k的表达式是式一:
广义连续体力学方法计算格林-拉格朗日应变张量ε=(JTJ-I)/2,这里J是位置矢量斜率矩阵,它在节点k的表达式是式二:
左柯西-格林变形张量表示为这里下标e、p分别表示弹性、塑形;对于剑桥粘土地面模型,其弹性柯西-格林变形张量表示为,
b)建立地面有限单元或变形体动力方程:
对于地面有限单元或变形体,实际功原理可以表示为式三:
这里V是单元体积,ρ是质量密度,r任意点的位置向量,fb是体力向量,式三中的第二项用广义内力表示为式四:
这里δe绝对节点坐标有限单元上节点坐标的变化值,Qs是广义内力,这样式三能够进一步导出为运动方程,即式五:
其中,M是固定不变的系统质量矩阵,Qs是广义内力矩阵,Qe是单元外加节点力矩阵;
c)建立一体化车辆-地面耦合模型的运动方程:
建立一体化车辆-地面耦合运动方程的增量形式,并用方程表示为式六:方程中下标r,f和a分别表示相对坐标,弹性坐标和绝对节点坐标,Mrr,Mrf,Mfr,Mff是浮动坐标公式中的分惯性矩阵,Maa是绝对节点坐标中的固定系统质量矩阵,Cq是约束点雅各比矩阵,λ拉格朗日乘子矩阵,Qr,Qf,和Qa分别是相对坐标,弹性坐标和绝对节点坐标中的广义力矩阵,Qc是二次速度矩阵,在相对坐标方程中广义坐标qr和qf用来描述经历较小变形的刚体和柔性体运动,在绝对节点坐标中的矢量qa用来描述经历较大变形和塑形变形的柔性体运动,矢量qa包括所有ANCF单元的节点坐标;质量矩阵Maa包括绝对节点坐标中的地面单元和车辆部件的质量矩阵,质量矩阵Maa通过乔莱斯基坐标变换成统一的质量矩阵;使用乔莱斯基转换矩阵Bc和,节点坐标e表示成乔莱斯基坐标p形式的e=Bcp;广义力矩阵Qa包括车辆-地面耦合的相互作用中的内力矩阵Qs和节点力矩阵Qe;
d)求解所述式六的方程:车体能够确定加速度矢量和和拉格朗日乘子λ;地面坐标就是有限单元节点上的坐标,即qa=e,加速度矢量用来求地面的坐标e和速度r=S(x,y,z)e(t),地面的坐标根据式二:
来求地面的总应变量,ε=(JTJ-I)/2,总应变量包括弹性应变εe和塑形应变εp2个部分,它们分别对应于Je和Jp;对应于总变形、弹性变形和塑形变形的右柯西-格林变形张量分别表示为Cr=JTJ,因此弹性柯西-格林应变张量表示为对于各向同性材料来说,左柯西-格林变形张量Cl表示为Cl=JJT,其弹性柯西-格林变形张量根据公式计算出来;这样,体积应变值是偏应变矢量是其中δ=[1 1 1]T;偏应变值为
柯西应力向量σK的主方向与弹性左变形向量的主方向是一样的,柯西应力张量计算表示为式七:这里,
这里ψ是储能函数,即P0是屈服曲面上的硬化参数,是弹性压缩比,α,μ0是常数;同时,2阶皮奥拉-基尔霍夫应力张量可以表示为σP2=J-1σKJ-1T,该应力张量与应变张量一起去计算式四中的绝对节点坐标力矢量Qs。由于地面的特性、屈服方程、流动准则都已包含在上述的公式中,因此,上述方程,即式六所表示的结构允许将所有地面的模型系统地代入到复杂车辆使用的多体动力学模型中。
通过式三,将车辆的多体动力学模型与地面的有限单元模型都集成到一个方程内,它一方面从数学上涵盖了车辆与地面的力学关系,如惯性矩阵,质量矩阵,广义力矩阵,约束点雅各比矩阵等。另一方面,从求解过程中也反映出它们之间的耦合关系,即,通过先求解车辆的加速度矢量,再以此计算地面的坐标和总应变,然后再根据地面的力学特性、屈服方程、流动准则、地面的变形状态(弹性或塑形变形)去计算地面的弹性力矢量Qs等,最后再将这些弹性力矢量带回式三中进行新一轮的计算,通过这样的循环迭代去完成全时域的车辆和地面的所有变形、运动和受力状态等的计算过程。
由于地面的大变形或塑形变形,所以单元的描述不仅要有位移坐标还要有转动坐标,即位移斜率坐标,从这一点来说,绝对节点坐标方法都具备,而其他方法不具备。
本发明的有益效果在于:
1)提供了建立车辆与地面耦合的一体式动力学模型的方法,将车辆的多体动力学模型与地面的有限单元模型都集成到一个模型内,实现了与连续体地面力学模型进行一体化耦合。
2)真正解决了车辆与地面的动力学耦合问题。
3)适用于解决地面的大变形情形,如塑形变形问题。
具体实施方式
下面结合具体实施例对本发明进一步说明。
本发明的目的是这样实现的:
1.建立绝对节点坐标系
绝对节点坐标单元不是采用无限小的旋转或有限旋转作节点坐标,而是用节点处的绝对斜率和位移作节点坐标。单元j上任意点位置矢量rj可以表示成全局坐标系XYZ成rj=Sj(xj,yj,zj)ej(t),在该方程中xj,yj,和zj是单元的空间坐标,Sj是形状函数矩阵,ej是时刻t时单元节点坐标矢量。节点坐标矢量ejk在节点k的表达式是:
对于完全参数化的绝对节点坐标单元不仅不排除使用弹性力学的公式,而且允许使用广义连续体力学方法计算Green-Lagrange(格林-拉格朗日)应变张量ε=(JTJ-I)/2,这里J是位置矢量斜率矩阵,它在节点k的表达式是:
根据绝对节点坐标方法的描述,左Cauchy-Green(柯西-格林)变形张量Cl可以表示为Cl=JJT。在大应变的塑性变形公式中,可以使用乘法分解式J=JeJp,这里Je相应于弹性变形的位置矢量斜率矩阵,Jp相应于塑性变形的位置矢量斜率矩阵。同样,左Cauchy-Green(柯西-格林)变形张量可表示为这里下标e、p分别表示弹性、塑形。对于Cam-Clay地面模型,其弹性Cauchy-Green(柯西-格林)变形张量可表示为,完全参数化的绝对节点坐标板单元和体单元能保证在节点处位置斜率的连续性以及车体碾压而产生的地面几何形状的变化。
2.建立地面有限单元或变形体动力方程
对于地面有限单元或变形体,实际功原理可以表示为:
这里V单元体积,ρ是质量密度,r任意点的位置向量,fb是体力向量。方程(3)中的第二项可以用广义内力表示为
这里δe绝对节点坐标有限单元上节点坐标的变化值,Qs是广义内力。这样方程(3)可以进一步导出为下列运动方程:
式中,M是固定不变的系统质量矩阵,Qs是广义内力矩阵,Qe是单元外加节点力矩阵。
3.建立一体化车辆-地面耦合模型的运动方程
由于ANCF有限单元可以将包含约束节点和大变形在内的所有单元和约束的统一在一个方程内,这样就能建立一个统一的求解车辆多体动力学模型和地面力学模型相互耦合的计算程序,而不是过去的两个独立的计算程序。根据节点方程和运动方程,一体化车辆-地面耦合运动方程的增量形式可表示如下:
式中下标r,f和a分别表示相对坐标,弹性坐标和绝对节点坐标,Mrr,Mrf,Mfr,Mff是浮动坐标公式中的分惯性矩阵,Maa是绝对节点坐标中的固定系统质量矩阵,Cq是约束点雅各比矩阵,λ拉格朗日乘子矩阵,Qr,Qf,和Qa分别是相对坐标,弹性坐标和绝对节点坐标中的广义力矩阵,Qc是二次速度矩阵,也就是,在相对坐标方程中广义坐标qr和qf用来描述经历小变形的刚体和柔性体运动,在绝对节点坐标中的矢量qa用来描述经历大变形和塑形变形的柔性体运动,矢量qa包括所有ANCF单元的节点坐标。质量矩阵Maa包括绝对节点坐标中的地面单元和车辆部件的质量矩阵,该矩阵可以通过Cholesky(乔莱斯基)坐标变换成统一的质量矩阵,它是一个优化的稀疏矩阵。使用Cholesky(乔莱斯基)转换矩阵Bc和,节点坐标e可以表示成Cholesky(乔莱斯基)坐标p形式的e=Bcp。广义力矩阵Qa包括车辆-地面耦合的相互作用中的内力矩阵Qs和节点力矩阵Qe。由于地面的特性、屈服方程、流动准则都已包含在上述的公式中,因此,方程(6)的结构允许将所有地面的模型系统地代入到复杂车辆使用的多体动力学模型中。
通过式三,将车辆的多体动力学模型与地面的有限单元模型都集成到一个方程内,它一方面从数学上涵盖了车辆与地面的力学关系,如惯性矩阵,质量矩阵,广义力矩阵,约束点雅各比矩阵等。另一方面,从求解过程中也反映出它们之间的耦合关系,即,通过先求解车辆的加速度矢量,再以此计算地面的坐标和总应变,然后再根据地面的力学特性、屈服方程、流动准则、地面的变形状态(弹性或塑形变形)去计算地面的弹性力矢量Qs等,最后再将这些弹性力矢量带回式三中进行新一轮的计算,通过这样的循环迭代去完成全时域的车辆和地面的所有变形、运动和受力状态等的计算过程。
由于地面的大变形或塑形变形,所以单元的描述不仅要有位移坐标还要有转动坐标,即位移斜率坐标,从这一点来说,绝对节点坐标方法都具备,而其他方法不具备。
上述实施例仅为本发明的优选实施例,并不用来限制本发明的保护范围,本领域的普通技术人员在本实施例的启发下,可以做出修改和变化,均在本发明的保护范围之内。
Claims (1)
1.一种建立车辆与地面耦合的一体式动力学模型的方法,其特征在于,包括以下步骤:
a)建立绝对节点坐标系:
单元j上任意点位置矢量rj表示成全局坐标系xyz为rj=Sj(xj,yj,zj)ej(t),在rj=Sj(xj,yj,zj)ej(t)中xj,yj,和zj是单元的空间坐标,Sj是形状函数矩阵,ej是时刻t时单元节点坐标矢量,节点坐标矢量ejk在节点k的表达式是式一:
采用广义连续体力学方法计算格林-拉格朗日应变张量ε=(JTJ-I)/2,这里J是位置矢量斜率矩阵,它在节点k的表达式是式二:
左柯西-格林变形张量表示为其中上标e、p分别表示弹性、塑形;对于剑桥粘土地面模型,其弹性柯西-格林变形张量表示为,
b)建立地面有限单元或变形体动力方程:
对于地面有限单元或变形体,实际功原理可以表示为式三:
这里V是单元体积,ρ是质量密度,r’是任意点的位置向量,fb是体力向量,式三中的第二项用广义内力表示为式四:
这里δe是绝对节点坐标有限单元上节点坐标的变化值,Qs是广义内力矩阵,这样式三能够进一步导出为运动方程,即式五:
其中,M是固定不变的系统质量矩阵,Qs是广义内力矩阵,Qe是单元外加节点力矩阵;
c)建立一体化车辆-地面耦合模型的运动方程:
建立一体化车辆-地面耦合运动方程的增量形式,并用方程表示为式六:方程中下标r,f和a分别表示相对坐标,弹性坐标和绝对节点坐标,Mrr,Mrf,Mfr,Mff是浮动坐标公式中的分惯性矩阵,Maa是绝对节点坐标中的固定系统质量矩阵,Cq是约束点雅各比矩阵,λ是拉格朗日乘子矩阵,Qr,Qf,和Qa分别是相对坐标,弹性坐标和绝对节点坐标中的广义力矩阵,Qc是二次速度矩阵,在相对坐标方程中广义坐标qr和qf用来描述经历较小变形的刚体和柔性体运动,在绝对节点坐标中的矢量qa用来描述经历较大变形和塑形变形的柔性体运动,矢量qa包括所有ANCF单元的节点坐标;质量矩阵Maa包括绝对节点坐标中的地面单元和车辆部件的质量矩阵,质量矩阵Maa通过乔莱斯基坐标变换成统一的质量矩阵;使用乔莱斯基转换矩阵Bc和节点地面坐标e’表示成乔莱斯基坐标p形式的e′=Bcp;广义力矩阵Qa包括车辆-地面耦合的相互作用中的广义内力矩阵Qs和单元外加节点力矩阵Qe;
d)求解所述式六的方程:车体能够确定加速度矢量和和拉格朗日乘子矩阵λ;地面坐标就是有限单元节点上的坐标,即qa=e’,加速度矢量用来求地面的坐标e’和速度r=S(x,y,z)e(t),地面的坐标根据式二:
来求地面的格林-拉格朗日应变张量,ε=(JTJ-I)/2,格林-拉格朗日应变张量包括弹性应变εe和塑形应变εp2个部分,它们分别对应于Je和Jp;对应于总变形、弹性变形和塑形变形的右柯西-格林变形张量分别表示为Cr=JTJ,因此弹性柯西-格林应变张量表示为对于各同性材料来说,左柯西-格林变形张量Cl表示为Cl=JJT,其弹性柯西-格林变形张量根据公式计算出来;这样,体积应变值是偏应变矢量是其中δ=[1 1 1]T;偏应变值为
柯西应力向量σK的主方向与弹性左变形向量的主方向是一样的,柯西应力张量计算表示为式七:这里,
这里ψ是储能函数,即P0是屈服曲面上的硬化参数,是弹性压缩比,α,μ0是常数;同时,2阶皮奥拉-基尔霍夫应力张量可以表示为σP2=J-1σKJ-1T,使用2阶皮奥拉-基尔霍夫应力张量σP2=J-1σKJ-1T和格林-拉格朗日应变张量ε=(JTJ-I)/2计算式四中的广义内力矩阵Qs。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310257524.4A CN103366048B (zh) | 2013-06-25 | 2013-06-25 | 建立车辆与地面耦合的一体式动力学模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310257524.4A CN103366048B (zh) | 2013-06-25 | 2013-06-25 | 建立车辆与地面耦合的一体式动力学模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103366048A CN103366048A (zh) | 2013-10-23 |
CN103366048B true CN103366048B (zh) | 2017-06-06 |
Family
ID=49367379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310257524.4A Expired - Fee Related CN103366048B (zh) | 2013-06-25 | 2013-06-25 | 建立车辆与地面耦合的一体式动力学模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103366048B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292016A (zh) * | 2017-06-15 | 2017-10-24 | 中车唐山机车车辆有限公司 | 仿真数据处理方法及装置 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107220421A (zh) * | 2017-05-18 | 2017-09-29 | 北京理工大学 | 一种空间复杂柔性结构多体系统动力学建模与计算方法 |
CN116306043B (zh) * | 2023-05-23 | 2023-07-28 | 北京云道智造科技有限公司 | 一种多体系统的求解方法、装置、电子设备及仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727527A (zh) * | 2009-12-28 | 2010-06-09 | 中国农业大学 | 动力学分析软件与有限元分析软件间数据自动转换的方法 |
CN102799740A (zh) * | 2012-07-30 | 2012-11-28 | 株洲南车时代电气股份有限公司 | 基于有限元的屏蔽门导靴三维滑动摩擦特性仿真分析方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7096166B2 (en) * | 2002-01-22 | 2006-08-22 | The Hong Kong Polytechnic University | Method of modeling muscular tissue with active finite elements |
JP2013057638A (ja) * | 2011-09-09 | 2013-03-28 | Sumitomo Rubber Ind Ltd | ゴム材料のシミュレーション方法 |
-
2013
- 2013-06-25 CN CN201310257524.4A patent/CN103366048B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727527A (zh) * | 2009-12-28 | 2010-06-09 | 中国农业大学 | 动力学分析软件与有限元分析软件间数据自动转换的方法 |
CN102799740A (zh) * | 2012-07-30 | 2012-11-28 | 株洲南车时代电气股份有限公司 | 基于有限元的屏蔽门导靴三维滑动摩擦特性仿真分析方法 |
Non-Patent Citations (4)
Title |
---|
3D Road Surface Digital Modeling in Time Domain for Virtual Proving Ground of Vehicles;Yonglin Zhang等;《Fuzzy System and Knowledge Discovery (FSKD),2010 Seventh International Conference on》;20100812;第6卷;全文 * |
Simulation Research for Quarter Vehicle ABS on Complex Surface Based on PID Control;Qiang Fu等;《Consumer Electronics, Communications and Networks(CECNet),2012 2nd International Conference on》;20120423;全文 * |
考虑地面变形特性的车辆地面耦合系统的建模与仿真;张晓阳等;《机械工程学报》;20091231;第45卷(第12期);全文 * |
车辆与路基路面耦合的动力学分析及其仿真;董慧;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(月刊)》;20080215;第2008年卷(第02期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292016A (zh) * | 2017-06-15 | 2017-10-24 | 中车唐山机车车辆有限公司 | 仿真数据处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103366048A (zh) | 2013-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Manko | A general model of legged locomotion on natural terrain | |
Sun et al. | Dynamic modeling of cable towed body using nodal position finite element method | |
Peiret et al. | Multibody system dynamics interface modelling for stable multirate co-simulation of multiphysics systems | |
Jain et al. | Roams: Planetary surface rover simulation environment | |
CN111730605B (zh) | 一种机器人姿态控制方法、装置、可读存储介质及机器人 | |
CN103366048B (zh) | 建立车辆与地面耦合的一体式动力学模型的方法 | |
Wu et al. | Design and analysis of a sixteen-legged vehicle with reconfigurable close-chain leg mechanisms | |
Ames et al. | First steps toward formal controller synthesis for bipedal robots with experimental implementation | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
Gasoto et al. | A validated physical model for real-time simulation of soft robotic snakes | |
Hwang | Recursive Newton–Euler formulation for flexible dynamic manufacturing analysis of open-loop robotic systems | |
CN101635063A (zh) | 一种椭圆截面弹性柱体自由扭转/变形的模拟方法 | |
Ghotbi et al. | Vehicle-terrain interaction models for analysis and performance evaluation of wheeled rovers | |
Lee et al. | The development of a sliding joint for very flexible multibody dynamics using absolute nodal coordinate formulation | |
Abo-Shanab et al. | Tip-over stability of manipulator-like mobile hydraulic machines | |
Recuero et al. | ANCF continuum-based soil plasticity for wheeled vehicle off-road mobility | |
Pappalardo et al. | Design and development of a virtual model of an electric vehicle of category L7 | |
CN115310339A (zh) | 基于物质点法的具有表面张力效应的固液耦合模拟方法 | |
Csomay-Shanklin et al. | Nonlinear Model Predictive Control of a 3D Hopping Robot: Leveraging Lie Group Integrators for Dynamically Stable Behaviors | |
Jeong et al. | Analysis of mobile robot stability through 3D dynamics and lumped parameter tire modeling | |
Guo et al. | An improved 3D multi-sphere DE-FE contact algorithm for interactions between an off-road pneumatic tire and irregular gravel terrain | |
Roy | Experimental validation of non-cohesive soil using discrete element method | |
CN105740485A (zh) | 一种机械系统力学的运动仿真方法及系统 | |
Rubtsov et al. | Cooperative simulation process of unmanned track vehicle in RecurDyn and Simulink | |
Poulakis et al. | Port-based modeling and simulation of planetary rover locomotion on rough terrain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170606 Termination date: 20210625 |
|
CF01 | Termination of patent right due to non-payment of annual fee |