均匀设计法与有限元网格划分结合可行性.pdfVIP

均匀设计法与有限元网格划分结合的可行性 庞毅，张旭明，张慧 河海大学工程力学系，南京 (210098) E-mail ：pangyi995@ 摘 要：有限元在划分网格的过程中,有很多都是凭个人经验来划分的,疏密的程度没有统一 的标准哪里应该加密哪里应该划疏，或者是知道疏密的程度但苦于计算机软件与程序的 实现或精度的准确。自适应划分网格是一种比较有效的方法,但本文从均匀设计的理论出 发借助MATLAB 强大的程序模拟,可以得出网格疏密的划分重点,相信在用ANSYS 划分网 格时可以提供强有力的帮助 。具体思路是在求解同一个工程问题时,分别用均匀设计的方 法由MATLAB 计算得出结果(应力最优点),再与有限元方法由 ANSYS 计算得出结果(应力 梯度) 比较,应力变化图如果基本吻合的话,就给均匀设计与有限元网格划分结合的可行性 提供了强大的佐证,亦即当ANSYS 在划分网格前用 MATLAB 优化的结果的程序结合提供 了可行性。因为结果说明均匀设计得到的优化点等价于受力大或变化梯度大的点，加密梯 度变化大的点当然精度更高。 关键词： 均匀设计；有限元；网格划分；结构优化设计；梯度变化 中图分类号：O 39 1．引言 均匀设计是考虑试验点在试验范围内均匀散布的一种试验方法，其优点是很大数量的因 素和水平可以用一个中等数量的实验来容纳[1]。1978 年王元和方开泰两位数学家提出了均匀 设计，采用均匀设计安排试验时，每个因素的每个水平只做一次试验，当水平数增加时试验 数随水平数的增加而增加,实验数的增加会带来更细微的结果[2],这与有限元划分网格的方法 是一致的。 均匀设计是通过一套精心设计的表来进行试验设计的，这就是均匀设计表[3]。每个表 有一个代号 UL (q n ) ，其中“U”表示均匀设计，“L”表示试验数，“q”表示每个因素的水平 数“n”表示列数。 王元和方开泰先生用均匀设计思想来做的配方均匀设计，生成的表称为配方均匀表， 记为 n UML (L ) ，其设计方案如下： 1给定L ，由它生成的向量为H= （ h ,, h ），取H 前n-1 项（h ,, h ）为生成向量 1 n 1 n−1 HH ，并由此生成向量产生均匀设计表UL (Ln−1 ) ，用{qki }记UL (Ln−1 ) 中元素。 2对每个i，计算 2qki −1 Cki = ， k=1,…,L 2L 3计算 - 1 - 1 1 ⎧ ⎛ ⎞i −1 ⎜ n−i ⎟ n−j − x n 1−C C i 1,2, L, n 1 ⎪ ki ⎜ ki ⎟∏ kj ⎪ ⎝ ⎠j 1 ⎨ 1 ⎪ n−1