ANSYS数值计算在土木工程中的应用及进展上.ppt

土木工程结构数值分析 概述 1 现代土木工程结构正向大型化 复杂化方向发展 需要高精度的结构分析与精确的过程仿真 作为结构设计 施工的基础 利用数字技术替代结构整体试验了解结构受力性能已成为经济 可靠的捷径 计算力学的发展为结构数值分析提供了理论基础 有限差分法 有限元法 边界元法 新算法研究 并行计算 结构优化 反问题算法 本构理论的发展等 计算机科学的发展为高性能计算提供了条件 三维可视化建模技术 并行大容量计算系统 网络计算机 概述 2 1991年 美国 国家关键技术委员会 提出计算机仿真分析与建模是美国新时期应优先发展的关键技术之一 工程研究进入新时期 美国自然科学基金建立机构 GeorgeE BrownJr NetworkforEarthquakeEngineeringResearch NEES 专门从事地震工程研究 集实验 计算为一体 利用实验室和宽带网构建研究平台 15个装备先进的大型实验室 网络计算系统NEESgrid 数值分析的实现需要操作者具备良好的计算力学基础 工程结构知识 分析经验 熟练建模 选择计算方法 对结果作出正确的解释是任一项分析所要求的 结构分析对软件的实用性 准确性和适用性等方面的要求越来越高 工程项目的有限元计算 成都市市政中心抗震性能分析某雷达天线塔考虑土 桩 上部结构共同工作自振特性研究福堂水电站调压井开挖稳定性分析高技术厂房结构微振响应分析 成都市政中心4号搂 典型平面图 结构总高度72 4m 总宽52 5m 单肢宽13 5m 结构长度134 368m 结构沿x轴对称 沿y轴严重不对称 竖向布置上采取了一边退台 一边挑出的方案 使重心逐渐向左偏移 各层质量 刚度相差较大 层平面图 第一层 第十一层 第十五层 第四层 有限元建模 杆件单元 BEAM188 楼板 墙单元 SHELL181 整个模型共为69243个单元 55907个节点 SHELL181 混凝土本构关系 Saenz公式表示的混凝土应力 应变关系 计算中简化的混凝土本构关系 钢材用两线型本构模型 计算工况 模态分析谱分析中震下弹性时程分析大震下弹塑性时程分析 Mises屈服准则 随动强化模型 输入地震波 人工模拟波 美国SanFernandoM6 5地震记录二个分量 结构自振周期及振型 第一周期1 69s 第三周期1 12s 第二周期1 54s 位移反应 最大水平位移和同时刻平均位移比值 层间位移 变形形式 头部弯剪型 尾部弯曲型 剪力墙应力反应 柱内力反应 TH1波作用下柱最大剪力及剪压比 柱等效应力 楼板应力反应 大震下位移反应 4号楼抗震性能分析结论 该结构具有较好的抗震性能 中震作用下 结构基本处于弹性状态 大震作用下 有足够的抵御变形的能力 保证不倒塌 中震作用下 结构变形呈弯剪型 无位移薄弱层 剪力墙承担了85 的总剪力 说明剪力墙布置合理 结构左端剪力墙形成的筒体为薄弱区域 在中震作用下 局部达到屈服状态 应予以加强 柱剪压比能满足规范规定的小震下的限值 表明柱的抗剪承载力较高 轴压比均不大于0 7 应力分析表明 柱未见屈服 梁在圆弧端部和与连廊连接的部位出现明显的轴力 设计时应按偏压 拉 构件计算 雷达天线塔考虑土 桩 上部结构共同工作自振特性研究 某雷达天线塔 高9m 半径3 2m 采用人工挖孔扩底灌注桩基 共8根桩 桩径1m 承台厚1 3m 位于较软弱的地基上 因工艺要求需要了解结构的自振频率 有限元模型 采用三维实体元 上部结构单元数32872 节点113921 地基 承台和桩单元数17865 节点数46257 地基范围取2 5倍承台直径 四周为对称边界 底部固定 分析工况 上部结构 土 桩全系统自振频率分析 用弹簧替代土模型自振频率分析 弹簧刚度系数由地基土的静力分析获得 自振频率影响因素研究 工况1分析结果 第一振型 f1 1 6093Hz 第二振型 f2 1 6098Hz 第三振型 f3 1 9612Hz 如不考虑土的影响 按结构基顶固定考虑第一频率为32Hz 工况2分析结果 f1 3 35Hz f2 8 45Hz f3 16 95Hz f4 21 20Hz f5 32 74Hz 工况3分析结果 模拟不同加固地基情况时结构的自振频率 注 第一次换土指 0 3 4 3 m的土换为混凝土材料 第二次换土指桩底高度1m以上的原土换为C15混凝土 小结 当土较软弱时 土 桩 上部结构共同作用效应显著 计算结构的自振特性应考虑土 桩 上部结构的共同工作 当采用全模型计算结构的特征值时 应正确判断结构的模态和土局部模态 利用土弹簧模型计算特征值是可行的 但要考虑合理的方法确定弹簧刚度 福堂水电站调压井开挖稳定性分析 建立三维有限元模型 研究在自重荷载作用下围岩与衬砌的受力和变形情况 根据现场监测数据 反演初始地应力和材料参数 利用ANSYS软件的生死单元功能模拟调压井的施工开挖状况 以确定合理的施工顺序和施工方法 分析施工过程围岩的稳定性 研究调压井在支护 钢筋混凝土衬砌及围岩共同作用下的工作性能 工程概况 福堂水电站位于四川省阿坝藏族羌族自治州汶川县境内的岷江干流上 装机容量360MW 由首部枢纽 引水系统和厂区枢纽组成 引水隧洞全长约20km 最大引用流量251m3 s 调压井设计为开敞差动式 圆形断面 内径31m 深110m 为亚洲第一深井 有限元模型 围岩部分采用八节点六面体单元Solid45及其锐化的十节点四面体单元Solid92 衬砌部分采用混凝土单元Solid65 模型边界条件 四周侧立面采用对称边界条件 底部采用固定约束 材料分两层 上部为四类岩体 下部为三类岩体 单元数69316 节点数36803 反演分析 实测位移与计算位移对比表 根据勘测资料建议范围取不同的参数组合 利用现场监测数据 反复试算 反演 反演条件主要考虑不同高程后边坡的沉降为控制值 确定弹性模量为1260m高程以上E 2 5GPa 1260m高程下部为E 7GPa 自重荷载作用下结构变形情况 整体结构变形图 Z向位移分布图 最大竖向位移为10 8mm 发生在后边坡 井周最大竖向位移6 5mm 自重荷载作用下衬砌变形情况 井筒径向位移最大值6 19mm 发生在后边坡侧1260m高程 井筒环向位移最大值约5mm 发生在福堂沟侧 结构等效应力分布图 衬砌等效应力分布图 开挖稳定性分析 研究随开挖 衬砌的进行 调压井围岩的应力和位移的变化规律 通过对比分层开挖围岩应力与位移的变化 研究不同开挖施工方法对围岩稳定的影响 一次开挖 一次衬砌 一次衬砌 两层开挖 两次衬砌 未开挖 一次开挖 衬砌1 未开挖 开挖1 衬砌2 开挖2 应用生死单元 开挖分析应力位移结果 1260m高程截面控制点的应力位移值表 一次性开挖时1260m高程截面控制点的应力位移值表 结果分析 调压井有限元分析结果与监测成果基本相符 表明模型合理地模拟了复杂地形对结构体的影响 衬砌的分析表明其内力分布不均匀 局部应力较大 计算显示后边坡侧的变形较大 应该加强锚固 进行变形监测 以防边坡滑移和塌落影响调压井的正常使用 通过 二次开挖 衬砌 与 一次性开挖 衬砌 的计算结果对比分析 表明分层开挖时围岩应力 应变变化幅度相对较小 有利于围岩稳定 根据现场监测数据 进行的材料弹性模量的反演计算 所确定的参数能基本反映材料的性能 有一定的可靠度 为下一步多参数反演分析打下了基础 工作展望 结构数值分析 利用学术 技术和软件优势 为实际工程服务 解决各类结构计算问题与过程模拟仿真 结构类型 建筑结构 桥梁结构 岩土工程 道路工程 地下工程 水工结构等分析项目 非线性静 动力分析 弹塑性分析 地震反应分析 温度应力分析 上下部共同作用分析 工程安全稳定分析 施工模拟分析分析手段