垂直微结构剖面仪主体结构优化.docx

1、垂直微结构剖面仪主体结构优化姜欣刘玉红宋诗军王延辉天津大学，天津,300072摘要：垂直微结构剖面仪是获取海洋微结构湍流数据最有效的测量平台，其下降过程中产生的流致振动是影响剖面仪测量材度的主要因素。分析了影响剖面仪流致振动的因素，基于ANSYS-Workbench数值仿真平台，采用双向流固耦合算法，得到了割面仪流致振动与其结构外形及参数之间的美系。在此基础上，以减,小流致振动为目标，采用单目标优化方法分别对剖面仪主体结构外形及结构参数进行了优化.优化后的剖面仪测量耕度提高了57%。研究成果为垂直微结构割面仪及类似站构的合理设计提供了理论指导。关键词：垂直微结构剖面仪；测址精度；流致振动；流固

2、耦合；结构优化中图分类号:THU3DOI：10.3969/j.issn.1004-132X.2015.01.017StructureOptimizationforMainBodyofVerticalMicrostructureProfilerJiangXinLiuYuhongSongShijunWangYanhuiTianjinUniversity,Tianjin,300072Abstract：Verticalmicrostructureprofiler(VMP)isamostcommonandeffectiveplatformtoobtaindataofmicrostructureturbu

3、lenceintheocean.TheFIVisthemainfactoraffectingthemeasurementprecisionoftheVMP.FactorsinfluencingtheFIVoftheVMPwereanalyzedherein.BasedonthenumericalsimulationplatformofANSYS-Workbench,theFIVoftheVMPwithvariousshapesandstructureparameterswascalculatedusingalgorithmoftwo-wayfluid-structurecoupling.Aim

4、edtodecreasetheFIV,themainbodyshapeandthestructureparametersofthemainbodywereoptimizedrespectively.Themeasurementprecisionoftheoptimizedprofilerisimprovedby57%.ThepresentinvestigationprovidesscientificguidancefordesigningtheVMP.Keywords：verticalmicrostructureprofilerVMP)；measurementprecision；flow-in

5、ducedvibra-tion(FIV)?fluid-structurecoupling；structureoptimization0引言海洋微结构湍流对海水的动披、热量和质量输运都有重要贡献，对海水运动速度、温盐特性及水中溶解态、颗粒态物质的分布有显著影响，是产生海洋宏观现象的原动力。搭载翼型剪切流传收稿日期；2013-03-21基金项目，国家高技术研究发腰计划(863计划)资助项目(2010AA09Z102)i国家自然科学基金资助项目(51105268)YangWangyucHuAnmin.SunZuqing.EffectofAusteniteGrainSizeonStrainEnhanc

6、edTransformationinaLowCarbonSteelJJ.ActaMetallurgiesSinica.200036(10):1055-1060.12尹云洋.杨王珂，李龙飞，等.基于过冷奥氏体动态相变的热轧TRIP钢组织控制I.原始奥氏体品粒尺寸J.金属学报，2010,46(2),155-160.YinYunyang.YangWangyue,LiLongfei*etal.MicrostruciureControlofHotRolledTripSteelBasedonDynamicTransformationofUndercooledAusteniteI.PriorAustenit

7、eGrainSizeJ.Acta感器的垂直微结构剖面仪是精确测量海洋微结构湍流数据的最有效和普遍的测量平台。垂直微结构剖面仪（以下简称剖面仪）的研制始于20世纪70年代，由Osborn研制成功，并搭载了翼型剪切流传感器，成功地对海洋微结构湍流进行了测虽。目前，加拿大、美国、日本及欧洲部分国家对剖面仪的研制处于领先水平J国内仅天津大学最早于2004年开始进行剖面仪MetallurgiesSinica,2010.46(2)：155-160.(编辑陈勇)作者简介：以男.1978年生.燕山大学机械工程学院期教授.主要研究方向为大型锻件工艺及优化.畋，女.1988年生.燕山大学机械工程学院硕士研究生。金

8、森，界.1968年生.燕山大学机械工程学院教授、博士研究生导师。!玉邺.男，1978年生.燕山大学国家冷轧板带装备及工艺工程技术研究中心实验师.簇潸华.舅.1963年生.中国第二重型机械集团公司高级工程师，的研制刃，历经原理样机、试验样机，于2011年研制成功我国第一台剖面仪工程样机目前剖面仪不仅可以测量:海洋微结构湍流，而且可以搭载CTD、浊度计等进行海洋温度、盐度、浊度等海洋环境参数的测量.剖面仪下降过程中产生的流致振动是影响其测量精度的主要因素。流致振动现象是剖面仪下降过程中不可避免的一种有害现象，会影响剖面仪的运动稳定性从而影响湍流测量数据的准确性。为了减小剖面仪下降过程中产生的流致振

9、动影响,本文以天津大学研制的第二代剖面仪工程样机为对象，分析影响剖面仪流致振动的结构因素，以减小流致振动为目标，对剖面仪外形及结构参数进行优化设计。1剖面仪的流致振动及计算方法1.1剖面仪的流致振动剖面仪在下降过程中受到流体力的作用，会改变剖面仪的运动状态，剖面仪运动状态的改变反过来乂会影响剖面仪周围的流场，从而使剖面仪受到的流体力发生改变。剖面仪与流体的相互作用最终会达到稳定，使剖面仪呈现出规律的摆动，这就是剖面仪的流致振动。Mourn等对Chameleon剖面仪的实验数据进行分析后指出，流致振动决定着剖面仪湍动能耗散率的测量下限。Osborn在使用第一台垂直剖面仪测筮时发现，低频率的凹切流

10、数据不是实际的流速脉动，而是剖面仪本体由于流致振动而产生的低频钟摆运动。Miller等在AMP系列剖面仪研究过程中发现，当剖面仪的流致振动频率接近固有频率时会引起剖面仪的共振，严重影响剖面仪的测量精度，他还指出平滑的剖面仪秃体（如减少凸起）可改善剖面仪的流致振动情况。以上研究结果表明，剖面仪的流致振动与其结构密切相关.因此，可通过优化剖面仪的结构形状及结构参数来降低剖面仪的流致振动效应。1.2流致振动计算方法基于ANSYS-Workbench数值仿真平台，将流体计算软件CFX与结构计算软件ANSYS结合进行流固双向耦合计算，由CFX进行流场计算，ANSYS进行结构计算。采用双向分离式顺序计算方

11、法.即流固耦合系统的流体域和固体域按各自划分的网格与设定的参数分别进行计算，然后在流固耦合面上进行计算数据的传输，实现耦合计算，其中CFX向ANSYS传输的是剖面仪的受力情况，ANSYS向CFX传输的是剖面仪的位移情况。文献13证明了利用该方法进行剖面仪流致振动计算的准确性与可靠性。流固耦合计算区域计算区域为长方体，如图1所示。为保证流体充分发展，剖面仪前端距流体入曰的距离L】=2L,剖面仪后端距流体出口的距离L2=3L,其中L为剖面仪主体氏度（图1中点A、O之间距离）；由于剖面仪流致振动的计算涉及动网格，剖面仪在横向进行摆动，所以水域的宽度W要取得足够大，在此取H=W=16D,其中D为剖面仪

12、的直径。图1中O点为摆动中心的位曾,F,为剖面仪浮心的位置,G,为剖面仪重心的位置,A点为观测点。计算时，记录观测点A的位移与速度，从而获得剖面仪的流致振动数据。图I剖面仪计算区域示意图流体区域与结构区域均采用四面体进行网格划分。为了保证计算精度，经网格独立性检验，流场区域剖面仪最大面网格设为16mm,增K率设为1.1,结构区域剖面仪面网格最大尺寸为30mm。流体区域的入口设为速度入口，流速取剖面仪实际工作速度0.65m/s,流体区域出口设为压力出口，流体四周的区域边界设为对称边界条件，剖面仪的壁面设置为流固耦合面。流场计算采用k-e湍流模型，时间步长为0.1s,总计算时间为80s。计算时，记

13、录观测点A处的速度及位移变化曲线。2剖面仪结构外形优化最初OsbornF的垂直剖面仪采用的是圆球状.但在海试试验中发现圆球状的剖面仪在下降过程中会产生较大的旋转和振动，从而产生较大的测量误差。经过诸多研究者的不断研究、试骚与改进，发现长圆柱体的外形相较于其他形状具有更好的下降稳定性，而且下降过程中不易倾斜和偏转，所以如今的剖面仪均采用长圆柱的主体外形。虽然整体的形状可以确定为长圆柱形，但局部形状的不同也会导致流致振动的不同。因此，本文针对天津大学研制的第二代长圆柱体形剖面仪.对M局部结构（头部及主体）形状进行优化.2.1头部形状优化天津大学研制的第二代剖面仪的原始头部端面设计为一平面，如图2a

14、所示，这种结构易使剖面仪下降过程中在头部流场产生边界层分离，从而增大流致振动效应。为改善头部流场悄况，将剖面仪头部设计为半椭球形，如图2b所示.可以减小头部末端的边界层分离.为优化头部椭球形状.令椭球头部轴线方向尺寸L在60120mm之间变化，径向尺寸保持不变（为剖面仪主体宜径尺寸160mm）,分别建立模型进行流致振动计算。计算结果如图3所示。（苗惊始头部形状（b）m球形头部图2剖面仪头部形状图3剖面仪头部椭球轴向尺寸对流致振动的影响图3中，左侧纵轴为剖面仪前端监测点A的y向加速度，右侧纵轴为剖面仪流致振动频率。由图3知，剖面仪头部形状尺寸对流致振动的频率影响很小，频率在0.16Hz左右变化；

15、桶球体头部轴向尺寸的增大会使流致振动的加速度幅值减小，两者近似成指数关系；当轴向尺寸大于80mm（即半球形）时.随柠尺寸的增大，流致振动的减小幅度较小、程度变级。同时考虑到安装传感器的需求，剖面仪头部前端曲率不能过大因此剖面仪前端设计成半球形较为合适.2.2主体结构外形优化原始设计中，剖面仪前端具有凹阶结构如图2所示。研究发现.剖面仪下降时在凹阶处会发生边界层分离现象.如图4所示。由图4以看出.凹阶处发生了涡的脱落.并且脱落的涡又会在后部再次附若.造成后面涡结构的混乱.导致流致振动程度的增强。因此，对有无凹阶结构的剖面仪进行流致振动对比计算，得到加速度的幅值和振动频率，见&1。由表】W以石出.

16、去掉凹阶结构使剖面仪的流致振动加速度幅值减小r31.17%.明显减缓流致振动幅度.因此在设计中应尽ht避免凹阶结构的出现。图4凹阶处涡脱落图表I有无凹阶对流致振动的影响主体外形加速度帼fiVmm/s：）顿率（Hz）无凹阶2.650.147有凹阶3.850.1423剖面仪结构参数优化由上文分析知，剖面仪的最优外形结构如图5所示.头部为半球形，主体为长圆柱.而IL无凹阶结构。除外形影响流致振动外.剖血仪的结构参数（如长径比和重浮心距离）对流致振动也具有取要影响,同样需要:优化以减小剖面仪的流致振动。津=4图5割面仪主体外形示意图3.1长径比按图5所示外形建立流固糟合计算模型。设计时.当剖面仪质址确

17、定后.剖面仪的体积也为一定值，因此.计算中改变长径比时应保uE体积不变。令剖面仪长度L分别为】.72.1m,对应算出剖血仪的住径D，从而得到剖面仪的K径比L/D.图6所示为监测点A的y向加速度和振动频率随长径比的变化曲线。由图6知.随着长径比的增大.流致振动的频率越来越低；削面仪的振动加速度幅值与长径比成作线性关系随齐长径比的增大，剖面仪流致振动幅值越来越小。根据计算结果，流致振动的频率，与K径比成线性关系.其关系式可表述为1.59.50.1750.1700.1650.1600.1550.1500.14511.012.5长径比与55555432(S.0.14014.0图6长&比对流致振动的影响

18、/=-8.14X10-刈+0.2522(1)其中，代表剖面仪的长径比。振动加速度幅值。随长径比的变化按三次多项式进行似合，其关系式为Q=一0.0262工j+1.1775犬-17.555.+88.72(2)拟合关系式的R平方值为0.9896,说明拟合公式可信。剖面仪下降过程中，要求振动的幅值越小越好，而且振动频率应远离剖面仪的固有频率，因此，优化模型中选取振动加速为优化目标，以剖面仪的固有频率为约束条件。优化模型的目标函数为min(a)10与115.6Hz(5)求解式(3)和式(5)组成的优化模型，得到长径比的最优解为12。3.2重浮心距离计算时保持模型形状不变，因此浮心位置不变，改变重心位置，

19、使重浮心距离在150450mm之间变化。图7所示为削面仪流致振动加速度和振动频率随重浮心距离的变化曲线。由图7可看出，重浮心距离与剖面仪流致振动频率成正比关系，其关系式可表述为/=0.376x24-0.0365(6)其中，了2代表重浮心距离。而流致振动加速度幅值随重浮心距离的变化则较为笈杂。本文采用五次多项式进行拟合。五次多项式拟合流致振动加1.5L，0.08150250350450政浮心距腐ixi/mm图7重浮心距离对流致振动的赃响速度与剖面仪重浮心距离的关系式为a=-9708葛+11686工；-4843葛+753.3x1-17.95x2-0.1028(7)拟合关系式的R平方值为0.9359

20、.说明拟合公式可用。实际中重浮心距离变化的范围为200400mm。仍然以剖面仪流致振动的加速度幅值最小作为优化目标，以剖面仪的固有频率作为约束条件，建立优化目标模型。目标函数为min(a)200mm以400mm(8)约束条件为式(5)。联立式(8)和式(5),得到重:浮心距离的最优解为322mm。4优化效果分析优化后的剖面仪结构外形及相关尺寸如下:剖面仪头部形状为半球形，主体为长圆柱体，无凹阶结构.K径比为12,重浮心距离为322mm。优化前剖面仪结构外形及相关尺寸如下：剖面仪头部形状为平面，主体为长圆柱体，有凹阶结构(图2a),长径比为11.重浮心距离为360mm。分别对优化前后的剖面仪进行流致振动计算，得到优化前后剖面仪流致振动加速度幅值分别为3.85mm/s2和2.52mm/s2剖面仪的测觉精度主要指剖面仪对海洋湍流动能耗散率的测量范围，即测量下限和测量匕限。其中剖面仪测定上限主要由其搭载的剪切流传感器决定，而测段下限则主要由剖面仪的流致振动决定。为验证结构优化后剖面仪测量:精度改善情况，可对结构改进前后剖面仪的测量下限进行比较。剖面仪的测杭下限估算公式为j=7.5Z.)，LiuYuhong.WangYanhui,ctal.VibrationSourceAnalysisandVibrat