基于响应面法ANSYS灵敏度分

基于应面法基于应面法ANSYSANSYS灵敏度分析灵敏度分析 肖禧成 东南大学机械工程学院 120176 摘摘 要要 基于响应面法的ANSYS灵敏度分析是一种可靠有效的分析方法 本文详细阐述了 ANSYS中基于响应面法的灵敏度分析基本原理 为基于ANSYS的结构灵敏度分析提供参考 应 用ANSYS对一个流固耦合模型中固体接触面内应力受冷 热水流体温度的影响进行了灵敏度 分析 并寻求一组最佳参数值 使当冷 热水入口的速度 温度在一定范围内变化的情况下 使中心块的内应力最小 并通过对计算结果的分析 验证了基于响应面法ANSYS灵敏度分析 的高效性和可靠性 关键词关键词 灵敏度分析 响应面法 流固耦合 ANSYS Abstract Abstract The sensitivity analysis SA based on response surface method in ANSYS is a reliable and effective method the basic SA theory based on response surface method in ANSYS is presented in detail which can be a reference of structural SA using ANSYS A Fluid solid coupling model is analyzed by using ANSYS and the sensitivity values of intra stress on the contract surfaces of the solid responses to different velocity and temperature of the inlet and outlet are calculated which can be used to compute a fittest parameter making the intra stress minimum The analysis of the calculated results verifies the reliability and effectiveness of SA based on experimental data KeyKey words words sensitivity analysis SA response surface methods Fluid solid coupling ANSYS 1 目录目录 1 响应面法 2 2 ANSYS 灵敏度分析的基本原理 4 2 1 ANSYS 概率设计系统 PDS 4 2 2 基于 Spearman 秩相关系数 9 的 ANSYS 灵敏度分析 5 3 流场腔内固定块的应力场的灵敏度分析 5 3 1 ANSYS CFX 分析 5 3 2 求解壁面应力灵敏度 11 4 结论 17 参考文献 18 2 1 响应面法响应面法 响应面方法是进行灵敏度分析的一种有效方法 其思想是通过一系列确定 性试验拟合一个响应面来模拟真实极限状态曲面 假设与系统随机参量Z 的关系可用式描述 通过随机抽样法得到随机参量的N个样 12 r QQ QQ 1 本值 对这个样本值数值计算得到系统响应的一组样本值利用N 12 s z zz 最小二乘法得到该系统函数 用该响应面方程代替有限元模型进行失效模式分 析 在结构响应未知的情况下 用响应面函数代替结构的真实响应 将大大Z 节约计算时间 1 0 11 RRR iiijij iij i ZaaQa QQ 1 式中为待定系数 共个 0 1 iij a a a iR jiR 1 1 2nn n 采用矩阵法对每个随机变量取三个水平点 按照某种法则得出中心所在点和 边中点作为样本值点 图1 1表示三变量 样本值点 123 Q Q Q 图1 1 变量的样本点值 Fig 1 1 Swatch of Variable 3 符合某分布的变量可用前法确定变量水平 s x s q n f q dqp 1 2 3n 2 式中为变量的概率密度函数 为水平点 选取 f q n p 正态分布变量 123 0 01 0 5 0 99 ppp 1 sn qp 3 式中为平均值 为标准差 为标准正态分布函数 可查表 1 n p 对参变量的个样本点值数值模拟 得出个输出点 对数据SS 12 s z zz 进行回归分析得 2 4 2 0 111 sRRR iiiijij iiij i szaa qa q q 4 令 0 0 0 1 2 0 1 2 i ij s a s iR a s iR jiR a 对该式进行求解计算 确定式系数的估测值 进而确定系统的函数关系式 1 近年来 响应面法不仅仅在化学工业 生物学 医学 生物制药 食品学 以及生态学领域都得到了广泛应用 同时 在工程学方面也涉及到了响应面法 的应用 5 在工程学方面 响应面法主要用来做结构的优化设计和可靠性分 析 而对于计算数控机床精度可靠性尚属首次 对于大型复杂机械系统来说 如何有效地计算其可靠度具有重要的现实意义 灵敏度分析 Sensitivity Analysis 是评价因设计变量或参数的改变而引 起结构响应特性变化率的方法 结构灵敏度的研究是一个很特别的领域 它是当 前计算力学和结构工程领域的主要研究方向之一 结构灵敏度分析在确认系统 的前后关系 优化算法 响应量测的仪器精度 系统性能的可靠性评估以及结 4 构冗余 度 5 研究等方面发挥了重要作用 实际上 在确立结构优化 可靠性评估和参 数识别时 结构灵敏度分析是一个主要的先决条件 6 文 6 7 对各种灵敏 度分析方法进行了总结回顾 图1 2为不同类型的灵敏度分析方法 6 其中 基 于响应面法的灵敏度分析方法是一种可靠有效的分析方法 而在结构分析领域广 泛应用的大型通用软件ANSYS即提供了此功能 8 但国内鲜有基于ANSYS的结 构灵敏度分析文献 本文对ANSYS灵敏度分析原理的详尽阐述以及算例可为基于 ANSYS的结构灵敏度分析提供借鉴和参考 本文采用基于响应面法的分析方法 应用有限元软件ANSYS对一个流固耦合模型中固体接触面内应力受冷 热水流体 温度的影响进行了灵敏度分析 并寻求一组最佳参数值 使当冷 热水入口的速 度 温度在一定范围内变化的情况下 使中心块的内应力最小 并通过对计算 结果的分析 验证了基于响应面法ANSYS灵敏度分析的高效性和可靠性 图1 2 灵敏度分析方法分类 Fig 1 2 The classification of SA methods 2 2 ANSYSANSYS 灵敏度分析的基本原理灵敏度分析的基本原理 2 12 1 ANSYSANSYS 概率设计系统 概率设计系统 PDSPDS ANSYS灵敏度分析包含于其概率设计系统 8 PDS 模块中 ANSYS概率设 计系统是一种基于有限元的概率设计 用来评估输入参数的不确定性对于系统输 灵敏度分析方法灵敏度分析方法 基于实验数据基于实验数据基于计算策略基于计算策略 回归法回归法概率法概率法变分法变分法离散法离散法 直接法直接法伴随法伴随法基于有限元法基于有限元法基于理论公式基于理论公式 半解析法半解析法有限差分法有限差分法解析法解析法直接法直接法伴随法伴随法 5 出的影响行为及其特性 其概率设计方法包括蒙特卡罗 Monte Carlo 法和响 应面法 本文采用概率分析中比较常用的响应面法 响应面法的使用参数设计和 系数估计 参数设计是在参数空间中通过选择合理的采样点 构造出符要求的 简单 比较接近真实的响应面 使其在有限的区域内能够有效的逼近真实响应 系数估计的任务是利用有限的采样点数据求解响应面函数中的相关系数 目前 常采用最小二乘法或加权最小二乘法 2 22 2 基于基于 SpearmanSpearman 秩相关系数秩相关系数 9 9 的 的 ANSYSANSYS 灵敏度分析灵敏度分析 假设考察结构响应 比如最大位移 对结构参数 比如m 个杆件y 12 m x xx 的截面面积 的灵敏度 首先指定每个结构参数的概率分布函数此后 12 m x xx 的随机模拟计算中每个结构参数的取值分布即基于其概率分布函数 假设进行 次随机模拟运算 得到个结构响应数值 其中 第 次模拟运算的nn 12 n y yy i 结构响应为 结构参数为 i y 12 1 2 iimi xxxin 由第个结构参数n次随机模拟运算的样本值和j 1 2 jm 12 jjjn xxx 次随机模拟运算结构响应样本值构成如下 个数据对 n 12 n y yy m 12 12 1 2 jnjj n xxx jm yyy 对于每个数据对 利用式 1 求得其Spearman秩相关系数 sj r 1 2 jm 反应了结构响应y和第个结构参数之间的相关性 为正 说明结构应y sj rj sj r 随结构参数的增大而增大 减小而减小 反之相反 而 反应了结构响应y和 j x sj r 结构参数之间的相关程度 即本文所指的灵敏度 越大 说明结构响应y对 j x sj r 于结构参数越敏感 即结构参数的变化对于结构响应y的影响越显著 反之相 j x j x 反 6 3 流场腔内固定块的应力场的灵敏度分析流场腔内固定块的应力场的灵敏度分析 3 13 1 ANSYSANSYS CFXCFX 分析分析 如图3 1 1是模型图 冷热流体从两个不同端口进入一腔体内 腔体中央是 一中心块 其底部固定 由于受到冷 热水流体温度的影响 中央块内会产 生内应力 现要分析进 出水的水流速度及温度对中心块表面应力的影响 以及应力场对个参数的灵敏度 图3 1 1 模型图 Figure 3 1 1 Model graphic 对上图流体模型在ANSYS Workbench mesh模块中划分网格如图3 1 2 冷水进口 温度为 15 35 速度为 1 2 5m s 热水进口 温 度为 79 95 速度为 0 5 1 5m s 水出口端 相对压力为 0Pa 中间固定块 7 图3 1 2 流体有限元模型 Figure 3 1 2 FEA model of fluid 在ANSYS Workbench中 采用ANSYS ICEM CFD划分网格生成流体的有限元模 型后 导入到CFX Pre中 网格定义了我们需要求解的流体域 然后按照以下步 骤执行 定义相关参数 1 定义流体域属性 通用选项界面中基本设定下的流体列表可选择参与计算的流体种类 流体 采用不可压缩流体 流体类型为水 域模型下设定参考压强为1 atm 其它选 项保持默认值 求解方法采用标准k 模型 总仿真时间为0 02s 仿真时间步 长取为0 0002s 总仿真步数为100 步 如图3 1 3所示 8 图3 1 3 定义流体域属性 Figure 3 1 3 Definition of fluid attributes 2 定义边界条件 想要完整定义一个问题 需要设定边界条件 进行流体动力学分析时边 界条件和载荷包括入口 出口 液体和管壁的接触三部分 本文模型需要设 置3 种边界条件 首先创建一个入口 inlet 边界条件 只允许流体流入的计 算域 在边界细节中 定义4个表达式 冷水入口速度coldinletvel 1 75m s 温 度coldinlettemp 22 热水入口速度hotinletvel 1 0m s 温度hotinlet temp 85 如图3 1 4 所示 设定完入口边界条件后 其次创建出口 outlet 边界条件 允许流体流出流入计算域 设定为Opening 即开口边界 在边界 细节中 设置出口相对静压为0Pa 出口不施加任何载荷和约束 最后创建壁 面 central block 边界条件 不允许流体穿越 并将 coldinletvel coldinlettemp hotinletvel和hotinlettemp均设为输入参 数 9 图3 1 4 流体边界条件设置 Fig 3 1 4 Fluid boundary conditions set 3 求解器控制 设定求解器控制表中的参数 如可设定最大循环步数 时间尺度选择等 这些数值会对收敛精度有影响 参数设定如图3 1 5 所示 冷水流入口 coldinlet 热水流入口 hotinlet 流出口 outlet 壁面 central block 10 图3 1 5 求解控制设置 Fig 3 1 5 Solution control establishment 5 数值求解 求解结束后 得到central block壁面处得温度 如图3 1 6 图3 1 6 固体壁面温度 11 Figure 3 1 6 Temperature of central block 3 23 2 求解壁面应力灵敏度求解壁面应力灵敏度 经过 CFX 流体分析得到接触面温度分布函数后 便要建立 Steady State Thermal ANSYS 项 再导入 CFX 的温度值 如图 3 2 1 所示 图3 2 1 导入CFX的温度值 然后插入Equivalent Stress并使其Maximum为输出参数 便可以建立结构分 析项了 之后选定中心块底面为固定端约束 求解后 再建立Response Surface项 如图3 2 2所示 然后设置输入参数的范围 设置冷水入口温度为 10 35 283 308K 速度为1 0 2 5m s 热水入口温度为75 97 348 368 速度设为0 5 1 5m s 如图3 2 3所示 12 图3 2 2 建立Response Surface项 13 图3 2 3 设置参数范围 设置输入参数后 得到响应面与个输入参变量间的二维曲线图 如图3 2 4 14 15 图3 2 4 不同参数间的二维曲线图 最后得到响应面对个输入参数的灵敏度柱状图和曲线图 如图3 2 5 图中 TempspreadTlow为central block接触面处得最大温差 为设置的CFX输出参数 16 灵敏度柱状图 17 灵敏度曲线图 图3 2 5 个输入参数灵敏度 4 结论结论 从灵敏度图中可以看出 响应面 即中心块的6个接触面 对冷水进水口的 温度变化 coldinlettemp 最灵敏 热水口水流速度其次 对冷水口进水速度 coldinletvel 变化灵敏度最低 出现这种现象的原因 可能跟流场腔体的形 状及中心块在腔体中的位置有关 同时也与中心块的材料有关 18 参考文献参考文献 1 赵联春 球轴承振动的研究 D 杭州 浙江大学 2003 Zhao Lianchun Research on vibration of ball bearing D Hangzhou Zhejiang University 2003 in Chinese 2 Harsha S P Kankar P K Stability analysis of a rotor bearing system due to surface waviness and number of balls J International Journal of Mechanical Sciences 2004 46 7 1057 1081 3