基于响应面法ANSYS灵敏度分

基于应面法ANSYS灵敏度分析肖禧成〔东南大学机械工程学院,120176〕摘要：基于响应面法的ANSYS灵敏度分析是一种可靠有效的分析方法,本文详细阐述了ANSYS中基于响应面法的灵敏度分析基本原理,为基于ANSYS的结构灵敏度分析提供参考;应用ANSYS对一个流固耦合模型中固体接触面内应力受冷、热水流体温度的影响进行了灵敏度分析,并寻求一组最正确参数值，使当冷、热水入口的速度、温度在一定范围内变化的情况下，使中心块的内应力最小。并通过对计算结果的分析，验证了基于响应面法ANSYS灵敏度分析的高效性和可靠性。关键词：灵敏度分析;响应面法;流固耦合;ANSYSAbstract:Thesensitivityanalysis(SA)basedonresponsesurfacemethodinANSYSisareliableandeffectivemethod,thebasicSAtheorybasedonresponsesurfacemethodinANSYSispresentedindetail,whichcanbeareferenceofstructuralSAusingANSYS.AFluid-solidcouplingmodelisanalyzedbyusingANSYS,andthesensitivityvaluesofintra-stressonthecontractsurfacesofthesolidresponsestodifferentvelocityandtemperatureoftheinletandoutletarecalculated,whichcanbeusedtocomputeafittestparametermakingtheintra-stressminimum.TheanalysisofthecalculatedresultsverifiesthereliabilityandeffectivenessofSAbasedonexperimentaldata.Keywords：sensitivityanalysis(SA);responsesurfacemethods;Fluid-solidcoupling;ANSYS学习文档仅供参考学习文档仅供参考目录TOC\o"1-5"\h\z\o"CurrentDocument"响应面法 2\o"CurrentDocument"ANSYS灵敏度分析的基本原理 4\o"CurrentDocument"ANSYS概率设计系统〔PDS〕 4\o"CurrentDocument"基于Spearman秩相关系数⑼的ANSYS灵敏度分析 5\o"CurrentDocument"流场腔内固定块的应力场的灵敏度分析 5\o"CurrentDocument"ANSYSCFX分析 5\o"CurrentDocument"求解壁面应力灵敏度 11\o"CurrentDocument"结论 18\o"CurrentDocument"参考文献 191.响应面法响应面方法是进行灵敏度分析的一种有效方法，其思想是通过一系列确定性试验拟合一个响应面来模拟真实极限状态曲面.假设Z与系统随机参量Q=[Q,Q，…,Q]的关系可用式⑴描述，通过随机抽样法得到随机参量的N个样1 2本值，对这N个样本值数值计算得到系统响应的一组样本值(工,工,・・・,z)利用最1 2小二乘法得到该系统函数；用该响应面方程代替有限元模型进行失效模式分析，在结构响应Z未知的情况下，用响应面函数代替结构的真实响应，将大大节约计算时间【1】．aQ+iii=aQ+iii=1£R£RaQQijiji=1 j=i(1)式中a,a,a(i=1,…,R;j=i,・..，R)为待定系数，共n+1+n(n+1)/2个.0iij采用矩阵法对每个随机变量取三个水平点，按照某种法则得出中心所在点和边中点作为样本值点.图1.1表示三变量(Q,Q,Q)样本值点。123内◎中心点・边中点图1.1变量的样本点值.1SwatchofVariable

符合某分布的变量可用前法确定变量水平XsJaf(q)dq=p,n=1,2,3 (2)n一g式中f(q)为变量的概率密度函数；p为水平点，选取p=0.01,p=0.5,p=0.99.n 1 23(3)正态分布变量(3)q=N+gvt(p)

sn式中n为平均值，o为标准差，V(•)为标准正态分布函数，W-1(p)可查表.n对参变量的S个样本点值数值模拟，得出S个输出点(z,z,・..,z)，对数据1 2s进行回归分析得【2-4】s二[s二[zii=1i=1aq+工工aqq)]2

ii ijiji=1j=i(4)=0i=1,2,…小daiji=1,2,...R;j=i,...，Rdaij对该式进行求解计算，确定式(1)系数的估测值，进而确定系统的函数关系式。近年来，响应面法不仅仅在化学工业，生物学、医学、生物制药、食品学以及生态学领域都得到了广泛应用。同时，在工程学方面也涉及到了响应面法的应用【5】。在工程学方面，响应面法主要用来做结构的优化设计和可靠性分析，而对于计算数控机床精度可靠性尚属首次，对于大型复杂机械系统来说，如何有效地计算其可靠度具有重要的现实意义。灵敏度分析(SensitivityAnalysis)是评价因设计变量或参数的改变而引起结构响应特性变化率的方法。结构灵敏度的研究是一个很特别的领域,它是当前计算力学和结构工程领域的主要研究方向之一。结构灵敏度分析在确认系统的前后关系、优化算法、响应量测的仪器精度、系统性能的可靠性评估以及结构冗余度［5］研究等方面发挥了重要作用。实际上,在确立结构优化、可靠性评估和参数识别时,结构灵敏度分析是一个主要的先决条件【6】。文【6,7】对各种灵敏度分析方法进行了总结回忆,图1.2为不同类型的灵敏度分析方法【6】。其中,基于响应面法的灵敏度分析方法是一种可靠有效的分析方法,而在结构分析领域广泛应用的大型通用软件ANSYS即提供了此功能⑻，但国内鲜有基于ANSYS的结构灵敏度分析文献,本文对ANSYS灵敏度分析原理的详尽阐述以及算例可为基于^5丫5的结构灵敏度分析提供借鉴和参考。本文采用基于响应面法的分析方法,应用有限元软件ANSYS对一个流固耦合模型中固体接触面内应力受冷、热水流体温度的影响进行了灵敏度分析,并寻求一组最正确参数值，使当冷、热水入口的速度、温度在一定范围内变化的情况下，使中心块的内应力最小。并通过对计算结果的分析，验证了基于响应面法/ANSYS灵敏度分析的高效性和可靠性。图1.2灵敏度分析方法分类Fig.1.2TheclassificationofSAmethodsANSYS灵敏度分析的基本原理ANSYS概率设计系统〔PDS〕ANSYS灵敏度分析包含于其概率设计系统⑻〔PDS〕模块中。ANSYS概率设计系统是一种基于有限元的概率设计,用来评估输入参数的不确定性对于系统输出的影响行为及其特性。其概率设计方法包括蒙特卡罗(Monte-Carlo)法和响应面法,本文采用概率分析中比较常用的响应面法。响应面法的使用参数设计和系

数估计。参数设计是在参数空间中通过选择合理的采样点，构造出符要求的简单，比较接近真实的响应面，使其在有限的区域内能够有效的逼近真实响应。系数估计的任务是利用有限的采样点数据求解响应面函数中的相关系数。目前常采用最小二乘法或加权最小二乘法。基于Spearman秩相关系数【9】的ANSYS灵敏度分析TOC\o"1-5"\h\z假设考察结构响应丁（比方最大位移）对结构参数x,x,・・.,x（比方m个杆件1 2m的截面面积）的灵敏度。首先指定每个结构参数的概率分布函数x,x,・・・,x此后1 2m的随机模拟计算中每个结构参数的取值分布即基于其概率分布函数。假设进行九次随机模拟运算，得到n个结构响应数值y,y,・・.,y，其中，第%•次模拟运算的结1 2n构响应为y，结构参数为x,x，…,x,i=l,2,.・・,〃。i 1i2i mi由第j（j=1,2,...,㈤个结构参数n次随机模拟运算的样本值x,x,・・.x和nj1j2jn次随机模拟运算结构响应样本值y,y,・・.,y构成如下m个数据对：1 2nxj1xj1y1x)j2y2)jnn对于每个数据对，利用式⑴求得其Spearman秩相关系数r（j=1,2,・・.,m）。sjTOC\o"1-5"\h\zr反应了结构响应y和第j个结构参数之间的相关性，r为正，说明结构应ysj sj随结构参数x的增大而增大，减小而减小;反之相反。而卜1反应了结构响应y和j sj结构参数x之间的相关程度，即本文所指的灵敏度。|r|越大,说明结构响应y对j sj于结构参数x越敏感，即结构参数x的变化对于结构响应y的影响越显著；反之相jj反。流场腔内固定块的应力场的灵敏度分析ANSYSCFX分析如图3.1.1是模型图，冷热流体从两个不同端口进入一腔体内，腔体中央是一中心块，其底部固定。由于受到冷、热水流体温度的影响，中央块内会产生

内应力。现要分析进、出水的水流速度及温度对中心块外表应力的影响，以及应力场对个参数的灵敏度。对上图流体模型在ANSYSWorkbenchmesh模块中划分网格如图3.1.2.

图3.1.2流体有限元模型Figure.3.1.2FEAmodeloffluid在ANSYSWorkbench中，采用ANSYSICEMCFD划分网格生成流体的有限元模型后，导入到CFX-Pre中，网格定义了我们需要求解的流体域。然后按照以下步骤执行，定义相关参数。定义流体域属性通用选项界面中基本设定下的流体列表可选择参与计算的流体种类，流体采用不可压缩流体，流体类型为水。域模型下设定参考压强为1[atm]；其它选项保持默认值。求解方法采用标准he模型，总仿真时间为，仿真时间步长取为，

总仿真步数为100步。如图总仿真步数为100步。如图3.1.3所示。De后也MDctaiiltDcnuinii'iFlawAnalysis1OutheDoman:DefedtDomainIDetaibofDefaultDomaininFlowAnalysis：1日血SehngsAudModels Iri的闻ionHeatTrarsferTOC\o"1-5"\h\zOption TnermalEnergy *二|[nd.ViscoiEC^ipatonlurtufenceOption k-^psfcn t|WalFincfan Scalable -]TiibJentFluxCtereHeatTransfer EAdvancedTlibden'SCnntrd ECcfflbustai□phon rtone *"[haemalRsdiaban□ptian Hone ・□Efedramwre怔Model图3.1.3定义流体域属性Figure.3.1.3Definitionoffluidattributes定义边界条件想要完整定义一个问题，需要设定边界条件。进行流体动力学分析时边界条件和载荷包括入口、出口、液体和管壁的接触三部分，本文模型需要设置3种边界条件。首先创建一个入口（皿1©力边界条件，只允许流体流入的计算域。在边界细节中，定义4个表达式，冷水入口速度coldinletvel=1.75m/s,温度00^1口10廿0山口=22℃,热水入口速度hotinletvel=1.0m/s,温度hotinlettemp=85℃,如图3.1.4所示。设定完入口边界条件后，其次创建出口（outlet）边界条件，允许流体流出流入计算域，设定为Opening，即开口边界，在边界细节中，设置出口相对静压为0Pa，出口不施加任何载荷和约束。最后创建壁面（centralblock）边界条件，不允许流体穿越。并将coldinletvel、coldinlettemp、hotinletvel和hotinlettemp均设为输入参数。OutlineDomain:DefeutDomainExpressicnsExpressonsJ向Expressionsgcoldinlettemp骂5n5将吗|csldinletvel 入箝师&人1]hotinlettemp3我石岗吗|Htinle加I图流体边界条件设置Fig.Fluidboundaryconditionsset求解器控制设定求解器控制表中的参数，如可设定最大循环步数、时间尺度选择等，这些数值会对收敛精度有影响，参数设定如图所示。图3.1.5求解控制设置Fig.Solutioncontrolestablishment〔5〕数值求解求解结束后，得到centralblock壁面处得温度，如图.

Figure.Temperatureofcentralblock求解壁面应力灵敏度经过CFX流体分析得到接触面温度分布函数后，便要建立Steady-StateThermal(ANSYS)项，再导入CFX的温度值，如图所示。IMultiple-Sj^Em's-MechanicalMuhiph^rsics]FileEdrtViewUnitsToolsHelpSake▼ShowErran乎lifc才置).因图,.'/c-xsheei*FileEdrtViewUnitsToolsHelp十府厅卜附喻痘啕用♦I3牛食既虱闻&Q半谏闻旨?ShowVertwe初附refnam® IIEdgeColoring* XwXHhlTiicksnAmcrtatiDn^口向Sho*M匹h点HFEnvironmen的hh- Conditions*Outline]FitterNan心・91mtedt&囹Hodri(乜E4)田Outline]FitterNan心・91mtedt&囹Hodri(乜E4)田•怎5m啊底■/■点CoaiirflteSystEns二短心Uli10.0淞dWitaheThwimal[D4J/*号1IritdTarctfaliireAnal*p3Kfilings(3-ytilWgfedUsd[泌血十;由■ CmparsedTemperainreSohtkm[Dfi)建「ISahbonInfbrmabankstrudijral[E51AnalyseSdtrnjpfixed汕DflrlImportwlLead(EflUrni)日二.喇SolutkrtfE6)R：Stg^dy-期r£金Tlh啪mlImpartedIemp^raiureUnit:叱2013/B/2720；026里州四号帕热fiijcug59JIMS5T1O3T55103153JO3B51JO249J014.juooa4£lW2fti1hDJJ05%, PrintPnevim%R^tHirtPrevisw/SourceTimeSaurDC-TimeStepAnaly&kTime(i)Scal«Offsel固T]10.级110图导入CFX的温度值然后插入EquivalentStress并使其Maximum为输出参数，便可以建立结构分析项了；之后选定中心块底面为固定端约束，求解后，再建立ResponseSurface项，如图所示，然后设置输入参数的范围，设置冷水入口温度为10〜35℃〔283〜308K〕，速度为1.0〜2.5m/s,热水入口温度为75〜97℃〔348〜368〕，速度设为0.5〜1.5m/s,如图所示。

图建立ResponseSurface项OutlineofSchematicF2:DesignofExperimentsAB1Enabled2El"DesignofExperiments3日InputParameters4日FluidFlow；CFX)仁邛5CpP5-coldinlettemp底16CpP4-hotinletvel回7举P2■coldinletvel回@--PP3-hotinlettemp国g日OutputParameters10□FluidFlow；CFX):C0iipmP6■TempspreadTlcw12.日国StaticSfe-uctural(E3J13-pJP7-EquivalentStressMaxinnumid.[=1

|PropertiesofSchematicF2;DesignofExpe--PXPropertiesofSdiematicF2:DesignofExpe▼4XABTAB■¥1PropertyValue1PropertyValue工sGeneral3□Gener^3UnitsK3Unitsms.A-14TypeDesignVariable4TypeDesignVariable5ChssifiotionContinuous '5ClarificationCantnuQus '6SValues6aValues7Value358.157Value1SLowerBound3^88LowerBound0.59UpperBound3689UpperBoundl.,510UseManufacturableLA.1..一一■10UsehlanufactLiratoleu_L ■P「口PE「P「口PE「tiRqof5±ematicF2:DesignofExpeAB41Property-ValueNsGeneral3UnitsIC4TypeDesignVariable5ChssificatiQnContinuaus二16sValues7Value235.158LowerBound2839UpperBound30810UseManufacturableET~AB1PropertyValue2HGeneral3Unitsms^-14TypeDesigny/ariable5ClassificationContinuQusJ三6sValues7Value1.758LowerBound1gUpperBaund2.510UseManufacturable□PropertiesofSchematicF2:DesignofE:图设置参数范围设置输入参数后，得到响应面与个输入参变量间的二维曲线图，如图。

3.44012.94013.4451ResponseChartfarP7-EquivalentStressMaximumResponseChartidrP7-EquivalentStressMaximum(or)EnExBs««whmauB-Enb山.GaResponseChartforP73.44012.94013.4451ResponseChartfarP7-EquivalentStressMaximumResponseChartidrP7-EquivalentStressMaximum(or)EnExBs««whmauB-Enb山.GaResponseChartforP7-EquivalentStressMaxinnurn■DMEnEHgEKMUh^luu-qentlUJ—《CL3.78853.6SS52.58853.48853.38853.2885b.less3.08S52.0885

R&gpanseChartftsrP7-EquivalentStressMswirriumj3.69113.59113.49113.39112.29113.19113.09112.09112.89112.79112.69110.5 0.6 0.7 0.8 0.9 1L.L1.2 1.3 1.4 1.5ulrlEMEW蛆值0.5 0.6 0.7 0.8 0.9 1L.L1.2 1.3 1.4 1.5ulrlEMEW蛆值3」4。4=3而>'5仃山1id图3.2.4不同参数间的二维曲线图最后得到响应面对个输入参数的灵敏度柱状图和曲线图，如图。图中TempspreadTlow为centralblock接触面处得最大温差，为设置的CFX输出参数。LocalSensitivityd?p«uucn-Em765.4..3.21-JQo.o.o..2o.4o.P5-coldinlettempP4-hotinletvelP2-coldinletvelP3-hotinlettempP6-TempEpraadTlowP7-EquivalentStressMaximumOutputParameters灵敏度柱状图LocalSensitivityCurvesforP6-TempspreadTlowP5PRseponsePcoldinlethP4-hotinle2-coldinlebint■smptveltvelP3-hotinletttsmpQQ.l0.2 0.3 0.4 0.5 0.6 Q.7 Q.8 0.9 1灵敏度曲线图图3.2.5个输入参数灵敏度结论从灵敏度图中可以看出，响应面〔即中心块的6个接触面〕对冷水进水口的温度变化[coldinlettem〕p最灵敏,热水口水流速度其次，对冷水口进水速度[coldinletvei变化灵敏度最低。出现这种现象的原因，可能跟流场腔体的形状及中心块在腔体中的位置有关，同时也与中心块的材料有关。参考文献[1]赵联春.球轴承振动的研究⑪