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

土木工程结构数值分析概述（1）现代土木工程结构正向大型化、复杂化方向发展，需要高精度的结构分析与精确的过程仿真，作为结构设计、施工的基础；利用数字技术替代结构整体试验了解结构受力性能已成为经济、可靠的捷径；计算力学的发展为结构数值分析提供了理论基础，有限差分法、有限元法、边界元法…，新算法研究（并行计算、结构优化、反问题算法），本构理论的发展等；计算机科学的发展为高性能计算提供了条件，三维可视化建模技术、并行大容量计算系统、网络计算机。概述（2）

1991年，美国“国家关键技术委员会”提出计算机仿真分析与建模是美国新时期应优先发展的关键技术之一；工程研究进入新时期，美国自然科学基金建立机构“GeorgeE.BrownJr.

NetworkforEarthquakeEngineeringResearch(NEES)”，专门从事地震工程研究，集实验、计算为一体，利用实验室和宽带网构建研究平台。15个装备先进的大型实验室+网络计算系统NEESgrid；数值分析的实现需要操作者具备良好的计算力学基础+工程结构知识+分析经验，熟练建模、选择计算方法、对结果作出正确的解释是任一项分析所要求的。结构分析对软件的实用性、准确性和适用性等方面的要求越来越高。成都市政中心4号搂典型平面图结构总高度72.4m，总宽52.5m，单肢宽13.5m。结构长度134.368m；

结构沿x轴对称，沿y轴严重不对称，竖向布置上采取了一边退台、一边挑出的方案，使重心逐渐向左偏移，各层质量、刚度相差较大。层平面图第一层第十一层第十五层第四层有限元建模杆件单元：BEAM188楼板、墙单元：SHELL181整个模型共为69243个单元、55907个节点混凝土本构关系Saenz公式表示的混凝土应力-应变关系计算中简化的混凝土本构关系钢材用两线型本构模型计算工况模态分析谱分析中震下弹性时程分析大震下弹塑性时程分析Mises屈服准则+随动强化模型输入地震波人工模拟波美国SanFernandoM6.5地震记录二个分量位移反应最大水平位移和同时刻平均位移比值楼层A1波TH1波TH3波Umax(mm)U0Umax/U0Umax(mm)U0Umax/U0Umax(mm)U0Umax/U013.081.621.913.842.001.923.701.951.90320.110.91.8424.413.61.8023.513.21.78431.717.41.8238.023.61.6137.121.41.73659.833.11.8171.045.11.5770.041.01.71775.841.21.8989.455.41.6188.650.61.751012965.71.9715584.61.8415178.21.931216789.31.872031171.741931021.8814205120.01.702491521.632341391.6815223140.01.592711761.542541611.57层间位移变形形式：头部弯剪型，尾部弯曲型柱内力反应楼层1区2区3区V（kN）V/fcAV（kN）V/fcAV（kN）V/fcA15260.0174780.0273560.02622810.0091350.0082800.02034500.0144460.02610930.07944100.0134460.0267980.05854470.0146050.0353720.035TH1波作用下柱最大剪力及剪压比

柱等效应力楼板应力反应大震下位移反应雷达天线塔考虑土-桩-上部结构共同工作自振特性研究

某雷达天线塔，高9m，半径3.2m,采用人工挖孔扩底灌注桩基，共8根桩，桩径1m，承台厚1.3m，位于较软弱的地基上。因工艺要求需要了解结构的自振频率有限元模型

采用三维实体元，上部结构单元数32872，节点113921，地基、承台和桩单元数17865，节点数46257。地基范围取2.5倍承台直径，四周为对称边界，底部固定。分析工况上部结构-土-桩全系统自振频率分析；用弹簧替代土模型自振频率分析，弹簧刚度系数由地基土的静力分析获得；自振频率影响因素研究。工况2分析结果f1=3.35Hzf2=8.45Hzf3=16.95Hzf4=21.20Hzf5=32.74Hz工况3分析结果12345678910未加固3.358.4516.921.232.734.266.767.5110.7110.8换第一层土3.613.708.7935.950.450.585.788.5118.4125.3桩周加混凝土4.1110.722.341.858.959.793.496.5120.1125.9第二次换土4.454.5610.843.559.159.193.997.4118.5125.5模拟不同加固地基情况时结构的自振频率注：第一次换土指（-0.3—4.3）m的土换为混凝土材料，第二次换土指桩底高度1m以上的原土换为C15混凝土。工程概况福堂水电站位于四川省阿坝藏族羌族自治州汶川县境内的岷江干流上，装机容量360MW，由首部枢纽、引水系统和厂区枢纽组成。引水隧洞全长约20km，最大引用流量251m3/s；调压井设计为开敞差动式，圆形断面，内径31m，深110m,为亚洲第一深井。

有限元模型围岩部分采用八节点六面体单元Solid45及其锐化的十节点四面体单元Solid92，衬砌部分采用混凝土单元Solid65。模型边界条件：四周侧立面采用对称边界条件，底部采用固定约束。材料分两层：上部为四类岩体，下部为三类岩体。单元数69316、节点数36803

。反演分析高程实测应力（MPa）实测位移（mm）计算应力（MPa）计算位移（mm）12303.263.9412608.165.429.696.19127010.627.6012803.484.91实测位移与计算位移对比表根据勘测资料建议范围取不同的参数组合，利用现场监测数据，反复试算、反演；反演条件主要考虑不同高程后边坡的沉降为控制值；确定弹性模量为1260m高程以上E=2.5GPa，1260m高程下部为E=7GPa。自重荷载作用下结构变形情况整体结构变形图Z向位移分布图最大竖向位移为10.8mm，发生在后边坡，井周最大竖向位移6.5mm自重荷载作用下衬砌变形情况井筒径向位移最大值6.19mm，发生在后边坡侧1260m高程；井筒环向位移最大值约5mm，发生在福堂沟侧结构等效应力分布图衬砌等效应力分布图

开挖稳定性分析研究随开挖、衬砌的进行，调压井围岩的应力和位移的变化规律。通过对比分层开挖围岩应力与位移的变化，研究不同开挖施工方法对围岩稳定的影响。一次开挖、一次衬砌一次衬砌两层开挖、两次衬砌未开挖一次开挖衬砌1未开挖开挖1衬砌2开挖2应用生死单元开挖分析应力位移结果载荷步控制点径向位移(mm)环向位移(mm)径向应力(Mpa)环向应力(Mpa)未开挖1-3.970.72-0.52-0.4320.70-3.46-0.35-1.6633.28-0.24-0.41-0.2241.343.57-0.29-1.55开挖11-7.260.31-0.04-0.362-0.45-3.37-0.05-0.8630.71-0.39-0.04-0.3540.143.71-0.06-0.94衬砌11-6.920.31-0.06-0.382-0.26-3.31-0.08-0.8830.87-0.39-0.05-0.3740.353.65-0.09-0.96开挖21-7.740.38-0.07-0.392-0.56-3.32-0.09-0.9030.65-0.42-0.09-0.3940.313.51-0.07-0.99衬砌21-7.780.39-0.07-0.412-0.66-3.49-0.10-0.9130.56-0.44-0.09-0.3940.253.48-0.09-0.991260m高程截面控制点的应力位移值表

一次性开挖时1260m高程截面控制点的应力位移值表

载荷步控制点径向位移(mm)环向位移(mm)径向应力(Mpa)环向应力(Mpa)未开挖1-4.000.57-0.52-0.4520.85-3.42-0.32-1.2633.31-0.29-0.41-0.2941.633.44-0.30-1.40开挖11-8.13-1.10-0.04-0.302-0.99-3.55-0.06-0.9030.37-0.39-0.02-0.424-0.403.71-0.05-1.01衬砌11-7.89-1.06-0.06-0.442-0.91-3.47-0.09-0.9530.56-0.50-0.04-0.4740.453.45-0.08-1.01结果分析调压井有限元分析结果与监测成果基本相符，表明模型合理地模拟了复杂地形对结构体的影响。衬砌的分析表明其内力分布不均匀，局部应力较大。计算显示后边坡侧的变形较大，应该加强锚固，进行变形监测，以防边坡滑移和塌落影响调压井的正常使用。通过“二次开挖、衬砌”与“一次性开挖、衬砌”的计算结果对比分析,表明分层开挖时围岩应力、应变变化幅度相对较小,有利于围岩稳定。根据现场监测数据，进行的材料弹性模量的反演计算，所确定的参数能基本反映材料的性能。有一定的可靠度。为下一步多参数反演分析打下了基础。工作展望结构数值分析：利用学术、技术和软件优势，为实际工程服务，解决各类结构计算问题与过程模拟仿真。结构类型：建筑结构、桥梁结构、岩土工程、道路工程、地下工程、水工结构等