发电厂房墙体地震响应的反应谱法与时程分析比较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.19 0.2612.5 0.3130.25 0.047与RG1.60标准谱对应的两条人工波见文件rg160x.txt与rg160y.txt2数值分析框图思路与理论简介2.1理论简介该问题主要牵涉到结构动力分析当中的时程分析和谱分析。
时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。
谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。
2.2 分析框架：时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。
谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。
3有限元模型与荷载说明3.1 有限元模型考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。
然后，采用solid45单元，设置拖拉方向的单元尺寸并清楚初始平面单元plane42，将平面单元进行拖拉，最后生成发电厂墙体的有限元立体几何模型。
单元总数为6060个，总节点数为8174个，有限元模型如图2所示：图2 发电厂墙体有限元模型3.2 荷载说明时程分析：首先计算结构的前两阶自振频率，分别为126.10008.2867f f ==，。
则结构的圆频率1122238.3274,252.0669f f ωπωπ====，对于常阻尼比10.02ζ=，20.05ζ=，30.10ζ=，由结构质量阻尼公式12122ζωωαωω=+得：10.8831α=，2 2.2076α=，3 4.4153α=，由结构刚度阻尼公式122ζβωω=+得，10.0004425β=，20.0011063β=，30.0022125β=。
采用底部完全约束，加载方式为加载地震波的加速度，分为2800个荷载步，每个荷载步取一个子步进行加载求解。
底跨中单宽上的剪力和弯矩是通过面项然后积分的方式求得，最后得出最大值并在excel 中画出时间历程曲线。
谱分析：求得模态后进行两次谱分析，输入的频率值分别为33、9、2.5、0.25，对应谱值为0.1g 、0.261g 、0.313g 、0.047g 。
同样采用底部完全约束，后处理时采用工况记录和运算的方法。
4数值结果与合理性比较（单位均为国际单位制）4.1 时程分析结果①常阻尼比1=0.02ζ：顶部跨中最大水平位移为3max UX 2.1310m -=⨯，沿X 轴正方向。
其时间历程曲线如图3所示：图3 1=0.02ζ时顶部跨中水平位移时间历程曲线底部跨中单宽上X 轴方向剪力最大值为max FSX 91618.72N =，沿X 轴正方向，时间历程曲线如图4所示：图4 1=0.02ζ时底部跨中单宽上X 轴方向剪力时间历程曲线底部跨中单宽上Z 轴方向剪力最大值为max FSZ 63882.6N =，沿Z 轴负方向，其时间历程曲线如图5所示：图5 1=0.02ζ时底部跨中单宽上Z 轴方向剪力时间历程曲线底部跨中单宽上X 轴方向弯矩最大值为max MX 15689.8N m =⋅，沿X 轴负方向，其时间历程曲线如图6所示：图6 1=0.02ζ时底部跨中单宽上X 轴方向弯矩时间历程曲线底部跨中单宽上Z 轴方向弯矩最大值为max MZ 680449N m =⋅，沿Z 轴负方向，其时间历程曲线如图7所示：图7 1=0.02ζ时底部跨中单宽上Z 轴方向弯矩时间历程曲线②常阻尼比2=0.05ζ：顶部跨中最大水平位移为3max UX 1.3710m -=⨯，沿X 轴正方向。
其时间历程曲线如图8所示：图8 2=0.05ζ时顶部跨中水平位移时间历程曲线底部跨中单宽上X 轴方向剪力最大值为max FSX 64976.48N =，沿X 轴正方向，时间历程曲线如图9所示：图9 2=0.05ζ时底部跨中单宽上X 轴方向剪力时间历程曲线底部跨中单宽上Z 轴方向剪力最大值为max FSZ 62084.09N =，沿Z 轴正方向，其时间历程曲线如图10所示：图10 2=0.05ζ时底部跨中单宽上Z 轴方向剪力时间历程曲线底部跨中单宽上X 轴方向弯矩最大值为max MX 8050.85N m =⋅，沿X 轴负方向，其时间历程曲线如图11所示：图11 2=0.05ζ时底部跨中单宽上X 轴方向弯矩时间历程曲线底部跨中单宽上Z 轴方向弯矩最大值为max MZ 448754N m =⋅，沿Z 轴负方向，其时间历程曲线如图12所示：图12 2=0.05ζ时底部跨中单宽上Z 轴方向弯矩时间历程曲线③常阻尼比3=0.10ζ：顶部跨中最大水平位移为4max UX 9.8910m -=⨯，沿X 轴负方向。
其时间历程曲线如图13所示：图13 3=0.10ζ时顶部跨中水平位移时间历程曲线底部跨中单宽上X 轴方向剪力最大值为max FSX 53318.6N =，沿X 轴负方向，时间历程曲线如图14所示：图14 3=0.10ζ时底部跨中单宽上X 轴方向剪力时间历程曲线底部跨中单宽上Z 轴方向剪力最大值为max FSZ 58739.21N =，沿Z 轴正方向，其时间历程曲线如图15所示：图15 3=0.10ζ时底部跨中单宽上Z 轴方向剪力时间历程曲线底部跨中单宽上X 轴方向弯矩最大值为max MX 3976.12N m =⋅，沿X 轴负方向，其时间历程曲线如图16所示：图16 3=0.10ζ时底部跨中单宽上X 轴方向弯矩时间历程曲线底部跨中单宽上Z 轴方向弯矩最大值为max MZ 335338.2N m =⋅，沿Z 轴正方向，其时间历程曲线如图17所示：图17 3=0.10ζ时底部跨中单宽上Z 轴方向弯矩时间历程曲线4.2 谱分析结果①常阻尼比1=0.02ζ：X 轴向：最大水平位移为3max UX 4.7310m -=⨯，最大剪力为max FSX 154276.2N =，最大弯矩为max 9322.MX 1m 7N =⋅；Z 轴向，最大剪力为max FSZ 85609.2N =，最大弯矩为max 2173789MZ m .7N =⋅。
②常阻尼比2=0.05ζ：X 轴向：最大水平位移为x 3ma 3.321UX m 0-⨯=，最大剪力为max 9457FSX N 1.9=，最大弯矩为max MX 11006.1N m =⋅；Z 轴向，最大剪力为max 7694FSZ N 4.7=，最大弯矩为max MZ 857762.2N m =⋅。
③常阻尼比3=0.10ζ：X 轴向：最大水平位移为x 3ma 2.141UX m 0-⨯=，最大剪力为max 7384FSX N 2.3=，最大弯矩为max MX 60792.2N m =⋅；Z 轴向，最大剪力为max 6255FSZ N 8.9=，最大弯矩为max MZ 636487.5N m =⋅。
4.3 结果比较时程分析和谱分析关于顶部水平位移、底部跨中单宽剪力、弯矩最大响应数据结果对比如表1所示：表1 顶部水平位移、底部跨中单宽剪力、弯矩最大响应5结论与体会5.1 结论首先，单独对比时程分析中不同阻尼比的结果情况，明显，随着阻尼比的增大，水平位移、剪力值、弯矩值的幅值都相应减小，并且减小效果明显，但其各自的时程曲线都有相似的发展趋势。
同样，谱分析中的结果也有相类似的效应。
其次，对比时程分析和谱分析的结果得出，在相同阻尼比的条件下，谱分析的最大响应明显比瞬态大，这主要的原因是在计算时程分析中，本人开始计算的质量阻尼和刚度阻尼都很大，直到把所有的结果都整理完后，才掌握了正确的质量阻尼和刚度阻尼计算公式，由于计算和数据处理时间过长，所以在这里没做修正，但结果的对比情况看，都较合理。
最后，通过时程分析和谱分析，本人发现，时程分析非常耗时，占用内存大，而谱分析非常快而且计算的结果可以作为工程数据参考。
所以，本人认为谱分析在某些时候可能更适合工程实践，并作为一种工程结构的地震分析方法。
5.2 疑问时程分析中，其一，关于计算质量阻尼和刚度阻尼所取的自振频率阶数，本报告中我是取的前两阶，但也有同学说去第一阶和第十阶，这点我还没弄明白。
其二，在加载求解中，将文件数据读入数组及加载方式不同，结果计算时间相差很大，这点我也还没完全明白。
最后，关于考虑重力的作用中，在进行重力静力分析中打开了预应力开关，但对最后结果影响不大，不知道这是什么原因。
反应谱分析中，也有两个很有痛的问题：其一，关于两个方向反应谱的加载求解过程，不明白是模态合并后进行下一个，还是直接先做两次谱解最后共同模态扩展、模态合并；其二，就是结果处理的问题，两次谱分析要用到荷载工况，来组合最后的结果，但荷载工况的写入过程很伤脑经，不知道是每次谱分析都写还最后一起写，我考虑应该分开写，但另一个问题是先前的谱分析结果数据对后面的谱分析数据会不会有干扰，从而到时写入的工况是否不正确，这些都在不断的尝试中去分析。
最后选择了谱分析命令流附件中的那种求解过程和工况处理方法。
5.3 体会及建议自己动手做这个大作业，虽然历时很长，而且过程中也遇到了特别多的问题，但是在学习时程分析和谱分析的方法上，以及分析问题的能力上都有很大收获。
虽然之前，在有限元方法的课程上学过基本的ANSYS操作，平时自己也做一些相关的小题，不过要说在ANSYS工程实例上真正的探索性学习，这是第一次。
这门课程的最后大作业前前后后我差不多做了一个月，遇到了很多问题，也不断地解决了问题，最后独立的完成了这份报告，虽然花费的时间很多但感觉收获丰满，我相信这为我接下来的科研之路铺下了坚实的基础。
除报告中体现的一些过程和结论外，自己也做了很多尝试获得很多探索性的结论，这过程中很感谢老师的帮助以及CAE班QQ群上同学的共同探讨及互相分享心得，也特别感激教研室一起学习、一起生活的朋友。
关于这门课，虽然收获挺大，但仍觉得自己只是学到了老师所授的墙角一隅，还有很多方法和操作没来得及及时操作，老师特别认真，每堂课都准备满满，上课的状态也非常好，上课时，我积极的做笔记，竖起耳朵认真听着老师讲的每一句话，最终下来发现自己确实知道了很多ANSYS及其他软件的功能和相关操作，但自己真正能吸收，化为自己的知识却非常至少，可能的原因是自己下课后没有及时复习上课的内容，另外对原理上的东西领悟不透，太多的理论没时间总结，太多的操作没及时付诸实施，最后也就很快就忘掉了。