土木与环境工程学报  2019, Vol. 41 Issue (2): 99-105   PDF    
随机人群运动荷载作用下大跨度连廊的振动响应
张琼 1, 南娜娜 1, 朱前坤 1,2, 杜永峰 1     
1. 兰州理工大学 防震减灾研究所, 兰州 730050;
2. 大连理工大学 建设工程学部, 辽宁 大连 116024
摘要:基于社会力模型,研究了随机人群运动荷载作用下大跨连廊结构的振动舒适度问题。通过结合交通学领域、生物力学领域对人群运动特性及人体运动参数的统计分析,对随机人群运动进行仿真,得到不同密度工况下,任意行人在任意时刻的步速、落点位置等参数,建立随机人群荷载模型,计算不同工况密度下结构的加速度峰值。研究表明,加速度峰值与行人步速关系很大,随着行人密度增大,行人速度减小,结构加速度响应则先增大后减小。根据得到的结构加速度峰值与行人密度关系曲线,可以找出出现结构振动舒适度问题的行人密度区间。
关键词社会力模型    随机人群荷载    步速    峰值加速度    振动舒适度    
Vibration response of long-span corridor under stochastic crowd load
Zhang Qiong 1, Nan Nana 1, Zhu Qiankun 1,2, Du Yongfeng 1     
1. Institute of Earthquake Protection and Disaster Mitigation, Lanzhou University of Technology, Lanzhou 730050, P. R. China;
2. Faculty of Infrastructure Engineering; Dalian University of Technology, Dalian 116024, Liaoning, P. R. China
Abstract: Based on social force model which is combined statistics and analysis of crowd motion characters and body motion parameters in the field of transportation and biodynamic, vibration comfort of long-span corridor under stochastic crowd load is studied and stochastic crowd motion is simulated.The critical parameters such as pace and landing position are obtained in different densities, the stochastic crowd model is built and the peak acceleration is calculated.The study shows, peak acceleration is closely related to the pacing frequency, and pedestrian velocity becomes large along with increase of pedestrian density while the acceleration response of structure increased initially and then decreased. By the relationship curve of structure peak acceleration and pedestrian density the pedestrian density interval in which the structure suffering vibration serviceability problems can be found.
Keywords: social force model    stochastic crowd load    walking speed    peak acceleration    vibration serviceability    

近年来,随着新型轻质材料的发展,以及人们对结构新颖、美观的追求,使得轻质、低频结构的使用越来越广泛。人群荷载作为连廊、人行桥等结构的主要荷载,是引起结构振动的主要原因,当结构自振频率较低时,易与行人步频接近而发生共振,结构振动响应超过人体对振动舒适度的要求,可能引起行人的恐慌甚至导致结构的安全问题[1-3]

人群荷载模型是研究人致结构振动的基础,目前对单个行人对结构的响应研究较多,但对人群荷载的研究较少,研究表明,人群对结构的作用并不是单个行人的简单叠加[4]。行人荷载与行人步频有较大关系[5],人群平均步速与人群密度有关,行人在不同人群密度下受到其他行人的影响不同,每个行人的移动方式及移动特征的表现也不尽相同,鉴于人群移动的复杂性以及随机性,目前,人群对结构的作用大都采用简化的计算方法,主要包括固定多人步行力法[6]、移动多人步行力法、等效同步人数法[7]等,这些模型都在一定程度上提供了人群对结构动力响应的计算方法。以上假设均没有考虑行人个体的差异性以及个体间的随机性。丁国等[8]通过引入步频变异系数反映个体内的随机性,引入体重及均值步频的变化反映个体间的随机性,研究了行人荷载随机性对楼盖结构振动响应的影响,结果表明,荷载随机性对结构响应影响显著,在舒适度评价中不可忽略。杨娜等[9]根据生物力学以及交通领域内对人体质量、行人步频、步速等参数的统计分析,得到各参数的分布形式,再利用蒙特卡洛法产生相应的伪随机数,并结合人群密度的影响作用生成考虑人群集度的随机荷载模型。上述方法考虑了行人之间步态参数分布的随机性以及个体在时间域内的差异性,但是忽略了行人之间以及行人与环境的交互作用。因此,得到的行人落点、步速等并不符合行人的实际情况,计算结果与实际结构响应可能会有较大偏差。

由Helbing等[10]提出的社会力模型,将人的主观愿望、行人与行人之间、行人与环境之间的相互作用用社会力的概念表示,成功模拟了“快即是慢”“瓶颈摆动”等行人现象。笔者结合交通学以及生物力学领域内对行人步速、体重以及行人尺寸等参数的统计,考虑行人个体差异性以及个体之间的随机性,利用社会力模型模拟人群的运动过程,得到任意时刻所有的行人落点、步速等;建立基于社会力模型的人群荷载随机模型,并对不同行人密度下连廊结构动力响应进行计算,得到结构峰值加速度与行人密度的关系曲线,对随机人群荷载作用下结构振动舒适度进行评价。

1 单人步行荷载模型

人行走时两条腿交替移动,步行荷载对结构的作用具有周期性,其对连廊等结构产生的激励主要为竖向激励,时域内步行荷载广泛采用傅里叶级数的形式[11]

$ F(t)=G+G \sum\limits_{i=1}^{n} \alpha_{i} \sin \left(2 i \pi f_{\mathrm{p}} t+\varphi_{i}\right) $ (1)

式中:αii阶谐波的动力荷载因子,定义为αi=Ai/Gφi为第i阶谐波的相位角;G为行人重量;fp为行人步频;n为模型中考虑的阶数;相位角φi=0。考虑前3阶的影响,动载因子为

$ \begin{aligned} \alpha_{1}=&-0.2649 f_{\mathrm{p}}^{3}+1.3206 f_{\mathrm{p}}^{2}-\\ & 1.7597 f_{\mathrm{p}}+0.7613 \end{aligned} $ (2.1)
$ \alpha_{2}=0.07 $ (2.2)
$ \alpha_{3}=0.06 $ (2.3)
2 随机人群荷载模型的建立
2.1 基于社会力模型随机人群运动仿真

社会力模型认为行人的运动在社会力的作用下发生,其中包括自驱力、行人之间的作用力、行人与周边环境之间的作用力。笔者在此基础上引入了减速避让机制[12]以及自停止机制[13]。当行人的动态空间被其他行人占据时,当前行人会采取减速机制以避免与其他行人发生碰撞,而当周围行人对当前行人斥力过大,使行人的速度小于0,则赋与当前行人速度为0。自驱力

$ F_{\mathrm{d}}=m_{i} \frac{{\mathit{\boldsymbol{v}}}_{\mathrm{e}} {\mathit{\boldsymbol{e}}}_{i}-{\mathit{\boldsymbol{v}}}_{i}}{T_{\alpha}} $ (3)

式中:mi为行人质量;ve为行人的期望速度;ei为行人期望运动方向;vi为行人的运行速度;Tα为松弛时间。行人之间的相互作用力包括行人之间的心理排斥力、减速避让力和物理作用力。心理排斥力为

$ \begin{array}{c}{{\mathit{\boldsymbol{f}}}_{i j}^{\mathrm{soc}}=A \exp \left(\frac{r_{i j}-d_{i j}}{B}\right)} \cdot \\ {\left(\lambda_{i}+\left(1-\lambda_{i}\right) \frac{1+\cos \varphi_{i j}}{2}\right) \boldsymbol{n}_{i j}}\end{array} $ (4)

式中:A为行人之间的作用强度;B为行人之间的作用力范围;(rijdij)为行人之间的距离;λi为各向异性参数;φij为行人间的斥力与期望运动方向的夹角。减速避让力

$ \begin{array}{c}{{\mathit{\boldsymbol{f}}}_{i}^{\mathrm{avo}}=-\delta_{i}(t) {\mathit{\boldsymbol{v}}}_{i}(t) m_{i}} \cdot \\ {\left(\lambda_{i}+\left(1-\lambda_{i}\right) \frac{1+\cos \left(\varphi_{i k}\right)}{2}\right) \boldsymbol{n}_{i j}}\end{array} $ (5)

当(rijdij)<dα(t)时,行人会受到减速避让力。式中:di=|(a+b|vi(t)|)ei(tnib|为按原速运动后的行人之间的位置差;$ - {\delta _i}{\mathit{\boldsymbol{v}}_i}(t) = - \frac{{{\delta _i}(t)}}{{{T_{\rm{r}}}}}$为行人i由速度vi(t)减至0而采取的加速度;Tr为行人反应时间;φik为行人i实际步行速度与行人k作用于行人i的排斥力的反方向的夹角。

$ {\mathit{\boldsymbol{F}}_{ij}} = \left( {K\mathit{\Theta }\left( {{r_{ij}} - {d_{ij}}} \right)} \right){\mathit{\boldsymbol{n}}_{ij}} + k\mathit{\Theta }\left( {{r_{ij}} - {d_{ij}}} \right)\Delta {\mathit{\boldsymbol{v}}_{ij}}{\mathit{\boldsymbol{t}}_{ij}} $ (6)
$ \mathit{\Theta } = \left\{ {\begin{array}{*{20}{c}} {\left( {{r_{ij}} - {d_{ij}}} \right),}&{\left( {{r_{ij}} - {d_{ij}}} \right) \le 0}\\ {0,}&{\left( {{r_{ij}} - {d_{ij}}} \right) > 0} \end{array}} \right. $ (7)

式中:nij为行人j指向行人i的单位向量;tij由单位法向量nij逆时针旋转90°得到;Kk为常数;Δvij为两行人速度的矢量差。

行人与障碍物的相互作用与行人之间的相互作用原理一致,包括行人对障碍物的心理排斥力、减速避让力以及物理作用力。

$ \begin{array}{*{20}{c}} {F_{i{\rm{w}}}^{{\rm{soc}}} = \left( {{A_{\rm{w}}}\exp \left( {\frac{{{r_i} - {d_{i{\rm{w}}}}}}{{{B_{\rm{w}}}}}} \right) + K\mathit{\Theta }\left( {{r_i} - {d_{i{\rm{w}}}}} \right)} \right){\mathit{\boldsymbol{n}}_{i{\rm{w}}}} + }\\ {k\mathit{\Theta }\left( {{r_i} - {d_{i{\rm{w}}}}} \right)\left( {{v_i} \cdot {\mathit{\boldsymbol{t}}_{i{\rm{w}}}}} \right){\mathit{\boldsymbol{t}}_{i{\rm{w}}}}} \end{array} $ (8)

式中:Aw为行人与障碍物作用力强度;Bw为行人与障碍物作用范围;(ridiw)为行人半径与行人到障碍物的法向距离差;niw为由障碍物指向行人的单位法向量;tiw为平行于障碍物的单位切向力;[vi· tiw]为行人速度在障碍物方向上的投影。

利用MATLAB软件对社会力模型进行仿真,仿真环境为一个长22.0 m、宽3.0 m的连廊。假定行人质量在50~80 kg的范围内服从均值为65、方差为5的正态分布,行人期望速度与行人中老人、妇女的所占比重有关,假定在0.7~1.4 m/s的范围内服从均值为1.29、方差为0.29的正态分布,半径在0.2~0.3 m的范围内服从均匀分布。其他参数的取值参考文献[13],如表 1,仿真模拟如图 1所示。模拟结果表明,行人能够有效地对其他行人以及障碍物进行避让,将得到的密度-速度曲线与其他人实测的或模拟结果进行对比,结果吻合度较高。

表 1 模型参数 Table 1 Parameters of model

图 1 人群模拟 Fig. 1 Crowd simulated

2.2 行人速度-密度关系

利用社会力模型分别对2~90人范围内选取几种不同密度工况进行仿真,行人位置分布会影响每个行人的受力,进而影响行人速度以及结构的动力响应,考虑到行人的位置分布、期望速度等的随机性以及仿真过程的波动性,每种密度工况下产生7组数据,计算每种工况下速度的平均值,进而得到行人速度-密度曲线,见图 2。将均值速度作为行人在某一密度下的移动速度,并与其他学者的结果进行对比,证明该社会力模型的可靠性。将得到的点进行拟合得到了密度-速度关系式

图 2 速度-密度曲线 Fig. 2 Speed-density curve

$ v=1.28-0.33 k $ (9)

式中:k为行人密度,即单位面积的行人数。由图 2速度-密度曲线可知,在行人密度小于0.5时行人速度的变化并不稳定,原因是,当密度小于0.5时,行人处于自由运动的状态,速度具有较大的波动性,与实际情况符合;当密度大于0.5时,处于中密度状态,行人之间的相互作用开始变得明显,速度随着密度的增大呈稳定下降趋势。由文献[14],速度与行人步频的关系可表示为

$ f_{\mathrm{p}}=0.35 v_{\mathrm{p}}^{3}-1.59 v_{\mathrm{p}}^{2}+2.96 v_{\mathrm{p}} $ (10)
2.3 随机人群荷载

随机人群荷载是单个行人荷载的组合。社会力模型中考虑到个体内的随机性,质量在50~80 kg的范围内服从均值为65、方差为5的正态分布,期望速度在0.7~1.4 m/s的范围内服从均值为1.29、方差为0.29的正态分布。而行人真实速度则由社会力模型中任意时刻行人受力以及前一步速度共同决定。上述过程中,不仅在时间域内考虑了行人个体内以及个体间的随机性,并且在空间域内考虑了行人运动过程中周围环境以及与周围行人之间的相互作用。根据社会力模型,行人下一时刻的速度以及位置由当前行人受到的来自周围行人以及障碍物的力共同决定,不同密度下行人受到其他行人影响不同。计算不同密度下任意行人任意时刻的位置分布、速度等,由式(10)可得到任意时刻行人的步频,再由式(11)对随机人群中所有荷载进行叠加则可得到人群随机荷载,而人体自重引起的结构动力响应相较于简谐荷载而言可忽略不计[5]

$ F P(t)=\sum\limits_{p=1}^{m}\left(G+G \sum\limits_{i=1}^{n} \alpha_{i} \sin \left(2 i \pi f_{\mathrm{p}} t+\varphi_{i}\right)\right) $ (11)

式中:m为行人数;n为在计算中考虑的总阶数;αii阶谐波的动力荷载因子,定义为步行荷载傅里叶幅值谱峰值与人体体重之比;φi为第i阶谐波的相位角;G为行人体重;fp为行人步频。

图 3图 4FP为所有行人步行荷载之和,G为结构上所有行人的体重之和。图 3为单个行人荷载时程曲线图,从图 3可以明显看出人行荷载具有周期性。图 4对比了行人同步行走荷载与随机人群荷载时程曲线,其中人数为10人,同步行走中将10人排成3列,前2列均为4人,第3列为2人,行、列间距一致均为0.75 m。行走过程中行人步速均为1.3 m/s, 步频均为1.93 Hz。经过对比可知,同步人群荷载时程与单个行人荷载时程类似,具有很明显的周期性,而随机人群荷载由于步频并不统一,无明显周期,行人行走过程中出现同步情况较少,因此,荷载峰值也比同步行走时的小。

图 3 单人步行荷载曲线 Fig. 3 Single walking load curve

图 4 人群荷载时程曲线 Fig. 4 Crowd load time-history curve

3 连廊的振动控制方程

将连廊简化为简支梁,单位长度质量为m、刚度为k、阻尼为cI为惯性矩、E为弹性模量,连廊长度为L=21.8 m, 宽为3 m。

3.1 随机人群荷载下连廊的振动控制方程

在人群荷载的共同作用下,结构的振动控制方程

$ m \frac{\partial^{2} u}{\partial t^{2}}+c \frac{\partial u}{\partial t}+E I \frac{\partial^{4} u}{\partial x^{4}}=\sum\limits_{i=1}^{n} F_{i}(t) \delta\left(x_{i}-v_{i} t_{i}\right) $ (12)

式中:下标i表示为第i行人;δ(xvt)为Dirac函数,其性质为

$ \delta(x-v t)=0, \quad x \neq v t $ (13.1)
$ \int_{0}^{L} \delta(x-v t) \mathrm{d} x=1 $ (13.2)
$ \int_{0}^{L} f(x) \delta(x-w t) \mathrm{d} x=f(v t) $ (13.3)
3.2 控制方程的求解

结构的位移变形可表示为

$ u(x, t)=\sum\limits_{i=1}^{\infty} \varphi_{i}(x) z_{i}(t) $ (14.1)

式中:φi(x)为振型; zi(t)为广义坐标。结构的振动过程,以一阶模态为主,故忽略其他阶的影响,振型函数为$\varphi_{i}(x)=\sin \frac{n \pi x}{L}$n为振型的阶数,本文中只考虑一阶振型。

$ u(x, t)=\sum\limits_{i=1}^{\infty} \sin \frac{\pi x}{L} z_{i}(t) $ (14.2)

将式(13.2)带入结构振动方程,其中,等式右边广义荷载

$ \sum\limits_{i=1}^{n} F_{i}=\int_{0}^{L} \delta\left(x_{i}-v_{i} t\right) F_{\mathrm{p}_{i}} \sin \frac{\pi x_{i}}{L} \mathrm{d} x=\sum\limits_{i=1}^{n} F_{\mathrm{p}} \sin \left(\frac{\pi v_{i} t}{L}\right) $ (15.1)

等式左边广义质量

$ \widetilde{m}=\int_{0}^{L} m \sin ^{2} \frac{\pi x}{L} \mathrm{d} x=\frac{m L}{2} $ (15.2)

广义刚度数

$ \tilde{k}=\int_{0}^{L} E I\left(\frac{\pi^{2}}{L^{2}}\right) \sin ^{2} \frac{\pi x}{L} \mathrm{d} x=\frac{\pi^{4} E I}{2 L^{3}} $ (15.3)

则结构的振动控制方程可化为

$ \ddot{z}(t)+2 \zeta w_{n} \dot{z}(t)+w_{n}^{2} z(t)=\frac{m L}{2} \sum\limits_{i}^{n} {F_{{\rm{p}}i}} \sin \left(\frac{\pi v_{i} t}{L}\right) $ (16)

利用MATLAB平台,通过基于Newmar k-β法自编的程序对控制方程进行求解,计算流程如图 5所示。

图 5 结构动力响应计算流程图 Fig. 5 Flow chart of structural dynamic response calculation

4 算例

某商贸城的钢结构连廊跨度L=21.8 m,宽度为3.0 m,结构的阻尼比ξ为0.01。该钢结构连廊截面竖向等效抗弯刚度为EI=3.268×109 N·m2,单位长度质量m=1 063.5 kg,单位长度阻尼系数c为950.8 N/(m·s-1),连廊上基于社会力模型的人群随机荷载分布如图 6所示。《建筑楼盖结构振动舒适度技术规范(送审稿)》中规定,连廊的第一阶竖向自振频率不宜小于3 Hz,不封闭连廊竖向峰值加速度限制为0.5 m/s2。该结构竖向基频为4.72 Hz,满足频率限制的要求。当单个行人步频为结构基频的1/3通过该连廊,由步频和速度的关系式得行人的速度为0.838 m/s结构加速度峰值如图 7所示为0.12 m/s2,小于加速度限制。

图 6 行人随机分布图 Fig. 6 Pedestrian random distribution

图 7 单个行人加速度时程曲线 Fig. 7 Acceleration time history curve of single pedestrian

《建筑楼盖结构振动舒适度技术规范(送审稿)》中指出,当行人密度为0.1人/m2是属于完全自由行走状态。当为0.5人/m2时行人较多,行人仍然能够自由选择行走,当行人密度等于1.0人/m2时,行人行走有不舒适感,人群拥挤。分别对行人密度为0.167、0.333、0.5、1.0人/m2共4种行人密度工况下的人群荷载引起的结构响应进行计算。得到结果如图 8所示。从图 8可以看出当行人密度为0.167、0.333、0.5、1.0人/m2时结构加速度最大响应分别为0.066、0.143、0.255、0.535 m/s2,当行人密度为1.0人/m2时,结构加速度响应大于0.5 m/s2,超过连廊振动舒适度限制。由于不同行人以及期望速度等对行人的运动特性影响较大,而在同一种行人密度工况下,行人位置分布随机,期望速度取值服从正态分布。考虑到行人的随机性以及仿真过程的波动性,同一行人密度工况下引起的结构响应不同,针对10种行人密度工况各产生了8组不同行人模拟。表 2为在各行人情况下8组不同行人分布下计算得到的结构加速度峰值,其计算流程如图 5所示,取8组结构加速度峰值的平均值作为对应行人密度下加速度响应的代表值,用以评估结构的振动舒适度。为了得到更加明显的变化规律,将得到的峰值加速度与人群密度进行函数关系拟和,该拟合只针对该结构。拟合曲线如图 9R2=0.922拟合度较高,拟合结果为

图 8 不同行人密度下加速度时程曲线 Fig. 8 Acceleration time history curve in different pedestrian density

表 2 不同人群密度下结构加速度峰值 Table 2 Peak acceleration of structure in different pedestrian density

图 9 人群密度-加速度时程曲线 Fig. 9 Pedestrian density-acceleration time history curve

$ \ddot{y}(t)=-0.314 p^{2}+0.836 p-0.0282 $ (17)

为了验证曲线的可靠性,分别计算0.25、0.583、0.75、0.917人/m2几种行人密度工况下结构加速度峰值,与公式拟合的值进行对比,结果见表 3。通过对比实际峰值加速度与拟合公式计算得到的加速度峰值,发现结果较为吻合。由图 9密度-加速度峰值关系曲线可知,当密度在1.15~1.30人/m2之间时,速度在0.89~0.84 m/s范围内变化,相应的由拟合公式计算得到的速度平均值为0.841 m/s,与之对应的步频在1.57~1.63 Hz范围内变化,接近结构基频的1/3。故而结构在1.15~1.30人/m2的密度范围内取最大值。当行人密度在1.0~1.4人/m2时,结构加速度响应超出了舒适度限制,应该采取一定措施控制人群密度出现在这个范围的概率,或者采取减振措施。

表 3 实际与拟合的加速度峰值 Table 3 Peak acceleration of actual value and fitted

5 结论

1) 以社会力模型为基础,结合生物力学以及交通领域内对于行人特性的研究,考虑了行人时间和空间内的随机性以及行人之间、行人与环境之间的交互性,建立了更符合实际的随机人群荷载模型。

2) 人群随机荷载引起的结构振动与行人步速有较大关系,行人步速与行人密度有直接关系,可以通过行人密度与速度的关系曲线来预测不同密度下人群引起的结构加速度响应,找出人群荷载引起结构过量振动的行人密度范围,采取措施进行减振控制。

3) 人群随机荷载引起的结构动力响应,并不会一直呈现增大的趋势,当人群密度过大时,人群平均步频远小于结构基频的1/3,此时可将行人荷载看作是静载,引起的动力响应很微弱。

参考文献
[1]
樊健生, 李泉, 李全旺, 等. 考虑人-结构相互作用的楼盖振动控制研究[J]. 振动与冲击, 2010, 29(11): 230-236.
FAN J S, LI Q, LI Q W, et al. Floor vibration control considering human- structure interaction[J]. Journal of Vibration and Shock, 2010, 29(11): 230-236. DOI:10.3969/j.issn.1000-3835.2010.11.050
[2]
简方梁, 吴定俊, 李奇. 上海虹桥车站人行走廊人致振动分析[J]. 振动与冲击, 2010, 29(8): 136-140.
JIAN F L, WU D J, LI Q. Human-induced vibration analysis of pedestrian suspension corridor in Shanghai Hongqiao Station[J]. Journal of Vibration and Shock, 2010, 29(8): 136-140. DOI:10.3969/j.issn.1000-3835.2010.08.031
[3]
张家华, 高承勇, 秦卫红, 等. 大跨度楼面人群荷载作用响应分析[J]. 建筑结构学报, 2011, 32(11): 177-182.
ZHANG J H, GAO C Y, QIN W H, et al. Analysis of vibration responses of long-span floor under crowd loads[J]. Journal of Building Structures, 2011, 32(11): 177-182.
[4]
VENUTI F, RACIC V, CORBETTA A. Pedestrian-structure interaction in the vertical direction: coupled oscillator-force model for vibration serviceability assessment[C]//Ⅸ International Conference on Structural Dynamics Eurodyn, Porto, Portugal, 2014.
[5]
朱前坤, 李宏男, 杜永峰, 等. 不同行走步速下人行桥振动舒适度定量化评估[J]. 工程力学, 2016, 33(10): 97-104.
ZHU Q K, LI H N, DU Y F, et al. Quantitative evaluation of vibration serviceability of pedestrian bridge under different walking speed[J]. Engineering Mechanics, 2016, 33(10): 97-104.
[6]
甄伟, 盛平, 张力, 等. 财富中心悬吊楼盖结构舒适度设计及检测[J]. 建筑结构, 2011, 41(1): 52-54, 34.
ZHEN W, SHENG P, ZHANG L, et al. Vibration comfort design and experiment of the suspended floor of prosper center[J]. Building Structure, 2011, 41(1): 52-54, 34.
[7]
潘宁, 杜文博, 倪志伟, 等. 计算人群行走楼板响应的随机振动方法[J]. 建筑科学, 2016, 32(5): 39-43.
PAN N, DU W B, NI Z W, et al. Random vibration method for calculating floor vibration response under crowd walking[J]. Building Science, 2016, 32(5): 39-43.
[8]
丁国, 陈隽. 行人荷载随机性对楼盖振动响应的影响研究[J]. 振动工程学报, 2016, 29(1): 123-131.
DING G, CHEN J. Influences of walking load randomness on vibration responses of long-span floors[J]. Journal of Vibration Engineering, 2016, 29(1): 123-131.
[9]
杨娜, 邱雪. 人群集度荷载模型及人致振动模拟分析[J]. 北京交通大学学报, 2016, 40(3): 88-96.
YANG N, QIU X. Crowd density load model and analysis of crowd-induced vibration[J]. Journal of Beijing Jiaotong University, 2016, 40(3): 88-96. DOI:10.11860/j.issn.1673-0291.2016.03.015
[10]
HELBING D, BUZNA L, JOHANSSON A, et al. Self-organized pedestrian crowd dynamics:experiments, simulations, and design solutions[J]. Transportation Science, 2005, 39(1): 1-24.
[11]
PFEIL M, AMADOR N, PIMENTEL R, et al. Analytic-numerical model for walking person-footbridge structure interaction[C]//International Conference on Structural Dynamics, Eurodyn, 2014.
[12]
李珊珊, 钱大琳, 王九州. 考虑行人减速避让的改进社会力模型[J]. 吉林大学学报(工学版), 2012, 42(3): 623-628.
LI S S, QIAN D L, WANG J Z. Improved social force model considering pedestrian deceleration to avoid collision[J]. Journal of Jilin University (Engineering and Technology Edition), 2012, 42(3): 623-628.
[13]
PARISI D R, GILMAN M, MOLDOVAN H. A modification of the Social Force Model can reproduce experimental data of pedestrian flows in normal conditions[J]. Physica A Statistical Mechanics and Its Applications, 2009, 388(17): 3600-3608. DOI:10.1016/j.physa.2009.05.027
[14]
BRUNO L, VENUTI F. Crowd-structure interaction in footbridges:Modelling, application to a real case-study and sensitivity analyses[J]. Journal of Sound and Vibration, 2009, 323(1): 475-493.