自然单元法在水工结构数值分析中的应用.pdfVIP

第 30 卷第 1期 华 　北 　水 　利 　水 　电　学 　院 　学 　报 V ol130 No11 200 9年 2 月 Jou rna l o fN orth Ch in a In stitu te of W ate r Con servancy and Hydroe lec tric Powe r Fe 1 2009 文章编号 : 1002 - 5634 ( 200 9) 0 1 - 00 31 - 0 3 自 然 单 元 法 在 水 工 结 构 数 值 分 析 中 的 应 用 于国辉 , 赵继伟 (华北水利水 电学 院 ,河南 郑州 450011) 摘 　要 : 自然单元法既具有无网格方法和经典有限元方法的优点又克服了两者的一些缺点 ,为水工 结构分析 开拓了新途径 . 编制了弹性 自然单元法 C + +程序 SGN EM 1. 0 ,利 用其对几个水工 结构物进行分析 , 通过和有 限元分析结果进行比较 分析 ,结果表明自然单元法在水工结构数值分析上 的应 用是可行的. 关 键 词 : 自然单元法 ;重力坝 ;无网格方法 ;有限元方法 ; C + +语言 中图分类号 : TV 3　　　文献标识码 : A 　　自然单元法基于弗洛诺伊 (Vorono i) 图和德劳 为张 、压情况下模量没有改变的理想情况 [ 4 - 6 ] . 材料 ( ) ( 内 D e lau nay 三 角形 几何 结构 , 采用 希布 逊 Si 2 属性见表 1. ) ( ) son 或非希布逊 N on—Si son 插 值方法全域构造 插值函数 , 既具有无 网格方法和经典有限元方法 的 优点又克服了两者 的一些缺点 ,它可 以像有限元法 一样直接施加本质边界条件 ,克服 了基于移动最小 二乘拟合 的无 网格方法在这一点的不足 ,而且 自然 单元法是无网格方法 ,可 以方便地处理有限元方法 较难处理 的一些问题 ,例如移动边界和大变形等 问 题 [ 1 - 3 ] . 笔者将 自然单元法 的理论和水工结构 的特 点结合 ,利用 C ++语言编制程序进行计算 ,验证该 方法在水工结构分析中的可行性. 1　编 　程 自然单元法的程序实现需要为积分准备背景积分 网格的 Delaunay三角化、自然邻接形函数 的计算以及 获得和求解平衡方程等 ,这里主要列出基于 C ++语言 的 SGN EM 1. 0编程的主要步骤 ,如图 1所示. 2　重力坝弹性分析 2. 1　模型简介 如图 2 是重力坝连 同地基 的结构剖面 图. 为简 单起见 ,认为坝体和坝基分别 由一种材料构成 ,同时 也不考虑材料在张 、压情况下模量不同的效应 ,而认 收稿 日期 : 2 008 - 09 - 2 7 作者简介 :于国辉 ( 197 4 —) ,男 ,黑龙江兰西人 ,助教 ,硕士 ,主要从事水工结构与数值 仿真方面 的研究. 32 　华 　北 　水 　利 　水 　电　学 　院 　学 　报 　　　　　　　　　　　2 009 年 2 月 表 1　重力坝材料参数表 2. 2　模型的离散和背景积分网格 材料 模量 M / GPa 泊松比 μ 重度 γ/ (N ·m - 3 ) 图 3 是坝体 的离散点的分布图 ,为 了清晰地加 坝体材料 40 0. 167 26 52 0 以显示 , 以线交点表示离散点. 图 4是整体模型的背 坝基材料 30 0. 2 00 25 000 景积分网格 图 ,图 5为指定取值点的示意 图. 图 3　坝体计算模 型的节点离散 图 4　坝体坝基的背景积分网格 较. 在图 6 中 ,纵坐标为 目标值 ,横坐标为相对位置 ,