爆破地震波特性研究

爆破地震波特性研究
爆破地震波特性研究

爆破地震波特性研究3

张义平,吴桂义

(贵州大学矿业学院, 贵州贵阳 550003)

摘 要:结合现场爆破震动信号,从爆破地震波的传播形式、传播方式、波的特征、波的衰减吸收及传播介质的力学模型等方面分析了爆破地震波特性。结果表明:爆破地震波是一种与自然地震波相似但又相区别的非常复杂的随机过程,它是不同幅值、不同频率与不同相位的各种波型叠加而成的复合波。爆破地震波在传播过程中会发生多次反射、折射、绕射、衍射、波型转换甚至波导、层间波等复杂现象,传播过程中波的有关参数和时频特征常与爆源条件、传播介质的物理性质、场地特征及地形等因素紧密相关。地震波在发生几何衰减的同时,还因粘弹性介质的内摩擦和热传导导致能量耗散,使得波能不断衰减。

关键词:爆破地震波;波特性;衰减与吸收

中图分类号:T D235.1 文献标识码:A

文章编号:1005-2763(2007)06-0068-05

Study on Character isti cs of Bl a sti n g-Caused Se is m i c W ave

Zhang Y iping,W u Guiyi

(College ofM ining,Guizhou University,Guiyang,

Guizhou550003,China)

Abstract:Combined with the data collected fr om the in-site monit oring of blasting vibrati on,the characteristics of blasting-caused seis m ic wave are analyzed comp rehensively fr om its p r op2 agati on for m,p r opagati on mode,p r operties,attenuati on,ab2 s or p ti on and the mechanics model of p r opagati on mediu m s.The results show that blasting-caused seis m ic wave,which is a very comp lex random p r ocess rese mbling t o be diffence fr om seis m ic wave,is a composite wave composed of kinds of waves with dif2 ferent ranges,frequencies and phases.I n the p r opagati on p r ocesses of blasting-caused seis m ic wave,comp lex phenome2 na such as many ti m es of wave reflecti on,refracti on,diffracti on and wave type diversi on even wave-guide and layer wave will happen,and relati onal para meters and ti m e-frequency charac2 teristics of waves are cl osely related t o the conditi ons of exp l osi on s ource,physical p r operties of p r opagati on mediu m,field charac2

teristics and terrain.The dissi pati on of wave energy caused by the inner fricti on and heat exchange of viscous-elastic mediu m s accompanied with its geometry attenuati on induces the gradual attenuati on of wave energy.

Key W ords:B lasting-caused seis m ic wave,Characteristics of wave,A ttenuati on and abs or p ti on

爆破是矿山开采中的一个重要环节。当炸药在岩体中爆炸时,一部分能量使炸药周围的介质引起扰动,并以波的形式向外传播。在爆破近区、中区传播的依次是冲击波、应力波,地震波由应力波在传播远区到达界面产生反射和折射叠加而形成[1],它是一种由爆源附近的应力波转换而来在岩土介质中传播的一种能量逐渐衰减的扰动,尽管只占爆炸所释放能量中的一小部分[2],但爆破地震波的特性对研究爆破地震波的传播机理、衰减规律及危害控制都具有重要意义。为此,人们从不同侧面对爆破地震波特性进行了大量研究[3~10]。

本文结合现场爆破震动监测信号,在查阅大量文献基础上,从爆破地震波的传播形式、传播方式、波的特征、波的衰减吸收及传播介质的力学模型等方面进行分析,探讨爆破地震波的特性。

1 爆破地震波的分类

爆破地震波包括在地层内部传播的体波和在地层表面或介质体表面传播的面波。体波可分为纵波(P波)、横波(S波);面波主要有Rayleigh波(R波)和Love波(L波)。

体波中的纵波指质点的振动方向与波的前进方向一致,使介质产生压缩和膨胀,因此又称为压缩波、疏密波、无旋转波或P(p re m ier)波。体波中的横波指质点振动方向与波的前进方向垂直,使介质被

I SS N1005-2763 CN43-1215/T D 矿业研究与开发第27卷第6期

M I N I N G R&D,Vol.27,No.6

2007年12月

Dec.2007

3收稿日期:2007-04-16

基金项目:国家自然科学基金(50764001);贵州省科技攻关项目(黔科合GY字(2007)3015);教育部“春晖计划”资助项目;贵州省优秀青年科技人才资助项目.

作者简介:张义平(1970-),男,湖南邵东人,博士,副教授,主要从事岩体力学、矿山开采及安全技术教学与科研工作,Email:c me.

ypzhang@https://www.360docs.net/doc/be1356892.html,。

剪切,故又称为剪切波、等体积波、旋转波或S (sec 2ond )波。S 波又分为S V 波与SH 波,其中S V 波是在入射平面内且与传播方向垂直振动的S 波,SH 波是在垂直于入射平面的方向上振动射线和平行的S 波。与S 波相比,P 波具有传播速度快、周期短、振幅小的特点。

面波是体波经地层界面的多次反射形成的次生波,是在地表或结构体表面以及结构层面传播的波,已发现存在R 波和L 波两种形式。半无限介质表面附近的R 波是P 、S 波由于边界作用产生的,其波

速V R 通过Rayleigh 方程可以确定[11]

。R 波传播时,质点在波的传播方向和自由面法线组成的平面内作椭圆运动,而在与该平面垂直的水平方向没有振动,其振动随深度呈指数衰减。

半无限弹性介质中P 波、S 波和R 波的传播速度之比和μ的关系见图1。可知V R

R 波不存在频散现象。

图1 半无限介质中P 波,S 波及R 波的

传播速度之比和泊松比μ的关系

只有在半无限空间上至少覆盖一低速地表层

时,L 波才会出现。L 波传播时,质点作与波的传播方向垂直的水平横向剪切振动,而无垂直分量运动,其传播速度介于最上层横波速度与最下层横波速度之间。

体波具有周期短、振幅小、衰减快的特点;面波其特点是周期长、振幅大、传播速度慢、衰减慢和携带的能量大。体积波特别是其中的P 波能使岩石产生压缩和拉伸变形,它是爆破时造成岩石破裂的主要原因,体波在爆破近区起主要作用;表面波特别是其中的R 波,由于它的频率低、衰减慢、携带较多的能量,是造成地震破坏的主要原因,面波在爆破远区起主要作用。地基土泊松比在1/4~1/2范围之

内,其纵波波速V P ≥1.75V S ,R 波波速V R =(0.9194~0.9554)V S 。所以从振源向外传播时,由于传播速度的差别,在远区体波与面波在时空上相互分开,P 波最先到达接收点,S 波次之,R 波最后到达。在弹性波的总输入能量中,R 波约占67%,剪切波占26%,压缩波占7%,即总能量的2/3是以R 波的形式向外传播。因此,在研究波的特性和抗震时,R 波的研究具有重要意义。

2 地震波传播速度

由于地震波在传播过程中遇到不同的地质条件,它将影响地震波的传播速度。

波速指波在介质中的传播速度,是波在介质中的能量衰减过程。质点振动速度指在外力作用下质点相对平衡位置作反复微幅运动时的速度,是质点运动能量的衰减过程。波速通常是波作用下的质点振动速度的几个数量级。在爆破地震效应分析中,通常关注质点振动速度而不太重视波的传播速度。实际传播介质并非理想的弹性介质而是粘弹性介质,波速不但和介质的成分、弹性、密度有关;还和介质的孔隙度、孔隙中所含流体的种类、相态有关。此外,它还和介质的埋藏深度、地质年代、经受地质

构造运动的历史等因素有关[12]

3 爆破地震波传播的特性

3.1 地震波的特征

由于爆源的复杂性(炸药、装药结构、爆破参数的多样性),传播介质的物理力学特性和地形地貌

的多变性,使得爆破地震波具有随时间作复杂变化

的随机不可重复的特性[13]

。不同条件下的现场实测爆破波形和频谱反映出,不同条件下的爆破所产生的震动波形是有明显区别的,不但在震动幅值上变化复杂,而且波的频率和持续时间也与爆源特性、爆心距、爆破规模及介质的不同显出明显的差异性。

爆破能量传播是一个衰减的过程,在各种不同介质中爆破地震波所包含的能量仅占爆炸总能量的3%~20%,其作用时间也较短,具有瞬态冲击振动的特性,其在爆破远区的危害性往往为人们所忽视。爆破地震波含有各种频率成份,是一种宽频带波。在传播过程中,由于介质的滤波作用,爆破地震波在离爆源较近时高频成份较丰富,随着波向远处的传播,高频逐渐吸收,而低频能传播到较远的距离。爆破地震波包含一个或几个主要的频率成份,

9

6 张义平,等: 爆破地震波特性研究

不同的频率成份对结构、设备及人员的影响也明显不同。大部分的爆破地震波频率主要集中在低频段,如果与结构的固有频率接近,就会产生共振现象,从而加大对结构体的破损影响,所以爆破地震波的频率特性不容忽视[14]。

3.2 地震波传播的方式

地震波在传播过程中将发生反射、透射与绕射、衍射、波型转换、波导、层间波等复杂现象,使地震波传播方向与途径发生变化。地震波在各种界面处垂直入射时要发生多次反射、透射与绕射。斜入射条件下还要发生P波与S波的波型转换,如入射的P 波或S V波常会产生反射与透射的P波和S V波,而SH波就只产生反射与透射的SH波。所以层面介质对地震波的传播速度和特性有重要影响。爆炸引起的地面振动是非常复杂的随机过程,测到的波形既有体波也有面波,它是不同幅值、不同频率与不同相位的各种波型叠加而成的复合波[15]。

3.3 波的衰减与吸收[16]

波按其波阵面的几何形状可以简单地分成球面波、柱面波和平面波。地震波在介质中传播时,通常以球面波或柱面波的形式传播,随着传播距离的增大,波阵面的几何形状将发生扩展,其波能在波阵面上均匀分布,从而使波能随传播距离R的增大而分别以R-2和R-1的速率减小,即发生几何衰减。若在理想弹性介质中,能量在波的传播过程中不会发生衰减(不考虑几何衰减),只在介质的分界面上产生透射和反射,仅使波能的传播方向发生改变,而不会出现能量形式的转换。而实际爆破地震波的传播介质是非理想的粘弹性介质,地震波在传播过程中要克服介质质点之间的内摩擦或粘滞作用,同时伴随波传播过程中的热传导及驰豫效应等现象,从而使波的能量在介质中传播时发生衰减,即能量被介质吸收,产生声能转换为热能的不可逆过程。即在粘弹性介质中波的衰减是几何衰减和传播介质吸收共同作用的结果,其中介质吸收又包括粘滞吸收和热传导吸收。

(1)粘滞吸收。介质的粘滞性是能量衰减的一个主要原因。在地震波传播时,相邻的体积元之间存在相对运动,并产生内摩擦或粘滞力,而地震波传播即介质形变(形状或体积的改变)传播过程中,需要能量来克服内摩擦或粘滞力,使得这种形变在传播过程中逐渐减小,即能量在传播过程中逐渐衰减。

(2)热传导吸收。热传导是导致介质中波能被吸收的另一主要原因。当波在介质中传播时,介质会产生体积压缩和膨胀的变化,压缩部分体积变小、温度升高,而膨胀部分体积变大、温度降低。在理想弹性介质中,温度的变化和体积的变化是同相的,即体积膨胀到极大值,温度为极小值;反之亦然,并且温度与体积间的变化是可逆的。而对于实际的非理想的粘弹介质,温度的变化和体积的变化总要滞后一段时间,另外在高温的压缩区和低温的膨胀区之间,不可避免地发生热能的转移,并且这种转移是不可逆的。最终使传播介质中产生热能从而导致介质温度升高。

(3)介质的品质因数。岩石既具粘滞性又具弹性,在应力波传播时既有能量损耗又有弹性动能与弹性位能的相互转换。岩石在一个周期性的波通过时其能量损耗可用下式描述:

1

Q

=

1

△E

E

=

Ca

πf

(1)

式中,Q为岩石(或介质)的品质因数,1

Q

为岩石(或介质)对应力波能量损耗的量度,C是岩石的声波速度,△E、E分别设为一个周期内损耗的能量和储存的弹性能。从式(1)可知:波能的损耗与a或△E

E

(一个周期内岩石或介质能量的相对损耗)成正比,同时与频率f成反比。实际工程上常用品质因数Q作为评价岩体是否完善及风化程度、裂隙发育情况的指标。

(4)频率的选择性吸收。波在介质中传播时,除高频比低频衰减更大外,介质还会对频率进行有选择性的吸收,例如实际地震记录上观察到的频率显著低于高能爆炸短促冲击可望产生的频率,还有军事上挖掘坑道而产生的噪音传播经过地层后听起来均属低频,接近听觉的下限,这些都说明地层对频率的吸收是有选择性的。描述这种类型的吸收频带最简单的方式是透射因子:

exp-f

f1

q x

x1

(2)式中:f———频率;

f1———恒定参考频率;

x1———常数;

x———波在介质中的传播距离;

q———决定吸收频带边缘陡度的参量。

如果确定x

1

,则f1就表示振幅在x=x1的距离上衰减至原振幅1/e时的频率,其中e为自然对数的

07矿业研究与开发 2007,27(6) 

底。描出函数exp -(f/f

1

)q」为代表的各种q值的一簇曲线(见图2),q可取任何实整数值,对于任意q

值,都形成高频在其中受到抑制而甚低频则透过的一种吸收谱。这些曲线表示可能被吸收的谱的频带。当q增大时,吸收频带边缘陡度增加;而当q无限增

大时,则所有低于f

1

的频率均自由透过,而所有高于f1的频率均完全被吸收掉。当然这里假设不存在频散,即所有频率成分的波均以相同速度传播通过大地[17]。

3.4 粘弹性介质的力学模型

由于传播介质的粘弹性特性,实际上爆破地震波的传播随着时间的增加而衰减,其能量耗散主要来自散射和吸收。传统经典弹性固体理论并没有考虑介质的粘弹性特征。1845年,St okes首先建立了St okes方程,从而使内摩擦作用在波动中有所反映。粘弹性介质的力学模型主要有3种:Max well模型, Kelvin模型,标准线性固体模型。

Max well模型同时具有固体及粘性流体性质,其变形速率随时间增长而趋于某一常数,属于粘性流体的范畴,但该模型不具备蠕变特征。Kelvin模型同时具有弹性固体及粘性流体性质,属于弹性固体的范畴,但该模型不能考虑应力作用下应变的突然变化,也不能表示应力消失后的剩余应变。即两者都不能充分描述大多数粘弹性介质的特性。而标准线性固体模型同时考虑了应变突然变化和剩余应变及蠕变特征,它是由Kelvin模型再串联一个弹簧组成的三元固体,如图3所示。此模型中固体既有瞬时弹性也有逼近弹性,其本构方程的标准形式为:

σ+p

1dσ

d t

=q0ε+q1

d t

(3)

其中:p

1=

η

E1+E2

,q0=

E1E2

E1+E2

,q1=

E1η

E1+E2

式中:σ———应力;

ε———应变;

E———弹性模量;

η———粘性系数。

在粘弹性介质中,波动能量的耗散不同于完全弹性介质中能量的耗散,它主要与介质的特性和振动频率等因素相关。若将粘弹性纵波和横波的衰减

系数分别记为a

p 、a

s

,经推导有[18]:

a p=-ωI m

ρ

λ3(iω)+2G3(iω)

(4)

a s=-ωI m

ρ

G3(iω)

(5)

式中:I

m

———复数的虚部;

ρ———岩石的体积密度;

ω———波的圆频率;

λ3、G3——

—粘弹性介质与弹性介质中λ、G

对应的量。

由上看出:粘弹性介质中波的衰减与频率有关,

频率越高其衰减程度越大。

图2 被选择吸收的频谱曲线

图3 标准线性固体模型

4 爆破地震波与自然地震波

爆破地震和自然地震都属于能量释放引起地表

振动的现象,它们所引起的振动有如下相似之处:

(1)两者突然释放的能量均以波的形式通过介

质从振源向外传播,并引起强烈的地表或建(构)筑

物的振动;

(2)两者的质点振动强度与振源能量和振源距

离紧密相关;

(3)质点的振动参数都明显地受地质、地形等

因素的影响;

(4)两者对结构体的破坏机制是相同的。

由于两者多方面的相似性,人们常常将自然地

震领域相对完善的分析理论和方法应用于爆破领

域,如信号分析技术和反应谱理论等。但在工程实

际中二者又有明显的差异,具体表现在:

(1)自然地震频率都很低,一般低于5Hz。爆

破地震频率复杂且较高,常在0~200Hz的范围内。

与建(构)筑物结构的固有频率相比,前者与之接

近,而后者却高得多;

17

 张义平,等: 爆破地震波特性研究

(2)目前世界上测到的地震最大加速度约为1.3g,而在大爆破附近测得的地表振动加速度高达25.3g。尽管爆破地震波振动幅值高,但由于频率高和能量小其衰减很快,而自然地震波衰减缓慢得多;

(3)爆破地震持续时间一般不超过0.5s,如果是雷管段数增加和接力,也不会超过数秒。而自然地震常持续达数分钟甚至更久;

(4)由于自然地震频率低、衰减慢、持续时间长和携带能量巨大,其所造成的损失远远超过爆破地震所带来的危害;

(5)爆破震源的大小和位置以及作用方向可以控制,爆破震动的延续时间可以知道,从而可以控制爆破震动对结构的效应。

5 结 语

(1)爆破地震波是爆炸应力波转换而来在岩土介质中传播的一种能量逐渐衰减的扰动,它以体波和面波的形式向外传播,爆炸引起的地面振动是非常复杂的随机过程,它是不同幅值、不同频率与不同相位的各种波型叠加而成的复合波。体波在爆破近区起主要作用;爆破远区,表面波特别是R波是造成地震破坏的主要原因。

(2)爆破地震波富含各种频率成份,具有瞬态性、随机性和危害性的特征。地震波在传播过程中将发生反射、透射与绕射、衍射、波型转换、波导、层间波等复杂现象,使地震波传播方向与途径发生变化。传播过程中波的有关参数和时频特征常与爆源条件、传播介质的物理性质、结构特征及地形地貌等因素密切相关。

(3)爆破地震波的传播介质是非理想的粘弹性介质,在传播过程中波的能量衰减是几何衰减、传播介质粘滞吸收和热传导吸收共同作用的结果,标准线性固体模型是以St oke方程为基础、同时考虑了应变突然变化和剩余应变及蠕变特征的力学模型。

(4)爆破地震波与自然地震波具有诸多共同之处,地震信号的分析方法及相应的抗震理论非常类似,但信号的时频特征及其危害性又各不相同。

参考文献:

[1]孟吉复,惠鸿斌.爆破振动测试技术[M].北京:冶金工业出版

社,1992:80~138.

[2]张 奇,白春华,刘庆明.爆炸地震波频谱特性研究[J].北京

理工大学学报,1999,19(3):306~308.

[3]李 铮,朱瑞赓.爆炸地震波震速的特征系数与衰减指数的研

究[J].爆炸与冲击,1986,6(3):221~229.

[4]朱瑞赓,李 铮.爆炸地震波的现场测量及其安全距离[J].

爆炸与冲击,1982,2(1):60~67.

[5]宋光明.爆破震动小波包分析理论与应用研究[D].长沙:中

南大学,2001.

[6]W ada Yuji,Ma Gui Chen,Nakaji m a Toshihide,et al.Study on the

contr ol of vibrati on caused by a blasting[J].Kayaku Gakkaishi/

Journal of the Japan Exp l osives Society,1994:55(4):174~180.

[7]Rodgers Jay.The effect of confine ment on gr ound vibrati on a mp li2

tude[J].Pr oceedings of the Annual Conference on Exp l osives and

B lasting Technique,2003,(II):131~144.

[8]Dowding C H.I n Engineering Handbook[M].Prentice Hall,

1985:746~760.

[9]李夕兵,张义平,刘志祥,等.爆破震动信号的小波分析与

HHT变换[J].爆炸与冲击,2005,25(6):528~535.

[10]凌同华,李夕兵.地下工程爆破震动信号能量分布的小波包分

析[J].爆炸与冲击,2004,24(1):63~68.

[11]徐 建.建筑振动工程手册[M].北京:中国建筑工业出版社,

2002:21~28.

[12]张雪亮,黄树棠.爆破地震效应[M].北京:地震出版社,1981:

1~211.

[13]林大超,施惠基,白春华,等.爆炸地震效应的时频分析[J].

爆炸与冲击,2003,23(1):31~35.

[14]张义平.爆破震动信号的HHT分析与应用研究[D].长沙:中

南大学,2006.

[15]熊代余,顾毅成.岩石爆破理论与技术新进展[M].北京:冶金

工业出版社,2002:13~165.

[16]楚泽涵.声波测井原理[M].北京:石油工业出版社,1987:49

~58.

[17](美)N H瑞克.粘弹性介质中的地震波[M].北京:地质出版

社,1981:7~51.

[18]吴爱祥,孙业志,刘湘平.散体动力学理论及其应用[M].北

京:冶金工业出版社,2002:163~169.

(上接第41页)

锌硫矿物的有效抑制剂,S N9对铅矿物也能选择性捕收,闭路试验中较好地实现了微细粒铅锌的浮选分离。

参考文献:

[1]黄传兵,陈兴华,等.选择性絮凝技术及其在矿物分选中的应

用[J].矿业工程,2005,6(3):27.[2]岳铁兵,魏德洲,曹进成,等.细粒金红石矿浮选工艺研究

[J].化工矿物与加工,2005,(3):1~3.

[3]魏宗武,叶雪均,陈建华,等.难选铅锑锌分离新工艺研究

[J].中国矿业,2006,9(9):88~90.

[4]刘智林,魏宗武,叶雪均,等.某选矿厂矿泥系统铅锌分离探索

研究[J].矿产保护与利用,2005,6(3):34~36.

[5]刘红缨,朱书全,等.提高选择性絮凝脱硫效果的途径[J].化

工矿物与加工,2001,7(7):15~17.

27矿业研究与开发 2007,27(6) 

喷管流动特性与管道截面变化规律的关系

喷管流动特性与管道截面变化规律的关系 摘要:针对管内流动规律的一般应用中存在的问题,着重讨论了喷管内工质流动特性与管道截面变化规律的关系,从而更准确更完整地反映了喷管内工质流动规律。 关键词:喷管;流动特性;变化规律 通常在研究喷管内工质流动特性时,只着重于对喷管外形的确定,所以总是以状态参数变化为前提,去探讨工质流动截面(即管道截面)的相应变化。这时由可逆绝热流动的基本方程组,即连续性方程、能量方程和过程方程,整理出如下两个关系式: 很明显,式(1)、(2)反映了工质流速c、压力P、截面A之间的变化关系。从数学角度而言,这几个量是可以互为变化前提的。但对具体的管内流动来说,究竟谁是其中的决定性因素,从而控制着(导致)其它两个量的相应变化,这自然是一个非常重要的问题。但这一问题在很多文献[1~3]中并无明确地阐述。 显然,要揭示清楚喷管内工质的流动规律,必须揭示清楚上式中各个量的决定与被决定关系,不然问题的实质就不会充分地显现出来,所得结论也是不完整的,也就无法满足实际应用的需要。特别是个别文献还错误地强调了这种关系,从而让人产生各种疑惑甚至是误解。这也是许多人在学习了喷管内流动特性之后,对一些管内流动现象还仍然解释不清,甚至出现概念上的错误的根本原因。 1对喷管内流动特性与管道截面变化规律关系的分析 任何一种流动都是在一定的外部条件作用下产生的。随流动条件的不同,管内流动现象才是多种多样的。就喷管流动而言,其流动条件应包括如下两个方面:(一)力学条件:即喷管前后的压差;(二)几何条件:即喷管长度L和喷管流动方向(设为x方向)的截面变化规律A=f(x)。 工质降压升速、升压减速等流动特性,即工质压力P、比容v、流速c包括流动截面A的相互变化关系,应属流体自身属性,这种属性不会自发地表现出来,它是从属于流动的外部条件而存在的。这里的力学条件是工质流动和膨胀的动力,几何条件是工质连续降压增速的保证。在流动产生前和流动过程中,其力学条件和几何条件都是客观的,两者共同确定了相应的流动特性,缺一不可。比如,即使在力学条件完全具备的情况下,若没有几何条件的保证,流体降压升速等属性也不会自发地表现出来。对此还可以用一个简单的例子来加以说明:设流动的 力学条件为初压P 1与背压P b ,在流动产生之前,只有P 1 、P b 是客观存在的,P 1 与P b 之间的其它压力以及其它参数都不是客观的。只有在流动产生之后才在各

在 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 实验目的 (1)深入了解减压器工作原理及其工作特性。 (2)研究减压器的静态特性,掌握测定减压器静态特性的方法,掌握减压器静态特性的一般规律。 (3)了解减压器的过渡过程压力曲线测定方法,增加对减压器动态特性的感性认识。 2 实验背景 2.1减压器的应用 减压器不仅广泛应用于油、气工业、化工行业、能源工业、基础设施建设等行业,在航空航天领域也发挥着重要作用。在航天行业中,减压器可应用于地面设备(包括地面试验设备)、导弹/运载火箭和卫星航天器。具体而言,减压器可用于: (1)地面试验吹除系统。受系统工作压力的限制,此类减压器出口压力较低,精度要求也不是很高,但质量流量大,要求有较好的启动稳定性。 (2)地面试验或弹箭体供气系统。对于使用气体推进剂的地面发动机试验系统或弹箭体而言,其供气系统中都必须使用到减压器,以保证稳定的压力和流量供应,对减压器的精度!动态特性要求较高。 (3)地面试验或弹箭体液体推进剂输运系统。减压器为推进剂储箱提供恒定的压力,进而为发动机提供需要的推进剂,其出口压力影响到发动机的工作状态,直接关系到整个系统推进剂供应的准确性与安全性,是影响整个发动机推力稳定性的一个重要因素,因此对减压器精度要求较高。 (4)航天器的姿态和轨道控制。在卫星、探空火箭、宇航控制系统、空间站对接操纵系统中以及弹体姿态控制系统中的的冷气推进系统中,减压器出口的气体直接送至喷管进行姿态或轨道控制,具有开启次数频繁,流量变化大的特点,对动态特性、工作范围、控制精度、可靠性和寿命都有较高的要求。 (5)提供基准压力或控制其它调节器。利用减压器出口压力稳定的特点,

地震波的概念种类特点及地表影响

关于地震波 摘要:地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面,将地球内部分为地壳、地幔和地核三个圈层。 关键词:地震波辐射地球内部 一:背景 ①2008年5月12日14时28分04秒,四川汶川、北川,8级强震猝然袭来,大地颤抖, 山河移位,满目疮痍,生离死别……西南处,国有殇。这是新中国成立以来破坏性最强、波及范围最大的一次地震。此次地震重创约50万平方公里的中国大地!为表达全国各族人民对四川汶川大地震遇难同胞的深切哀悼,国务院决定,2008年5月19日至21日为全国哀悼日。自2009年起,每年5月12日为全国防灾减灾日。 ②1976年7月28日北京时间03时42分53.8秒,在中国河北省唐山、丰南一带(东经118.2°, 北纬39.6°)发生了强度里氏7.8级(矩震级7.5级),震中烈度Ⅺ度,震源深度23千米的地震。地震持续约12秒。有感范围广达14个省、市、自治区,其中北京市和天津市受到严重波及。强震产生的能量相当于400颗广岛原子弹爆炸。整个唐山市顷刻间夷为平地,全市交通、通讯、供水、供电中断。唐山地震没有小规模前震,而且发生于凌晨人们熟睡之时,使得绝大部分人毫无防备,造成24.2万人死亡,重伤16.4万人,名列20世纪世界地震史死亡人数第一。 ③邢台地震由两个大地震组成:1966年3月8日5时29分14秒,河北省邢台专区隆尧县 (北纬37度21分,东经114度55分)发生震级为6.8级的大地震,震中烈度9度强; 1966年3月22日16时19分46秒,河北省邢台专区宁晋县(北纬37度32分,东经115度03分)发生震级为7.2级的大地震,震中烈度10度。两次地震共死亡8064人,伤38000人,经济损失10亿元。这是一次久旱之后的大震。

热工学实验

实验十 渐缩(缩放)喷管内压力分布和流量测定 一、实验目的 1.验证并加深对喷管中的气流基本规律的理解,树立临界压力,临界流速,最大流量等喷管临界参数的概念,把理性认识和感性认识结合起来。 2.对喷管中气流的实际复杂过程有概略的了解。 3.通过渐缩喷管气流特性的观测,要明确:在渐缩喷管中压力不可能低于临界压力,流速不可能高于音速,流量仍不能大于最大流量。 4.根据实验条件,计算喷管(最大)流量的理论值,并与实侧值进行对比。 二、实验设备 本设备由2x 型真空泵,PG -Ⅲ型喷管(见图10-1)和计算机(控制与显示设备)构成。由于真空泵的抽吸,空气自吸气口2进入进气管1,流过孔板流量计3,流量的大小可以从U 型管压差计4读出。喷管5用有机玻璃制成,有渐缩、缩放两种型式(见图10-2、10-3),可根据实验要求,松开夹持法兰上的螺丝,向右推开进气管的三轮支架6,更换所需的喷管。喷管各截面上的压力是由插在其中,外径0.2mm 的测压探针连至可移动真空表8测得,探针的顶封死,中段开有测压小孔,摇动手轮——螺杆机构9,即可移动探针,从而改变测压小孔在喷管中的位置,实现对喷管不同截面的压力测量。在喷管的排气管上装有背压真空表10,排气管的下方为真空罐12,起稳定背压的作用,背压的高低用调节阀11调节。罐前的调节阀用作急速调节,罐后的调节阀作缓慢调节,为减少震动,真空罐与真空泵之间用软管13连接。 在实验中必须观测四个变量:(1)测压孔所在截面至喷管进口的距离x ;(2)气流在该截面上压力P ;(3)背压P b ;(4)流量m 。这些变量除可分别用位移指针的位置、移动真空表,背压真空表及 U 形管压差计的读数来显示读出外,还可分别用位移电位器、负压传感器、压差传感器把它们转换为电信号,由计算机显示并绘出实验曲线。位移电位器将在螺杆之旁,它实际上是一只滑杆变阻器。负压传感器和压差传感器分别装在真空表和U 形管压差计附近,其内部结构为一直流电桥,压力和压差改变时将改变电桥中两臂的电阻,从而获得电桥的不平衡电压输出。为了使这些传感器可靠而稳定地工作,都由直流稳压电源供电。 三、实验原理 1.喷管中气流的基本规律 气流在喷管中稳定流动后,喷管任何截面上的质量流量m 均相等,有连续性方程: M= 2 2 21 1 1C A C A AC υυυ = = =定值,[kg/s] (10-1) 式中:A —— 截面积[m 2] C —— 气体流速[m/ s] υ —— 气体比容[m 3/kg] 下标1—— 喷管进口 下标2——喷管出口 气体在喷管中作绝热膨胀,C 1<C 2,工质为理想流体时,喷管的理论流量可按下式计算: ])()[(121 1 22 12112 2 2 2k k k p p p p p k k A C A m +-?-== υυ (10-2) 式中: k —— 绝热指数,对于空气k=1.4 P 1 —— 喷管进口压力(初压) [N/ m 2] P 2 —— 喷管出口压力 [N/ m 2] 喷管中气体状态参数P 、υ和流动参数C 的变化规律和流通截面积A 的变化以及喷管

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

时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系,能真实地反应结构地震相应随时间变化的全过程,是抗震分析的一种重要方法[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. 地震波动力学:研究地震波在传播过程中波形、振幅、频率、相位等特征的及其变化规律,以及这些变化规律与地下的地层结构,岩石性质及流体性质之间存在的联系。 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.验证喷管中气流的基本规律,加深对临界压力、临界流速和最大流量等喷管临界参数的理解。 2.比较熟练地掌握压力、压差及流量的测量方法。 3.重要概念1的理解:应明确在渐缩喷管中,其出口处的压力不可能低于临界压力,流速不可能高于音速,流量不可能大于最大流量。 4.重要概念2的理解:应明确在缩放喷管中,其出口处的压力可以低于临界压力,流速可高于音速,而流量不可能大于最大流量。 二、实验装置 整个实验装置包括实验台、真空泵(规格为1401型,排气量3200L/min)。实验台由进气管、孔板流量计、喷管、测压探针、真空表及其移动机构、调节阀、真空罐等几部分组成,如图6-4所示。 图6-4 喷管实验台 1-进气管;2-空气吸气口;3-孔板流量计;4-U形管压差计;5-喷管; 6-三轮支架; 7- 测压探针; 8-可移动真空表; 9-位移螺杆机构及位移传感器; 10-背压真空表; 11-背压用调节阀;12-真空罐;13-软管接头;14-仪表箱;15-差压传感器;16-被压传感器;17-移动压力传感器 进气管为φ57×3.5无缝钢管,内径φ50。空气从吸气口入进气管,流过孔板流量计。孔板孔径φ7,采用角接环室取压。流量的大小可从U形管压差计或微

压传感器读出。喷管用有机玻璃制成,配有渐缩喷管和缩放喷管各一只。根据实验的要求,可松开夹持法兰上的固紧螺丝,向左推开进气管的三轮支架,更换所需的喷管。喷管各截面上的压力是由插入喷管内的测压探针(外径φ1.2)连至“可移动真空表”测得,由于喷管是透明的,测压探针上的测压孔(φ0.5)在喷管内的位置可从喷管外部看出,它们的移动通过螺杆机构移动,标尺或位移传感器实现测量读数。喷管的排气管上还装有“背压真空表”,其压力大小用背压调节阀进行调节。真空罐直径φ400,起稳定压力的作用。罐的底部有排污口,供必要时排除积水和污物之用。为减小震动,真空罐与真空泵之间用软管连接。 在实验中必须测量四个变量,即测压孔在喷管内的不同截面位置X 、气流在该截面上的压力P 、背压P b 、流量m ,这些量可分别用位移指针的位置、可移动真 空表、背压真空表以及U 形管压差计的读数来显示。 实验装置特点: 1.可方便地装上渐缩喷管或缩放喷管,观察气流沿喷管各截面的压力变化。 2.可在各种不同工况下(初压不变,改变背压),观察压力曲线的变化和流量的变化,从中着重观察临界压力和最大流量现象。 3.除供定性观察外,还可作初步的定量实验。压力测量采用精密真空表,精度0.4级。流量测量采用低雷诺数锥形孔板流量计,适用的流量范围宽,可从流量接近为零到喷管的最大流量,精度优于2级。 4.采用真空泵为动力,大气为气源。具有初压初温稳定,操作安全,功耗和噪声较小,试验气流不受压缩机械的污染等优点。喷管用有机玻璃制作,形象直观。 5.采用一台真空泵,可同时带两台实验台对配给的渐缩、缩放喷管做全工况观测。因装卸喷管方便,本实验台还可用作其他各种流道喷管和扩压管的实验。 三、实验原理 1、喷管中气流的基本规律 (1)由能量方程: 221dc dh dq += 及 dp dh dq ν-= 可得 cdc dp =-ν 可见,当气体流经喷管速度增加时,压力必然下降。 (2)由连续性方程: 有 及过程方程 常数=k p ν 常数=?=??????=?=?νννc A c A c A 222111c dc d A dA -=νν

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 ………………

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

工程热力学喷管特性实验

实 验 报 告 评分 实验题目:喷管特性实验 实验目的:验证并进一步加深对喷管中气流基本规律的理解,建立临界压力、临界流速 和最大流量等喷管临界参数的概念;比较熟练地掌握用热工仪表测量压力(负压)、压差及流量的方法;明确渐缩喷管出口处的压力不可能低于临界压力,流速不可能高于音速,流量不可能大于最大流量;明确缩放喷管中的压力可以低于临界压力,流速可高于当地音速,而流量不可能大于最大流量;对喷管中气流的实际复杂过程有所了解,能定性解释激波产生的原因。 实验原理: 1.喷管中气流的基本原理 由连续方程、能量方程和状态方程结合声速公式KPV a =得: c dc M A dA ? ?? ? ?-=12 马赫数M=c/a 显然,要使喷管中气流加速,当M<1时,喷管应为渐缩型(dA<0);当气流M>1时, 喷管应为渐扩型(dA>0)。 2.气体流动的临界概念 喷管中气流的特征是dp<0,dc>0,dv>0,三者之间互相制约。当某一截面的速度达到当地音速时,气流处于从亚音速变为超音速的转折点,通常称为临界状态。 临界压力比112-? ?? ??+=K K K ν ,对于空气,ν=0.528 当渐缩喷管出口处气流速度达到音速或缩放喷管喉部达到音速时,通过喷管的气体流量 便达到了最大值,或成临界流量。可由下式确定: 1112 1212m i n m a x V P K K K K A m ?-??? ??++= 式中: min A —最小截面积(对于渐缩喷管即为出口处的流通截面积;对于缩放喷管即为喉部的面 积。本实验台的两种喷管最小截面积均为11.44)。 3.气体在喷管中的流动 (1)渐缩喷管 渐缩喷管因受几何条件(dA<0)的限制。有式(4)可知:气体流速只能等于或低于音速(a C ≤);出口截面的压力只能高于或等于临界压力(c P P ≥2);通过喷管的流量只能等于或小于最大流量(max m m =)。 (2)缩放喷管

工程热力学喷管特性实验

工程热力学喷管特性实验 实验报告评分 实验题目:喷管特性实验 实验目的:验证并进一步加深对喷管中气流基本规律的理解,建立临界压力、临界流速 和最大流量等喷管临界参数的概念;比较熟练地掌握用热工仪表测量压力 (负压)、压差及流量的方法;明确渐缩喷管出口处的压力不可能低于临界 压力,流速不可能高于音速,流量不可能大于最大流量;明确缩放喷管中的压力可以低于临界压力,流速可高于当地音速,而流量不可能大于最大流量; 对喷管中气流的实际复杂过程有所了解,能定性解释激波产生的原因。实验原理: 1(喷管中气流的基本原理 a,KPV由连续方程、能量方程和状态方程结合声速公式得: dAdc2,,,M,1,,,,Ac 马赫数M=c/a 显然,要使喷管中气流加速,当M<1时,喷管应为渐缩型(dA<0);当气流M>1时,喷管应为渐扩型(dA>0)。 2(气体流动的临界概念 喷管中气流的特征是dp<0,dc>0,dv>0,三者之间互相制约。当某一截面的速度达到当地音速时,气流处于从亚音速变为超音速的转折点,通常称为临界状态。 K 2,,K,1,,,,K,1,, 临界压力比,对于空气,,=0.528 当渐缩喷管出口处气流速度达到音速或缩放喷管喉部达到音速时,通过喷管的气体流量便达到了最大值,或成临界流量。可由下式确定:

2P2K2,,K,11,m,A,,,maxminK,1K,1V,,1 式中: A—最小截面积(对于渐缩喷管即为出口处的流通截面积;对于缩放喷管即为喉部的面min 积。本实验台的两种喷管最小截面积均为11.44)。 3(气体在喷管中的流动 (1)渐缩喷管 渐缩喷管因受几何条件(dA<0)的限制。有式(4)可知:气体流速只能等于或低于音 P,P2cC,a速();出口截面的压力只能高于或等于临界压力();通过喷管的流量只能等 ,,m,mmax于或小于最大流量()。 (2)缩放喷管 缩放喷管的喉部dA=0,因而气流可达到音速(c=a);扩大段dA>0,出口截面处的流速可超音速(c>a),其压力可低于临界压力(P2

燃烧学实验报告1

燃烧学 实验报告 学院 专业 学号 学生姓名 指导教师 2017 年 01 月

目录 实验原理系统图、实验仪器仪表型号规格及燃料物理化学性质 (3) 实验一Bensun火焰及Smithell法火焰分离 (6) 实验二预混火焰稳定浓度界限测定 (8) 实验三气体燃料的射流燃烧、火焰长度及火焰温度的测定 (11) 实验四本生灯法层流火焰传播速度的测定 (15) 燃烧喷管及石英玻璃管说明 燃烧喷管共4根,分别标记为: I号长喷管—细的长喷管(喷口内径7.18mm) II号长喷管—粗的长喷管(相配的冷却器出口直径10.0mm) I号短喷管—细的短喷管(喷口内径5.10mm) II号短喷管—粗的短喷管(喷口内径7.32mm) 石英玻璃套管共3个,分别标记为: I号玻璃管—最细的石英玻璃管(本生灯火焰内外锥分离用) II号玻璃管—中间直径的石英玻璃管(观察Burk-Schumann火焰现象及测定射流火焰长度用) III号玻璃管—最粗的石英玻璃管(测定射流火焰温度用)

燃烧学实验注意事项 1.实验台上的玻璃管须轻拿轻放,用完后横放在实验台里侧,以防坠落。 2.燃烧火焰的温度很高,切勿用手或身体接触火焰及有关器件。 3.燃烧完后的喷嘴口、水平石英管的温度仍很高,勿碰触,以防烫伤。 4.在更换燃烧管时,手应握在下端,尽量远离喷嘴口。 本生灯燃烧实验系统介绍 一、本生灯燃烧实验系统 本生灯燃烧实验系统如图1所示。压缩空气(蓝色导管)通过减压阀门调节,进入浮子流量计调节测量。甲烷(桔红导管)由气瓶开关、减压阀提供,接通单向电磁阀进入浮子流量计调节测量。主控制面板上设计了一个主控单向电磁阀开关,只有当电源接通,开关按下时,燃料才能供应。空气和燃气调节好当量比被分别输送进入一个混合管(混合管的功用是缓和气流的脉动,并使甲烷、压缩空气两股气体充分混合,以保证在本生灯管的进口处获得稳定、均匀的流动),该混合管连接一个内径为25mm的本生灯管,用于火焰结构演示实验、钝体火焰稳定实验。根据火焰传播的余弦定理,已知管口直径,测量出火焰的高度,就能算出火焰锥角。再通过浮子流量计的读数,测出气流流出管口的速度,就可以算出火焰的传播速度。 图 1 本生灯燃烧实验系统

喷管中气体流动基本特性实验报告

喷管中气体流动基本特性实验报告 一、实验目的 1. 验证并进一步对喷管中气流基本规律的理解。牢固树立临界压力、临界流速和最大流量等喷管临界参数的概念。 2. 掌握喷管实验装置的实验原理、实验方法和操作步骤,比较熟练地用热工仪表测量压力(负压)、压差及流量。 3. 测量并绘制喷管内的压力分布曲线及流量曲线,做出定性的解释。 二、实验原理 喷管是一些热工设备的重要部件,这些设备的工作过程和喷管中气体的流动过程有密切的关系。实验观察气流完全膨胀时沿喷管各界面的压力变化,测定流量曲线和临界压力比,可以帮助了解喷管中气体流动现象的基本特性,并且通过观察渐缩渐扩喷管中膨胀不足和膨胀过度的现象,还可进一步了解工作条件对喷管中流动过程的影响。 气体在喷管的流动过程中,气体的状态参数P 、V ,流速C 和喷管截面积f 之间的基本关系可用下面三个方程表示: c dc f df v dv f df c dc vdp cdc M )1(02 -==-+-= (4-1) 式中:M 为马赫数,是表示气体流动特性的一个重要特性值。M<1时,表明气体流速小于当 地音速,M>1时,气体流速大于当地音速,气体作超音速流动。 方程指出:气体流经喷管时,压力降低,流速增大,喷管的截面积亦随之变化,而喷管的截面变化情况则取决于M值. 1) 当气流流速小于音速(即M<1)时,欲使流速增大,喷管截面应该是收缩的; 2) 当气流流速大于音速(即M>1)时,喷管截面应该是扩放的; 3)当流速等于音速时,喷管截面最小,此处正是气流流速由亚音速过渡到超音速,喷管由收缩形过渡到扩放形的转折点。这点的参数称为喷管的临界参数,用脚码C 表示,如临界压力P C 、临界流速C C 等等。 1.渐缩喷管 气体流经喷管的膨胀程度可以用喷管的背压P 2与进口压力P 1之比β表示。P P 1 2= β称 为压力比。而气体在渐缩喷管中膨胀所能达到的最低压力,是使喷管出口的气流速度达到当地音速的压力,即临界压力P C 。所以,气流在渐缩喷管中流动时最大膨胀程度决定于临界压

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后的),然后将一组人

相关文档
最新文档