[发明专利]一种预应力束与混凝土间粘结性能的有限元模拟方法

摘要

本发明公开了一种预应力束与混凝土间粘结性能有限元模拟方法，提出了预应力束粘结滑移等效模拟方法。首先，给出了预应力束整体滑移面及有效工作面确定原则，建立了预应力束界面微单元受力平衡关系，提出了预应力束极限粘结力计算方法；其次，给出了预应力束极限粘结力的等效原则，并据此推导了单根预应力筋等效最大粘结应力计算公式，更新了BPE粘结滑移模型，明确了模拟粘结滑移的弹簧单元的属性，实现了单根预应力筋等效粘结滑移有限元模拟；在此基础上，基于单根预应力筋等效粘结力叠加获得预应力束粘结力、单根预应力筋平均滑移作为预应力束滑移的原则，实现了预应力束粘结滑移的有限元模拟。最后，结合现有试验结果验证了本发明方法的精度。

主权项

一种预应力束与混凝土粘结性能的有限元模拟方法，其特征在于包括以下步骤：步骤1：建立几何模型；根据实际工程结构，确定模型几何参数，混凝土与预应力孔道中的压浆体的几何模型通过采用分离式建模方法获得；预应力束几何模型的建立可通过两步实现，首先采用实体切分法在钢绞线的位置将混凝土切分，然后在切分处重新创建多根与钢绞线几何参数相同的线作为预应力束几何模型，普通钢筋的几何模型采用同样的方法得到；模型中不考虑混凝土与压浆体间的粘结滑移，视两者为完全粘结，忽略普通钢筋与混凝土间的粘结滑移，在两者间采取建立约束方程的方法连接；步骤2：确定预应力束整体滑移面及有效工作面的构成；具体包括：钢绞线与周围混凝土接触面、钢绞线间混凝土浆体的有效工作面积和粘结力构成；步骤3：确定预应力束极限粘结力及预应力束中单根预应力筋等效最大粘结应力；建立预应力束界面微单元受力平衡关系，提出预应力束极限粘结力计算方法；在此基础上，给出预应力束极限粘结力的等效原则，并据此推导单根预应力筋等效最大粘结应力；步骤4：定义单元类型与材料本构关系；1)混凝土与压浆体材料采用三维实体单元Solid65模拟，并采用Hongnestad材料本构模型；2)预应力束中钢绞线采用Link8单元模拟，普通钢筋同样采用Link8单元模拟，钢绞线与普通钢筋材料采用简化的Menegotto本构模型；3)预应力束中钢绞线与混凝土间的粘结滑移关系采用COMBIN39非线性弹簧单元模拟，此单元属性需通过钢绞线与混凝土单元间粘结滑移(F‑S)关系定义；步骤5：定义预应力束中单根钢绞线与混凝土间弹簧单元属性；首先，基于预应力束中单根钢绞线等效最大粘结应力τmax计算公式，更新BPE粘结滑移模型，推导出钢绞线单元与混凝土单元间的粘结滑移(F‑S)关系，获得模拟单根钢绞线与混凝土间粘结滑移的弹簧单元的属性，进而实现预应力束中单根钢绞线等效粘结滑移有限元模拟；单根钢绞线与混凝土间的粘结滑移采用COMBIN39弹簧单元进行模拟，具体是在钢绞线与混凝土单元的重合节点间采用三根长度为零的弹簧单元连接；在三根弹簧中，法向与横切向的弹簧变形相对于纵切向的变形可忽略不计，因此，可以将这两个方向弹簧的刚度系数K取为无穷大；钢绞线与混凝土间的粘结滑移主要通过纵向弹簧模拟，其单元属性由单根钢绞线与混凝土单元间粘结滑移(F‑S)关系获得，其中单根钢绞线与混凝土单元间粘结力F通过下式计算：F＝Ar·τ＝n·br·lr·τ      (1)式中，Ar为预应力筋单元表面积，br为预应力筋单元直径，lr为预应力筋单元长度，τ为预应力筋与混凝土单元间粘结应力，可根据BPE粘结滑移模型确定：τ=τmax×(ss1)α0≤s≤s1τmaxs1≤s≤s2τmax-(τmax-τ3)(s2-s)s2-s3s2≤s≤s3τ3s>s3---(2)]]>式中，τmax为预应力束中单根钢绞线等效最大粘结应力，S1、S2、S3、S4、τ2、α具体参考文献“Analysis of the bond behavior between prestressed strands and concrete in fire,Jamal Khalaf,Zhaohui Huang,Construction and Building Materials,128(2016)12‑23”。步骤6：网格划分与加载方式的定义；①网格划分，网格的质量将影响到计算精度与收敛性，尽量采用规则六面体单元，对于非线性响应强烈的区域，为避免应力集中而提前开裂，单元尺寸不宜过小；②定义模型的加载方式，模型中采用分级加载，荷载施加于预应力筋端部节点；模型中应避免将约束直接施加在混凝土节点上，故在模型两侧各添加一块刚性垫板用于施加约束；步骤7：有限元求解及后处理；①求解，定义分析类型与求解控制选项，设置平衡迭代最大次数、收敛准则等参数进行有限元求解计算；②提取结果及后处理，提取各级荷载模型下的剖面应力云图，提取每根钢绞线的荷载滑移曲线，并取其平均值作为预应力束的荷载滑移曲线；步骤8：确定预应力束粘结滑移曲线并验证其精度。预应力束粘结力通过单根预应力筋等效粘结力叠加获得，预应力束滑移取单根预应力筋平均滑移值；结合已有实验数据验证本发明方法的精度。