【JY】ETABS纤维铰的使用
ETABS中框架单元的非线性行为可以通过塑性铰和纤维铰两种方式实现，塑性铰是通过建立力与变形关系，从整体上描述截面的非线性行为，而纤维铰则是通过将截面离散为有限根纤维，并假定每根纤维处于单轴应力状态，在平截面假定的基础上依据纤维材料的本构关系确定截面的非线性行为，这在模拟轴力（特别是变轴力）与双向弯矩耦合作用下的非线性行为非常有效。
图1 结构受力形态与四个刚度法单元模拟效果
图2 常规框架集中铰设置
柔度法则是以单元截面力插值为出发点，形函数可以由力平衡条件直接确定，单元内任意截面的内力都可以由形函数得到，利用虚力原理，得到单元的柔度矩阵，求逆后得到单元刚度矩阵。柔度法虽然计算量有所增加，但是一个单元即可较精确地模拟构件的非线性行为，因而总体上计算效率仍然很高。
纤维铰从积分方法上还可以分为集中铰和分布铰两种形式，由于单元刚度矩阵需要截面刚度沿杆件长度进行积分得到，因此需要进行数值积分，积分点位置选取不同，计算结果也有差异。当选取塑性铰区域中某一有代表性的位置（通常取为塑性铰区的中间位置）的截面，用此截面刚度代表整个塑性铰区域的截面刚度，进而获得单元的刚度矩阵的方法叫做特征值截面法，也就是集中铰法。分布铰法是在构件中，依据积分法则选取若干截面，计算截面刚度后，再依据积分法则积分得到整个构件的单元刚度。
2. 设置纤维铰
2.1 纤维铰的定义
用户可以点击命令【定义】>【截面属性】>【铰属性】>【添加铰属性】，铰属性类型选择为FiberP-M2-M3，程序将弹出铰属性数据对话框，如图3所示。对话框中纤维定义选项即纤维的划分，当用户选择纤维定义来自于框架截面的默认值，程序会对钢截面、混凝土柱截面以及SD截面自动进行纤维的划分。图4为程序默认的混凝土柱的纤维划分，默认的划分原则是将纤维分为保护层、核心区以及钢筋三个部分。其中保护层采用的是默认的混凝土本构，可以认为是素混凝土本构；核心区的混凝土则可能不同，在定义混凝土本构时，如果参数化定义中选择了Mander，那么核心区混凝土将会自动变更为约束混凝土，否则将采用默认的混凝土本构；而钢筋纤维的位置和面积则来源于截面定义中对钢筋的布置。
图3 纤维铰的铰属性数据对话框
图4 柱铰纤维自动划分
当用户在纤维定义时选择用户自定义时，则可以处理更加复杂的情况，或者更加灵活地定义纤维及其本构曲线。如图5中用户自定义的梁纤维铰，与柱类似，分别考虑了保护层、核心区以及钢筋的分布，但是由于梁主要是绕3轴的弯曲行为，因此将纤维简化为一列，这样做减少了纤维数量会大大提高分析的效率而不会造成精度损失。
铰属性对话框中另一项为铰长度。铰长度的输入有两种方式，一种为相对长度（默认选项），一种为绝对长度。相对长度表示的是纤维铰区域占整个构件长度的比例，取值在0到1之间。而绝对长度则为纤维铰区域的绝对长度，默认单位为mm，如果认为铰长度与柱截面高度有关时，可以选择这种定义方式。
铰长度的取值应该慎重，铰长度不但对分析结果有影响，对可接受准则的影响更大。对于常规的梁柱，工程中常简化取铰长度为构件长度的0.1倍，并将铰布置于杆件的两端，这通常是符合实际破坏状态的。
但是对于一些特殊情况则不然，例如端柱，因为端柱本身应作为剪力墙的一部分，并不会像普通柱那样，有反弯点，两端发生弯曲破坏，而是作为墙的边缘构件受拉或受压，破坏状态也是整个构件受拉或受压破坏，此时铰长度取值应该与墙的铰长度相同，甚至直接取为1。
纤维铰中各个纤维的行为是通过材料定义的。用户可以通过【定义】>【材料】>【高级属性】>【非线性属性】进行定义，材料非线性属性定义包括三个部分：应力-应变曲线、滞回曲线和可接受准则。
应力-应变曲线用户可以通过两种方式进行定义，一种为参数化定义，一种为用户自定义。参数化定义是程序内置了相关材料的骨架曲线计算公式，这些公式可以通过【帮助】>【文档】>【Technical Notes】>【Material Stress-Strain Curves】进行查询，用户只需指定一些控制点即可，如图6红框部分所示，用户只需指定混凝土的峰值强度对应应变、极限应变以及极限应变后的下降段斜率即可，其中峰值强度取值为设计强度参数中的标准值（钢筋材料使用的是期望的屈服强度）；下降段的斜率为负值，其绝对值越大，下降越陡峭。
材料的滞回曲线共有8种（如图7），分别是弹性、随动、各向同性、退化、BRB硬化、Takeda、Pivot和素混凝土，可通过图6中紫色框的滞回类型进行设置，其中对于混凝土材料建议使用Takeda或素混凝土模型，而对于钢材或钢筋则宜使用随动或退化模型，使用不同的滞回曲线，纤维铰的行为也不尽相同。
图7 ETABS中滞回曲线
可接受准则用于判断构件的破坏程度，对于纤维铰的可接受准则是使用应变，先确定同类材料纤维的最大应变，再依据材料的可接受准则进行判断，可通过图6中蓝色框区域对可接受准则进行设置。
2.3 纤维铰的布置
用户选择需要指定纤维铰的框架对象后，可以通过【指定】>【框架】>【铰】，将纤维铰指定给框架对象，通常将铰指定于铰长度区域的中心。如图8示例，该杆件两端纤维铰的中心位于相对距离为0.05和0.95的位置。
纤维铰的布置也不必拘泥于放置于在构件两端，也不一定只放置一个纤维铰。如果需要布置多个纤维铰，用户可通过插入多个铰来近似模拟沿单元长度分布的塑性。例如，用户在一个杆件单元上布置10个铰，相对位置为0.05，0.15，0.25，...，0.95，假定每个铰的长度为杆件长度的 1/10 ，则可以模拟整根杆件任意位置的塑性行为。显然，添加更多的铰会增加计算量，因此只有在需要的情况下才这样做。
3.铰的选择
图9 ASCE41-17钢筋混凝土梁的塑性铰属性考虑了剪力影响
图11 ETABS模拟结果
图12 默认铰P铰定义
同样，塑性铰也存在一些缺陷，例如PMM铰只能采用各向同性滞回，这大大高估了柱的耗能能力，并且PMM铰在处理变轴力下的非线性行为时，很难提供足够的精度；塑性铰的骨架曲线也难以考虑加载制度对骨架曲线的影响。因此用户应该有足够的经验来判断选择使用哪种铰进行模拟。
如果用户做pushover分析：
由于不需考虑构件的滞回行为，更加推荐使用塑性铰，因为塑性铰考虑的因素更多一些，并且性能状态的判断更加直观，也有规范的支持，计算效率也更高一些。
如果用户做非线性时程分析：
对于梁、支撑等以弯曲或轴力破坏为主的构件，可以使用塑性铰（M铰、V铰、P铰等）；对于框架柱等以PMM行为为主的构件，更加推荐使用纤维铰，因为纤维铰对柱的压弯行为以及耗能能力描述的更加准确，这对时程分析而言是很重要的。
本文介绍了纤维铰的分类，按有限元实现的方法分为刚度法和柔度法，按积分方法分为集中铰法和分布铰法，其中ETABS分别采用的是刚度法和集中铰法。用户需要依据构件可能发生的破坏形式选择合理的铰长度以及铰位置，做到有的放矢。ETABS中有不同的铰类型可供选择，用户应结合构件的行为以及分析类型合理地选择不同类型的铰。
