基于广义乘子法的有约束的自由曲面形态创建.docxVIP

基于广义乘子法的有约束的自由曲面形态创建 自由曲线的结构以其优秀的建筑表现力成为当前空间结构发展的主要趋势之一。所谓的“自由曲线”是指不能通过分析和表达简单曲线的形状。 2006年崔昌禹 上述研究推动了自由曲面形态创建的发展, 但大多数创建没有引入约束条件.而实际结构设计时, 往往会受到一些因素的约束和限制, 因此在自由曲面形状优化算法中有必要引入约束条件以使曲面创建结果更为合理. 本文采用NURBS (非均匀有理B样条) 方法生成自由曲面, 以应变能最小为优化目标, 以体积为约束条件进行自由曲面形态创建, 并将广义乘子法与梯度法相结合对有约束的优化问题进行求解, 通过Fortran语言编制了有约束的优化算法程序.通过对一个正高斯自由曲面和一个负高斯自由曲面进行有约束的形态创建, 并与无约束的创建结果进行对比分析, 验证本文提出的方法的有效性和实用性. 1 优化问题的描述 形状优化是寻找结构的最优几何形状来满足所需结构性能的方法.本文的自由曲面形状优化问题表述如下 其中:f是优化目标, 本文采用能够反映结构刚度的应变能作为优化目标;g 2 数值优化方法 2.1 自由曲线几何建模方法 本文采用NURBS生成自由曲面.表达式如下 其中:P 2.2 基于广义验证方法的优化算法 广义乘子法 则公式 (1) 转换为: 之后采用基于有限差分法 其中:γ为优化步长, 通过黄金分割法确定 拉格朗日乘子λ的迭代方程为: 3 基于有限单元法的模型设计 将上述NURBS方法、广义乘子法和梯度法相结合, 采用Fortran语言编制相应程序, 从而实现自由曲面有约束的形状优化.其中, 结构分析采用有限单元法完成, 单元采用三角形薄壳单元 3.1 约束和有约束的最优尺寸 本算例模型为一方形正高斯自由曲面壳体, 采用混凝土材料, 两边均为24 m, 四边简支约束, 壳体厚度为0.1 m, 材料弹性模量为3 GPa, 泊松比为0.167, 壳体承受竖向均布荷载作用, 大小为8 kN/m 优化的目标是整体应变能最小, 约束为壳体体积小于等于76 m 图2为有约束条件优化过程中的应变能变化曲线和体积变化曲线 (C 进一步比较有约束和无约束的最优形状, 如图4所示, 观察分析可见, 加入约束条件后最优曲面的高度从14.24 m下降到了11.22 m. 3.2 算例2:表面活性剂表面改性 本算例模型为负高斯自由曲面壳体, 两方向的跨度分别为24 m和16 m, 两长边简支约束, 壳体厚度为0.1 m, 材料特性和上一算例相同, 壳体承受8 kN/m 优化的目标是整体应变能最小, 约束条件为体积小于等于34 m 将两种情况下的曲面优化结果进行比较 (如图7所示) .可见, 加入约束条件后, 最优曲面有变平的趋势, 最大高度从7.68 m下降到了4.43 m. 4 优化问题的描述及验证 本文考虑实际工程情况, 在曲面形态创建中引入约束条件, 采用NURBS (非均匀有理B样条) 方法生成自由曲面, 以应变能最小为优化目标, 以体积为约束条件进行自由曲面形态创建, 并将广义乘子法与梯度法相结合对该有约束的优化问题进行求解.最后, 通过对两个典型算例进行分析, 优化结果具有良好的收敛性, 验证了本文方法的适用性和有效性, 实现了实际工程在满足结构性能的同时对材料用量的要求与限制, 更具有实际应用价值