三向地震波的合理选取和人工定义

18.1。时程分析与三向地震波18.2。三向地震波的合理选取18.3。如何人工定义地震波18。三向地震波的合理选取和人工定义18.1。时程分析与地震波弹性、弹塑性时程分析均与地震波相关。TAT、SATWE、PMSAP、EPDA等软件时程分析时均需选取地震波。旧版软件采用的是按照场地土区分的单向地震波库；新版软件采用的是按照特征周期区分的三向地震波库。三向地震波可以退化为单向地震波进行计算。可以通过填写文本文件的方式增加用户地震波。新抗震规范5.1.2条规定，“特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑，应采用时程分析法进行多遇地震下的补充计算”，“采用时程分析法时，应按建筑场地类别和设计地震分组选用不少于两组的实际强震记录和一组人工模拟的加速度时程曲线”。地震波与反应谱应在“统计意义上相符”。时程分析法单波和平均值的底部剪力应不小于按反应谱方法得到的底部剪力的“65％”和“80％”等限值。新抗震规范5.5.3条规定，除可以采用简化方法计算外的建筑结构，可采用静力弹塑性分析方法或弹塑性时程分析方法。18.2。三向地震波的合理选取按照规范的要求，至少应该选择三条地震波进行地震时程反应的分析，并规定了最小基底剪力。当计算的基底剪力不满足规范要求，则应认为该地震波不合格，应重新选择分析，直至选到合适的地震波为止。而实际上，只有在建筑物所在地的地震波才有可能有意义。但是大多数地区不具备这个条件，则可以用实测的人工波来代替。目前重要建筑物的场地波都是通过实测和人工模拟产生的，即实测人工波。选取地震波实测地震波——特征参数实测地震波——反应谱周期(秒)规准加速度谱1234563.69实测地震波——东南向40度作用时间(秒)加速度方向：S50E，记录时长：40.00秒实测地震波——西南向50度作用时间(秒)加速度方向：S40W，记录时长：40.00秒实测地震波——竖向作用时间(秒)加速度方向：VERT，记录时长：40.00秒保留的旧版地震波库18.3。如何人工定义地震波在当前的工程目录下建立相应的地震波文件。文件名应采用“USER”加上“1”或“2”或其他阿拉伯数字。使用“.X”、“.Y”和“.Z”文件后缀给出主方向、次方向和竖向所对应的地震波波形。如果用户给出了无后缀的文件，则认为该文件中的内容为主方向的地震波波形。例如“USER1”、“USER2.X”、“USER2.Y”、“USER2.Z”等文件名都是合法的。文件中第一行输入用户地震波步数N；在第2~第N+1行写入地震波加速度值，单位任意，但要一致。人工输入地震波选择19.1。弹塑性分析目的、意义19.2。弹塑性分析的规范规定19.3。简化弹塑性分析方法及适用范围19.4。静力弹塑性分析方法19.5。动力弹塑性分析方法19。罕遇地震下三种薄弱层弹塑性变形验算方法及其适用范围19.1。弹塑性分析目的、意义三水准设防中的“大震不倒”。两阶段设计中的“第二阶段弹塑性变形验算”。强震下变形验算的基本问题：计算和确定薄弱层位移反应和变形能力；通过改善结构均匀性、加强薄弱层和薄弱部位使得层间位移角满足弹塑性变形验算限值要求。19.2。弹塑性分析的规范规定《建筑抗震设计规范》GB50011-2001《高层混凝土结构技术规程》JGJ3-2002《高层民用建筑钢结构技术规程》JGJ99-98《建筑抗震设计规范》3.4.3条竖向不规则结构应（宜）进行弹塑性变形分析3.6.2条弹塑性分析可以根据具体情况采用弹塑性静力、时程、简化方法5.5.2条何种结构需要进行弹塑性变形验算5.5.3条弹塑性变形验算方法5.5.4条弹塑性分析的简化方法5.5.5条弹塑性层间位移角限值《高层混凝土结构技术规程》4.6.4条,4.6.5条,5.1.13条，4.6.4条有具体规定基本遵从于《建筑抗震设计规范》《高层民用建筑钢结构技术规程》5.3.6条～5.3.10条、5.4.4条、5.5.3条，有具体规定，有层间侧移延性比规定19.3。简化弹塑性分析方法及应用范围新抗震规范5.5.3条规定，罕遇地震下薄弱层（部位）弹塑性变形验算可采用下列方法：“不超过12层且层刚度无突变的钢筋混凝土框架结构、单层钢筋混凝土柱厂房可采用5.5.4条的简化计算方法。”新抗震规范5.5.4条规定的简化弹塑性分析方法包含两方面内容：薄弱层按照楼层区分强度系数确定。弹塑性层间位移角由罕遇地震弹性层间位移角折减得到。可以看出，简化的弹塑性分析方法：有明确的适用范围，超出此范围不能采用。薄弱层的判断和相应弹塑性层间位移角的确定均是估算结果。19.4。静力弹塑性分析方法抗倒塌分析图静力弹塑性分析方法的特点静力弹塑性分析方法是将动力地震作用静力化的一种罕遇地震分析方法。考虑结构的弹塑性性质。较动力弹塑性分析方法能一定程度上节省计算时间。通过静力推覆分析过程可以了解结构的抗倒塌能力。通过能力谱方法可以得到结构的罕遇地震下最大弹塑性位移角。能力谱方法存在“以第一振型振动为主、结构可以等效为单自由度体系”等前提假定，能否适用于超高层结构仍然需要探讨；但推覆分析过程有一定的普适性。19.5。动力弹塑性分析方法动力弹塑性分析方法的特点将罕遇地震作用以较为真实的加速度时程方式进行输入。考虑结构的弹塑性性质。对结构没有过多限制其应用范围的基本假定，适用范围广泛，可以认为是一种仿真分析方法。多条地震波分析时，计算时间相对较长。选取不同的地震波进行分析时，计算结果可能差别较大，需要使用者进行合理的判断。20.1。弹塑性分析软件整体功能简介20.2。弹塑性静力分析软件PUSH简介20.3。弹塑性静力分析软件PUSH工程实例20.4。弹塑性静力分析软件PUSH验证20。弹塑性静力分析的正确应用和普及20.1。弹塑性分析软件整体功能简介20.2。弹塑性静力分析软件PUSH简介较为先进的单元类型。先进的弧长法加载策略。非线性方程叠代方法的多种选择。波前法解线性方程。病态方程的特殊解法处理。接力SATWE、PMSAP程序，适应的结构类型广泛。PUSH软件主要参数说明荷载类型：有倒三角形和矩形两种选择，通常可以选择倒三角形。基底剪力与总重量的比值：通过该参数定义侧向荷载的总和，比如填0.5意思是侧向荷载总量最大可以施加到50%的结构总重量。荷载方向与X轴的夹角：一次静力弹塑性分析只在一个方向上施加侧向荷载，该荷载的方向通过荷载正向与X轴正向的夹角决定。单位度。从头运行和接力运行：PUSH软件具有重启动功能。通过该功能可以接力原来的计算结果进行连续计算。停机控制：配合重启动功能进行计算步数选择。材料参数调整：通过该系列参数的调整，可以按照用户实际情况调整混凝土和钢的单轴本构关系曲线形状，考虑混凝土破坏、受到约束等情况。铰的相对刚度界限：为杆件塑性铰的判断条件，例如填入0.8的含义是杆端截面的刚度与初始刚度相比退化了80％时认为出现了塑性铰。墙高斯点破坏条件：剪力墙通过用户给出破坏应变界限来判断破坏情况。杆元细分、墙元细分：用户可以人为的细分杆件单元。楼板考虑为梁翼缘的相对宽度：通过该参数，用户可以考虑楼板对于梁的刚度贡献。20.3。弹塑性静力分析软件PUSH工程实例——高层混凝土结构高层钢结构20.4。弹塑性静力分析软件PUSH工程实例单层钢框架模型11.21.41.61.822.2050100150200250300350abaqusPUSH荷载因子-位移曲线，PUSH与ABAQUS的比较9层钢框架模型11.051.11.151.21.251.30100200300400500600700abaqus-0.01PUSH-0.01abaqus-1PUSH-1荷载因子-位移曲线，PUSH与ABAQUS的比较00.20.40.60.811.20102030405060abaqusepsa混凝土框架模型荷载因子-位移曲线，PUSH与ABAQUS的比较00.20.40.60.811.20100200300400500荷载因子-位移曲线，PUSH与ABAQUS的比较11.11.21.31.41.51.61.71.8050100150abaqusPUSH荷载因子-位移曲线，PUSH与ABAQUS的比较11.522.533.544.5050100150200250abaqus－rebar×1abaqus－rebar×2push－rebar×1push－rebar×2push－rebar×5abaqus－rebar×5荷载因子-位移曲线，PUSH与ABAQUS的比较21.1。弹塑性动力分析软件EPDA简介21.2。弹塑性动力分析软件EPDA工程实例21.3。弹塑性动力分析软件EPDA验证21。弹塑性动力分析的正确应用和讨论21.1。弹塑性动力分析软件EPDA简介单元模型•梁、杆、柱、撑采用纤维束模型•剪力墙采用弹塑性壳单元方程解法•PCG解线性方程•多种解动力微分方程方法•多种解非线性方程方法接力SATWE、PMSAP程序，适用的结构类型广泛弹塑性动力时程分析参数选择EPDA软件主要参数说明“地震波作用方向角(度)”：地震波主方向与结构X轴夹角，如用户希望地震波主方向作用沿着Y轴方向，此处应添“90”。“主分量峰值加速度(cm/s2)”：地震烈度对应的罕遇地震主方向峰值加速度。“次分量峰值加速度(cm/s2)”：地震烈度对应的罕遇地震次方向峰值加速度。“竖直分量峰值加速度(cm/s2)”：地震烈度对应的罕遇地震竖直方向峰值加速度。“混凝土本构关系类型”：用户可以选择“双线性模型”和“三线性模型”两种混凝土本构关系。ε(压)双线性模型ε(压)σ(压)εcE0σtεu强化弱化σ(压)εcσtεuσc退化斜率退化起始界三线性模型混凝土本构关系模型“塑性铰判断方法”：给出了“弹性积分点比例”和“截面刚度退化比例”两种判断杆系构件的塑性铰的方法。“弹性积分点比例”方法是按照构件截面的积分点仍然保持弹性的比例来判断构件的端部是否出现塑性铰。“截面刚度退化比例”是按照结构进入弹塑性状态后的杆端截面刚度与初始截面刚度的比值来判断构件是否形成塑性铰。“塑性铰判断参数”：该参数与“塑性铰判断方法”相对应，填入0.0～1.0之间的一个数值。当通过“弹性积分点比例”判断塑性铰时，如果填入“0.3”表示“只有30％的端截面积分点保持弹性时出现塑性铰”。当选择当通过“截面刚度退化比例”判断塑性铰时，如果填入“0.3”表示“截面刚度退化为初始截面刚度的30％时出现塑性铰”。“动力微分方程组解法”：目前程序提供给用户两种求解动力微分方程组的方法，Wilson-θ法和Newmark-β法。这两种方法的计算结果差别不大，用户根据需要选择。“非线性方程组解法”：程序提供了两种求解非线性方程组的迭代方法，Newton-Raphson迭代和modifiedNewton-Raphson迭代。这两种方法的迭代次数和适用条件是不同的。对于混凝土结构一般建议采用Newton-Raphson迭代进行计算。“非线性迭代步数限值”：该限值规定了非线性迭代的最多次数，当达到该步数限值时，如果还没有收敛，需要缩短步长进行计算。该值不宜取的过大，“10”左右比较合适，否则会明显增加计算时间。“非线性迭代收敛精度”：EPDA程序衡量非线性迭代是否收敛的依据是“不平衡力向量范数”，一般认为0.01～0.001左右的精度是可以满足工程要求的。该值不宜取的过小，否则程序将难以收敛。“最佳收敛次数”：EPDA程序中提供了计算步长自动伸缩求解策略。所谓“最佳收敛次数”是指用户认为多少次迭代收敛是比较合理的。当程序的迭代次数小于该值时，表示计算步长过小，程序会自动增加计算步长；当程序的迭代次数大于该值时，表示计算步长过大，程序会自动缩减计算步长。通常认为一个计算时间步内，结构应该在3～5次迭代后收敛。“不收敛时步长缩减次数限值”：指当结构计算不收敛时，程序自动缩减步长的最多次数。该参数