概述:以Koyna混凝土坝为对象进行地震响应计算。将自编的八节点UEL和二十节点UEL应用到计算中。分别进行了混凝土坝模态计算和地震时程计算。
其中,在模态计算中共设置四种计算工况,分别为:ABAQUS-C3D8、UEL-C3D8、ABAQUS-C3D20、UEL-C3D20。
在地震时程计算中设置两种计算工况,分别为:ABAQUS-C3D8、UEL-C3D8。
计算结果表明,自编UEL与ABAQUS自带单元结果一致。
()模型信息
Koyna混凝土重力坝位于印度孟买东南200 km处,1967年12月11日,Koyna混凝土重力坝遭受里氏6.5级的地震(Koyna地震),该地震给大坝和水电站等水工建筑带来了巨大的损坏,给下游数十万居民的人身生命财产安全造成了巨大损失,该大坝地震案例事后成为诸多学者进行地震作用下坝基动力相互作用、混凝土材料的动态力学性能等领域研究的对象。
Koyna混凝土重力坝的横断面尺寸如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图1](static/picture/f54cc0b2fcc34af59f7df8080148314e-formatwebp.jpg)
该坝坝高113 m,坝顶宽度14.8 m,坝底宽度70 m,正常蓄水位91.75 m,坝基尺寸分别取向上游、下游延伸两倍坝高,向地基方向同样延伸两倍坝高,向横河向延伸20 m,最终的三维几何模型示意图如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图2](static/picture/f49b971bf5b74bd9a4c8318c8f6c3637-formatwebp.jpg)
采用六面体单元离散,有限元计算模型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图3](static/picture/08079aec62d2470b8ff2754961a3a04f-formatwebpresizew_760.jpg)
该模型的坝体和坝基共计用17950个六面体单元,其中坝体5000个,坝基12950个,坝基的网格采用疏密渐进过渡的方式避免计算结果在网格突变处不连续。
模型的材料属性如下:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图4](static/picture/aa129a5717b642c292ba523a011aacc6-formatwebpresizew_760.jpg)
注:这里地基的密度为零是因为没有考虑土-结构动力相互作用。
()模态分析
四种工况的前100阶模态具体数值为:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图5](static/picture/ebdaa7daf9cf495eae777d915548bbb9-formatwebp.jpg)
数值对比为:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图6](static/picture/15fb67e9adae4b41b20614b3b45c84fc.png)
相对误差为:(注:这里的标签应为:ABAQUS-C3D8——UEL-C3D8 ABAQUS-C3D20——UEL-C3D20)
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图7](static/picture/06b659123c2f463bbd2ff0f349296da7-formatwebpresizew_760.png)
自编UEL与ABAQUS自带单元比较,相对误差保持在0.1%以下。
若干振型图如下:
第1阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图8](static/picture/ff75cf77cd6e4212a4529a4eee2cbcab-formatwebpresizew_760.jpg)
第2阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图9](static/picture/9cb9c28c0caa4fdf89d78cb64ac38f9e-formatwebpresizew_760.jpg)
第3阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图10](static/picture/d2dbe3eb08454fcf95729f5a56b7c38e-formatwebpresizew_760.jpg)
第10阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图11](static/picture/9897d403e7a849b6aabce65c89141519-formatwebpresizew_760.jpg)
第50阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图12](static/picture/49f2f243fe2c44d18af035b762695cf0-formatwebpresizew_760.jpg)
第100阶振型如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图13](static/picture/935d3d0de67840b99de6aeb6020391fd-formatwebpresizew_760.jpg)
()地震时程分析
坝顶坝锺顺河向相对位移如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图14](static/picture/aae459bafb3341cda0876c5643738fc4-formatwebpresizew_760.jpg)
坝顶坝锺竖直向相对位移如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图15](static/picture/6584857bb3944c3b9a24bde0f1e7f623-formatwebpresizew_760.jpg)
坝顶坝锺横河向相对位移如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图16](static/picture/00632de3bf7445e6969f3afa45fa2cda-formatwebpresizew_760.jpg)
坝顶坝锺顺河向相对速度如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图17](static/picture/735139249a4447b8b873d906402d484e-formatwebpresizew_760.jpg)
坝顶坝锺竖直向相对速度如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图18](static/picture/f465556e1e024ddca120a68cd5f99df0-formatwebpresizew_760.jpg)
坝顶坝锺横河向相对速度如下图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图19](static/picture/036d51c69bdd44b59730525845f2c34c-formatwebpresizew_760.jpg)
顺河向速度变化云图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图20](static/picture/youku-case-formatwebpresizew_760.png)
加速度变化云图:
![ABAQUS UEL 二次开发(Koyna混凝土坝地震响应计算)的图21](static/picture/youku-case-formatwebpresizew_760.png)
()附件
----------------------- ----------------------- KOYNA HISTORY ABAQUS-C3D8 JOB-1.INP KOYNA_HACCEL.INP KOYNA_VACCEL.INP RUN.BAT UEL-C3D8 BBAR.OBJ JOB-1.INP KOYNA_HACCEL.INP KOYNA_VACCEL.INP RUNBBAR.BAT DATA.OPJU MODEL ABAQUS-C3C8 JOB-1.INP RUN.BAT ABAQUS-C3C20 JOB-1.INP RUN.BAT UEL-C3C8 BBAR.OBJ JOB-1.INP RUNBBAR.BAT UEL-C3C20 A.OBJ JOB-1.INP RUNBBAR.BAT DATA ----------------------- -----------------------