CN104252566B - 一种箱体结构的简化及装夹变形仿真分析方法 - Google Patents
一种箱体结构的简化及装夹变形仿真分析方法 Download PDFInfo
- Publication number
- CN104252566B CN104252566B CN201410508770.7A CN201410508770A CN104252566B CN 104252566 B CN104252566 B CN 104252566B CN 201410508770 A CN201410508770 A CN 201410508770A CN 104252566 B CN104252566 B CN 104252566B
- Authority
- CN
- China
- Prior art keywords
- clamping
- mtd
- msub
- casing
- deformation
- 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
Abstract
本发明公开了一种箱体结构的简化及装夹变形仿真分析方法,该方法在发动机上曲轴箱体结合面的铣削加工中,合理的夹紧布局方案和夹紧力对于加工变形的控制具有重要的意义。首先对装夹的约束条件进行了探讨,其次分析了箱体加工过程的装夹工艺以及装夹载荷的计算,然后建立了发动机箱体铣削加工的有限元模型,获得箱体的装夹变形云图,最后分析了工件在不同夹紧力和布局下的加工变形,为零件加工精度的提高提供依据。
Description
技术领域
本发明属于机加工领域,尤其涉及一种箱体结构的简化及装夹变形仿真分析方法。
背景技术
薄壁零件刚性差,在加工中因受切削力、夹紧力以及切削热等因素,影响极易产生变形。因此控制加工变形是保证薄壁零件加工质量的关键。发动机箱体属于薄壁多孔类箱体零件。有生产实践证明,对于所有的复杂箱体和薄壁类零件,装夹变形都是造成加工后平面度、位置度达不到精度要求的重要原因之一,工件的加工误差20%-60%都是由装夹引起的。研究表明,装夹导致的加工误差来自三方面:a)薄壁零件刚性差,加工时装夹引起弹性变形影响表面尺寸精度和形状、位置精度;b)施加装夹过程中,装夹顺序不当或者装夹力施加不均匀,工件定位不准,定位后的工件和理想状态存在位置和角度误差;c)装夹力不足,工件在加工过程中会出现滑移和振动,导致加工误差乃至刀具和工件的损坏。缸体和缸盖是发动机的主要机加对象,在设备调试和加工过程中,由于夹具夹紧力的影响,加工尺寸会发生超差现象。
发明内容
本发明的目的在于提供一种箱体结构的简化及装夹变形仿真分析方法,旨在分析夹紧力和夹紧力布局对缸体加工精度的影响,从而为控制发动机缸体中夹具的布局和夹紧力的确定提供参考。
本发明是这样实现的,一种箱体结构的简化及装夹变形仿真分析方法包括:
步骤一、节点四面体二次单元,采用10节点四面体二次单元对发动机箱体进行有限元网格划分,10节点四面体二次单元即为在原4节点四面体单元的每一条棱边上再增加一个位于中点位置的内部节点,采用体积坐标应用Lagrange插值公式计算形状函数矩阵;
步骤二、有限元模型的建立及其结果,包括几何模型及网格划分、材料属性及接触对定义、约束定义及载荷施加;
步骤三、不同夹紧位置下装夹变形分析,包括不同夹紧位置下装夹变形分析、不同夹紧力下装夹变形分析。
进一步,有限元模型的建立及其结果,具体方法为:
第一步、几何模型及网格划分:
对发动机箱体进行有限元建模时,采用虚拟拓扑方法,对整体力学性能影响小的几何细节进行整合;
第二步、材料属性及接触对定义,将夹具和工件定义为弹性体以提高分析精度,由于箱体尺寸大,采用弹性无接触模型来进行有限元分析;
第三步、约束定义及载荷施加,限制箱体一端面的移动自由度来限制z向移动自由度,限制箱体左右两个缸面除z方向外的5个自由度;
刀具对箱体的加工过程是连续的,此处基于以下假设:刀具在某工况下其铣削力对箱体的作用是准静态,即可将连续铣削过程离散;刀具相对于铸铝合金箱体为刚性,铣削力在箱体上均匀分布;由于考虑铣削过程材料去除量相对于整个箱体体积极其小,不考虑加工过程中箱体刚度的变化;为计算方便,假设在铣削过程中,铣削力随着刀具的移动,只发生位置和方向的变化,其大小保持不变;
进一步、不同夹紧位置下装夹变形分析,具体方法为:
第一步、计算待夹紧面面积,按照夹紧力的计算方法,得到夹紧压强并保持不变,夹紧力位置分别为箱体左右两侧的缸孔上的不规则面分析不同工况下夹具对加工变形的产生的影响;
第二步、不同夹紧力下装夹变形分析,采用不同的三种装夹方案,各个方案下夹紧压强不同,铣削力保持不变,计算箱体的加工变形。
进一步,体积坐标的定义方式如下:对于给定的四面体,其体内任何一点G的体积坐标:
等价于
Li=hi/Hi
其中,Hi为顶点i到对面的距离;hi为G点到i点对应面的距离;体积Gjmp为四面体Gjmp的体积值;Δ为给定四面体的体积;
10节点四面体单元共有30个节点位移自由度(DOF),其节点位移列阵为:
由于该单元有4个角节点,在每个棱边上的中点有一个中间节点(共6个),总共10个节点,由函数构造的Pascal三角形,取单元位移场模式为
u(x,y,x)=a1+a2x+a3y+a4z+a5xy+a6yz+a7xz+a8x2+a9y2+a10z2
以自然坐标表示,那么单位位移场模式为:
u(x,y,z)=a′1L1+a′2L2+a′3L3+a′4L4+a′5L1L2+a′6L2L3
+a′7L3L4+a′8L1L3+a′9L1L4+a′10L2L4
=N1u1+N2u2+N3u3+…+N10u10
对于角节点,相应的形状函数为:
Ni=(2Li-1)Li(i=1,2,3,4)
对于棱边上的中间节点,相应的形状函数为
N5=4L1L2,N6=4LaL2,N7=4L1L3
N8=4L1L4,N9=4L2L4,N10=4L3L4
得到形状函数后,按照有限元分析的标准过程即可推得相应的几何矩阵、刚度矩阵、节点等效载荷矩阵以及刚度方程。
效果汇总
本发明在发动机上曲轴箱体结合面的铣削加工中,合理的夹紧布局方案和夹紧力对于加工变形的控制具有重要的意义。首先对装夹的约束条件进行了探讨,其次分析了箱体加工过程的装夹工艺以及装夹载荷的计算,然后建立了发动机箱体铣削加工的有限元模型,获得箱体的装夹变形云图,最后分析了工件在不同夹紧力和布局下的加工变形,为零件加工精度的提高提供依据。
附图说明
图1是本发明实施例提供的节点四面体二次单元的结构示意图;
图2是本发明实施例提供的箱体结构的简化及装夹变形仿真分析方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明是这样实现的,一种箱体结构的简化及装夹变形仿真分析方法包括:
步骤一、节点四面体二次单元,采用10节点四面体二次单元对发动机箱体进行有限元网格划分,10节点四面体二次单元即为在原4节点四面体单元的每一条棱边上再增加一个位于中点位置的内部节点,如图1所示,采用体积坐标应用Lagrange插值公式计算形状函数矩阵,图中4个角节点上的坐标是基于体积坐标下的坐标值;
体积坐标的定义方式如下:对于给定的四面体,其体内任何一点G的体积坐标:
等价于
Li=hi/Hi
其中,Hi为顶点i到对面的距离;hi为G点到i点对应面的距离;体积qmp为四面体qmp的体积值;△为给定四面体的体积;
10节点四面体单元共有30个节点位移自由度(DOF),其节点位移列阵为:
由于该单元有4个角节点,在每个棱边上的中点有一个中间节点(共6个),总共10个节点,由函数构造的Pascal三角形,取单元位移场模式为
u(x,y,z)=a1+a2x+a3y+a4z+a5xy+a6yz+a7xz+a8x2+a9y2+a10z2
以自然坐标表示,那么单位位移场模式为:
u(x,y,z)=a′1L1+a′2L2+a′3L3+a′4L4+a′5L1L2+a′6L2L3
+a′7L3L4+a′8L1L3+a′9L1L4+a′10L2L4
=N1u1+N2u2+N3u3+…+N10u10
对于角节点,相应的形状函数为:
Ni=(2Li-1)Li(i=1,2,3,4)
对于棱边上的中间节点,相应的形状函数为
N5=4L1L2,N6=4L1L2,N7=4L1L3
N8=4L1L4,N9=4L2L4,N10=4L3L4
得到形状函数后,按照有限元分析的标准过程即可推得相应的几何矩阵、刚度矩阵、节点等效载荷矩阵以及刚度方程;
步骤二、有限元模型的建立及其结果,具体方法为:
第一步、几何模型及网格划分:
对发动机箱体进行有限元建模时,采用虚拟拓扑方法,对整体力学性能影响小的几何细节进行整合;
第二步、材料属性及接触对定义,将夹具和工件定义为弹性体以提高分析精度,由于箱体尺寸大,采用弹性无接触模型来进行有限元分析;
有限元模型中,箱体与V型夹具都定义为弹性体,箱体的材料为铸铝合金ZL702A,弹性模量为70Gpa,泊松比为0.33,密度为2.7×103kg/m3,V型夹具材料为合金钢,弹性模型为210Gpa,泊松比为0.3。箱体和夹具均采用10节点四面体二次单元C3D10,其中V型夹具有82657个单元,1个警告单元(0.0012098%);箱体有382194个单元,2212个警告单元(0.5787%)。
第三步、约束定义及载荷施加,箱体具备6个自由度,其中5个自由度通过与V型夹具接触被限制住,余下为z方向的移动自由度,实际上它并不影响结合面的加工,但是为了有限元计算有解,限制箱体一端面的移动自由度来限制z向移动自由度,限制箱体左右两个缸面除z方向外的5个自由度。
刀具对箱体的加工过程是连续的,此处基于以下假设:刀具在某工况下其铣削力对箱体的作用是准静态,即可将连续铣削过程离散;刀具相对于铸铝合金箱体为刚性,铣削力在箱体上均匀分布;由于考虑铣削过程材料去除量相对于整个箱体体积极其小,不考虑加工过程中箱体刚度的变化;为计算方便,假设在铣削过程中,铣削力随着刀具的移动,只发生位置和方向的变化,其大小保持不变;
施加铣削力预测值中的最大值到有限元模型中并保持恒定,在X方向施加66N载荷,在Y方向上施加156.75N载荷,在Z方向上施加90.75N载荷。铣削力可施加在结合面的中间和两端位置。
步骤三、不同夹紧位置下装夹变形分析,具体方法为:
第一步、计算待夹紧面面积(计算结果为:2127.036mm2),按照夹紧力的计算方法,得到夹紧压强为4MPa并保持不变,夹紧力位置分别为箱体左右两侧的缸孔上的不规则面,分析不同工况下夹具对加工变形的产生的影响;
采用方案1进行装夹时,计算箱体在不同铣削力作用下的加工变形,结果,当铣削力施加于加工面右下端时最大变形量为22.3μm,当铣削力施加于加工面中间位置时其最大变形量为6.2μm,当铣削力施加于加工面右上端位置时其最大变形量为9.0μm。
改变夹紧力布局,计算得到工件加工表面不同的加工变形量,见表1。
由表可知,当铣削力位于加工面右侧中间位置时,加工变形量最小,当铣削力位于右侧下端位置时,加工变形量最大,右上端其次。在铣削力和装夹不变情况下,加工变形量随着加工位置的改变而变化。
表1不同装夹位置下的加工变形
在加工位置不变的情况下,随着夹紧位置的变化,加工变形量也发生着变化。加工面右下端处的变形量随着夹紧位置的变化先变大后减少,中间位置处变形量先减小后增大,右上端处变形量先增大后减少。根据计算结果,方案1在右下端及右上端处,相比另外两个方案变形量小,因此,方案1的装夹方案最优。
第二步、不同夹紧力下装夹变形分析,采用不同的三种装夹方案,各个方案下夹紧压强分别取3、4、5、6MPa,铣削力保持不变,计算箱体的加工变形,结果见表2-表4,分析表中数据,可得:当夹紧布局方案不变的情况下,当铣削力位于右下端和右上端时,工件的加工变形随着夹紧力的增加而减小。除了方案3外,当铣削力位于右中间时,加工变形量也随夹紧力增加而减小。统计各个工况变形量最小时对应的方案,得到装夹布局方案1最好,方案2其次,方案3不好。另外,方案1在6MPa夹紧载荷下,工件变形量最小。
表2方案1下工件加工变形
表3方案1下工件加工变形
表4方案1下工件加工变形
本发明在发动机上曲轴箱体结合面的铣削加工中,合理的夹紧布局方案和夹紧力对于加工变形的控制具有重要的意义。首先对装夹的约束条件进行了探讨,其次分析了箱体加工过程的装夹工艺以及装夹载荷的计算,然后建立了发动机箱体铣削加工的有限元模型,获得箱体的装夹变形云图,最后分析了工件在不同夹紧力和布局下的加工变形,为零件加工精度的提高提供依据。仿真结果表明:夹紧位置选择两端四缸孔方案最佳(方案1),夹紧载荷为6MPa时,加工变形量最小。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性的劳动即可做出的各种修改或变形仍在本发明的保护范围之内。
Claims (3)
1.一种箱体结构的简化及装夹变形仿真分析方法,其特征在于,所述的箱体结构的简化及装夹变形仿真分析方法包括:
步骤一、节点四面体二次单元,采用10节点四面体二次单元对发动机箱体进行有限元网格划分,10节点四面体二次单元即为在原4节点四面体单元的每一条棱边上再增加一个位于中点位置的内部节点,采用体积坐标应用Lagrange插值公式计算形状函数矩阵;
步骤二、有限元模型的建立及结果,包括几何模型及网格划分、材料属性及接触对定义、约束定义及载荷施加;
步骤三、不同夹紧位置下装夹变形分析,包括不同夹紧位置下装夹变形分析、不同夹紧力下装夹变形分析;
体积坐标的定义方式如下:对于给定的四面体,体内任何一点G的体积坐标:
等价于
Li=hi/Hi
其中,Hi为顶点i到对面的距离;hi为G点到i点对应面的距离;体积Gjmp为四面体Gjmp的体积值;Δ为给定四面体的体积;
10节点四面体单元共有30个节点位移自由度,节点位移列阵为:
<mrow>
<munder>
<msup>
<mrow>
<mo>{</mo>
<mi>q</mi>
<mo>}</mo>
</mrow>
<mi>e</mi>
</msup>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>&times;</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</munder>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>u</mi>
<mn>10</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>v</mi>
<mn>10</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>w</mi>
<mn>10</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
由于单元有4个角节点,在每个棱边上的中点有一个中间节点,共6个,总共10个节点,由函数构造的Pascal三角形,取单元位移场模式为:
u(x,y,z)=a1+a2x+a3y+a4z+a5xy+a6yz+a7xz+a8x2+a9y2+a10z2
以自然坐标表示,那么单位位移场模式为:
u(x,y,z)=a′1L1+a′2L2+a′3L3+a′4L4+a′5L1L2+a′6L2L3
+a′7L3L4+a′8L1L3+a′9L1L4+a′10L2L4
=N1u1+N2u2+N3u3+…+N10u10
对于角节点,相应的形状函数为:
Ni=(2Li-1)Li
i=1,2,3,4,对于棱边上的中间节点,相应的形状函数为:
N5=4L1L2, N6=4L1L2, N7=4L1L3
N8=4L1L4, N9=4L2L4, N10=4L3L4
得到形状函数后,按照有限元分析的标准过程即推得相应的几何矩阵、刚度矩阵、节点等效载荷矩阵以及刚度方程。
2.如权利要求1所述的箱体结构的简化及装夹变形仿真分析方法,其特征在于,有限元模型的建立及结果,具体方法为:
第一步、几何模型及网格划分:
对发动机箱体进行有限元建模时,采用虚拟拓扑方法,对整体力学性能影响小的几何细节进行整合;
第二步、材料属性及接触对定义,将夹具和工件定义为弹性体以提高分析精度,由于箱体尺寸大,采用弹性无接触模型来进行有限元分析;
第三步、约束定义及载荷施加,限制箱体一端面的移动自由度来限制z向移动自由度,限制箱体左右两个缸面除z方向外的5个自由度;
刀具对箱体的加工过程是连续的,刀具在某工况下其铣削力对箱体的作用是准静态,即将连续铣削过程离散;刀具相对于铸铝合金箱体为刚性,铣削力在箱体上均匀分布;由于考虑铣削过程材料去除量相对于整个箱体体积极其小,不考虑加工过程中箱体刚度的变化;在铣削过程中,铣削力随着刀具的移动,只发生位置和方向的变化,大小保持不变。
3.如权利要求1所述的箱体结构的简化及装夹变形仿真分析方法,其特征在于,不同夹紧位置下装夹变形分析,具体方法为:
第一步、计算待夹紧面面积,按照夹紧力的计算方法,得到夹紧压强并保持不变,夹紧力位置分别为箱体左右两侧的缸孔上的不规则面,分析不同工况下夹具对加工变形的产生的影响;
第二步、不同夹紧力下装夹变形分析,采用不同的三种装夹方案,各个方案下夹紧压强不同,铣削力保持不变,计算箱体的加工变形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410508770.7A CN104252566B (zh) | 2014-09-28 | 2014-09-28 | 一种箱体结构的简化及装夹变形仿真分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410508770.7A CN104252566B (zh) | 2014-09-28 | 2014-09-28 | 一种箱体结构的简化及装夹变形仿真分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104252566A CN104252566A (zh) | 2014-12-31 |
CN104252566B true CN104252566B (zh) | 2017-12-19 |
Family
ID=52187454
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410508770.7A Active CN104252566B (zh) | 2014-09-28 | 2014-09-28 | 一种箱体结构的简化及装夹变形仿真分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104252566B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016170552A1 (en) | 2015-04-22 | 2016-10-27 | Alesamonti S.R.L. | Method to minimize workpiece deformations due to its clamping on a machine tool |
CN106484960A (zh) * | 2016-09-20 | 2017-03-08 | 江苏科技大学 | 一种抑制柴油机关键孔系加工变形的工艺参数优化方法 |
CN107992638B (zh) * | 2016-10-26 | 2021-07-09 | 北京汽车动力总成有限公司 | 一种发动机曲轴箱结构模型的建立方法及装置 |
CN107451371B (zh) * | 2017-08-11 | 2021-01-26 | 四川大学 | 提高三角单元计算精度的方法 |
CN109214051B (zh) * | 2018-08-03 | 2020-11-24 | 上海交通大学 | 基于直列发动机缸体的装夹优化方法及系统 |
CN113703395B (zh) * | 2021-07-07 | 2022-10-21 | 西北工业大学 | 一种面向加工变形控制的变工况铣削加工夹紧力预测方法 |
CN113779707B (zh) * | 2021-11-11 | 2022-02-01 | 成都航空职业技术学院 | 基于位移载荷仿真分析的超大型薄壁零件变形控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510229A (zh) * | 2009-03-20 | 2009-08-19 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
EP2237175A1 (en) * | 2009-04-01 | 2010-10-06 | Livermore Software Technology Corporation | Explosion simulation in finite element analysis |
-
2014
- 2014-09-28 CN CN201410508770.7A patent/CN104252566B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510229A (zh) * | 2009-03-20 | 2009-08-19 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
EP2237175A1 (en) * | 2009-04-01 | 2010-10-06 | Livermore Software Technology Corporation | Explosion simulation in finite element analysis |
Non-Patent Citations (2)
Title |
---|
汽车发动机缸体加工变形分析与精度控制研究;雷锋杰;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120715;摘要,第4-5、16-17、24、29-35、44页 * |
矿用隔爆变压器箱体的有限元分析及结构优化;黄向明;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20051115;第15、70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104252566A (zh) | 2014-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104252566B (zh) | 一种箱体结构的简化及装夹变形仿真分析方法 | |
Hong et al. | Static structural analysis of great five-axis turning–milling complex CNC machine | |
Cheung et al. | Finite element implementation | |
US10838403B2 (en) | Simulation method for milling by use of dynamic position error | |
CN106484960A (zh) | 一种抑制柴油机关键孔系加工变形的工艺参数优化方法 | |
Ahmadi | Finite strip modeling of the varying dynamics of thin-walled pocket structures during machining | |
CN106334972A (zh) | 一种球头铣刀平面加工的刀刃接触判据方法 | |
CN106294977A (zh) | 一种机器人铣削加工中工件装夹位置优划方法 | |
Ghorbani et al. | An analytical model for cutter-workpiece engagement calculation in ball-end finish milling of doubly curved surfaces | |
CN111008493B (zh) | 一种砂轮磨削的仿真方法 | |
CN107322008A (zh) | 一种车外圆工艺中已加工表面的三维建模方法 | |
CN108021756B (zh) | 一种基于d-h表示法的可重构机床运动学误差模型建立方法 | |
Sun et al. | Iso-planar feed vector-fields-based streamline tool path generation for five-axis compound surface machining with torus-end cutters | |
Ozlu et al. | Comparison of one-dimensional and multi-dimensional models in stability analysis of turning operations | |
CN106096146A (zh) | 切削过程中薄壁件动力学参数的预测方法 | |
CN108416153A (zh) | 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 | |
Rai et al. | FEM-based prediction of workpiece transient temperature distribution and deformations during milling | |
JP4346630B2 (ja) | 加工条件取得装置およびそのプログラム | |
Cho et al. | Volumetric error analysis of a multi-axis machine tool machining a sculptured surface workpiece | |
CN107016160A (zh) | 三维椭圆振动辅助切削装置设计方法 | |
Liu et al. | Final machining of large-scale engine block with modularized fixture and virtual manufacturing technologies | |
CN111061218B (zh) | 一种基于acis的复杂曲面加工切削力预测方法及设备 | |
Zhang et al. | Tool path regeneration for mold design modification | |
Tsai et al. | The mathematical models of the basic entities of multi-axis serial orthogonal machine tools using a modified Denavit–Hartenberg notation | |
CN108334688B (zh) | 基于matlab的旋转功能梯度厚板动力学行为计算方法 |
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 |