基于Matlab求解建筑结构地震响应的时程分析法_孟宪萍 (1)

2008年第6期总第120期

福 建 建 筑

F u j i a nA r c h i t e c t u r e ＆C o n s t r u c t i o n

N o 6·2008

V o l ·120

基于M a t l a b 求解建筑结构地震响应的时程分析法

孟宪萍

(开封市供水总公司　475004)

摘　要:本文基于m a t l a b 阐述了我国《建筑抗震设计规范》(G B 50011-2001)规定的求解建筑结构地震响应的时程分析法,应用m a t l a b 语言编制了时程分析法求解建筑结构地震响应的计算程序,并以一三层钢筋混凝土结构为工程算例,应用基于m a t l a b 的时程分析法进行结构的地震响应计算。结果表明,基于m a t l a b 的时程分析计算效率较高。关键词:M A T L A B 地震响应　时程分析法

中图分类号:T U 312+.1 文献标识码:A 文章编号:1004-6135(2008)06-0038-03

T h e t i m e -h i s t o r y m e t h o db a s e d o nm a t l a b o f r e s o l v i n g t h e e a r t h q u a k e r e s p o n s e o f t h e s t r u c t u r e s

M e n g X i a n p i n g

(K a i f e n g Wa t e r S u p p l y C o m p a n y 475004)

A b s t r a c t :I nt h i s p a p e r ,t h e t i m e -h i s t o r y m e t h o dw h i c h i s m e n t i o n e d i n t h e c o d e f o r s e i s m i c d e s i g n o f b u i l d i n g s (G

B 50011-2001)t o

r e s o l v e t h e e a r t h q u a k e r e s p o n s e o f t h e s t r u c t u r e s i s d i s c u s s e d b a s e d o nm a t l a b .T h e c a l c u l a t i o np r o g r a m s o f t h et i m e -h i s t o r y m e t h o da r e w o r k e do u t u s i n g t h e l a n g u a g e m a t l a b .T a k i n g a t h r e es t o r y r e i n f o r c e d c o n c r e t e f r a m e s t r u c t u r e a s a ne x a m p l e ,t h e e a r t h q u a k e r e s p o n s e o f t h e s t r u c t u r e i s r e s o l v e d b y u s i n g t h e c a l c u l a t i o n p r o g r a m s o f t h e t i m e -h i s t o r y m e t h o d .T h e r e s u l t i n d i c a t e s t h a t T h e t i m e -h i s t o r y m e t h o d b a s e do n m a t l a bo f r e s o l v i n gt h e e a r t h q u a k e r e s p o n s e o f t h e s t r u c t u r e s i s e f f i c i e n t .K e y w o r d s :

M A T L A B e a r t h q u a k e r e s p o n s e t i m e -h i s t o r y a n a l y s i s m e t h o

d 作者简介:孟宪萍,女,1966年出生,主要从事建筑结构设

计及建筑咨询。

收稿日期:2008-03-25

1　引言

我国《建筑抗震设计规范》(G B 50011-2001)第5章对时程分析法的使用情况作出了规定。时程分析法又称为直接动力法或逐步积分法。采用时程分析法可以计算出结构在地震过程中每一瞬时的反应,可用来求解建筑结构的几何及物理线性与非线性动力响应。与经典的反应谱方法相比,有很多的优点,但是它也存在许多不足,主要有计算模型的合理选择困难;地震波输入的不确定性;在计算过程中要进行刚度矩阵等的不断修正,每一时刻的结果都受到此刻之前的结果的影响等,导致计算分析工作量较大。虽然目前在结构弹塑性时程分析时,结构动力增量微分方程已有较为成熟的算法以及相关的大型分析软件可以利用,但是其计算分析工作量仍然十分繁重,不但耗费机时,结果处理复杂,而且同计算者本身的经验和对结构在地震作用下的损伤形态和破坏顺序的

假定相关,这些都带有一定的主观性。但是随着计算机的普及,时程分析法正逐步被抗震规范接受。本文在详细阐述了时程分析法基本原理基础上,结合m a t -l a b 语言编制了时程分析法求解建筑结构地震响应的计算程序,并以一三层钢筋混凝土结构为例进行验证。

2　时程分析法基本原理

2.1　结构在地震作用下的动力微分方程

多自由度体系建筑结构在地震作用下的运动运动微分方程为

[M ]{x ··

}+[C ]{x ·

}+[K ]{x }=-[M ]{x ··

g }

(1)

其中,[M ],[C ],[K ]分别为建筑结构质量、阻尼和刚度矩阵,{x ··

g }为地面运动加速度。2.2　建筑结构的计算模型

建筑结构计算模型一般应根据结构形式及构造特点、分析精度以及计算机容量等情况确定。用时程分析法求解时,由于计算工作量大,在尽量真实再现结构动力反应特点的前提下,尽可能对结构予以简化。对于传统的多层房屋结构,最简单且应用最广的模型是层间剪切模型,如图1所示,在这种模型中,房

屋的质量集中在各楼层,在振动过程中各楼层式中保持为水平,结构的变形表现为层间错动,各层的层间位移互不影响。对于以剪切变形为主的结构,一般都可以采用这种模型,如多层砖房以及横梁线刚度远比柱线刚度大的强梁弱柱型框架结构.但是对于强柱弱梁型的框架结构,用这种模型计算时误差较大,但有时为了简化计算,对于各跨相等的底层框架和建筑物宽度远大于高度的多层框架结构亦可以近似地应用

。

图1　层间剪切型模型计算简图　图2　杆系模型计算简图

为了精确计算可以采用如图2所示的杆系模型。该模型假定楼板在自身平面内为绝对刚性,以构件作为基本杆件单元而将梁柱简化为以中性轴表示的无质量杆,将质量集中于各结点,利用构件连接处的变形协调条件建立各构件变形关系,利用构件的恢复力特性集成整个结构的弹塑性刚度

[1]

,然后采用数值积

分方法对结构进行地震反应分析。这种模型一般适用于强柱弱梁的框架结构,可以用它来求出地震过程中各杆逐渐开裂并进入塑性阶段的过程以及整个结构的影响,但是计算较繁。我国《建筑抗震设计规范》(G B 50011-2001)规定,规则结构可以使用平面杆系模型计算罕遇烈度地震作用下的结构弹塑性变形,对于不规则结构应采用空间杆系结构模型。2.3　输入地震波的选择

2.3.1　选波意义,结构的地震反应是输入地震和结构自身的动力特性共同作用的结果。在采用时程分析法求解结构地震反应时,需要首先确定输入结构的地震地面运动加速度。而观测结果表明,地震波具有强烈的随机性,结构的弹塑性时程分析结果也表明,结构的地震反应随输入地震波的不同会有很大不同,误差有时会高达几倍至几十倍,所以要保证时程分析结果的可靠性,必须正确合理地选择地震波。2.3.2　选择方法及选用数量,目前关于地震波的选择有三种方法:(1)拟建场地的实际地震记录。(2)过去典型的强震记录。(3)人工合成场地地震波。

过去典型强震记录常用的主要有E l c e n t r o 波、T a f t 波、天津波等。在选择时应充分考虑地震动三要

素。即对于振幅,首先要考虑该被选择的地震与所考虑设防的地震有较好的相似性,有足够大的幅值及其峰值加速度应与设计烈度相当;其次,记录地点的地基土质条件与拟建场地土质条件宜相近。对于频谱选择,应该使所选地震波的卓越周期尽可能与拟建场地的地震周期一致,并且所选地震波的震中距尽可能与拟建场地震中距一致。对于强震持时,工程实践中确定持时的原则是应使地震记录最强烈部分包括在持时之内、若对结构仅进行弹性最大地震反应分析,持续时间可取短些;若对结构进行弹塑性最大地震反应分析,持时应取长些、一般可以取持时为结构基本周期的5～10倍。

人工合成场地地震波是根据随机振动理论产生的符合所需统计特征的地震波。对于人工地震波的

选择主要有三点要求,(1)人工地震波加速度峰值与设计烈度相当.(2)人工地震波的反应谱尽可能与规范规定的相应场地上的反应谱特征相近.(3)人工地震波的总持续时间和强震的平稳段持续时间应与按建筑场地的地震地质背景资料确定的相符。

研究表明,在充分考虑地震动三要素情况下,采用3～5条地震波可基本上保证时程分析结果的合理性。我国《建筑抗震设计规范》(G B 50011-2001)规定至少选择两条实际地震波和一条人工波[2]

。

2.4　恢复力模型

国内外已提出许多结构及构件的恢复力模型,其中最常用的恢复力模型主要有双线性模型(兰伯格-奥斯古德模型)和克拉夫、武腾青等提出的三线型刚度退化模型。根据构件的不同材料及不同力学特性,可以选用不同的更接近构件实际受力特性的计算模型。如对钢结构时程分析时可以采用最简单的双线型恢复力模型(见图3);对以受弯为主的钢筋混凝土构件则可以采用三线型刚度退化恢复力模型(见图4)

。

图3　双线型恢复力模型 图4　退化三线型恢复力模型

2.5　地震反应的数值模拟

在地震反应运动方程⑴中,地震地面运动加速度是一系列复杂的随机函数,同时在弹塑性反应中刚度

图5　逐步积分法地震动输入

矩阵与阻尼矩阵亦随时间变化,因此不可能求出解析解,而只能采取数值分析方法求解。所以常将式⑴转变成下面的增量方程式。

[M ]{■x ··

}+[C ]{■x ·

}+[K ]{■x }=-[M ]{■x ··g }

⑵

图6　弹性时程分析法程序设计框图

对于上式⑵需用逐步积分法求解。即将时间转化成一系列微小时间段■t ,在■t 时间内可采取一些假设,从而对增量方程⑵直接积分得出地震反应增量。以该步的终态值作为下一时间段的初始值,逐步积分,可求得结构在地震作用下振动反应的全过程。根据对于时间段内加速度变化所作的假设不同,就形成了不同的逐步积分法。常用的主要有线性加速度法、W i l s o n -θ法、N e w m a r k -β法等。3　时程分析法的M a t l a b 编程实现

建筑结构是一个庞大而又复杂的非线性结构,要想对其进行精确的动力反应分析,必须依靠计算机来处理。然而,不同的计算机语言直接影响着编程的繁琐和解决问题的快慢程度。与传统的F O R T R A N 和C 等高级语言相比,M A T L A B 语言用直观的、符合人们思维习惯的代码,代替了C 语言和F O R -T R A N 语言的冗长代码,使编程人员从繁琐的程序代码中解放出来,给用户带来了最直观最简洁的程序开发环境。将M A T L A B 语言编

写用于求解地震反应的弹性时程分析法的程序,其M A T L A B 程序框图,如图6所示。4　工程实例

工程概况,某三层钢筋混凝土结构,结构的特性参数为:第一至第三层质量m =2762k g ,2760k g ,2300k g ,第一至第三层刚度k=2.485×104

N /m ,

1.920×104N /m ,1.522×104

N /m 。地震波采用加速度峰值为200g a l 的E l -C e n t r o 波,采样周期为0.02s 。算模型及分析方法:采用上述层间剪切模型,利用弹性时程分析法,地震反应的求解过程采用W i l s o n -u 法,按照前述的程序设计框图,编制了M A T L A B 程序,计算结果见表1及图7。从计算耗时角度看,在保证计算结果可靠性的基础上,基于m a t l a b 语言求解结构地震反应的动力时程分析法所用时间比用其它语言编制的计算程序所用时间短,这主要是因为m a t l a b 语言本身具有强大的以矩阵计算为基础的数值计算功能,因而更适合求解结构特征值问题。能从根本上缩短计算时间,提高了时程分析法的计算效率。

表1　钢筋混凝土结构的自振周期振型号123自振频率

4.1041

10.4906〗

14.

9514

图7　顶层位移反应

5　结语

本文在简要分析时程分析法基本原理的基础上,阐述了用时程分析法求解结构地震反应的几个关键环节,然后利用简洁高效的M a t l a b 语言对时程分析法进行了编程实现,并以一三层钢筋混凝土结构为例运用编制的程序进行了计算。结果表明,基于m a t l a b 语言求解结构地震反应的时程分析法所用时间比用其它语言编制的计算程序所用时间短,这主要是因为m a t l a b 语言更适合求解结构特征值问题,因而从根本上缩短了计算时间,提高了时程分析法的计算效率。

参考文献

[1]刘晶波,杜修力.结构动力学[M ].北京:机械工业出版社,2006.

[2]胡聿贤著.地震工程学[M ].北京:地震出版社,2006.

层次分析法及matlab程序

层次分析法建模 层次分析法（AHP－Analytic Hierachy process）---- 多目标决策方法 70 年代由美国运筹学家T·L·Satty提出的，是一种定性与定量分析相结合的多目标决策分析方法论。吸收利用行为科学的特点，是将决策者的经验判断给予量化，对目标（因素）结构复杂而且缺乏必要的数据情况下，採用此方法较为实用，是一种系统科学中，常用的一种系统分析方法，因而成为系统分析的数学工具之一。 传统的常用的研究自然科学和社会科学的方法有： 机理分析方法：利用经典的数学工具分析观察的因果关系； 统计分析方法：利用大量观测数据寻求统计规律，用随机数学方法描述（自然现象、 社会现象）现象的规律。 基本内容：（1）多目标决策问题举例AHP建模方法 （2）AHP建模方法基本步骤 （3）AHP建模方法基本算法 （3）AHP建模方法理论算法应用的若干问题。 参考书： 1、姜启源，数学模型（第二版，第9章；第三版，第8章），高等教育出版社 2、程理民等，运筹学模型与方法教程，（第10章），清华大学出版社 3、《运筹学》编写组，运筹学（修订版），第11章，第7节，清华大学出版社 一、问题举例： A．大学毕业生就业选择问题 获得大学毕业学位的毕业生，“双向选择”时，用人单位与毕业生都有各自的选择标准和要求。就毕业生来说选择单位的标准和要求是多方面的，例如： ①能发挥自己的才干为国家作出较好贡献（即工作岗位适合发挥专长）； ②工作收入较好（待遇好）； ③生活环境好（大城市、气候等工作条件等）； ④单位名声好（声誉-Reputation）； ⑤工作环境好（人际关系和谐等） ⑥发展晋升（promote, promotion）机会多（如新单位或单位发展有后劲）等。 问题：现在有多个用人单位可供他选择，因此，他面临多种选择和决策，问题是他将如何作出决策和选择？——或者说他将用什么方法将可供选择的工作单位排序？

地震响应的反应谱法与时程分析比较 (1)

发电厂房墙体地震响应的反应谱法与时程分析比较 1问题描述 发电厂房墙体的基本模型如图1所示： 图1 发电厂墙体几何模型 基本要求：依据class 9_10.pdf的最后一页的作业建立ansys模型,考虑两个水平向地震波的共同作用(地震载荷按RG1.60标准谱缩放,谱值如下)，主要计算底部跨中单宽上的剪力与弯矩最大值,及顶部水平位移。要求详细的ansys反应谱法命令流与手算验证过程。以时程法结果进行比较。分析不同阻尼值(0.02,0.05,0.10)的影响。 RG1.60标准谱 (1g=9.81m/s2) (设计地震动值为0.1g) 频率谱值(g) 33 0.1 9 0.261 2.5 0.313 0.25 0.047 与RG1.60标准谱对应的两条人工波见文件rg160x.txt与rg160y.txt 2数值分析框图思路与理论简介 2.1理论简介 该问题主要牵涉到结构动力分析当中的时程分析和谱分析。时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。 2.2 分析框架： 时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。 谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。 3有限元模型与荷载说明 3.1 有限元模型 考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。然后，采用solid45

框架结构地震响应时程分析的计算模型

框架结构地震响应时程分析的计算模型 摘要：在结构进行地震响应时程分析时，必须首先确定结构的计算模型，以便确立结构的层间刚度。在地震作用下，结构计算模型是结构进行地震响应时分析的主体，由几何模型和物理模型两部分组成。其中几何模型反映了结构计算模型的几何构成，物理模型反映了材料或构件的力学性能。目前在工程上常用的计算模型主要有层间模型、杆系模型和杆系—层间模型。本文针对这三种模型进行全面的分析，并对它们的优缺点展开论述。 1前言 在求解结构在地震作用下的运动方程时，必须要计算结构的刚度矩阵[k]，而要计算结构的刚度矩阵[k]，就得确定结构的计算模型。因此，确定结构的计算模型是结构进行动力分析时必不可少的内容。对于多层框架结构，目前应用最广泛的模型是层间模型、杆系模型和杆系—层间模型。 2 层间模型 层间模型是在假定建筑各层楼板在其自身平面内刚度无穷大，水平地震作用下同层各竖向位移相同，以及建筑结构刚度中心和质量中心相重合，水平地震作用下没有绕竖轴扭转发生的基础上建立起来的。在这种模型中，将结构视为一根竖向杆，结构的质量集中于各楼层处，如图1（a）所示。 (a) (b) (c) (d) 图1 层间模型 （a）层间模型一般形式；（b）层间剪切模型；（c）层间弯曲模型；(d) 层间弯剪模型计算时，层间模型取各层为基本计算单元，采用层恢复力模型来表示地震作用过程中层刚度随层剪力的变化关系，而不考虑弹塑性阶段层刚度沿层高的变化。其几何模型相当于串联质点模型，物理模型的重要参数是层间刚度及其非线性变化规律。根据结构形式、构造特点以及结构侧向变形情况不同，层间模型又分为层间剪切模型、层间弯曲模型及层间弯剪模型，如图1（b）—（d）所示。其中，层间弯曲模型主要用于结构侧向变形以弯曲为主的剪力墙结构中。 而在进行框架结构动力分析时，常用的层间模型是层间剪切模型和层间弯剪模型。当框架横梁与柱的线刚度之比较大时，即“强梁弱住”型框架结构，在振动过程中各楼层始终保持水平，结构的变形表现为层间的错动，其侧向变形主要是层间剪切变形，那么应该采用层间剪切模型。 当框架梁对柱的约束相对较弱时，如一些高层框架，即“强柱弱梁”型结构，其侧向变形包含有层间弯曲和剪切两种成分，层间剪切模型已不能完全反映其变形特点，那么应该采用层间弯剪模型。 层间模型的优点在于自由度数较少，动力方程逐步积分所耗时也较少，但方法比较粗糙，计算精度较差，无法求出结构各杆件的时程反应，也不能确定结构各杆单元的内力和变形。因此，在工程实践中，层间模型主要是用于确定结构的层间剪力和层间侧移，以校核结构在地震作用下层间剪力是否超过层间极限承载力和检验结构在地震作用下的薄弱层位置。 3 杆系模型 杆系模型是较为精确的计算模型，它是在假定楼板在其自身平面内为绝对刚性的基础上建立起来的。这种模型将整个框架结构的梁柱构件离散为杆元，以结构的各杆件作为基本计算单元，将结构的质量集中于框架的各个节点，如图2所示。

地下结构地震破坏形式与抗震分析方法综述

地下结构地震破坏形式与抗震分析方法综述 摘要：随着人口的在激增以及经济的发展，人们的需求也开始狂飙式的增长。然而，城市的空间有限，地面空间已经被充分利用，人们的视线开始转为地下，地下结构的开发缓解了城市的地面压力。然而，由于地下结构的抗震技术的发展还并不成熟，在地震后，往往会造成地下结构的损坏甚至直接丧失继续工作的能力，给人们的财产安全带来威胁，影响人们的正常生活。因此在此文中对地下结构的震害形式以及近年来地下结构抗震分析的研究成果进行展示。以加深对地下结构震害的了解，并引起人们对地下结构抗震减震的重视。 关键词：地下结构抗震，震害形式，抗震分析，抗震减震 0 引言 地震是自然界自然界一种常见的自然灾害，地球上每年约发生500多万次地震，即每天要发生上万次地震。其中绝大多数太小或太远以至于人们感觉不到。真正能对人类造成严重危害的地震大约有一二十次，能造成特别严重灾害的地震大约有一两次。然而，这种地震不仅仅会给损害人们的财产安全，更有甚者会威胁到生命安全。 以往的抗震研究主要集中在地上建筑。认为地下结构受到的外界环境较少，各方向约束较多，刚度较大，且高度较小，加之过去地下结构的建设规模相对较少，地下结构受地震作用引起的结构的严重破坏的相关资料也较少，因此地下结构的工程抗震研究及设计长期未得到足够的重视。 1923年日本关东大地震（M8.2），震区内116座铁路隧道，有82座受到破坏；1952 年美国加州克恩郡地震（M7.6），造成南太平洋铁路的四座隧道损坏严重；1976年唐山地震（M7.8），唐山市给水系统完全瘫痪，秦京输油管道发生五处破坏；1978年日本伊豆尾岛地震（M7.0）震后出现了横贯隧道的断裂，隧道衬砌出现了一系列的破坏；特别是1995年日本阪神大地震（M7.2）中，神户市及阪神地区几座城市的供水系统和污水排放系统受到严重破坏，其中神户市供系统完全破坏，并基本丧失功能。神户市部分地铁车站和区间隧道受到不同程度的破坏，其中大开站最为严重，一半以上的中柱完全倒塌，导致顶板坍塌和上覆土层大量沉降，最大沉降量达2.5m。 地震对地下结构造成大规模破坏的同时，地震对地下结构的安全性构成的威胁也开始引起了人们的重视，地下结构工程抗震从业者在震后获取了大量的地震动作用在地下结构上产生的动力特性及影响结构动 力响应的影响因素等宝贵资料，对地下结构工程抗震减震领域的发展具有极大的推动作用。 近年来，关于地下结构的工程抗震分析方法的文献大量涌现。学者从不同角度对地下结构抗震进行阐述，并且有不少理论转化为工程技术，在工程实践中得到了论证。笔者试图综合前人的研究成果，在本文中简要介绍地下结构在地震作用下的破坏形式以及地下结构抗震分析方法，以便加深对地下结构工程抗震的了解，也可增加人们对地下结构工程抗震的重视程度。 1 地下结构震害 由于所处环境、约束情况等的差异，地下结构的破坏形式与结构破坏的影响因素与地上结构有很多不同之处。 1.1 地下结构震害形式 以下以日本阪神地震为主要对象，结合其他地震造成的震害，总结了地铁车站、地下管道、地下隧道的主要震害形式。

层次分析法matlab程序

disp('请输入判断矩阵A(n阶)'); A=input('A='); [n,n]=size(A); x=ones(n,100); y=ones(n,100); m=zeros(1,100); m(1)=max(x(:,1)); y(:,1)=x(:,1); x(:,2)=A*y(:,1); m(2)=max(x(:,2)); y(:,2)=x(:,2)/m(2); p=0.0001;i=2;k=abs(m(2)-m(1)); while k>p i=i+1; x(:,i)=A*y(:,i-1); m(i)=max(x(:,i)); y(:,i)=x(:,i)/m(i); k=abs(m(i)-m(i-1)); end a=sum(y(:,i)); w=y(:,i)/a; t=m(i); disp(w);disp(t); %以下是一致性检验 CI=(t-n)/(n-1);RI=[0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56

1.58 1.59]; CR=CI/RI(n); if CR<0.10 disp('此矩阵的一致性可以接受!'); disp('CI=');disp(CI); disp('CR=');disp(CR); end function AHPInit1(x,y) %层次分析的初始化 %默认只有两层x为准则数，y为方案数 %CToT为准则对目标生成的比较阵 %EigOfCri为准则层的特征向量 %EigOfOpt为选项层的特征向量 EigOfCri=zeros(x,1);%准则层的特征向量 EigOfOpt=zeros(y,x); dim=x;%维度 RI=[0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51];%RI标准%生成成对比较阵 for i=1:dim CToT(i,:)=input('请输入数据:'); end CToT %输出 pause, tempmatrix=zeros(x+1);

地震反应谱分析实例

结构地震反应谱分析实例 在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0

!进行模态求解 ANTYPE,MODAL MODOPT,LANB,30 SOLVE FINISH !进行谱分析 /SOLU ANTYPE,SPECTR SPOPT,SPRS,30,YES SVTYP,2 !加速度反应谱 SED,1,1 !X与Y向 FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167 FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000 SV,0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716 SV,0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.4853,0.4853,0.4853,0.4853 SV,0.05,0.2588,0.2167 SOLVE FINISH !进行模态求解(模态扩展) /SOLU ANTYPE,MODAL EXPASS,ON MXPAND,30,,,YES,0.005 SOLVE FINISH !进行谱分析(合并模态) /SOLU ANTYPE,SPECTR SRSS,0.15,disp SOLVE FINISH /POST1 SET,LIST !结果1 /INP,,mcom

大型地下结构三维地震响应特点研究

第43卷第3期2003年5月 大连理工大学学报 Jour nal of Dalian University of Technology Vol .43,No .3May 2003 文章编号:1000-8608(2003)03-0344-05 收稿日期:2002-04-01;　修回日期:2003-03-25. 基金项目:国家自然科学基金资助项目(50209002);辽宁省自然科学基金资助项目(20022130). 作者简介:陈健云*(1968-),男,副教授;林　皋(1929-),男,教授,博士生导师,中国科学院院士. 大型地下结构三维地震响应特点研究 陈健云*,　胡志强,　林　皋 (大连理工大学土木水利学院,辽宁大连　116024) 摘要:采用阻尼影响抽取法分析了地下结构无限围岩介质的动刚度特性,建立了岩石地下 结构抗震分析的实用相互作用分析时域模型,比较研究了地下结构-围岩动力相互作用分析中地震动输入机制、无限围岩动刚度及结构特性等各种主要因素对地下结构地震响应的影响程度.指出几种常用地下结构地震响应近似分析方法只在一定条件下适用,无限介质的阻尼特性对结构响应起着重要的作用. 关键词:地下洞室;地震反应分析;动刚度;优化;阻尼影响抽取法中图分类号:T U 35;TU 9;TV3 文献标识码:A 0　引　言 随着国民经济的发展,地下空间得到了越来越广泛的使用.然而近几年世界范围内发生了一 系列大地震,造成了巨大的灾难,不少地下结构遭受破坏.由于与围岩的相互作用,地下结构的动力特性十分复杂,其响应特点与地面结构有明显的差别.研究表明[1] ,对地下结构采用施加惯性力的地震响应分析,即使采用几倍于结构尺寸的地基离散模型,施加不同的边界条件对地震位移响应的影响可达10倍,应力差别达5～6倍. 目前各种实际地下结构的动力响应分析仍以各种近似方法为主.包括各种拟静力方法,如位 移响应法[2、3] ,地基影响参数通常根据简化假定采 用经验参数.动力近似分析通常将结构简化为二维问题处理[4],对于地下管线等结构形式具有一定的适用性.对于处于比较复杂地质、地形条件下的地下结构,或者形式较复杂的大型地下空间结构,要合理地反映地下结构的地震响应,则必须进行三维动力响应分析. 当前常用的地下结构三维地震分析方法,主要有在模型外边界施加各种人工透射边界解决能量向无限远处辐射[5]的波动分析方法;以地下结构为主体,围岩的作用通过相互作用力来求解的相互作用分析方法[6] ,通常采用有限元、边界元、 解析法或半解析法等耦合求解;以及在外边界施 加粘性阻尼器的惯性力方法.前两种方法属于较精确的数值方法,后一种方法则为近似方法. 由于围岩介质对结构的动力影响在时间与空 间都是耦合的,较精确的地下结构地震响应分析具有一定难度,时域求解复杂且求解代价很大. 本文采用相互作用分析方法,结合溪洛渡超 大型地下洞室群的地震响应分析,研究动力相互 作用运动方程中各主要因素对地下结构地震响应的影响程度,为地下结构的简化分析提供依据. 1　地下结构地震响应的相互作用分 析方法 地下结构的相互作用分析主要采用各种耦合 方法,如有限元与边界元的耦合分析.本文则采 用阻尼影响抽取法得到地基刚度与有限元进行耦合分析. 1.1　阻尼影响抽取法的基本概念 [7] 将无限地基截取有限区域,其刚度阵为S t (X )=K -X 2 M (1) 式中:K 和M 分别为有限域的刚度阵与质量阵. 引入量纲一的频率a 0=X ?r 0/c s 及刚度阵K 与质量阵M ,则式(1)可表达为 　S t (X )=Gr s -2 0(K -a 20M )=Gr s -2 0S (a 0) (2)

层次分析法实现代码(MATLAB)

%% AHP weight calculation %%data input clc clear all A =[1 3 5 7 9 5;1/3 1 3 9 3 3;1/5 1/3 1 3 3 1/3;1/7 1/9 1/3 1 5 1/3;1/9 1/3 1/3 1/5 1 1/3;1/5 1/3 1 3 3 1]; %%Consistency calculation and weight vector calculation [n,n] = size(A); [v,d] = eig(A); r = d(1,1); CI = (r-n)/(n-1); RI = [0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.52 1.54 1.56 1.58 1.59]; CR = CI/RI(n); if CR<0.10 CR_Result = 'pass'; else CR_Result = 'no pass'; end % % Weight vector calculation w = v(:,1)/sum(v(:,1));

w = w'; % % output disp('The judgment matrix weight vector calculation report：'); disp('coincidence indicator：');disp(num2str(CI)); disp('Consistency ratio：');disp(num2str(CR)); disp(' Consistency test results：');disp(CR_Result); disp('eigenvalue:');disp(num2str(r)); disp('weight vector:');disp(num2str(w));

基于Matlab的层次分析法及其运用浅析

基于Matlab的层次分析法及其运用浅析 本文通过使用Matlab软件进行编程,在满足同一层次中各指标对所有的下级指标均产生影响的假定条件下,实现了层次分析法的分析运算。本程序允许用户自由设定指标层次结构内的层次数以及各层次内的指标数,通过程序的循环,用户只需输入判断矩阵的部分数据,程序可依据层次分析法的计算流程进行计算并作出判断。本程序可以方便地处理层次分析法下较大的运算量,解决层次分析法的效率问题,提高计算机辅助决策的时效性。 标签：Matlab层次分析法判断矩阵决策 在当前信息化、全球化的大背景下,传统的手工计算已不能满足人们高效率、高准确度的决策需求。因此计算机辅助决策当仁不让地成为了管理决策的新工具、新方法。基于此,本文在充分发挥计算机强大运算功能的基础上,选用美国MathWorks公司的集成数学建模環境Matlab R2009a作为开发平台,使用M语言进行编程,对计算机辅助决策在层次分析法中的运用进行讨论。试图通过程序实现层次分析法在计算机系统上的运用,为管理决策探索出新的道路。 1 层次分析法的计算流程 根据层次分析法的相关理论,层次分析法的基本思想是将复杂的决策问题进行分解,得到若干个下层指标,再对下层指标进行分解,得到若干个再下层指标,如此建立层次结构模型,然后根据结构模型构造判断矩阵,进行单排序,最后,求出各指标对应的权重系数,进行层次总排序。 1.1 构造层次结构模型在进行层次分析法的分析时,最主要的步骤是建立指标的层次结构模型,根据结构模型构造判断矩阵,只有判断矩阵通过了一致性检验后,方可进行分析和计算。其中,结构模型可以设计成三个层次,最高层为目标层,是决策的目的和要解决的问题,中间层为决策需考虑的因素,是决策的准则,最低层则是决策时的备选方案。一般来讲,准则层中各个指标的下级指标数没有限制,但在本文中设计的程序尚且只能在各指标具有相同数量的下级指标的假定下,完成层次分析法的分析,故本文后文选取的案例也满足这一假定。 1.2 建立判断矩阵判断矩阵是表示本层所有因素针对上一层某一个因素的相对重要性的比较给判断矩阵的要素赋值时,常采用九级标度法(即用数字1到9及其倒数表示指标间的相对重要程度),具体标度方法如表1所示。 1.3 检验判断矩阵的一致性由于多阶判断的复杂性,往往使得判断矩阵中某些数值具有前后矛盾的可能性,即各判断矩阵并不能保证完全协调一致。当判断矩阵不能保证具有完全一致性时,相应判断矩阵的特征根也将发生变化,于是就可以用判断矩阵特征根的变化来检验判断的一致性程度。在层次分析法中,令判断矩阵最大的特征值为λmax,阶数为n,则判断矩阵的一致性检验的指标记为:

反应谱与时程理论对比

反应谱是在给定的地震加速度作用期间内，单质点体系的最大位移反应、速度反应和加速度反应随质点自振周期变化的曲线。用作计算在地震作用下结构的内力和变形。更直观的定义为：一组具有相同阻尼、不同自振周期的单质点体系，在某一地震动时程作用下的最大反应，为该地震动的反应谱。 反应谱理论考虑了结构动力特性与地震动特性之间的动力关系，通过反应谱来计算由结构动力特性(自振周期、振型和阻尼)所产生的共振效应，但其计算公式仍保留了早期静 力理论的形式。地震时结构所受的最大水平基底剪力，即总水平地震作用为: FEK= αG 其中α为地震影响系数，即单质点弹性体系在地震时最大反应加速度。另一方面地震影响系数也可视为作用在质点上的地震作用与结构重力荷载代表值之比。 目前，反应谱分析法比较成熟，一些主要国家的抗震规范均将它作为基本设计方法。不过，它主要适合用于规则结构。对于不规则结构以及高层建筑，各国规范多要求采用时程分析法进行补充计算。 地震作用反应谱分析本质上是一种拟动力分析，它首先使用动力法计算质点地震响应，并使用统计的方法形成反应谱曲线，然后使用静力法进行结构分析。但它并不是结构真实的动力响应分析，只是对于结构动力响应最大值进行估算的近似方法，在线弹性范围内，反应谱分析法被认为是高效而且合理的方法。反应谱分为加速度反应谱、速度反应谱和位移反应谱。基于不同周期结构相应峰值的大小，我们可以绘制结构速度及加速度的反应谱曲线。一般情况下，随着周期的延长，位移反应谱为上升曲线，速度反应谱为平直曲线，加速度反应谱为下降曲线，目前结构设计主要依据加速度反应谱。 加速度反应谱在短周期部分为快速上升曲线，并且在结构周期与场地特征周期接近时出现峰值，后面更大范围为逐渐下降阶段。峰值出现的时间与对应的结构周期和场地特征周期有关。一般来说结构自振周期的延长，地震作用将减小。当结构自振周期接近场地特征周期时，地震作用最大。 反应谱分析方法需要先求解一个方向地震作用响应，再基于三个正交方向的分量考虑结构总响应，即基于振型组合求解一个方向的地震响应，再基于方向组合求解结构总响应。 振型组合方法有SRSS法，CQC法。 1.SRSS法 SRSS法是平方和平方根法，这种方法假定所有最大模态值在统计上都是相互独立的，通过求各参与阵型的平方和平方根来进行组合。该法不考虑各振型间的藕联作用，实际上结构模态都是相互关联的，不可避免的存在藕联效应，对那些相邻周期几乎相等的结构，或者不规则结构不适用此法。《抗规》GB50011-2010规定的SRSS法为如下所示：

地震工程学心得体会

精心整理《地震工程学》课程总结? 1.对所学内容的综述? 1.1结构地震反应分析的方法? 结构地震反应分析的方法很多，下面主要介绍反应谱理论和时程反应分析法? 绍。 也并不是一次地震动作用下的反应谱，而是不同地震反应的包线。 1.1.2?? 时程分析法? 时程分析法又称作动态分析法。它是将地震波段按时段进行数值化后，输入结构体系的振动微分方程，采用逐步积分法进行结构弹塑性动力反应分析，计算出结构在整个强震时域中的振动状态过程，给出各个时刻各杆件的内力和变形以及各杆

件出现塑性铰的顺序。? 时程分析法计算地震反应需要输入地震动参数，该参数具有概率含义的加速度时程曲线、结构和构件的动力模型考虑了结构的非线性恢复力特性，更接近实际情况，因而时程分析方法具有很多优点。它全面地考虑了强震三要素；比较确切地、具体地和细致地给出了结构弹塑性地震反应。? 1.1.3地震信号频域分析? ???? X(f）， 1.2? 1.2.1 （1) ??(2 （3 ?（4 性和有效性；? ?? （5）验证抗震理论、结构地震反应分析方法、结构振动控制算法等的可靠性和适用性。? 1.2.2? 结构抗震试验的实施程序? ??

（1）确定研究目标和试验方法，含试验目的、试验设备和试件的采用、需要测量的物理量等；? ?? （2）荷载施加，含与试验设备相关的荷载施加方式和加载规则等；? ?（3）测点布置和数据采集，含各类传感器和数采设备的采用、测点数量的选择；? ??（4）数据分析，含测试数据的常规处理和特殊分析。? （1 ? （2 ????旨在 （3 ?? 入下结构或构件的地震反应，研究和验证结构地震破坏机理、破坏特征、抗震能力和抗震薄弱环节。 ?（4）振动台试验? ?????振动台试验是利用振动台装置进行的结构强迫振动试验，是地震工程研究中最重要的实验手段之一。?

ANSYS地震反应谱SRSS分析共24页

ANSYS地震反应谱SRSS分析 我在ANSYS中作地震分解反应谱分析，一次X方向，一次Y 方向，他们要求是独立互不干扰的，可是采用直进行一次模态分析的话，他生成的*.mcom文件好像是包含了前面的计算 结果，命令流如下： !进入PREP7并建模 /PREP7 B=15 !基本尺寸 A1=1000 !第一个面积 A2=1000 !第二个面积 A3=1000 !第三个面积 ET,1,beam4 !二维杆单元 R,1,0.25,0.0052,0.0052,0.5,0.5 !以参数形式的实参 MP,EX,1,2.0E11 !杨氏模量 mp,PRXY,1,,0.3 mp,dens,1,7.8e3 N,1,-B,0,0 !定义结点 N,2,0,0,0 N,3,-B,0,b

N,4,0,0,b N,5,-B,0,2*b N,6,0,0,2*b N,7,-B,0,3*b N,8,0,0,3*b E,1,3 !定义单元 E,2,4 E,3,5 E,4,6 E,3,4 E,5,6 e,5,7 e,6,8 e,7,8 D,1,ALL,0,,2 FINISH ! !进入求解器，定义载荷和求解 /SOLU D,1,ALL,0,,2 !结点UX=UY=0

sfbeam,1,1,PRES,100000, sfbeam,3,1,PRES,100000, sfbeam,7,1,PRES,100000, SOLVE FINISH allsel NMODE=10 /SOL !* ANTYPE,2 !* MSAVE,0 !* MODOPT,LANB,NMODE EQSLV,SPAR MXPAND,NMODE , , ,1 LUMPM,0 PSTRES,0 !* MODOPT,LANB,NMODE ,0,0, ,OFF

结构地震反应谱分析实例

在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X 与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0<T<=0.04 秒 0.4853*(0.10/T)^(-0.686) 0.04<T<=0.1 秒 0.4853 0.1<T<=1.2 秒 0.4853*(1.2/T)^1.5 1.2<T<=4 秒 以下是命令流程序 ---------------------------------------------------------------------------------------------------- /filname,SPEC,1 /PREP7 !定义单元类型及材料特性 ET,1,45 MP,EX,1,2.8E10 MP,DENS,1,2.4E3 MP,NUXY,1,0.18 !建立模型 BLOCK,0,1,0,1,0,5 !网格剖分 ESIZE,0.5 VMESH,all /VIEW,,-0.3,-1,1 EPLOT FINISH /SOLU !施加底部约束 ASEL,,LOC,Z,0 DA,ALL,ALL ALLSEL !施加自重荷载 ACEL,0,0,10 !进行模态求解

ANTYPE,MODAL MODOPT,LANB,30 SOLVE FINISH !进行谱分析 /SOLU ANTYPE,SPECTR SPOPT,SPRS,30,YES SVTYP,2 !加速度反应谱 SED,1,1 !X与Y向 FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167 FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000 SV,0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716 SV,0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.4853,0.4853,0.4853,0.4853 SV,0.05,0.2588,0.2167 SOLVE FINISH !进行模态求解(模态扩展) /SOLU ANTYPE,MODAL EXPASS,ON MXPAND,30,,,YES,0.005 SOLVE FINISH !进行谱分析(合并模态) /SOLU ANTYPE,SPECTR SRSS,0.15,disp SOLVE FINISH /POST1 SET,LIST !结果1 /INP,,mcom lcwrite,11

层次分析法matlab程序举例

层次分析法程序举例： A=[1 1/7 1/5 2 4 1/3;7 1 3 5 5 3;5 1/3 1 5 5 3;1/2 1/3 1/5 1 2 1/3;1/4 1/5 1/5 1/2 1 1/5;3 1/3 1/3 3 5 1]; [v,d]=eig(A); eigenvalue=diag(d); lamda=max(eigenvalue); cil=(lamda-6)/5; crl=cil/1.26; w1=v(:,1)/sum(v(:,1)) 挑选合适的工作。经双方恳谈，已有三个单位表示愿意录用某毕业生。该生根据已有信息建立了一个层次结构模型，如下图所示。 程序： A=[1 1/7 1/5 2 4 1/3;7 1 3 5 5 3;5 1/3 1 5 5 3;1/2 1/3 1/5 1 2 1/3;1/4 1/5 1/5 1/2 1 1/5;3 1/3 1/3 3 5 1]; [v,d]=eig(A); eigenvalue=diag(d); lamda=max(eigenvalue); ci=(lamda-6)/5

cr=ci/1.26 w1=v(:,1)/sum(v(:,1)) B1=[1 1/4 1/2;4 1 3;2 1/3 1]; [v,d]=eig(B1); eigenvalue=diag(d); lamda=max(eigenvalue); cil1=(lamda-3)/2 cr1=cil1/0.52 b1w=v(:,1)/sum(v(:,1)) B2=[1 1/4 1/5;4 1 1/2;5 2 1]; [v,d]=eig(B2); eigenvalue=diag(d); lamda=max(eigenvalue); cil2=(lamda-3)/2 cr2=cil2/0.52 b2w=v(:,1)/sum(v(:,1)) B3=[1 1/2 2; 2 1 3;1/2 1/3 1]; [v,d]=eig(B3); eigenvalue=diag(d);

三 设计地震动反应谱确定的规范方法

三设计地震动反应谱确定的规范方法 设计地震动是通过对地震环境和场地环境的分析判断和分类方法确定。工程勘察单位至少提供： 设计基本地震加速度和设计特征周期 场地环境：覆盖层厚度、剪切波速、土层钻孔资料 1.设计基本地震加速度和设计特征周期 根据场地在中国地震动参数区划图上的位置判断确定。

土层剪切波速的测量应符合下列要求: 1 在场地初步勘察阶段对大面积的同一地质单元测量土层剪切波速的钻孔数量不宜少于3。 2 在场地详细勘察阶段对单幢建筑测量土层剪切波速的钻孔数量不宜少于2 个数据变化较大时可适量增加对小区中处于同一地质单元的密集高层建筑群测量土层剪切波速的钻孔数量可适量减少但每幢高层建筑下不得少于一个。 3 对丁类建筑及层数不超过10 层且高度不超过30m 的丙类建筑当无实测剪切波速时可根据岩土名称和性状按表 4.1.3 划分土的类型再利用当地经验在下表的剪切波速范围内估计各土层的剪切波速.

建筑场地覆盖层厚度的确定应符合下列要求: 1 一般情况下应按地面至剪切波速大于500m/s 的土层顶面的距离确定（且其下卧层沿途的剪切波速均不小于500m/s）。 2 当地面5m 以下存在剪切波速大于（其上部各土层）相邻上层土剪切波速2.5 倍的土层且其下卧岩土的剪切波速均不小于400m/s 时可按地面至该土层顶面的距离确定 3 剪切波速大于500m/s 的孤石、透镜体应视同周围土层 4.土层中的火山岩硬夹层应视为刚体其厚度应从覆盖土层中扣除

例题：某类建筑场地位于7度烈度区，设计地震分组为第一组，设计基本地震加速度为0.1g，建筑结构自振周期T=1.4s，阻尼比为0.08，该场地在建筑多遇地震条件下地震影响系数a为多少。 同一个场地上甲乙两座建筑物的结构自震周期分别为T甲=0.25sT乙=0.60s，一建筑场地类别为Ⅱ类，设计地震分组为第一组，若两座建筑的阻尼比都取0.05，问在抗震验算时甲、乙两座建筑的地震影响系数之比最接近下列那个选项。 A 1.6 B 1.2 C 0.6 D 条件不足无法计算 例题：吉林省松原市某民用建筑场地地质资料如下： （1）0-5m粉土，=150 =180m/s (2) 5-12m中砂土=200 =240m/s （3）12-24m粗砂土=230 =310m/s (4) 24-45m硬塑粘土=260 =300m/s （5）45-60m泥岩=500 =520m/s 建筑物采用浅基础，埋深2m，地下水位2.0m，阻尼比为0.05，自震周期为1.8s该建筑进行抗震设计时 （1）进行第一阶段设计时，地震影响系数应取多少 （2）进行第二阶段设计时，地震影响系数应取多少 例题：吉林省松原市某民用建筑场地地质资料如下： （1）0-5m粉土，=150 =180m/s (2) 5-12m中砂土=200 =240m/s

结构抗震课后习题答案解析

《建筑结构抗震设计》课后习题解答建筑结构抗震设计》第 1 章绪论 1、震级和烈度有什么区别和联系？震级是表示地震大小的一种度量，只跟地震释放能量的多少有关，而烈度则表示某一区域的地表和建筑物受一次地震影响的平均强烈的程度。烈度不仅跟震级有关，同时还跟震源深度、距离震中的远近以及地震波通过的介质条件等多种因素有关。一次地震只有一个震级，但不同的地点有不同的烈度。 2.如何考虑不同类型建筑的抗震设防？规范将建筑物按其用途分为四类：甲类（特殊设防类）、乙类（重点设防类）、丙类（标准设防类）、丁类（适度设防类）。 1 ）标准设防类，应按本地区抗震设防烈度确定其抗震措施和地震作用，达到在遭遇高于当地抗震设防烈度的预估罕遇地震影响时不致倒塌或发生危及生命安全的严重破坏的抗震设防目标。 2 ）重点设防类，应按高于本地区抗震设防烈度一度的要求加强其抗震措施；但抗震设防烈度为 9 度时应按比 9 度更高的要求采取抗震措施；地基基础的抗震措施，应符合有关规定。同时，应按本地区抗震设防烈度确定其地震作用。 3 ）特殊设防类，应按高于本地区抗震设防烈度提高一度的要求加强其抗震措施；但抗震设防烈度为 9 度时应按比 9 度更高的要求采取抗震措施。同时，应按批准的地震安全性评价的结果且高于本地区抗震设防烈度的要求确定其地震作用。 4 ）适度设防类，允许比本地区抗震设防烈度的要求适当降低其抗震措施，但抗震设防烈度为 6 度时不应降低。一般情况下，仍应按本地区抗震设防烈度确定其地震作用。 3.怎样理解小震、中震与大震？ 小震就是发生机会较多的地震，50 年年限，被超越概率为63.2%；中震，10%；大震是罕遇的地震，2%。 4、概念设计、抗震计算、构造措施三者之间的关系？ 建筑抗震设计包括三个层次：概念设计、抗震计算、构造措施。概念设计在总体上把握抗震设计的基本原则；抗震计算为建筑抗震设计提供定量手段；构造措施则可以在保证结构整体性、加强局部薄弱环节等意义上保证抗震计算结果的有效性。他们是一个不可割裂的整体。

地层地震反应对地下结构的影响

地层地震反应对地下结构的影响 隧道二班谭坤（07011227） 地震对地下工程影响的一般规律 地震对地下工程的影响规律总体上有以下的特点: 1) 地下结构的振动变形受周围地基土壤的约束作用显著,结构的动力反应 一般不明显表现出自振特性的影响。 2) 地下结构的振动形态受地震波入射方向的影响很大,地震波的入射方向 发生不大的变化,地下结构各点的变形和应力可以发生较大的变化,相位差别也 十分明显。但主要应变一般与地震加速度大小的联系不很明显,随埋深发生的变化也不很明显。 3) 地下结构地基的相互作用都对它的动力反应产生重要影响,对结构动力 反应起主要作用的因素是地基的运动特性,一般来说,结构形状的改变对动力反 应的影响相对较小,只引起量的变化。而地下结构的存在对周围地基震动的影响一般很小(指地下结构的尺寸相对于地震波长的比例较小的情况) 。 岩体隧道震害的形式主要有裂纹、剥落、底部隆起或倾斜,破坏程度主要取决于地震作用力方向及现场地质条件,一般发生于存在破碎带的地层中。 对于土体隧道,土体对地震的响应要明显强于岩体,所以隧道破坏的可能性 也更大。又由于土体隧道多用于城市地铁,车站较多,整体结构形式不均一,容易产生应力集中,使破坏多集中在车站上。 1) 并行隧道距离越小, 其地震内力反应越大, 当距离小于隧道断面外径D , 尤其是小于0. 5D 时, 抗震设计应给予足够的重视; 2) 地震引起的地基变形是影响盾构隧道地震反应的决定性因素, 因而在抗震设计时需要合理考虑盾构隧道应承受的地基变形, 因此相对于地震系数法, 反应位移法的设计思想更为合理; 3) 相对于目前广泛采用的设计基本地震加速度, 对地铁区间隧道等地下结构进行抗震分析及设计时采用地面峰值相对位移作为设计地震动参数更为合理。 上述结论是基于三类建筑场地条件得出的, 可供地铁盾构区间隧道等地下 结构抗震设计参考。对于其他场地条件, 还有待进一步研究。

层次分析法计算权重在matlab中的实现

信息系统分析与设计作业 层次分析法确定绩效评价权重在matlab中的实现 小组成员：孙高茹、王靖、李春梅、郭荣1 程序简要概述 编写程序一步实现评价指标特征值lam、特征向量w以及一致性比率CR的求解。 具体的操作步骤是：首先构造评价指标，用专家评定法对指标两两打分，构建比较矩阵，继而运用编写程序实现层次分析法在MATLAB中的应用。 通过编写MATLAB程序一步实现问题求解，可以简化权重计算方法与步骤，减少工作量，从而提高人力资源管理中绩效考核的科学化电算化。 2 程序在matlab中实现的具体步骤 function [w,lam,CR] = ccfx(A) %A为成对比较矩阵，返回值w为近似特征向量 % lam为近似最大特征值λmax，CR为一致性比率 n=length(A(:,1)); a=sum(A); B=A %用B代替A做计算 for j=1:n %将A的列向量归一化 B(:,j)=B(:,j)./a(j); end s=B(:,1); for j=2:n s=s+B(:,j); end c=sum(s);%计算近似最大特征值λmax w=s./c; d=A*w lam=1/n*sum((d./w)); CI=(lam-n)/(n-1);%一致性指标 RI=[0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45,1.49,1.51];%RI为随机一致

性指标 CR=CI/RI(n);%求一致性比率 if CR>0.1 disp('没有通过一致性检验'); else disp('通过一致性检验'); end end 3 案例应用 我们拟构建公司员工绩效评价分析权重，完整操作步骤如下： 3.1构建的评价指标体系 我们将影响员工绩效评定的指标因素分为：打卡、业绩、创新、态度与品德。 3.2专家打分，构建两两比较矩阵 A = 1.0000 0.5000 3.0000 4.0000 2.0000 1.0000 5.0000 3.0000 0.3333 0.2000 1.0000 2.0000 0.2500 0.3333 0.5000 1.0000 3.3在MATLAB中运用编写好的程序实现 直接在MATLAB命令窗口中输入 [w,lam,CR]=ccfx(A) 继而直接得出 d = 1.3035 2.0000 0.5145 0.3926 w = 0.3102 0.4691 0.1242 0.0966 lam =4.1687