发明名称 基于Ansys Workbench的潜油电机温度场模拟方法
摘要 基于Ansys Workbench的潜油电机温度场模拟方法,属于潜油电机温度场模拟技术领域。本发明是为了解决现有潜油电机的温度场分析方法不能对潜油电机进行实时监控,以进行过热保护的问题。它在SolidWorks系统中对潜油电机的定子和转子进行三维建模,获得潜油电机模型;将潜油电机模型导入到Ansys Workbench系统中,进行自动网格剖分;采用潜油电机的热力学微分方程和热力学微分方程边界条件,确定潜油电机模型求解域的温度场;进行潜油电机仿真,获得潜油电机在不同负载下定子、转子和导条的温度场分布云图。本发明用于潜油电机温度场的仿真。
申请公布号 CN103246764B 申请公布日期 2016.02.03
申请号 CN201310138051.6 申请日期 2013.04.19
申请人 哈尔滨工业大学 发明人 王立国;徐殿国;安天琪;韩宇泽;吴松霖;胡东;张世博
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 一种基于Ansys Workbench的潜油电机温度场模拟方法,其特征在于,它包括以下步骤:步骤一:在SolidWorks系统中对潜油电机的定子和转子进行三维建模,获得潜油电机模型;步骤二:将潜油电机模型导入到Ansys Workbench系统中,采用Ansys Workbench系统对潜油电机模型进行自动网格剖分;步骤三:对自动网格剖分后的潜油电机模型的温度场进行分析,采用潜油电机的热力学微分方程和热力学微分方程边界条件,确定潜油电机模型求解域的温度场,该温度场为由对流散热产生的温度场;步骤四:根据潜油电机的热力学微分方程和热力学微分方程边界条件设置潜油电机模型的模型材料、初始边界条件和负载,进行潜油电机仿真,获得潜油电机在不同负载下定子、转子和导条的温度场分布云图;所述步骤三中采用潜油电机的热力学微分方程和热力学微分方程边界条件,确定潜油电机模型求解域的温度场的具体方法为:所述潜油电机的热力学微分方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><mi>T</mi></mrow><mrow><mo>&part;</mo><mi>&tau;</mi></mrow></mfrac><mo>=</mo><mi>a</mi><mrow><mo>(</mo><mfrac><mrow><msup><mo>&part;</mo><mn>2</mn></msup><mi>T</mi></mrow><mrow><mo>&part;</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&part;</mo><mn>2</mn></msup><mi>T</mi></mrow><mrow><mo>&part;</mo><msup><mi>y</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&part;</mo><mn>2</mn></msup><mi>T</mi></mrow><mrow><mo>&part;</mo><msup><mi>z</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>q</mi><mi>v</mi></msub><mrow><mi>&rho;</mi><mi>c</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000011.GIF" wi="692" he="180" /></maths>式中T为温度,τ为时间微元,a为散热系数,x为空间坐标系中X轴方向的坐标值,y为空间坐标系中Y轴方向的坐标值,z为空间坐标系中Z轴方向的坐标值,q<sub>v</sub>为单位体积生热率,ρ为密度,c为质量定容热容,其中<img file="FDA0000812440170000012.GIF" wi="186" he="148" />式中λ为导热系数;热力学微分方程的热力学第一边界条件为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>T</mi><msub><mo>|</mo><msub><mi>S</mi><mn>1</mn></msub></msub><mo>=</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><mi>T</mi><msub><mo>|</mo><msub><mi>S</mi><mn>1</mn></msub></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000013.GIF" wi="730" he="124" /></maths>式中,T<sub>0</sub>为已知温度;f(x,y,z,t)为已知温度函数,t为时间,S<sub>1</sub>为热力学第一类边界面;热力学第一边界条件用于潜油电机模型的热力学第一类边界面S<sub>1</sub>上温度为已知的情况;热力学微分方程的热力学第二边界条件为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>-</mo><mi>&lambda;</mi><mfrac><mrow><mo>&part;</mo><mi>T</mi></mrow><mrow><mo>&part;</mo><mi>n</mi></mrow></mfrac><msub><mo>|</mo><msub><mi>S</mi><mn>2</mn></msub></msub><mo>=</mo><mi>q</mi><mo>,</mo><mo>-</mo><mi>&lambda;</mi><mfrac><mrow><mo>&part;</mo><mi>T</mi></mrow><mrow><mo>&part;</mo><mi>n</mi></mrow></mfrac><msub><mo>|</mo><msub><mi>S</mi><mn>2</mn></msub></msub><mo>=</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000021.GIF" wi="808" he="170" /></maths>式中n为微元点,S<sub>2</sub>为热力学第二类边界面,q为热流密度,g(x,y,z,t)为热流密度函数;热力学第二边界条件用于热力学第二类边界面S<sub>2</sub>上的热流密度为已知的情况;热力学微分方程的热力学第三边界条件为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mo>-</mo><mi>&lambda;</mi><mfrac><mrow><mo>&part;</mo><mi>T</mi></mrow><mrow><mo>&part;</mo><mi>n</mi></mrow></mfrac><msub><mo>|</mo><msub><mi>S</mi><mn>3</mn></msub></msub><mo>=</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000022.GIF" wi="580" he="201" /></maths>式中S<sub>3</sub>为热力学第三类边界面,T<sub>f</sub>为热力学第三类边界面S<sub>3</sub>周围流体的温度,根据上述潜油电机的热力学微分方程和热力学微分方程边界条件,确定潜油电机模型求解域的温度场;所述步骤四中获得潜油电机在不同负载下定子、转子和导条的温度场分布云图的具体方法为:根据潜油电机的热力学微分方程和热力学微分方程边界条件确定潜油电机单位体积生热率q<sub>v</sub>;然后根据潜油电机的定子阻值、转子阻值、额定工作电流和额定工作电压确定转子绕组的铜损p<sub>r_copper</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>r</mi><mo>_</mo><mi>c</mi><mi>o</mi><mi>p</mi><mi>p</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>=</mo><mn>3</mn><msup><msub><mi>I</mi><mi>s</mi></msub><mn>2</mn></msup><mfrac><mrow><msup><mi>s</mi><mn>2</mn></msup><msup><msub><mi>&omega;</mi><mi>e</mi></msub><mn>2</mn></msup><msup><msub><mi>L</mi><mi>m</mi></msub><mn>2</mn></msup></mrow><mrow><msup><msub><mi>R</mi><mi>r</mi></msub><mn>2</mn></msup><mo>+</mo><msup><mi>s</mi><mn>2</mn></msup><msup><msub><mi>&omega;</mi><mi>e</mi></msub><mn>2</mn></msup><msup><msub><mi>L</mi><mi>r</mi></msub><mn>2</mn></msup></mrow></mfrac><msub><mi>R</mi><mi>r</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000023.GIF" wi="887" he="217" /></maths>式中I<sub>s</sub>为定子线圈的电流有效值,<sub>s</sub>为转差率,ω<sub>e</sub>为同步角速度,L<sub>m</sub>为定子与转子之间的互感,R<sub>r</sub>为转子阻值,L<sub>r</sub>为转子自感;根据定子线圈的电流有效值及定子线圈的阻值,定子铜损P<sub>s_copper</sub>为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>s</mi><mo>_</mo><mi>c</mi><mi>o</mi><mi>p</mi><mi>p</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>=</mo><mn>3</mn><msubsup><mi>I</mi><mi>s</mi><mn>2</mn></msubsup><msub><mi>R</mi><mi>s</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000812440170000024.GIF" wi="453" he="122" /></maths>R<sub>s</sub>为定子线圈的阻值;利用热力学微分方程的热力学第二边界条件求解潜油电机相互接触的各个部分之间由于热传导和对流换热产生的温度场;利用热力学微分方程的热力学第二边界条件求解潜油电机内部由于对流换热产生的温度场,该对流换热发生在两部分,一为定子内表面及转子外表面与气隙空气的对流换热;二为机壳外表面与原油的对流换热;通过设定已知温度T<sub>0</sub>、热力学第三类边界面S<sub>3</sub>周围流体的温度T<sub>f</sub>、导热系数λ、密度ρ和质量定容热容c,确定热力学微分方程的热力学第三边界条件,再对潜没电机的每个微元求解热力学微分方程,获得潜油电机在不同负载下定子、转子和导条的温度场分布云图。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号