1、陕西理工学院毕业设计第 页 共 页外文资料原文陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页陕西理工学院毕业设计第 页 共 页译文火灾时一般框架的简化弹塑性分析摘要：针对高温下的三维框架的分析，本文提出了一种可供选择的简化方法。在塑性发展过程中，考虑垂直力及弯矩作用，则一般塑性铰原理包括了温度的影响。利用刚度的折减系数，来改善附加的允许的单元逐步屈服。同时讨论了如何确定内力矢量，其考虑了非线性热效应。利用实例反映的不同情况，对本文建立模型
2、与实验结果进行了比较。由此可知，本文建议方法可用于分析实际火灾中的二维及三维框架，同时其精度合理，计算成本低。关键词：弹塑性分析；钢结构；三维结构；塑性铰；非线性分析；SAFIR；有限位移1. 引言在发生偶然火灾时，为了确保生命安全和减少经济损失，在一定时间内，结构必须保持稳定性。如果是钢结构，细长构件的高导热性及低热惯性促进了温度的升高，而温度的膨胀和高温下材料性质的变化将导致塑性变形和大位移，最终导致整个结构破坏。因此，对钢结构的防火安全来讲，一个重要部分是研究火灾时结构的行为分析。由于它的复杂性，不能利用手工计算来估算整体结构的行为。另一方面，实验室试验受到操作性和经济方面的限制作为自然
3、影响因素，促进了基于优先元法的数值模型的发展。为了研究火灾时的钢结构，一般采用塑性区法。关于此模型的发展研究可以参阅文献 28及Friedmen 的研究。关于这方面，Friedmen 曾给出了综述，而 Olenick 和 Carpenter21进行了发展。这些研究中的一个现代化的计算规范是 SAFIR，比利时国王大学 Sanad 等人 26的研究表明，一些特殊实例可以采用非线性结构分析的商业化程序，如 ABAQUS1。本文研究的目标是提出一种可供选择的方法。习惯上将塑性铰的基本理论用来分析高温状态下的非弹性结构，但在本文中被延伸到了分析热效应。附加的改进之处包括模拟应力及屈服的分布。我们已经提
4、出了框架的塑性铰模型。Liew 等人 17已经进行了分析，但是没有给出具体的公式。Li 和 Jiang15给出了在单元内力变化为线性时的热膨胀效应模型。Toh 等人 30对塑性铰法和塑性区法的分析结果进行了比较。Lu 光和 Chan10给出了基于塑性铰的模型，此模型在刚度矩阵中考虑了热效应。在下部分中我们将介绍包括冷却相变的改进模型，所有这些模型都是二维的。在本文中，公式均延伸到三维的情况。同时，传统的塑性铰模型是精确的，它通过单元模拟逐步屈服过程。这种方法的发展是基于以前的公式 5，19 ，它的运动理论是基于 Oran22的研究基础之上的，Oran 容许在小应变的情况下，产生大位移。Argy
5、ris 2假设几何刚度矩阵包括切线弯矩，尤其是对于空间结构。本文在 2.7 节讨论了内力的计算，此内力考虑了温度效应。在动力荷载作用下，利用Yang33的研究方法，在临界状态后，随着荷载的降低，位移增大。Newton-Raphson 判断了在恒定荷载作用下，热量增加的位置。通过与精确模型以及实验结果相比较，证明了此数值方法的可行性。2. 公式2.1 基本假设数值程序是建立在以下假设基础之上的：（1）结构承受具有小应变的有限位移。（2）最初单元是垂直、等截面、双轴对称的横截面。（3）根据平截面仍保持平面以及它仍同纵轴正交，进行了伯努力-欧拉假设，忽略了剪切变形和翘曲变形。（4）假设随着周围环境温
6、度以及高温作用，单元可提高其整体塑性弯矩，而不是局部屈服的性能，陕西理工学院毕业设计第 页 共 页同时单元具有适量的转动能力。（5）作用在单元结点的荷载（径向荷载）成比例地增加，同时形状独立。在热增量分析中，荷载保持不变，热量增加。（7）一般塑性铰理论中 18包括正应力的影响，它包括港督的逐渐减小。（8）采用塑性流变规则。对于材料的性能，当热量比在 2-500C/分钟时，采用欧洲规范 34的应力-应变关系是可以的，它包括蠕变效应。对于塑性铰，如图 1 所示，近似采用双线性模型，同时采用欧洲规范推荐的高温时的屈服强度极限及弹性模量的折减息书，热膨胀系数恒等于 14 10-6/0C。图 1 高温下
7、钢材的应力-应变关系2.2 运动方程Oran2提出在运动时，当小应变发生较大位移时，即使在粗大网格中也可以提供较好的结果。为了描述运动学原理，随后利用 3 3 阶矩阵进行了解释，如图 2 所示。=结点修正矩阵：柱子时坐标轴 和 轴的方向余弦，其与端点 i（i=1，2）刚接，同时在开始时与世界坐标系 X，Y ，Z 平行。=端部修正矩阵：柱子是结点 i 方向余弦的垂线，同时沿着截面的主要方向。=单元转动矩阵：行是局部坐标轴 x，y，z 的余弦。X 轴通过截面的重力形心，y 轴及 z 轴平行于中部垂直方向的主要方向。在结点转动中，适量考虑几何非线性的影响时重要的。可以通过一般迭代长度（L t）和初始
8、长度（L 0）来确定轴向相对位移 Ux ，通过相对广义位移来定义单元变形运动学的优点是可以反映杆件端截面的转动。在计算正应力及剪应力合力方面，利用一般纵轴的弦长可能引起错误的结果，一般它在屈曲状态方面出现错误。利用本文的运动学，它的每个单元的端部在参考系统中是独立的。因此即使一个单元仅仅划分若干个区格，也可以提供精确的结果，文献 19， 22，28 给出了更为详细的解释。2.3 塑性铰模型利用梁-柱有限元，钢框架的非弹性分析常用二种形式：塑性区和塑性铰。塑性区分析又称为分布或分散塑性，同时可以模拟通过横截面和沿着杆件长度方向的逐步屈服。这个模型也很容易考虑热量梯度。另一方面，与塑性铰法相比，塑
9、性区法计算复杂，成本高。在传统的塑性铰或集中的塑性铰模型中（通常用于钢结构的塑性设计中） ，在单元中要保持塑性，仅在端部考虑塑性效应。一旦弯矩达到塑性铰，塑性铰理论就认为完全塑性。在许多情况下，这样的模型完成的很好，与精确塑性区法相比，可以为整体框架的非线性刚度陕西理工学院毕业设计第 页 共 页和强度提供预测 14。对于单个独立的单元和框架，结构失效试验的最主要方面是沿着长度方向的屈服分布，因此有必要改进弹-塑性铰模型 9。在弹- 塑性铰模型中，在预测非弹性的荷载-位移曲线和强度方面进行修正，可以引起明显的改进。图 2 （a ）单元的世界及局部坐标轴和初始参考轴；（b）各截面的法向及基本方向2
10、.4 屈服函数在近期的研究中，采用传统的塑性铰理论（考虑屈服弯矩和垂直力之间的耦合作用） ，对矩形截面的梁柱，沿长轴方向的塑性弯矩 Mpz 以及轴向屈服荷载 Np，可用文献 19的屈服函数表示式中 和 ，分别表示轴向荷载与屈服荷载的比值以及长轴屈服极限弯矩与塑性弯矩的比值。对于平面框架的 I 字型截面，可利用文献 12表示屈服函数对三维框架，屈服面为 23式中 my=My/Mpy 表示短轴屈服极限弯矩与塑性弯矩的比值。同时，下一部分将进行更为详尽的描述，目的是获得更实际的反应，它是沿着单元建立内应力与屈曲分布的反应。2.5 残余应力以及屈服分布利用修正系数 和 来考虑这二方面的效应。为了考虑应
11、力-应变关系中热轧型钢梁中的残余应力效应（如文献 3） ，可用有效模量代替弹性模量式中 为材料刚度的这算系数，计算见文献 6考虑到塑性发展沿着单元的塑性分布，文献 12给出了附加的刚度折减系数用刚度折减系数 ，切向弯曲刚度从最大值逐渐降低到 0，在 0.5-1.0 之间变化。同时，因为弯矩沿着单元并不相同（参见文献 34） ，也就是说，计算每个单元端部的折减系数陕西理工学院毕业设计第 页 共 页，同时通过折减系数 、 来计算扩大刚度矩阵。至于切线模量，在 Np、M py 和 Mpz 中利用屈曲强度折减系 Ky,T 考虑。2.6 刚度矩阵为了得到适用于三维框架的表达式，Argris 确定弯矩是准
12、切线，同时扭矩是半切线，并提出了他们之间的关系。这种关系是用来在维持节点转角和弯矩一致的情况下，推导采用 Argris 提出的单元刚度矩阵。这个刚度矩阵包括直接推导的几何非线性部分，因此可通过数值积分来提高计算效果。2.7 内力矢量作为非线性分析，常采用增加塑性（参见文献 18的例 2）和以下假设来进行（与弹塑性及温度作用有关）：（1）弹性增量 ，塑性增量 和热位移增量（2）存在与应力合力的屈曲标准有关的温度式中 是应力合力矢量，T 是温度。（3）协调条件：在塑性变形时，合力必须保持在屈服面，（4）内力和位移之间的弹性关系为式中 是总位移中的弹性部分， 是切线刚度矩阵（利用 Argris2几何
13、非线性效应的函数） ，内力增量可表示为（5）根据常规条件确定的塑性位移增量为式中 是一个正的标量和 是屈服函数的梯度。将式（13）代入式（8） ，整理后将 代入式（12） ，得将式（10）代入即得陕西理工学院毕业设计第 页 共 页因此，将式（14）代入即式（17）可写成式（19）表达了单元内力的增量，它考虑了几何非线性效应、热膨胀和屈服温度，若仅考虑动力荷载，式（19）可简化为它与文献 5的结果相符。在反复荷载作用下，位移增量迭代为3. 计算方法这个计算方法采用一个与迭代增量解答有关的修正拉格朗日表达式 22，28 。采用控制荷载锃亮的分析方法，在参考文献 24中可以看到有关这方面的更为详细地
14、介绍，它提供了超过临界状态，甚至是软化行为的荷载-位移路径。为了分析热增量，采用 Newton-Raphon 方法。在这种情况下，开始施加动力荷载在建筑物上，同时找到其相应的平衡位置。然后，在荷载不变的情况下，通过结构单元来介绍温度增量。通过式（19）确定每个单元的内力增量，同时组合整个结构的内力到矢量。通过内力与外力的不同，给出不平衡力。根据相应的不平衡合资阿，建立结构位移增量方程这些位移增量被增加到总位移，同时修正建筑物的形状，利用式（22）和（23）逐步迭代，直到内力及外力之间达到平衡。然后，在结构上施加一个新的热增量，重新开始计算。4. 实例4.1 实例 1图 3 给出了出自文献 23
15、的门式刚架。考虑到材料的弹塑性，取弹性模量为 200000MPa，屈服强度为 250 MPa。梁-柱横截面为矩形截面，面积为 31612mm2 和惯性矩等于 8.328 107mm4。式（2）表示矩形截面，塑性铰得到的屈服函数如图所示。框架分为 5 个单元。利用 ABAQUS1应用塑性区法计算。图 4 给出了结果。有图中可知，与粗分法相比，本文公式结果吻合很好。为了与 ABAQUS 结陕西理工学院毕业设计第 页 共 页果相比，采用 32 个单元（每个柱 10 个，梁 12 个单元） 。4.2 实例 2Rubert 和 Schaumann 对一系列平面框架进行了试验，其中之一是如图 5 所示的门
16、式刚架。它是利用电力装置进行均匀加热，所有截面均为 IPE80，弹性模量及屈服强度的初始值分别为210000MPa 和 382MPa。同时利用式（3）的屈服标准，图 6 给出了试验和图 5 所示的计算结果的比较。图 3 例 1 中用到的 5 个单元的框架图 4 室温时框架的荷载-位移曲线图 5 例 2 中用到的单层单跨框架4.3 实例 3同以前实例25一样，进行了同样系列试验。图 7 给出了单层双跨框架局部加热的情况。在前陕西理工学院毕业设计第 页 共 页面实例中，所有截面均为 IPE80，在室温时材料屈服强度为 355 MPa， ，弹性模量为 210000MPa。如图 8 所示。模拟及试验结
17、果吻合很好。图 6 单层单跨框架的位移-温度曲线图 7 例 3 中用到的单层双跨框架图 8 单层双跨框架的位移-温度曲4.4 实例 4Mcnamee 和 Lu 对图 9 所示的三层框架进行了试验。为了研究温度作用，考虑动力侧向荷载，即 P=30KN，在二层模拟火灾时的热曲线。同样，式（3）作为屈服标准，图 10 给出了本文研究和SAFIR 的比较。尽管低估了位移，但是临界温度仅相差 0.8%（本文 693.70C，SAFIR 为 689.00C） 。4.5 实例 5图 11 给出了出自文献 13的三维框架。荷载 P=250KN，材料弹性模量初始值为 210000MPa，屈服强度为 325MPa
18、，所有截面均为 H150 150 7 10。通过式（14）的屈服条件进行计算，图 12 陕西理工学院毕业设计第 页 共 页给出了在荷载 H1、H 2 点之间的水平位移，同时与 SAFIR 结果进行了比较。同前述实例一样，高估了塑性铰模型的荷载路经，破坏温度之间仅仅相差 4.1%（本文 5120C，SAFIR 为 5330C） 。图 9 例 4 中用到的三层单跨框架图 10 三层单跨框架的位移-温度曲线图 11 例 5 中的三维框架陕西理工学院毕业设计第 页 共 页图 12 三维框架的位移-温度曲线5. 结论本文提出了一种适用于火灾中的三维钢框架分析的改进简化方法。这种模型是建立在一般塑性铰原理
19、基础之上的。它包括在法向内力和弯矩共同作用下，通过横截面的塑性铰函数。它主要在以下几个方面进行了改进：（1）利用刚度折减系数考虑了单元内应力和屈服的逐步分布，同时在计算效应方面有所增加。（2）在每个单元端部，采用单独的系统描述运动特性，同时采用端部近似转动来表示，如图 1 所示。在有限位移效应方面，Argris 提出的刚度矩阵可以提供精确解答。（3）内力矢量（应仔细详细地进行确定）考虑屈服级温度变化的组合效应。（4）用于分析荷载增量的控制方法，在超过临界状态时可以确定平衡路径，同时在荷载降低时，位移增加。通过试验结果与先进规范的比较，如 SAFIR 给出的结果，二者吻合很好。由结论可知，本文建议改进的一般塑性铰方法，在温度升高时可以用于分析平面及空间框架，同时其精确度合理，计算效率较高。
