《一种重型数控落地铣镗床热特性数值仿真方法.pdf》由会员分享，可在线阅读，更多相关《一种重型数控落地铣镗床热特性数值仿真方法.pdf（12页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104200019A43申请公布日20141210CN104200019A21申请号201410427702822申请日20140827G06F17/5020060171申请人电子科技大学地址611731四川省成都市高新区（西区）西源大道2006号72发明人黄智李俊杰杜丽陈学尚衡凤琴许可王正杰74专利代理机构成都顶峰专利事务所普通合伙51224代理人李崧岩54发明名称一种重型数控落地铣镗床热特性数值仿真方法57摘要本发明公开了一种重型数控落地铣镗床热特性数值仿真方法，解决了现有重型数控机床热态特性有限元数值分析方法的分析结果不准确、精度低的问题。该数字仿真方法将机床作为一个整。
2、体进行热特性数值仿真，再结合对各零部件结合面的接触热阻和静压导轨油膜的产热，有效地提高了分析结果的准确性和精度。本发明分析结果准确性和精度高，为其大范围的推广应用，奠定了坚实的基础。51INTCL权利要求书3页说明书7页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书3页说明书7页附图1页10申请公布号CN104200019ACN104200019A1/3页21一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，包括以下步骤1确定重型数控落地铣镗床的热源分布，分析机床与环境之间的热交换；2在三维建模软件PRO/E中建立重型数控落地铣镗床装配体模型，并对装配体模型进行简化；3。
3、通过PRO/E和有限元软件WORKBENCH的无缝接口，将简化后的装配体模型导入有限元软件WORKBENCH中，并进行网络划分；4根据重型数控落地铣镗床中各个轴承的径向和轴向刚度计算等效圆柱环的等效径向弹性模量和轴向弹性模量，并将等效径向弹性模量和轴向弹性模量导入WORKBENCH中；5根据重型数控落地铣镗床的材料设置零部件物理性能参数，计算步骤1中所确定的各热源大小和对流传热系数，将计算结果作为边界条件施加在机床上，然后利用有限元软件WORKBENCH对重型数控落地铣镗床分别进行稳态温度场分析和瞬态温度场分析，得出重型数控落地铣镗床温度分布和温升的变化情况，为重型数控落地铣镗床的热位移分析奠。
4、定基础；6根据重型数控落地铣镗床的实际安装位置，在床身底面施加固定约束，将重型机床温度场分析的温度作为载荷，采用间接法，基于热结构耦合分析重型数控落地铣镗床的变形场，通过在有限元软件WORKBENCH中求解计算，得到重型数控落地铣镗床的热变形场云图。2根据权利要求1所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，所述步骤3中装配体模型的简化内容如下31去除零件中的孔，去掉螺钉和螺栓孔；32去除大型铸件的倒角、阶梯孔，忽略结构中的台阶；33将轴承和滚珠丝杠副简化成一个相同尺寸的等截面圆柱体；34去除主轴部件内部的拉刀机构；35去除机床各进给传动电机。3根据权利要求1所述的一种重型数控。
5、落地铣镗床热特性数值仿真方法，其特征在于，所述网络划分的方法为自由分网。4根据权利要求1所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，所述等效圆柱环的弹性模量的确定方法为径向等效弹性模量的确定首先，在WORKBENCH结构分析模块中，对轴承模型外圈施加全约束，并在轴与轴承内圈定义接触类型；其次，对轴施加径向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承的刚度变化；再其次，根据刚度计算公式，获得轴承有限元计算径向刚度；当轴承有限元计算径向刚度与简化前轴承刚度之差最小时，轴承对应的弹性模量即为轴承的径向等效弹性模量；轴向等效弹性模量的确定首先，在。
6、WORKBENCH结构分析模块中，对轴承模型的一侧约束自由度，轴承模型的外圈约束出轴向以外的自由度；其次，对轴承模型的另外一侧施加轴向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承轴向的刚度变化；再其次，根据刚度计算公式，获得轴承有限元计算轴向刚度；当滚动轴承有限元计算轴向刚度权利要求书CN104200019A2/3页3与简化前轴承刚度之差最小时，轴承对应的弹性模量即为轴承的轴向等效弹性模量。5根据权利要求1所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，所述步骤5中机床承受的热载荷还包括由于重型数控落地铣镗床各零部件结合面不完全接触而产生的热。
7、量损失。6根据权利要求5所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，仿真热量损失的方法如下通过分析重型数控落地铣镗床的各零件的结合面的加工和装配工艺，确定各结合面的换热系数，然后通过采用在有限元软件中对装配体模型的接触属性设置换热系数的方法仿真热量损失。7根据权利要求6所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，重型数控落地铣镗床零部件之间的接触类型定义如下主轴轴承与主轴及轴套之间的接触为绑定，立柱与滑座之间的接触类型为绑定，其余零部件之间接触类型均为不分离。8根据权利要求1所述的一种重型数控落地铣镗床热特性数值仿真方法，其特征在于，所述热源包括滚动轴承产热、滚。
8、珠丝杠副产热、电机产热、静压导轨油膜发热，热源和对流传热系数的计算方式如下A滚动轴承产热计算Q11047MN1式中N为轴承转速，单位为R/MIN；M为滚动轴承摩擦力矩，单位为NM，其中，MM0M1，M0反映和轴承所受载荷无关的损失的摩擦力矩；M1反映与轴承所受载荷滑动摩擦相关的损失的摩擦力矩；B滚珠丝杠副产热计算Q212NM2式中N为滚珠丝杠转速，单位为R/MIN；M为螺母的摩擦力矩，单位为NM；C电机产热计算式中NM是驱动电机的功率，单位为W；是驱动电机的效率；MM为驱动电机输出力矩，单位为NM；N为驱动电机转速，单位为R/MIN；D静压导轨油膜发热计算支承滑动面上的剪切应力为/H4支承滑动。
9、面上的摩擦力为FFAA/H5摩擦功率损失为PFFFA2/H6式中为油液的动力粘度，单位为PAS；A为密封带部分的支承面积，单位为M2；E对流传热系数计算式中NU为努谢尔特准则数；为对流传热系数，单位为W/M2；为流体的导热权利要求书CN104200019A3/3页4系数，单位为W/M；D为放热壁面的定型尺寸，单位为M。权利要求书CN104200019A1/7页5一种重型数控落地铣镗床热特性数值仿真方法技术领域0001本发明涉及一种分析方法，具体的说，是涉及一种重型数控落地铣镗床热特性数值仿真方法。背景技术0002重型数控机床是我国能源、航天航空、大型舰船、汽车等行业发展的急需，是衡量一个国家技。
10、术水平和综合国力的战略物资。随着精密机械加工制造技术的发展，关于重型机床机械加工精度的要求越来越高。然而，仅就重型机床的加工精度而言，国内的产品与国际企业的产品还存在较大差距，重型机床精度差距制约了我国制造业的发展，提高重型机床加工精度已经成为亟待解决的我国机械制造业关键问题之一。经过研究分析表明热变形最终导致刀尖点偏移误差占重型机床加工误差的2575。因此，控制热变形误差是机床的加工精度的保证，也是提高其加工精度的关键保障技术。0003目前，利用有限元方法对于中小型数控机床的热变形特性分析研究较多，但针对重型数控落地铣镗床整机的的热特性分析研究较少，主要是重型机床所处的工况较为恶劣及多变，外。
11、界环境条件与内部发热状况与中小型机床相比更复杂，现有技术并没有专门针对重型机床热变化特性的准确定量计算的成熟、可行的方法。目前的重型数控机床热态特性有限元数值分析方法主要有00041关键零件热特性数值分析方法，主要是针对关键零件如主轴、主轴箱、滑枕、立柱进行温度场和变形场数值仿真，或将关键部件的若干零件通过布尔运算粘成一个零件，不考虑零件之间的配合进行数值仿真。该方法仅仅能分析单个零件的温度场和变形场，并且忽略了其他零件对它的影响，数值分析结果精度低，无法预知整体部件或者整机的温度场与变形场情况。00052关键部件热特性数值分析方法，主要计算重型数控落地铣镗床的关键部件如主轴部件、滚珠丝杠部件。
12、的温度场和热变形。该方法首先将零件模型进行装配，并对装配体进行整体求解，其结果精度可大为提高。所以这种是重型数控落地铣镗床热特性数值仿真的常用方法。00063整机热特性数值分析方法，主要针对重型数控落地铣镗床整机的温度场和热变形进行数值仿真。该方法首先建立整机的装配体模型，并对三维模型进行简化求解运算，充分考虑零件之间的热传递，各构件之间的约束与配合，确保数值仿真的准确度。这种方法对重型数控落地铣镗床热特性数值仿真也更符合实际情况。0007现有技术存在的缺陷主要包括以下三个方面00081忽略静压导轨油膜产热。对于重型数控落地铣镗床，滑座在床身上的移动、主轴箱沿立柱移动和滑枕在主轴箱内移动，均装。
13、配有静压导轨。在重型数控落地铣镗床热源分析和计算中，目前上述方法均忽略了静压导轨油膜产热。但在工程实践过程中，以上部件在静压导轨上移动会因为油膜粘性阻力产生热量，由此引发在静压导轨配合处产生温度升高，进而导致静压导轨处的热变形。忽略静压导轨油膜发热，会导致对重型数控落地铣镗床说明书CN104200019A2/7页6热特性的数值仿真与实际不符，影响其温度场和热变形结果的准确性。00092忽略构件结合面换热的热量损失。在重型数控落地铣镗床温度场计算中，装配体的结合面都认为是充分接触没有热量损失的。然而，在亚微观的角度上，任何金属物质的表面都不可能是绝对光滑的。在结合面的热传递方面，实际接触面积仅占。
14、理论接触面积的很小一部分约为1/1000。因为结合面的非充分接触，实际结合面的接触是通过离散点完成的，在结合面会有明显的温度下降。忽略结合面的热量损失，影响温度场的数值分析结果的精度。00103缺乏针对重型数控落地铣镗床整机热特性数值仿真。针对重型数控落地铣镗床，目前常用的研究方法都是针对于机床中关键零部件如主轴部件、滑枕、立柱，但是缺少对重型数控落地铣镗床整机的热特性数值仿真的准确方法。对于整机的热特性的数值仿真，考虑的热源分布众多，模型庞大且复杂，数值仿真的效率，但是整机的热特性数值仿真更符合工程实践的实际状况，数值仿真结果更准确和完整，其分析结果可以作为数控系统的实时热变形补偿产生直接参。
15、考和借鉴作用。发明内容0011本发明的目的在于克服上述缺陷，提供一种分析结果准确性及精度高的重型数控落地铣镗床热特性数值仿真方法。0012本发明针对现有技术的缺陷，主要从以下三方面进行了改进0013一、考虑静压导轨油膜发热，主要针对滑枕、立柱和滑座处装配的静压导轨油膜产热。以上构件在静压导轨上移动会因为油膜粘性阻力产生热量，根据经验公式，可以计算出静压导轨摩擦副在实际工作过程中摩擦功率损失。将静压导轨的油膜产热看作相应部位的热源，进行重型数控落地铣镗床热特性数值仿真，提高数值仿真结果的准确性。0014二、考虑结合面的传热系数，主要针对装配体结合面的非充分接触导致热量损失的状况。影响结合面的接触。
16、热传导过程的因素有很多，通常结合面的传热系数通常由实验确定，进一步的，可根据重型数控落地铣镗床关键零部件的加工工艺和装配过程的公差状况，同时参考工艺系统常用表面间的结合面换热系数取值表，可确定机床中各结合面的换热系数经验取值。换热系数经验取值确定后，可通过有限元软件进行接触热阻的模拟，并在模型接触属性中设置换热系数值，其结果较为真实地反映了模型配合处的热流量传导状况。0015三、以重型数控落地铣镗床整机作为研究对象，主要是将重型数控落地铣镗床作为整体进行热特性数值仿真，而非现有技术中，仅仅单独考虑某个关键零部件的温度场和变形场。通过充分考虑各构件的相互配合和约束，考虑各热源之间的相互影响以及热。
17、量在各构件之间的传递，对于重型数控落地铣镗床热特性的数值仿真指标上考虑更加全面和完整，其温度场和热变形的数值仿真结果准确度更高。0016为了实现上述目的，本发明采用了技术方案如下0017一种重型数控落地铣镗床热特性数值仿真方法，包括以下步骤00181确定重型数控落地铣镗床的热源分布，分析机床与环境之间的热交换；00192在三维建模软件PRO/E中建立重型数控落地铣镗床装配体模型，并对装配体模型进行简化；说明书CN104200019A3/7页700203通过PRO/E和有限元软件WORKBENCH的无缝接口，将简化后的装配体模型导入有限元软件WORKBENCH中，并进行网络划分；00214根据重。
18、型数控落地铣镗床中各个轴承的径向和轴向刚度计算等效圆柱环的等效径向弹性模量和轴向弹性模量，并将等效径向弹性模量和轴向弹性模量导入WORKBENCH中；00225根据重型数控落地铣镗床的材料设置零部件物理性能参数，计算步骤1中所确定的各热源大小和对流传热系数，将计算结果作为边界条件施加在机床上，然后利用有限元软件WORKBENCH对重型数控落地铣镗床分别进行稳态温度场分析和瞬态温度场分析，得出重型数控落地铣镗床温度分布和温升的变化情况，为重型数控落地铣镗床的热位移分析奠定基础；00236根据重型数控落地铣镗床的实际安装位置，在床身底面施加固定约束，将重型机床温度场分析的温度作为载荷，采用间接法，。
19、基于热结构耦合分析重型数控落地铣镗床的变形场，通过在有限元软件WORKBENCH中求解计算，得到重型数控落地铣镗床的热变形场云图。0024其中，所述步骤3中装配体模型的简化内容如下002531去除零件中对分析结果影响不大的孔，去掉螺钉和螺栓孔；002632只考虑机床的热变形，而不考虑结构中热应力影响，故去除大型铸件的倒角和小的阶梯孔，忽略结构中影响小的台阶；002733轴承和滚珠丝杠副结构复杂，工作原理类似，向有限元结构转化时，将轴承和滚珠丝杠副简化成一个相同尺寸的等截面圆柱体；有限元分析中材料属性的弹性模量用等效弹性模量代替，以此降低模型简化对计算精度的影响；002834去除主轴部件内部的拉。
20、刀机构；002935去除机床各进给传动电机。0030考虑到重型数控落地镗铣床整机装配体模型庞大结构复杂，故网络划分采用自由分网的方式，对于关键的零件如主轴，滑枕控制网格大小，提高网格质量。0031所述等效圆柱环的弹性模量的确定方法为0032径向等效弹性模量的确定0033首先，在WORKBENCH结构分析模块中，对轴承模型外圈施加全约束，并在轴与轴承内圈定义接触类型；其次，对轴施加径向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承的刚度变化；再其次，根据刚度计算公式获得轴承有限元计算径向刚度；当轴承有限元计算径向刚度与简化前轴承刚度之差最小时，轴承对应的弹性。
21、模量即为轴承的径向等效弹性模量；0034轴向等效弹性模量的确定0035首先，在WORKBENCH结构分析模块中，对轴承模型的一侧约束自由度，轴承模型的外圈约束出轴向以外的自由度；其次，对轴承模型的另外一侧施加轴向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承轴向的刚度变化；再其次，根据刚度计算公式，获得轴承有限元计算轴向刚度；当滚动轴承有限元计算轴向刚度与简化前轴承刚度之差最小时，轴承对应的弹性模量即为轴承的轴向等效弹性模说明书CN104200019A4/7页8量。0036进一步的，所述步骤5中机床承受的热载荷还包括由于重型数控落地铣镗床各零部件结合面不完。
22、全接触而产生的热量损失。具体的说，通过分析重型数控落地铣镗床的各零件的结合面的加工和装配工艺，确定各结合面的换热系数，然后通过采用在有限元软件中对装配体模型的接触属性设置换热系数的方法仿真热量损失。0037更进一步的，重型数控落地铣镗床零部件之间的接触类型定义如下主轴轴承与主轴及轴套之间的接触为绑定，立柱与滑座之间的接触类型为绑定，其余零部件之间接触类型均为不分离。0038本发明中，所确定的热源分布包括滚动轴承产热、滚珠丝杠副产热、电机产热、静压导轨油膜发热，进一步的，热源和对流传热系数的计算方式如下0039A滚动轴承产热计算0040Q11047MN10041式中N为轴承转速，单位为R/MIN。
23、；M为滚动轴承摩擦力矩，单位为NM，其中，MM0M1，M0反映和轴承所受载荷无关的损失的摩擦力矩；M1反映与轴承所受载荷滑动摩擦相关的损失的摩擦力矩；0042B滚珠丝杠副产热计算0043Q212NM20044式中N为滚珠丝杠转速，单位为R/MIN；M为螺母的摩擦力矩，单位为NM；在实际工程应用中，滚珠丝杠副近一般被看作仅有轴向负载下的向心推力球轴承。对于滚珠丝杠副，螺母的摩擦力矩M由丝杠的驱动力矩MD与滚珠螺旋的阻力矩MP组成，即MMDMP；0045C电机产热计算00460047式中NM是驱动电机的功率，单位为W；是驱动电机的效率；MM为驱动电机输出力矩，单位为NM；N为驱动电机转速，单位为R。
24、/MIN；0048D静压导轨油膜发热计算0049支承滑动面上的剪切应力为0050/H40051支承滑动面上的摩擦力为0052FFAA/H50053摩擦功率损失为0054PFFFA2/H60055式中为油液的动力粘度，单位为PAS；A为密封带部分的支承面积，单位为M2；滑枕、立柱和滑座处的静压导轨油膜产热的大小均由式4至6计算；0056E对流传热系数计算00570058式中NU为努谢尔特准则数；为对流传热系数，单位为W/M2；为流体的导热系数，单位为W/M；D为放热壁面的定型尺寸，单位为M。说明书CN104200019A5/7页90059与现有技术相比，本发明的有益效果在于00601本发明将重型。
25、数控落地铣镗床作为一个整体进行热特性数值仿真，避免了单独针对关键零部件热特性数值仿真的局限性，充分考虑机床复杂的构建配合和热源分布，更准确的预测重型数控落地铣镗床整体的温度及热变形趋势，可为控制重型数控机床的热变形误差比如结构优化改进、关键部位温度控制等可靠的技术数据。00612本发明针对重型数控落地铣镗床温度场和热变形数值仿真时，通过计算滑枕、立柱和滑座处的静压导轨油膜产热，充分考虑重型数控落地铣镗床在工作工程中的热源，减少理论分析中不合理简化热源产生的对重型数控落地铣镗床温度场和热变形数值仿真结果分析的影响，提高了重型数控落地铣镗床热特性数值仿真的准确性。00623本发明通过分析重型数控落。
26、地铣镗床的各零件的结合面的加工和装配工艺，确定各结合面的换热系数；然后，采用在有限元软件中对装配体模型的接触属性设置换热系数的方法，其仿真结果能直接反映模型配合处的热流量传导状况，有效地克服了忽略结合面的接触热阻对重型数控落地铣镗床的热特性数值分析结果不符合工程实际的缺点，进而提高了温度场数值分析结果的精度。附图说明0063图1为TK6916DA热源分布示意图。0064图2为对流换热类型示意图。0065附图中所对应的附图标记名称如下1主轴箱进给驱动系统，2主传动系统，3镗轴进给系统，4滑枕进给系统，5立柱进给系统。具体实施方式0066下面结合实施例对本发明作进一步说明。本发明的实施方式包括但不。
27、限于下列实施例。0067实施例0068本发明提供了一种重型数控落地铣镗床热特性数值仿真方法，该数字仿真方法将机床作为一个整体进行热特性数值仿真，再结合对各零部件结合面的接触热阻和静压导轨油膜的产热，有效地提高了分析结果的准确性和精度。下面以TK6916DA型重型数控落地铣镗床为例，结合该方法进行详解0069首先，分析TK6916DA的热源分布，如图1所示，其热源分布包括主轴轴承摩擦产热、滑座进给传动摩擦产热、立柱静压导轨产热和主轴箱部件产热。分析TK6916DA与环境之间的热交换，如图2所示。0070其次，根据以下公式计算各热源分布的产热0071A滚动轴承产热计算0072Q11047MN100。
28、73B滚珠丝杠副产热计算0074Q212NM20075C电机产热计算说明书CN104200019A6/7页1000760077D静压导轨油膜发热计算0078PFFFA2/H60079E对流传热系数计算00800081根据TK6916DA的二维图纸，在三维建模软件PRO/E中建立其装配体模型，并对装配体模型进行合理的简化；其简化内容包括1、去除零件中的小孔和对分析结果影响不大的孔，去掉螺钉和螺栓孔；2、只考虑机床的热变形，而不考虑结构中热应力影响，故去除大型铸件的倒角，小的阶梯孔，忽略结构中影响小的台阶；3、轴承和滚珠丝杠副结构复杂，工作原理类似，向有限元结构转化时，将轴承和滚珠丝杠副简化成一个。
29、相同尺寸的等截面圆柱体；有限元分析中材料属性的弹性模量用等效弹性模量代替，以此降低模型简化对计算精度的影响；4、去除主轴部件内部的拉刀机构；5、去除机床各进给传动电机。0082通过PRO/E和有限元软件WORKBENCH的无缝接口，将简化后的装配体模型导入有限元软件WORKBENCH中，并进行自由分网。0083然后，根据各个轴承的径向和轴向刚度计算等效圆环的等效径向和轴向弹性模量，并在WORKBENCH中修改各个等效圆环的弹性模量。其中，圆柱环的弹性模量的确定方法为0084径向等效弹性模量的确定0085首先，在WORKBENCH结构分析模块中，对轴承模型外圈施加全约束，并在轴与轴承内圈定义接触。
30、类型；其次，对轴施加径向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承的刚度变化；再其次，根据刚度计算公式，获得轴承有限元计算径向刚度；当轴承有限元计算径向刚度与简化前轴承刚度之差最小时，轴承对应的弹性模量即为轴承的径向等效弹性模量；0086轴向等效弹性模量的确定0087首先，在WORKBENCH结构分析模块中，对轴承模型的一侧约束自由度，轴承模型的外圈约束出轴向以外的自由度；其次，对轴承模型的另外一侧施加轴向集中载荷，以材料参数中的弹性模量为输入参数，轴截面质心的位移为输出参数，进而确定轴承轴向的刚度变化；再其次，根据刚度计算公式，获得轴承有限元计算轴向。
31、刚度；当滚动轴承有限元计算轴向刚度与简化前轴承刚度之差最小时，轴承对应的弹性模量即为轴承的轴向等效弹性模量。0088根据TK6916DA的材料设置零部件物理性能参数，根据计算的热源大小和对流传热系数，确定机床承受的热载荷，求解运算，对TK6916DA分别进行稳态温度场分析和瞬态温度场分析。0089根据TK6916DA的实际安装位置，在床身底面施加固定约束，采用间接法，基于热结构耦合分析TK6916DA的热变形，将热载荷施加在机床上，通过求解计算，可得到机床热变形场云图。0090按照上述实施例，便可很好地实现本发明。值得说明的是，基于上述设计原理的前说明书CN104200019A107/7页11。
32、提下，为解决同样的技术问题，即使在本发明所公开的结构基础上做出的一些无实质性的改动或润色，所采用的技术方案的实质仍然与本发明一样，故其也应当在本发明的保护范围内。0091本发明公开了一种重型数控落地铣镗床热特性数值仿真方法，解决了现有重型数控机床热态特性有限元数值分析方法的分析结果不准确、精度低的问题。该数字仿真方法将机床作为一个整体进行热特性数值仿真，再结合对各零部件结合面的接触热阻和静压导轨油膜的产热，有效地提高了分析结果的准确性和精度。本发明分析结果准确性和精度高，为其大范围的推广应用，奠定了坚实的基础。说明书CN104200019A111/1页12图1图2说明书附图CN104200019A12。