工程抗震概述.doc

硕士研究生非笔试课程考核报告（以论文或调研报告等形式考核用） 2023 至 2023 学年 第 2 学期考核课程： 工程抗震 提交日期： 2023 年 7 月 15 日报告题目： 工程抗震综述 考核成绩考核人姓 名 学 号 年 级 2023级 专 业 建筑与土木工程 所在学院 土木工程学院 山东建筑大学研究生处制工程抗震综述研究生第一年已经接近尾声，这学期学习了工程抗震这门课，让我对抗震这门课有了更进一步的理解，以及对结构设计来说，抗震设防的重要性和必要性。与此同时，也学习了常见结构的抗震措施和解决方式，从理论上有了宏观的把握，对以后结构设计打下了基础。工程抗震这门课内容很多，从“地震动特性”到“结构消能、减震”等等。老师结合建筑抗震设计规范,强调了结构概念设计，推导了规范最基本的相关公式，如底部剪力法、振兴分解反映谱法等，介绍了目前国内外关于抗震研究现状，加深了对强柱弱梁、强剪弱弯和强节点若杆件等结构设计概念的理解。现在结合工程抗震老师上课所讲和自己的笔记，简朴谈谈对这门课程的结识和理解。一，地震动及其特性 随着经济建设的迅速发展和社会财富的快速积累，灾害的防治越来越受到重视。灾害是指那些由于自然的、人为的或人与自然的因素，对人类的生存和社会发展导致损害的各种现象。在国家自然科学基金委员会工程与材料学部的重大工程灾害与防治论证报告中给出了灾害的一种更科学的定义：灾害是指自然环境或工程系统演变过程中失去固有平衡或稳定期导致人类赖以生存的基础破坏或功能失效的突发事件。联合国成立了国际减灾十年委员会，希望通过国际间的合作与努力，力求在十年中有效地减少和防治自然灾害。国际减灾十年委员会（19902023）提供的资料表白：1963年至1992年的三十年内，全球自然灾害的受灾总人数30.08亿人，死亡总人数36万人，直接经济损失3400亿美元，其中地震引起的灾害占很大的比例。1.地震的成因（Cause of EQS） 根据地震成因来划分地震的类型： 塌陷地震（Collapse EQS）：50年来仅1954和1965年四川自贡发生两次。 火山地震（Volcanic EQS）：火山喷发引起的大地振动。 构造地震（Tectonic EQS）：地壳构造变动引起的地震，占90以上。破坏性地震一般均为构造地震。以后讨论的地震均为构造性地震。 引发构造地震的因素： 外因地震成因的宏观背景：板块构造学，用板块构造学上升到理论模型 内因地震成因的局部机制：弹性回跳理论理论模型 板块构造运动（二十世纪60年代提出）： 全球地壳可以分为六大板块（太平洋、欧亚、美洲、非洲、印度澳洲、南极），各大板块内部可以提成较小的板块。板块厚80到100千米（地壳30到40千米加上岩石层40到70千米）。板块在上地幔推动下发生移动，约25cm/年，各板块运动不均匀，同一板块的运动也不均匀。板块运动引起岩层变形积累，达成一定限度忽然破裂释放的巨大能量以波的形式向外传播，引起地面运动，迫使结构振动，产生破坏。在大板块边沿地震最多，这些地震称为板缘地震（浅源）。从世界地震的震中分布图上可以看到两条大地震带： 1）环太平洋地震带占世界地震总数75以上。2）欧亚地震带（至印度与环太平洋地震带相遇）占地震总算22。这两个地震带属于板缘地震，大小板块边沿都是地震集中处。 有些地震不发生在边沿，而是在板块内，成为板内地震（或称内陆地震），约占15分布零散，危害性大（震源浅，震级大，地震区人口集中）；机制复杂，如唐山地震。中国地震：中国位于欧亚板块东南端，东为太平洋板块，南为印度洋板块，欧亚板向东、太平洋板块向西、印度洋向北挤压中国大陆，因此中国是地震多发国家。除少数省外，都发生过6级以上破坏地震，基本烈度大于6度，占国土总面积的60。全国450个城市，74.5位于地震区，约一半城市位于基本烈度7度及其以上地区。28个百万以上人口的大城市有85.7位于地震区。本世纪全世界发生的7级以上强震中，中国占三分之一。2.地震的描述震源（Hypocenter，Source，Focus，Center）：地壳深处岩石发生断裂、错动的地方，即一次地震时应变能释放区。震中（Epicenter，Epifocus）：震源在地表的投影点。断层：地震形成的地壳岩石中的大断裂面（可达成几十几百公里长）。地震时破裂开始于断层面上的一个点或一个局部，然后向四周扩展。 事实上震源不是一个点，而是具有相称尺度（大小）的破裂面，与地震的大小有关。地震学在研究中采用三种震源模型（理论模型）： 点源：中小震（M6），断层长度：几十到几百公里； 面源：中大震（M6），断层长度：几十到几百公里。对一次地震震源（的位置）（Hypocenter）认为是断层上一方面破裂的那个点（的位置）。能量中心：能量释放中心，一般是断层破裂面的几何中心。小震：震源能量中心。 大震：震源和能量中心不重合。3.地震的三要素 发震时间时间 发震地点空间：震中位置（经度，纬度）震源深度 地震大小强度 对已发生的地震，实际的震中、震源位置的拟定会有几公里到几十公里的误差，极个别的也许存在百公里的误差。一个地震由下列参数定义：时间、空间、强度、简称时、空、强三要素。对一个已发生地震的拟定需要给出三要素，在地震预报中预报一个地震也需要拟定这三个要素。根据震源深度的划分，地震可以分为： 浅源地震（shallow focus earthquake）：震源深度，h70km，占总地震的72 Elcentro 1940 EQ M=7.1，h=24Km 唐山 1976 EQ M=7.8，h=16Km 特点：峰值加速度大，卓越频率520HZ，高频成份多。中源地震（Intermediate focus EQ）：70Kmh300Km，占总地震的4 特点：峰值加速度小，卓越频率0.23HZ（高频被滤掉）。4， 地震动（1）地震波分为弹性波面波（次声波）；其中以面波对人类导致的危害较大。（2）地震动的反映取决于三要素：峰值、频谱和连续时间。对于三要素的影响因素: 震源机制、震中距、传播途径及场地条件。5.、地震动的估计及模拟（1） 烈度估计：得到烈度和参数之间的关系。（2） 衰减规律。（3） 震源机制（4） 人造地震动分为地震学方法和工程方法。具体方法为比例法和级数法。6.、地震的宏观评价“震害指数评价法”：破坏作用分级，总分为：、。二、 减轻震害的措施1、 能较及时的预报地震发生的时间、地点和强度。但是地震具有群发性、迁移性、等距性、复发性、前兆性和大小年，增长了预测的难度。2、 转移分散。3、 地震工程方法：地震危险性分析+结构抗震设计；其中结构抗震设计是比较常见的传统的方法。三、 结构抗震中的不拟定因素1、 地面运动的不拟定性。非平稳随机过程，地震的记录资料少。短时间的作用。有选择的破坏作用。2、 结构自身动力特性的不拟定性：涉及自重、刚度、自振周期。结构自重：基础与上层建筑结构质量分布不均匀。结点刚度：导致力的作用的偏心、扭转和P-效应。非结构构件：弹塑性性质。3、 结构构件的材料对地震作用的影响不同：涉及时间、速率、弹塑性能和反复加载的影响。4、 阻尼的影响，弹性状态下的阻尼系数2-3%。5、 结构分析模型的不拟定性根据结构的特点以及研究的构件等因素可以给结构建立不同的模型，例如：平面模型、平面协同模型、空间模型、杆系模型、墙元模型、壳元模型等。6、 抗震规范的局限性。7、 抗震的基本需求,其指标有：刚度、承载力和稳定性。刚度是指有足够的吸能能力，承载力指的是足够的耗能能力，稳定性指的是足够的变形能力。四、 结构抗震设计理论的演变1、 静力法阶段（20世纪2030年代） 惯性力F=ma=a=KG其中：a为地面运动的加速度；k为震设系数，可取0.1. K=，Sn地面运动最大加速度。2、 反映谱阶段（1）40年代初：计算结构的反映谱（2）1952年将该方法列入规范，即反映谱理论。反映谱理论结构多自由度体系弹性阶段分析。用振型分解（组合）法总反映。反映谱法单自由度最大反映。（3）1970年，新西兰的Parly pank提出能力设计法承载力变形能力：耗能能力(i)、合理的破坏机制。(ii)、防止脆性破坏。(iii)、变形能力设计3、 基于性能的抗震设计（1） 概念：工程结构在预定的使用年限内，在不同强度水平地震作用下达成预定的不同的性能目的。 生命安全 性能目的 结构的损伤控制 正常使用（2） 基于位移的性能化设计（以位移为控制指标）a:能力谱法 :将多自由度体系等效为单自由度体系得到“”曲线，然后转化为“”加速度与位移曲线。 多自由度体系的 单自由度体系的 加速度与位移曲线“MD” 加速度与位移曲线“SD”其中第一阵型的模态质量为，阵型参与系数为。则存在下列公式：= =b:求需求谱i：由建筑抗震设计规范中的方法可得“”的曲线，其中=ii：特定的地震动得到“”的曲线。 性能点 T iii:找性能点，与的交点即为性能点，若无交点，则说明抗震能力局限性。其中性能点所相应的最大位移为，而实际结构性能点位移=。（3） 基于预案拟定设计方案由政府主管部门拟定性能目的、指标由业主来拟定 性能目的的选用是否用现行规范方式设计？建筑结构的方案实行目的的措施新技术 新方法具体的设计、必要的实验专家认证通过目的评价经济分析 不满足性 能目的 否 五、 地震危害性评价1、概述：地震危害性的评价是对地震的活动性、分布规律的研究，涉及地震的参数、强度的研究。目的：提供设防的依据 对重大工程项目进行场地分析2.、评价的内容（1） 地震活动性分析（150km以内的地震活动） （2） 地震构造的调查（2）近场区域地震构造活动的调查分析（3） 烈度及地震动衰减关系（4） 危害性分析（5） 场地工程地震条件（6） 场地的地震灾害评价（7） 场地的地震反映分析（8） 提出场地的地震动参数六、 现行抗震规范所体现的设计原则1.设防的标准建筑工程抗震设防分类标准中对设防烈度分为甲、乙、丙、丁四个等级，可以在以下情况设防(1)地震作用的计算（2）采用抗震措施时设防（一） 概率为基准的设计方法1、 三水准设防。涉及第一水准设防，（小震）为众值烈度，重现期为63.2%，为弹性状态；第二水准设防，（中震）为基本烈度，重现期为10%，部分进入弹塑性状态；第三水准设防，（大震）罕遇地震，重现期23%.地震发生的概率符合概率公式： 密度函数： f= 分布函数： F= f 63.2% 10% 23% 基本烈度 罕遇地震 F 1/600 1/200 1/50 地震作用引起的变形区域代表结构所能抵抗地震的能力区域代表能修阶段所能抵抗地震的能力区域代表继续消耗能量不超过1/50的层间变形，即为可充足变形来抵抗地震作用的能力。2、 二阶段设计方法第一阶段设计：第一水准强度规定要满足； 基于概念设计的内力调整和验算 位移线性（小震下的弹性位移） 延性设计第二阶段设计 ：对于易倒塌的结构。（三）基本框架1、 结构抗震分析构件设计。2、 设防烈度+结构特性=构件的抗震等级。3、 三水准设防+二阶段设防七、 关于概念设计概念设计是指一些在计算中或在规范中难以作出具体规定的问题，必须由工程师运用“概念”进行分析，作出判断，以便采用相应的措施。例如结构破坏机理的概念，力学概念以及由震害实验现象等总结提供的各种宏观和具体的经验等，这些概念及经验要贯穿在方案拟定及结构布置过程中，也体现在计算简图或计算结果的解决中。建筑结构的抗震设计，是以现有科学水平和经济条件为前提的。目前地震及结构所受地震作用尚有许多规律未被结识，人们在总结历次大地震灾害经验中结识到：一个合理的抗震设计，在很大限度上取决于良好的“概念设计”。抗震概念设计重要有如下几点：1建筑的体型力求简朴、规则、对称、质量和刚度变化均匀。2抗震结构体系，应符合以下规定：（1）具有明确的计算简图和合理的地震作用传递途径；（2）具有多道抗震防线，避免因部分结构或构件破坏而导致整个体系丧失抗震能力或对重力荷载的承载能力；（3）应具有必要的强度、良好的变形能力和耗能能力；（4）具有合理的刚度和强度分布，避免因局部削弱或突变形成薄弱部位，产生过大的应力集中或塑性变形集中；对也许出现的薄弱部位，应采用措施提高抗震能力。3抗震结构的各类构件应具有必要的强度和变形能力（或延性）。4抗震结构的各类构件之间应具有可靠的连接。5抗震结构的支撑系统应能保证地震时结构稳定。6非结构构件（围护墙、隔墙、填充墙）要合理设立。建筑抗震概念设计是根据地震灾害和工程经验等所形成的基本设计原则和设计思想，进行建筑和结构总体布置并拟定细部构造的过程，是一种基于震害经验建立的抗震基本设计原则和思想。八、地震作用的分析（一）分析方法等效荷载法涉及静力法和反映谱法动力方程直接积分：可以得到地震反映时间关系（及时程关系）。（二） 单自由度体系的反映（右图为单自由度体系的模型）惯性力：m+ m 恢复力：k 阻尼力：c k由达朗贝尔原理得 m+c+k=-m关于阻尼：粘滞阻尼与c成正比。 摩擦阻尼与常数成正比。 滞变阻尼：与成正比。 流体阻尼：与成正比。齐次通解为= 当时，不衰减； 当时，衰减； 当 即时，不振，此时，无阻尼。 单自由度体系在地震作用下的动能反映式为 即杜哈默积分，包含位移反映、速度反映、加速度反映。（三）地震作用 可求出“时间历程”的最大值。绝对加速度 的影响因素有、频谱、强度（烈度）、震中距、场地特性、刚度、质量、阻尼、固有特性。对于单自由度体系：所以=的物理意义：最大惯性力与重力的比值，考虑了动力特性及结构自身特性。下图为谱： 0.45 T（四）多质点（多自由度系统）弹性体系地震反映。 1、一般运动微分方程度自由度体系的运动方程为2 、地震反映 3 、固有振动 自由振动的一般解为 4、主阵型矩阵 其中， 是主质量； 是主刚度主阵型矩阵为 5 、多自由度体系的地震作用下的运动方程。 整理得： 得最终多自由度体系的反映时程： 6 、阵型的叠加原理 低阶阵型的奉献较大（1） 对大量低矮的房屋，自由度少，一般取13阵型。（2） 高层建筑，自由度多，一般取915阵型。（3） 桥梁、高耸建筑去更多阵型。阵型的参与质量要保证大于等于90%。线性系统为所以，地振反映时程为“反映时程”为， 此为随机振动理论。7、阵型分解反映谱理论 结构选用多自由度系统。多自由系统的反映是采用振型分解法由单自由度系统的反映叠加组合而成。单自由度系统的反映有反映谱法拟定。基本假定 反映为线弹性反映，可以使用叠加原理。 基础是刚座，所有支撑点地震作用相同。 最大反映最不利。 地震动属于平稳随机过程。（1）反映谱法原理通过坐标变换，可得到最大反映值为 随机抗震理论（CQC）（考虑地震动相关性，是建立在地震效应的基础上。）地震反映最大值为： -阵型相关系数当 时， 此为平动和扭动的耦合作用。CQC-SRSS法 （2）地震作用 最大值为 在此基础上，将地震作用作用在物体上产生的地震效应，再进行叠加。（五）规范的阵型分解反映谱法1、阵型最大的地震作用 第j振型的第i个质点 第j个振子单自由度体系的振型影响系数，其中重力代表值 计算地震作用效应 ，n个振型的地震作用效应。2、阵型组合（叠加） 结构比较柔，即自振频率慢，自振周期长，高阶阵型影响。（六）底部剪力法对阵型分解反映谱法多阵型的组合。对于结构较刚，,质量和刚度分布均匀，受剪切变形，高度在40m以内比较合用。地震反映以第一阵型为主。 等效重力系数一般取0.85，其计算公式为其中，等效重荷载 由 得到 当 时，此时对女儿墙等边端部位，作用在该处的作用再扩大三倍解决。（七 ）、结构竖向地震作用1、竖向地震作用的影响受竖向地震作用影响的结构重要是距离震源200km以内， 在1/21.6的长悬臂结构；大跨结构；高耸、高层结构；柱、悬挂结构等以轴向力为主的结构构件；砌体等脆性材料结构。2 计算方法（1）静力法（2）水平地震作用折减法（3）竖向地震反映谱 （4）直接动力法（验算手段）3、反映谱法反映谱法重要合用于高耸、高层建筑结构。， 4、屋盖、屋架屋盖、屋架结构的就算一般使用重力荷载代表值乘以地震作用系数。5 、长悬臂和其他大跨度结构。（八）建筑结构的扭转地震效应1、因素（1）结构偏心、不规则引起的平面作用与扭转作用的偶联。（2）地面的扭转分量的影响，此时不偶联。2、方法（1）规则结构不考虑扭转偶联计算，对于周边结构，地震作用增大。（2）不规则结构需要考虑平-扭偶联计算。基本假设 楼盖平面内绝对刚性，即发生刚体位移。 结构系统是由多个平面抗侧力结构构成。 结构的抗扭刚度是由多个平面抗侧力结构的侧向恢复力提供。 质量集中于楼盖。第j阵型第i层的水平地震作用： 单项水平地震作用效应为： (CQC) 双向地震作用影响两者取不利值。（九）结构基本周期的近似计算1、能量法 得到 2、顶点位移法弯曲振动：剪切振动：弯-剪型振动：3、等效质量法4、关于对周期（刚度）的修正 时程反映，仿真分析。九、地震反映分析的直接动力法 其中，-考虑阻尼模型的比例系数； -数值方法的离散函数。初始状态： 对此初值问题在 进行时域分析、品与分析、时程分析。采用数值积分法作为近似计算方法。其特点是： 能全程满足，在有限时间节点上满足。节点个数为 个，其中 ，。误差分析：重要有阶段误差和舍去误差，这些误差影响收敛；积累误差回导致数值稳定问题。（一）弹性体系的分析1、 （全量分析）假设“”的关系，采用常量加速度方法 对 积分，得 进行代换，最终带入动力方程可得环节分析：（1） 准备工作 求出 、 拟定初值 、 在 范围内取增值 ，拟定、的取值。 形成等效刚度矩阵 （2） 分析： 时刻的等效荷载 时刻的等效荷载 由此求出 求得 、 ，拟定 的取值 2、线性加速度法（增量） 推导： 为常量使用泰勒级数展开，进行多步代入求解，可得到拟刚度矩阵 和拟荷载向量 可以得到 ， 的选取依据： 结构的自振周期 地震波的卓越周期 地震波的加速度变化率 非线性刚度、阻尼当时失效，常用3、法用代替各式中的，得 （卓越周期）时，无条件稳定。（二）恢复力特性 1、弹塑性恢复力曲线（1）恢复力变形曲线（2）非线性分析的基础（3）三个层面，分别是：材料层面、构件层面、结构层面（4）两个要素，分别是骨架曲线和滞回曲线，它们所包含的信息有能量吸取与消散，延性，强度、刚度及蜕化性能，阻尼性能等。2、钢筋混凝土构件（1）受弯构件 滞回环饱满（2）压弯构件（N-M）（3）节点 梁柱节点钢筋滑移导致节点滞回曲线收缩。模型分为直线型模型和曲线型模型两种，直线型模型涉及双线性和三线型两种。3、重要模型及参数 （1）双线型模型钢结构的双线型模型重要有： 抱负弹塑性双线型，参数为:， 硬化双线型，参数为:， 退化双线型，参数为:，（2） 三线型模型退化三线型的参数为：，（三）动力分析模型及刚度矩阵1、层模型（1） 层间剪切型模型 P 进行弹性分析，第j层层间剪切刚度，框架柱为， 考虑梁柱转角，则对结果进行分析可得到层间剪力和层间位移。关于薄弱层的分析，可对三线型模型进行弹塑性分析 层间开裂剪力 开裂位移对于j层 层间侧移刚度 屈服剪力，对j层进行状况分析，可发现有三种情况：强柱弱梁、强梁弱柱和混合型。对于强梁弱柱型 对于强柱弱梁型 对于混合型：应具体分析对待。若是强梁弱柱型，按实际计算，若是强柱弱梁型，采用等效计算。 屈服位移其中，-弹塑性阶段割线刚度系数（2）层间弯剪模型此时刚度矩阵为满阵。 柔度矩阵求逆法静力法方程： 对每一层施加单位水平力，则各层的水平位移为 ： ： 等效刚度矩阵法等效刚度矩阵法的求解方法是建立静力法方程。对下一层施加荷载，计算各层的剪力，对于j层 2、杆模型通常以梁、柱为单元。杆单元的特点： 质量集中于节点，即质点。 质点的自由度有水平方向、竖直方向和自由转动；刚性楼盖的自由度只有水平方向的。 刚度矩阵的建立，采用聚缩法对于弹塑性分析：单分量模型：塑性变形所有集中于杆端塑性铰区。 完全弹性双分量模型：将杆件提成两部分，如下 弹性部分 塑性部分三分量模型： 弹性部分 开裂后，屈服前 屈服后十、结构耗能、减震。1、 概述 隔震装置特性：（1）竖向承载力（2）隔震性能（3）复位特性（4）阻尼耗能2、 动力分析模型令代入上式得：抗震装置的自振频率阻尼系数阻尼比为：隔震效果分析： =0.10.3当时，可减少2.55度。3、 隔震结构的计算设计反映谱法采用隔震结构的地震影响系数：=所以= 工程抗震程序分析十一.时程分析法计算结构在地震作用下的响应1.1.使用软件为平面结构弹塑性地震响应分析软件NDAS2D。1.2.单位：质量：吨（t） 时间：秒（sec） 长度：米（m） 力：千牛（KN）1.3.分析结构为一跨两层钢筋混凝土框架结构（图10）：图10 框架计算简图层高4.2m;柱距3.9m，柱尺寸：0.4mX0.4m ;梁尺寸：0.25mX0.40m;混凝土采用：C301.4单元属性 柱 截面积：0.16 ; 惯性矩：0.002278 弹性摸量：30000 梁 截面积：0.1; 惯性矩：0.0013333弹性摸量：300001.5结点质量梁自重：0.25x0.4x3.9x2.5=0.975t柱自重：0.4x0.4x4.2x2.5=1.68t转化为结点自重：3,4为2.2t;5,6为1.3t;1.6阻尼参数采用瑞雷阻尼，即阻尼矩阵为初始刚度矩阵和质量矩阵的线性组合 ： 1.7采用的地震波采用埃而森特罗波，地震时间间隔0.02秒。 -1.4 -10.8 -10.1 -8.8 -9.5 -12.0 -14.2 -12.8 -11.0 -8.5 -8.5 -13.1 -17.6 -19.4 -16.2 -14.4 -10.8 -8.2 -4.2 -6.6 -13.1 -19.0 -19.6 -6.6 3.0 14.1 -4.9 -12.8 -14.4 -20.3 -26.0 -32.5 -30.6 -17.2 -19.7 -16.3 -16.4 -6.7 2.5 15.0 23.6 25.2 33.6 46.3 49.2 41.9 35.9 27.1 23.5 33.9 41.2 53.0 63.9 73.2 65.2 59.9 40.0 40.0 6.3 -51.5 -78.7 -60.3 -48.4 -25.0 -5.9 13.4 30.8 49.9 71.0 99.5 121.9 152.9 144.9 115.5 93.5 89.2 92.6 83.9 90.1 99.3 120.9 32.8 -147.5 -206.6 -198.9 -203.4 -181.6 -172.5 -175.2 -175.3 -180.5 -163.0 -134.7 -108.7 -78.2 -42.9 -1.7 36.0 78.5 116.4 159.8 196.0 241.2 272.9 303.6 320.0 341.7 282.1 232.4 -119.8 -237.3 -164.0 -186.5 -109.5 -75.3 -17.3 11.3 53.3 89.5 118.6 175.7 57.6 -263.1 -154.7 -172.9 -101.2 -57.9 23.7 -67.0 -198.0 -164.1 -168.5 -148.1 -123.1 -100.1 -75.1 -52.3 -27.1 -4.4 18.8 -9.5 -43.3 -83.8 -95.1 -71.6 -59.9 -33.4 -10.8 18.5 42.0 67.3 -9.7 -37.2 -4.0 1.1 34.4 56.5 88.3 113.0 136.3 21.9 24.1 68.3 68.9 131.8 135.3 204.0 -93.1 -130.8 -69.2 -54.6 7.2 67.5 -106.7 -148.8 -107.1 -116.2 -76.2 -55.9 -21.5 -12.6 -67.4 -32.4 -33.7 -10.9 1.7 29.9 48.8 60.8 22.2 -3.2 -24.5 7.7 21.1 56.8 82.6 120.6 147.8 173.7 42.1 2.9 25.9 29.3 -5.5 -14.7 14.3 20.6 49.9 64.5 95.7 112.8 144.7 162.9 194.5 185.6 198.4 176.9 125.0 -120.7 -54.2 -38.4 -31.1 -111.8 -166.1 -246.4 -202.5 -183.5 -131.7 -96.0 -32.5 15.4 81.6 131.9 181.8 -5.8 -16.9 28.5 44.7 98.3 142.4 185.3 245.6 168.5 -138.0 -99.9 -108.9 -90.7 -46.9 -125.0 -211.1 -161.7 -169.2 -130.6 -111.1 -77.3 -51.0 -54.4 -120.0 -120.9 -115.8 -114.5 -71.7 -54.6 6.4 -80.4 -163.4 -85.9 -96.1 -39.6 -14.7 31.9 64.8 87.6 47.2 19.8 -