基于有限元的振动筛拓扑优化分析.pdfVIP

学兔兔 石油天然气学报 (江汉石油学院学报)2014年8月第36卷 第8期 JournalofOilandGasTechnology (J．JPI)Aug．2014Vo1．36No．8 基于有限元的振动筛拓扑优化分析 闵庆新 (华北油田天成实业有限公司飞达公司，河北~)／xll061ooo) 邓杰文 (华北电力大学能源动力与机械工程学院，北京1ooooo) [摘要]大型直线型振动筛广泛应用于石油、矿产等工业生产过程 中，振动筛的物料处理量直接影响生产 效率 ，大型振动筛虽处理量大，但存在设备 自重大、振幅低等 问题 。采用有限元方法对振动筛结构进行 拓扑优化设计 ，并 以对优化结果进行强度分析和模态分析 ，分析结果表 明优化后 的结构 自重降低 ，提高 了固有频率 ，其强度符合实际生产 需求，该优化过程为大型振动筛的设计工作提供可靠 的理论依据。 [关键词]振动筛 ；有限元法；拓扑优化；结构优化 [中图分类号]TE95l [文献标志码]A [文章编号]1000—9752(2014)08—0165—04 为了适应石油工业生产发展需求 ，必须发展大型生产机械 。振动筛是石油固控 系统的主要设备之 一 ， 采用大型振动筛可以提高生产效率 ，降低生产成本 川。然而，大型振动筛存在 自身重量较大 、振幅 较低等问题 ，直接影响实际生产效率 。笔者采用有限元方法对结构进行拓扑分析 ，并对优化后的结构进 行静态分析，对振动筛的优化设计结果进行验证 ，为大型振动筛的设计和优化提供可靠的理论依据 。 1 拓扑数学模型 在进行连续结构的实际优化计算时，一般采用变密度法作为其理论依据 。所谓变密度法[2就是假设 一 种密度存在变化且不是均匀的理想材料 ，采用有限元网格划分的方法将结构体离散为有限数量的网格 单元 ，然后对各个网格单元的密度进行数值附加 ，将附加到网格单元上的密度值作为计算变量 ，则有 ： p— x po (1) 式 中：x 代表各个 网格单元 的相对密度 ；p。代表在计算域里的每个 网格单元的附加密度值 ；p代表拓扑计 算变量。当各个 网格单元的相对密度 x 一 1时 ，说明该单元 内存在材料，因此将对该单元不做删除只 做加减处理；当X 一 。时，说 明该单元不存在材料 ，则对该单元进行删除处理。通过上述的判断方法 对 网格单元进行加减和删除操作就完成了对结构的拓扑过程 ，此外在进行拓扑优化时，应该尽可能使材 料 的相对密度为 。或 1并且分布在计算域 内。 由于材料的相对密度只有 。或 1两种数值 ，这就造成拓扑计算变量p不存在连续性 ，因此对应的导 数方程在连续变量优化计算时是无法求解的。为使求解方程得以继续 ，须对方程进行连续化_3]。 K 一 (X )K。 (2) 式 中：K 为单元刚度，N／m；K。为单元 固有刚度 ，N／m； 为惩罚 因子 。 ‘ C ===F D — D KD (3) 式 中：C为柔顺度 ；F为载荷矢量 ；K 为位移矢量；D为变形位移矢量 。 对连续结构进行拓扑优化时通常的计算 目标为柔顺度最小 ，其中，柔顺度数值为单元总应变能数值 大小 ，根据上述 的计算 目标 ，再将体积 比作为约束条件 ，则有优化计算 的数学模型如式 (4)：求 X 一 (X】，X2，… ，X ) 。 min： C — F D [收稿 日期