CN110222463B - 基于vcct确定不同裂纹长度的能量释放率的优化方法 - Google Patents
基于vcct确定不同裂纹长度的能量释放率的优化方法 Download PDFInfo
- Publication number
- CN110222463B CN110222463B CN201910530776.7A CN201910530776A CN110222463B CN 110222463 B CN110222463 B CN 110222463B CN 201910530776 A CN201910530776 A CN 201910530776A CN 110222463 B CN110222463 B CN 110222463B
- Authority
- CN
- China
- Prior art keywords
- crack
- matrix
- node
- energy release
- len
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于VCCT确定不同裂纹长度的能量释放率的优化方法,包括以下步骤:S1:在ABAQUS中建立复合式路面结构模型;在复合式路面模型上预设裂纹和裂纹参数,然后将所述复合式路面模型进行网格划分,其中,裂纹参数包括裂纹扩展方向、时间增量步和牵引开裂值;S2:从ABAQUS中提取复合式路面结构模型的裂纹扩展面的粘结节点并将所述粘结节点按照裂纹扩展路径重新排列,形成裂纹扩展节点矩阵;S3:根据所述裂纹扩展节点矩阵,在ABAQUS中计算不同裂纹长度对应的能量释放率。本发明实现在ABAQUS的单个模型中对不同裂纹长度对应的能量释放率的计算,具有建模工作量小、计算量小和精度高的有益技术效果。
Description
技术领域
本发明属于道路工程数据处理技术领域,尤其涉及一种基于VCCT确定不同裂纹长度的能量释放率的优化方法。
背景技术
虚拟裂纹扩展分析技术(VCCT)是基于裂纹扩展时释放的能量等于将裂纹闭合一个扩展增量所需的功的假设的一种断裂力学分析方法,该方法假设虚拟裂纹尖端后面的展开位移和实际裂纹尖端后面的张开位移近似相等,通过虚拟裂纹扩展分析技术可获得应力或应变的精确解,但基于虚拟裂纹扩展分析技术现有的ABAQUS中的一个模型只能获得单一裂纹长度对应的能量释放率,无法在单个模型中获得不用裂纹长度对应的能量释放率,若要获得不同裂纹长度对应的能量释放率只能对不同裂纹长度进行建模,如计算n中裂纹长度对应的能量释放率,需要建立n个模型,这样大大增大了建模工作量和计算量,且对计算机的配置要求高。
因此,亟需一种基于虚拟裂纹扩展分析技术确定不同裂纹长度的能量释放率的优化方法。
发明内容
有鉴于此,本发明提供一种基于VCCT确定不同裂纹长度的能量释放率的优化方法以解决现有技术的不足。
本发明提供一种基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:包括以下步骤:
S1:在ABAQUS中建立复合式路面结构模型;在复合式路面模型上预设裂纹和裂纹参数,然后将所述复合式路面模型进行网格划分,其中,所述裂纹参数包括裂纹扩展方向、时间增量步和牵引开裂值;
S2:从ABAQUS中提取所述裂纹扩展面的粘结节点并将所述粘结节点按照裂纹扩展路径重新排列,形成裂纹扩展节点矩阵;
S3:根据所述裂纹扩展节点矩阵,在ABAQUS中计算不同裂纹长度对应的能量释放率。
进一步,所述步骤S2包括以下步骤:
S21:将位于预设裂纹两侧的粘结节点分别放置于L矩阵和R矩阵;所述粘结节点至少包括节点坐标和节点编号;
S22:将所述L矩阵按照列进行升序排列形成新的矩阵L1;将所述R矩阵按照列进行升序排列形成新的矩阵R1;
S23:将矩阵L1中的第三列赋值给矩阵NL,将矩阵R1的第三列赋值给矩阵NR;即将L1矩阵的节点编号赋值给矩阵NL,将R1矩阵的节点编号赋值给矩阵NR,即矩阵NL和矩阵NR为节点编号;
S24:定义列向量N1和N2:
S241:判断W与len的大小,其中,W为中间变量,其初始值为1,len表示矩阵NL中的元素个数或NR中的元素个数,若W小于等于len,进一步判断W除row的余数是否是1,其中,row=len/col,col表示节点列数,若余数不是1,进入步骤S242,若余数是1,则进入步骤S243,若W大于len,则退出;
S242:N1[K]=NL[W],N2[K]=NR[W],其中,K为列向量N1和N2的行数,W表示阵列NL和阵列NR的行数;
S243:W=W+1,K=K+1;
S244:重复步骤S241至步骤S243,直至W大于len退出;
S25:定义列向量N3和N4;
S251:判断j与len的大小,其中,j为中间变量,其初始值为1,len表示矩阵NL中的元素个数或NR中的元素个数,若j小于等于len,进一步判断j除row的余数是否是0,其中,row=len/col,col表示节点列数,若余数不是0,进入步骤S252,若余数是0,则进入步骤S253,若j大于len,则退出;
S252:N3[n]=NL[j],N2[n]=NR[j],其中,n为列向量N3和N4的行数,j表示阵列NL和阵列NR的行数;
S253:j=j+1,n=n+1;
S254:重复步骤S251至步骤S253,直至W大于len退出;
S26:输出裂纹扩展节点矩阵NODE,NODE=[N1 N2 N3 N4];
S27:将裂纹扩展节点矩阵NODE作为自定义单元写入AQAQUS的inp文件中,形成修改后的inp文件。
进一步,所述裂纹扩展节点矩阵的一行的四个节点表示一个哑结点断裂单元的四个节点。
进一步,其中所述L矩阵和R矩阵均为多行三列矩阵,从左至右依次为,第一列为节点的x轴坐标,第二列为节点的y轴坐标,第三列为节点的编号。
进一步,所述能量释放率包括Ⅰ型能量释放率和Ⅱ型能量释放率。
进一步,所述不同裂纹长度对应的能量释放率在ABAQUS中采用如下方法计算:
S31:从inp文件读入自定义单元的第i行的节点,即裂纹扩展节点矩阵NODE的第i行的节点,其中,i为中间变量,其初始值为1,形成弹性刚度矩阵AMATRX,并更新右手端矢量RHS;
S32:计算哑结点断裂单元的能量释放率,并赋值给状态变量矩阵SVARS,其中,所述状态变量矩阵SVARS包括SVARS(1)和SVARS(2);
S33:判断分析步时间是否小于等于分析步一,即Step-1的结束时间,若是,则弹性刚度矩阵AMATRX元素数值不变,裂纹不发生扩展,进入步骤S34;若否,进一步判断状态变量SVARS的元素SVARS(2)或SVARS(1)是否大于预设的牵引开裂临界值,若是,则弹性刚度矩阵AMATRX元素全部赋值为0,裂纹扩展,裂纹长度增长△a,进入步骤S34;若否,弹性刚度矩阵AMATRX元素数值不变,裂纹不发生扩展,进入步骤S34;
S34:判断当前分析步时间是否大于Step-2结束时间,其中,step-2结束时间表示裂纹扩展的计算总时间,若否,则分析步时间增加一个时间增量步,进入步骤S31;若是,退出。
进一步,所述步骤S33还包括将所述状态变量矩阵SVARS赋值给ABAQUS状态变量SDV1和SDV2,并输出SDV1和SDV2。
本发明的有益效果:本发明将裂纹扩展面粘结节点按预设裂纹扩展路径进行排序,以及裂纹扩展的条件判断和每个裂纹扩展增量的能量释放率的计算,实现在ABAQUS的单个模型中对不同裂纹长度对应的能量释放率的计算,具有建模工作量小、计算量小和精度高的有益技术效果。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1为生成自定义单元Matlab程序流程图。
图2为计算不同裂纹长度对应的能量释放率的流程图。
具体实施方式
以下结合说明书附图对本发明做出进一步的说明:
本发明提供一种基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:包括以下步骤:
S1:在ABAQUS中建立复合式路面结构模型;在复合式路面模型上预设裂纹和裂纹参数,然后将所述复合式路面模型进行网格划分,其中,所述裂纹参数包括裂纹扩展方向、时间增量步和牵引开裂值;
S2:从ABAQUS中提取所述裂纹扩展面的粘结节点并将所述粘结节点按照裂纹扩展路径重新排列,形成裂纹扩展节点矩阵;
S3:根据所述裂纹扩展节点矩阵,在ABAQUS中计算不同裂纹长度对应的能量释放率。
其中,VCCT表示Virtual cracK closure technique的简称,即虚拟裂纹扩展分析技术或虚拟裂纹闭合法的简称。
通过上述技术方案,对裂纹扩展面粘结节点按预设裂纹扩展路径进行排序,以及裂纹扩展的条件判断和每个裂纹扩展增量的能量释放率的计算,实现在ABAQUS的单个模型中对不同裂纹长度对应的能量释放率的计算,具有建模工作量小、计算量小和精度高的有益技术效果。解决了现有ABAQUS中,计算不同裂纹长度对用个的能量释放率需要重复建立多个模型,以及计算量大的问题,实现了在单个模型中计算不同裂纹长度对应的能量释放率,即实现了在ABAQUS中模拟预设裂纹的连续扩展。
在本实施例中,所述步骤S2包括以下步骤:
S21:将位于预设裂纹两侧的粘结节点分别放置于L矩阵和R矩阵;所述粘结节点至少包括节点坐标和节点编号;所述节点坐标可以是平面二维坐标也可以是立体的三维坐标,所述节点编号是ABAQUS程序对模型内的节点的编号;其中所述L矩阵和R矩阵均为多行三列矩阵,其中所述L矩阵和R矩阵均为多行三列矩阵,从左至右依次为,第一列为节点的x轴坐标,第二列为节点的y轴坐标,第三列为节点的编号。
S22:将所述L矩阵按照列进行升序排列形成新的矩阵L1;将所述R矩阵按照列进行升序排列形成新的矩阵R1;
S23:将矩阵L1中的第三列赋值给矩阵NL,将矩阵R1的第三列赋值给矩阵NR;
S24:定义列向量N1和N2:
S241:判断W与len的大小,其中,W为中间变量,其初始值为1,len表示矩阵NL中的元素个数或NR中的元素个数,若W小于等于len,进一步判断W除row的余数是否是1,其中,row=len/col,col表示节点列数,若余数是1,进入步骤S242,若余数不是1,则进入步骤S243,若W大于len,则退出;
S242:N1[K]=NL[W],N2[K]=NR[W],其中,K为列向量N1和N2的行数,W表示阵列NL和阵列NR的行数;
S243:W=W+1,K=K+1;
S244:重复步骤S241至步骤S243,直至W大于len退出;
S25:定义列向量N3和N4;
S251:判断j与len的大小,其中,j为中间变量,其初始值为1,len表示矩阵NL中的元素个数或NR中的元素个数,若j小于等于len,进一步判断j除row的余数是否是0,其中,row=len/col,col表示节点列数,若余数是0,进入步骤S252,若余数不是0,则进入步骤S253,若j大于len,则退出;
S252:N3[n]=NL[j],N4[n]=NR[j],其中n为列向量N3和N4的行数,j表示阵列NL和阵列NR的行数;
S253:j=j+1,n=n+1;
S254:重复步骤S251至步骤S253,直至W大于len退出;
S26:输出裂纹扩展节点矩阵NODE,NODE=[N1 N2 N3 N4];
S27:将裂纹扩展节点矩阵NODE作为自定义单元写入AQAQUS的inp文件中,形成修改后的inp文件。
其中,本步骤的计算本领域技术人员可根据实际选择实现方式,如手动计算、或通过计算机实现,在本实施例中,通过MATLAB软件来实现。
通过上述技术方案,将裂纹面的粘结节点从复合式道路模型中读取出来,并根据预设的裂纹扩展路径重新排序,使重新排序后矩阵的符合裂纹尖端的哑结点断裂单元的节点排序。
在本实施例中,所述裂纹扩展节点矩阵的一行的四个节点表示一个哑结点断裂单元的四个节点。通过上述技术方案可为后续裂纹的连续扩展提供哑结点断裂单元的扩展路径。
在本实施例中,所述能量释放率包括Ⅰ型能量释放率和Ⅱ型能量释放率。所述Ⅰ型能量释放率表示垂直方向的拉伸裂纹驱动力,所述Ⅱ型能量释放率表示水平方向的剪切裂纹驱动力。
在本实施例中,所述不同裂纹长度对应的能量释放率在ABAQUS中采用如下方法计算:
S31:从inp文件读入自定义单元的第i行的节点,即裂纹扩展节点矩阵NODE的第i行的节点,其中,i为中间变量,其初始值为1,形成弹性刚度矩阵AMATRX,并更新右手端矢量RHS;
S32:计算哑结点断裂单元的能量释放率,并赋值给状态变量矩阵SVARS,其中,所述状态变量矩阵SVARS包括SVARS(1)和SVARS(2);
所述Ⅰ型能量释放率GⅠ和Ⅱ型能量释放率GⅡ采用如下方法计算:
其中,Fx、Fy分别为节点x轴和y轴方向的节点力,ΔV1、ΔU1分别为相邻节点编号x轴和y轴方向的位移差,Δa为裂纹扩展长度,B表示裂纹厚度;
S33:判断时间增量步是否小于等于分析步一,即Step-1的结束时间,此步骤即是判断负载是否加载完成,若是,即负载尚未加载完成,则弹性刚度矩阵AMATRX元素数值不变,裂纹不扩展,进入步骤S34;若否,进一步判断状态变量SVARS的元素SVARS(2)或SVARS(1)是否大于预设的牵引开裂临界值,若是,则弹性刚度矩阵AMATRX元素全部赋值为0,即弹性刚度矩阵之间以没有牵引力,裂纹扩展,裂纹长度增长△a,进入步骤S34;若否,弹性刚度矩阵AMATRX元素数值不变,裂纹不发生扩展,进入步骤S34;
S34:判断当前分析步时间是否大于Step-2结束时间,其中,step-2结束时间表示裂纹扩展的计算总时间,若否,则分析步时间增加一个时间增量步,进入步骤S31;若是,退出。
进一步,所述步骤S33还包括将所述状态变量矩阵SVARS赋值给ABAQUS状态变量SDV1和SDV2,并输出SDV1和SDV2。通过上述技术方案,将机器可读但用户无法获得的能量释放率转换为机器和用户均可识别的能量释放率。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:包括以下步骤:
S1:在ABAQUS中建立复合式路面结构模型;在复合式路面模型上预设裂纹和裂纹参数,然后将所述复合式路面模型进行网格划分,其中,所述裂纹参数包括裂纹扩展方向、时间增量步和牵引开裂值;
S2:从ABAQUS中提取复合式路面结构模型的裂纹扩展面的粘结节点并将所述粘结节点按照裂纹扩展路径重新排列,形成裂纹扩展节点矩阵;
所述步骤S2包括以下步骤:
S21:将位于预设裂纹两侧的粘结节点分别放置于L矩阵和R矩阵;所述粘结节点至少包括节点坐标和节点编号;
S22:将所述L矩阵按照列进行升序排列形成新的矩阵L1;将所述R矩阵按照列进行升序排列形成新的矩阵R1;
S23:将矩阵L1中的第三列赋值给矩阵NL,将矩阵R1的第三列赋值给矩阵NR;
S24:定义列向量N1和N2:
S241:判断W与len的大小,其中,W为中间变量,其初始值为1,len表示矩阵NL中的元素个数,若W小于等于len,进一步判断W除row的余数是否是1,其中,row=len/col,col表示节点列数,若余数是1,则进入步骤S243,若余数不是1,则进入步骤S242,若W大于len,则退出;
S242:N1[K] =NL[W], N2[K] =NR[W],其中,K为列向量N1和N2的行数,W表示阵列NL和阵列NR的行数;
S243:W=W+1,K=K+1;
S244:重复步骤S241至步骤S243,直至W大于len退出;
S25:定义列向量N3和N4;
S251:判断j与len的大小,其中,j为中间变量,其初始值为1,len表示矩阵NL中的元素个数或NR中的元素个数,若j小于等于len,进一步判断j除row的余数是否是0,其中,row=len/col,col表示节点列数,若余数是0,进入步骤S253,若余数不是0,则进入步骤S252,若j大于len,则退出;
S252:N3[n] =NL[j], N4[n] =NR[j],其中,n为列向量N3和N4的行数,j表示阵列NL和阵列NR的行数;
S253:j=j+1,n=n+1;
S254:重复步骤S251至步骤S253,直至j大于len退出;
S26:输出裂纹扩展节点矩阵NODE,NODE=[N1 N2 N3 N4];
S27:将裂纹扩展节点矩阵NODE作为自定义单元写入AB AQUS的inp文件中,形成修改后的inp文件;
S3:根据所述裂纹扩展节点矩阵,在ABAQUS中计算不同裂纹长度对应的能量释放率;所述不同裂纹长度对应的能量释放率在ABAQUS中采用如下方法计算:
S31:从inp文件读入自定义单元的第i行的节点,即裂纹扩展节点矩阵NODE的第i行的节点,其中,i为中间变量,其初始值为1,形成弹性刚度矩阵AMATRX,并更新右手端矢量RHS;
S32:计算哑结点断裂单元的能量释放率,并赋值给状态变量矩阵SVARS,其中,所述状态变量矩阵SVARS包括SVARS(1)和SVARS(2);
S33:判断分析步时间是否小于等于分析步一,即Step-1的结束时间,若是,则弹性刚度矩阵AMATRX元素数值不变,裂纹不发生扩展,进入步骤S34;若否,进一步判断状态变量SVARS的元素SVARS(2)或SVARS(1)是否大于预设的牵引开裂临界值,若是,则弹性刚度矩阵AMATRX元素全部赋值为0,裂纹扩展,裂纹长度增长△a,进入步骤S34;若否,弹性刚度矩阵AMATRX元素数值不变,裂纹不发生扩展,进入步骤S34;
S34:判断当前分析步时间是否大于Step-2结束时间,其中,step-2结束时间表示裂纹扩展的计算总时间,若否,则分析步时间增加一个时间增量步,进入步骤S31;若是,退出。
2.根据权利要求1所述基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:所述裂纹扩展节点矩阵的一行的四个节点表示一个哑结点断裂单元的四个节点。
3.根据权利要求1所述基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:其中所述L矩阵和R矩阵均为多行三列矩阵,从左至右依次为,第一列为节点的x轴坐标,第二列为节点的y轴坐标,第三列为节点的编号。
4.根据权利要求1所述基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:所述能量释放率包括Ⅰ型能量释放率和Ⅱ型能量释放率。
5.根据权利要求1所述基于VCCT确定不同裂纹长度的能量释放率的优化方法,其特征在于:所述步骤S33还包括将所述状态变量矩阵SVARS赋值给ABAQUS状态变量SDV1和SDV2,并输出SDV1和SDV2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910530776.7A CN110222463B (zh) | 2019-06-19 | 2019-06-19 | 基于vcct确定不同裂纹长度的能量释放率的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910530776.7A CN110222463B (zh) | 2019-06-19 | 2019-06-19 | 基于vcct确定不同裂纹长度的能量释放率的优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110222463A CN110222463A (zh) | 2019-09-10 |
CN110222463B true CN110222463B (zh) | 2023-01-24 |
Family
ID=67817811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910530776.7A Active CN110222463B (zh) | 2019-06-19 | 2019-06-19 | 基于vcct确定不同裂纹长度的能量释放率的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110222463B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116030920B (zh) * | 2023-01-17 | 2023-07-21 | 中国人民解放军陆军工程大学 | 一种裂纹尖端应变能释放率分析方法、装置、电子设备及存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844941A (zh) * | 2017-01-18 | 2017-06-13 | 西安建筑科技大学 | 一种基于虚拟裂纹闭合法获取模拟声发射信号的方法 |
-
2019
- 2019-06-19 CN CN201910530776.7A patent/CN110222463B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110222463A (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109165404B (zh) | 一种扩展光滑无网格伽辽金法 | |
CN110188423A (zh) | 一种基于有限元网格划分的线性工程结构快速bim建模方法 | |
CN110222463B (zh) | 基于vcct确定不同裂纹长度的能量释放率的优化方法 | |
CN105354388A (zh) | 一种加强筋的生长式拓扑优化设计方法 | |
CN111177861B (zh) | 适用于增材制造成形技术的常平环结构轻量化设计方法 | |
CN112116128A (zh) | 一种用于冲击载荷作用下结构层裂与多重层裂模拟预测方法 | |
CN111125963A (zh) | 基于拉格朗日积分点有限元的数值仿真系统及方法 | |
CN112487673B (zh) | 基于机床工作状态的关键主机组件结构优化设计方法 | |
CN109063401B (zh) | 采用等效静态位移实现结构接触力优化的方法 | |
CN112942837B (zh) | 一种悬挑结构混凝土3d打印方法及系统 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Šešok et al. | Global optimization of trusses with a modified genetic algorithm | |
CN109271655A (zh) | 一种基于非对称有限元算法的材料尺度效应分析方法 | |
CN107545097A (zh) | 一种高效铝合金榫卯连接结构的cae模拟优化分析方法 | |
KR100911167B1 (ko) | 내부에 보강체가 포함된 대상체의 가상절점을 이용한 유한요소 모델링 방법 | |
CN108776743B (zh) | 用于增材制造中针对悬垂结构基于正态分布布置支撑位的生成方法 | |
CN102637223A (zh) | 基于伪密度排序并考虑拔模制造约束的拓扑优化设计方法 | |
CN107273558A (zh) | 用于三维模型的机械性能校验的方法和装置 | |
CN112685936B (zh) | 一种用于贝壳珍珠母微结构有限元分析的建模方法 | |
Zhao et al. | Simulation of the dynamic cracking of brittle materials using a nonlocal damage model with an effective strain rate effect | |
CN111259589B (zh) | 一种考虑破损-安全的连续体频率约束拓扑优化设计方法 | |
Lv et al. | Gyroid minimal surface-based composite porous structure with superior mechanical property via triangular grid design and stress distribution optimization | |
CN107832537B (zh) | 一种复杂形式残余应力的数值模拟方法 | |
De Filippo et al. | Concept of a pseudo-lower bound solution for reinforced concrete slabs | |
CN116030920B (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 |