本发明专利技术公开了一种基于灵敏度分析的定量式静压转台优化设计方法，该方法建立静压转台台面的参数化模型，根据流体力学相关理论计算转台油膜的刚度、阻尼及泵功率；使用Ansys计算转台台面各设计参数对于转台台面的综合变形、质量、及固有频率的灵敏度并找到关键设计参数。使用Matlab计算各设计参数对刚度、阻尼及泵功率的灵敏度并找到关键设计参数；使用Isight连接ansys,catia及matlab，选用粒子群算法对找出来的关键设计参数进行重点优化；通过本发明专利技术方法的优化，对于提升转台的设计速度、设计质量有重要作用。
【技术实现步骤摘要】
一种基于灵敏度分析的定量式静压转台优化设计方法
本专利技术属于静压转台设计领域，涉及一种双圈油垫支承的定量式静压转台的优化设计方法，更具体是一种基于灵敏度分析的定量式静压转台优化设计方法，该方法首先运用灵敏度分析找出对转台各项性能影响最大的设计参数，而后对其进行有针对性的优化。
技术介绍
静压转台(HydrostaticRotaryTable)用有压力的流体使有相对运动的两个表面分开并借助流体静压来承载。由于运动副之间完全被油膜隔开，所以运动副间的摩擦力大大减小，同时其承载能力、运动精度与寿命却大大提高。正因为液体静压支承的诸多优点，所以它广泛的应用于重型机床并成为其关键部件之一。静压转台的各项性能对于机床的精度及加工质量都有重要的影响，找到转台的最优设计参数对于全面提升转台的性能有重要作用。
技术实现思路
本专利技术的目的是提供一种基于灵敏度分析的定量式静压转台优化设计方法，首先建立转台静压部分及结构部分的性能模型，之后对设计参数进行灵敏度分析找出影响最大的设计参数，最后用Isight软件对转台进行优化设计。为实现上述目的，本专利技术是采用以下技术手段实现的：1、建立静压转台台面的参数化模型，根据流体力学相关理论计算转台油膜的刚度、阻尼及泵功率；2、使用Ansys计算转台台面各设计参数对于转台台面的综合变形、质量、及固有频率的灵敏度并找到关键设计参数。使用Matlab计算各设计参数对刚度、阻尼及泵功率的灵敏度并找到关键设计参数；3、使用Isight连接ansys,catia及matlab，选用粒子群算法对找出来的关键设计参数进行重点优化；本专利技术的特点在于基于灵敏度分析及优化理论，建立了转台台面及静压部分的优化目标函数并根据粒子群算法对转台的性能进行优化，
技术实现思路
包括三部分，在第一部分中建立静压转台优化设计的优化目标函数，在第二部分中计算各设计参数的局部灵敏度并加以分析找到最关键的设计参数；在第三部分中对关键参数进行优化设计。最后通过实例论证本专利技术所述方法的有效新性与可行性。与现有技术相比，本专利技术具有如下有益效果。静压转台的设计参数众多它们对于静压转台性能的影响也各不相同，在进行转台的优化设计时如果对每一个设计参数都进行优化那工作量将十分巨大，需要对众多设计参数进行筛选找到对转台性能影响最明显的参数，然后有针对性的进行优化，这对于提升转台的设计速度，设计质量有重要作用。附图说明图1静压转台结构简图。图2静压转台台面结构简图。图3转台台面内部结构图。图4转台台面三维模型图。图5支承油垫结构简图。图6预压油垫结构简图。图7优化流程图。图8各个设计参数的灵敏度。具体实施方式本专利技术实施一种双圈油垫支承的静压转台的优化设计方法，下面结合附图，对本专利技术的实施进行具体说明。图1为双圈支承的静压转台台面结构简图，转台由转台台面、基座、支承油垫、预压油垫组成，转台的自重为G，各油垫均由定量泵供油其中支承油垫的供油量为Q0，预压油垫的供油量为Q1。图2、图3、图4分别为转台台面的结构简图及三维图，图5图6为支承油垫与预压油垫的结构简图，预压油垫是环形油垫而支承油垫是圆形油垫。一种基于灵敏度分析的定量式静压转台优化设计方法，其特征在于：该静压转台由转台台面、基座、支承油垫、预压油垫组成，转台的自重为G，各油垫均由定量泵供油其中支承油垫的供油量为Q0，预压油垫的供油量为Q1；该设计方法流程如下，步骤1目标函数的建立1.1根据静压转台台面结构简图绘制转台台面的三维参数化模型，完成建模后的转台的质量、固有频率以及最大设计载荷作用下的变形量，这些性能与每一个设计参数都有关系，于是性能与设计参数的关系写成函数如下所示：freq＝fr(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(1)mass＝fm(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(2)deform＝fd(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(3)上式中freq表示转台台面的固有频率，mass表示转台台面的质量，deform表示转台台面在一固定位置施加固定载荷时转台的最大变形量；而a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b都为转台台面的设计参数，那么转台结构部分的目标函数写为f1＝min(freq,-mass,deform)(4)1.2根据量纲分析将连续性方程及N-S方程简化并求解得到一维雷诺方程和流出封油边的油液的流量为：1.2.1支承油垫承载力计算支承油垫的结构简为圆形阶梯结构，R1为油垫内径R2为油垫的外径hi为油垫的油膜厚度，则对于圆形定量补偿的油垫有边界条件：将边界条件(7)带入方程(5)和(6)解得各支承油垫的油腔压力p0i和压力分布pi(r)进而得到各支承油垫的承载力Fi为：于是支承油垫的刚度阻尼与泵功率分别为：其中KSi为支承油垫的刚度，CSi为支承油垫的阻尼，NTi为支承油垫的泵功率；1.2.2预压油垫承载力的计算预压油垫为环形油垫，对于环形油垫有边界条件为：将(12)式代如(5)、(6)式得到预压油垫的油腔压力p0y封油边压力分布p1y(r)、p2y(r)为：进而计算出预压油垫的承载力Fy为：所以预压油垫的刚度Ky阻尼Cy与泵功率NTy为：1.2.3转台的整体刚度阻尼与泵功率的计算转台的支承系统由各支承油垫与预压油垫并联而成那么转台的整体刚度、阻尼与泵功率为：上式中KZ为转台的轴向刚度，Kt为倾覆刚度，CZ为轴向阻尼，Ct为倾覆阻尼，NT为总泵功率，n1为第一圈支承油垫的个数，n2为第二圈支承油垫的个数根据经验n2＝2n1；所以静压部分的目标函数写为:f2＝min[KZ,Kt,CZ,Ct,-NT]T(19)所以将静压部分的目标函数与结构部分的目标函数组合起来则转台的整体目标函数为f＝min(f1|f2)(20)步骤2灵敏度计算转台的各设计参数对结构特性影响程度不一；对转台进行优化设计时需要找出对结构影响程度小的参数；而灵敏度是关注指标对某设计参数的变化梯度，灵敏度的计算如下式所示，对转台进行灵敏度分析帮助找出对结构影响程度小的参数；由上式计算出各设计参数对于各设计性能的灵敏度值，由于静压部分的灵敏度值差别太大所以需要对其进行取对数处理而结构部分的灵敏度值不需要进行这样的处理，对于每一项性能从各参数的灵敏度值中挑出灵敏度值最大一项，然后其余各参数的灵敏度值均除以最大的灵敏度值，这样就对各参数的灵敏度值进行了无量纲化处理；之后对于结构部分选取无量纲灵敏度值大于0.60或小于-0.6的参数作为关键设计参数，而对于静压部分要选择无量纲灵敏度值大于0.7或小于-0.7的参数作为关键设计参数，每一设计参数只要对一项性能的灵敏度值满足上述条件就认为此参数为关键参数；之后将对这些关键参数进行重点优化设计；步骤3优化设计用Isight连接Catia，Ansys和Matlab，设置由灵敏度分析得到的关键参数为被优化的参数，给出这些关键参数的取值范围，各参数的取值范围统一设置为在各参数初始值左右变动，变动量为初始值的一个固定百分比，取20％；而后设置目标函数及限制条件，目标函数就是使固有频率、刚度、阻尼最大而转台重量、变形及泵功率最小，限制条件为油垫的内径小于外径；本文档来自技高网...
【技术保护点】
一种基于灵敏度分析的定量式静压转台优化设计方法，其特征在于：该静压转台由转台台面、基座、支承油垫、预压油垫组成，转台的自重为G，各油垫均由定量泵供油其中支承油垫的供油量为Q0，预压油垫的供油量为Q1；该设计方法流程如下，步骤1目标函数的建立1.1根据静压转台台面结构简图绘制转台台面的三维参数化模型，完成建模后的转台的质量、固有频率以及最大设计载荷作用下的变形量，这些性能与每一个设计参数都有关系，于是性能与设计参数的关系可以写成函数如下所示：freq＝fr(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b) (1)mass＝fm(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b) (2)deform＝fd(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b) (3)上式中freq表示转台台面的固有频率，mass表示转台台面的质量，deform表示转台台面在一固定位置施加固定载荷时转台的最大变形量；而a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b都为转台台面的设计参数，那么转台结构部分的目标函数可写为f1＝min(freq,‑mass,deform) (4)1.2根据量纲分析将连续性方程及N‑S方程简化并求解的到一维雷诺方程和流出封油边的油液的流量为：1r∂∂r(rh312η∂p∂r)=∂h∂t---(5)]]>Q(r)=-πrh36η∂p∂r---(6)]]>1.2.1支承油垫承载力计算支承油垫的结构简为圆形阶梯结构主要尺寸在图中标明，R1为油垫内径R2为油垫的外径hi为油垫的油膜厚度，则对于圆形定量补偿的油垫有边界条件：将边界条件(7)带入方程(5)和(6)可以解得各支承油垫的油腔压力p0i和压力分布pi(r)p0i=6ηln(R2R1)πhi3(Q0+πR12∂hi∂t-π(R22-R12)2ln(R2R1)∂hi∂t)---(8)]]>pi(r)=r23ηπhi3∂hi∂t-(p0iln(rR2)+3η(R22ln(rR1)-R12ln(rR2))hi3∂hi∂t)ln(R2R1)---(9)]]>进而可以得到各支承油垫的承载力Fi为：Fi=πR12p0i+2π∫R1R2rpi(r)dr---(10)]]>于是支承油垫的刚度阻尼与泵功率分别为：KSi=9Q0η(R22-R12)hi4,CSi=3πη(R22-R12)22hi3,NTi=p0iQ0=6ηln(R2R1)Q02πhi3---(11)]]>其中KSi为支承油垫的刚度，CSi为支承油垫的阻尼，NTi为支承油垫的泵功率；1.2.2预压油垫承载力的计算预压油垫为环形油垫，对于环形油垫有边界条件为：将(12)式代如(5)、(6)式可得到预压油垫的油腔压力p0y封油边压力分布p1y(r)、p2y(r)为：p0y=6ηln(RC4RC3)ln(RC2RC1)πhy3ln(RC4RC2RC3RC1)(Q1+π(RC32-RC22)∂hy∂t-π(ln(RC2RC1)(RC42-RC32)+ln(RC4RC3)(RC12-RC22))2ln(RC2RC1)ln(RC4RC3)∂hy∂t)---(13)]]>p1y(r)=r23ηπhy3∂hy∂t-(p0yln(rRC1)+3η(RC12ln(rRC2)-RC22ln(rRC1))hy3∂hy∂t)ln(RC1RC2)---(14)]]>P2y(r)=r23ηπhy3∂hy∂t-(p0yln(rRC4)+3η(RC42ln(rRC3)-RC32ln(rRC4))hy3∂hy∂t)ln(RC4RC3)---(15)]]>进而可以计算出预压油垫的承载力Fy为：Fy=π(RC32-RC22)p0y+2π∫Rc1Rc2rp1y(r)dr+2π∫Rc3Rc3rp2y(r)dr---(16)]]>所以预压油垫的刚...
【技术特征摘要】
1.一种基于灵敏度分析的定量式静压转台优化设计方法，其特征在于：该静压转台由转台台面、基座、支承油垫、预压油垫组成，转台的自重为G，各油垫均由定量泵供油其中支承油垫的供油量为Q0，预压油垫的供油量为Q1；该设计方法流程如下，步骤1目标函数的建立1.1根据静压转台台面结构简图绘制转台台面的三维参数化模型，完成建模后的转台的质量、固有频率以及最大设计载荷作用下的变形量，这些性能与每一个设计参数都有关系，于是性能与设计参数的关系写成函数如下所示：freq＝fr(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(1)mass＝fm(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(2)deform＝fd(a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b)(3)上式中freq表示转台台面的固有频率，mass表示转台台面的质量，deform表示转台台面在一固定位置施加固定载荷时转台的最大变形量；而a1,a2,a3,a4,a5,a6,a7,RL,RS,Dy,Dr,b都为转台台面的设计参数，那么转台结构部分的目标函数写为f1＝min(freq,-mass,deform)(4)1.2根据量纲分析将连续性方程及N-S方程简化并求解得到一维雷诺方程和流出封油边的油液的流量为：1.2.1支承油垫承载力计算支承油垫的结构简为圆形阶梯结构，R1为油垫内径R2为油垫的外径hi为油垫的油膜厚度，则对于圆形定量补偿的油垫有边界条件：将边界条件(7)带入方程(5)和(6)解得各支承油垫的油腔压力p0i和压力分布pi(r)进而得到各支承油垫的承载力Fi为：于是支承油垫的刚度阻尼与泵功率分别为：其中KSi为支承油垫的刚度，CSi为支承油垫的阻尼，NTi为支承油垫的泵功率；1.2.2预压油垫承载力的计算预压油垫为环形油垫，对于环形油垫有边界条件为：将(12)式代如(5)、(6)式得到预压油垫的油腔压力p0y封油边压力分布p1y(r)、p2y(r)为：
【专利技术属性】
技术研发人员：刘志峰，湛承鹏，李小燕，罗兵，张敬莹，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。