CN105868446B - 一种高速冲压装备关键部件热态特性的跨尺度仿真方法 - Google Patents
一种高速冲压装备关键部件热态特性的跨尺度仿真方法 Download PDFInfo
- Publication number
- CN105868446B CN105868446B CN201610172556.8A CN201610172556A CN105868446B CN 105868446 B CN105868446 B CN 105868446B CN 201610172556 A CN201610172556 A CN 201610172556A CN 105868446 B CN105868446 B CN 105868446B
- Authority
- CN
- China
- Prior art keywords
- micro
- contact
- mrow
- bulge
- critical component
- 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
- 239000007787 solid Substances 0.000 claims abstract description 40
- 238000012546 transfer Methods 0.000 claims abstract description 33
- 230000007246 mechanism Effects 0.000 claims abstract description 25
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 16
- 238000010438 heat treatment Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 238000013459 approach Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 10
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 125000000174 L-prolyl group Chemical group [H]N1C([H])([H])C([H])([H])C([H])([H])[C@@]1([H])C(*)=O 0.000 claims description 4
- 238000004080 punching Methods 0.000 claims description 4
- 230000001788 irregular Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 4
- 238000010998 test method Methods 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 239000010687 lubricating oil Substances 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 206010037660 Pyrexia Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000314 lubricant Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种高速冲压装备关键部件热态特性的跨尺度仿真方法。该方法首先确定高速冲压装备关键部件在工作过程中的主要热源及其发热量。然后确定影响关键部件温度场分布的热传递途径与热传递方式，计算对流换热系数。获取关键部件固体粗糙表面的分形维数和特征长度，模拟固体粗糙表面，综合考虑弹性、塑性、弹塑性三种不同变形机制，计算其接触热阻，从而得到各类固体接触面的接触传热系数。最后对高速冲压装备关键部件进行热态特性有限元仿真分析。本发明考虑多变形机制，在考虑动摩擦因素情况下，将变形分为三种情况分别计算基体热阻，使得仿真方法对高速冲压装备关键部件热态特性的分析更加全面、精确，更加符合实际情况。
Description
技术领域
本发明涉及高速冲压技术领域，尤其涉及一种高速冲压装备关键部件热态特性的跨尺度仿真方法。
背景技术
高速冲压装备关键部件热态特性分析主要有试验法和仿真方法。试验法需进行大量实验来完成整机温度场的测量，温度信息采集时，无法整机都布置传感器，而通常是在冲压装备某些部位进行布置，传感器的位置直接影响温度信息的准确釆集，最终影响热分析模型的精度。此外试验法需要的设备较多，成本较昂贵，不具备通用性。仿真方法借助现有的有限元软件对高速冲压装备关键部件进行仿真建模与分析计算，能够高效低成本地获得高速冲压装备关键部件内部温度场的变化规律与瞬态热平衡时间，在工程中得到广泛应用。
要保证高速冲压装备关键部件热态特性仿真分析结果的有效性，其关键在于准确计算热态特性仿真分析时各类热环境参数，构建与实际相符的热态特性仿真分析模型，其中，固体接触面接触传热参数即接触热阻的准确计算是难点。接触热阻作为决定固体接触面间热传递的一个重要参数，对高速冲压装备关键部件的热态特性有着直接而重要的影响，因此，高速冲压装备关键部件的热态特性仿真分析时必须充分考虑固体接触面接触热阻的影响。
高速冲压装备运行过程中，滚动轴承、滑动轴承连接处作为主要热源部分，其接触热阻对温度场的截流作用更为明显，需充分考虑其接触面的微观形貌特征和弹性、塑性、弹塑性三种形变方式，建立与实际相符的固体接触面的接触热阻计算模型，准确计算各接触面的接触传热参数，从而保证高速冲压装备关键部件热态特性仿真分析结果的准确性。
发明内容
本发明的目的在于针对现有技术的不足，提供一种高速冲压装备关键部件热态特性的跨尺度仿真方法。首先根据高速冲压装备关键部件的具体运行情况确定该关键部件在工作过程中的主要热源，并计算其发热量。然后根据各部件接触面的具体接触类型确定热传递方式，计算强制对流换热系数和自然对流换热系数。在考虑固体接触面真实微观形貌与多边形机制的情况下，给出各类固体接触面的接触热阻计算公式，从而得到相应的接触传热系数。最后利用Pro/E建立关键部件三维模型，导入ANSYS workbench中，划分网格，设置材料属性和计算得到的热环境参数，对高速冲压装备关键部件的热态特性进行有限元仿真分析。
为实现上述目的，本发明采用的技术方案是：一种高速冲压装备关键部件热态特性的跨尺度仿真方法，该方法包括以下步骤：
1)确定影响高速冲压装备关键部件温度场分布的各主要热源，并计算各主要热源的发热量。
2)确定影响高速冲压装备关键部件温度场分布的热传递途径与热传递方式，计算出对流换热系数。
3)细观尺度下考虑多变形机制的固体接触粗糙表面接触传热系数计算，具体包括以下子步骤：
3.1)利用功率谱法与W-M分形函数获取高速冲压装备关键部件固体接触粗糙表面的分形维数D和特征长度G，利用分形-蒙特卡罗法模拟冲压装备关键部件的固体接触粗糙表面；所述粗糙表面由微凸峰不规则排列构成，所述微凸峰由微凸体经自仿射变换构造新微凸体堆叠，再由新生成微凸体经自仿射变换构造新微凸体堆叠，如此往复，直到新生成的微凸体的微接触直径小于设定值；
3.2)计算关键部件固体接触粗糙表面上构成微凸峰的微凸体的收缩热阻和该微凸峰的不同面积接触点间热阻；
3.3)将关键部件固体接触粗糙表面上构成微凸峰的微凸体，根据其微接触半径的大小，分成弹性变形微凸体、弹塑性变形微凸体和塑性变形微凸体；具体分类方法如下：计算弹性变形的临界微接触半径和塑性变形的临界微接触半径；将两个临界微接触半径作为分段点，若微凸体的微接触半径小于弹性变形的临界微接触半径，则将其归类为弹性变形微凸体；若微凸体的微接触半径大于塑性变形的临界微接触半径，则将其归类为塑性变形微凸体；其余归为弹塑性变形微凸体；
3.4)针对发生不同形变的微凸体，分别计算其基体热阻；
3.5)将固体接触粗糙表面上构成微凸峰的微凸体的基体热阻、收缩热阻和该微凸峰的不同面积接触点间热阻串联，得到高速冲压装备关键部件固体接触粗糙表面上各微凸峰的接触热阻；
3.6)将固体接触粗糙表面上所有的微凸峰的接触热阻并联，计算得到关键部件固体接触粗糙表面的接触热阻。
3.7)将关键部件固体接触粗糙表面的接触热阻计算结果取倒数，得到固体接触粗糙表面的接触传热系数。
4)宏观尺度下高速冲压装备关键部件热态特性仿真分析：利用Pro/E构建关键部件的三维模型，导入ANSYS workbench，并划分网格，设置部件中各零件的材料属性。根据步骤1)、2)、3)的计算结果，设置关键部件热态特性仿真所需的发热量、对流换热系数和接触传热系数，在宏观尺度下对高速冲压装备关键部件进行热态特性仿真分析。
进一步地，基于多变形机制的弹性形变的临界微接触半径R1和塑性变形的临界微接触半径R2的计算公式如下：
式中：σy为材料屈服强度；K是两接触材料各自硬度与屈服强度的比值的较小值；l为微凸体变形前截面横向宽度；E为有效弹性模量；G为特征长度；D为分形维数。
本发明的有益效果是：在考虑动摩擦因素的情况下，综合考虑了微凸体的弹性、弹塑性、完全塑性三种变形情况下基体热阻应有不同的计算方式。同时，接触面上的热阻网络模型综合考虑了基体热阻、收缩热阻以及不同面积接触点之间的热阻，使计算模型更全面，使得所提出仿真方法对高速冲压装备关键部件热态特性的分析更加全面、精确，更加符合实际情况。
附图说明
图1为高速冲压装备关键部件热态特性的跨尺度仿真流程图；
图2为高速冲压装备施力机构1/2简化对称结构图；
图3为用分形-蒙特卡罗法产生粗糙表面的流程图；
图4为跨尺度仿真方法所得的高速冲压装备施力机构的温度场变化趋势图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步说明。
高速冲压装备关键部件热态特性的跨尺度仿真流程如图1所示。以高速冲压装备施力机构热态特性仿真分析为例，说明高速冲压装备关键部件热态特性跨尺度仿真的主要步骤。高速冲压装备施力机构主要由主轴、连杆、上横梁、滑块、轴承等构成。图2是本发明相应的高速冲压装备施力机构1/2简化对称结构图。仿真方法具体如下：
1)高速冲压装备施力机构热源分析与发热量计算；
根据高速冲压装备的具体运行情况确定影响高速冲压装备施力机构温度场分布的主要热源来自滑动轴承发热量和滚动轴承发热量。
上横梁与主轴的轴瓦、连杆与主轴的轴瓦实际均是滑动轴承，与滚动轴承相似，滑动轴承工作时，轴承的热量由摩擦功转变而来。轴承每秒产生的热量可利用轴承的摩擦因数，轴承所受载荷和轴承圆周线速度计算求得。
滚动轴承发热由外加载荷摩擦力矩引起的滚动体与滚道的摩擦生热，由润滑剂粘度力矩引起的克服滚动体粘性阻力而产生的粘性摩擦生热和由滚动体在与内外圈滚道接触区上发生自旋产生的自旋摩擦生热所构成，可分别计算各摩擦热从而得到滚动轴承的发热量。
高速冲压装备施力机构各热源发热量计算结果如表1所示，W为单位——瓦特。
表1
2)高速冲压装备施力机构热传递分析与对流换热系数计算；
2.1)高速冲压装备施力机构热传递分析
确定影响高速冲压装备施力机构温度场分布的热传递途径与热传递方式，计算出对流换热系数。
热量的传递是由温差引起的，在温差作用下，热量自发地从高温物体传向低温物体，实现热量的传递。高速冲压装备施力机构主要热源的热传递方式如表2所示。
表2
2.2)高速冲压装备施力机构对流换热系数计算
高速冲压装备施力机构对流换热主要包括：主轴上轴承接触部分模拟润滑油的强制热对流、主轴与空气的强制热对流、高速冲压装备外壳与空气的自然热对流。热对流遵循牛顿冷却公式，其大小由对流换热系数表示。由雷诺数判定冷却液流动状态，并依据努塞尔数可计算出主轴不同部件的对流换热系数。
强制流换热系数可以通过部件周围流体的导热系数、努塞尔数、轴承所在轴颈的直径或者主轴当量直径进行求解。
自然流换热系数可以通过格拉晓夫数和普朗特数进行求解。
高速冲压装备施力机构对流换热系数计算结果如表3所示，W、m和K为单位——瓦特、米和开尔文。
表3
| 热对流 | (W/(m2·K)) |
| 模拟润滑油强制热对流 | 204.9 |
| 主轴旋转强制热对流 | 85.2 |
| 压力机外壳自然热对流 | 5 |
3)细观尺度下考虑多变形机制的固体接触粗糙表面接触传热系数计算，具体包括以下子步骤：
3.1)利用功率谱法与W-M分形函数获取高速冲压装备关键部件固体接触粗糙表面的分形维数D和特征长度G，利用分形-蒙特卡罗法模拟冲压装备关键部件的固体接触粗糙表面；所述粗糙表面由微凸峰不规则排列构成，所述微凸峰由微凸体经自仿射变换构造新微凸体堆叠，再由新生成微凸体经自仿射变换构造新微凸体堆叠，如此往复，直到新生成的微凸体的微接触直径小于设定值；
3.2)计算关键部件固体接触粗糙表面上构成微凸峰的微凸体的收缩热阻和该微凸峰的不同面积接触点间热阻；
3.3)将关键部件固体接触粗糙表面上构成微凸峰的微凸体，根据其微接触半径的大小，分成弹性变形微凸体、弹塑性变形微凸体和塑性变形微凸体；具体分类方法如下：计算弹性变形的临界微接触半径R1和塑性变形的临界微接触半径R2，计算公式如下：
式中：σy为材料屈服强度；K是两接触材料各自硬度与屈服强度的比值的较小值；l为微凸体变形前截面横向宽度；E为有效弹性模量；G为特征长度；D为分形维数。
将两个临界微接触半径作为分段点，若微凸体的微接触半径小于弹性变形的临界微接触半径，则将其归类为弹性变形微凸体；若微凸体的微接触半径大于塑性变形的临界微接触半径，则将其归类为塑性变形微凸体；其余归为弹塑性变形微凸体；
3.4)针对发生不同形变的微凸体，分别计算其基体热阻；
3.5)将固体接触粗糙表面上构成微凸峰的微凸体的基体热阻、收缩热阻和该微凸峰的不同面积接触点间热阻串联，得到高速冲压装备关键部件固体接触粗糙表面上各微凸峰的接触热阻；
3.6)将固体接触粗糙表面上所有的微凸峰的接触热阻并联，计算得到关键部件固体接触粗糙表面的接触热阻。
3.7)将关键部件固体接触粗糙表面的接触热阻计算结果取倒数，得到固体接触粗糙表面的接触传热系数，计算结果如表4所示，W、m和K为单位——瓦特、米和开尔文。
表4
4)宏观尺度下高速冲压装备施力机构热态特性仿真分析：利用Pro/E构建高速冲压装备施力机构的三维模型，将其导入ANSYS workbench中，根据需要确定网格大小和划分方式。网格划分后，设置材料属性，并根据步骤2)和3)计算结果，设置该热态特性仿真模型各热源的发热量、对流换热系数和接触传热系数，对高速冲压装备施力机构热态特性进行仿真分析，仿真结果如图4所示。
Claims (1)
1.一种高速冲压装备关键部件热态特性的跨尺度仿真方法，其特征在于，该方法包括以下步骤：
1)确定影响高速冲压装备关键部件温度场分布的各主要热源，并计算各主要热源的发热量；
2)确定影响高速冲压装备关键部件温度场分布的热传递途径与热传递方式，计算出对流换热系数；
3)细观尺度下考虑多变形机制的固体接触粗糙表面接触传热系数计算，具体包括以下子步骤：
3.1)利用功率谱法与W-M分形函数获取高速冲压装备关键部件固体接触粗糙表面的分形维数D和特征长度G，利用分形-蒙特卡罗法模拟冲压装备关键部件的固体接触粗糙表面；所述粗糙表面由微凸峰不规则排列构成，所述微凸峰由微凸体经自仿射变换构造新微凸体堆叠，再由新生成微凸体经自仿射变换构造新微凸体堆叠，如此往复，直到新生成的微凸体的微接触直径小于设定值；
3.2)计算关键部件固体接触粗糙表面上构成微凸峰的微凸体的收缩热阻和该微凸峰的不同面积接触点间热阻；
3.3)将关键部件固体接触粗糙表面上构成微凸峰的微凸体，根据其微接触半径的大小，分成弹性变形微凸体、弹塑性变形微凸体和塑性变形微凸体；具体分类方法如下：计算弹性变形的临界微接触半径和塑性变形的临界微接触半径；基于多变形机制的弹性变形的临界微接触半径R1和塑性变形的临界微接触半径R2的计算公式如下：
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K&sigma;</mi>
<mi>y</mi>
</msub>
<msup>
<mi>l</mi>
<mi>D</mi>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&pi;EG</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>110</mn>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</msup>
<msub>
<mi>K&sigma;</mi>
<mi>y</mi>
</msub>
<msup>
<mi>l</mi>
<mi>D</mi>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&pi;EG</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mfrac>
</mrow>
式中：σy为材料屈服强度；K为两接触材料各自硬度与屈服强度的比值的较小值；l为微凸体变形前截面横向宽度；E为有效弹性模量；G为特征长度；D为分形维数；
将两个临界微接触半径作为分段点，若微凸体的微接触半径小于弹性变形的临界微接触半径，则将其归类为弹性变形微凸体；若微凸体的微接触半径大于塑性变形的临界微接触半径，则将其归类为塑性变形微凸体；其余归为弹塑性变形微凸体；
3.4)针对发生不同形变的微凸体，分别计算其基体热阻；
3.5)将固体接触粗糙表面上构成微凸峰的微凸体的基体热阻、收缩热阻和该微凸峰的不同面积接触点间热阻串联，得到高速冲压装备关键部件固体接触粗糙表面上各微凸峰的接触热阻；
3.6)将固体接触粗糙表面上所有的微凸峰的接触热阻并联，计算得到关键部件固体接触粗糙表面的接触热阻；
3.7)将关键部件固体接触粗糙表面的接触热阻计算结果取倒数，得到固体接触粗糙表面的接触传热系数；
4)宏观尺度下高速冲压装备关键部件热态特性仿真分析：利用Pro/E构建关键部件的三维模型，导入ANSYS workbench，并划分网格，设置部件中各零件的材料属性；根据步骤1)、2)、3)的计算结果，设置关键部件热态特性仿真所需的发热量、对流换热系数和接触传热系数，在宏观尺度下对高速冲压装备关键部件进行热态特性仿真分析。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610172556.8A CN105868446B (zh) | 2016-03-24 | 2016-03-24 | 一种高速冲压装备关键部件热态特性的跨尺度仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610172556.8A CN105868446B (zh) | 2016-03-24 | 2016-03-24 | 一种高速冲压装备关键部件热态特性的跨尺度仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105868446A CN105868446A (zh) | 2016-08-17 |
| CN105868446B true CN105868446B (zh) | 2018-05-29 |
Family
ID=56625555
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610172556.8A Active CN105868446B (zh) | 2016-03-24 | 2016-03-24 | 一种高速冲压装备关键部件热态特性的跨尺度仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105868446B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109084975B (zh) * | 2018-09-29 | 2020-02-14 | 南京理工大学 | 一种厚壁圆筒内壁磨损定量性研究方法 |
| CN110705147B (zh) * | 2019-09-18 | 2023-08-04 | 北京工业大学 | 一种数控机床主轴热态特性综合性理论建模与分析方法 |
| CN111177842B (zh) * | 2020-01-03 | 2022-09-09 | 中铁第四勘察设计院集团有限公司 | 一种盾构隧道管片接触面的动力特性确定方法、装置及存储介质 |
| CN114580222B (zh) * | 2022-05-09 | 2022-07-22 | 山东大学 | 用于电子器件的高速瞬态热仿真方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011065238A (ja) * | 2009-09-15 | 2011-03-31 | Ecogold Co Ltd | 構造物内温度計算方法及び当該計算方法を実施するプログラム |
| CN104820748A (zh) * | 2015-05-07 | 2015-08-05 | 北京宇航***工程研究所 | 一种运载火箭大气层内飞行段舱段温度场分布确定方法 |
| CN104834760A (zh) * | 2014-02-12 | 2015-08-12 | 上海通用汽车有限公司 | 提取发动机缸孔各阶次变形的处理方法 |
| CN105033073A (zh) * | 2015-08-07 | 2015-11-11 | 苏州明和行新材料科技有限公司 | 低温润滑供给装置、具有其的高速冲压模具及其控制方法 |
-
2016
- 2016-03-24 CN CN201610172556.8A patent/CN105868446B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011065238A (ja) * | 2009-09-15 | 2011-03-31 | Ecogold Co Ltd | 構造物内温度計算方法及び当該計算方法を実施するプログラム |
| CN104834760A (zh) * | 2014-02-12 | 2015-08-12 | 上海通用汽车有限公司 | 提取发动机缸孔各阶次变形的处理方法 |
| CN104820748A (zh) * | 2015-05-07 | 2015-08-05 | 北京宇航***工程研究所 | 一种运载火箭大气层内飞行段舱段温度场分布确定方法 |
| CN105033073A (zh) * | 2015-08-07 | 2015-11-11 | 苏州明和行新材料科技有限公司 | 低温润滑供给装置、具有其的高速冲压模具及其控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105868446A (zh) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105868446B (zh) | 一种高速冲压装备关键部件热态特性的跨尺度仿真方法 | |
| Abbasi et al. | Temperature and thermoelastic instability at tread braking using cast iron friction material | |
| Haitao et al. | Simulation of thermal behavior of a CNC machine tool spindle | |
| CN103366050B (zh) | 一种高速电主轴力-热耦合建模方法 | |
| CN105022900B (zh) | 基于热固耦合分析的重型数控立车静压转台结构优化方法 | |
| Wang et al. | Investigation into the meshing friction heat generation and transient thermal characteristics of spiral bevel gears | |
| CN104978465B (zh) | 一种考虑微凸体的弹塑性变形和空气介质热阻的接触热阻建模方法 | |
| Meng et al. | Research on thermal resistance network modeling of motorized spindle based on the influence of various fractal parameters | |
| Li et al. | Calculation method of convective heat transfer coefficients for thermal simulation of a spindle system based on RBF neural network | |
| Li et al. | Assessment of six turbulence models for modeling and predicting narrow passage flows, part 1: Impingement jets | |
| CN108414105A (zh) | 一种用于干式摩擦副的仿真与试验耦合的动态测温方法 | |
| Luo et al. | Flow and forced-convection characteristics of turbulent flow through parallel plates with periodic transverse ribs | |
| Liu et al. | Experimental and numerical studies on the excitation loads and vortex structures of four circular section cylinders in a square configuration | |
| CN105302964A (zh) | 一种用于芯片结构的热分布分析方法 | |
| CN107862130B (zh) | 一种导致继电器贮存退化的内部关键因素的分析方法 | |
| Dai et al. | Temperature measurement point optimization and experimental research for bi-rotary milling head of five-axis CNC machine tool | |
| CN107526868A (zh) | 一种用于雷达电子机柜***的热设计方法 | |
| CN109726477A (zh) | 一种间隙转动副有限元建模及分析方法 | |
| Ozkol et al. | Entropy generation in laminar natural convection from a constant temperature vertical plate in an infinite fluid | |
| Tang et al. | Fluid analysis and improved structure of an ATEG heat exchanger based on computational fluid dynamics | |
| Iftikhar et al. | Large eddy simulation of the flow past a soccer ball | |
| CN113884266A (zh) | 一种等效式高速轴承流场参数的测量方法 | |
| Fan et al. | Prediction of heat dissipation effect by simulating radiator group of mining dump-truck | |
| Wu et al. | Transient mixed convective heat transfer predictions around three heated cylinders in a horizontal channel | |
| Brecher et al. | Metrological investigation and simulation of thermo-mechanical interactions in externally driven spindles |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20160817 Assignee: Huzhou Machine Tool Factory Assignor: Zhejiang University Contract record no.: 2019330000026 Denomination of invention: Cross-scale simulation method for thermal-state characteristic of key part of high-speed stamping device Granted publication date: 20180529 License type: Common License Record date: 20190307 |
|
| EE01 | Entry into force of recordation of patent licensing contract |