CN105447316B - 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 - Google Patents

一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 Download PDF

Info

Publication number
CN105447316B
CN105447316B CN201510860878.7A CN201510860878A CN105447316B CN 105447316 B CN105447316 B CN 105447316B CN 201510860878 A CN201510860878 A CN 201510860878A CN 105447316 B CN105447316 B CN 105447316B
Authority
CN
China
Prior art keywords
fatigue
unit
blade
composite material
formula
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
CN201510860878.7A
Other languages
English (en)
Other versions
CN105447316A (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.)
Zhengzhou Foamtech Nano Material Co Ltd
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201510860878.7A priority Critical patent/CN105447316B/zh
Publication of CN105447316A publication Critical patent/CN105447316A/zh
Application granted granted Critical
Publication of CN105447316B publication Critical patent/CN105447316B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种预测含初始缺陷复合材料桨叶结构疲劳寿命的方法,该方法有三大步骤:步骤一、建立含不同初始缺陷复合材料桨叶结构有限元模型;步骤二、含不同初始缺陷复合材料桨叶结构渐进损伤分析;步骤三、含不同初始缺陷复合材料桨叶结构剩余寿命累积损伤计算。本发明简单实用、需要的模型参数少、计算精度高,并能合理给出桨叶结构在疲劳载荷作用下的裂纹扩展路径等优点。

Description

一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法
技术领域
本发明提供一种预测含初始缺陷复合材料桨叶剩余寿命寿命的渐进损伤方法,属于复合材料结构疲劳可靠性技术领域。
背景技术
复合材料以质量轻、强度高等优点被广泛地应用于航空航天领域,尤其是在直升机结构的设计中,复合材料已占据了直升机的主要部分。而桨叶作为直升机的主要承力构件,在一些恶劣工作环境下使用时,很容易受到损伤,且其疲劳性能关系到整个直升机的飞行安全。现今对含初始缺陷复合材料结构疲劳寿命预测的方法主要是疲劳试验和有限元计算,试验方法需试样和试验成本高,且不方便设计使用;而一些有限元计算方法,需要通过大量的试验确定模型参数,导致预测成本高,计算复杂,误差大,计算精度和效率低;目前尚缺乏预测含初始缺陷复合材料桨叶结构疲劳寿命的简便方法,为此,本发明提出一种预测含初始缺陷复合材料桨叶结构疲劳寿命的方法,该方法简单实用、需要的模型参数少、计算精度高,并能合理给出桨叶结构在疲劳载荷作用下的裂纹扩展路径等优点,具有重要的学术意义和工程应用价值。
发明内容
1、目的:本发明的目的是提供一种预测含初始缺陷复合材料桨叶剩余寿命寿命的渐进损伤方法,该方法简单实用、需要的模型参数少、计算精度高,并能合理给出桨叶结构在疲劳载荷作用下的裂纹扩展路径等优点。
2、技术方案:本发明提供一种预测含初始缺陷复合材料桨叶剩余寿命寿命的渐进损伤方法,该方法具体步骤如下:
步骤一、建立含不同初始缺陷复合材料桨叶结构有限元模型
复合材料桨叶结构主要由碳纤维复合材料蒙皮和两根复合材料单向带大梁组成。为提高计算精度,有限元建模时,桨叶蒙皮和大梁分别采用solid46层单元和solid45三维实体单元进行模拟。初始缺陷包括中心圆孔、两侧缺口以及冲击损伤的三种不同缺陷类型。复合材料桨叶结构承受载荷以压强的形式施加到桨叶的上下翼面,复合材料桨叶结构边界条件为桨叶端面固支。
步骤二、含不同初始缺陷复合材料桨叶结构渐进损伤分析
含不同初始缺陷复合材料桨叶结构渐进损伤分析主要包括:(1)疲劳载荷下结构有限元应力分布计算(2)单元剩余强度失效判别准则(3)失效单元材料性能退化方式。具体为:
(1)疲劳载荷下结构有限元应力分布计算。采用步骤一中的三维有限元模型,可以准确得到桨叶结构的三维应力分布状态。由于桨叶属于柔性结构,使用过程中会有几何大变形产生,因此应力分析时需考虑几何大变形对应力分布的影响。
(2)单元剩余强度失效判别准则。复合材料结构受疲劳载荷加载时,其剩余强度是随着循环次数的增加而逐渐减小的。本文采用应力控制剩余强度的方程计算各单元的剩余强度值,具体表达式为
nelm=C(s-S0)m[R0-R(n)]b (1)
式(1)中nelm为单元疲劳寿命,s为疲劳应力,S0为拟合疲劳极限,R0为初始剩余强度值,R(n)为经过n次循环后材料的剩余强度,m、C和b为由疲劳试验拟合参数。应力控制剩余强度的方程可以由复合材料层合板疲劳试验得到。进行渐进损伤分析时,若单元应力小于其剩余强度,则认为单元失效。
(3)失效单元材料性能退化方式。若单元未发生失效,而疲劳载荷循环次数不断增加,则按照式(1)对未失效单元剩余强度进行退化;若单元发生失效,即单元各方向应力小于其剩余强度,则用生死单元技术杀死失效单元,即单元刚度降为0。
渐进损伤分析具体流程算法:输入桨叶单元类型、单元属性、几何尺寸、边界条件、疲劳应力,建立桨叶有限元模型;然后计算所有复合材料单元各主方向应力分布;使用剩余失效准则,对每个桨叶单元进行失效判定。如果桨叶单元无任何失效发生,则增加循环增量Ni,用疲劳剩余寿命模型式(1)计算该单元各主方向剩余强度,并采用材料性能逐渐退化方法进行材料剩余强度退化;如果桨叶单元有失效发生,则杀死失效单元,令其刚度将为0,并判定桨叶结构是否发生彻底破坏,若桨叶结构没有彻底破坏,则在此循环次数下进行单元应力分析,继续判定有无新的单元失效发生,直到该循环下桨叶结构无新的失效单元出现为止;增加循环增量Ni并按照式(1)对复合材料单元进行剩余强度退化,进行下一次循环;当判定桨叶结构发生彻底破坏,则将当前循环数作为桨叶结构疲劳寿命。
步骤三、含不同初始缺陷复合材料桨叶结构剩余寿命累计损伤计算
根据步骤一和二,计算不同恒幅疲劳载荷水平下的疲劳寿命。桨叶结构的载荷-寿命表达式如式(2)
式(2)中,F为疲劳载荷最大值,N为计算得到的桨叶疲劳寿命,m1和C1为拟合参数。式(2)中的拟合参数可以通过最小二乘法拟合得到。
由Miner累积损伤理论
式(3)中ni为疲劳载荷谱中第i级载荷的循环次数,Ni为桨叶结构在第i级疲劳载荷下发生失效的循环次数,j为疲劳载荷谱中载荷级数,T为载荷谱块数。
将式(2)代入式(3)可以得到复合材料桨叶结构谱载下的疲劳寿命为
式(4)中Fi表示疲劳载荷谱中第i级疲劳载荷值。
3、优点及功效:本发明提供一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法,其特点是简单实用、操作方便、计算精度高。
附图说明
图1为复合材料桨叶结构示意图。
图2为复合材料桨叶结构载荷边界条件示意图。
图3为本发明所述流程框图。
图中符号说明如下:
图3中的Ni为桨叶有限元模型中第i个单元的循环次数,F为疲劳载荷,N为复合材料桨叶结构不同恒幅疲劳载荷寿命。
具体实施方式
图3为本发明所述方法的流程框图,本发明分三大步骤实现,具体为:
步骤一、建立含不同初始缺陷复合材料桨叶结构有限元模型
复合材料桨叶结构主要由碳纤维复合材料蒙皮和两根复合材料单向带大梁组成。为提高计算精度,有限元建模时,桨叶蒙皮和大梁分别采用solid46层单元和solid45三维实体单元进行模拟。初始缺陷包括中心圆孔、两侧缺口以及冲击损伤的三种不同缺陷类型。复合材料桨叶结构承受载荷以压强的形式施加到桨叶的上下翼面,复合材料桨叶结构边界条件为桨叶端面固支。
步骤二、含不同初始缺陷复合材料桨叶结构渐进损伤分析
含不同初始缺陷复合材料桨叶结构渐进损伤分析主要包括:(1)疲劳载荷下结构有限元应力分布计算(2)单元剩余强度失效判别准则(3)失效单元材料性能退化方式。具体为:
(1)疲劳载荷下结构有限元应力分布计算。采用步骤一中的三维有限元模型,可以准确得到桨叶结构的三维应力分布状态。由于桨叶属于柔性结构,使用过程中会有几何大变形产生,因此应力分析时需考虑几何大变形对应力分布的影响。
(2)单元剩余强度失效判别准则。复合材料结构受疲劳载荷加载时,其剩余强度是随着循环次数的增加而逐渐减小的。本文采用应力控制剩余强度的方程计算各单元的剩余强度值,具体表达式为
nelm=C(s-S0)m[R0-R(n)]b (5)
式(5)中nelm为单元疲劳寿命,s为疲劳应力,S0为拟合疲劳极限,R0为初始剩余强度值,R(n)为经过n次循环后材料的剩余强度,m、C和b为由疲劳试验拟合参数。应力控制剩余强度的方程可以由复合材料层合板疲劳试验得到。进行渐进损伤分析时,若单元应力小于其剩余强度,则认为单元失效。
(3)失效单元材料性能退化方式。若单元未发生失效,而疲劳载荷循环次数不断增加,则按照式(5)对未失效单元剩余强度进行退化;若单元发生失效,即单元各方向应力小于其剩余强度,则用生死单元技术杀死失效单元,即单元刚度降为0。
渐进损伤分析具体流程算法:输入桨叶单元类型、单元属性、几何尺寸、边界条件、疲劳应力,建立桨叶有限元模型;然后计算所有复合材料单元各主方向应力分布;使用剩余失效准则,对每个桨叶单元进行失效判定。如果桨叶单元无任何失效发生,则增加循环增量Ni,用疲劳剩余寿命模型式(5)计算该单元各主方向剩余强度,并采用材料性能逐渐退化方法进行材料剩余强度退化;如果桨叶单元有失效发生,则杀死失效单元,令其刚度将为0,并判定桨叶结构是否发生彻底破坏,若桨叶结构没有彻底破坏,则在此循环次数下进行单元应力分析,继续判定有无新的单元失效发生,直到该循环下桨叶结构无新的失效单元出现为止;增加循环增量Ni并按照式(5)对复合材料单元进行剩余强度退化,进行下一次循环;当判定桨叶结构发生彻底破坏,则将当前循环数作为桨叶结构疲劳寿命。
步骤三、含不同初始缺陷复合材料桨叶结构剩余寿命累计损伤计算
根据步骤一和二,计算不同恒幅疲劳载荷水平下的疲劳寿命。桨叶结构的载荷-寿命表达式如式(6)
式(6)中,F为疲劳载荷最大值,N为计算得到的桨叶疲劳寿命,m1和C1为拟合参数。式(6)中的拟合参数可以通过最小二乘法拟合得到。
由Miner累积损伤理论
式(7)中ni为疲劳载荷谱中第i级载荷的循环次数,Ni为桨叶结构在第i级疲劳载荷下发生失效的循环次数,j为疲劳载荷谱中载荷级数,T为载荷谱块数。
将式(6)代入式(7)可以得到复合材料桨叶结构谱载下的疲劳寿命为
式(8)中Fi表示疲劳载荷谱中第i级疲劳载荷值。

Claims (1)

1.一种预测含初始缺陷复合材料桨叶结构疲劳寿命的方法,其特征在于:该方法具体步骤如下:
步骤一、建立含不同初始缺陷复合材料桨叶结构有限元模型
复合材料桨叶结构主要由碳纤维复合材料蒙皮和两根复合材料单向带大梁组成;为提高计算精度,有限元建模时,桨叶蒙皮和大梁分别采用solid46层单元和solid45三维实体单元进行模拟;初始缺陷包括中心圆孔、两侧缺口以及冲击损伤的三种不同缺陷类型;复合材料桨叶结构承受载荷以压强的形式施加到桨叶的上下翼面,复合材料桨叶结构边界条件为桨叶端面固支;
步骤二、含不同初始缺陷复合材料桨叶结构渐进损伤分析
含不同初始缺陷复合材料桨叶结构渐进损伤分析主要包括:(1)疲劳载荷下结构有限元应力分布计算(2)单元剩余强度失效判别准则(3)失效单元材料性能退化方式;具体为:
(1)疲劳载荷下结构有限元应力分布计算;采用步骤一中的三维有限元模型,可以准确得到桨叶结构的三维应力分布状态;由于桨叶属于柔性结构,使用过程中会有几何大变形产生,因此应力分析时需考虑几何大变形对应力分布的影响;
(2)单元剩余强度失效判别准则;复合材料结构受疲劳载荷加载时,其剩余强度是随着循环次数的增加而逐渐减小的;采用应力控制剩余强度的方程计算各单元的剩余强度值,具体表达式为
nelm=C(s-S0)m[R0-R(n)]b (1)
式(1)中nelm为单元疲劳寿命,s为疲劳应力,S0为拟合疲劳极限,R0为初始剩余强度值,R(n)为经过n次循环后材料的剩余强度,m、C和b为由疲劳试验拟合参数;应力控制剩余强度的方程可以由复合材料层合板疲劳试验得到;进行渐进损伤分析时,若单元应力小于其剩余强度,则认为单元失效;
(3)失效单元材料性能退化方式;若单元未发生失效,而疲劳载荷循环次数不断增加,则按照式(1)对未失效单元剩余强度进行退化;若单元发生失效,即单元各方向应力小于其剩余强度,则用生死单元技术杀死失效单元,即单元刚度降为0;
渐进损伤分析具体流程算法:输入桨叶单元类型、单元属性、几何尺寸、边界条件、疲劳应力,建立桨叶有限元模型;然后计算所有复合材料单元各主方向应力分布;使用剩余失效准则,对每个桨叶单元进行失效判定;如果桨叶单元无任何失效发生,则增加循环增量Ni,用疲劳剩余寿命模型式(1)计算该单元各主方向剩余强度,并采用材料性能逐渐退化方法进行材料剩余强度退化;如果桨叶单元有失效发生,则杀死失效单元,令其刚度将为0,并判定桨叶结构是否发生彻底破坏,若桨叶结构没有彻底破坏,则在此循环次数下进行单元应力分析,继续判定有无新的单元失效发生,直到该循环下桨叶结构无新的失效单元出现为止;增加循环增量Ni并按照式(1)对复合材料单元进行剩余强度退化,进行下一次循环;当判定桨叶结构发生彻底破坏,则将当前循环数作为桨叶疲劳寿命;
步骤三、含不同初始缺陷复合材料桨叶结构剩余寿命累计损伤计算
根据步骤一和二,计算不同恒幅疲劳载荷水平下的疲劳寿命;桨叶结构的载荷-寿命表达式如式(2)
式(2)中,F为疲劳载荷最大值,N为计算得到的桨叶疲劳寿命,m1和C1为拟合参数;式(2)中的拟合参数可以通过最小二乘法拟合得到;
由Miner累积损伤理论
式(3)中ni为疲劳载荷谱中第i级载荷的循环次数,Ni为桨叶结构在第i级疲劳载荷下发生失效的循环次数,j为疲劳载荷谱中载荷级数,T为载荷谱块数;
将式(2)代入式(3)可以得到复合材料桨叶结构谱载下的疲劳寿命为
式(4)中Fi表示疲劳载荷谱中第i级疲劳载荷值。
CN201510860878.7A 2015-12-01 2015-12-01 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 Active CN105447316B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510860878.7A CN105447316B (zh) 2015-12-01 2015-12-01 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510860878.7A CN105447316B (zh) 2015-12-01 2015-12-01 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法

Publications (2)

Publication Number Publication Date
CN105447316A CN105447316A (zh) 2016-03-30
CN105447316B true CN105447316B (zh) 2018-02-13

Family

ID=55557486

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510860878.7A Active CN105447316B (zh) 2015-12-01 2015-12-01 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法

Country Status (1)

Country Link
CN (1) CN105447316B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404649B (zh) * 2016-09-06 2020-11-03 中国航空工业集团公司沈阳飞机设计研究所 一种橡胶材料产品寿命预测方法
CN106650119B (zh) * 2016-12-27 2020-08-14 沈阳航空航天大学 Cfrp与钛合金叠层结构钻孔刀具寿命的预测方法
CN106768980B (zh) * 2017-01-19 2018-10-12 上海理工大学 零件无限寿命设计下的可回收性质量评价方法
CN107884290B (zh) * 2017-10-16 2019-11-19 北京理工大学 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法
CN108681641B (zh) * 2018-05-18 2021-01-08 北京航空航天大学 一种考虑缺陷尺寸效应的复合材料缺陷容限性能表征方法
CN108693054B (zh) * 2018-05-18 2020-05-12 北京航空航天大学 一种复合材料结构谱载疲劳寿命的渐进损伤算法
CN110929316A (zh) * 2019-09-16 2020-03-27 江苏卓然智能重工有限公司 一种大型塔吊缆索疲劳评估中的可靠性分析模型
CN110987388B (zh) * 2019-11-27 2021-06-22 南京航空航天大学 一种基于缺口疲劳强度的等效机械加工缺口的方法
CN111046547B (zh) * 2019-12-04 2022-09-16 中国直升机设计研究所 一种直升机矩形复合材料柔性梁高周寿命快速计算方法
CN110884684B (zh) * 2019-12-04 2022-04-12 中国直升机设计研究所 一种直升机无轴承尾桨叶弹击后的强度试验设计方法
CN111597716B (zh) * 2020-05-19 2022-04-29 北京航空航天大学 一种含分层损伤复合材料层合板的疲劳寿命预测方法
CN112710539B (zh) * 2020-12-22 2024-05-17 沈阳工业大学 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法
CN113449375B (zh) * 2021-04-20 2023-03-24 中国直升机设计研究所 一种复合材料桨叶疲劳寿命的半解析计算方法
CN117672436A (zh) * 2024-01-31 2024-03-08 湖南大学 一种考虑制造工艺缺陷的复合材料结构性能预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103698199A (zh) * 2013-09-26 2014-04-02 北京航空航天大学 一种基于细观力学退化模型的复合材料结构失效分析方法
CN103927422A (zh) * 2014-04-24 2014-07-16 北京航空航天大学 一种基于渐进损伤模型的预测复合材料螺栓连接失效的强度包线法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103698199A (zh) * 2013-09-26 2014-04-02 北京航空航天大学 一种基于细观力学退化模型的复合材料结构失效分析方法
CN103927422A (zh) * 2014-04-24 2014-07-16 北京航空航天大学 一种基于渐进损伤模型的预测复合材料螺栓连接失效的强度包线法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Effects of matrix ductility and progressive damage on fatigue strengths of unnotched and notched carbon fibre plain woven roving fabric laminates;M.Kawai et al;《Composites Part A: Applied Science and Manufacturing》;19961231;第27卷(第6期);493-502页 *
Fatigue damage behaviors of carbon fiber-reinforced epoxy composites containing nanoclay;Shafi Ullah Khan et al;《Composites Science and Technology》;20101130;第70卷(第14期);2077-2085页 *
Simulation of fatigue failure in a full composite wind turbine blade;Mahmood M.Shokrieh et al;《Composite Structures》;20060831;第74卷(第3期);332-342页 *
复合材料疲劳可靠性研究;熊峻江等;《材料工程》;20001231(第5期);40-44页 *
稳态循环载荷下疲劳/断裂可靠性寿命估算;熊峻江等;《应用力学学报》;19970930;第14卷(第3期);14-18页 *

Also Published As

Publication number Publication date
CN105447316A (zh) 2016-03-30

Similar Documents

Publication Publication Date Title
CN105447316B (zh) 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法
CN103983467B (zh) 一种基于服役状态的单机服役使用寿命监控方法
JP6023882B2 (ja) 機械部品の寿命消費の信頼できる予測
CN107506529B (zh) 一种复合材料加筋壁板轴压稳定性计算方法
CN112710539B (zh) 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法
CN110083916A (zh) 基于自构隶属函数的机械结构模糊疲劳可靠度优化方法
CN104820892A (zh) 一种基于数据传递的航空发电系统定量危害性分析方法
CN104007351B (zh) 一种航天器电子组件热循环试验方案确定方法
CN104462700A (zh) 基于故障物理的电子产品可靠性仿真试验方法
CN105335582A (zh) 一种飞机复材壁板重量分析建模方法
CN108693054B (zh) 一种复合材料结构谱载疲劳寿命的渐进损伤算法
Bacharoudis et al. A probabilistic approach for strength and stability evaluation of wind turbine rotor blades in ultimate loading
Chen et al. A critical review of damage and failure of composite wind turbine blade structures
US20150227658A1 (en) Generating a simplified calculation model and predicting life consumption of a component
CN103324800A (zh) 悬臂式掘进机截割头计算机辅助设计方法及cad系统
CN105138850B (zh) 一种民用飞机系统维修时间间隔计算方法
CN114169144A (zh) 一种力热振复合环境下随机载荷谱编制方法
CN104123458A (zh) 一种基于应变能理论的横断型斜裂纹转子变刚度特性计算方法
CN105334032B (zh) 一种机型结构寿命的计算方法
Muyan et al. Strength analysis of a 5-m composite wind turbine blade under static and fatigue loading conditions
CN104361193A (zh) 一种产品寿命主机理确定方法
Lekou Probabilistic design of wind turbine blades
CN111339627B (zh) 计算流体动力学分析异常症候预测系统及方法
CN104331788B (zh) 发电系统可靠性评估的严重故障集快速筛选方法
CN112115626A (zh) 一种风电机组叶片疲劳寿命预测方法和装置

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210325

Address after: 100191 building 66-111, 14 Huayuan North Road, Haidian District, Beijing

Patentee after: Beijing mubangren Technology Co.,Ltd.

Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road

Patentee before: BEIHANG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240201

Address after: 451100 north side of Yanhe Road, Hezhuang Town, Xinzheng City, Zhengzhou City, Henan Province

Patentee after: ZHENGZHOU FOAMTECH NANO MATERIAL Co.,Ltd.

Country or region after: China

Address before: 100191 building 66-111, 14 Huayuan North Road, Haidian District, Beijing

Patentee before: Beijing mubangren Technology Co.,Ltd.

Country or region before: China