结构地震反应谱分析实例
地震响应的反应谱法与时程分析比较
地震响应的反应谱法与时程分析比较(总13页)-CAL-FENGHAI.-(YICAI)-Company One1-CAL-本页仅作为文档封面,使用请直接删除发电厂房墙体地震响应的反应谱法与时程分析比较1问题描述发电厂房墙体的基本模型如图1所示:图1 发电厂墙体几何模型基本要求:依据class 的最后一页的作业建立ansys模型,考虑两个水平向地震波的共同作用(地震载荷按标准谱缩放,谱值如下),主要计算底部跨中单宽上的剪力与弯矩最大值,及顶部水平位移。
要求详细的ansys反应谱法命令流与手算验证过程。
以时程法结果进行比较。
分析不同阻尼值,,的影响。
标准谱 (1g=s2) (设计地震动值为频率谱值(g)339与标准谱对应的两条人工波见文件与2数值分析框图思路与理论简介理论简介该问题主要牵涉到结构动力分析当中的时程分析和谱分析。
时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。
谱分析是模态分析的扩展,是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。
分析框架:时程分析:在X和Z两个水平方向地震波作用下,提取底部跨中单宽上的剪力、弯矩值和顶部水平位移,并求出最大响应。
谱分析:先做模态分析,再求谱解,由于X和Z两个方向的单点谱激励,因此需进行两次谱分析,分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。
3有限元模型与荷载说明有限元模型考虑结构的几何特性建立有限元模型,首先建立平面几何模型,并将模型进行合理的切割,采用plane42单元,使用映射划分网格的方法生产平面单元(XOY平面)。
然后,采用solid45单元,设置拖拉方向的单元尺寸并清楚初始平面单元plane42,将平面单元进行拖拉,最后生成发电厂墙体的有限元立体几何模型。
单元总数为6060个,总节点数为8174个,有限元模型如图2所示:图2 发电厂墙体有限元模型荷载说明时程分析:首先计算结构的前两阶自振频率,分别为126.10008.2867f f ==,。
某连续梁桥地震反应谱分析
纵向剪力
0 0 0
横 向剪力 F
10 12 3 . 249 6 1. 326 8
纵 向剪 力
O 0 O
用条件下 , 制动 墩主要 承担桥 梁上部结 构所 引起 的 向剪力 ; 3 由计 算结 果可以看 到 : 向剪力各墩 比较平 均 , ) F 这是 由于 桥梁支座为单 向滑动 的 支座 , 拟 时要约 束 住 F 模 向的 自由度 。 因此各墩共 同承担这个方 向的剪力 ;
本文采用 “ m法 ” 模拟桩土相互作用 , 此法是一 种常用 的模拟
f 模型 『
结构桩一土相互之间作用力
输入地震力 模式
桩基静力 的方法 , 在设计单位 经常采用 , 在建模 时 , 据本桥 结构 依
的特点 , 采用三 维梁 单元 离散桥 梁 的主题 结构 , 括桥 墩 、 梁、 包 主
向的组合 , 不模 拟地震 的竖 向作 用力 , 分别对 两种模 型地 震反 应
图 2 箱 梁横 断面 I 单位 :i ) tl 1
进行模拟 , 见表 1 。
表 1 反 应 谱 分 析 模 型 和 输 入地 震 力模 式
2 桩一 土 结构 相 互作 用分 析方 法和 结构 模型 建 立
1. 8 6 21. 1
1 . 44
32. 1 49. 7
58 .
1 2 0. l 1 6.
86 .
2 9 1. 3 6 3.
模拟 , 向活动支 座具有 竖 向转 动和单 一滑 移性 能 ; 单 固定 支座 仅 具有竖 向转 动性能 。本文采用两种模 型来进行 模拟 , 不考 虑桩 土 相互作用 的模 型 A, 墩第固结 , 图 3所示 , 一种采 用需要 考 将 如 另
大底盘双塔连体结构的地震反应谱分析
第一作者简介 : 佞 (9 l , , 肖 18 一) 男 昆明理工大学硕士研究生 。
1 7期
肖 佞, : 等 大底 盘双塔连体结构 的地震反应谱分析
4 0 03
7 0 24 0
… … … … … — 一
0 0 3 5 0 5 0 04 5 0 0 0 0 0
摘
要
大底 盘双 塔结构是 指将底部 几层公共空间设置为大底 盘, 在楼层上部用 两个 塔楼作 为主体 的结构 , 而在塔楼 间的某
些楼层用连体连接起来 , 就成为 了大底盘双塔楼连体 高层建筑 结构。在近 些年 , 此类结构 在 国 内外 的各 种工程 中均有应 用。 运用 A S S 限元程序采用振 型分 解反应谱法对大底 盘双塔连体 结构进行 分析 , NY 有 得到不 同的抗震性能结果。 关键词 大底盘 双塔 连 体结构 连接方式 A 地 震响应 中图法分类号 T 32 1 U 5 .; 文献 标志码
4
鬻
雾
霪
中
睾 拿
一
丰
+ ÷ 牛
串 孛
串 +
串
誓一
毒一 誊 簪
榔
孛 案
拳 一
案
一
魂 霪
礴 莲 。} l ; l 6 , j
图 1 大 底 盘 平 面 布 置 图
一 一 一
地震 反 应是 指结 构 由地震 引起 的振 动 , 它包 括 地 震在 结构 中 引起 的 内 力 、 形 和 位 移 , 般 都 是 变 一 通 过结 构 动力学 的分析 方法 来进 行求 解 。
~ 一
} ¨
i 一 梆
—
…
一
0 i
考虑土_结构相互作用体系地震反应谱分析_徐源
( m·s - 1 ) 357. 2 160. 0 131. 0 60. 0
场地 类别
Ⅱ Ⅱ Ⅲ Ⅲ
3 计算结果
3. 1 自振周期的计算 高度小于 40 m 的以剪切变形为主的规则框
架结构,第一振型为主要位移振型,表 2 列出了各 工况的基本周期和竖向一阶周期的计算结果.
表 2 考虑和不考虑相互作用的周期对比 Table 2 Comparison of cycle considering
58. 32 151. 632
0. 3
6
180. 00 432. 000
0. 2
6
352. 80 846. 720
0. 2
18
6. 48 18. 144
0. 4
剪切波速 等效剪切波速
( m·s - 1 )
250 300 360 420 450 480 160 60 90 160 180 300 420 60
Structure Interaction
XU Yuan,CONG Long,HU Ying-peng,LIANG Yan
( School of Resources & Civil Engineering,Northeastern University,Shenyang 110819,China. Corresponding author: XU Yuan,E-mail: xuyuanneu@ gmail. com)
土层 编号
1 2 3 4 5 6 1 1 2 3 4 5 6 1
表 1 地基土参数 Table 1 Foundation soil parameters
土层厚度 剪切模量 弹性模量
m
MPa
MPa
ABAQUS时程分析法计算地震反应的简单实例
ABAQUS时程分析法计算地震反应的简单实例ABAQUS时程分析法计算地震反应的简单实例(在原反应谱模型上修改)问题描述:悬臂柱高12m,工字型截面(图1),密度7800kg/m3,EX=2、1e11Pa,泊松比0、3,所有振型的阻尼比为2%,在3m高处有一集中质量160kg,在6m、9m、12m处分别有120kg的集中质量。
反应谱按7度多遇地震,取地震影响系数为0、08,第一组,III类场地,卓越周期Tg=0、45s。
图1 计算对象第一部分:反应谱法几点说明:λ本例建模过程使用CAE;λ添加反应谱必须在inp中加关键词实现,CAE不支持反应谱;λ *Spectrum不可以在keyword editor中添加,keyword editor不支持此关键词读入。
λ ABAQUS的反应谱法计算过程以及后处理要比ANSYS方便的多。
操作过程为:(1)打开ABAQUS/CAE,点击create model database。
(2)进入Part模块,点击create part,命名为column,3D、deformation、wire。
continue(3) Create lines,在分别输入0,0回车;0,3回车;0,6回车;0,9回车;0,12回车。
(4)进入property模块,create material,name:steel,general-->>density,mass density:7800mechanical-->>elasticity-->>elastic,young‘s modulus:2、1e11,poisson’s ratio:0、3、(5) Create section,name:Section-1,category:beam,type:beam,Continuecreate profile, name:Profile-1, shape:I,按图1尺寸输入界面尺寸,ok。
第三章2 工程结构地震反应分析与抗震验算.ppt
h 1 ---直线下降段的斜率调整系数;按下式确定
h1 = 0.02 + (0.05 - z ) / 8 当h1 < 0时,取h1 = 0
h2 - -阻尼调整系数,h2 < 0.55时,取h2 = 0.55
h2
=1+
0.05 - z 0.06 +1.7z
Tg : 特征周期,见表3.2
max:水平地震系数的最大值 α max = kβ max ,β max= 2.25
结构在地震持续过程中经受的最大地震作用为
F
=
F (t ) max
= m &x&(t) + &x&g (t) max
= mSa
= mg Sa
&x&g (t) max = Gk = G
&x&g (t) max
g
G ---集中于质点处的重力荷载代表值;
g ---重力加速度
= Sa
&x&g (t) max
地震特征周期分组的特征周期值(s)
场地类别
Ⅰ
Ⅱ
Ⅲ
Ⅳ
第一组 0.25
0.35
0.45 0.65
查表确定 Tg Tg = 0.3
第二组 0.30
0.40
第三组 0.35
0.45
0.55 0.75 0.65 0.90
例:单层单跨框架。屋盖刚度为无穷大,质量集中于屋 盖处。已知设防烈度为8度,设计地震分组为二组,Ⅰ类 场地;屋盖处的重力荷载代表值G=700kN,框架柱线刚 度 ic = EIc / h = 2.6104 kN m ,阻尼比为0.05。试求该结构多 遇地震时的水平地震作用。
ABAQUS地震反应谱分析
ABAQUS反应谱法计算地震反应的简单实例Fan.hj2010年4月4日清明小长假,琢磨了下ABAQUS如何进行地震反应谱计算。
现通过一小算例说明。
问题描述:(本例的问题引用《有限元法及其应用》一书中陆新征博士ANSYS算例的问题)悬臂柱高12m,工字型截面(图1),密度7800kg/m3,EX=2.1e11Pa,泊松比0.3,所有振型的阻尼比为2%,在3m高处有一集中质量160kg,在6m、9m、12m处分别有120kg的集中质量。
反应谱按7度多遇地震,取地震影响系数为0.08,第一组,III类场地,卓越周期Tg=0.45s。
图1 计算对象几点说明:●本例建模过程使用CAE;●添加反应谱必须在inp中加关键词实现,CAE不支持反应谱;●*Spectrum不可以在keyword editor中添加,keyword editor不支持此关键词读入;●ABAQUS的反应谱法计算过程以及后处理要比ANSYS方便的多。
操作过程为:(1)打开ABAQUS/CAE,点击create model database。
(2)进入Part模块,点击create part,命名为column,3D、deformation、wire。
OK (3)Create lines:connected,分别输入0,0;0,3;0,6;0,9;0,12。
OK。
退出sketch。
(4)进入property模块,create material,name:steel,general-->>density,mass density:7800,mechanical-->>elasticity-->>elastic,young‘s modulus:2.1e11,poisson’s ratio:0.3.OK(5)Create section,name:I,category:beam,type:beam, Continue, create profile, name:I, shape:I, 按图1尺寸输入界面尺寸,ok。
振型反应谱分析法和底部剪力法例题
m1 270t K1 245 MN/m
T1 0.467s T2 0.208s T3 0.134s
(2)计算各振型的地震影响系数
查表得max 0.16
Tg 0.4s
第一振型Tg T1 5Tg
1
(Tg T
) 2 max
0.139
2 m ax 0.45 m ax
(Tg T
)
2 m ax
振型分解反应谱法计算图示框架多遇地震时的层间剪
(震1设解)防:求烈体度系为的8自度振,周Ⅱ期类和场振地型,设计地震m3 分180t组K为3 第98MN二/m 组
0.334
0.667
4.019
X 1 0.667 X 2 0.666 X 3 3.035 m2 270t K2 195MN/m
第二振型V21 120.9 120.7 120.8 120.8kN
80.9kN
V22 120.7 120.8 0.1kN
107.2kN
V23 120.8
第三振型
振型分解反应谱法计算图示框架多遇地震时的层间剪
(震1设解)防:求烈体度系为的8自度振,周Ⅱ期类和场振地型,设计地震m3 分180t组K为3 第98MN二/m 组
0.063
振型分解反应谱法计算图示框架多遇地震时的层间剪
(震1设解)防:求烈体度系为的8自度振,周Ⅱ期类和场振地型,设计地震m3 分180t组K为3 第98MN二/m 组
0.334
0.667
4.019
X 1 0.667 X 2 0.666 X 3 3.035 m2 270t K2 195MN/m
120.8kN
120.7kN
120.9kN
第二振型
17.8kN 80.9kN 107.2kN
ABAQUS时程分析法计算地震反应的简单实例
时程分析法计算地震反应的简单实例时程分析法计算地震反应的简单实例(在原反应谱模型上修改)问题描述:悬臂柱高12m,工字型截面(图1),密度78003,2.1e11,泊松比0.3,所有振型的阻尼比为2%,在3m高处有一集中质量160,在6m、9m、12m处分别有120的集中质量。
反应谱按7度多遇地震,取地震影响系数为0.08,第一组,类场地,卓越周期0.45s。
图1 计算对象第一部分:反应谱法几点说明:本例建模过程使用;添加反应谱必须在中加关键词实现,不支持反应谱;*不可以在中添加,不支持此关键词读入。
的反应谱法计算过程以及后处理要比方便的多。
操作过程为:(1)打开,点击。
(2)进入模块,点击,命名为,3D、、。
(3),在分别输入0,0回车;0,3回车;0,6回车;0,9回车;0,12回车。
(4)进入模块,,:,>>,:7800>>>>,‘s :2.1e11,’s :0.3.(5),:1,:,:,, : 1, ,按图1尺寸输入界面尺寸,。
在选择I,选择。
(6),选择全部,,弹出的对话框选择:1,。
(7),选择全部,默认值确定。
(8) >> ,在弹出的对话框里勾选,,以可视化梁截面形状。
(9)添加集中质量,>>>>,:1,:,,选择(0,3)位置点,:160,。
,:2,:,,选择0,6;0,9;0,12位置点(按多选),,:120,,。
(10) >> ,选(),。
(11) >> ,:1,选,在选项卡中,选择频率提取方法,本例选用法,,选,输入10。
再,,:2,选,在选项卡中,选择单向,选择()法,:(反应谱的,后面再中添加),方向余弦(0,0,1),:1.进入选项卡,阻尼使用直接模态(),勾选,:1,:8,:0.02。
(12)进入模块,>> ,:,选择,选择,选择,选择0,0点,,勾选u13所有6个自由度。
钢管混凝土拱桥空间地震反应谱分析
两 种材 料组合 而成 , 计 算 时一 般 可有 两 种 方 法来 在 处 理 : 方法是 将 钢 管 和混 凝 土作 为 两 根 杆件 模 一种 拟 计算 , 证其 节点 坐标相 同 , 保 也就 是在两 个节 点之 问利 用两 个单元 , 分别模 拟钢 管和混 凝土 , 同受 力 共
计算 ; 另一种方法是将钢管和混凝土换算 成一种材 料来计算。本研究 中采用第二种方法 。
如果应用ansys的主从节点方法则会使刚体区域的6个自由度完全相同实际上这是不准确的本模型采取了刚性材料法即把这个区域的单元单独进行剖分通过增大材料的弹性模量来增加构件的刚度材料的弹性模量按照真实值的15倍输入这时构件交汇区的刚度相对普通构件就是一个刚性的区域
第 1 期 1
北 方 交 通
中图分类号 :4 82 U 4.2 5 文献标识码 : B 文章编 号:63— 0 2 20 ) 1— 0 1 o 17 6 5 (0 8 1 05 一 3
桥 梁振动 特性 是研 究 动 力 反应 一抗 震 、 风设 抗 计和车 桥振动 分析 的基 础 ¨ , 是进 行 其 他 动力 学 J也 分析 , 如反应 谱分析 、 响应分 析 的起 点 。虽 然我 诸 谐 国的公路规 范 当中并 没有对 桥梁 的动力特 性做 出特 殊 要求 , 是随着 桥梁 跨度 的增 大 , 但 体系 的几何非 线 性 突出 , 有必 要 研究 动力 特 性 , 因此 为桥 梁 在设 计 、 施工 、 营中的管理 提供 依据 。随着计 算机 的发展 , 运 应用 数值方法 , 即有 限元 法 把 结构 离 散 为杆 、 、 梁 壳 等各 种单元 , 不但容 易 而 且 准确 的求 出结 构 振 型 和 频率 , 还可 以考 虑结构 的细节变 化 。 运 用有 限元法 , 结 构 简化 成 梁 、 将 杆单 元 体 系 , 应 用大 型通用 程序 A S S对简 化 模 型进 行 了模 态 NY 分析 。大桥建模 技术 可 以为 同类桥 梁动力 分析 提供 有益 的思路 , 寻找应 用 A S S分 析桥梁 结构 的有 效 NY
单跨飞机库结构在多向地震下的反应谱分析
o tocos rm w r os ut n on ,h nls eu s h w dta ai o pes nrt os co i s ntes c r fh — fw r — a e okcnt ci it te ayi r l o e t a cm r i i o cnt t n o t o r t eo eO sf r oj s a s s ts h l x so ao f u r i jn h t u u t
Z
4o 0 30 o 20 0 10 0 0 2 4 6 8 1 1 1 1 1 2 2 2 0 2 4 6 8 0 2 4
机库所 在上海地 区设 防烈度 7度 , 地类别为 Ⅳ类 , 场 设计 地
震分组第 一组 , 根据抗震规范及上海抗震规程 的规定 ’ , 3 多遇地 ] 震下地震影响系数最大值 = . 8设计 特征周 期 :0 9 。 00 , .08
作用下 , 施工缝对结构整体位移的影响是随轴压 比增大而增大的。
表 4 有施工缝框架与无施工缝框架在不同轴压 比下 X方 向位移表
轴压 比
低
由
应严格分层布料 、 细致振 捣并综合 运用 二次 振捣 、 模外 振等方 法
将施工缝附近新浇筑的混凝土振捣密实 , 使其 紧密结合 。
参考文献 : [ ] G 0 1.0 1 建筑抗震设计规 范[ ] 1 B50 120 , s.
[ ] 山西建筑 ,0 1 3 ( 6 :14 . J. 2 1 ,7 1 )4 —2
2 结语
通过 A S S的模 拟证 明 , 于钢筋 混凝 土框架柱 上 的施 工 NY 对 缝, 对混凝土框架结构抗震性 能的影响是不可 忽视 的。施工缝 对
ห้องสมุดไป่ตู้
长周期结构的地震反应和位移反应谱(方小丹7.21)
• 基本周期4~5s以上的长周期结构通常为超 高层建筑,承受的风荷载较大;更高的建 筑,如基本周期6~7s以上,结构体系和抗 侧力构件断面尺寸常为避免明显的重力二 阶效应的刚度需求(也即刚重比)所控制, 结构构件有较高的承载力安全度储备。
• 由于周期长,刚度小,地震作用相对较小, 结构的“保有耐力”大,往往无需增加额外投 资即可做到中震甚至大震弹性或不屈服。因 此,对于长周期结构,通常没有结构的安全 问题。但结构的位移可能比较大。
算例三:珠海仁恒滨海中心(H=190.400m)
结构顶点位移计算值对比
αmax=0.23
Tg=0.35
Tg=0.6
Tg=0.9
Tg=1.2
Tg=1.5
XYXYXYXYXY
加速度谱计算 值Δ
0.063 0.071 0.174 0.197 0.277 0.424 0.360 0.552 0.441 0.676
• 算例一:广州西塔(H=432m)
振型 1
2
3
4
5
6
7
8
9 10 11 12
周期
7.707 7.673 2.660 2.190 2.163 1.165 1.147 0.678 0.577 0.557 0.501 0.262
结构顶点位移计算值对比
Tg=0.35
αmax=0.23
Tg=0.6
Tg=0.9
结构顶点位移计算值对比
Tg=0.35
αmax=0.23
Tg=0.6
Tg=0.9
Tg=1.2
Tg=1.5
X
加速度谱 计算值Δ
位移谱估 计值Δ
0.056 0.053
误差 -4.7%
Y 0.077 0.073 -5.5%
ABAQUS反应谱法计算地震反应实例
ABAQUS反应谱法计算地震反应实例问题描述:悬臂柱高12m,工字型截面(图1),密度7800kg/m3,EX=2.1e11Pa,泊松比0.3,所有振型的阻尼比为2%,在3m高处有一集中质量160kg,在6m、9m、12m处分别有120kg 的集中质量。
反应谱按7度多遇地震,取地震影响系数为0.08,第一组,III类场地,卓越周期Tg=0.45s。
图1 计算对象第一部分:反应谱法几点说明:λ 本例建模过程使用CAE;λ 添加反应谱必须在inp中加关键词实现,CAE不支持反应谱;λ *Spectrum不可以在keyword editor中添加,keyword editor不支持此关键词读入。
λ ABAQUS的反应谱法计算过程以及后处理要比ANSYS方便的多。
操作过程为:(1)打开ABAQUS/CAE,点击create model database。
(2)进入Part模块,点击create part,命名为column,3D、deformation、wire。
continue(3) Create lines,在分别输入0,0回车;0,3回车;0,6回车;0,9回车;0,12回车。
(4)进入property模块,create material,name:steel,general-->>density,mass density:7800mechanical-->>elasticity-->>elastic,young‘s modulus:2.1e11,poisson’s ratio:0.3.(5) Create section,name:Section-1,category:beam,type:beam,Continuecreate profile, name: Profile-1, shape:I,按图1尺寸输入界面尺寸,ok。
在profile name选择I,material name 选择steel。
地震动反应谱方法
70年代
国内
我国对反应谱研究始于60年代,王前信等人做出了大量工作, 并制定了适合于我国地质条件的设计反应谱
华南理工大学
2.2.1 反应谱的定义
有阻尼单自由度弹性体系在地震作用下的运动方程:
t 2 x t 2 x t g t x x
(1)
(5)
t , x t , x t 分别为结构相对于基底的加速度、速度和位移列向量 x
I
为单位位移时候结构各节点产生的位移
用振型分解法对(5)式解耦,将 x t 表示为阵型叠加的形式:
x t j Y i t
j 1 n
(5)
由(3)、(4)和(5)式:
Sd x t Sv x t
相对位移反应谱
m ax
相对速度反应谱
m ax
S a t g t x x
相对加速度反应谱
m ax
简单来说,反应谱是指单质点体系地震最大反应与结构自振周期 之间的关系。它是跟阻尼比和周期有关的函数。
Tg为场地特征周期
0 . 1 max0
地震影响系数曲线
水平地震影响系数最大值
烈度
max
6 0.12
7 0.23
8 0.45
9 0.90
CHONGQING UNIVERSITY
华南理工大学
2.3 振型分解反应谱法
自由度为N的弹性体系在地面运动加速度为
ut 其运动方程为
M t C x t K x t M I ut x
j max
E max
j
t
D j ,
ANSYS地震分析实例
ANSYS地震分析实例土木工程中除了常见的静力分析以外,动力分析,特别是结构在地震荷载作用下的受力分析,也是土木工程中经常碰到的题目。
结构的地震分析根据现行抗震规范要求,一般分为以下两类:基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。
本算例将以一个4质点的弹簧-质点体系来说明如何使用有限元软件进行地震分析。
更复杂结构的分析其基本过程也与之类似。
关键知识点:(a) 模态分析(b) 谱分析(c) 地震反应谱输进(d) 地震时程输进(e) 时程动力分析(1) 在ANSYS窗口顶部静态菜单,进进Parameters菜单,选择Scalar Parameters选项,在输进窗口中填进DAMPRATIO=0.02,即所有振型的阻尼比为2%(2) ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete,添加Beam 188单元(3) 在Element Types窗口中,选择Beam 188单元,选择Options,进进Beam 188的选项窗口,将第7个和第8个选项,Stress/Strain (Sect Points) K7, Stress/Strain (Sect Nods) K8,从None 改为Max and Min Only。
即要求Beam 188单元输出积分点和节点上的最大、最小应力和应变(4) 在Element Types 窗口中,继续添加Mass 21集中质量单元(5) 下面输进材料参数,进进ANSYS主菜单Preprocessor->Material Props-> Material Models菜单,在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic 属性,输进材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。
结构地震反应分析详解演示文稿
D 1 2
由初始条件: x0 x0 x0 x0
xt et[x0 cosDt
x0
c1 x0
c2
x0
x0 D
x0 D
sin
Dt]
第八页,共94页。
xt
et [x0
cos Dt
x0
x0 D
sin Dt]
有阻尼单自由度体系自由振动位移时程
D 1 2
有阻尼自振频率
当c=0,无阻尼,即ξ=0
强迫振动
当地面运动为冲击运动时的特解
xg ( )
冲击力: P mxg 0
0 dt
dt
xg
质点0~dt时间内的加速度为
a
p m
xg
dt时刻的速度和位移分别为:
v xg dt
d
1 2
xg dt2
0
dt
冲击运动使质点产生初速度,然后进行自由振动
第十页,共94页。
将 得
v xg dt 代入 xt et[x0
第二十页,共94页。
3.设计反应谱
由反应谱计算地震作用: ?
F
m
xg
x
max
m
xg
x max
mSa T
需要设计反应谱
F mSa T mg
xg max g
Sa T
xg max
Gk T
G k T
G ---作用在质点处的重力荷载代表值
k ---地震系数
T ---动力系数
第二十一页,共94页。
度 ic EIc / h 2.6104 kN m ,阻尼比为0.05。试求该结构多 遇地震时的水平地震作用。
解:(1)求结构体系的自振周期
地震反应谱分析实例
结构地震反应谱分析实例在多位朋友的大力帮助下,经过半个多月的努力,鄙人终于对结构地震反应谱分析有了一定的了解,现将其求解步骤整理出来,以便各位参阅,同时,尚有一些问题,欢迎各位讨论!为叙述方便,举一简单实例:在侧水压与顶部集中力作用下的柱子的地震反应谱分析,谱值为加速度反应谱,考虑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,45MP,EX,1,2.8E10MP,DENS,1,2.4E3MP,NUXY,1,0.18!建立模型BLOCK,0,1,0,1,0,5!网格剖分ESIZE,0.5VMESH,all/VIEW,,-0.3,-1,1EPLOTFINISH/SOLU!施加底部约束ASEL,,LOC,Z,0DA,ALL,ALLALLSEL!施加自重荷载ACEL,0,0,10!进行模态求解ANTYPE,MODALMODOPT,LANB,30SOLVEFINISH!进行谱分析/SOLUANTYPE,SPECTRSPOPT,SPRS,30,YESSVTYP,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.0000SV,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.2167SOLVEFINISH!进行模态求解(模态扩展)/SOLUANTYPE,MODALEXPASS,ONMXPAND,30,,,YES,0.005SOLVEFINISH!进行谱分析(合并模态)/SOLUANTYPE,SPECTRSRSS,0.15,dispSOLVEFINISH/POST1SET,LIST !结果1/INP,,mcomLCASE,11PRRSOL, !结果2SET,FIRSTPRRSOL, !结果3SET,NEXTPRRSOL, !结果4SET,NEXTPRRSOL, !结果5SET,NEXTPRRSOL, !结果6 FINISH!静力分析/SOLUANTYPE,STATIC!施加水压荷载NSEL,,LOC,Y,0NSEL,R,LOC,Z,0,5 SFGRAD,PRES,0,Z,0,-10000 SF,ALL,PRES,50000!施加集中荷载NSEL,,LOC,Y,0NSEL,R,LOC,z,5F,ALL,FY,10000ALLSELEPLOTSOLVEFINISH/POST1set,lastlcwrite,12Lcase,11Lcoper,add,12Lcwrite,13LCASE,12PRRSOL, !结果7PRRSOL, !结果8FINISH---------------------------------------------------------------------------------------------------以下是计算的结果---------------------------------------------------------------------------------------------------结果1:(Results Summary)1 21.6472 21.6473 121.514 121.51结果2:(单独谱分析反力LCASE,11)VALUE 2467.9 2290.1 18384.结果3:(单独谱分析反力SET,FIRST)VALUE 0.13334E+06-0.15785E+07-0.18819E-06结果4:(单独谱分析反力SET,NEXT)VALUE -0.15785E+07-0.13334E+06-0.48918E-06结果5:(单独谱分析反力SET,NEXT)VALUE -0.87805E+07 0.27008E+08 0.86846E-07结果6:(单独谱分析反力SET,NEXT)VALUE 0.27008E+08 0.87805E+07 0.79325E-06结果7:(单独静力分析反力LCASE,12)VALUE 0.22901E-08-0.15500E+06 0.12000E+06结果8:(谱分析与静力分析叠加反力LCASE,13)VALUE 2467.9 -0.15271E+06 0.13838E+06---------------------------------------------------------------------------------------------------以下是问题的讨论---------------------------------------------------------------------------------------------------1、模态提取数为30,即取前30阶振型数,但在谱分析时得到的是4阶,这4阶是什么意思?2、在单独谱分析时,为何结果2、3、4、5、6会相差如此之大?(其应力和位移也是如此。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3、在进行谱分析的合并模态步骤中,模态合并方法应该选用SRSS,0.15,DISP还是SRSS,0.15,ACEL(其对结果影响很大)?
4、在谱分析中,似乎应该将水压力折算成各节点处的附加质量MASS21,但对于非竖值方向(如本例为Y方向),MASS21单元如何加上去?
3 121.51
4 121.51
结果2:(单独谱分析反力LCASE,11)
VALUE 2467.9 2290.1 18384.
结果3:(单独谱分析反力SET,FIRST)
VALUE 0.13334E+06-0.15785E+07-0.18819E-06
结果4:(单独谱分析反力SET,NEXT)
在多位朋友的大力帮助下,经过半个多月的努力,鄙人终于对结构地震反应谱分析有了一定的了解,现将其求解步骤整理出来,以便各位参阅,同时,尚有一些问题,欢迎各位讨论!
为叙述方便,举一简单实例:
在侧水压与顶部集中力作用下的柱子的地震反应谱分析,谱值为加速度反应谱,考虑X与Y向地震效应作用。已知地震影响系数a与周期T的关系:
LCASE,11
PRRSOL, !结果2
SET,FIRST
PRRSOL, !结果3
SET,NEXT
PRRSOL, !结果4
SET,NEXT
PRRSOL, !结果5
SET,NEXT
PRRSOL, !结果6
FINISH
!静力分析
/SOLU
ANTYPE,STATIC
以下是问题的讨论
---------------------------------------------------------------------------------------------------
1、模态提取数为30,即取前30阶振型数,但在谱分析时得到的是4阶,这4阶是什么意思?
/VIEW,,-0.3,-1,1
EPLOT
FINISH
/SOLU
!施加底部约束
ASEL,,LOC,Z,0 D Nhomakorabea,ALL,ALL
ALLSEL
!施加自重荷载
ACEL,0,0,10
!进行模态求解
ANTYPE,MODAL
MODOPT,LANB,30
8、有朋友说可以“做模态分析之前先做打开预应力的静力分析,然后再做打开预应力的模态分析,接着和以前一样”,但我试了试,得到的结果相差很大,不知何故,如果哪位朋友曾经作过这方面的例子,能否再试一下本例,以便进行比较!
!施加水压荷载
NSEL,,LOC,Y,0
NSEL,R,LOC,Z,0,5
SFGRAD,PRES,0,Z,0,-10000
SF,ALL,PRES,50000
!施加集中荷载
NSEL,,LOC,Y,0
NSEL,R,LOC,z,5
F,ALL,FY,10000
ALLSEL
结果7:(单独静力分析反力LCASE,12)
VALUE 0.22901E-08-0.15500E+06 0.12000E+06
结果8:(谱分析与静力分析叠加反力LCASE,13)
VALUE 2467.9 -0.15271E+06 0.13838E+06
---------------------------------------------------------------------------------------------------
5、在选取(离散)频率和谱值时,似乎没有一定的标准,而其基本原则是须选取具有代表性的点。
6、模态分析中采用Subspace法与Block Lanczos法所得结果相差2%~4%左右,用后一种方法计算更快,但占用内存大50%左右。
7、在网格剖分时采用VMESH和VSWEEP两种方法,对静力分析结果没有影响,但对谱分析结果却有一定影响,何也?
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
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
EPLOT
SOLVE
FINISH
/POST1
set,last
lcwrite,12
Lcase,11
Lcoper,add,12
Lcwrite,13
LCASE,12
PRRSOL, !结果7
LCASE,13
PRRSOL, !结果8
FINISH
---------------------------------------------------------------------------------------------------
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
/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
MXPAND,30,,,YES,0.005
SOLVE
FINISH
!进行谱分析(合并模态)
/SOLU
ANTYPE,SPECTR
SRSS,0.15,disp
SOLVE
FINISH
/POST1
SET,LIST !结果1
/INP,,mcom
lcwrite,11
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 秒
以下是命令流程序
----------------------------------------------------------------------------------------------------
以下是计算的结果
---------------------------------------------------------------------------------------------------
结果1:(Results Summary)
1 21.647
2 21.647
VALUE -0.15785E+07-0.13334E+06-0.48918E-06
结果5:(单独谱分析反力SET,NEXT)
VALUE -0.87805E+07 0.27008E+08 0.86846E-07
结果6:(单独谱分析反力SET,NEXT)
VALUE 0.27008E+08 0.87805E+07 0.79325E-06