液体静压支承转台系统优化目标函数的响应面拟合.pdfVIP

北京力学会第18届学术年会论文集：计算力学 液体静压支承转台系统优化目标函数的响应面拟合· 文聘 叶红玲 沈静娴 (北京工业大学机电学院100124) 摘要：本文对液体静压转台系统固有振动频率进行了有限元分析，采用最大差值极小化方 法建立频率与油腔外半径和初始油膜厚度之间的函数关系，为进一步对转台进行动态优化 奠定了理论基础。 关键词：液体静压转台，频率，最大差值最小化，响应面拟合 一、引言 响应面方法是基于一组评估函数的实验点采用统计学的综合实验技术解决复杂系 统输入输出关系的方法。本文采用的最大差值极小化拟合响应面函数方法借助精确映射 反演方法【l】，将样本点的一次响应面函数和二次带交叉项响应面函数的拟合问题转换为 求解一类线性规划问题。变量试验点的选取采用中心扩展设计的原则。为了使响应面方 法拟合出的方程与实际物理量之间的误差在可以接受的范围内，选取合适的置信区间半 径，并对响应逼近函数的预估性能进行评价。本文对液体静压支承系统中两个设计变量 与优化目标函数进行了拟合。 二、液体静压转台系统优化目标函数 液体静压支承转台系统振动频率的大小及其振动形式不仅影响液体静压支承转台 系统的工作性能，而且直接影响到数控机床的加工精度。因此优化目标是使系统低阶频 率落在最优Ⅸ域内。一般运动部件的低阶频率要避免靠近相应工作转速，部件才能更好 的工作f2J，此时低阶频率即落在最优区域内，能有效的防止共振。 系统P阶固有频率为f(f-l…p)，部件工作转速为玎，其多倍转速分别为 吩=ixn(i=2…p+1)，令q=(乃+％)／幻=1··一，则 F1=(1-砰f1)2、F2=(1-瓦f2。)2,oo鼻=(1一虿ft2)2，㈨”p。 1L f2 (1) 则目标函数可以表示为：F=∑(1一箫)2一min(i=l…p) 在工程中液体静压转台振动检测中通常取前四阶频率值，因此设定P=4。目标函 数与设计变量(油腔外半径、初始油膜厚度)关系是隐式的，因此需要进行响应面拟合。 ‘北京市教委项目(蹦2009lO005005)，北京市人才强教计划中青年骨干项目(PHR20l1080)支持． VI．109 北京力学会第18届学术年会论文集：计算力学 二次响应面函数拟合 3．1目标函数拟合 我们选择转台外半径、初始油膜厚度为设计变量，采用二次最大差值极小化拟合响 应面函数方法将目标函数式(1)显式化，形成的显式目标函数为 (2) F(五功=ai+呸五+a3x2+％彳+呜《+％．_恐=ai+呸R+呸％+％露+呜聍+％冗％ 采用中心扩展设计的原则，初始中心点选取(R‘叭，‰‘o’)，尺(o)：0．03525， 对5个试验点进行数值模拟计算，得到这五个试验点的数值如表l所示。 表l初始试验点数值 根据试验点的相应点数值，利用最大差值极小化拟合响应面方法，对上述试验点及 其响应值进行二次响应面拟合，得到函数关系为 (3) 3．2目标函数响应面拟合精度检验 选取设计变量的初始取值范围的边界点进行特殊点检验，从表2可以看出：检验点 的拟合值与数值模拟值之间的误差很小，拟合精度为0．0476；这些数据表明由最大差值 最小化拟合二次响应面得到的目标函数表达式非常接近于真实的目标函数表达式，拟合 程度较高，因此采用最大差值极小化拟合二次响应面法拟合目标函数是准确可靠的。 表2响应面拟合精度 捡熊盛 塑金鱼塾笪 查匿丞送型焦