闵浦二桥索塔锚固区计算报告

闵浦二桥索塔锚固区计算报告北京迈达斯技术有限公司2010．02midasFEA闵浦二桥索塔锚固区计算报告1一、模型介绍采用有限元软件midasFEA分析，有限元模型按实际尺寸1：1建立，模型全部采用8节点六面体网格划分，保证网格节点的耦合。在底面上施加竖向约束，横桥向约束中间部位个别节点，顺桥向约束中间部位个别节点，保证约束处的横桥向和顺桥向反力为0，以接近真实情况。模型网格如图1：图1模型网格图预应力筋布置与模型设计图一致，采用植入式钢筋单元的方式，每束预应力筋按设计施加预应力，不考虑自重，索力按照面压力均布加载到钢垫板上，如图2：midasFEA闵浦二桥索塔锚固区计算报告2图2预应力布置示意图二、预应力和顶推荷载共同作用下结果（a）midasFEA闵浦二桥索塔锚固区计算报告3（b）图3索孔下出口处横桥向应力分布图（a）FEA结果（b）ansys结果从图3可以看到，测点附近区域的横桥向应力在-1.38～-8.20MPa之间，而实测横向应力在-2.15～-9.67之间，两者基本上一致。相应ansys计算结果为0.01～-12.1MPa。（a）midasFEA闵浦二桥索塔锚固区计算报告4（b）图4索孔下出口处竖向应力分布图（a）FEA结果（b）ansys结果从图4可以看到，测点附近区域的竖向应力在1.12～-1.82MPa之间，而实测竖向应力在1.72～-1.78之间，两者基本一致。相应ansys计算结果为1.34～-2.29MPa。图5顶面横桥向应力分布图从图5可以看到，测点附近区域横桥向应力计算值在-5.50～-11.19MPa，与实验实测值-7.01～-12.49MPa相比基本一致。相应ansys计算结果为-6.29～-12.8MPa。midasFEA闵浦二桥索塔锚固区计算报告5图6索孔下出口处主拉应力分布图（a）FEA结果（b）ansys结果从图6可以看到索孔下出口附近的主拉应力范围在1.65～-0.46MPa，与实验实测值1.74～-2.97MPa相比基本一致。相应ansys计算结果为1.59～-0.52MPa。三、仅在预应力作用下结果midasFEA闵浦二桥索塔锚固区计算报告6（a）（b）图7索孔下出口处横桥向应力分布（a）FEA结果（b）ansys结果从图7可以看到，测点附近区域的横桥向应力在-4～-12.58MPa之间，而实测横向应力在-4.01～-14.85之间，两者基本上接近。相应ansys计算结果为-1.06～-10.28MPa。midasFEA闵浦二桥索塔锚固区计算报告7（a）（b）图8索孔下出口处竖向应力分布（a）FEA结果（b）ansys结果从图8可以看到，测点附近区域竖向应力计算值在0.74～-2.88MPa，与实验实测值0.93～-2.40MPa相比基本一致。相应ansys计算结果为-0.04～-2.84MPa。midasFEA闵浦二桥索塔锚固区计算报告8图9顶面横桥向应力分布图从图9可以看到，测点附近区域横桥向应力计算值在-4.56～-13.36MPa，与实验实测值-9.14～-14.61MPa相比基本一致。相应ansys计算结果为-2.65～-16.47MPa。四、仅在顶推作用下结果图10索孔下出口处横桥向应力分布从图10可以看到仅在顶推荷载作用下测点附近横桥向应力分布范围在midasFEA闵浦二桥索塔锚固区计算报告91.49～5.57MPa，与实测范围1.28～5.56MPa基本一致。相应ansys计算结果为0.77～5.55MPa。图11索孔下出口处竖向应力分布从图11可以看到仅在顶推荷载作用下测点附近竖向应力分布范围在-0.40～1.57MPa，与实测范围-0.38～2.94MPa基本一致。相应ansys计算结果为-0.57～1.73MPa。图12顶面横桥向应力分布图从图12可以看到仅在顶推荷载作用下顶面测点附近横桥向应力分布范围在midasFEA闵浦二桥索塔锚固区计算报告10-3.57～3.05MPa，与实测范围0.10～2.12MPa基本一致。相应ansys计算结果为-4.88～3.33MPa。五、裂缝开展模拟结果对该模型进行非线性加载模拟，本构模型采用FEA专为混凝土材料提供的总应变裂缝模型。总应变裂缝模型考虑了裂缝的产生对刚度的影响，并且考虑了三维横向裂缝和横向约束的影响，广泛应用于正常使用极限状态和承载能力极限状态的计算，可以良好地模拟裂缝的影响和发展。本模型先一次性加载预应力后，对模型加载顶推荷载，荷载大小为3P，荷载步分为20步，程序在算到第18步，也就是顶推荷载加载到0.9*3P=2.7P时由于局部裂缝的影响很难收敛了，以下是每个荷载步裂缝的发展，为避免预应力张拉端和沿途应力集中的影响，仅截取顶推面查看裂缝的分布。1.05P步1.8P步1.95P步2.1P步midasFEA闵浦二桥索塔锚固区计算报告112.25P步2.4P步2.55P步2.7P步图13各个荷载步裂缝应变纵向视图1.95P步横向1.95P步斜向midasFEA闵浦二桥索塔锚固区计算报告122.25P步横向2.25P步斜向2.55P步横向2.55P步斜向2.85P步横向2.85P步斜向图14各个荷载步裂缝应变横向视图及斜视图从裂缝区域的发展可以看到，加载至1.95P之前，裂缝仅在索孔管道两侧和midasFEA闵浦二桥索塔锚固区计算报告13预应力张拉沿线有零星的分布，法向裂缝应变在10e-5量级，可视为无法观测到的微裂缝。加载到1.95P的时候，索孔下出口上缘处（图中红圈标注）出现密集的裂缝分布并向上方扩展，法向裂缝应变达到10-4量级（颜色明显比其他区域深），而索孔管道两侧微裂缝应变大小和范围没有发展，根据法向裂缝应变的大小，结合单元在裂缝法向的宽度，裂缝宽度约在10e-2mm量级，可以认为在索孔下出口上缘裂缝开始向上发展。这与试验测得在2.0P出现裂缝基本一致。加载到2.1P时，索孔下出口下缘出现了密集的裂缝分布并向下方扩展，法向裂缝应变在10-4量级（颜色明显比其他区域深），可以认为出现了裂缝。两处裂缝都是略朝Y轴正方向发展，这也与试验观察到的裂缝2，裂缝3的发展方向一致。图中蓝色圆片代表的裂缝分布由于法向裂缝应变均在10-5量级以下，虽然区域较大，但均为混凝土内部微小裂纹，不能视为宏观裂缝。此外索孔区域以外的局部裂缝分布或是由预应力张拉带来的应力集中造成的，或是荷载加载部位的应力集中造成的，是有限元程序难免，均应忽略。实验所得的A、B两侧裂缝发展不对称，是由实验的误差造成的，而有限元网格和边界荷载均对称，所以两侧裂缝的分布和发展呈相同状态。六、结论经过与ansys计算结果对比，可以发现FEA由于采用植入式钢筋方式模拟钢筋，预应力张拉造成的应力集中现象相比ansys小通过对索孔下出口和顶面的横桥向应力、竖向应力、主拉应力对比，两种模拟的应力分布范围和趋势基本一致。与实测结果相比，FEA计算的应力范围与实测结果基本一致，而且相比ansys结果更接近实测结果一些，由于网格点位置和实际测点位置可能有些差别，读取的计算结果直接和实测点结果相比会有部分误差。另外裂缝产生的部位和裂缝随加载发展的状态也与实验结果比较一致，可见总应变裂缝模型能对裂缝的发展进行有效的模拟，得到裂缝分布状态和发展方向，具有现实的指导意义。midasFEA“AdvancedNonlinearandDetailAnalysisSystem”