采用Hermite矩形单元建立类桁架连续体.docVIP

 PAGE 8  MACROBUTTON MTEditEquationSection2 Equation Chapter 1 Section 1 SEQ MTEqn \r \h \* MERGEFORMAT  SEQ MTSec \r 1 \h \* MERGEFORMAT  SEQ MTChap \r 1 \h \* MERGEFORMAT 采用Hermite矩形单元建立类桁架连续体 周克民，李霞 (华侨大学土木工程学院，福建362021) E-mail:zhoukm@hqu.edu.cn Tel: 0595Fax:595摘要：采用先建立类桁架连续体，再转化为离散杆系结构的方法研究了单工况应力约束下的结构拓扑优化方法. 在用有限元方法建立类桁架连续体时，以结点位置的杆件密度和方向为设计变量可以容易地构造C0阶连续的杆件分布场. 杆件在结点位置的密度和方向根据结点位置的应力或应变计算. 一般有限元分析方法直接得到结点的位移，由位移场求导得到应变和应力场. 所以应变和应力的精度比位移低. 结点位置的应变需要由围绕该结点的单元在结点位置的应变值平均. 本文采用Hermite矩形单元，将结点位置的应变直接作为基本变量，无需对位移求导. 1 引言 以连续体为基结构的拓扑优化方法中，均匀化方法[1-4]，进化结构优化方法[5-8]、ICM方法[9,10]是通过优化有限单元的“有”和“无”来实现结构拓扑优化，一般存在数值计算不稳定问题[11,12]. 、水平集方法和泡泡法通过优化连续体的边界（包括产生新边界）实现拓扑优化. 除均匀化方法以外，都是采用各向同性材料模型. 自由材料优化方法[13]把材料的弹性张量的所有分量作为优化设计变量，具有最大的设计空间，且不存在数值计算不稳定问题. 但优化结果与真实材料缺乏明确的对应关系，将其转化为真实材料仍是一件十分困难的工作[14,15]. Michell理论揭示了拓扑优化结构的一个重要基本性质：类桁架连续体. 本文通过有限元方法建立这种类桁架连续体. 类桁架连续体难以制造，结果不便于在工程上直接应用. 但一些研究表明，离散化的Michell桁架与理论上的最优解非常接近[16,17]. 所以，我们首先建立类桁架连续体，实现理论上的最优化；再将其转化为杆系结构，实现工程上的可用性. 作为一种新的拓扑优化方法[18,19]，具有重要研究价值. 它改变了目前各种结构拓扑优化方法中结构拓扑大多都是用有限元网格来描述的方法. 2 类桁架材料模型 本文仅限于研究单工况应力约束平面应力状态的拓扑优化问题，即Micehll桁架问题，或称最小重量桁架. 研究表明，它也是体积约束下的最小柔度设计，还是拓扑优化的带孔板在体积比趋于零时极限状态的结构. Michell桁架一般是由无间隙的无限细的杆件构成的连续体，又称为类桁架连续体（truss-like）。在Michell桁架中共有5种杆件分布状态，如图1b所示.图1c描述了对应应变状态的应变圆，图1a和图1d分别是各种杆件分布状态的符号和名称. 本文采用正交异性材料模型描述上述杆件分布，即假设在任意点位置杆件正交分布. 这种材料模型可以容易地描述正交两向拉压区域. 对于单向拉伸和压缩的区域，只要令一个方向的杆件密度为零即可. 但对于“任意方向杆件是优化的”的2种区域则不能描述. 但在实际优化问题中，这种区域一般仅限于在奇异点出现. 对主要区图1 Michell桁架中杆件的各种分布状态 a) b) c) d) 单向拉伸 单向压缩 各向均匀拉伸 两向正交拉压 各向均匀压缩 域，这种材料模型仍然是有效的. 假设两个正交方向的杆件密度为，应变、应力和单位宽度上的内力分别记作和(i=1, 2). 假设材料遵守虎克定律，即， ,  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 和， .  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 假设剪切力和剪应变的本构关系 .  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 由于Michell桁架中平行杆件之间