多点地震动时程人工合成
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
* " + * ) , & + 比较 有 代表性的方法 有 ! 等将 ) tr e 和屈铁军 * + , ’ + ) 是一致的 2 -. r {r u r \ 等推导出了一种不需要 做g 通过调整求和项 p e | _ b & w分解的渐进合成方法 ,
-! ^ , c ‘ ^ ~^ b ‘ _ v d ^ ‘ _ u c pc p _ \ d {_ ‘ ^ m r | b c d u wf e ‘ c p _ _ ^ ] _ \ k r | d _ b r \ u_ ^ ] _ \ k _ m c e ‘ b e f | r ] ] _ um e p _ ‘ _ \ m w{r c ‘ ^ ~
F .
导出一种在任一频率点处的时程幅值与迟滞相干矩 阵 特征值和特征 向 量 有 关 的 方 法 ! 研究了相干矩阵 特 征值和特征向 量 的 变 化 规 律 ! 提出了对其进行简 化的方法 ! 用算例验证了简化计算的可行性 "
0@ JK HK 2 ’ R ) 2I HK @ J HO @ J
# 平稳地震动时程合成公式推导
!x ,rb 5 y z 7 6 ; : 7 r b _ ue \c p _b c r c ^ e \ r ‘ w‘ r \ u e {h ‘ e m _ b bc p _ e ‘ w ^ {d | r c ^ e \{_ c p e uf e ‘ b h r c ^ r | m e ‘ ‘ _ | r c ^ k _c ^ {_
摘
要 !分析大跨度结构需要一种效率 1 精度都较高的多点地震动时程合成方法 2 基于平稳随机过程理论推
导了多点地震动时程合成公式 , 每一频率点处的时程幅值与迟滞相干矩阵特征值和特征向量有关 2 用数值计 算研究了迟滞相干矩阵特征对的变化规律 , 提出了对其进行简化的方法 , 用算例验证了简化计算的可行性 2 关键词 !多点输入 / 地震动时程 / 振型分解 中图分类号 !3 & " + 文献标识码 !4
) ,x # # 0 ’ " ,g * o p _t_ r u v d r ‘ c _ ‘ be f c p _3 i 4 qr k w _ ^ s ^ \ ]" p ^ \ r
" " )
( " ,3 -e -n o ,qr " # # # . ,g / [ \ ] ^ \ _ _ ‘ ^ \ ]a \ b c ^ c d c _e f [ \ ] ^ \ _ _ ‘ ^ \ ]g e ‘ h b i 4 j\ ^ k f l m ^ _ m p \ s ^ \ ]) p ^ \ r
’ . )
’ S ) 3 $ -3 ? = C 8 5 6 7 ’ 4 ( ( ) ;* 2 1 $9 : -1 0 4 1 ?0 : = 1 @ 当 -*= 时%
其中 ! 是第 $ 点时程在频率 4 ) 3 4 $ -’ 1 1 处的第 - 个分 量! 点 的时 间 ! 它是 ( $是地 震波从 某一参 考点 到达 $ 只与参考点到 $ 点的距离有关 ! 是 一个确定的量 ! : -1 当 -<= 或1 时! 和: 相互独立 " 随机相位角 ! <> : -1 = > 现 在通过确定 3 的值使合成的时程满足给定 ) 4 $ -’ 1 的谱密度矩阵 " 第$ 点的时程 & 的第 1个窄 频 带 信号 & ’ ) ’ ) ( ( $ $ 1 与第 ? 点的时程 & 的第 1个窄频带信号 & 的 ’ ) ’ ) ( ( ? ? 1 表达式为 %
,
’ . . )
’ 2 ) *C 8+ 3 ) ’ 94 9 B 4 5 6 7 4 ( ( $ ? 1 $ -’ 1 1 1 $
-* . ,
& ’ ( )* $ 1 & ’ ( )* ? 1
4) 5 6 7 ’ 4( 9 4(9 : ) ’ @ ) +3 ’
$ 1 1 1$ -1 -* . ,
)+ 3 ’ ) ’ 9: ) ; * 94 4 5 6 7 4 : ( ( -1 ? = 1 1 1 ? = 1
’ + 用 频 率 波 数 谱 表 示 的 合 成 公 式* 其本质上与方法 ,
关性符 合预 定值 , 随机场变界法也只能保证在频率 而其余各 频率 各 2 3 处选取点间的相关性得到保证 , 间 距点的 相关 性 与 给 定 的 相 关 值 不 符 , 地震动相干 函数不能精确满足 2 上 述 几 种 方 法 中, 方法)和方法,合成精度较 高, 但 效率 不 高 / 方 法 / 和 方 法 1 效 率 较 高, 但精度 不高 / 方法 - 精度和效率都较好 , 但只适用于相干函 数比较简单的情况 2 本文基于平稳随机过程理论推
,
’ . 2 )
3 $ -3 ? = 5 6 7 ’ @ 4 ( 94 ( ( ) ;9 C 8 1 1 $9 : -1 9 4 1 ?9 : = 1 @ 3 $ -3 ? = C 8 5 6 7 ’ 4 ( ( ) ;* 1 $9 : -1 0 4 1 ?0 : = 1 @ 3 $ -3 ? 5 6 7 ’ 4 ( ( ) 1 $0 4 1 ? @
% 5 6 7 8 9 8 : 8 ; < = > ? > 6 ; 7 8 @ ?@ 9 AB < 7 8 C @ 8 ? 7 D ; 6 7 E F B ; G >= 6 @ B ? HA@ 7 8 @ ?
% , RKSTSQ % % IJ KL M N O P Q P U V Q P M, L KS W X P U Y Z X N
-* . ,
,
它们的协方差为 % B ’ 2 )* C 8 & ’ ( ) & ’ ( ) ; $ ? 1 $ 1 ? 1 将式 ’ 式’ 代入式 ’ 得到 % @ ) E A ) D )
,
’ . @ )
’ D )
B ’ 2 )* @ W 4 T U 8 X ’ V 4 ) ; $ ? 1 $ ? 1 ’ 4 ) X ’ V 4 ) Z MX . . . @ ’ ) * Y V 4 X ’ V 4 ) X ’ 4 ) Z @ . @ @ [ [ \ X ’ V 4 ) ] . , X ’ V 4 ) @ , [
’ . F )
其中 % 8 U
0V 4 _
; *a 8 V b c U
0V 4 ( .
! U
0 , $ + 数 21 地震动随机场的投影展开法 * 任意选取 S , 个随 机变量 确定 投 影 基 , 只能保证选取点之间的相
地震动加速度谱密度矩阵做 g 得到某 p e | _ b & w分解 , 点地震动时程不同频率点处幅值和相位与加速度谱 密 度矩阵分解值 的 关 系 , 通过对有正态分布随机相 位角的余弦函数双求和 , 得到各点的地震动时程 2 , 将地 震 动场视为 均 值 为 零 的 空 间 时 间 随 机 场 , 得到
. ; 度函数分别为 8 %
’ V 4 ) X ’ V 4 ) Z X ’ 4 ) NX ^ , . , @ 百度文库 , 假 定各支点处的点谱相同为 X 则 式 可表 ’ ) ! ’ . D ) 4 2 示为 %
0V 4 _ ‘ 0V 4 _ Y ’ V 4 )* X ’ 4 ) 8 U ;B 8 U ; 2
= *.
’ . D )
式’ 中任一项的集合平均为 % F ) ) ’ 94 9: ) ’ ) ’ 94 9: ) ; * C8 3 4 5 6 7 4 ( ( 3 4 5 6 7 4 ( ( $ -’ 1 1 1 $ -1 ? = 1 1 1 ? = 1 3 $ -3 ? = 8 ’ @ 94 9: 94 9: ) ; 9 C 5 6 7 4 ( ( ( 1 1 $ -1 1 ? = 1 @ 3 $ -3 ? = 8 ’ 9: 04 0: ) ; ’ G ) 5 6 7 4 ( ( C 1 $ -1 1 ? = 1 @ 将: 和: 看 成 两 个 随 机 变 量! 它们的密 9: 0: -1 = 1 -1 = 1
= *.
4) 5 6 7 ’ 4( 9 4(9 : ) +3 ’
? = 1 1 1? = 1 = *.
’ A )
3 $ -3 ? ’ 04 ) * 5 6 7 4 ( ( 1 $ 1 ? + @ -* . 3 $ -3 ? - V ’ 4( 0 4 () U 1$ 1? ; @ 根据互功率谱与协方差的关系可得到 % 8+ T U
得到了合成 的幅值 来满足 对 互 谱 密 度 矩 阵 的 要 求 , 公式中各系数的解析表达式 , 合成效率大大提高 2 /
* . + 用 4. 但它 04 方 法 合 成 多 点 平 稳 地 震 动 时 程 , 的不 足在于 需要 预 先 知 道 地 震 动 的 互 协 方 差 函 数 ,
而对 地震动 而言 , 通常给出的是自功率谱和相干函
第 &卷 第 &期 ) # # )年 ’月
解 放 军 理 工 大 学 学 报 ( 自 然 科 学 版* 4 e d ‘ \ r | e f 3 i 4 j\ ^ k _ ‘ b ^ c we f l m ^ _ \ m _r \ uo _ m p \ e | e ] w
5e | &qe & 4 d \ ) # # )
’ . A )
而 ,个点处地震动时程的互功率谱密度矩阵为 %
B ’ 2 )* C 8 4 ) 5 6 7 ’ 4 ( 94 ( $ ? 1 $ -’ 1 1 1 $9 +3
-* . ,
: ) ’ 4 ) 5 6 7 ’ 4 ( 94 ( ) ;’ F ) -1 ? = 1 1 1 ?9 : = 1 +3
将第 $ 点的时程表示成如下形式 %
, /0 .
3 $ -3 ? = 5 6 7 ’ @ 4 ( 94 ( ( ) ;* 2 C 8 1 1 $9 : -1 9 4 1 ?9 : = 1 @
-1
& ’ ( )* $
4) 5 6 7 8 4( 9 4(9 : ; ++3 ’
$ 1 1 1$ -* . 1 *2
% -}\ p ^ b c e ‘ ^ _ be f {d | c ^ h e ^ \ c ] ‘ e d \ u{e c ^ e \^ b ] ^ k _ \ | wc p _ {e u r | u _ m e {h e b ^ c ^ e \e f c p _ | r ] ] _ um e p _ ‘ _ \ m w{r % -4\_ rb ^ {h | ^ f ^ _ u{_ c p e u^ bh ‘ e h e b _ u ~ r {h | _^ | | d b c ‘ r c _ bc p _f _ r b ^ " ^ | ^ c we f c p _b ^ {h | ^ f ^ m r c ^ e \ !{d % /_ /{e # > $%@ 6 H z | c ^ h e ^ \ c ^ \ h d c r ‘ c p v d r & _] ‘ e d \ u{e c ^ e \ u r | u _ m e {h e b ^ c ^ e \ 多点输入的研究对象是平面尺寸较大的结构 物, 担负着重要的社会职能 , 设计这些结构时一般要 采用直接积分法计算结构的反应 2 由于台阵记录的 数据有限 , 且结构物的平面尺寸多种多样 , 现有的差 动 台 阵记录不可 能 满 足 多 点 输 入 抗 震 设 计 要 求 , 需 要有比较合理的多点地震动时程合成方法 2 人工合成多点地震动时程 , 不少学者做过研究 ,
收稿日期 万方数据 ! ) # # " % # ’ % ) ’ 作者简介 ! 夏友柏 ( 男, 博士生 " $ ’ ’ ( * ,
第 A期
夏友柏 ! 等 %多点地震动时程人工合成 M2 . ’ H9 @ J ) @ D J : ’ H )* L @ . ’ @ Q0 H ) @ D J N2 当 -<= 时% HI0 @ J
文章编号 !" # # $ % & ’ ’ & ( ) # # ) * # & % # # + # % # ’
多点地震动时程人工合成
夏友柏 " , 王年桥 " , 鄢常舒 )
解放军理工大学 工程兵工程学院 , 江苏 南京 ) 解放军海军司令部 , 北京 " ( " " # # # . /) # # 0 ’ " *
-! ^ , c ‘ ^ ~^ b ‘ _ v d ^ ‘ _ u c pc p _ \ d {_ ‘ ^ m r | b c d u wf e ‘ c p _ _ ^ ] _ \ k r | d _ b r \ u_ ^ ] _ \ k _ m c e ‘ b e f | r ] ] _ um e p _ ‘ _ \ m w{r c ‘ ^ ~
F .
导出一种在任一频率点处的时程幅值与迟滞相干矩 阵 特征值和特征 向 量 有 关 的 方 法 ! 研究了相干矩阵 特 征值和特征向 量 的 变 化 规 律 ! 提出了对其进行简 化的方法 ! 用算例验证了简化计算的可行性 "
0@ JK HK 2 ’ R ) 2I HK @ J HO @ J
# 平稳地震动时程合成公式推导
!x ,rb 5 y z 7 6 ; : 7 r b _ ue \c p _b c r c ^ e \ r ‘ w‘ r \ u e {h ‘ e m _ b bc p _ e ‘ w ^ {d | r c ^ e \{_ c p e uf e ‘ b h r c ^ r | m e ‘ ‘ _ | r c ^ k _c ^ {_
摘
要 !分析大跨度结构需要一种效率 1 精度都较高的多点地震动时程合成方法 2 基于平稳随机过程理论推
导了多点地震动时程合成公式 , 每一频率点处的时程幅值与迟滞相干矩阵特征值和特征向量有关 2 用数值计 算研究了迟滞相干矩阵特征对的变化规律 , 提出了对其进行简化的方法 , 用算例验证了简化计算的可行性 2 关键词 !多点输入 / 地震动时程 / 振型分解 中图分类号 !3 & " + 文献标识码 !4
) ,x # # 0 ’ " ,g * o p _t_ r u v d r ‘ c _ ‘ be f c p _3 i 4 qr k w _ ^ s ^ \ ]" p ^ \ r
" " )
( " ,3 -e -n o ,qr " # # # . ,g / [ \ ] ^ \ _ _ ‘ ^ \ ]a \ b c ^ c d c _e f [ \ ] ^ \ _ _ ‘ ^ \ ]g e ‘ h b i 4 j\ ^ k f l m ^ _ m p \ s ^ \ ]) p ^ \ r
’ . )
’ S ) 3 $ -3 ? = C 8 5 6 7 ’ 4 ( ( ) ;* 2 1 $9 : -1 0 4 1 ?0 : = 1 @ 当 -*= 时%
其中 ! 是第 $ 点时程在频率 4 ) 3 4 $ -’ 1 1 处的第 - 个分 量! 点 的时 间 ! 它是 ( $是地 震波从 某一参 考点 到达 $ 只与参考点到 $ 点的距离有关 ! 是 一个确定的量 ! : -1 当 -<= 或1 时! 和: 相互独立 " 随机相位角 ! <> : -1 = > 现 在通过确定 3 的值使合成的时程满足给定 ) 4 $ -’ 1 的谱密度矩阵 " 第$ 点的时程 & 的第 1个窄 频 带 信号 & ’ ) ’ ) ( ( $ $ 1 与第 ? 点的时程 & 的第 1个窄频带信号 & 的 ’ ) ’ ) ( ( ? ? 1 表达式为 %
,
’ . . )
’ 2 ) *C 8+ 3 ) ’ 94 9 B 4 5 6 7 4 ( ( $ ? 1 $ -’ 1 1 1 $
-* . ,
& ’ ( )* $ 1 & ’ ( )* ? 1
4) 5 6 7 ’ 4( 9 4(9 : ) ’ @ ) +3 ’
$ 1 1 1$ -1 -* . ,
)+ 3 ’ ) ’ 9: ) ; * 94 4 5 6 7 4 : ( ( -1 ? = 1 1 1 ? = 1
’ + 用 频 率 波 数 谱 表 示 的 合 成 公 式* 其本质上与方法 ,
关性符 合预 定值 , 随机场变界法也只能保证在频率 而其余各 频率 各 2 3 处选取点间的相关性得到保证 , 间 距点的 相关 性 与 给 定 的 相 关 值 不 符 , 地震动相干 函数不能精确满足 2 上 述 几 种 方 法 中, 方法)和方法,合成精度较 高, 但 效率 不 高 / 方 法 / 和 方 法 1 效 率 较 高, 但精度 不高 / 方法 - 精度和效率都较好 , 但只适用于相干函 数比较简单的情况 2 本文基于平稳随机过程理论推
,
’ . 2 )
3 $ -3 ? = 5 6 7 ’ @ 4 ( 94 ( ( ) ;9 C 8 1 1 $9 : -1 9 4 1 ?9 : = 1 @ 3 $ -3 ? = C 8 5 6 7 ’ 4 ( ( ) ;* 1 $9 : -1 0 4 1 ?0 : = 1 @ 3 $ -3 ? 5 6 7 ’ 4 ( ( ) 1 $0 4 1 ? @
% 5 6 7 8 9 8 : 8 ; < = > ? > 6 ; 7 8 @ ?@ 9 AB < 7 8 C @ 8 ? 7 D ; 6 7 E F B ; G >= 6 @ B ? HA@ 7 8 @ ?
% , RKSTSQ % % IJ KL M N O P Q P U V Q P M, L KS W X P U Y Z X N
-* . ,
,
它们的协方差为 % B ’ 2 )* C 8 & ’ ( ) & ’ ( ) ; $ ? 1 $ 1 ? 1 将式 ’ 式’ 代入式 ’ 得到 % @ ) E A ) D )
,
’ . @ )
’ D )
B ’ 2 )* @ W 4 T U 8 X ’ V 4 ) ; $ ? 1 $ ? 1 ’ 4 ) X ’ V 4 ) Z MX . . . @ ’ ) * Y V 4 X ’ V 4 ) X ’ 4 ) Z @ . @ @ [ [ \ X ’ V 4 ) ] . , X ’ V 4 ) @ , [
’ . F )
其中 % 8 U
0V 4 _
; *a 8 V b c U
0V 4 ( .
! U
0 , $ + 数 21 地震动随机场的投影展开法 * 任意选取 S , 个随 机变量 确定 投 影 基 , 只能保证选取点之间的相
地震动加速度谱密度矩阵做 g 得到某 p e | _ b & w分解 , 点地震动时程不同频率点处幅值和相位与加速度谱 密 度矩阵分解值 的 关 系 , 通过对有正态分布随机相 位角的余弦函数双求和 , 得到各点的地震动时程 2 , 将地 震 动场视为 均 值 为 零 的 空 间 时 间 随 机 场 , 得到
. ; 度函数分别为 8 %
’ V 4 ) X ’ V 4 ) Z X ’ 4 ) NX ^ , . , @ 百度文库 , 假 定各支点处的点谱相同为 X 则 式 可表 ’ ) ! ’ . D ) 4 2 示为 %
0V 4 _ ‘ 0V 4 _ Y ’ V 4 )* X ’ 4 ) 8 U ;B 8 U ; 2
= *.
’ . D )
式’ 中任一项的集合平均为 % F ) ) ’ 94 9: ) ’ ) ’ 94 9: ) ; * C8 3 4 5 6 7 4 ( ( 3 4 5 6 7 4 ( ( $ -’ 1 1 1 $ -1 ? = 1 1 1 ? = 1 3 $ -3 ? = 8 ’ @ 94 9: 94 9: ) ; 9 C 5 6 7 4 ( ( ( 1 1 $ -1 1 ? = 1 @ 3 $ -3 ? = 8 ’ 9: 04 0: ) ; ’ G ) 5 6 7 4 ( ( C 1 $ -1 1 ? = 1 @ 将: 和: 看 成 两 个 随 机 变 量! 它们的密 9: 0: -1 = 1 -1 = 1
= *.
4) 5 6 7 ’ 4( 9 4(9 : ) +3 ’
? = 1 1 1? = 1 = *.
’ A )
3 $ -3 ? ’ 04 ) * 5 6 7 4 ( ( 1 $ 1 ? + @ -* . 3 $ -3 ? - V ’ 4( 0 4 () U 1$ 1? ; @ 根据互功率谱与协方差的关系可得到 % 8+ T U
得到了合成 的幅值 来满足 对 互 谱 密 度 矩 阵 的 要 求 , 公式中各系数的解析表达式 , 合成效率大大提高 2 /
* . + 用 4. 但它 04 方 法 合 成 多 点 平 稳 地 震 动 时 程 , 的不 足在于 需要 预 先 知 道 地 震 动 的 互 协 方 差 函 数 ,
而对 地震动 而言 , 通常给出的是自功率谱和相干函
第 &卷 第 &期 ) # # )年 ’月
解 放 军 理 工 大 学 学 报 ( 自 然 科 学 版* 4 e d ‘ \ r | e f 3 i 4 j\ ^ k _ ‘ b ^ c we f l m ^ _ \ m _r \ uo _ m p \ e | e ] w
5e | &qe & 4 d \ ) # # )
’ . A )
而 ,个点处地震动时程的互功率谱密度矩阵为 %
B ’ 2 )* C 8 4 ) 5 6 7 ’ 4 ( 94 ( $ ? 1 $ -’ 1 1 1 $9 +3
-* . ,
: ) ’ 4 ) 5 6 7 ’ 4 ( 94 ( ) ;’ F ) -1 ? = 1 1 1 ?9 : = 1 +3
将第 $ 点的时程表示成如下形式 %
, /0 .
3 $ -3 ? = 5 6 7 ’ @ 4 ( 94 ( ( ) ;* 2 C 8 1 1 $9 : -1 9 4 1 ?9 : = 1 @
-1
& ’ ( )* $
4) 5 6 7 8 4( 9 4(9 : ; ++3 ’
$ 1 1 1$ -* . 1 *2
% -}\ p ^ b c e ‘ ^ _ be f {d | c ^ h e ^ \ c ] ‘ e d \ u{e c ^ e \^ b ] ^ k _ \ | wc p _ {e u r | u _ m e {h e b ^ c ^ e \e f c p _ | r ] ] _ um e p _ ‘ _ \ m w{r % -4\_ rb ^ {h | ^ f ^ _ u{_ c p e u^ bh ‘ e h e b _ u ~ r {h | _^ | | d b c ‘ r c _ bc p _f _ r b ^ " ^ | ^ c we f c p _b ^ {h | ^ f ^ m r c ^ e \ !{d % /_ /{e # > $%@ 6 H z | c ^ h e ^ \ c ^ \ h d c r ‘ c p v d r & _] ‘ e d \ u{e c ^ e \ u r | u _ m e {h e b ^ c ^ e \ 多点输入的研究对象是平面尺寸较大的结构 物, 担负着重要的社会职能 , 设计这些结构时一般要 采用直接积分法计算结构的反应 2 由于台阵记录的 数据有限 , 且结构物的平面尺寸多种多样 , 现有的差 动 台 阵记录不可 能 满 足 多 点 输 入 抗 震 设 计 要 求 , 需 要有比较合理的多点地震动时程合成方法 2 人工合成多点地震动时程 , 不少学者做过研究 ,
收稿日期 万方数据 ! ) # # " % # ’ % ) ’ 作者简介 ! 夏友柏 ( 男, 博士生 " $ ’ ’ ( * ,
第 A期
夏友柏 ! 等 %多点地震动时程人工合成 M2 . ’ H9 @ J ) @ D J : ’ H )* L @ . ’ @ Q0 H ) @ D J N2 当 -<= 时% HI0 @ J
文章编号 !" # # $ % & ’ ’ & ( ) # # ) * # & % # # + # % # ’
多点地震动时程人工合成
夏友柏 " , 王年桥 " , 鄢常舒 )
解放军理工大学 工程兵工程学院 , 江苏 南京 ) 解放军海军司令部 , 北京 " ( " " # # # . /) # # 0 ’ " *