采用ABAQUS数值模拟方法,建立了有力筋块体自由摇摆振动的数值模型,以块体的高宽比、力筋刚度、初始位移为分析参数,研究了不同参数水平下块体自由摇摆振动的恢复系数。研究发现恢复系数具有一定的随机离散性;高宽比对恢复系数有显著影响,力筋的施加以及初始位移的大小对恢复系数几乎没有影响;恢复系数的理论值,可以直接用于高瘦型块体的计算,但对矮胖型块体的计算则误差较大。给出了不同情况下恢复系数的试验平均值,可供计算时参考。 The values of coefficient of restitution (COR) of free rocking blocks with tendons were investigated by establishing numerical models using ABAQUS. The analysis parameters were the aspect ratio, tendon stiffness, and initial displacement. The results reveal that the COR showed some random nature. The aspect ratio had a significant effect on the COR, however the tendon stiffness and the initial displacement had little effect on the COR. The theoretical COR formula was accurate enough for slender blocks, but showed larger errors for dumpy blocks. The average values of COR for cases covered in this paper were given for reference use in related calculations.
梁艳苹,苏小卒*
同济大学土木工程学院建筑工程系,上海
收稿日期:2018年8月23日;录用日期:2018年9月10日;发布日期:2018年9月17日
采用ABAQUS数值模拟方法,建立了有力筋块体自由摇摆振动的数值模型,以块体的高宽比、力筋刚度、初始位移为分析参数,研究了不同参数水平下块体自由摇摆振动的恢复系数。研究发现恢复系数具有一定的随机离散性;高宽比对恢复系数有显著影响,力筋的施加以及初始位移的大小对恢复系数几乎没有影响;恢复系数的理论值,可以直接用于高瘦型块体的计算,但对矮胖型块体的计算则误差较大。给出了不同情况下恢复系数的试验平均值,可供计算时参考。
关键词 :块体,自由摇摆振动,恢复系数,力筋,ABAQUS
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
目前,除了采用消震与隔震措施,一般的建筑结构在地震作用下是通过自身的不可恢复变形和损伤来耗散地震能量,在大震作用下,即使结构没有发生倒塌,因耗散地震能量而产生的变形和损伤也常常会导致结构难以修复或加固费用昂贵,甚至需要推倒重建。而摇摆结构在地震作用下以刚体运动为主,变形和损伤均很小,因此受到越来越多研究者的关注 [
摇摆结构的概念由Housner [
Schau [
目前,虽然对摇摆振动已有一定的理论与试验研究,但理论研究基本都是基于块体和基础都为刚性且无弹跳与滑移的假定,同时相应的试验研究成果较少,如何建立合理有效的摇摆振动模型仍处于探索之中,各参数对振动响应的影响也有待进一步分析与讨论。
本文提出一种模拟摇摆振动的ABAQUS建模方法,通过与已有试验结果对比验证方法的正确性。模型根据实际情况定义材料特性和接触特性,避免了理论研究的刚体假定和无弹跳、滑移的假定,之后研究了不同参数水平下块体的碰撞恢复系数。
本文以Kalliontzis [
图1. 摇摆块体 [
图2. 有力筋的摇摆块体 [
图3. 试验中的摇摆系统(单位:m)
图4. 试验中摇摆系统的ABAQUS模型
对摇摆振动选用显式动力学方法建立接触条件更为容易,因此选用ABAQUS 6.13/Explicit作为分析平台。模型中各部分尺寸及材料参数按照试验取定,材料为弹性,力筋线膨胀系数根据钢结构设计规范取1.2 × 10−5。完成的ABAQUS模型如图4所示,对建模中关键内容说明如下。
在系统摇摆过程中,整体以刚体运动为主,只有柱底部与水泥层有碰撞,为减小模型的体量,将柱沿高度方向分段划分网格。上部和下部采用C3D8R单元,增强沙漏控制(enhanced),中间过渡部分选择C3D4单元。力筋选用T3D2,基础、质量块、水泥层以及锚板的单元类型与柱上部相同。
柱与质量块、基础与水泥层分别以一个部件建模;柱顶部和锚板、基础底部和锚板之间用束缚约束(Tie)连接在一起。以柱底面和水泥层顶面作为接触面定义通用接触,切向特性中摩擦公式选罚函数(Penalty),摩擦系数取0.6,法向特性中选硬接触(“Hard” Contact),其余取默认值。
试验通过五个分析步来实现。初始分析步,加z向重力g = −9.8 m/s2,约束基础底面三个方向的平动位移。第二个分析步,以降温法给力筋施加初始预应力,分析步时长0.01 s。第三个分析步,给柱体及质量块施加绕y轴初始角位移0.03 rad,时长0.2 s;为防止出现三维振动,约束柱体和质量块的y向平动位移直至分析结束。第四个分析步,使柱体及质量块在初始位移状态下保持稳定,时长0.05 s。第五个分析步,释放初始位移,使结构摇摆,时长5 s。
在ABAQUS/Explicit分析中,可以通过各种能量之间的比较,帮助评估一个分析是否得到了合理及准确的响应(ABAQUS 6.13帮助文件),通常总能量(total energy, ETOTAL)应保持为一个常数或接近于一个常数,伪应变能(artificial energy, ALLAE)与实际的能量如动能(kinetic energy, ALLKE)、内能(internal energy, ALLIE)相比应是可以忽略的。图5给出摇摆过程中的能量图,可以看到,总能量的数值基本保持在2~5之内,伪应变能接近于零,与动能和内能相比,总能量接近于一个常数,伪应变能可以忽略,因此认为模型结果计算有效。
通过GetData Graph Digitizer软件对试验结果(文献 [
由图6可以看到,在振动前4 s,试验结果与数值模拟所得结果在振动周期上吻合度很好,振幅间的差值随着振动时间的增加逐渐增加。造成振幅有差值的原因主要是数值模拟中没有考虑力筋与管道之间的摩擦产生的能量损失。实际试验中,力筋与管道之间不可避免的存在着摩擦,在结构摇摆的过程中,力筋反复伸长与回缩,与管道之间产生摩擦进而产生能量损耗,导致振幅减小。由于目前的建模方法还不能准确的模拟这部分耗能对振幅的影响,同时考虑到本文主要探讨不同参数对恢复系数的影响,因而本文数值试验中均不考虑此项摩擦。
图7为有力筋块体摇摆振动的计算简图。块体以角点O或O’为转动轴作刚体运动,当块体运动到 θ = 0 时,与基础发生碰撞,随之转动轴发生改变,因而摇摆振动的方程是分段连续的。根据块体绕角点O的动量平衡推导得到块体绕一个角点的运动方程如式(1):
I O d 2 θ d t 2 + W R sin ( α − θ ) + ( F 0 + F e ) b = 0 (1)
式中 I O = 4 MR 2 / 3 是矩形块体对角点O的动量矩,W是块体重力, F e = E s A s ε = E s A s tan α | sin θ / 2 | 是由力筋伸长产生的恢复力, F 0 = σ 0 A s 是由初始预应力产生的恢复力。由(1)式可以看到,在振动过程中,有力
图5. 系统能量图
图6. 角位移时程对比图
图7. 有力筋块体计算简图
筋的块体有三部分恢复力,W、 F e 及 F 0 ,当 F 0 = 0 和 F e = 0 同时满足时,(1)式即为无力筋的摇摆振动系统的运动方程,只有重力恢复力,同Housner [
理论方程中一般假定能量损耗集中在碰撞过程,并以Housner [
数值试验模型包括块体、基础、力筋及锚板(图8)。基础和块体的材料选C30混凝土,弹性模量为3 × 1010 Pa,密度为2400 kg/m3。块体与基础间的摩擦系数取0.6。所有模型中的基础尺寸均为1.5 m × 1.5 m × 0.4 m (长 × 宽 × 高)。力筋采用钢绞线,弹性模量为1.95 × 1011 Pa,密度为7850 kg/m3,线膨胀系数 α = 1.2 × 10 − 5 / ℃ ;锚板选用钢板,弹性模量2.06 × 1011 Pa,密度7850 kg/m3。为防止三维振动,约束块体y方向中面(即ABC面)的y向位
根据定义,恢复系数由碰撞后与碰移,分析步设置、边界条件以及荷载施加方式都与标定试验相同。
定义力筋刚度为 ∏ E = E s A s / W ,角位移为 Θ = θ / α ,数值试验中取块体宽度为0.4 m,初始预应力为0,以块体高宽比h/b、力筋刚度 ∏ E 及初始角位移 Θ 0 为变量设计了具有不同水平的12个数值模型(见表1),可以看到,这三个参数均是无量纲参数,当 ∏ E = 0 时即为无力筋的情况。模型名称按照“h/b- ∏ E - Θ 0 ”的方式编写,如4-300-0.1表示块体高宽比为4,力筋刚度 ∏ E = 300 ,初始角位移 Θ 0 =0 .1 。
撞前的动能之比求得,在本文计算中,由于在摇摆过程中块体无明显损伤,故认为块体质量在整个摇摆过程中保持不变,由碰撞后与碰撞前速度的平方之比求得各数值模型的恢复系数。图9(a)~图9(d)给出了12个数值试验的恢复系数随碰撞次数n的变化图,图中实线为恢复系数理论值(以下简称理论值),高宽比为4、6、10时,理论值分别为0.83132、0.92056、0.97052。
由图9可以看到,因为考虑了块体和基础的实际材质以及接触条件,数值试验所得的恢复系数(以下简称试验值)并不是一个固定值,而是在振动过程中以理论值为基准线随机波动,整体呈现离散状态,由文献 [
对块体在整个振动过程中的恢复系数取平均值,并将试验平均值r1与理论值r0进行比较,以e = (r1 − r0)/r0表示试验平均值与理论值的相对差值,见表2。
模型名称 | h/b | ∏ E | Θ 0 |
---|---|---|---|
4-0-0.1 | 4 | 0 | 0.1 |
4-300-0.1 | 4 | 300 | 0.1 |
4-650-0.1 | 4 | 650 | 0.1 |
6-0-0.1 | 6 | 0 | 0.1 |
6-300-0.1 | 6 | 300 | 0.1 |
6-650-0.1 | 6 | 650 | 0.1 |
10-0-0.1 | 10 | 0 | 0.1 |
10-650-0.1 | 10 | 650 | 0.1 |
10-6500-0.1 | 10 | 6500 | 0.1 |
6-0-0.5 | 6 | 0 | 0.5 |
6-70-0.5 | 6 | 70 | 0.5 |
6-300-0.5 | 6 | 300 | 0.5 |
表1. 模型试件参数
高宽比(初始位移) | 无力筋情况 | 有力筋情况1 | 有力筋情况2 | 理论解 | |||
---|---|---|---|---|---|---|---|
r1 | e (%) | r1 | e (%) | r1 | e (%) | r0 | |
4 (0.1) | 0.69588 | −16.292 | 0.76454 | −8.033 | 0.77014 | −7.360 | 0.83132 |
6 (0.1) | 0.95589 | 3.838 | 0.95761 | 4.025 | 0.96437 | 4.759 | 0.92056 |
10 (0.1) | 0.96489 | −0.580 | 0.9394 | −3.207 | 0.97978 | 0.954 | 0.97052 |
6 (0.5) | 0.94899 | 3.088 | 0.95207 | 3.423 | 0.94718 | 2.892 | 0.92056 |
表2. 恢复系数平均值
注:有力筋情况1包括4-300-0.1、6-300-0.1、10-650-0.1、6-70-0.5,有力筋情况2包括4-650-0.1、6-650-0.1、10-6500-0.1、6-300-0.5。
图8. 数值分析的ABAQUS模型
图9. 模型试件的恢复系数
从表2看,对高宽比为4的块体,力筋的施加对恢复系数r的平均值有较为明显的影响,但由图9(a)可以看到,因为r本身具有极大的离散性,所以实际的r并不因力筋的施加而具有显著变化;对高宽比为6和10的块体,无论是从表2中r的试验平均值方面还是图9(b)~图9(d)中r随碰撞次数变化方面,力筋的施加对恢复系数都不具有显著影响。由图9(b)、图9(d)以及表2可以看到,初始位移对恢复系数没有显著影响。
1) 在振动过程中,恢复系数r具有一定的随机离散性,并不像理论所得为一恒定值,且矮胖型块体的恢复系数离散性更大。
2) 对瘦高型的块体,其恢复系数可以用理论计算值表示,相对误差很小;但是对于矮胖型块体,理论值存在较大误差。本文给出的各工况下的恢复系数试验平均值,可供相关计算参考。
3) 块体是否加力筋以及所加力筋的大小,对恢复系数似无显著影响。
4) 初始位移的大小对恢复系数无明显影响。
本文的建模方法未考虑力筋与管道之间的摩擦,同时因为计算时间原因对高宽比为6和10的块体并未分析其运动全过程,只分析了前10 s的运动,之后的研究可以在这些方面进一步改进与探讨。
国家自然科学基金项目(51178328)。
梁艳苹,苏小卒. 有力筋块体自由摇摆振动恢复系数有限元分析 FEM Analysis of Restitution Coefficient of Free Rocking Block with Tendon[J]. 土木工程, 2018, 07(05): 725-733. https://doi.org/10.12677/HJCE.2018.75086