使用塑性铰做桥梁的动力弹塑性分析.docx
《使用塑性铰做桥梁的动力弹塑性分析.docx》由会员分享，可在线阅读，更多相关《使用塑性铰做桥梁的动力弹塑性分析.docx（16页珍藏版）》请在冰豆网上搜索。
使用塑性铰做桥梁的动力弹塑性分析
使用塑性铰做桥梁的动力弹塑性分析
北京迈达斯技术有限公司
2008年7月
目录
1．概要3
2．midasCivil中的塑性铰4
3．桥梁资料6
4．输入质量6
5．修改边界条件8
6．结构的非线性特性9
7．定义时程分析数据11
8．运行结构分析11
9．定义分析结果函数12
10．查看分析结果13
1．概要
结构抗震设计根据设防的目标不同有两种不同形式：
一种是弹性设计法，另一种是弹塑性设计法。
弹性设计法主要适合在较小地震作用下的结构抗震设计，它是以结构在设计地震作用下截面的应力保持在线弹性范围内为目标，用结构的弹性强度抵抗地震荷载。
与弹性设计法不同，弹塑性设计法是允许截面应力在地震时进入塑性范围的抗震设计方法，主要是通过提高结构极限变形能力的途径改善它的抗震性能，而不是简单地增加截面尺寸、提高截面强度来加强结构的抗震能力。
《公路桥梁抗震设计规范》（报批稿）6.3.6条，根据抗震设防的原则，E2地震作用下，允许结构出现塑性，发生损伤；即在E2地震作用下，桥梁已经进入非线性工作范围，因此只有进行结构非线性时程地震反应分析才能比较真实地模拟结构实际反应。
梁柱单元的弹塑性可以采用Bresler建议的屈服面来表示，也可采用非线性梁柱纤维单元模拟。
《公路桥梁抗震设计规范》（报批稿）7.4.1条，E2地震作用下，一般情况下，应按式7.4.2验算潜在塑性铰区域沿纵桥向和横桥向的塑性转动能力，但是对于规则性桥梁，可按式7.4.6验算桥墩墩顶位移，对于矮墩（高宽比小于2.5）的桥墩，可不验算桥墩的变形，但应按7.3.2条验算强度。
（7.4.2）
式中，
：
在E2地震作用下，潜在塑性铰区域的塑性转角；
：
塑性铰区域的最大容许转角。
（7.4.6）
式中，
：
在E2地震作用下墩顶的位移反应；
：
桥墩容许位移。
《铁路工程抗震设计规范》GB50111-2006中的7.3.3条，钢筋混凝土桥墩在罕遇地震作用下的弹塑性变形分析，宜采用非线性时程反应分析法，延性验算应满足下式的要求：
（7.3.3）
式中，
：
非线性位移延性比；
：
允许位移延性比，取值为4.8；
：
桥墩的非线性响应最大位移；
桥墩的屈服位移。
2．midasCivil中的塑性铰
屈服强度（面）的计算方法有用户输入和自动计算，在没有相关试验结果的情况下，一般采用自动计算。
按类型区分有“梁柱”、“弹簧”及“桁架”，在梁柱中又分为“集中铰”和“分布铰”。
“集中铰”通过转动和平移弹簧把结构的非弹性性能集中在单元的两端和中心，结构的其它位置假定为弹性，集中非弹性铰通过力矩与转角或者力与位移之间的关系定义，输出的时程分析结果非弹性铰的变形RX/RY/RZ都是转角。
“分布铰”假定整个构件均为非弹性，分布非弹性铰通过在积分点处力矩与转角或者力与位移之间的关系定义，输出的时程分析结果非弹性铰的变形RX/RY/RZ都是曲率。
作用类型有“无”、“强度P-M”及“状态P-M-M”。
“无”是不考虑轴力与弯矩的相互作用。
“强度P-M”即PM铰，考虑轴力对铰的弯曲屈服强度的影响，但对于两个方向弯矩间的相互作用是不考虑的。
PM铰只能考虑初始轴力P（初始重力荷载或用户输入的初始轴力）。
“状态P-M-M”即PMM铰，反映轴力和两个方向上弯矩的相互作用，P值是可变的，即可以考虑变化的轴力对屈服面的影响。
所谓变化的轴力是指在地震作用下产生的附加轴力引起的轴力的变化。
选择PMM铰后，一般只要定义FX方向的铰特性值，但这时只能选择kinematichardening滞回模型（随动强化模型），My、Mz方向的铰和FX方向的是相关联的，所以不用定义。
当轴力的变化对屈服的影响较大时（例如跨度较大的情况），或者需要考虑M-M之间的影响时（例如弯桥的情况），可以采用PMM铰；当轴力的变化不大且对屈服的影响比较小时，如跨度较小且桥墩的长细比不是很大时（不需要考虑P-delta效应），采用PM铰更好一些。
铰特性值，包括成分、铰位置（铰数量）、滞回模型等。
成分一般输入一个轴向和两个弯距成分的数据。
铰位置，选择集中类型时被激活，一般轴力成分选择单元中间，弯矩成分选择I端、j端或者两端。
铰数量，选择分布类型时被激活，输入积分点的数量，最多可输入到20个，根据输入的个数计算各个截面的力－位移或变形的关系。
kinematichardening滞回模型（随动强化型），初期加载时的效应点是在3条骨架曲线上移动的，卸载刚度与弹性刚度相同，随着荷载的增加强度有增加的趋势，适用于金属材料，它考虑了金属材料的包辛格效应，对于混凝土材料这样会过高的评价其耗能能力。
对于钢筋混凝土材料，考虑刚度退化影响的双直线模型是一种比较常用的模型，其中的Takeda（武田）模型和Clough（克拉夫）模型是两个比较典型的双直线弹塑性模型，但两者在考虑刚度退化时的变形规则是不一样的。
钢筋混凝土受弯构件的破坏过程常用混凝土开裂点、钢筋初始屈服点和极限强度三个特殊点来描述，如将原点、开裂点、屈服点、破坏点之间的力和变形关系用直线连接起来，就得到三直线模型骨架曲线，像Takeda（武田）模型也有对应的三直线模型和四直线模型。
特性值中的“延性系数”，是选择计算延性的基础。
当选择D/D1时延性系数是当前变形除以第一屈服变形，当选择D/D2时延性系数是当前变形除以第二屈服变形。
程序中采用集中铰时，D/D1是桥墩的非线性响应最大转角除以桥墩的第一屈服转角，不是直接和《公路桥梁抗震设计规范》（报批稿）中式（7.4.2）对应的，但是可以直接用于按《铁路工程抗震设计规范》GB50111-2006中的7.3.3条进行验算。
特性值中的“铰状态”，输入参考延性系数，即将铰的状态划分为5个不同的状态。
对于非对称铰程序会取各时间步骤中铰在正（+）、负（-）两个方向上较大的状态值来确定其状态。
Lever-1（0.5）表示铰还处于弹性阶段，Lever-2
（1）表示铰已达到屈服状态，Lever-3
（2）、Lever-4（4）、Lever-5（8）表示各构件不同的延性。
在分析结果中对于上述5个状态分别以蓝色、深绿色、浅绿色、黄色和红色来表示。
3．桥梁资料
该桥梁是高速公路上的一座大桥，属于新规范中的B类桥梁，所处区域地震烈度是7度，因此需要按8度进行抗震设防，按照新规范报批稿7.3与7.4节需要验算大震E2作用下桥墩的强度、潜在塑性铰区域转角、墩顶位移等。
桥墩高18m，截面是2m×2m的矩形截面，采用C50混凝土，桩基础，上面有盖梁，静力分析模型如图1所示。
图1.静力分析模型
4．输入质量
在midasCivil中输入质量有两种类型。
一个是将所建结构模型的自重转换为质量，还有一个是将输入的其它恒荷载（铺装及护栏荷载等）转换为质量。
对于结构的自重不需另行输入，即可在模型>结构类型对话框中完成转换。
而二期恒载一般是以外部荷载（梁单元荷载、楼面荷载、压力荷载、节点荷载等）的形式输入的，可使用模型>质量>荷载转换为质量功能来转换。
本例题也使用上述两种方法来输入质量。
首先将所输入的二期恒载（梁单元荷载）转换为质量。
模型/质量/将荷载转换成质量
质量方向>X,Y,Z
转换的荷载种类>梁单元荷载（开）
重力加速度（9.806）;荷载工况>上部结构自重
组合值系数
（1）;添加
图2.将梁单元荷载转换为质量
下面将单元的自重转换为质量。
模型/结构类型
将结构的自重转换为质量>按集中质量法转换
转换到X,Y,Z
图3.将结构的自重自动转换为质量
5．修改边界条件
建立桥梁抗震分析模型时，一般对扩大基础、沉井基础、锚碇等的处理比较简单，可视为固端。
而桩基础的处理比较复杂，常用的有两种方法：
（1）在墩底作用六个方向的弹簧等代群桩的作用，这六个弹簧刚度是竖向刚度、顺桥向和横桥向的抗推刚度、绕竖轴的抗扭刚度和绕两个水平轴的抗弯刚度。
它们的计算方法与静力法相同，只是考虑到在瞬间荷载作用下的土抗力应比持续荷载作用的大，一般取
。
（2）用三维梁单元模拟实际的桩基础，用土弹簧单元模拟桩周围土抗力的影响，地震波从桩端或者土弹簧处输入。
这里我们采用第一种方法来模拟桩基础，根据相关资料得到六个方向的刚度如表1所示。
表1桩基弹簧刚度
K1
竖向
K2
纵向
K3
横向
K4
扭转
K5
绕纵轴转
K6
绕横轴转
桩基弹簧
5.2e6
7.4e5
6.6e5
1.3e7
3.1e8
1.4e8
在midasCivil中，将静力分析模型中的墩底固定约束修改为弹簧约束，用模型>边界条件>节点弹性支承来模拟。
模型/边界条件/节点弹性支承
类型>线性
SDx（5.2e6）;SDy（7.4e5）;SDz（6.6e5）;
SRx（1.3e7）;SDx（3.1e8）;SDx（1.4e8）;
适用
图4.节点弹性支承
6．结构的非线性特性
柱截面设计
为了输入塑性铰的非线性特性，首先要在设计>钢筋混凝土构件设计参数>柱截面验算数据中输入柱截面钢筋。
然后在设计>钢筋混凝土构件设计参数>编辑钢筋混凝土材料特性中编辑混凝土、主筋以及箍筋的特性值。
当用户要确认塑性铰的延性系数时，要定义构件的屈服强度。
配筋截面的屈服强度，可使用程序中的设计功能自动计算。
如果是简单的截面可不经过本步骤，由用户直接输入也可。
图5.柱截面验算数据
图6.钢筋混凝土材料特性
定义非弹性铰特性值
下面定义非弹性铰特性值。
模型/材料和截面特性/非弹性铰特性值>添加
名称（PM）
屈服强度（面）计算方法>自动计算
类型>集中
定义>骨架
作用类型>强度P-M
材料>类型:
钢筋混凝土,规范:
ACI,名称:
1:
C50
构件>柱
截面>名称:
1:
Column
特性值>
Fx>铰的位置>Center
滞回模型>Clough
特性值:
自动计算（默认）
My>铰的位置>I
滞回模型>Clough
特性值:
自动计算（默认）
Mz的特性输入同My
屈服面特性值:
自动计算
确认
Fx
My
Mz
图7.定义非弹性铰特性值
分配非弹性铰
下面将前面定义的铰特性值分配给单元。
模型/材料和截面特性/分配非弹性铰
选项:
添加/替换
单元类型:
梁
非弹性铰特性值:
PM
选择单元:
1to9确认
图8.分配非弹性铰
7．定义时程分析数据
时程分析具体步骤可以参考技术资料“使用纤维模型做桥梁的动力弹塑性分析”，本文省略。
这里仅对需要注意的几点特别说明一下。
定义的塑性铰作用类型是PM铰，也就是考虑轴力对结构几何刚度的影响，所以将桥墩自重和上部结构自重都定义为时变静力荷载，在此基础上做地震动作用下的非线性时程分析。
该区域地震烈度是7度，E2地震作用（大震、罕遇地震）对应的加速度峰值是220cm/s2，因此在定义地震波时程函数时，对峰值加速度是341.7cm/s2要乘上0.644的折减系数。
这里在定义地震作用的非线性时程荷载工况时，阻尼的计算方法采用瑞利阻尼，取结构的第一阶和第三阶振型，它们的周期分别是2.25s和0.63s，阻尼比都取0.05。
8．运行结构分析
输入完上述数据后，运行结构分析。
9．定义分析结果函数
查看时程分析内力、位移等结果的方法参考“桥梁抗震分析与设计例题”，这里主要查看非弹性铰的地震响应。
为了以图形方式输出时程分析结果，定义结果输出函数。
荷载/时程分析数据/时程结果函数
定义函数>图形函数:
非弹性铰变形/内力>添加新的函数
名称:
1-My
单元类型>梁；单元号>1
结果类型>内力；成分>My；位置>I-Pos
时程荷载工况>Quake
确认↵
荷载/时程分析数据/时程结果函数
定义函数>图形函数:
非弹性铰变形/内力>添加新的函数
名称:
1-Ry
单元类型>梁；单元号>1
结果类型>变形；成分>Ry；位置>I-Pos
时程荷载工况>Quake
确认↵
图9.定义时程结果函数
10．查看分析结果
下面查看定义了塑性铰的柱的时程分析结果。
查看墩底处潜在塑性铰区域弯矩转角关系曲线。
结果/时程分析结果/时程分析图形
函数列表>1-My（开）
竖轴>从列表中添加;
横轴>1-Ry
图形标题（墩底处截面的弯矩转角关系）
X轴小数点位数（4）;Y轴小数点位数
（1）
类型>时程图表
图表↵
图10查看弯矩－转角关系
查看桥墩塑性铰的发展情况，这里通过查看“延性系数D/D1”和“铰状态”判断结构的状态。
结果/时程分析结果/非弹性铰状态
时程荷载工况名称>Quake
步骤>0.02
时间函数>Elcent_h
结果类型>延性系数
成分>Ry；数值类型>正
显示类型:
图例（开）;数值（开）;
小数点以下位数
（1）
动画（开）;开始时间（0.02）;结束时间（0.02）;增幅（0.2）
适用↵
图11查看塑性铰延性系数
除了动画的方式查看结果，还可以用表格的方式输出塑性铰的情况。
结果/分析结果表格/非弹性铰/延性系数（D/D1）
时程荷载工况>Quake（TH:
最大）
确定↵
图12查看塑性铰延性系数表格
升级会员 