糖化程过中麦糟的研究和分析
郑翔鹏
（福建省燕京惠泉啤酒股份有限公司，362100）
摘要：通过对麦糟的定性老化淀粉的分析、外观分析以及定量碘值分析，找出麦糟有关指标的变化情况以及在生产过程中控制的范围，同时研究了麦糟与麦汁碘值之间的关系，表针了两者之间存在的关系，为生产过程的控制起到指导的作用。
关键词：麦糟；麦汁；碘值；老化淀粉
前言
利用麦芽所含的各种水解酶，在适宜的条件（温度、PH值、时间等）下，将麦芽辅助原料中的不溶性高分子物质（淀粉、蛋白质、半纤维素及其中间分解产物等），逐步分解为可溶性的低分子物质，这个分解过程称为糖化。糖化中最主要的任务是淀粉的分解，通过α-淀粉酶、β-淀粉酶、R酶、界限糊精酶、α-葡萄糖苷酶和麦芽糖酶等的作用，把淀粉不断降解为麦芽糖、麦芽三糖、葡萄糖等以及低分子糊精。麦芽淀粉在化学结构上分为两种，直链淀粉占麦芽总淀粉的含量的17%~24%，系由60~2000个葡萄糖残基形成α-1，4健链连接的螺旋形不分支长链，其相对分子量波动在1万~50万之间，易溶于水，直链淀粉遇碘液，碘从螺旋形中通过，与直链淀粉之间构成吸附化合物呈蓝色反应，其中呈色反应与其葡萄糖残基的聚合度有关（见表1）；支链淀粉占麦芽总淀粉的含量的76%~83%，其相对分子量在100万~600万之间，需加热才能溶解，支链淀粉遇到碘液，碘只能与支点外部的20几个葡萄糖基结合，故呈红到紫红色反应。在生产与实验室的分析检测中，通常对麦汁糖份组成的分析以及碘试分析来判断淀粉的分解情况，通常状况下，其反应是属正常现象。
但是，对于麦糟的分析往往给于忽略，很多企业由于生产实行全自动化控制，麦糟在经沉淀洗糟完后直接排放到麦糟回收厂。麦糟的组成直接取解于原料的组成、粉碎方法和粉碎度以及糖化方法，麦糟中含有未分解的小颗粒淀粉、分解不完全的胚乳细胞壁物质和蛋白质，假如这些物质含量过大，将构成了一种细腻的粘性沉淀物质，它们充斥麦糟中，将严重地影响麦层的渗透性，影响麦汁过滤的速度，也影响原料的消耗和成品酒的稳定性。一般来讲，麦糟的产率随麦芽谷皮含量、辅料添加比例和粉碎方式不同而有较大的出入，每吨麦芽可产麦糟1.8~2.0立方米，每吨辅料可产糟0.3~0.4立方米，从麦糟的组成我们可以判断糖化过程的各参数的执行的情况，同时，对麦糟中的老化淀粉的分析可进一步的完善糖化的工艺参数和对在原料选择上的优化。
本文对糖化过程中麦糟的一些指标进行定性定量分析，其中选择优化可行的化学分析方法，利用数学统计表明差异性，最终达到指导生产过程控制的目的。
一、实验方法介绍
在本研究中引入定量碘值检测方法，应用于麦汁、糊化醪、麦糟、啤酒。其检测的原理是在经处理的样品中加入酒精，淀粉或高分子糊精会在酒精的作用下沉淀下来，经离心分离后，重新溶解于磷酸盐缓冲溶液中，加入碘液，根据显红糊精和淀粉的分子量和支链度的不同而形成由红到蓝的颜色，通过分
光光度计测定相应的吸光度，以差值表示碘值的大小。
1、仪器
离心机、离心管、震荡器、移液管、加样管、分光光度计、比色皿（3、4cm）、玻璃搅拌棒
2、试剂
95%酒精、0.02M碘液、磷酸缓冲溶液
3、步骤
取10.0（8.0）ml 经处理的样品于离心管中，加入40.0（32）ml的酒精，机械震荡10分钟，后离心5分钟，保留离心的沉淀物，倒去上清夜，加入20.0ml 的缓冲溶液，再离心后取2.0ml 的上清液加8.0ml的缓冲液于比色皿中，以缓冲液为参比在578nm下测定吸光度，即E2；同时E1为空白样，E3为样品处理后加碘液的测定值，（具体说明略）。
4、公式
△E=（E3—E2—0.952E1）×5
根据公式计算的碘值判断各样品情况。
二、麦糟处理分析
对于麦糟的处理分析中，我们引入麦糟中老化淀粉的分析和研究，取一定量的麦糟在酸碱的作用下，通过盐析的作用迫害麦糟中淀粉系统，通过碘试法定性的判断颜色情况来分析是否存在着老化淀粉，再通过定量碘值的方法分析处理完的麦糟醪液的碘值进一步判断老化淀粉的情况；其中加入麦糟外观分析以判断粉碎情况。
1、麦糟外观分析
取同种品种不同批次的麦糟20批，称取定量的麦糟，观察粉碎情况，具体的分析如下表2所示。
通过对上不同批次的麦糟外观情况分析可见，大部分是粉碎正常，针对麦芽颗粒度大小问题，可以
与麦芽供应商协调以优化，其中，在适当过程中调节粉碎度，使粉碎更为均匀，也可在糖化过程中添加一定量的酶试剂以改善麦芽的分解情况，同时提高得率。
2、麦糟中老化淀粉定性分析
在酸碱的作用下使细胞壁破裂，通过碘试定性分析残留的老化淀粉情况，以判断糖化的分解程度和生产中的得率问题。依上方法取不同批次的麦糟20批进行分析，具体结果如下说明。
表3 麦糟老化淀粉分析
为了进一步说明问题，我们运用数学统计方法分析100批次麦糟的老化淀粉碘试验结果，统计点为：蓝紫~紫、紫~红色、蓝紫~紫~红色、紫红~红~浅红、红~浅红，具体的结果如下面积法图表示。
通过上图1的分析可见，由于在糖化过程中麦芽质量的不均匀性以及粉碎的不均匀性，造成了麦糟
的差异性，在麦糟老化淀粉的分析上，颜色介于紫~红色~浅红与紫~红色占了61%，加上红~浅红色段的7%，总共68%，也就是说只有68%左右的麦糟在老化淀粉上残留较少，糖化较完全，而剩下的32%糖化较不完全，或说使用率较低。因此，对于麦糟老化淀粉的定性分析在一定程度上可以对生产过程中的糖化程度、粉碎情况等进行一定程度上的分析与考核。
3、麦糟定量碘值分析
麦糟的定量分析是建立在定性分析的基础上，对经中和的后麦糟乳液进行处理后依碘值方法测定，其中我们随机分析了100批次麦糟的情况，根据检测的结果以及麦糟外观情况，我们把碘值分为几个分段进行讨论分析，碘值为2＞，2~3，3~4，4~5，5＜，具体统计如图2说明。
由上图2分析统计结果可见，麦糟的碘值在3~4之间占了48%，而碘值大于4的占了27%，碘值小于3的占了25%。根据我们的经验以及其他一些参数条件，认为：麦糟的碘值在3~4之间表示糖化较完全（相对于麦糟的浸出情况），碘值小于3的表示糖化程度好，麦糟的碘值在4~5之间的表示糖化较不完全，而麦糟的碘值大于5的表示糖化不好。通过上的分析表明了在实际生产过程中的控制均匀性和稳定性只有80%左右。因此根据麦糟的碘值的具体情况，我们可以判断糖化过程的完全度情况，从而达到指导生产与优化工艺参数、稳定控制的目的，其中具体的参数说明略。
三、麦糟碘值与麦汁碘值之间的关系分析
麦糟中存在着些粗蛋白、脂肪、纤维素与一些灰分成分，这些物质在一定程度上表针了麦糟与糖化参数的关系，其中与麦汁碘值之间的关系在一定程度上说明了同个问题，麦汁碘值定量的反应了糖化的程度（碘值＜0.3表示糖化完全，麦糟碘值定量的反应了糖化残留的程度。下面我们分析统计了100批次的麦汁和相应的麦糟的碘值情况，找出相应的关系。
备注：图3的各点表示5个批次的组合。
从上图3可见，麦汁碘值与麦糟碘值图象相符合，图中麦汁碘值最小为0.06，相对应麦糟的碘值是2.48，也是该分析组中麦糟的最小碘值，在较大的麦汁碘值以及相应麦糟碘值分析中（麦汁麦糟碘值组0.35： 4.52；0.34：3.97；）也大体体现相应的关系，而其他分析组在数值上大体也体现麦汁碘值与麦糟碘值的正比关系，其中在图中出现的差异点是实验或生产工艺参数差异性而引起的。
麦汁碘值的分析可在一定程度上对生产过程出现的异常进行防范，如在发酵液中补加淀粉分解酶等，有关具体的分析已有很多文献阐述过，这里不在加以说明。
四、其他因素
对于麦糟中其他因素的分析，如麦糟碘值与啤酒碘值的关系、麦糟老化淀粉的详细分解等等，这里就不加以一一的详细说明。
结束语
麦糟作为一种在啤酒生产中最大量的废弃物，对其的研究在一定程度上能对生产过程中中的一些参数控制起到很好的表征作用，同时其有很高的营养价值，可以通过一定的处理后作为高蛋白饲料，这是外话了。总之，对于麦糟的研究不止上阐述中所涉及的问题，还有其他的一些项目，如蛋白成分、糟皮成分等等，通过对研究麦糟有关指标研究，达到均匀有效的控制生产的一些参数，为稳定产品的质量提供定性定量的指标。
参考材料
1、啤酒工业手册
2、酿造酒工艺学
3、内参
计算书：课程设计计算书（题一） 根据加速度调幅公式：m i a t a a a /)(max ,00*= )/(29002902s mm Gal a m == 得：58/)(72900/)(3500i i t a t a a =*= )(i t a =[0 600 1100 150021002500 2900350 2050
15001000600200 -700 -1300-1700 -2000 -1800-1500 -700-250200 -100 0 0 0]; 所以经调幅后为0a =[0 72.6 133.1 181.5 254.1 302.5 350.9 42.4 248.1 181.5 121 72.6 24.2 -84.7 -157.3 -205.7 -242 -217.8 -181.5 -84.7 -30.3 24.2-12.1 0 0 0 ] 6.7206.72''1''2=-=-U U 5.60 6.721.133''2''3=-=-U U 依次类推可以求出地面运动加速度的差值。 因为km c 2=ζ,08.0=ζ , m kN k /9000=, m s kN m /2502?= 代入可以算得m s kN c /240?= 一、表格第一行数据计算： t c t m k K i i /3/62++=*, t=0.05s 代入得m N K i /623400 =* ）△△2 /3()3/6(''''''''t U U c U t U U m P i i g i *++---=* N 18150-6.72250-=*= **=i i P U K △△ mm K P U i i 03.0623400/18150 /-=-==**△△ 起始时刻时：0=U 0'=U 0''=U 因为'''2''3/6/6i i U t U t U U -*-*=△△ 所以7205.0/)03.0(62''1 -=-*=U △
型钢框架混凝土核心筒和钢框架支撑核心筒结构弹塑性时程分析
第十届中日建筑结构技术交流会南京 型钢框架一混凝土核心简和钢框架一支撑核心 简 结构弹塑性时程分析 王斌张翠强吕西林 同济大学土木工程防灾国家重点实验室同济大学结构工程与防灾研究所 AbstraCt Currently noIllinear time llistory amlySis of seismic analysis of mgh-rise buildings has b een widely use 也 but itS amlysis methodS still rleed deVel 叩ment and improVement ． Sino-Japanese S 饥Jctural Engine 甜ng Con6毒rence decided to iIlitiate a nonlinear time histo 呵analysis conlpmtive study in 20 l 2，and t11en organized eight corplofatio 璐at home and abroad for the same case study ．In this paper ，the two cases ， steel reinf ．orced coIlcrete 丘arr 圮-concrete tIl_be smlcture and steel 台arIle-braced n|be s 仃uctllre ，were analyzed based on so 胁are NosaCAD20 l 0 and Midas Building respectiVely ．The nonliTlear time 11istoD ，analysis with 7孕ound motio 璐、Ⅳas 训ed out under me rarely ear 廿1quake with inteIlsit)，8．The def ．0丌】[】ation and damage deVel 叩ment of the s 虮lctllre we 陀stlldied ． Key 帅rds Hybrid stru 【c 咖； noIllinear ti 眦llisto 巧amlysis ； s eisIIlic perf ．0nmnce 1引言 2012年中日建筑结构技术交流会中日双方研究决定进行中日高层建筑结构弹塑性时程的算例对 比分析活动，组织了国内外8家单位对相同案例进行分析比较【l 】。本文针对此次分析活动中2个案例： 钢框架．混凝土核心筒和钢框架．支撑核心筒结构，分别采用NosaCAD2010和Midas Building 有限元分 析程序建立整体结构模型。其中压弯构件采用纤维模型，梁采用塑性铰模型，支撑采用塑性 铰模型， 墙体采用非线性平板壳单元，以反映构件非线性复杂受力情况。通过8度罕遇烈度下7条地震输入的 弹塑性时程分析，研究了该案例结构的变形和破坏情况，探讨了弹塑性时程分析在实际工程中的应用 要点。 第一部分：型钢框架一混凝土核心筒结构 2．1工程概况 钢框架．混凝土核心筒结构共32层，结构总高度129m ，平面基本尺寸为48m×48m，首层5．0m ， 其它层高均为4．Om 。楼板无大开洞，形成刚性横隔板，把核心筒与外框架联系在一起。核心筒采用普 通钢筋混凝土剪力墙，外框架由型钢混凝土柱和钢梁构成的组合结构框架，标准柱距为9．6、米，矩形 型钢混凝土柱直径从基底逐渐减少并延伸至屋顶，外框架梁采用焊接H 型钢梁与柱刚接，与核心筒墙 体铰接，其典型楼层布置和立面见图l 所示。 2．2计算分析程序和主要参数 采用No 鼢CAD20lO 分析程序对该结构进行弹塑性时程分析，对该结构抗震性能和抗震机理进行 研究。 2．2．1构件有限元模型 ：， 梁柱杆单元采用三段变刚度杆单元模型，由位于中部的线弹性区段和位于杆两端的弹塑性段组成。 以受弯为主的钢梁和混凝土梁单元截面的弹塑性段弯矩一曲率骨架曲线分别采用二折线和三折线模型。 由于柱受双向弯矩作用，并到受轴力变化影响，柱单元弹塑性段采用纤维模型，钢和钢筋纤维采用理 235
一、YJK转ABAQUS 1、YJK模型的合理简化 ⑴YJK的模型，如果存在次梁布置不规则、次梁与核心筒搭接不规则、次梁与核心筒开洞相交等情况，会造成模型转化失败，因此，转之前需对模型进行一些合理的简化，既要避免模型转化失败，同时尽可能保持原有模型的特性，防止简化过多，造成简化的模型与原模型在结构动力特性上差别较大，总之一句话，模型简化坚持“简单但不失真”的原则。 此过程不可能一蹴而就，需要反复尝试，简化从少入多，简化越少越好。 ⑵验证简化模型的有效性。 模型转过来以后并不是万事大吉，还需要对比模型进行检验。首先转成线弹性模型，此模型的目的就是采用ABAQUS分析模型的动力特性，查看YJK与ABAQUS两软件计算所得的质量与周期是否一致。若在误差允许范围内，则可进行下一步操作，反之，则需对简化的YJK模型就行修改。 ⑶模型验证有效后，下一步转成弹塑性时程分析模型。转弹塑性时程分析模型之前，有几个问题需要注意： ①关于楼板 楼板是采用刚性楼板还是采用弹性楼板，取决于楼板有没有缺失，若整层楼板开洞很小，且我们不关注楼板的应力状态，则分析时采用刚性楼板即可，后续abaqus弹塑性时程分析时不对楼板细分，会节约计算成本；反之，若楼板缺失严重，且楼板应力分布是重点关注的东西，则YJK要对板指定弹性板3或弹性板6或弹性模。后续ABAQUS分析时会对板就行细分。板内钢筋根据施工图进行确定，但目前导入ABAQUS却不能查看板内钢筋应力分布情况（此问题有待继续研究）。 ②关于梁柱 ABAQUS采用纤维单元进行模拟。梁柱内钢筋采用等效的矩形钢管进行模拟，后续可以查看钢筋的受压损伤因子与受拉损伤因子。梁柱单元细分数目可取2m。 ③关于材料强度 由于ABAQUS分析未考虑箍筋的作用。因此可通过取材料平均值来适当考虑箍筋对混凝土的约束作用。 ⑷参数设置成功以后即可计算，当然计算之前需对电脑进行设置，保证程序可以自动调入子程序。 ⑸ABAQUS分析结果查看，ABAQUS的默认历史时程输出只有能量的输出，我们关心的顶点时程位移曲线，层间位移角，基底剪力这些需要自己编写命令流输出，以供后续处理。 ⑹弹塑性时程分析报告编写 需要涵盖梁、柱、板、墙以及钢筋在大震下的应力分布情况。
弹性时程分析——YJK软件操作篇
操作菜单1 上部结构计算——弹性时程分析 2 常用活动菜单——计算参数+计算分析 3 结果菜单——WDYDA+层位移+层位移角+层剪力+层弯矩+反应谱对比
计算参数 根据《建筑抗震设计规范》（GB50011-2010）表5.1.2-2，多遇 地震，自动对主次方向的峰值加速度取值1第一级对话框——参数输入-弹性时程分析信息次方向的峰值加速度取值取为默认值时，CQC 法结果是考虑了主 次波组合情况下的计算结果。WZQ 中CQC 法的计算结果始终是单向地震下的分量计算结果，未考虑双向地震组合。所以两份文件的CQC 法计算结果只有在单向地震情况下，次方向的峰值加速度取值取为0时保持一致2只计算主方向地震效应：程序对结构地震波效应的计算结果分为0°与90°两种情况，每种情况又各自有主次两个方向分量的效应。 在后续对弹性时程结果的运用中，次方向的效应一般不会用到3
第二级对话框——地震波选择对话框 1 本级菜单一般条件下无需进行调整 2 查看反应谱——PGA、EPA、加速度谱、速度谱、位移谱
第三级对话框——自动筛选最优地震波组合 1 地震波组合晒选限制条件 ?单条地震波基底剪力满足规范要求——±35% ?地震波组合平均基地剪力满足规范要求——±20% ?平台与第一周期领域平均值筛选—— 《结构时程分析法输入地震波的选择控制指标》——仅供参考！ ①一是同欧洲规范，对地震记录加速度反应谱值在[0.1, Tg]平台段的均 值进行控制，要求所选地震记录加速度谱在该段的均值与设计反应谱 相差不超过10% ②二是对结构基本周期T1附近[T1-DT1，T1+DT2 ]段加速度反应谱均 值进行控制，要求与设计反应谱在该段的均值相差不超过10% ③由于实际结构在大震作用下常进入非线性状态，结构刚度发生退化， 结构基本周期随之不断延长，在选取DT1和DT2时，可使 DT2=0.5s>=DT1。Tol为限值
静力弹塑性分析方法与与动力弹塑性分析方法的优缺点 Pushover）分析法 1、静力弹塑性分析方法（Pushover）分析法优点： （1）作为一种简化的非线性分析方法，Pushover方法能够从整体上把握结构的抗侧力性能，可以对结构关键机构及单元进行评估，找到结构的薄弱环节，从而为设计改进提供参考。 （2）非线性静力分析可以获得较为稳定的分析结果，减小分析结果的偶然性，同时花费较少的时间和劳力，较之时程分析方法有较强的实际应用价值。 2、静力弹塑性分析方法（Pushover）分析法缺点： （1）它假定所有的多自由度体系均可简化为等效单自由度体系，这一理论假定没有十分严密的理论基础。 （2）对建筑物进行Pushover分析时首先要确定一个合理的目标位移和水平加载方式，其分析结果的精确度很大程度上依赖于这两者的选择。（3）只能从整体上考察结构的性能，得到的结果较为粗糙。且在过程中未考虑结构在反复加载过程中损伤的累积及刚度的变化。不能完全真实反应结构在地震作用下性状。 二、弹塑性时程分析法
1、时程分析法优点： （1）采用地震动加速度时程曲线作为输入，进行结构地震反应分析，从而全面考虑了强震三要素，也自然地考虑了地震动丰富的长周期分量对高层建筑的不利影响。 （2）采用结构弹塑性全过程恢复力特性曲线来表征结构的力学性质，从而比较确切地、具体地和细致地给出结构的弹塑性地震反应。 （3）能给出结构中各构件和杆件出现塑性铰的时刻和顺序，从而可以判明结构的屈服机制。 （4）对于非等强结构，能找出结构的薄弱环节，并能计算出柔弱楼层的塑性变形集中效应。 2、时程分析法缺点： （1）时程分析的最大缺点在于时程分析的结果与所选取的地震动输入有关，地震动时称所含频频成分对结构的模态n向应有选择放大作用，所以不同时称输入结果差异很大。 （2）时程分析法采用逐步积分的方法对动力方程进行直接积分，从而求得结构在地震过程中每一瞬时的位移、速度和加速度反应。所以此法的计算工作十分繁重，必须借助于计算机才能完成。而且对于大型复杂结构对计算机要求更高，耗时耗力。 （3）对工程技术人员素质要求较高，工程应用要求较高。从结构模型建立，材料本构的选取、地震波选取，到参数控制及庞大计算结果的整理及甄别都要求技术人员具有扎实的专业素质以及丰厚的工程经验。
目录 1 工程概况 (64) 1.1工程介绍 (64) 1.2进行罕遇地震弹塑性时程分析的目的 (64) 2分析方法及采用的计算软件 (65) 2.1分析方法 (65) 2.2分析软件 (65) 2.3材料模型 (65) 2.3.1 混凝土材料模型 (65) 2.3.2 钢材本构模型 (66) 2.4构件模型 (66) 2.4.1 梁单元 (66) 2.4.2 楼板模型 (67) 2.5分析步骤 (67) 2.6结构阻尼选取 (67) 3 结构抗震性能评价指标 (68) 3.1结构的总体变形 (68) 3.2构件性能评估指标 (68) 4 动力特性计算 (69) 5 施工加载过程计算 (69) 5.1施工阶段设置 (69) 5.2施工阶段计算结果 (69) 6 罕遇地震分析总体信息结果汇总 (70) 6.1地震波选取 (70) 6.2基底剪力 (72) 6.3层间位移角 (74) 6.3.1 左塔楼 (74) 6.3.2 右塔楼 (78) 6.4结构顶点水平位移 (82) 6.5柱底反力 (85) 6.8结构弹塑性整体计算指标评价 (86) 7构件性能分析 (87) 7.1钢管混凝土柱 (87) 7.2斜撑 (87) 7.3连梁 (88) 7.3主要剪力墙 (89) 7.4钢梁的塑性应变 (96) 7.5楼板应力及损伤 (96) 8 罕遇地震作用下结构性能评价 (99)
1 工程概况 1.1 工程介绍 上海中心，地下5层，地上33层，结构总高度为180m；主体结构采用框架-核心筒体系，外框架为圆钢管混凝土柱、钢框架梁。 钢管混凝土柱截面为Φ1200x1140～Φ900x860。核心筒采用钢筋混凝土剪力墙体系，外墙厚750mm～400mm，内墙厚500mm～300mm，部分墙体内配置10mm厚钢板。在32层以下，结构由左右两个塔楼构成，中间通过钢梁及6-7层、17-20层两道“人”字形斜撑连接，斜撑截面为BOX 560x1060x80x80。 上部主体结构分析时，以地下室顶板为嵌固端。 图1.1 工程整体效果图（中间一栋） 主要构件信息： （1）框架柱均采用圆钢管混凝土柱，混凝土强度等级为C60。钢管为Q390。 （2）核心筒内连梁： ?上下纵筋配筋率各为1.0%； ?SATWE模型中有钢板的连梁需要考虑内嵌钢板（钢板尺寸20x600）； ?核心筒内其他主梁:上下纵筋配筋率各为1.0%； （3）楼板（C40）:单向配筋率为0.3%。 （4）剪力墙（C60）： ?加强区（66m标高以下及巨型支撑层上下层（含支撑层））： ?暗柱纵筋配筋率为10%（含型钢）； ?墙体的竖向和水平分布筋配筋率均为0.6%； ?其他区域（66m标高以上）： ?角部及与巨型支撑连接处的暗柱纵筋配筋率为5%，其他暗柱1.6%； ?墙体的竖向和水平分布筋配筋率均为0.35%； 图1.2 标准层结构布置图 图1.3 abaqus整体模型图1.4 桁架层 图1.5 典型楼板单元剖分 1.2 进行罕遇地震弹塑性时程分析的目的 对此工程进行罕遇地震作用下的弹塑性时程分析，以期达到以下目的：
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：7800 mechanical-->>elasticity-->>elastic，young‘s modulus：2.1e11，poisson’s ratio：0.3.
弹性时程分析 上部结构计算包含了3个主菜单：前处理及计算、设计结果、弹性时程分析。图为弹性时程分析的各个菜单 一、弹性时程分析计算的目标 《抗规》5.1.2-3条：（及《高规》4.3.5条） 特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑，应采用时程分析法进行多遇地震下的补充计算；当取三组加速度时程曲线输入时，计算结果宜取时程法的包络值和振型分解反应谱法的较大值；当取七组及七组以上的时程曲线时，计算结果可取时程法的平均值和振型分解反应谱法的较大值。 时程分析法是针对特别不规则结构、特别重要结构、较高结构的补充计算，对于时程分析结果的应用，如《抗规》5.1.2条文说明：应把时程法计算结果的底部剪力、楼层剪力、层间位移和上部结构计算的振型分解反应谱法的结果进行比较，当时程分析法大于振型分解反应谱法时，相关部位的构件内力和配筋作相应的调整。 二、对用户选用的地震波提供符合规范要求的的检测 《抗规》5.1.2-3条：采用时程分析法时，应按建筑场地类别和设计地震分组选用实际强震记录和人工模拟的加速度时程曲线，其中实际强震记录的数量不应少于总数的2/3，多组时程曲线的平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符，其加速度时程的最大值可按表5.1.2-2采用。弹性时程分析时，每条时程分析曲线所得结构底部建立不应小于振型分解反应谱法的65%，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的的80%。 什么是“在统计意义上相符”，如《抗规》5.1.2条文说明：多组时程波的平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在对应于结构主要振型的周期点上相差不大于20%。计算结果在结构主方向的平均底部剪力一般不应小于振型分解反应谱法计算结果的80%，每条地震波输入的计算结果不应小于65%。但计算结果也不能太大，每条地震波输入计算不大于135%，平均不大于120%。 根据如上规范要求，软件将对用户选择的每一条地震波进行检测，如下图的超限提示是使用该波计算的基底剪力超过要求。 下图的作用是，比较多组时程波的平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，对应于结构主要振型的周期点上相差是否大于20%。如果大于20%则给出红色的超限提示。
弹性、弹塑性时程分析法在结构设计中的应用 杨志勇黄吉锋 (中国建筑科学研究院北京 100013 0 前言 地震作用是建筑结构可能遭遇的最主要灾害作用之一。几十年来,人们积累了大量的实测地震资料,这些资料多以位移、速度或者加速度时程的形式体现。与此相对应,时程分析方法也被认为是最直接的一种计算建筑结构地震响应的方法。但是,由于地震作用随机性导致计算结果的不确定性,弹性时程分析方法只是结构设计的一种辅助计算方法;虽然如此,抗震规范为了增强重要结构的抗震安全性,还是将弹性时程分析方法规定为常遇地震作用下振型分解反应谱法的一种补充计算方法;尤其是考虑了结构的弹塑性性能后,弹塑性时程分析方法更是被普遍认为是一种仿真的罕遇地震作用响应计算方法。 《建筑抗震设计规范》 (GB50011-2001第3.6.2,5.1.2, 5.5.1,5.5.2,5.5.3等条文规定了时程分析相关的内容。下面结合TAT,SATWE,PMSAP和EPDA等软件应用,探讨如何将弹性、弹塑性时程分析正确应用到结构设计中去。 1 弹性时程分析的正确应用 正确地在软件中应用弹性时程分析方法需要对规范的相关条文规定有正确的认识。以下几点是需要特别明确的: (1抗震规范第5.1.2条第3点规定,“可取多条时程曲线计算结果的平均值与振型分解反应谱法计算结果的较大值”。在设计过程中,如何实现“较大值”有不同的做法: 1设计采用弹性时程分析的构件内力响应包络值的多波平均值与振型分解反应谱法计算结果二者的较大值直接进行构件设计;2在实现振型分解反应谱方法时,放大地震力使得到的楼层响应曲线包住时程分析楼层响应曲线的平均值。
弹塑性时程分析方法将结构作为弹塑性振动体系加以分析，直接按照地震波数据输入地面运动，通过积分运算，求得在地面加速度随时间变化期间内，结构的内力和变形随时间变化的全过程，也称为弹塑性直接动力法。 基本原理 多自由度体系在地面运动作用下的振动方程为： 式中、、分别为体系的水平位移、速度、加速度向量；为地面运动水平加速度，、、 分别为体系的刚度矩阵、阻尼矩阵和质量矩阵。将强震记录下来的某水平分量加速度－时间曲线划分为很小的时段，然后依次对各个时段通过振动方程进行直接积分，从而求出体系在各时刻的位移、速度和加速度，进而计算结构的内力。 式中结构整体的刚度矩阵、阻尼矩阵和质量矩阵通过每个构件所赋予的单元和材料类型组装形成。动力弹塑性分析中对于材料需要考虑包括：在往复循环加载下，混凝土及钢材的滞回性能、混凝土从出现开裂直至完全压碎退出工作全过程中的刚度退化、混凝土拉压循环中强度恢复等大量非线性问题。 基本步骤 弹塑性动力分析包括以下几个步骤： (1) 建立结构的几何模型并划分网格； (2) 定义材料的本构关系，通过对各个构件指定相应的单元类型和材料类型确定结构的质量、刚度和阻尼矩阵； (3) 输入适合本场地的地震波并定义模型的边界条件，开始计算； (4) 计算完成后，对结果数据进行处理，对结构整体的可靠度做出评估。 计算模型 在常用的商业有限元软件中，ABAQUS、ADINA、ANSYS、MSC.MARC都内置了混凝土的本构模型，并提供了丰富的单元类型及相应的前后处理功能。在这些程序中一般都有专用的钢筋模型，可以建立组合式或整体式钢筋。 以ABAQUS为例，它提供了混凝土弹塑性断裂和混凝土损伤模型以及钢筋单元。其中弹塑性断裂和损伤的混凝土模型非常适合于钢筋混凝土结构的动力弹塑性分析。它的主要优
第九章时程分析法 第一节时程分析法的概念 振型分解法仅限于计算结构在地震作用下的弹性地震反应。时程分析法是用数值积 分求解运动微分方程的一种方法，在数学上称为逐步积分法。这种方法是从t=0时刻开始，一个时段接着一个时段地逐步计算，每一时段均利用前一时段的结果，而最初时段应根 据系统的初始条件来确定初始值。即是由初始状态开始逐步积分直至地震终止，求出结 构在地震作用下从静止到振动、直至振动终止整个过程的地震反应。 时程分析法是对结构动力方程直接进行逐步积分求解的一种动力分析方法。时程分 析法能给出结构地震反应的全过程，能给出地震过程中各构件进入弹塑性变形阶段的内 力和变形状态，因而能找出结构的薄弱环节。 时程分析法分为弹性时程分析法和弹塑性时程分析法两类。 第一阶段抗震计算“小震不坏”中，采用时程分析法进行补充计算，这时计算所采用 的结构刚度和阻尼在地震作用过程中保持不变，称为弹性时程分析。 在第二阶段抗震计算“大震不倒”中，采用时程分析法进行弹塑性变形计算，这时结 构刚度和阻尼随结构及其构件所处的非线性状态，在不同时刻可能取不同的数值，称为 弹塑性时程分析。弹塑性时程分析能够描述结构在强震作用下在弹性和非线性阶段的内力、变形，以及结构构件逐步开裂、屈服、破坏甚至倒塌的全过程。 第二节时程分析法的适用范围 一、时程分析法的适用范围 时程分析法是根据选定的地震波和结构恢复力特性曲线，对动力方程进行直接积分，采用逐步积分的方法计算地震过程中每一瞬时的结构位移、速度和加速度反应，从而可观察到结构在强震作用下弹性和非弹性阶段的内力变化以及构件开裂、损坏直至结构倒塌的全过程。但此法的计算工作十分繁重，须借助计算机，费用较高，且确定计算参数尚有许多困难，目前仅在一些重要的、特殊的、复杂的以及高层建筑结构的抗震设计中应用。《建筑抗震设计规范》对时程分析法的适用范围规定如下：
建 筑 结 构 学 报(增刊1) Journal of Building Structures (Supp le mentary Issue 1) 弹性与弹塑性动力时程分析方法中若干问题探讨 杨志勇,黄吉锋,邵 弘 (中国建筑科学研究院结构所,北京100013) 摘要:依据大量实际工程弹性、弹塑性动力时程分析经验,结合实际工程应用,探讨了弹性、弹塑性动力时程分析方法中的一些基本问题。针对性地分析了动力时程分析方法中地震波的离散性;地震波如何与反应谱曲线在统计意义上相符;人工模拟地震波方法及其工程应用;弹性、弹塑性时程分析法选取地震波的基本原则;弹性时程分析法地震波的选取数量;如何将反应谱分析结果与时程分析结果取较大值等方面的问题。通过大量的算例分析可以看出,正确地应用弹性、弹塑性动力时程分析方法需要从多个方面进行准确理解和把握,教条地应用很难发挥弹性、弹塑性动力时程分析应有的作用。关键词:弹性时程分析法;弹塑性时程分析法;地震波;反应谱中图分类号:T U31113 文献标识码:A D iscussi on on linear and nonlinear ti m e hist ory analysis method Y ANG Zhiyong,HUANG Jifeng,SHAO Hong (Building Structure Research I nstitute,China Academy of Building Research,Beijing 100013,China ) Abstract:This paper discussed linear and nonlinear ti m e hist ory analysis method,es pecially concerning with the following issues:the disperse of earthquake wave,scaling the earthquake wave t o fit the design res ponse s pectrum of China code,the earthquake wave si mulati on method,the basic p rinci p le of earthquake wave selection,the number of waves required in ti me hist ory analysis,and the maxi mum structural res ponse fro m s pectrum analysis and ti me hist ory analysis .A s sho wn in many examp les,linear and nonlinear ti m e hist ory analysis method should be used app r op riately t o obtain useful results . Keywords:linear ti me history analysis method;nonlinear ti me history analysis method;earthquake wave;res ponse s pectrum 基金项目:建设部软科学研究资助项目(062K9231)。 作者简介:杨志勇(1974— ),男,黑龙江齐齐哈尔人,工学博士,副研究员。收稿日期:2008年6月 0 前言 《建筑抗震设计规范》(G B 50011—2001)、《高层建 筑混凝土结构技术规程》(JGJ 3—2002)、《高层民用建筑钢结构技术规程》(JGJ 99—98)等对于弹性、弹塑性 动力时程分析方法进行了具体的规定,涉及到弹性、弹塑性时程分析方法适用范围,地震波的选取原则,变形 验算的限值规定等方面[123] 。随着复杂、超限结构的增多,弹性、弹塑性动力时程分析方法在实际建筑结构抗震设计中得到了越来越多的工程应用。通过对一定数量的实际工程弹性、弹塑性动力时程分析实例的参与,发现在实际应用中存在着较多方面的问题,对其中的一些重要问题做一总结和探讨,为弹性、弹塑性动力时程分析方法的进一步完善提供量化依据。 1 地震波的离散性 图1所示为一幢17层高层混凝土结构模型,该结构有2层地下室,抗震设防烈度为8度,Ⅱ类场地,多遇地震特征周期0135s 。图2、图3给出了该结构4条天然波和4条人工波的多遇地震弹性时程分析法和反应谱分析法计算得到的顶点位移、基底剪力响应结果,这些地震波来源于PKP M 软件的地震波数据库。表1为多遇地震时8条地震波的弹性时程分析与反应谱分析响应的对比结果,对比曲线见图2。表2为多遇地震弹性时程分析法中地震波离散性的分析结果,对比曲线见图3。 3 12
时程分析法 定义：由结构基本运动方程沿时间历程进行积分求解结构振动响应的方法。 概述：时程分析法是世纪60年代逐步发展起来的抗震分析方法。用以进行超高层建筑的抗震分析和工程抗震研究等。至80年代,已成为多数国家抗震设计规范或规程的分析方法之一。 原理：时程分析法在数学上称步步积分法,抗震设计中也称为“动态设计”。由结构基本运动方程输入地面加速度记录进行积分求解，以求得整个时间历程的地震反应的方法。此法输入与结构所在场地相应的地震波作为地震作用,由初始状态开始, 一步一步地逐步积分,直至地震作用终了。 是对工程的基本运动方程，输入对应于工程场地的若干条地震加速度记录或人工加速度时程曲线，通过积分运算求得在地面加速度随时间变化期间结构的内力和变形状态随时间变化的全过程，并以此进行结构构件的界面抗震承载力验算和变形验算。 时程分析法是世纪60年代逐步发展起来的抗震分析方法。用以进行超高层建筑的抗震分析和工程抗震研究等。至80年代,已成为多数国家抗震设计规范或规程的分析方法之一。 “时程分析法”是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的一种结构动力计算方法，也为国际通用的动力分析方法。 “时程分析法”常作为计算高层或超高层的一种（补充计算）方法，也就是说满足了规范要求的时候是可以不用它计算结构的。规范规定：对于特别不规则的建筑、甲类建筑及超过一定高度的高层建筑，宜采用时程分析法进行补充计算。所以有较多设计人员对应用时程分析法进行抗震设计感到生疏。近年来,随着高层建筑和复杂结构的发展,时程分析在工程中的应用也越来越广泛了。 地震动输入对结构的地震反应影响非常大。目前的现状是，输入地震动的选择大多选择为数不多的几条典型记录（如：1940年的El Centro（NS）记录或1952年的Taft记录），国内外进行结构时程分析时所经常采用的几条实际强震记录主要有适用于I类场地的滦河波、适用于II、III类场地的El-Centrol波(1940，N-S)和Taft波（1952，E-w）、适用于IV 类场地的宁河波等。
80 第40卷 增刊 建 筑 结 构 2010年6月 北京某超高层商住楼动力弹塑性时程分析 徐晓龙，高德志，桂满树，姜毅荣，何四祥，王 侃 （北京迈达斯技术有限公司，北京 100044） [摘要] 基于梁柱塑性铰和剪力墙纤维模型，利用MIDAS Building 软件实现了超高层建筑结构的弹塑性时程分析。结合该结构研究了在大震作用下结构将出现的破坏模式、塑性发展特点等，并与弹性分析进行了对比，说明弹塑性分析更能反映实际情况，能对结构的抗震性能给出较为合理全面的评价，并对工程设计给出指导。 [关键词] 动力弹塑性时程分析；MIDAS Building ；纤维模型 Elastic-plastic time-history analysis on the super-high business-living building in Beijing Xu Xiaolong, Gao Dezhi, Gui Manshu, Jiang Yirong, He Sixiang, Wang Kan （Beijing MIDAS Technology Information Co.,Ltd,. Beijing 100044,China ） Abstract: Based on the theory of plastic hinges （beams and columns ） and fiber model （walls ）, elastic-plastic time-history analysis is performed on the super-high business-living building in Beijing by MIDAS Building software under the scarce earthquake load. Failure Modes and plastic zone development are researched according to the feature of the structure. Through the comparison with the elastic analysis, it is considered that evaluation on the structure can be deduced from the elastic-plastic analysis more reasonably and comprehensively, and there will be better instruction to the projects. Keywords: dynamic elastic-plastic analysis; MIDAS Building; fiber model 1 结构特点 某50层的超高层商住两用建筑，地上50层，结构高度达到236.3m ，采用钢骨混凝土柱框筒结构形式，平面尺寸64.8m ×43.8m （轴线尺寸）。结构已经超过型钢混凝土框架-钢筋混凝土筒体结构8度（0.2g ）抗震设防下的最大适用高度（150m ），该结构为抗震超限结构，故有必要对结构进行动力弹塑性时程分析，以考察其在罕遇地震作用下的响应、薄弱环节、破坏模式等。结构整体模型及首层平面见图1，2。 2 动力弹塑性时程分析 图1 结构模型图 图2 首层平面图 时程分析法[1]被认为是目前结构弹塑性分析的最可靠和最精确的方法，它不仅能对结构进行定性分析，同时又可给出结构在罕遇地震下的量化性能指标，并且得到结构在各个时刻的真实地震反应。弹塑性时程分析方法将结构作为弹塑性振动体系加以分析，直接按照地震波数据输入地面运动，通过逐步积分运算，求得在地面加速度随时间变化期间内，结构的内力和变形随时间变化的全过程，也称为弹塑性直接积分法。 弹塑性动力时程分析有如下优点：1）输入的是罕遇地震波的整个过程，可以真实反映各个时刻地震作用引起的结构响应，包括变形、内力、损伤状态（开裂和破坏）等；2）有些程序通过定义材料的本构关系来考虑结构的弹塑性性能，故可以准确模拟任何结构，计算模型简化较少；3）该方法基于塑性区的概念，对带剪力墙的结构，结果更为准确可靠。 基于MIDAS Building 动力弹塑性分析平台，对北京某超高层商住楼进行了罕遇地震作用下的动力时程分析，研究其各个抗震性能指标以及破坏模式。 2.1 弹塑性动力分析的基本方法 弹塑性动力分析包括以下几个步骤：1）建立结构
动力弹塑性时程分析的方法及其应用 彪仿俊1 阎晓铭1 陈志强1王传甲1王庆扬1，2张劲2 （1 深圳市电子院设计有限公司；2 中国石油大学） 摘要:本文对现有的弹塑性分析方法进行了概述，重点介绍了动力弹塑性时程分析的理论、优点和基本方法，及该方法在东莞一实际工程中的成功应用，对于动力弹塑性时程分析方法在高层、特别是超限高层分析中的推广应用提供了有益的参考和借鉴。 关键词: 静力弹塑性分析动力弹塑性时程分析 ABAQUS 混凝土塑性损伤模型 1.引言 《建筑抗震设计规范》5.5.2条规定，对于特别不规则的结构、板柱－抗震墙、底部框架砖房以及高度不大于150m的高层钢结构、7度三、四类场地和8度乙类建筑中的钢筋混凝土结构和钢结构宜进行弹塑性变形验算。对于高度大于150m的钢结构、甲类建筑等结构应进行弹塑性变形验算。《高层建筑混凝土结构技术规程》5.1.13条也规定，对于B级高度的高层建筑结构和复杂高层建筑结构，如带转换层、加强层及错层、连体、多塔结构等，宜采用弹塑性静力或动力分析方法验算薄弱层弹塑性变形。 历史上的多次震害也证明了弹塑性分析的必要性：1968年日本的十橳冲地震中不少按等效静力方法进行抗震设防的多层钢筋混凝土结构遭到了严重破坏，1971年美国San Fernando地震、1975年日本大分地震也出现了类似的情况。相反，1957年墨西哥城地震中11～16层的许多建筑物遭到破坏，而首次采用了动力弹塑性分析的一座44层建筑物却安然无恙，1985年该建筑又经历了一次8.1级地震依然完好无损。 可以看出，随着建筑高度迅速增长，复杂程度日益提高，完全采用弹性理论进行结构分析计算和设计已经难以满足需要，弹塑性分析方法也就显得越来越重要。2.现有弹塑性分析方法综述 2.1 静力弹塑性分析 静力弹塑性分析（PUSH-OVER ANAL YSIS，以下简称POA）方法也称为推覆法，它基于美国的FEMA-273抗震评估方法和ATC-40报告，是一种介于弹性分析和动力弹塑性分析之间的方法，其理论核心是“目标位移法”和“承载力谱法”。 1.计算方法 (1)建立结构的计算模型、构件的物理 参数和恢复力模型等； (2)计算结构在竖向荷载作用下的内 力； (3)建立侧向荷载作用下的荷载分布形 式，将地震力等效为倒三角或与第 一振型等效的水平荷载模式。在结 构各层的质心处，沿高度施加以上 形式的水平荷载。确定其大小的原 则是：水平力产生的内力与前一步 计算的内力叠加后，恰好使一个或 一批杆件开裂或屈服； (4)对于开裂或屈服的杆件，对其刚度 进行修改后，再增加一级荷载，又 使得一个或一批杆件开裂或屈服； (5)不断重复步骤（3）、（4），直至结构 达到某一目标位移或发生破坏，将 此时的结构的变形和承载力与允许
从传统的观点来看，底部剪力法，反应谱法和时程分析法是三大最常用的结构地震响应分析方法。那么正确的认识它们的一些关键概念，对于建筑结构的抗震设计具有非常重要的意义。HiStruct在此简单的总结一些，全当抛砖引玉。 1. 底部剪力法 高规规定：高度不超过40m、以剪切变形为主且质量和刚度沿高度分布比较均匀的高层建筑结构，可采用底部剪力法。底部剪力法适用于基本振型主导的规则和高宽比很小的结构，此时结构的高阶振型对于结构剪力的影响有限，而对于倾覆弯矩则几乎没有什么影响，因此采用简化的方式也可满足工程设计精度的要求。底部剪力法尚有一个重要的意义就是我们可以用它的理念，简化的估算建筑结构的地震响应，从而至少在静力的概念上把握结构的抗震能力，它还是很有用的。 2. 反应谱方法 高规规定：高层建筑结构宜采用振型分解反应谱法。对质量和刚度不对称、不均匀的结构以及高度超过100m的高层建筑结构应采用考虑扭转耦联振动影响的振型分解反应谱法。反应谱的振型分解组合法常用的有两种：SRSS和CQC。虽然说反应谱法是将并非同一时刻发生的地震峰值响应做组合，仅作为一个随机振动理论意义上的精确，但是从实际上它对于结构峰值响应的捕捉效果还是很不错的。一般而言，对于那些对结构反应起重要作用的振型所对应频率稀疏的结构，并且地震此时长，阻尼不太小（工程上一般都可以满足）时，SRSS是精确的，频率稀疏表面上的反应就是结构的振型周期拉的比较开；而对于那些结构反应起重要作用的振型所对应的频率密集的结果(高振型的影响较大，或者考虑扭转振型的条件下)，CQC是精确的。这是因为对于建筑工程上常用的阻尼而言，振型相关系数（见高规3.3.11-6）在很窄的范围内才有显著的数值。 3.反应谱分析的精确性 对于采用平均意义上的光滑反应谱进行分析而言，其峰值估计与相应的时程分析的平均值相比误差很小，一般只有百分之几，因此可以很好的满足工程精度的要求，正是在这个平均（普遍性）意义上，我们认为反应谱分析方法是精确的。但是对于单个锯齿形的反应谱而言，其分析结果与单个波的时程分析，误差可以达到10-30%之间，因此在个别（特殊性）意义上而言，反应谱分析结果是有误差的，因此，规范规定对于复杂的或者高层建筑需要采用时程分析进行补充计算和验证。 4.反应谱分析与时程分析对于高阶振型计算的不同之处 一般反应谱的高频段是采用平台段来表达的，实际上对于高阶振型反应不显著的结构而言，反应谱适用性很好，也足够准确。但是对于高柔结构而言，一般高阶振型的影响比较显著，采用时程分析的时候，等于其高频段的峰值并未被人为削成平台段，因此采用时程分析的时候此频段的地震响应可能很大，一般表现为高层建筑的顶部或者对其他结构对高阶振型影响显著部位，其地震响应峰值比反应谱分析结果要大（但是总体的剪力和弯矩差别则没这么明显）。 5.时程分析 理论上时程分析是最准确的结构地震响应分析方法，但是由于其分析的复杂性，且地震波的随机性，因此一般只是把它作为反应谱的验证方法而不是直接的设计方法使用。高规规定：3 7～9度抗震设防的高层建筑，下列情况应采用弹性时程分析法进行多遇地震下的补充计算：
PKPM软件园地 建筑结构.技术通讯 2007年1月 弹性、弹塑性时程分析法在结构设计中的应用 杨志勇 黄吉锋 （中国建筑科学研究院 北京 100013） 0 前言 地震作用是建筑结构可能遭遇的最主要灾害作用之一。几十年来，人们积累了大量的实测地震资料，这些资料多以位移、速度或者加速度时程的形式体现。与此相对应，时程分析方法也被认为是最直接的一种计算建筑结构地震响应的方法。但是，由于地震作用随机性导致计算结果的不确定性，弹性时程分析方法只是结构设计的一种辅助计算方法；虽然如此，抗震规范为了增强重要结构的抗震安全性，还是将弹性时程分析方法规定为常遇地震作用下振型分解反应谱法的一种补充计算方法；尤其是考虑了结构的弹塑性性能后，弹塑性时程分析方法更是被普遍认为是一种仿真的罕遇地震作用响应计算方法。 《建筑抗震设计规范》（GB50011－2001）第3.6.2，5.1.2，5.5.1，5.5.2，5.5.3等条文规定了时程分析相关的内容。下面结合TAT ，SATWE ，PMSAP 和EPDA 等软件应用，探讨如何将弹性、弹塑性时程分析正确应用到结构设计中去。 1 弹性时程分析的正确应用 11 正确地在软件中应用弹性时程分析方法需要对规范的相关条文规定有正确的认识。以下几点是需要特别明确的： （1）抗震规范第5.1.2条第3点规定，“可取多条时程曲线计算结果的平均值与振型分解反应谱法计算结果的较大值”。在设计过程中，如何实现“较大值”有不同的做法：1)设计采用弹性时程分析的构件内力响应包络值的多波平均值与振型分解反应谱法计算结果二者的较大值直接进行构件设计；2)在实现振型分解反应谱方法时，放大地震力使得到的楼层响应曲线包住时程分析楼层响应曲线的平均值。 图1 SATWE 地震作用放大系数 前一种做法可能使得构件配筋较大，因为在时程分析过程中，构件内力的最大响应具有不同时性，采用包络值进行设计会使得构件内力，尤其是压弯构件内力偏于保守。因此， TAT ，SATWE ，PMSAP 等软件均提供了地震力放大功能。SATWE 地震作用放大系数见图1，可以通过适当地放大振型分解反应谱法的地震作用来满足相应的规范要求。TAT 软件给出了一种折中的做法，如果设计者进行了弹性时程分析，则程序会将弹性时程分析结果作为一种地震荷载工况进行组合、设计。但是为了避免设计结果过于保守，程序会进行构件弹性时程分析内力的预组合。 （2）“采用时程分析方法时，应按建筑场地类别和设计地震分组选用不少于两组的实际强震记录和一组人工模拟的加速度时程曲线”。建筑结构在不同地震波作用下的响应差别可能较大，选用多条地震波的平均值可在一定程度上避免离散性。人工模拟地震波一般是以规范反应谱为基础，通过蒙特卡罗方法来得到，更加贴近规范反应谱或反映场地土的当地特征。TAT ，SATWE ，PMSAP ，EPDA 等软件按照结构的特征周期给出多组天然波和人工波，见图2。无论是进行弹性还是进行弹塑性时程分析，均要选取足够数量的地震波进行计算，以得到有代表意义的结果。 图2 按照特征周期区分的地震波库 （3）“多波平均地震响应系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符”。其条文说明解释为二者在各个周期点上相差不大于20％。对于人工波来说，这一规定一般是天然满足的，因为人工波是拟合规范反应谱得到的。对于天然波来讲则较难满足，因为规范的反应谱是依据众多实际采集的地震动时程曲线通过平滑化后的概率平均意义上的结果。图3所示为上述软件地震波库0.45s 特征周期中的2条天然波的动力放大系数谱曲线，可见与规范反应谱的差异还是明显的。那么如何执行规范的这条规定呢？其实规范的规定从概念上讲是合理的，因为频谱特征是地震波的最重要特征之一，一定程度上会影响时程分析结果的合理性。一种可行的做法是判断某条实际地震波 第一作者简介：杨志勇，男，1974.6出生，工学博士，副研究员。