基于有限元分析的钢筋锈胀时间与钢筋锈蚀率的计算.docxVIP

基于有限元分析的钢筋锈胀时间与钢筋锈蚀率的计算 混凝土表面的锈蚀期通常被认为是混凝土结构构件正常使用的极限状态的象征之一。许多国内外学者对此进行了研究。alono等人对混凝土层的锈腐烂率进行了试验分析。在此基础上，使用有限元法研究了不同实验部分中不同元素的混凝土层的腐蚀效果。屈文军等人利用温度膨胀环模拟了钢筋的腐蚀产物，并建立了有限分析模型混凝土腐蚀断裂。尹飞等人也进行了相关实验分析。根据试验结果和理论分析，确定了混凝土层的均匀氮侵蚀率，但缺乏广泛适用性。赵玉伟等人使用力学方法分析了混凝土层的裂缝过程。力学分析是基于材料的均匀性，因此有一定的局限性。 数值分析方法可以解决许多试验无法或难以解决的问题,是分析和研究混凝土结构性能的有力工具.随着对材料性能的深入了解,材料本构关系的正确揭示,以及计算机软件和硬件设备的发展,数值分析方法已经成为混凝土结构分析的重要方法. 本文采用数值分析方法确定混凝土构件锈蚀胀裂时刻的钢筋锈胀压力,包括均匀钢筋锈胀压力和非均匀钢筋锈胀压力的情况;然后根据作者先前建立的混凝土保护层胀裂时刻钢筋锈胀压力和钢筋锈蚀率的相关关系,以及本文数值分析的结果,提出了混凝土保护层胀裂时刻非均匀钢筋锈蚀率的计算方法,这对混凝土结构耐久性研究是非常有意义的. 1 均匀钢筋支护压力的有限分析 1.1 有限元模型建立的背景 本文的分析采用ANSYS有限元计算软件,对混凝土构件截面受钢筋锈胀压力的情况进行分析计算.在ANSYS有限元程序的基本单元中,SOLID65单元可以描述混凝土单元开裂的情况,因此SOLID65经常会被应用在混凝土构件受力并开裂的情况分析中.由于本研究关注混凝土在受锈胀压力而开裂的现象,为了模拟混凝土受锈胀压力开裂的情况,本文采用了SOLID65单元进行分析.根据SOLID65单元的特点,模型分析中本构关系和破坏准则分别采用Drucke-Prager模型和WillamWarnke五参数破坏准则,裂缝为弥散固定裂缝模型(smeared fixed crack model),并采用位移收敛准则. 本文在建模时,网格划分由软件自动生成,但在钢筋边缘加密了混凝土的单元网格的划分,以便更好地分析混凝土保护层锈胀开裂的情况;构件远端由于对分析结果影响甚小,采用较大网格划分,以节省计算所需空间与时间.有限元分析计算的网格划分如图1、2所示.图1为混凝土梁截面角部混凝土受钢筋锈胀压力的分析模型,采用C25混凝土,其有限元分析的输入参数分别为:受拉强度ftk=1.75 MPa, 受压强度fck=17.5 MPa,弹性模量Ec=11.9 MPa,泊松比νc=0.2.有限元分析模型的具体尺寸为:200 mm×400 mm,钢筋保护层厚度为25 mm,钢筋直径为20 mm.由于本模型采用二维分析,未考虑混凝土梁的跨度取值问题;约束条件为梁顶双向约束,即X与Y方向均为固端支座.图2为混凝土梁截面中部混凝土受钢筋锈胀压力的分析模型,其他条件同图1模型.混凝土受均匀钢筋锈胀压力的情况如图3所示,图中的q为混凝土受到的均匀钢筋锈胀压力. 1.2 混凝土开裂情况 根据有限元计算分析,当角部钢筋锈蚀产生的均匀钢筋锈胀压力达到1.45 MPa时,混凝土构件承受的钢筋锈胀压力达到极限值,即环向的钢筋锈胀压力已经加到最大值而无法继续加载,外围混凝土的开裂情况如图4所示;当中部钢筋锈蚀产生的均匀钢筋锈胀压力达到1.52 MPa时,混凝土构件承受的钢筋锈胀压力达到极限值,外围混凝土的开裂情况如图5所示.中部钢筋锈蚀情况的计算结果大于角部情况的计算结果,这是因为中部钢筋锈蚀产生钢筋锈胀压力有更多外围混凝土约束的缘故. 以上对角部和中部的钢筋锈蚀情况分别作了分析,其分析方法和计算结果对于混凝土胀裂时刻钢筋锈胀压力的确定是有意义的;而在对实际混凝土构件进行锈胀压力分析时,可以按照混凝土构件的实际几何尺寸建模进行钢筋锈胀压力作用的分析,这里举一例示意. 某混凝土梁,梁截面尺寸为200 mm×400 mm,采用C25混凝土,其分析参数同上述,梁底部设有3根直径为20 mm的钢筋,钢筋保护层厚度为25 mm.该梁在均匀钢筋锈胀压力为1.35 MPa时,达到钢筋锈胀压力极限值,计算结果如图6所示.按混凝土构件实际情况进行分析时,其钢筋锈胀压力极限值(1.35 MPa)略小于单独考虑角部(1.45 MPa)和中部(1.52 MPa)钢筋锈蚀的情况,这是因为本模型钢筋之间的距离较小,钢筋之间的混凝土裂缝贯通连接而使钢筋锈胀压力达到极限值的缘故. 2 非均匀钢筋混凝土结构力的有限分析 2.1 钢筋非均匀治理机理 在实际工程中,由于有害介质的单向侵蚀路径,大量的钢筋锈蚀情况属于非均匀锈蚀,如图7所示,其非均匀锈蚀钢筋净截面Art可按下式计算: Art=A?2∫α012(