地震地面运动具有强烈的随机性和不确定性.docx
- 文档编号:8306931
- 上传时间:2023-01-30
- 格式:DOCX
- 页数:11
- 大小:48.09KB
地震地面运动具有强烈的随机性和不确定性.docx
《地震地面运动具有强烈的随机性和不确定性.docx》由会员分享,可在线阅读,更多相关《地震地面运动具有强烈的随机性和不确定性.docx(11页珍藏版)》请在冰豆网上搜索。
地震地面运动具有强烈的随机性和不确定性
地震地面运动具有强烈的随机性和不确定性,地震作用下地面建筑物的动力响应同样具有强烈的随机性。
国外地震工程学者在40年代末就开始将随机振动理论引入到建筑结构抗震研究中来,国内有关研究工作始于60年代初期并在进入80年代后逐渐活跃。
但是应该看到,随机抗震理论的实际应用状况在国内外还不尽人意,运用随机抗震理论进行工程结构地震响应研究和分析的工作在国内外都开展得远远不够,与这种理论的先进性很不相称。
本文针对二维重力坝结构,详细论述结构随机地震动力分析方法,并对该方法实际运用过程中的一些问题提出作者的见解及具体处理方法,最后还给出一重力坝工程随机地震动力分析的实例。
1地震动随机输入模型
Housner在1947年引入白噪声地震动模型。
日本学者金井清在1957年提出了被称之为过滤白噪声模型的地震动随机输入模型,其形式较为简捷,物理概念明确,与大量地震动实测资料符合得较好。
Penzien和S.C.Liu曾经采用金井清模型产生50条人工地震波,以50条人工地震波样本得到的速度反应谱与Housner的设计反应谱曲线十分接近。
这一研究结果对金井清模型是有力的支持,不过金井清模型也存在一些缺陷,主要包括对地震动低频成分有所夸大,地震动位移、速度以及加速度导数过程的均方值因相应谱函数存在奇异点或广义积分不收敛而无法求得。
针对这些缺陷近二、三十年来国内外学者在金井清谱的基础上提出了不胜枚举的修正模型。
公正地说这些模型中不乏优者,但或许是对其合理性和可靠性的论证力度不够,至今尚未有哪个模型的影响程度真正能与金井清模型相匹敌。
下面的论述中我们不难看出金井清谱的上述缺陷对其实际运用并不形成大的障碍。
1.1KanaiTajimi地震动模型金井清地震动模型被认为是日本学者Kanai和Tajimi共同提出的,所以有时也称为KanaiTajimi地震动模型。
现有的强震地面运动加速度记录的傅里叶分析表明,即使是在频率的一个有限区段内傅里叶振幅谱也不是常数,这些振幅谱有点波动的性质,可能在一个或几个频率点上出现波峰,而波峰随频率的增大而衰减。
所有这些使我们可以假定只要适当选择过滤转换特性,一个平稳过滤白噪声随机过程模型更能代表实际地面运动。
Kanai和Tajimi建议的滤波函数为
|H(ω)|2=(ωg4+4ξg2ωg2ω2)/(ωg2-ω2)2+4ζg2ωg2ω2
(1)
式中ωg和ξg为地面特征频率和特征阻尼比。
根据这个滤波函数,地面加速度过程ag(t)的功率谱密度函数为
S(ω)=|H(ω)|3S0
(2)
式中S0为功率谱强度。
1.2谱强度与地震最大加速度关系目前使用最多的地震动统计参数是最大地震加速度am.我国国家地震局1980年修订的地震烈度表给出了对应各地震烈度的地面地震动最大加速度值。
然而在随机地震动力分析中需要用到的地震动参数是地面地震动功率谱强度S0.为更好地利用现已积累的地震动统计资料,建立地面地震动功率强度与地震最大加速度之间的关系十分重要。
将地面地震加速度过程视为随机过程,每一实测地震加速度记录即为随机过程的一次实现。
那么通常用到的地震最大加速度建议值可在统计意义上更严格地定义为地震动最大加速度期望值E[am]。
根据Davenport公式[2],有
E[am]=rpσa
(3)
式中rp称为地面地震动加速度峰值因子;σa为地面地震动加速度过程ag(t)均方根值。
式(3)建立了地震加速度均方值与最大地震加速度期望值之间的关系,通过该式可由地震最大加速度统计参数计算地震动功率谱强度S0.
地震加速度均方值
,式中Sa(ω)为地面地震动加速度过程的功率谱密度函数,采用金井清谱,将式
(2)和式
(1)代入上式,有
σa2=[(S0ωgπ)/4ξg](1+4ξg2)
(4)
峰值因子
rp=C+0.5572/C
(5)
其中
,T为强震持时,
为加速度导数过程的均方根值。
式(3)和式(4)即建立了谱强度与地震最大加速度之间的转换关系。
如果采用式(5)计算峰值因子,一定要计算地震加速度导数过程的均方根值
。
金井清地震动模型的
是没有收敛解的。
有一些金井清谱的修正模型正是以此为出发点提出对原模型的修正。
但实际上,式(3)是在统计意义上成立的一个关系式,通过大量实测地震记录的统计分析确定rp的取值,似为一更合理的途径。
在这个意义上来说,金井清谱的缺陷并不影响该模型的实际运用。
已有的研究工作表明rp的值是比较稳定的,受各参数的影响不敏感。
对岩石坝基rp约为3.0.
1.3竖向地震动功率谱及互相关谱对象重力坝一类的工程结构进行二维分析时,除了水平方向地震动还须考虑竖直向地震动作用。
此外,二个方向地震动之间的相关性也应加以考虑。
这里介绍一种简便的处理方法。
竖向地震加速度过程同样可采用金井清模型。
ωg和ξg取值与水平地震加速度相同,唯一待定的参数是竖向地震加速度功率谱强度Sgv。
大量实测地震记录的统计分析结果表明,竖向最大加速度约为水平最大加速度值的1/2~2/3,偏于安全通常取其比值为2/3。
于是有E[agVm]=(2/3)E[agHm],由式(3)、式(4)可得SgV=(4/9)SgH,即竖向地震加速度功率谱强度可取为水平向谱强度的4/9.
水平地震动与竖向地震动之间的相关关系可以通过大量实测地震记录的统计相关分析得到。
在缺乏统计研究结果的情况下,偏于安全可暂设3个方向地震动过程完全线性相关,即agV(t)=(2/3)agH(t).那么,二方向地震动过程互相关函数RgHv=E[agH(t)agV(t+τ)]=2/3RgH(τ),RgH(τ)为水平地震加速度的自相关函数。
因此,互功率谱函数
,式中SgH(ω)为水平地震加速度的功率谱密度函数。
按以上计算,显然有RHV(τ)=RVH(τ)和SgHV(ω)=SgVH(ω).
因此近似地处理,二方向地震加速度过程的互功率谱密度函数也取为金井清谱形式,其谱强度为水平向谱强度的2/3.
2二维结构随机地震动力响应分析
二维结构动力微分方程为
(6)
[M]、[C]、[K]分别为结构质量、阻尼和刚度矩阵;{u}、{
}、{ü}分别为位移、速度和加速度向量;{ag}为基础地震动加速度;[G]称为转换矩阵。
解藕后的运动方程为
(7)
Xj、
、
分别为广义坐标的位移、速度和加速度;ωj、βj、Pj为对应于j阶振型的频率、阻尼比及广义激励。
广义激励
Pj={φj}T[M][G]{ag}
(8)
式中,{φj}为第j阶振型向量。
作傅里叶变换
j阶振型广义坐标的功率谱密度函数
SXj(ω)=|Hj(ω)|2SPj(ω)
(9)
式中Hj(ω)为j阶振型的频率响应函数。
|Hj(ω)|2=1/[(ωj2-ω2)2+4βj2ωj2ω2].
2.1广义激励及广义坐标的功率谱密度函数从式(9)可见,结构随机地震动力分析首先需要计算广义激励的功率谱密度函数。
广义激励的功率谱密度函数
SPjK(ω)={φj}T[M][G][Sg(ω)][G][M]{φK}
(10)
式中
为基础地震动输入加速度功率谱矩阵,矩阵具有对称性。
设[F(ω)]=[M][G][Sg(ω)][G]T[M],则广义激励功率谱密度函数的计算公式化为十分简捷的形式
SPjK(ω)={φj}T[F(ω)]{φK}
(11)
[F(ω)]表示结构的地震动激励,可称为地震动激励功率谱矩阵。
通过矩阵运算可以得到结构地震动激励功率谱矩阵的具体形式如下
(12)
可以看到结构地震动激励功率谱矩阵的元素排列具有很强的规律性,且为对称矩阵。
广义坐标功率谱密度函数为
SXjk(ω)=Hj(-ω)Hk(ω)SPjk(ω)
(13)
广义坐标功率谱密度函数的计算公式具有十分简单的形式,这是结构随机振动理论的特点。
2.2结构随机地震动力响应结构位移分量
uj=[φj]{X}(j=1,2,…,n)
(14)
式中[φj]为振型矩阵的第j行[φj]=[φj1,φj2,…,φjn];{X}为广义坐标列阵。
结构位移响应uj(t)的功率谱密度函数
Suj(ω)=[φj][SX(ω)][φj]T
(15)
式中[SX(ω)]为广义坐标功率谱密度矩阵,其元素由式(13)计算。
结构位移响应均方值
(16)
如果忽略振型之间的相关性,即假设各振型是统计独立的,在这种情况下,结构位移响应的功率谱密度函数计算公式可简化为
(17)
结构位移响应均方值
.
由此,结构位移响应均方根值
(18)
式中σujk2是第k阶振型对应的振型位移响应(第j个位移分量)的均方值。
根据上式,结构位移响应的均方根值为各振型位移响应均方根值的平方和开平方,结构速度、加速度及动应力响应的均方根值也存在类似的振型组合关系。
结构动力响应最大值的期望值与响应均方值之间存在近似的比例关系,因此结构响应最大值可由各振型响应最大值的平方和再开平方获得。
这就是国内外广为使用的振型组合“平方和开平方法”,其理论基础来源于此。
根据结构速度响应、加速度响应与结构位移响应功率谱密度函数之间的转换关系,容易推求结构速度响应及加速度响应功率谱密度函数和响应均方值的计算公式。
3地震动力响应最大值
应用随机振动理论,将地震动输入过程模拟成均值为零的平稳随机过程,进行结构地震动力响应分析,结构响应也是均值为零的平稳过程。
因均值为零,结构随机响应过程的均方值是其最重要的统计特征量。
在平稳高斯过程的假设下,均方值给出了随机响应过程的完整描述。
根据结构响应过程及其导数过程的均方值可以计算结构最大响应的期望值及其方差。
对应于结构随机响应过程,结构地震最大响应为随机变量。
其期望值和方差是描述结构最大响应概率分布特征的主要统计参数。
结构地震最大响应的期望值和方差对结构抗震可靠度分析及结构工程设计具有重要参考价值。
随机过程最大值的计算方法通学采用Davevport1963年提出的一个近似算法。
结构地震响应过程y(t)超过限值b的事件,当b足够大时,可认为是稀有事件,它的概率分布符合Poission分布。
那么有
p(n,T)=[(NbT)n/n!
]exp(-NbT)
(19)
式中p(n,T)代表随机过程y(t)在T时段内n次超过限值b的概率;T为强震持时;Nb表示随机过程单位时间内超过限值b的平均次数。
根据Rice公式
(20)
式中σy和
分别为随机过程y(t)和随机导数过程(t)的均方根值。
由式(19),最大地震响应小于b的概率为
F(b)=p(0,T)=exp(-NbT)
(21)
上式即为最大响应的概率分布函数。
对其进行积分运算可得到最大响应概率密度函数p(b).再由概率统计理论期望值计算公式
,可以得到最大地震响应期望值计算公式为
E[b]=σy(C+0.5772/C)
(22)
式中系数
.根据方差计算公式σy2=E[b2]-(E[b])2,可得到最大地震响应均方差的计算公式
(23)
采用公式(22)和(23)可方便地计算结构最大地震响应的均值及方差,计算公式形式十分简捷。
4算例
选取某实际工程FM混凝土重力坝进行随机地震动力分析。
重力坝坝高77.5m,坝体混凝土质量密度2.35t/m3,动弹模30.0GPa,泊松比0.20;坝基动弹模11.0GPa,泊松比0.27.为消除坝基的放大作用,按无质量坝基进行分析。
库水的影响采用附加质量矩阵计算,库水附加质量按Wetergaard公式换算。
工程设计地震烈度7度,根据我国新的地震烈度表,7度地震对应的地面最大地震动加速度为125cm/s2。
坝基地震动加速度模拟为均值为0的高斯平稳随机过程,坝基地震加速度输入采用金井清过滤白噪声模型,取ωg=5πrad/s,ξg=0.60.按本文介绍的方法计算,水平向地震动加速度谱强度SgH=34.6cm2/s3,竖直向地震动加速度谱强度SgV=15.38cm2/s3,二方向地震动加速度互谱强度SgHV=23.07cm2/s3.
FM重力坝坝体特征部位地震动力响应均方根值及最大值的计算成果列于表1-表3,图1—图3绘出重力坝地震动力响应最大值的分布。
表1FM重力坝坝体地震动位移响应计算成果表
均方根响应
最大响应
特征部位
坝顶
坝踵
坝趾
下游折坡
特征部位
坝顶
坝踵
坝趾
下游折坡
水平
位移/mm
3.090
0.261
0.264
2.172
水平
位移均值/mm
9.942
0.847
0.854
6.975
速度/(mm/s)
63.075
5.872
5.711
43.502
均方差/mm
1.310
0.109
0.111
0.922
竖向
位移/mm
0.840
0.373
0.263
0.302
竖向
位移均值/mm
2.716
1.217
0.854
1.016
速度/(mm/s)
18.110
8.715
5.997
9.908
均方差/mm
0.354
0.156
0.110
0.122
表2FM重力坝坝体地震加速度响应计算成果表
均方根响应
最大响应
特征部位
坝顶
坝踵
坝趾
下游折坡
特征部位
坝顶
坝踵
坝趾
下游折坡
水平
加速度/(m/s2)
1.550
0.461
0.450
1.010
水平
均值/(m/s2)
2.258
1.535
1.478
3.292
导数/(m/s3)
56.00
13.60
11.40
23.20
均方差/(m/s2)
0.620
0.188
0.186
0.423
竖向
加速度/(m/s2)
0.662
0.509
0.462
0.616
竖向
均值/(m/s2)
2.243
1.697
1.522
2.106
导数/(m/s3)
23.80
15.20
12.20
24.60
均方差/(m/s2)
0.264
0.207
0.196
0.244
表3FM重力坝坝体地震动应力响应计算成果表
均方根响应
最大响应
特征部位
坝踵
坝趾
下游折坡处
坝体上部上游面应力集中处
特征部位
坝踵
坝趾
下游折坡处
坝体上部上游面应力集中处
竖向正应力
σy/MPa导数值
0.159
0.175
0.122
0.204
竖向正应力
均值/MPa
1.676
0.564
0.411
0.665
(MPa/s)
11.0
3.64
4.07
4.77
均方差/MPa
0.219
0.074
0.049
0.085
主应力
σ1/MPa导数值
0.549
0.312
0.190
0.205
主应力
均值/MPa
1.773
1.003
1.032
0.668
(MPa/s)
/
/
/
/
均方差/MPa
0.208
0.078
0.047
0.085
由随机地震动力分析得到的FM重力坝地震动位移分布规律,与常规动力分析方法所得规律是一致的。
最大动位移出现在坝顶上游侧,坝顶均方根水平动位移3.090mm,坝顶水平最大动位移均值9.942mm,均方差1.310mm,变异系数0.132.
根据FM重力坝坝体地震加速度随机响应分析结果,自坝基至坝顶,坝体加速度逐渐增大。
在重力坝随机地震动力分析成果中,重力坝的地震放到大效应也得到较好的体现。
坝顶水平均方根加速度1.550m/s2,坝顶水平最大加速度均值5.258m/s2,均方差0.620m/s2,变异系数0.118.根据坝体动应力分布图可以看到,坝体竖向地震动应力在坝踵附近存在明显的应力集中,下游坝面折坡处有一定程度应力集中,坝体上部上游面的应力集中并不太严重。
当然这与具体工程重力坝体型设计有关,坝体体型设计合理可使坝体上部的动应力集中现象得到一定程度的减轻。
竖向动应力在坝踵附近出现最大值,坝踵竖向应力均方根值为0.519MPa,其最大值期望值为1.676MPa,均方差0.219MPa,变异系数0.131.
图1最大水平动位移均值(mm)分布
图2最大水平加速度均值(m/s2)分布
图3最大竖向动应力均值(MPa)分布
5结语
本文论述了二维重力坝结构的随机地震动力分析方法,提出了一个考虑水平和竖向二个方向地震动的二维随机地震动输入实用模型。
模型使用方便,并能适当体现二方向地震动的相关关系。
针对二维结构情形,本文从实际应用角度详细推导了结构随机地震动力分析计算公式,给出了二维地震激励功率谱矩阵的具体形式,为实际运用提供了方便。
某实际工程重力坝随机地震动力响应分析成果表明,采用随机振动理论进行结构动力分析,其成果能够较好地体现结构地震响应特征,并能进一步地反映工程结构地震响应的统计规律。
参考文献:
[1](日)金井清著。
工程地震学[M]。
常宝琦、张虎男译,北京:
地震出版社,1987.
[2]DavenportAG.NoteontheDistributionoftheLargestValueofaRandomFunctionWithApplicationtoGustLoading[J]。
Proc.Inst.CivilEngr1963,28.
[3]何蕴龙。
混凝土重力坝随机地震响应及其动力可靠性研究[D]。
武汉:
武汉水利电力大学,1996.
网
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 地面 运动 具有 强烈 随机性 不确定性
![提示](https://static.bdocx.com/images/bang_tan.gif)