不等待
即关注
简介
在进行多点输入时程地震反应分析时,为了确定结构的起振时间,必须确定地震波的传播方向,一般多点输入的影响是复杂的,应沿不同的传播方向,根据地震波的传播速度确定各点的相位差,分别计算结构的反应。
同时应该明确的是,地震地面运动的方向与地震传播的方向是无关的,不能将地震波传播方向与地面运动方向混淆,在同一地震传播方向上,需要考虑不同的地面运动方向的影响。
那么分析多点激励应该怎么做较为合理呢?
Step1. 确定各点地震输入时程
首先应根据工程场地情况确定各点的地震输入时程,通常有以下两种方法:
场地波动分析方法。考虑土与结构的相互作用,需要根据结构所处场地,建立包括结构和土体的分析模型,在分析模型的边界输入地震波,进行场地波动分析,得到支承各点的地震输入时程。
相位差简化输入方法。当场地比较均匀,只需要考虑行波效应时,可采用该方法进行分析。利用场地波或以往记录的地震波在结构的基底进行输入,各点输入地震波的形状和峰值大小相同,考虑行波效应后,各点输入的地震波存在一个相位差,这个相位差需要由地震波的传播方向和传播速度共同确定。通常地震波的传播速度范围在100~1000m/s。
首先根据地震波的传播方向,确定结构的第一个起振点,以此点作为地震波时程的零相位点;假设地震波是沿着传播方向平行传播的,同一波面上的地震波相位相同,根据地震波的传播方向和波速,计算地震波传播到基底各点的时间,确定相应的相位差。
Step2.进行多点输入分析
建立结构的有限元模型,在基底各点输入前一步得到的地震时程,进行结构的整体分析。根据不同的地震波传播方向,分别确定各点的地震输入时程。同时,因为地震波的振动方向与地震波的传播方向是相互独立的,针对每一个地震波的传播方向,必要时还要考虑地面振动的不同方向,分别进行时程分析。
小结讨论:
对于桥梁结构,因为下部基础相对独立,基础振动不同步引发的拟静力反应的影响比较明显;而普通建筑结构因为地下室连接在一起,各部分基础的相对变形受到地下室底板的约束,一般不会产生基础的相对变形,拟静力反应的影响比较小。但多点地震输入的另一个影响是因为各点输入的相位不同,相当于输入了一个扭转分量,会产生结构整体的扭转反应,这一点对建筑结构的多点地震分析是主要因素。在考虑了多点地震输入后,对称的多塔楼结构就有可能产生扭转反应,这在单点地震输入分析中是无法实现的。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6c8346a8ea0ca145c2161fb258750127.gif)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6c8346a8ea0ca145c2161fb258750127.gif)
例题演示
建立一个四跨的梁式结构,地震传播到支座之间时间差为5s(此处为了体现多点激励与同步激励差别更明显,选择特别软弱的地基,剪切波速选择非常小),计算分析该结构的动力特性及动力结果。
(注:此结构仅作计算演示,无工程意义)
为了方便研究利用自主开发的工具将加速度地震输入转化为位移时程输入。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-5ec50eabd19a20d031502f5a9c80d084.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-ae99f60ac62bf340010d9370858e8d7d.png)
左图为平常所见的加速度地震波,转换为右图的位移加载时程。为了验证加速度地震波加载和位移时程加载的计算结果是否一致,我们先进行一致激励分析,验证结果。
Step1:计算结构的前三阶周期和前三阶振型
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-24b45a8ec2285690e73abc88785a5c27.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-e4ae106a701a4f60b8f5e8ce48b99bb6.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-013297474692ef092ab3b84c55512beb.png)
Step2:将加速度时程和位移时程的地震波分别加入结构进行计算分析(均采用直接积分法进行分析),取结构中点(如下图所示)的绝对加速度和绝对位移时程结果进行对比。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-a3a2bc9a1c23dc8a77433a302f9e102a.gif)
加速度加载结果时程显示
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-8fd010b35eb1557d69c3a45b07a7c836.gif)
位移加载结果时程显示
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-d542019ba5221b6f17b511b015d309c8.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-bb47b216e28dd63841b7980c8a82a200.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-62c60d193e50b4bb1bb62d4f7fda1919.png)
可见地震波通过加速度加载和位移加载结果基本一致,验证了做法有效性和可实施性,内力对比如下图所示,也呈现基本一致,说明加载方法可行。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-0edc9d7e85a5c973d0724a07ed0b17aa.png)
Step3:考虑行波效应的多点激励施加,施加方式如下图所示,且依旧提结构中点的绝对位移时程进行对比。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-37697f5c557810094cd80a5462b4d81b.gif)
多点激励加载结果时程显示
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-f14e73fe3f318d372a764840a3621270.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-52e0695ce15e3e80e9ab4f9bf0317c13.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-9e085e668368d7a12396d3d121aebcd6.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6a0e556a1f1ebc8c338a93d80b29dd49.png)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-84fc6b2f11256ca0a834a16f61974a6b.png)
可以发现,通过多点激励之后,结构的内力在各个支座不协同运动下增大了许多,而由于结构支座不协同运动,导致中点的位移反而偏小。
接下来将原理阐述下,感兴趣的看官可以继续一起探讨。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6c8346a8ea0ca145c2161fb258750127.gif)
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6c8346a8ea0ca145c2161fb258750127.gif)
理论补充
首先谈谈支座静位移时的动力方程,将结构自由度分为两类,即为内部自由度和支座自由度,于是平衡方程的分块矩阵为:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-76725a97125d5bef198c520fdf86ff11.png)
由于我们采用的是位移加载模式,因此一般支座位位移量,外荷载为已知,而结构内部自由度位移和反力为未知量为未知量,将上式展开可得:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-0e1127a082d6d6aafd6e9692edbd9ced.png)
上面第一个式子又可写成:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-9408540601bc5c90b0eb9ecd50b36f80.png)
若无外荷载激励而仅产生支座位移,可变形为:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-fd1ec846dab10cb6183718a7b5c26f09.png)
这个就是支撑节点发生位移时在非支撑节点产生的静位移。
考虑多点激励时,动力方程为:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-cd6ca99be0fd861339e8ba7723646859.png)
在结构承受随时间变化的支座移动时,结构的反应来源由两部分:支座移动引起的结构反应,称为拟静力反应;支座移动加速度导致的惯性力引起的结构反应,称为动力反应,因此将位移表示为:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-8722ff1d2a82a830e9034f5bfa2b001f.png)
由前半部分得到的可知:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-c25508f3ab0db9abf06cfc2377848774.png)
且求一二阶倒数:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-4b84ac743a808e3b04e10536f9b135e3.png)
将上式子带入动力方程可得:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-11f6e0d1609e0fe08687e5f443a1cd9c.png)
若阻尼力相对惯性力而言可忽略不计,则式子可变形为:
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-07eed3712ba8991b74a787b2b1c253e8.png)
上面的推导就是多点激励的动力平衡方程。由推导可以看出拟静力反应和动力反应均与结构刚度有关,对于线性结构,拟静力反应和动力反应可以分别单独求解。
![](http://jycmf.cn/wp-content/uploads/2021/04/wxsync-2021-04-6c8346a8ea0ca145c2161fb258750127.gif)
概念为先,机理为本,下期见!
建源之光——工程侠
往期推荐 ·
往期精彩
#性能分析
#其他
~路过别擦肩~
欢迎关注微信公众号: 建源学堂