火灾作用钢筋混凝土柱非线性有限元分析.docVIP

火灾作用钢筋混凝土柱非线性有限元分析摘要；利用ANSYS软件热-结构耦合功能，模拟分析了以往做过的试验，探索ANSYS软件进行结构火灾分析的技巧与优缺点。模拟分析了构件内部温度场和单调加载下钢筋混凝土柱的荷载-位移曲线，温度场数值分析结果与试验结果吻合较好，荷载-位移曲线与试验结果稍有差别，数值计算能够满足工程精度要求。 关键词；热-结构耦合分析;荷载-位移曲线;数值分析;配箍率;轴压比 火灾是各种灾害中发生频率最大且极具毁灭性的灾害之一，常见的火灾有建筑火灾、森林火灾、可燃材料堆场火灾、交通工具火灾等，火灾损失统计表明，发生次数最多、损失最严重者当属建筑火灾[1]。研究火灾作用下结构或构件的性能及火灾作用后结构或构件的损伤评估与修复是目前极为迫切的要求。由于受试验、经费等条件的限制，目前国内外常采用有限元数值分析的方法。本文用ANSYS有限元软件进行热-结构耦合分析，与相应的试验结果进行对比分析，探索ANSYS软件进行结构火灾分析的技巧与优缺点，并设计了数值试验试件，分析了配箍率、截面尺寸和轴压比对火灾作用后钢筋混凝土柱抗火性能的影响。 1 ANSYS热-结构耦合分析基本原理 热-结构耦合分析是指考虑温度和力两种不同物理状态荷载之间相互作用的分析，在ANSYS中通常用两种方法进行耦合分析：直接耦合和间接耦合（顺序耦合）。顺序耦合分析法包括两个或多个物理过程按一定顺序排列的分析，每一种分析属于某一个物理场分析，通过将前一个分析的结果作为荷载施加到下一个分析中的方式进行耦合，特点是各种分析相互独立。如在热-结构耦合分析中是先进行热分析，然后将热分析结果中得到的节点温度作为“体荷载”施加到随后的结构分析中去。直接耦合分析是指仅对耦合场进行一次分析，使用包含多场自由度的耦合单元，通过计算包含所需物理量的单元矩阵或载荷向量的方式进行耦合。如使用SOLID45，PLAND13或SOLID98单元的压电分析。 本文在进行热分析时混凝土和钢筋单元分别选用了热单元SOLID70和LINK33，在结构分析时要将热单元转换成结构单元，ANSYS程序默认状态下是将热单元SOLID70和LINK33分别转换成结构单元SOLID45和LINK8。 2 用ANSYS软件进行热-结构耦合分析可行性验证 本文选取文献[3]与[4]三面高温的偏心受压柱的试验研究为对象，针对其三面高温轴心受压的试验过程，用ANSYS9.0程序进行模拟分析。如图1，混凝土试件实测立方体抗压强度 =37.5MPa，棱柱体强度 =28.5MPa，纵筋常温下屈服强度 =340MPa，箍筋屈服强度 =256MPa。 2.1 ANSYS温度场数值计算值与试验值比较 按文献[3]试验升温曲线升温，考虑结果具有一般性，抽取升温时间为110分钟，加热温度800℃时的数值计算结果与试验数据进行比较，如图2和图3所示。 数值计算截面温度沿宽度和高度分布与试验实测温度分布吻合较好，温度分布的梯度也同样很接近，数值计算截面最高温度741.115℃与试验实测温度约为750℃相接近，非加热面数值计算最低温度为300℃与图3试验实测相同部位温度（约为280℃）略偏高，这主要是因为数值计算没有考虑试验时加热炉本身的散热和外部环境影响下其他方面的散热。 对比可见，两者结果接近，由于数值分析不能完全模拟试验情况，没有考虑钢筋与混凝土间的粘结滑移，数值计算的极限承载力比试验结果偏高，而极限挠度偏小。图4和图5将试验实测荷载-挠度曲线与数值计算的荷载-挠度曲线进行了比较，二者极限承载力和极限挠度差别不大，荷载-挠度曲线形式大致相同，自由升温阶段侧向挠度吻合较好。由于数值分析中钢筋与混凝土有限单元间没有考虑粘结滑移，致使混凝土数值分析曲线在屈服后比试验曲线较平缓。 通过以上对比分析可见，ANSYS软件模拟计算构件温度场精度较高，结构数值计算与试验数据比较吻合，能够满足工程精度要求。 3 总结 通过用ANSYS软件热-结构耦合模拟分析钢筋混凝土轴心受压试验，可以得出，ANSYS软件热分析计算出的构件温度场与试验值吻合很好。结构分析时，自由升温阶段变形与试验值吻合较好，在荷载作用下，所得出的极限荷载和变形与试验值相近，得出的荷载-位移骨架曲线下降段与试验值无明显差别。 利用ANSYS软件热-结构耦合模拟分析，可以应用到结构数值分析中。 注：本章论文的所有图表及公式以PDF形式查看 1