地震与扭波

地震与扭波
地震与扭波

地震波发生于地下,通过岩土质点的弹性振动进行传播,当地震波传到地面附近,受震物体有了扩展空间,则形成集团运动,使浅部地层发生旋扭,使地面上的物体产生三维扭动,地震仪所记录的只是地震仪所在地的物体集团运动的强度和周期变化,它既测不出集团运动的三维扭动轨迹,更测不到真正的地震波,由此得出一个错误理论,认为地震是先上下跳动(纵波所为),然后水平单向摇摆(横波所为),并据此采用了摇床作地震模拟实验,进行理论研究,皆与实际情况相悖。

在无法直接测出集团扭动三维轨迹之前,只能根据地震破坏遗迹进行分析,获取他们在破坏之前的运动信息。1976年7月28日唐山7.8级地震后,我们对唐山震害进行了大量实地观察,结合海城地震、溧阳地震和龙陵地震资料,总结出地震破坏的十大规律,搞清了地震破坏机制,并在此基础上,通过应力分析,发现了扭波,进而提出针对扭波设防抗震的八项措施。1979年发表了第一篇论文,题为"扭波与地震",此后又发表了多篇论文,有几十家报刊杂志作了报导。由于人所共知的原因,扭波理论至今仍未被权威专家所接受。

1995年1月17日日本神户发生了7.2级大地震,一家电视台的监视器录下了该办公室遭地震破坏的全过程,桌子、椅子都在旋扭,最大转角达170度。证明了我们在二十年前根据颓垣断壁分析出的扭波正确无误。可是在二十年后,日本和美国的地震

学家,在物体受震发生扭转运动的实况录像面前,仍然把地震破坏说成是由纵波和横波造成的。实在遗憾。为了改变这一不正常状况,现将1992年发表的一片论文"扭波的发现与设防抗震"重新发表,以正视听。

地震机制--扭波的发现与设防抗震原载于《新疆地质》增刊1992

目录

前言

一、地震破坏的规律和存在的问题

1、地震破坏的若干规律

(1)线性物体直立如初

(2)钢筋铁架抗震最佳

(3)裂开旋扭同时并存

(4)破裂形式如同老化

(5)间隔破裂轻重分带

(6)断裂两盘不受损伤

(7)地面平整,起伏罕见

(8)四面开花铅直坠落

(9)地下建筑保存良好

(10)人体抗震胜过厂房

2、存在的问题

(1)地震破坏的10条规律与纵波、横波及面波理论相悖(2)五个关键问题

二、弹性波的成生机制

1、弹性振动

2、弹性波

(1)弹性波的成生机制

(2)弹性波的传播

(3)弹性场与弹性波转向

(4)横扭波与纵扭波

三扭波与地震破坏

1 扭波的破坏作用

2 扭断类型

横扭裂

纵扭裂

3 地震破坏的分带性

4 扭波破坏与重力关系

5 地旋

四针对扭波设防抗震

1 截波

2 活端

3 底小

4 加筋

5 结牢

6 轻质

7 加箍

8 地下

参考文献

前言

中国是最早记录地震和记录古代地震最多的国家,也是首先发明地震仪和首先进行地震预报的国家。晋太康三年(282年)在魏襄王(公元前318-前296年)的墓中,出土了一部竹简,竹简上刻写着上至皇帝,下迄魏襄王二十二年(公元前297年)的历史,人称《竹书纪年》,里面记有4次地震,其中最早的一次是夏帝发七年(公元前1831年)泰山震。在这以后的许多有关古籍中,都有对地震的记载,从公元前1831-1955年3800年中,大约记录了9000次地震,其中约有1000次是破坏性地震。经整理不仅勾划出我国的地震区,同时还揭示出地震活跃期和平静期的周期性变化规律。

在国外,最早作地震记录的是日本,其记载是允恭天皇五年(416年)日本河内地震,比中国晚2247年。

东汉著名的科学家张衡(78-139年)首创了世界上第一台地震仪"候风地动仪",安放在洛阳城内,正确的报出了东汉永和三年(138年)二月初三在陇西所发生的地震,陇西位于今甘肃兰州、临洮一带,在洛阳以西700km。而欧洲在1880年才制成

地震仪,其原理和张衡地动仪基本相似,比中国晚1700多年。

在数千年以前,我国就力图从自然界本身去找地震原因。《竹书记年》在记述夏帝发十年的地震时写道"五星错行,夜中陨星如雨,地震。"春秋战国时期,阳伯父认为"天地之气,不失其序……阳伏而不能出,阴迫而不能蒸,于是有地震"。《庄子》则提出海水周流相薄说,认为:"海水三岁一周,流波相薄,故地动"。东汉杰出的思想家王允,最早提出了地震是地壳本身的"自动"现象,他不止一次的说:"地固将自动"。沈括认为,天地之变"率皆有法","物理有常,有变"。

根据地震前兆预报地震,在我国历史上也有许多记录,并总结出一套预报地震的经验。据山西省《虞乡县志》记载,清嘉庆二十年(1815)年,山西平陆地区曾根"淫雨后大热,宜防地震"的经验,预报出九月二十日午夜二时发生的一次强烈地震,据《清代地震档案史料》记载,清咸丰四年(1855年)十二月十一日,辽宁金县旗民根据地声,预报出即将发生的一次破坏性地震:"未震之时,先闻有声如雷,故该处旗民早已预防,俱各走避出屋,以未曾压毙多人,只伤男妇子女共七名"。宁夏《隆德县志》曾把地震前兆归纳为"地震六端":一、井水本湛静无波,倏忽挥如墨汁,泥渣上浮,势必地震。二、池沼之水,风吹成荇交萦,无端泡沫上腾,若沸煎茶,势必地震。三、海面遇雨,波浪高涌,奔腾洴汹,此常情;若日晴和,台飓不作,海水忽然绕起,汹涌异常,势必地震。四、夜半晦黑,天忽开朗,光明照耀,光异日

中,势必地震。五、天晴日暖,碧空清净,忽见黑云如缕,蜿入长蛇,横亘空际,久而不散,势必地震。六、时置盛夏,酷热蒸腾,挥汗如雨,蓦觉清凉,如受冰雪,冷气袭人,肌为之粟,势必地震。这6条对地震前的天气异常,海啸、地光、地震云等地震前兆都作了精辟的概括。日本人最近几年才开始研究地震云,所得结果与以上描述基本一致。

自1996年以后,我国已报准18次破坏性地震,如1975年的辽南海城地震(7.3级),1976年的云南龙陵,路西地震(7.5-7.6级)和四川松潘地震(7.2级)。这在世界上也是屈指可数的。从世界水平来看,地震机制问题还没有完全得到解决,根据表面现象作地震预报,只能对部分地震有效,不可能把所有地震都预报出来,1976年唐山地震未能报出,是很自然的,在地震研究方面,还有一大段路要走。

地震波有3种,除众所周知的纵波和横波以外还有扭波。扭波是地震破坏的主要动力,扭波的发现揭开了地震破坏的内幕,为制定抗震措施提供了理论依据。

一、地震破坏的规律和存在的问题

70年代中国接连发生了十几次强烈地震,生命财产遭到严重破坏,仅1976年唐山7.8级地震就死亡24.2万人,重伤16.4万人,从而引起人们对地震的强烈关注。唐山地震后,在观察中发现《地震烈度表》〔1〕所附的样板图样(图版 1)与事实有很大出入。例如,样板图上面的是电线杆和楼房随地震烈度的增加

逐渐倾斜,倾角逐渐增大XI-XII度时全部倒塌。实际上唐山地震并非如此,唐山路南区的破坏最严重,房屋全糟破坏,大多数塌平,按其它指标烈度定为XI度,可是电线杆仍然直立如初(图版2,4,5)〔2〕,另外无论是7层大楼还是平顶小屋坍塌后留下的断墙残壁基本上都是直的(图版 3)〔2〕,并不是先倾斜后破裂。这说明,人们对地震破坏的认识还有一定差距,有进一步研究的必要。通过对现场调查,结合溧阳地震、海城地震和龙陵地震,我们初步摸到了一些规律。

1、地震破坏的若干规律

?线性物体直立如初

唐山市和丰南县地震烈度为X-XI度,厂房和住宅几乎全部破坏,但其中有许多线性物体被完整地保存下来。例如:

①所有树木均安然无恙,均未直接受害(图版 2,4,5)。

②绝大多数电线杆直立如初(图版 2,5;图版 6,7)。

唐山河北矿治学院翻砂车间是钢筋水泥厂房。地震后柱断、梁弯、顶塌,但厂房门前的两根电线杆直立如初,房后的树木也依然如故(图版 2)。

唐山开滦医院原是新建的7层大楼,地震时大部分塌平,但紧靠它的树木和电杆仍安然无恙。残留的一角也是直立的(图版3)。

40cm宽的钢筋混凝土房柱被折断,而在同一地区的电杆和树木安然无恙(图版 7)。

唐山XI度区胜利桥墩倒梁落,而桥旁树木、电杆完好如初(图版 5)。

③唐山65m高的微波转播铁塔巍然屹立在大片废墟之中(图版 8);震后仍能使用(潜望镜的方向性要求很高不能相差几秒)。这说明地震对它几乎无丝毫影响。

④40m高的用钢筋混凝土建造的唐山车站水塔依然如故。

⑤天津化工厂厂房坍塌,而直立的反应塔和高大烟囱及电杆均保存完好。

在线性物体中也有破坏较重的,主要是没有钢筋的砖砌烟囱和砖砌水塔,沿个别砖缝破裂后引起整段脱落。

?钢筋铁架抗震最佳

①唐山两座铁路钢架天桥和开滦煤矿的十几对钢架竖井均完好无损(图版 9)。

②唐山路南区65m高钢架电视转播塔屹立如初(图版 8)。

③唐山启新水泥厂办公楼砖砌窗柱破碎散落,而柱内钢筋完好,支撑着上层建筑(图版 10)。

④唐山一些早年电杆,受震后水泥破裂脱落,钢筋无恙,电线杆仍然直立(图版 11)。

?裂开旋扭同时并存

无论是墙裂缝、地裂缝还是烟囱裂缝,凡是裂缝都有旋扭,以柱状体的断裂扭转最为明显。例如:

①唐山启新水泥厂的一栋三层库房,一楼和二楼基本完好,

三楼的所有窗柱全部断裂,旋转方向和角度各不相同,现存旋转角度最大的一个右转40o(图版 12);旋转角更大的当时即已脱落,无从查考。

②唐山河北矿冶学院翻砂车间,各墙柱均断裂旋转,旋转方向和旋转角度各个不同(图版 2)。

③唐山文化宫门柱断裂旋转(图版 13)。

④唐山十中地裂缝右行旋转1.2m(图版 14)。

⑤唐山文化宫烟囱水平旋转(图版 15)。

⑥辽宁海城县委会烟囱断成4节,最上一节右旋40o(图版16)。

⑦辽宁营口市砖砌烟囱水平断成5节,皆扭转错位(图版17)。

⑧河北宁河火车站局部铁轨发生水平扭转,其他铁轨和电杆均无损伤(图版 18)[2]。

⑨河北丰南火车站一条铁轨发生反复扭曲,右侧的铁轨和站台也发生变形,其它铁轨以及树木和电杆均保存完好(图版 19)[2]。

⑩云南龙陵镇安中学1976年7.3级地震后两个篮球架都发生了水平旋转。一个左旋14o,一个左旋18o,最大平移1.22m (图版 20)[2]。

(11)唐山地震京榆公路滦河大桥桥面塌落同时旋扭(图版21)[2]。

?破裂形式如同老化

地震造成的破坏主要是沿着物体本身的脆弱带发育,不论是地面上的建筑物,地表的道路还是地下的地层都是如此。例如:

①沿建筑物四角断裂(图版 12;图版 13)。

②沿门窗四角向外裂开(图版 22;图版 12)。

③沿支柱顶底面切断扭转(图版 13)。

④平行砖缝断裂(图版 15,16,17)。

⑤平行或垂直马路延长方向断裂(图版 23)。

⑥垂直水泥场地的边界线断裂。

⑦平行河边、塘边断裂。

地震所造成的断裂形式与非地震的"自然老化"和人为破坏所造成的断裂基本一致,地震只是加快了老化断裂的发展速度,是急性发作。它的破坏形式与震源方位无直接关系。地裂缝的发展也象建筑物裂缝一样只受控于接近地表的局部基底结构,沿其脆弱带产生地裂缝,其分布和延展方向也与震源和发震构造无直接关系。地裂缝都是比较浅的,通常只有几米,深十米以上的尚未见报导。

?间隔破裂轻重分带

①水平分带

唐山公安学校有3栋三层楼房。形状相同,相互间隔10m平行排列,南面一栋完全塌平,中间一栋保留了中段的部分骨架,北面一栋只前沿散落。

唐山水泥设计院新建的5栋楼房,均为砖混结构。同样设计和施工,基础条件也相近。但间隔坍塌,破坏轻重悬殊(图版 24)[3]。

唐山市路北区新华西路地震破坏轻重水平分带(图版 25)[2]。

1975年海城7.3级地震也有类似的情况。

②铅直分带

一栋楼房各层的破坏常是不一致的。在唐山地震中,有些楼房只一层破坏比较严重(图版 22),有的二层破坏严重(图版10);有的只第三层破坏严重(图版 12)。

唐山华新纺织厂的一个砖砌烟囱,高约40m,断为7节,间隔大致相等,每节都有扭转,但未坠落仍在使用。

海城县委会有一烟囱,震后断成4节,最上一节右旋40o[4](图版 16)。

江苏溧阳上兴农机厂的3个烟囱互不相连,却在同一高度产生同样断裂(图版 26),这决非巧合。

图版 26.江苏溧阳上兴农机厂,3个烟囱互不相连,震后在同一高度产生同样断裂,皆向四面胀裂。

?断裂两盘不受损伤

不论断裂多长,裂缝多大,断块两盘物体仍然坚固如初,断裂所影响到的范围一般很小,例如:

①唐山至乐亭公路被错断,右旋错位 1.2m,两侧树木直立

如初(图版 27)[2]。

②海城地震时,锦县八千公社冻土开裂成一条大型地裂缝,缝宽20cm左右,长达2km将100余棵树木劈开,其中一棵树杆下部被劈,一分为二,左旋8cm,上部仍然相连,附近其他树木安然无恙(图版 28)[2]。

③穿过唐山十中的一条地裂缝将围墙、厕所和道路切断,右旋平移 1.2m。但裂缝两侧的围墙、厕所隔墙、道路及其两侧的树木皆完好无恙。地下的管道也被错开,但两端并无损伤(图版14)。

?地面平整,起伏罕见

唐山地震后,除极个别地区受采空区塌陷或其它影响出现局部起伏外,绝大部分地面和路面仍同以前一样无明显变化,波浪起伏少见,地裂缝的两侧也多是平坦的(图版 23,27,28)。

?四面开花铅直坠落

建筑体的破坏尤其是砖石结构和水泥制件的破坏一般都是分段裂开,向四面八方开花崩塌,铅直坠落。如同折叠灯笼合在一起,很少有整体倾倒的。

①唐山机车车辆厂的烟囱分段向四面胀裂(图版 29)。

②河北矿冶学院水塔断块向四面散落(图版 30)。

③唐山工人医院砖砌方烟囱四角向外崩裂。

④1979年江苏溧阳6级地震,震中附近的上沛砖瓦厂,有两个砖窑共八个烟囱,全部向四面崩散(图版 31)。

⑤江苏溧阳上兴公社农机厂,3个烟囱并排,在同一高度向四面崩裂(图版 26)。

⑥唐山开滦医院7层大楼铅直塌落,上下楼板近于重合,余下的一角仍然直立(图版 3)。

⑦唐山铁路医院住院部塌落楼板重叠在一起(图版 32)。

⑧唐山河北矿冶学院阅览室三层楼房西段全部铅直坠落(图版 22),东段仍然直立。

?地下建筑保存良好

唐山的人防坑道除个别地方有些小裂纹外,其它均未受到破坏,开滦煤矿井下坑道与地面相比破坏极轻,吕家坨在地震时井下有600多名工人在工作,震后全部安全返回地面,马家沟3号井震后10天即恢复生产出煤。位于XI度区的铁路医院住院部,原是三层楼房,地上两层全部塌平,半地下的一层基本完好,只有一些小裂缝(图版 32),唐山有两条铁路隧道也未遭破坏。

?人体抗震胜过厂房

唐山地震死伤的人,绝大部分是由于建筑物的坍塌间接受害,没有人直接死于震动,人比钢筋水泥造的厂房还要抗震。人的抗震能力是很强的。

地震破坏的这些规律也普遍存在于国内外的各次大震中,如1970年1月5日云南通海7.7级地震;1923年9月1日日本关东的8.3级地震;1960年2月29日摩洛哥阿迪尔的5.8级地震等等。

2、存在的问题

人之躯体除大脑发达外,防御自然突变的能力一般比较差。冬天怕冷,夏天怕热;风吹站不稳,水淹立不住;唯独抗震能力极强,竟胜过钢筋水泥,赛过高楼大厦,奇妙得不可思议,简单到一目了然,无需佐证。人为什么能比厂房抗震呢?这是一个非常有趣的问题,能否正确地解释这种现象,是检验一切地震理论的试金石。要解决这个问题就必须首先弄清地震的破坏作用。过去认为地震的破坏是纵波、横波和面波所为,并以后两者为主,而实际上并非如此。前述的10条地震破坏规律绝大部分是不能用纵波和横波来解释的,也不能用所谓的面波来解释。

?地震破坏的10条规律与纵波、横波及面波理论相悖。

①地震的破坏若真是纵波和横波的作用,地面物体所受的振动应当是上下颤动(纵波)和水平摇晃(横波),其破坏作用应当与大风、海浪和原子弹爆炸产生的冲击波相似,首先被破坏的和受害严重的应该是人、树木、电杆、水塔和铁塔之类的头重脚轻,重心高的,直立的线状物体;而不应该是那些重心低,结构牢,承重大、强度高的钢筋混凝土厂房,而事实正相反。由此可见,地震破坏的主力不可能是纵波和横波;

②地震断裂皆具有旋转性,而纵波和横波的振动是没有旋转性的;

③若是纵波和横波的作用,建筑物应首先倾斜,然后断裂,向一侧倒塌。而实际上是向四面坍塌近铅直坠落,这再次证明地

震的破坏不是纵波和横波所为;

④不论是横波还是纵波,它们的传播都是连续的,强度是渐变的,从震中向外逐渐衰减。在同一个地震烈度区的几十到几百平方公里的范围内,同样的建筑物所受到的破坏应当是大致相同的,而实际情况并非如此,地震破坏具分带性,在短距离内常有很大差异,一栋楼房,一半塌平,另一半完好;两条并排的铁轨,一条扭曲,另一条完好如初。这些都不能用纵波和横波作用来解释。

⑤地震破坏也不是"面波"的作用(一般认为面波起主要作用)。所谓面波按目前的解释是"体波到达地表后激发的次生波,沿地表传播,它的震动方向有的与传播方向垂直,这和横波的性质一样,包括了两种类型。一类是在地面上来回振动(拉夫波),一类则在地面上滚动(瑞雷波)。面波能使地面出现波状起伏。"既然认为面波的性质与横波一致,当然也就不可能是地震破坏的主要因素,所谓"面波"是否存在恐怕还是个问题。

地震的破坏既不是纵波和横波,也不是面波又是什么波呢?是否可能还有某种尚未发现的未知波在起作用?寻找未知波有三种途径,一是进行实地调查,研究破坏应力的性质;二是重新分析地震仪的工作原理与图象;三是从理论上探讨弹性波的成生机制,弄清弹性振动与弹性波的关系。到目前为止,震源振动与弹性传播仍处于分别讨论的状态。还没有连接成为一个整体。

?五个关键问题

存在的问题很多,主要的可归纳为以下5点:

①为什么人比厂房抗震?

②地震如何破坏?

③面波是否存在?

④有无未知波?

⑤弹性波是怎样产生的,成生机制如何?

搞清弹性波的成生机制是解决5个问题的关键,其它问题迎刃而解。

二、弹性波的成生机制

1、弹性振动

一物体受外力作用,使体积或形态发生变化,而当外力取消后,该物体便恢复它原来体积和形状者,称为弹性体。

图33之圆(实线)代表一个圆球状弹性体,若沿圆球某一直径方向施加一对外压力FF′。球体受压产生弹性变形,成为椭球体(虚线),称为变形椭球体。变形椭球体中各点的具体变化(位移)各不相同,促使各点变化的局部力大小和方向也各不相同,这些实际上使球体各点产生变形(位移)的局部力即通常说的内力。与内力相平衡还有一个反作用力,称为反弹力或恢复力。当外力和内力消除后,变形椭球体即在反弹力的作用下恢复原状,复原成球体。反弹力和内力大小相等,方向相反,作用在同一点上。单位面积上所受的内力或反弹力称为应力,可用应力单

位度量之。

图34-A所示是一个球形弹性物体的原始状态。加一对外力FF′使之发生应变(图34-B),所产生的内力同时被大小相等方向相反的反弹力所平衡。一旦外力(FF′)取消,内力则随之消失,物体将在反弹力的作用下恢复到弹性体的原始状态(A)。若外力缓慢消失,能量将逐渐消耗,物体各点则逐渐恢复到原位,直接由(B)回到(A),称为"无振复原"或"静态复原"。若外力突然取消,物体各点则在反弹力的作用下产生加速度,由(B)向(A)呈加速度运动;愈接近原始位置(A),速度愈大,达到(A)时速度最大,停不下来,于是便以此为初速度继续向前运动,给物体一个反方向内力(对B而言),使物体向(C)应变,积蓄新的反弹力,动能逐渐转化为弹性位能(应变),速度逐渐减小,当动能全部转化为弹性能时,运动停止,如34-C所示。如在运动过程中无能量损耗,应变(C)的绝对值应与原应变(B)相等,只应变方向相反。到(C)后,立刻又在新反弹力的作用下向(A)转化,演化过程和以前一样,只是方向相反,由(C)过(A)最后回到(B);到(B)后,接着再由(B)过(A)到(C)。如此往返循环不已。弹性物体的应变和反弹力的这种周期性变化称为“弹性振动”。若介质为完全弹性体又无能量消耗时将一直振动下去,永不停息。实际上由于摩擦或其它阻力的作用,能量将逐渐消耗,每个周期的最大应变和最大反弹力的绝对值将逐渐减小,最后反弹力全部消失,振动停止,物体恢复到原始位置(A)。

弹性应变的这种恢复过程称为"振动复原"。

2、弹性波

(1)弹性波的成生机制

物体的振动可由弹性介质向外传播,形成弹性波。引起弹性波的物体振动称为"震源"。弹性波的波形可由图34所示的弹性振动过程分析得出。如果我们把图34之(A),(B),(C)重合在一起绘成图35的形式,任何方向上的弹性波波形皆可直接从图上求得。

图35之圆为应变前的原始球体,即振源,弹性波沿球半径方向向外传播,物体振动时各方向上的应变不同,传播出去的弹性波也不一样。弹性波的成生遵循以下法则:振源球面上任意点的振动沿径向传播出去就是该方向上的弹性波,其波形决定于该点的振动轨迹,这个法则叫做弹性波成生法则。振源体上的振动后半个周期与前半个周期相重,演化方向相反,但传播出去就称为一个完整的周期,连续不断的波列。

沿OA方向向外传播的弹性波,其波形由a点的振动轨迹所决定,a点的振动轨迹为直a1-a-a2,与传播方向(OA)方向一致,轨迹中任意点在相应椭球上的切线(切面)都垂直于径向Oa。所以,a点的振动沿OA方向传播出去,所生成的弹性波是纵波(P波)或称压缩波,传播质点沿传播方向往返振动(图35-A)。

b点的振动沿OB方向传播出去也是纵波(图35-B),波形

与OA方向是一样的,只是在时间上相差半个周期。在同一时间内一个向外运动,一个向内运动,但由于两个纵波的始点也是相差半个周期,所以在以O为圆心的各球面上,两者的波形仍然相同,但不同步。

在OC方向上应变椭球体为无伸缩面。在整个振动过程中C 点的位置保持不变,但C点在椭球体上的瞬时切线将随椭球体的形变而变,切线与OC的交角随时而异。C点的运动轨迹是原地旋转扭摆,故沿OC方向传播出去的波是扭波。它的传播方式与钟表摆轮相仿,传播质点偏离原始位置的扭摆角度,称为扭摆角,可用C点在椭球上的瞬时切线与在原始球体上的切线的交角α来表示。C 点在原始球体上的切线与OC垂直,扭摆角α为零,C点在椭球体上的各瞬间扭摆角皆大于零。振源体的振动幅度愈大,应变量愈大,扭摆角α的变化范围也愈大,扭摆角的最大值与扭波强度成正比。因此,可用最大扭摆角度来表示扭波的强度。除a、b、c等8个特殊点外,其它所有点的振动轨迹皆是S 型曲线。例如,d点的轨迹为d1-d-d2(图35-D).其振动沿OD 方向向外传播。为了便于分析,将d点的振动时轨迹放大如图36所示:图中dk为质点d向d2振动时轨迹上的任一位置,由d 至dk的运动可分解为3种分运动。

①沿OD方向的分运动(dk),沿OD方向传播出去形成纵波(p).

②垂直OD方向的分运动(Kdk),沿OD方向传播出去形成垂

直传播方向振动的横波(S)。

③扭摆分运动沿OD传播出去形成扭波(n).

由此可见,任何振动轨迹呈曲线的原始波都可以分解成为纵波、横波和扭波等3种基本体波(图35)。3种体波同源,同时生于同地。只是由于在传播过程中传播介质对3种波的传播速度不同,才逐渐拉开,单独运行。

实际的振源体大多不是球体,在这种情况下,弹性波是沿振源表面各点的曲率半径方向向外传播,振源表面各点的振动轨迹可以全部呈S形,都包含有3种基本成分,而不存在只产生一种波形的特殊方向。

(2)弹性波的传播

纵波是一种压缩波,凡是具有压缩弹性的物体都能传递纵波。固体以及受压气体和液体都有压缩弹性,都能传递纵波。真空无压缩弹性,不能传递纵波,如声波。纵波的振动方向与传播方向一致,传播质点在一维空间内传播,在直线上振动,故传播速度最快。

横波的振动垂直传播方向,只有当介质在垂直传播方向上有压缩弹性时才能传递横波。固体各方向都有压缩弹性,无论横波从什么方向穿过固体皆可以通过。液体和气体等流体则不然,在海洋等地表大型水体中,由于重力的作用,使水滴向上移动时须用力作功,具压缩弹性,所以沿水平方向可以传播横波,尤其在水面附近重力作用最显著传播横波的能力最大,如水波。在湖海

在 ansys 中如何 施加 地震波

三向输入简化后的单向输入 首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据: -0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据: 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件,并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行,同上 ACCEXYZ(0,3)=3 !08行,同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1 acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度 SOLVE

具有不同频谱特性的地震波(精)

具有不同频谱特性的地震波 对单塔悬索桥响应的影响分析 林瑞良(福州市建设委员会 350005) [提要]根据空间有限元计算模型,采用混合结构形式,以某市单塔悬索桥为研究对 象,运用时程分析法,探讨了具有不同频谱特性的地震波对单塔悬索桥响应的影响 问题。 [关键词]单塔悬索桥时程分析地震波 现行公路桥梁工程抗震设计规范《公路工程抗震设计规范》 (JTJ-004-89)是以反 应谱理论为基础的,针对这些问题,本文以某市悬索桥为工程实例,采用动力时程分 析法,探讨了不同频谱特性的地震波对单塔悬索桥横向、纵向和竖向地震响应的影响。 一、动力计算模型的基本假设 (1) 缆索在纵向分析中取水平位移和竖向位移两个自由度,横向分析中取水平位移 一个自由度,竖向分析中取竖向位移一个自由度;(2)吊杆为柔性索,考虑变形; (3) 主塔在纵向和横向分析中均取水平位移和转动两个自由度;(4)加劲桁架在纵向分析 中取水平位移、竖向位移和转动三个自由度,横向分析中取水平位移和转动两个自由 度,竖向分析中取竖向位移和转动两个自由度;(5)作用于全桥纵向、横向上的地震 输入波,均取与基础相垂直的水平方向;作用于全桥竖直方向上的输入波取水平向输

入波的65%加速度值[1]。 二、刚度矩阵与质量矩阵 由于悬索桥结构是由不同类型的构件组成,本文在有限元计算中采用混合结构 形式的三维有限元计算模型[2],将结构划分为如下三类单元:(1)空间梁单元,用 于加劲梁及塔架。(2)空间索单元,用于主缆。(3)杆面单元,由两根吊杆和一个虚 拟刚片组成,用来反映加劲梁与主缆之间的相互作用。单元质量矩阵采用集中(堆聚) 质量矩阵[2]。将单元刚度矩阵和单元质量矩阵经座标变换,组成总刚度矩阵和总质 量矩阵,再利用子空间迭代法计算出结构的特征值和特征向量,即可得到所需的各 阶频率和振型。 三、动力方程的建立和求解 当结构在地面运动加速度X¨g作用下,结构动力方程为 [M]*{U 1}+[C]*{U 1 }+[K]*{U 1 }=-[M]+*{I}X¨g(1) 式中:[M]*和[K]*分别为缩聚后的等效质量矩阵和等效刚度矩阵; U 1 有惯性力的位移;X¨g为输入地震加速度;[C]为阻尼矩阵,按瑞雷阻尼确定。 对于微分方程式(1),可采用逐步积分的数值解法,即求得各节点的位移量,本 文采用的是威尔逊θ法,用SAP5软件进行计算。 四、具有不同频谱特性的地震波对单塔悬索桥地震响应分析实例 某市悬索桥是福建省已建成跨径最大的钢筋砼加劲桁架单塔悬索桥(见图1所示),

地震波运动学理论

第二章地震波运动学理论 一、名词解释 1. 地震波运动学:研究在地震波传播过程中的地震波波前的空间位置与其传播时间的关系,即研究波的传播规律,以及这种时空关系与地下地质构造的关系。 2. 地震波动力学:研究地震波在传播过程中波形、振幅、频率、相位等特征的及其变化规律,以及这些变化规律与地下的地层结构,岩石性质及流体性质之间存在的联系。 3. 地震波:是一种在岩层中传播的,频率较低(与天然地震的频率相近)的波,弹性波在 岩层中传播的一种通俗说法。地震波由一个震源激发。 4. 地震子波:爆炸产生的是一个延续时间很短的尖脉冲,这一尖脉冲造成破坏圈、塑性带,最后使离震源较远的介质产生弹性形变,形成地震波,地震波向外传播一定距离后,波形逐渐稳定,成为一个具有2-3个相位(极值)、延续时间60-100毫秒的地震波,称为地震子波。地震子波看作组成一道地震记录的基本元素。 5.波前:振动刚开始与静止时的分界面,即刚要开始振动的那一时刻。 6.射线:是用来描述波的传播路线的一种表示。在一定条件下,认为波及其能量是沿着一条“路径”从波源传到所观测的一点P。这是一条假想的路径,也叫波线。射线总是与波阵面垂直,波动经过每一点都可以设想有这么一条波线。 7. 振动图和波剖面:某点振动随时间的变化的曲线称为振动曲线,也称振动图。地震勘探中,沿测线画出的波形曲线,也称波剖面。 8. 折射波:当入射波大于临界角时,出现滑行波和全反射。在分界面上的滑行波有另一种特性,即会影响第一界面,并激发新的波。在地震勘探中,由滑行波引起的波叫折射波,也叫做首波。入射波以临界角或大于临界角入射高速介质所产生的波 9.滑行波:由透射定律可知,如果V2>V1 ,即sinθ2 > sinθ1 ,θ2 > θ1。当θ1还没到90o时,θ2 到达90o,此时透射波在第二种介质中沿界面滑行,产生的波为滑行波。 10.同相轴和等相位面:同向轴是一组地震道上整齐排列的相位,表示一个新的地震波的到达,由地震记录上系统的相位或振幅变化表示。 11.地震视速度:当波的传播方向与观测方向不一致(夹角θ)时,观测到的速度并不是波前的真速度V,而是视速度Va。即波沿测线方向传播速度。 12 波阻抗:指的是介质(地层)的密度和波的速度的乘积(Zi=ρiVi,i为地层),在声学中称为声阻抗,在地震学中称波阻抗。波的反射和透射与分界面两边介质的波阻抗有关。只有在Z1≠Z2的条件下,地震波才会发生反射,差别越大,反射也越强。 13.纵波:质点振动方向与波的传播方向一致,传播速度最快。又称压缩波、膨胀波、纵波或P-波。 14.横波:质点振动方向与波的传播方向垂直,速度比纵波慢,也称剪切波、旋转波、横波或S-波,速度小于纵波约0.7倍。横波分为SV和SH波两种形式。 15.体波:波在无穷大均匀介质(固体)中传播时有两种类型的波(纵波和横波),它们在介质的整个立体空间中传播,合称体波。 16共炮点反射道集:在同一炮点激发,不同接收点上接收的反射波记录,称为共炮点道集。在野外的数据采集原始记录中,常以这种记录形式。可分单边放炮和中间放炮。 17.面波:波在自由表面或岩体分界面上传播的一种类型的波。 18.纵测线和非纵测线:激发点与接收点在同一条直线上,这样的测线称为纵测线。用纵测线进行观测得到的时距曲线称为纵时距曲线。激发点不在测线上,用非纵测线进行观测得到的时距曲线称为非纵时距曲线。

时程分析中地震波输入位置的讨论

时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系,能真实地反应结构地震相应随时间变化的全过程,是抗震分析的一种重要方法[1]。目前有限元软件可以实现结构的时程分析,但是在不同的软件中,其实现方式不同,主要区别在地震波的输入位置不同。本文通过有限元软件ABAQUS采用不同的地震波输入位置对同一结构进行时程分析分析,对比结构相同位置的时程位移曲线,结果表明结构在采用不同地震波输入位置的时程分析中,结构的地震响应基本一致。 关键词:时程分析、有限元软件、钢筋混凝土剪力墙 Abstract: The time history analysis method to analyze the available structure through direct power to the relationship between the corresponding changes over time, truly reflect the structure of earthquake corresponding to the whole process of change over time, is an important method of seismic analysis [1]. Finite element software can be time-history analysis of the structure, but in different software in different ways, the main difference between the different positions in the seismic wave input. In this paper the finite element software ABAQUS using different seismic wave input location on the same structure, process analysis analysis, contrast structure the same location of when the process displacement curve, the results show that the structure using different seismic waves enter the position time history analysis, the seismic response basically the same. Keywords: time history analysis, finite element software, reinforced concrete shear walls 一、引言 在时程分析等动力学问题中,地震力以加速度形式从基础固定处输入。由于结构的刚度不是无限大,在结构上的加速度反应与基础输入的加速度并不相同。在很多时候,结构的加速度比基础输入的加速度更大,即对输入的加速度有一个动力放大效应。在单自由度弹性体系中,体系最大绝对加速度与地面运动最大加速度的比值,即称为动力系数[2] (1) 动力系数与结构的动力学特性和输入的地震波的频率特性有关。它与地震系数k的乘积即为单自由度体系的地震影响系数。 因此,从原理上讲,时程分析是将地震波的加速度时程曲线作用到结构的基础约束处,得到上部结构的各种地震反应。但是在不同的软件中,其实现方

地震波使用说明

地震波使用说明 此目录下提供了四类场地土的地震波时程曲线和上海人工波。 按照场地土类型(1,2,3或4),选择时程曲线。在定义时程工况时,对于多遇或罕遇地震,按比例调整时程曲线的最大值。中国抗震规范规定,作为抗震计算中底部剪力法和振型分解反应谱法的补充方法,对于特别不规则,特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时,应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线,其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时,每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65%,多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80% 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析,可以是线性或非线性的。 此章对时程分析进行一般的描述,特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”,导入已定义的文本文件,即实测的时程曲线;也可使用程序内置的时程函数。

时程函数定义对话框 时程函数定义对话框中的条目解释如下: ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框,在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。

论地震勘探中几种主要地震波

论地震勘探中的几种主要地震波 论文提要 地震勘探,就是通过人工方法激发地震波,研究地震波在地层中传播的情况,以查明地下地质构造,为寻找油气田或其它勘探目的服务的一种方法。也可以理解为就是利用地震子波从地下地层界面反射回地面时带回来的旅行时间和形状变化的信息,用以推断地下的底层构造和岩性。地震勘探在勘探已有的各种物探方法中,是最有效地方法。在地震勘探中用炸药激发时,一声炮响之后会产生各种各样的地震波。按波在传播过程中质点震动的方向来区分,可以纵波和横波;根据波动所能传播的空间范围而言,地震波又可以分为体波和面波;按照波在传播过程中的传播路径的特点,又可以把地震波分为直达波、反射波、透射波、折射波,等等。地震勘探在石油勘探中除了能产生来自地层界面有用的反射波外,还会产生各种各样的干扰波。因此,我们要更好的了解各种波的产生、特点、用途,等等。下面简单介绍几种地震勘探中产生的地震波。 正文 一、反射波 (一)反射波的形成 1、几何地震学的观点 当炸药在井中爆炸激发地震波时,在雷管引爆几百微妙之内爆炸便完成了,在接近爆炸点的压强是一个延续时间很短的尖脉冲,爆炸脉冲向外传播,压强逐渐减少,地层开始产生弹性形变,形成地震波。地震波继续传播,由于介质对高频的吸收,地震波信号减小。当波入射到两种介质的分界面时(当上层介质波阻抗与下层介质波阻抗不等时,弹性地震波才会发生反射;上层介质波阻抗与下层介质波阻抗差别越大,反射波越强——反射波条件),一部分波回到第一种介质中,这就是所谓的反射波。如图所示 2、物理地震学观点 地震波从震源出发以球面波的方式向下传播,到达反射界面S,S可以就看成有许多

abaqus如何施加地震波

施加地震波: 1 *amplitude,name=amp,input=seismicdata.dat 输入地震波 2 *boundary,type=acceleration,amplitude=amp施加荷载 方法:module选load,在tools-----amplitude-----creat默认的continue在Edit A mplitude里面输入时间和加速度,点OK。点creat boundary condition,涌现对 话框creat boundary condition,选择acceleration/angular acceleration,continu e---选择要施加的边界---done----涌现对话框edit bondary condition对话框,在 amplitude里选择你所定义的时间和加速度。点ok就完工了。 在网上查了些方法: module选load,在tools-----amplitude-----creat默认的continue在Edit Amplitude 里面输入时间和加速度,点OK。点creat boundary condition,出现对话框creat boundary condition,选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框,在amplitude里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式,我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型,定义两个分析步。 第一个分析步是加自重,采用线性加载的方式。 (a) 加载方式:ABAQUS在施加Gravity时,默认为Instantaneous(瞬时加载),如果把结构自重以瞬间加载方式加到结构上,相当于对结构施加了一个脉冲荷载,会引起结构在竖向的振动,在不考虑结构阻尼的情况,这种振动会一直持续下去。如果是混凝土结构,这种竖向振动也会造成混凝土受拉损伤,所以这种加载方式不太合理。 (b)新建加载方式:创建一个新的Amplitude,Type=smooth tpye,0时刻Am=0,然后再选择一个0.5s~1s时刻,Am=1,在这个区间内线性插值,实现幅值从0到1。这种方式加载要优于上述瞬时加载,但是在起初的0.5s(或者1s,即smooth tpye中设置的终点时间)内计算结果是不准确的,所以要把这部分的计算结果剔除,剔除方法就是,创建2个step,第一个step主要分析自重作用,待自重稳定后开始第二个step地震时程反应分析。 第二个分析步就是加地震波。 输入地震波有两种方法: 1、在如下位置加入下面加黑的字体部分。格式如下:时间,地震波,时间,地震波,时间, 地震波,时间, 地震波…………每行8个数据(我下到的地震波文件是不带时间的,自己用C++处理了一下)。%%%%%%%%%%%%%%%%%%%%%% *End Assembly *Amplitude, name=Amp-1 0.005, -7.5e-08, 0.01, -3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-07 0.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-05 0.045, 3.54e-05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-05 ………………

1地震波

地震波 地震波 地震波是指从震源产生向四外辐射的弹性波。 分类 地震波按传播方式分为三种类型:纵波、横波和面波。 纵波、横波 纵波,又称P波,是推进波,地壳中传播速度为5.5~7千米/秒,最先到达震中,它使地面发生上下振动,破坏性较弱。 横波,又称S波,是剪切波:在地壳中的传播速度为3.2~4.0千米/秒,第二个到达震中,它使地面发生前后、左右抖动,破坏性较强。 面波又称L波,是由纵波与横波在地表相遇后激发产生的混合波。其波长大、振幅强,只能沿地表面传播,是造成建筑物强烈破坏的主要因素。 P波和S波的实际传播速度取决于岩石的密度和内在的弹性。对线弹性物质而言,当波与运行方向无关时,波速仅取决于两个弹性性质,称为弹性模量:岩石的体积模量E和剪切模量G。

面波 当P波和S波到达地球的自由面或位于层状地质构造的界面时,在一定条件下会产生其他类型地震波。这些波中最重要的是瑞利波和勒夫波。这两类波沿地球表面传播;岩石振动振幅随深度增加而逐渐减小至零。 勒夫波 勒夫波又称Q波。是地震面波中最简单的一种类型。它们是以1912年首次描述它们的勒夫的姓名命名的。是一种表面波通过切变波在表层内的多次内反射而传播。在半无限介质之上出现低速层的情况下,一种垂直于传播方向的在水平面内振动的波。 这个类型的波使岩石质点运动类似SH波,运动没有垂向位移。岩石运动在一垂直于传播方向上在水平面内从一边到另一边。虽然勒夫波不包括垂直地面运动的波,但它们在地震中可以成为最具破坏性的,因为它们常具有很大振幅,能在建筑物地基之下造成水平剪切。瑞利波 瑞利波Rayleighwave,瑞利面波具有相当不同的地面运动。 地震学中称其为R波或L波。于1885年首次由瑞利(LordRayleigh)描述,它们是地震波中最近似水波的。岩石质点向前、向上、向后和向下运动,沿波的传播方向作一垂直平面,质点在该平面内运动,描绘出一个椭圆。勒夫波和瑞利波的速度总比P波小,与S波的速度相

ANSYS地震波的输入

对于地震波的输入,可以把荷载记录作成文件,利用apdl的读取功能读入数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。 Fini /config,nres,1000 *dim,aceX,TABLE,3000,1 *dim,aceY,TABLE,3000,1 *dim,aceZ,TABLE,3000,1 *creat,ff *vread,aceX(1,1),acex,txt,,1 (e16.6) *vread,aceX(1,0),acexTT,txt,,,1 (e16.6) ACEX(0,1)=1 *end /input,ff *creat,ff *vread,aceY(1,1),txt,,1 (e16.6) *vread,aceY(1,1),ACETT,,,1 (e17.6) ACEY(0,1)=1 *end /input,ff *creat,ff *vread,aceZ(1,1),txt,,1 (e16.6) *vread,aceZ(1,0),ACETT,,,1 (e17.6) ACEZ(0,1)=1 *end /input,ff !地震波时程记录分成了3个文件,每个文件是一列。分别记录x,y,z方向的加速度。Accett是时间记录。 这样就可以把加速度记录读取到ansys数据库中作为数组。 也可以把加速度记录作成一个文件,这样程序就简单多了。 下面是计算部分语句: /SOLU ANTYPE,trans !求解其自己选了 TM_START=0.01 TM_END=15.00 TM_INCR=0.01 *DO,TM,TM_START,TM_END,TM_INCR TIME,tm

地震波的定义

地震波的定义

地震波的定义 地震是地壳的一切颤动,是一种自然现象。其主要能源来自地球的内部,是由地球内部自然力冲击引起的。地壳或地幔中发生振动的地方称为震源。震源在地面上的垂直投影称为震中。震中到震源的距离称为震源深度。地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面,将地球内部分为地壳、地幔和地核三个圈层。 发生原理 英文seismic wave.由地震震源发出的在地球介质中传播的弹性波。地球内 地震波 部存在着地震波速度突变的基干界面、莫霍面和古登堡面,将地球内部分为地壳、地幔和地核三个圈层。地震震源发出的在地球介质中传播的弹性波。地震发生时,震源区的介质发生急速的破裂和运动,这种扰动构成一个波源。由于地球介质的连续性,这种波动就向地球内部及表层各处传播开去,形成了连续介质中的弹性波。 概念介绍 地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面,将地球内部分为地壳、地幔和地核三个圈层。 传播方式 地震波按传播方式分为三种类型:纵波、横波和面波[1]。纵波是推进波,地壳中传播速度为5.5~7千米/秒,最先到达震中,又称P波,它使地面发生上下振动,破坏性较弱。横波是剪切波:在地壳中的传播速度为3.2~4.0千米/秒,第二个到达震中,又称S波,它使地面发生前后、左右抖动,破坏性较强。面波又称L波,是由纵波与横波在地表相遇后激发产生的混合波。其波长大、振幅强,只能沿地表面传播,是造成建筑物强烈破坏的主要因素。 纵波和横波

现象介绍 我们最熟悉的波动是观察到的水波。当向池塘里扔一块石头时水面被扰乱,以石头入水处为中心有波纹向外扩展。这个波列是水波附近的水的颗粒运动造成的。然而水并没有朝着水波传播的方向流;如果水面浮着一个软木塞,它将上下跳动,但并不会从原来位置移走。这个扰动由水粒的简单前后运动连续地传下去,从一个颗粒把运动传给更前面的颗粒。这样,水波携带石击打破的水面的能量向池边运移并在岸边激起浪花。地震运动与此相当类似。我们感受到的摇动就是由地震波的能量产生的弹性 岩石的震动。 假设一弹性体,如岩石,受到打击,会产生两类弹性波从源向外传播。第一类波的物理特性恰如声波。声波,乃至超声波,都是在空气里由交替的挤压(推)和扩张(拉)而传递。因为液体、气体和固体岩石一样能够被压缩,同样类型的波能在水体如海洋和湖泊及固体地球中穿过。在地震时,这种类型的波从断裂处以同等速度向所有方向外传,交替地挤压和拉张它们穿过的岩石,其颗粒在这些波传播的方向上向前和向后运动,换句话说,这些颗粒的运动是垂直于波前的。向前和向后的位移量称为振幅。在地震学中,这种类型的波叫P波,即纵波(图2.1),它是首先到达的波。 地震P波(纵波)和S波(横波)运行时弹性岩石运动的形态 弹性岩石与空气有所不同,空气可受压缩但不能剪切,而弹性物质通过使物体剪切和扭动,可以允许第二类波传播。地震产生这种第二个到达的波叫S 波,即横波。在S波通过时,岩石的表现与在P波传播过程中的表现相当不同。因为S波涉及剪切而不是挤压,使岩石颗粒的运动横过运移方向(图2.1)。这些岩石运动可在一垂直向或水平面里,它们与光波的横向运动相似。P和S波同时存在使地震波列成为具有独特的性质组合,使之不同于光波或声波的物理表现。因为液体或气体内不可能发生剪切运动,S波不能在它们中传播。P和S波这种截然不同的性质可被用来探测地球深部流体带的存在(见第6章)。 相关性质 带偏光眼镜以减弱散射光的人可能熟悉光的偏振现象,只有S波具有偏振现象。只有那些在某个特定平面里横向振动(上下、水平等)的那些光波能穿过偏光透镜。传过的光波称之为平面偏振光。太阳光穿过大气是没有偏振的,即没有光波振动的优选的横方向。然而晶体的折射或通过特殊制造的塑料如偏光眼镜,可使非偏振光成为平面偏振光。 当S波穿过地球时,他们遇到构造不连续界面时会发生折射或反射,并使其振动方向发生偏振。当发生偏振的S波的岩石颗粒仅在水平面中运动时,称为SH波。当岩石颗粒在包含波传播方向的垂直平面里运动时,这种S波称为SV

Midas地震波的选取方法

地震波的选取方法 建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定,正确选择输入的地震加速度时程曲线,要满足地震动三要素的要求,即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征,依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同,简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时T d的定义可分为两大类,一类是以地震动幅值的绝对值来定义的绝对持时,即指地震地面加速度值大于某值的时间总和,即绝对值|a(t)|>k*g的时间总和,k常取为0.05;另一类为以相对值定义的相对持时,即最先与最后一个k*a max之间的时段长度,k一般取0.3~0.5。不论实际的强震记录还是人工模拟波形,一般持续时间取结构基本周期的5~10倍。 说明: 有效峰值加速度 EPA=Sa/2.5 (1) 有效峰值速度 EPV=Sv/2.5 (2) 特征周期 Tg = 2π*EPV/EPA(3) 1978年美国ATC-3规范中将阻尼比为5%的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa,将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱),上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。 上述方法使用的是将频段固定的方法来求EPA和EPV,1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是:在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱,找出加速度反应谱平台段的起始周期T0和结束周期T1,然后在拟速度反应谱上选定平台段,其起始周期为T1(即加速度反应谱平台段的结束周期T1),结束周期为T2,将加速度反应谱在T0至T1之间的谱值求平均得Sa,拟速度反应谱在T1至T2之间的谱值求平均得Sv(注:生成谱的时候一定要用对数谱),加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5,按公式(1)、(2)、(3)求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震波数据生成器,生成后保存为SGS文件),用户可利用保存的SGS文件(文本格式文件)根据上面所述方法计算Sv、Sa、Tg=Sv/Sa。通过Tg值可判断该地震波是否适合当地场地和地震设计分组,然后将抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数(即放大系数),将其代入到地震波调整系数中。将地震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。 建筑抗震设计规范5.1.2条中规定,采用时程分析方法时,应按照场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线,其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。所谓“在统计意义上相符”指的是,其平均影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比,在各周期点上相差不大于20%。 在MIDAS程序中,可选取两组实际强震记录生成两个SGS文件(调整Sa后的),然后将一组人

结构抗震设计时程分地震波的选择

(1)设计用地震记录的选择和调整 用规范的确定性方法和地震危险性分析方法所确定的设计地震动参数,是选择天然地震加速度记录的依据。 (一)实际地震记录的选择方法 选择地震记录应考虑地震动三要素,即强度(峰值)、频谱和持续时间。对某一建筑的抗震设计,最好是选用该建筑所在场地曾经记录 到的地震加速度时间过程。但是,这种机会极少。为此,人们只能从现有的国内外常用的地震记录中去选择,尽可能挑选那些在震级、震中距和场地条件等方面都比较接近设计地震动参数的记录。他的文章给出了相应的地震数据的记录目录。 (二)实际地震记录的调整 1.强度调整。将地震记录的加速度值按适当的比例放大或缩小,使其峰值加速度等于事先所确定的设计地震加速度峰值。即令 其中a(为记录的加速度值为调整后的加速度值;A众为设计地震加速度峰值;。为记录的加速度峰值。这种调整只是针对原记录的强度进行的,基本上保留了实际地震记录的特征。也就是所说的(强度修正。将地震波的加速度峰值及所有的离散点都按比例放大或缩小以满足场地的烈度要求)

2.频率调整考虑到场地条件对地震地面运动的影响,原则上所选择的实际地震记录的富氏谱或功率谱的卓越周期乃至形状,应尽量与场地土相应的谱的特性一致。如果不一致,可以调整实际地震记录的时间步长,即将记录的时间轴“拉长”或“缩短”,以改变其卓越周期而加速度值不变也可以用数字滤波的方法滤去某些频率成分,改变谱的形状。另外,为了在计算中得到结构的最大反应,也可以根据建筑结构基本自振周期,调整实际地震记录的卓越周期,使二者接近。这种调整的结果,改变了实际地震记录的频率结构,从物理意义上分析是不合理的。 另外,在测定场地土和建筑结构的卓越周期时,运用不同的测试仪器和测试技术,往往得到不同的结果。即使是对同一个测试结果,在频谱上确定卓越周期时,不同的分析方法也会导致不同的结果。有的选取谱的第一个峰值所对应的周期作为卓越周期,有的选最大峰值时的,也有的取某一段周期等,很不一致。对如何确定地震加速度记录的卓越周期,也是各行其是,有的利用加速度反应谱,有的用伪速度谱,有的用富氏谱,结果当然是不一样的。上述各种作法在工程中引起了一些混乱。 王亚勇认为,用脉动测试方法测定场地土和结构的卓越周期及自振周期时,应采用速度摆型或加速度摆型的地震仪测定地运动和结构振动,然后计算其富氏谱或功率谱,以谱的最大峰值所对应的周期作为卓越周期和自振周期比较合适。反应而相应地根据记录的位移谱或速度谱。 这也就是所谓的滤波修正。可按要求设计滤波器,对地震波进行时域或频域的滤波修正。这样修正的地震资料不仅卓越周期满足要求,功率谱的形状和面积也可控制。卓越周期修正。将地震波的离散步长按人为比例改变,

地震波传播原理

菲涅尔体和透射波 摘要 在地震成像实验中,通常使用基于波动方程高频渐进解的几何射线理论,因此,通常假设地震波沿着空间中一条连接激发点和接受点的无限窄的线传播,称为射线。事实上,地震记录有非常多的频率成分。地震波频率的带限性就表明波的传播应该扩展到几何射线周围的有限空间。这一空间范围就成为菲涅尔体。在这片教案中,我们讲介绍关于菲涅尔体的物理理论,展示适用于带限地震波的波动方程的解。波动方程的有限频理论通过敏感核函数精确地描述了带限透射波和反射波的旅行时与振幅和地球介质中慢度扰动之间的线性关系。菲涅尔体和有限频敏感核函数可以通过地震波相长干涉的概念联系起来。波动方程的有限频理论引出了一个反直觉的结论-在三维几何射线上的点状速度扰动不会不会造成波长的相位扰动。因此,这说明在射线理论下的菲涅尔体理论是波动方程有限频理论在有限频下的一个特例。最后,我们还澄清了关于菲涅尔体宽度限制成像实验分辨率的误解。 引言 在地震成像技术中,射线理论通常在正演和反演中被用有构建正反演波长算子。射线理论之所以收到欢迎部分是由于计算机速度和内存的限制,因为射线理论具有较高的计算效率并且对于各种地震成像方法的应用也比较容易。而另一方面,地震成像实验清晰的表明,射线理论,由于他对波场传播的近似描述,对于散射效应严重的波场的成像是不完备的。Cerveny 给出了对于地震波射线理论的一个全面的理解。 在地震成像实验中,记录到的透射波和反射波信号都是由一个主要由低频信号组成的宽带震源激发产生的,因为地震波的高频信号在地层中很容易衰减。但是射线理论是基于高频近似的,这表明基于射线理论的成像技术和和测量波场这件之能会存在方法上的冲突。这个围绕射线且对带限地震波的传播起主要影响的空间范围就被叫做菲涅尔体。射线理论在地下构造尺度大于记录波场的第一菲涅尔带的介质中能够取得较好的效果。对于低频反射波(频率成分在10-70Hz之间)和透射波(频率成分在300-800Hz之间),第一菲涅尔体的宽度可以分别达到500m和50m的量级。这个宽度要大于我们在陆地和海洋的反射波地震勘探以及井间和垂直地震剖面中想要成像的地下地质特征。 在这篇教案中,我们将看到如何将地震分辨率扩展到识别体积小于第一菲涅尔带的不均匀体。我们将展示如把射线理论下的旅行时和振幅公式扩展到更精确的、可以应用与带限反射和透射地震信号波场近似理论。波动方程的有限频理论提出了反射和透射地震波的敏感核函数(也称作Frechet核函数)。这些有限频Frechet核函数将速度扰动和旅行时与振幅的扰动线性的联系起来。有限频波长近似被直接应用到各种地震成

新输入地震波

[结构分析] 地震波输入的问题 三向输入简化后的单向输入 首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据:-0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据: 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件,并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行 (3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行,同上 ACCEXYZ(0,3)=3 !08行,同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1

《地震波理论》复习最终版

《地震波理论》复习内容 一、弹性理论基础 1. 柯西公式的意义; 因此弹性体内一点的应力状态可以完全由作用于垂直坐标轴方向的三个截面上的应力向量或其分量所确定。 2. 应力与应变的关系; (为单位函数) 3. 杨氏模量E(纵向应力与纵向应变的比例常数就是材料的弹性模量E, 也叫杨氏模量) 泊松比ν(横向应变与纵向应变之比值称为泊松比,也叫横向变性系数,它是反映材料横向变形的弹性常数); 4. 拉梅常数λ、μ; 为引入均匀各向同性介质中应力与应变关系,引入λ、μ,μ表示剪切模量。 5. 运动的应力方程和位移方程; 运动应力方程:

运动位移方程: 6. 介质受应力作用产生位移由哪几部分组成; 由式上式可以看出处于应力应变状态下的物体其质点位移由三部分组成: ①平动: u,v,w,这是和参考点M一起作同样的运动,它不使物体形 状改变; ②弹性应变: eij,i,j=x,y,z 这是一种使物体形状和体积发生改变的 运动,称为弹性应变.应变有九个分量,考虑到它的对称性,只有其中六个分量独立的。exx,eyy,ezz称为正应变,exy,eyz,ezx称为切应变; ③旋转: ωx,ωy,ωz这是质点围绕参考点M的旋转运动,不使物体形 状和体积发生改变,不属弹性应变范畴. 7. 导出拉梅方程的前提条件; 在对空间求导时,只有λ、μ不随空间变化,即在均匀介质中才能导出拉梅方程。 8. 能流密度。 表示在单位时间内通过与它垂直的单位截面积的机械能。 二、弹性动力学中的基本波 1. 由拉梅方程导出纵波、横波方程;

拉梅方程 对上式进行散度运算,得到: 对上式进行旋度运算,得到: 2. 平面波、不均匀平面波; 平面波:等相位为平面,且与波的传播方向垂直的波动。 不均匀平面波:平面波传播的方向余弦为l 、m 、n 是复数,这样的波为 不均匀平面波。 3. 在什么情况下才能称为平面波; 离震源较远时可以将在局部等相位内,将点震源产生的球面波看成一个 平面。 4. 沿着x 方向传播的平面波的表示方法; 5. xoz 平面内的波剖面(图2-3);P 39 若该方程波函数为: ,则波剖面如下图所示 1(,,)()lx nz C f x z t f t + =-''''exp()(,,)exp()exp()Ksh z Kc f x z t A j x j t h ωθθ=--12(x,)()()x x C C f t f t f t =-+ +''' j θθθ=+

地震波

地震波 地震被按传播方式分为三种类型:纵波、横波和面波[1]。纵波是推进波,地壳中传播速度为5.5~7千米/秒,最先到达震中,又称P波,它使地面发生上下振动,破坏性较弱。横波是剪切波:在地壳中的传播速度为3.2~4.0千米/秒,第二个到达震中,又称S波,它使地面发生前后、左右抖动,破坏性较强。面波又称L波,是由纵波与横波在地表相遇后激发产生的混合波。其波长大、振幅强,只能沿地表面传播,是造成建筑物强烈破坏的主要因素。 [编辑本段] 地震纵波和横波 我们最熟悉的波动是观察到的水波。当向池塘里扔一块石头时水面被扰乱,以石头入水处为中心有波纹向外扩展。这个波列是水波附近的水的颗粒运动造成的。然而水并没有朝着水波传播的方向流;如果水面浮着一个软木塞,它将上下跳动,但并不会从原来位置移走。这个扰动由水粒的简单前后运动连续地传下去,从一个颗粒把运动传给更前面的颗粒。这样,水波携带石击打破的水面的能量向池边运移并在岸边激起浪花。地震运动与此相当类似。我们感受到的摇动就是由地震波的能量产生的弹性 地震波 岩石的震动。 假设一弹性体,如岩石,受到打击,会产生两类弹性波从源向外传播。第一类波的物理特性恰如声波。声波,乃至超声波,都是在空气里由交替的挤压(推)和扩张(拉)而传递。因为液体、气体和固体岩石一样能够被压缩,同样类型的波能在水体如海洋和湖泊及固体地球中穿过。在地震时,这种类型的波从断裂处以同等速度向所有方向外传,交替地挤压和拉张它们穿过的岩石,其颗粒在这些波传播的方向上向前和向后运动,换句话说,这些颗粒的运动是垂直于波前的。向前和向后的位移量称为振幅。在地震学中,这种类型的波叫P波,即纵波(图2.1),它是首先到达的波。 图2.1 地震P波(纵波)和S波(横波)运行时弹性岩石运动的形态 弹性岩石与空气有所不同,空气可受压缩但不能剪切,而弹性物质通过使物体剪切和扭动,可以允许第二类波传播。地震产生这种第二个到达的波叫S波,即横波。

时程分析时地震波的选取及地震波的反应谱化

时程分析时地震波的选取及地震波的反应谱化 发表时间:2017-12-29T15:40:37.810Z 来源:《防护工程》2017年第22期作者:金林飞 [导读] 目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。 常州市规划设计院江苏常州 213002 摘要:目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。时程分析法作为补充计算方法,在不规则、重要或较高建筑中采用。进行时程分析时,首先面临正确选择输入的地震加速度时程曲线的问题。时程曲线的选择是否满足规范的要求,则需要首先将时程曲线进行单自由度反应计算,得到其反应谱曲线,并按规范要求和规范反应谱进行对比和取舍。本文通过介绍常用的数值计算方法及计算步骤,实现将地震加速度时程曲线计算转化成反应谱曲线,从而为特定工程在时程分析时地震波的选取提供帮助。 关键词:时程分析,地震波,反应谱,动力计算 1 地震反应分析方法的发展过程 结构的地震反应取决于地震动和结构特性。因此,地震反应分析的水平也是随着人们对这两个方面认识的深入而提高的。结构地震反应分析的发展可以分为静力法、反应谱法、动力分析法这三个阶段。在动力分析法阶段中又可分为弹性和非弹性(或非线性)两个阶段。 [1] 目前,在我国和其他许多国家的抗震设计规范中,广泛采用反应谱法确定地震作用,其中以加速度反应谱应用得最多。反应谱是指:单自由度弹性体系在给定的地震作用下,某个最大反应量(如加速度、速度、位移等)与体系自振周期的关系曲线。反应谱理论是指:结构物可以简化为多自由度体系,多自由度体系的地震反应可以按振型分解为多个单自由度体系反应的组合,每个单自由度体系的最大反应可以从反应谱求得。其优点是物理概念清晰,计算方法较为简单,参数易于确定。 反应谱理论包括如下三个基本假定:1、结构物的地震反应是弹性的,可以采用叠加原理来进行振型组合;2、现有反应谱假定结构的所有支座处地震动完全相同;3、结构物最不利的地震反应为其最大地震反应,而与其他动力反应参数,如最大值附近的次数、概率、持时等无关。[1] 时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由于此法是对运动方程直接求解,又称直接动力分析法。可直接计算地震期间结构的位移、速度和加速度时程反应,从而描述结构在强地震作用下弹性和非弹性阶段的内力变化,以及结构构件逐步开裂、屈服、破坏甚至倒塌全过程。 根据我国《建筑抗震设计规范》(GB5011-2010)(以下简称《抗规》)第5.1.2-3条要求,特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑,应采用时程分析法进行多遇地震下的补充计算。此外《高层建筑混凝土结构技术规程》(JGJ3-2010) (以下简称《高规》)第4.3.4条也有相关要求。 2 时程分析时地震波的选取要求 在进行时程分析时,首先面临地震波选取的问题。所选的地震波需要符合场地条件、设防类别、震中距远近等因素。《抗规》对于地震波的选取主要有以下几点要求: 1、当取三组加速度时程曲线输入时,计算结果宜取时程法的包络值和振型分解反应谱法的较大值;当取七组及七组以上的时程曲线时,计算结果可取时程法的平均值和振型分解反应谱法的较大值(其中实际强震记录的数量不应少于总数的2/3)。 2、弹性时程分析时,每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65%,多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80%。 3、多组时程曲线的平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。根据规范条文说明,所谓“统计意义上相符”指的是,多组时程波的平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比,在对应于结构主要振型的周期点上相差不大于20%。但计算结果也不能太大,每条地震波输入计算不大于135%,平均不大于120%。 4、时程曲线要满足地震动三要素的要求,即频谱特性、有效峰值和持续时间均要符合规定。其中频谱特性可用地震影响系数曲线表征,依据所处的场地类别和设计地震分组确定;加速度的有效峰值按《抗规》表5.1.2-2中所列地震加速度最大值采用;输入的地震加速度时程曲线的有效持续时间,一般从首次达到该时程曲线最大峰值的10%那一点算起,到最后一点达到最大峰值的10%为止。有效持续时间一般为结构的基本周期的(5~10)倍,即结构顶点的位移可按基本周期往复(5~10)次。 根据以上《抗规》对地震波选取的要求,在工程中选取地震波时,首先需将地震波进行单自由度动力求解,得出相应的反应谱曲线,并与规范反应谱所用的地震影响系数曲线进行对比。对比频谱特性,查看是否在统计意义上相符,若相符则可将此地震波用于结构计算进行底部剪力验证,否则则需要更换地震波重新分析。 3 求解地震波单自由度体系反应的数值计算方法 对一般动力荷载反应的求解方法可分为叠加法和逐步法两类。 叠加法包括时域分析和频域分析,总反应计算采用独立反应贡献的组合。时域方法中,荷载被考虑为短暂持续时间的脉冲序列,由每个脉冲自由振动反应的独立贡献得到后续时间的总反应;频遇方法中,假设荷载为周期的,并用Fourier变换为离散的谐振分量Pn。再由这些荷载分量乘以结构的频率反应系数Hn,得到与其相应的结构谐振反应分量Vn。最后,由组合谐振反应分量(Fourier逆变换)获得结构的总反应。由于结果都使用了叠加,因此不适用于非线性反应分析。[2] 逐步法有很多种,但所有方法都是将荷载和反应历程分成一系列时间间隔或“步”。在每步期间均以此步开始时存在的初始条件(位移、速度和加速度)和该步期间的荷载历程来计算反应。因此每步反应是一个独立的分析问题。[2] 积分法是逐步法中的一般性的方法之一。积分法对每一时间步,从初始最终条件应用积分向前进一步,速度的变化依赖于加速度历程的积分,而位移的变化依赖于相应的速度积分。 用积分法进行分析时,首先需要假设在时间步内加速度是如何变化的。根据假设的不同,积分法可分为“基于常平均加速度”法(如

相关文档
最新文档