1结构抗震理论结构抗震理论的发展经历了由简单到复杂、由静力分析到动力分析的过程。早期的结构抗震设计主要采用静力分析法其核心思想是将地震作用等效为静力荷载进行计算因此形成了基于承载力的抗震设计方法。随着对地震作用和结构响应认识的深入反应谱分析法逐渐发展起来它能够考虑结构的动力特性并在一定程度上反映地震作用下结构的最大响应因此推动了基于承载力和延性的抗震设计方法的发展。进一步地动力分析法成为更精细的抗震分析手段它可以直接研究结构在地震波作用下随时间变化的响应过程并结合时程分析法和推覆分析法等具体方法对结构的弹塑性变形、损伤演化和抗震性能进行更全面的评价。由此抗震设计理念也从最初单纯基于承载力的设计逐步发展到兼顾承载力与延性的设计并最终向基于性能的抗震设计方向演进。1.1 静力分析法静力分析法也称为“底部剪力法”假设结构物与地震动具有相同的振动把结构物在地面运动加速度a作用下产生的惯性力视作静力作用于结构物上的抗震计算。惯性力计算公式为a为地震动最大水平加速度g为重力加速度Gmg为建筑物的重量m为建筑物的质量k为地面运动加速度峰值与重力加速度的比值称为地震系数其数值与结构动力特性无关。静力分析法忽略了结构的动力学特性这一重要因素把地震动加速度看作是结构地震破坏的单一因素因此有很大的局限性常导致对结构抗震能力的错误判断仅当结构物在振动时几乎不产生变形而被当做刚体的时候静力分析法才比较适用因此现在已较少使用。1.2 反应谱分析法反应谱分析法通过反应谱来计算结构动力特性自振周期、振型和阻尼所对应的共振效应但其计算公式仍保留了早期静力法理论的形式。地震时结构所受到的最大水平基底剪力即总水平地震作用为其中k为地震系数β(T)则是加速度反应谱与地震动最大加速度的比值它表示地震时结构振动加速度的放大倍数。反应谱理论尽管考虑了结构的动力特性然而在结构设计中它仍然把地震惯性力作为静力来对待所以它只能称为准动力理论。1.3 时程分析法时程分析法是利用实际地震地面运动加速度记录进行抗震分析与设计能考虑结构的弹性和弹塑性性态能反映地震动的三大要素频谱、振幅和持时因而在复杂工程结构抗震分析和设计中得到了广泛的应用。它具有如下特点1输入地震动参数要求采用符合场地情况具有概率意义的加速度时程对于复杂结构要求考虑地震动三个分量的时间过程及其空间相关性2要求采用的结构和构件的动力模型应接近实际情况应考虑结构和构件的非线性恢复力特性3该方法能计算出结构反应的全过程包括变形和能量损耗的积累。时程分析法是先进的但是在应用上还有很多困难和局限尤其是工程设计时的应用例如输入地震波的不确定性、结构性能的近似假定和模拟等使分析结果的可信度受到限制该方法需要专门的程序和应用知识输入输出数据量大计算技术复杂一般需要专门的技术人员进行分析。因此我国规范只要求少数重要超高或有薄弱位置的结构采用时程分析方法进行多发地震下的补充计算或罕见地震作用下薄弱层的弹塑性变形验算。2 模态分析理论模态分析理论有三个基本假设时不变性假设系统的动态特性不会随着时间的变化而相应的变化线性假设任意的激励引起的结构的响应等于各个激励分别作用于结构所引起响应的线性叠加可观测性假设全部的数据信息都可以观测得到即确定系统的动态特性所需要的所有数据为已知结构动力学基本方程通常结构自身的动力学特性与外载荷激励无关令{F(t)}0得到而阻尼对于结构自身的动力学特性影响不大工程上往往忽略阻尼项由此得到在有限元分析中常用的阻尼模型主要包括粘性阻尼、材料阻尼和摩擦阻尼三类。粘性阻尼是有限元法中最常使用的阻尼模型特点是阻尼力的大小与结构运动速度成正比。在工程分析中常采用线性阻尼模型表示粘性阻尼典型形式为瑞利阻尼即材料阻尼主要由材料内部能量耗散产生取决于材料的本构关系。在这种模型中阻尼应力通常与应变率成正比。例如Voigt模型可以表示为材料阻尼中粘弹性阻尼模型目前在工程中被广泛使用。摩擦阻尼则来源于结构或构件接触面之间的相对滑移和摩擦耗能。最简单的是库仑摩擦阻尼模型该模型通常属于非线性阻尼模型能够反映摩擦力与相对运动方向之间的关系。3 大坝未蓄水时的模态分析整体上可以分为以下步骤1读取网格定义单元LIRE_MAILIAGE读取网格格式medAFFE_MODELE定义有限元模型的单元类型定义拱坝为3D模型2定义材料赋予材料DEFI_MATERIAU定义拱坝相关的材料参数包括瑞利阻尼系数为后续抗震分析做准备AFFE_MATERIAU将定义好的材料属性赋予对应的有限元模型3边界条件AFFE_CHAR_MECA定义拱坝与地基交界面处的边界条件LIAISON_UNIF可以连接单个节点和一个节点组从而通过单个节点控制节点组的自由度load AFFE_CHAR_MECA(DDL_IMPO_F(DX0.0, DY0.0, DZ0.0, GROUP_NO(arch_bottom_1p, )), LIAISON_UNIF_F(DDL(DX, DY, DZ), GROUP_NO(arch_bottom, )), MODELEmodel)4计算质量、刚度矩阵ASSEMBLAGE计算拱坝结构的质量矩阵、刚度矩阵和阻尼矩阵并定义一个数据结构NUME_DDL这里的阻尼矩阵定义是为了后续的地震分析做准备ASSEMBLAGE(CHAM_MATERfieldmat, CHARGE(load, ), MATR_ASSE(_F(MATRICECO(M), OPTIONMASS_MECA), _F(MATRICECO(K), OPTIONRIGI_MECA), _F(MATRICECO(C), OPTIONAMOR_MECA)), MODELEmodel, NUME_DDLCO(NDDL))5模态计算CALC_MODES计算20阶模态6结果后处理及保存IMPR_RESU输出计算结果查看振型4 大坝蓄水后的模态分析蓄水后的分析步骤包括1矩阵投影PROJ_MATR_BASE / PROJ_BASE将质量矩阵、刚度矩阵和阻尼矩阵投影到前面计算的未蓄水时的模态空间中两个命令中前者一次只能完成一个投影操作后者可以同时进行多个矩阵投影操作PROJ_BASE(BASEmodes, MATR_ASSE_GENE(_F(MATRICECO(M1), MATR_ASSEM), _F(MATRICECO(K1), MATR_ASSEK), _F(MATRICECO(C1), MATR_ASSEC)), NUME_DDL_GENECO(NDDL1))NUME_DDL_GENE基于已计算的模态空间生成新的方程编号即新的模态特征值对应的方程从1~20进行编号不再与节点和物理自由度有关只与前面计算的特征模态有关。2库水附加质量MACRO_MATR_AJOU计算未蓄水大坝的库水附加质量并给库水赋予材料属性注意需要指定自由水面的压力值0MACRO_MATR_AJOU(DDL_IMPO_F(GROUP_NO(water_surface, ), PRES_FLUIDE0.0), FLUIDE_F(RHO1000.0, TOUTOUI), GROUP_MA_FLUIDE(water, ), GROUP_MA_INTERF(arch_face, ), MAILLAGEmesh, MATR_MASS_AJOUCO(WATER), MODELISATION3D, MODE_MECAmodes, NUME_DDL_GENENDDL1)3质量矩阵合并COMB_MATR_ASSE合并大坝结构质量矩阵和库水的附加质量矩阵M1 COMB_MATR_ASSE(reuseM1, COMB_R(_F(COEF_R1.0, MATR_ASSEWATER), _F(COEF_R1.0, MATR_ASSEM1)))4模态计算CALC_MODES运用原来的刚度矩阵和合并后的质量矩阵计算新的模态5转换模态空间基REST_GENE_PHYS用附加质量的方法计算的模态结果是基于合并后的模态空间坐标系该命令将计算结果恢复为广义的物理空间坐标系以便输出模态振型结果modes1 CALC_MODES(CALC_FREQ_F(NMAX_FREQ20), MATR_MASSM1, MATR_RIGIK1, OPTIONPLUS_PETITE)restran REST_GENE_PHYS(RESU_GENEmodes1)6结果后处理及保存IMPR_RESU输出模态振型结果5 大坝抗震分析1静模态计算MODE_STATIQUE计算坝基位移约束相关的静态模态modes0 MODE_STATIQUE(MATR_MASSM, MATR_RIGIK, MODE_STAT_F(AVEC_CMP(DX, DY, DZ), GROUP_NO(arch_bottom_1p, )))2计算地震谱CALC_CHAR_SEISME基于上面的静模态计算模型生成X/Y/Z三个方向地震载荷向量LIRE_FONCTION以函数的方式读入地震载荷谱数据X CALC_CHAR_SEISME(DIRECTION(1.0, 0.0, 0.0), GROUP_NO(arch_bottom_1p, ), MATR_MASSM, MODE_STATmodes0)ax LIRE_FONCTION(INDIC_PARA(1, 1), INDIC_RESU(1, 2), INTERPOL(LIN, ), NOM_PARAINST, PROL_DROITECONSTANT, PROL_GAUCHECONSTANT, TYPEFONCTION, UNITE2)3矩阵投影PROJ_BASE将质量矩阵、刚度矩阵和阻尼矩阵投影到蓄水模态计算空间中包括三个方向的地震载荷向量投影PROJ_BASE(BASErestran, MATR_ASSE_GENE(_F(MATRICECO(M3), MATR_ASSEM), _F(MATRICECO(K3), MATR_ASSEK), _F(MATRICECO(C3), MATR_ASSEC)), NUME_DDL_GENECO(DDL3), VECT_ASSE_GENE(_F(TYPE_VECTFORC, VECTEURCO(X1), VECT_ASSEX), _F(TYPE_VECTFORC, VECTEURCO(Y1), VECT_ASSEY), _F(TYPE_VECTFORC, VECTEURCO(Z1), VECT_ASSEZ)))4反应谱分析DYNA_VIBRA计算基于蓄水大坝模态空间的瞬态结构响应1.选择合适的时间步2.对比积分方案欧拉积分与龙格库塔积分3.施加地震载荷resharm DYNA_VIBRA(BASE_CALCULGENE, EXCIT(_F(FONC_MULTax, VECT_ASSE_GENEX1), _F(FONC_MULTay, VECT_ASSE_GENEY1), _F(FONC_MULTaz, VECT_ASSE_GENEZ1)), INCREMENT_F(INST_FIN4.0, INST_INIT0.0, PAS0.002), MATR_AMORC3, MATR_MASSM3, MATR_RIGIK3, TYPE_CALCULTRAN)5转换模态空间基REST_GENE_PHYS用附加质量的方法计算的模态结果是基于合并后的模态空间坐标系该命令将计算结果恢复为广义的物理空间坐标系以便输出模态振型结果6结果后处理及保存RECU_FONCTION后处理提取三个位移分量DX/DY/DZ生成时间历程曲线IMPR_RESU输出模态振型结果DX RECU_FONCTION(GROUP_NO(arch_edge_1p, ), NOM_CHAMDEPL, NOM_CMPDX, RESULTATrestran0)
开源CAE实战系列（十一）：Code_Aster应用实例之混凝土大坝的结构抗震分析
发布时间：2026/6/13 5:03:05
1结构抗震理论结构抗震理论的发展经历了由简单到复杂、由静力分析到动力分析的过程。早期的结构抗震设计主要采用静力分析法其核心思想是将地震作用等效为静力荷载进行计算因此形成了基于承载力的抗震设计方法。随着对地震作用和结构响应认识的深入反应谱分析法逐渐发展起来它能够考虑结构的动力特性并在一定程度上反映地震作用下结构的最大响应因此推动了基于承载力和延性的抗震设计方法的发展。进一步地动力分析法成为更精细的抗震分析手段它可以直接研究结构在地震波作用下随时间变化的响应过程并结合时程分析法和推覆分析法等具体方法对结构的弹塑性变形、损伤演化和抗震性能进行更全面的评价。由此抗震设计理念也从最初单纯基于承载力的设计逐步发展到兼顾承载力与延性的设计并最终向基于性能的抗震设计方向演进。1.1 静力分析法静力分析法也称为“底部剪力法”假设结构物与地震动具有相同的振动把结构物在地面运动加速度a作用下产生的惯性力视作静力作用于结构物上的抗震计算。惯性力计算公式为a为地震动最大水平加速度g为重力加速度Gmg为建筑物的重量m为建筑物的质量k为地面运动加速度峰值与重力加速度的比值称为地震系数其数值与结构动力特性无关。静力分析法忽略了结构的动力学特性这一重要因素把地震动加速度看作是结构地震破坏的单一因素因此有很大的局限性常导致对结构抗震能力的错误判断仅当结构物在振动时几乎不产生变形而被当做刚体的时候静力分析法才比较适用因此现在已较少使用。1.2 反应谱分析法反应谱分析法通过反应谱来计算结构动力特性自振周期、振型和阻尼所对应的共振效应但其计算公式仍保留了早期静力法理论的形式。地震时结构所受到的最大水平基底剪力即总水平地震作用为其中k为地震系数β(T)则是加速度反应谱与地震动最大加速度的比值它表示地震时结构振动加速度的放大倍数。反应谱理论尽管考虑了结构的动力特性然而在结构设计中它仍然把地震惯性力作为静力来对待所以它只能称为准动力理论。1.3 时程分析法时程分析法是利用实际地震地面运动加速度记录进行抗震分析与设计能考虑结构的弹性和弹塑性性态能反映地震动的三大要素频谱、振幅和持时因而在复杂工程结构抗震分析和设计中得到了广泛的应用。它具有如下特点1输入地震动参数要求采用符合场地情况具有概率意义的加速度时程对于复杂结构要求考虑地震动三个分量的时间过程及其空间相关性2要求采用的结构和构件的动力模型应接近实际情况应考虑结构和构件的非线性恢复力特性3该方法能计算出结构反应的全过程包括变形和能量损耗的积累。时程分析法是先进的但是在应用上还有很多困难和局限尤其是工程设计时的应用例如输入地震波的不确定性、结构性能的近似假定和模拟等使分析结果的可信度受到限制该方法需要专门的程序和应用知识输入输出数据量大计算技术复杂一般需要专门的技术人员进行分析。因此我国规范只要求少数重要超高或有薄弱位置的结构采用时程分析方法进行多发地震下的补充计算或罕见地震作用下薄弱层的弹塑性变形验算。2 模态分析理论模态分析理论有三个基本假设时不变性假设系统的动态特性不会随着时间的变化而相应的变化线性假设任意的激励引起的结构的响应等于各个激励分别作用于结构所引起响应的线性叠加可观测性假设全部的数据信息都可以观测得到即确定系统的动态特性所需要的所有数据为已知结构动力学基本方程通常结构自身的动力学特性与外载荷激励无关令{F(t)}0得到而阻尼对于结构自身的动力学特性影响不大工程上往往忽略阻尼项由此得到在有限元分析中常用的阻尼模型主要包括粘性阻尼、材料阻尼和摩擦阻尼三类。粘性阻尼是有限元法中最常使用的阻尼模型特点是阻尼力的大小与结构运动速度成正比。在工程分析中常采用线性阻尼模型表示粘性阻尼典型形式为瑞利阻尼即材料阻尼主要由材料内部能量耗散产生取决于材料的本构关系。在这种模型中阻尼应力通常与应变率成正比。例如Voigt模型可以表示为材料阻尼中粘弹性阻尼模型目前在工程中被广泛使用。摩擦阻尼则来源于结构或构件接触面之间的相对滑移和摩擦耗能。最简单的是库仑摩擦阻尼模型该模型通常属于非线性阻尼模型能够反映摩擦力与相对运动方向之间的关系。3 大坝未蓄水时的模态分析整体上可以分为以下步骤1读取网格定义单元LIRE_MAILIAGE读取网格格式medAFFE_MODELE定义有限元模型的单元类型定义拱坝为3D模型2定义材料赋予材料DEFI_MATERIAU定义拱坝相关的材料参数包括瑞利阻尼系数为后续抗震分析做准备AFFE_MATERIAU将定义好的材料属性赋予对应的有限元模型3边界条件AFFE_CHAR_MECA定义拱坝与地基交界面处的边界条件LIAISON_UNIF可以连接单个节点和一个节点组从而通过单个节点控制节点组的自由度load AFFE_CHAR_MECA(DDL_IMPO_F(DX0.0, DY0.0, DZ0.0, GROUP_NO(arch_bottom_1p, )), LIAISON_UNIF_F(DDL(DX, DY, DZ), GROUP_NO(arch_bottom, )), MODELEmodel)4计算质量、刚度矩阵ASSEMBLAGE计算拱坝结构的质量矩阵、刚度矩阵和阻尼矩阵并定义一个数据结构NUME_DDL这里的阻尼矩阵定义是为了后续的地震分析做准备ASSEMBLAGE(CHAM_MATERfieldmat, CHARGE(load, ), MATR_ASSE(_F(MATRICECO(M), OPTIONMASS_MECA), _F(MATRICECO(K), OPTIONRIGI_MECA), _F(MATRICECO(C), OPTIONAMOR_MECA)), MODELEmodel, NUME_DDLCO(NDDL))5模态计算CALC_MODES计算20阶模态6结果后处理及保存IMPR_RESU输出计算结果查看振型4 大坝蓄水后的模态分析蓄水后的分析步骤包括1矩阵投影PROJ_MATR_BASE / PROJ_BASE将质量矩阵、刚度矩阵和阻尼矩阵投影到前面计算的未蓄水时的模态空间中两个命令中前者一次只能完成一个投影操作后者可以同时进行多个矩阵投影操作PROJ_BASE(BASEmodes, MATR_ASSE_GENE(_F(MATRICECO(M1), MATR_ASSEM), _F(MATRICECO(K1), MATR_ASSEK), _F(MATRICECO(C1), MATR_ASSEC)), NUME_DDL_GENECO(NDDL1))NUME_DDL_GENE基于已计算的模态空间生成新的方程编号即新的模态特征值对应的方程从1~20进行编号不再与节点和物理自由度有关只与前面计算的特征模态有关。2库水附加质量MACRO_MATR_AJOU计算未蓄水大坝的库水附加质量并给库水赋予材料属性注意需要指定自由水面的压力值0MACRO_MATR_AJOU(DDL_IMPO_F(GROUP_NO(water_surface, ), PRES_FLUIDE0.0), FLUIDE_F(RHO1000.0, TOUTOUI), GROUP_MA_FLUIDE(water, ), GROUP_MA_INTERF(arch_face, ), MAILLAGEmesh, MATR_MASS_AJOUCO(WATER), MODELISATION3D, MODE_MECAmodes, NUME_DDL_GENENDDL1)3质量矩阵合并COMB_MATR_ASSE合并大坝结构质量矩阵和库水的附加质量矩阵M1 COMB_MATR_ASSE(reuseM1, COMB_R(_F(COEF_R1.0, MATR_ASSEWATER), _F(COEF_R1.0, MATR_ASSEM1)))4模态计算CALC_MODES运用原来的刚度矩阵和合并后的质量矩阵计算新的模态5转换模态空间基REST_GENE_PHYS用附加质量的方法计算的模态结果是基于合并后的模态空间坐标系该命令将计算结果恢复为广义的物理空间坐标系以便输出模态振型结果modes1 CALC_MODES(CALC_FREQ_F(NMAX_FREQ20), MATR_MASSM1, MATR_RIGIK1, OPTIONPLUS_PETITE)restran REST_GENE_PHYS(RESU_GENEmodes1)6结果后处理及保存IMPR_RESU输出模态振型结果5 大坝抗震分析1静模态计算MODE_STATIQUE计算坝基位移约束相关的静态模态modes0 MODE_STATIQUE(MATR_MASSM, MATR_RIGIK, MODE_STAT_F(AVEC_CMP(DX, DY, DZ), GROUP_NO(arch_bottom_1p, )))2计算地震谱CALC_CHAR_SEISME基于上面的静模态计算模型生成X/Y/Z三个方向地震载荷向量LIRE_FONCTION以函数的方式读入地震载荷谱数据X CALC_CHAR_SEISME(DIRECTION(1.0, 0.0, 0.0), GROUP_NO(arch_bottom_1p, ), MATR_MASSM, MODE_STATmodes0)ax LIRE_FONCTION(INDIC_PARA(1, 1), INDIC_RESU(1, 2), INTERPOL(LIN, ), NOM_PARAINST, PROL_DROITECONSTANT, PROL_GAUCHECONSTANT, TYPEFONCTION, UNITE2)3矩阵投影PROJ_BASE将质量矩阵、刚度矩阵和阻尼矩阵投影到蓄水模态计算空间中包括三个方向的地震载荷向量投影PROJ_BASE(BASErestran, MATR_ASSE_GENE(_F(MATRICECO(M3), MATR_ASSEM), _F(MATRICECO(K3), MATR_ASSEK), _F(MATRICECO(C3), MATR_ASSEC)), NUME_DDL_GENECO(DDL3), VECT_ASSE_GENE(_F(TYPE_VECTFORC, VECTEURCO(X1), VECT_ASSEX), _F(TYPE_VECTFORC, VECTEURCO(Y1), VECT_ASSEY), _F(TYPE_VECTFORC, VECTEURCO(Z1), VECT_ASSEZ)))4反应谱分析DYNA_VIBRA计算基于蓄水大坝模态空间的瞬态结构响应1.选择合适的时间步2.对比积分方案欧拉积分与龙格库塔积分3.施加地震载荷resharm DYNA_VIBRA(BASE_CALCULGENE, EXCIT(_F(FONC_MULTax, VECT_ASSE_GENEX1), _F(FONC_MULTay, VECT_ASSE_GENEY1), _F(FONC_MULTaz, VECT_ASSE_GENEZ1)), INCREMENT_F(INST_FIN4.0, INST_INIT0.0, PAS0.002), MATR_AMORC3, MATR_MASSM3, MATR_RIGIK3, TYPE_CALCULTRAN)5转换模态空间基REST_GENE_PHYS用附加质量的方法计算的模态结果是基于合并后的模态空间坐标系该命令将计算结果恢复为广义的物理空间坐标系以便输出模态振型结果6结果后处理及保存RECU_FONCTION后处理提取三个位移分量DX/DY/DZ生成时间历程曲线IMPR_RESU输出模态振型结果DX RECU_FONCTION(GROUP_NO(arch_edge_1p, ), NOM_CHAMDEPL, NOM_CMPDX, RESULTATrestran0)
的演进与通信机制)
