基于混凝土损伤塑性模型的桥墩地震响应分析

ABAQUS混凝土塑性损伤模型

4.5.2 混凝土和其它准脆性材料的塑性损伤模型 这部分介绍的是ABAQUS提供分析混凝土和其它准脆性材料的混凝土塑性损伤模型。ABAQUS 材料库中也包括分析混凝的其它模型如基于弥散裂纹方法的土本构模型。他们分别是在ABAQUS/Standard “An inelastic constitutive model for concrete,” Section 4.5.1, 中的弥散裂纹模型和在ABAQUS/Explicit, “A cracking model for concrete and other brittle materials,” Section 4.5.3中的脆性开裂模型。 混凝土塑性损伤模型主要是用来为分析混凝土结构在循环和动力荷载作用下的提供一个普遍分析模型。该模型也适用于其它准脆性材料如岩石、砂浆和陶瓷的分析；本节将以混凝土的力学行为来演示本模型的一些特点。在较低的围压下混凝土表现出脆性性质，主要的失效机制是拉力作用下的开裂失效和压力作用下的压碎。当围压足够大能够阻止裂纹开裂时脆性就不太明显了。这种情况下混凝土失效主要表现为微孔洞结构的聚集和坍塌，从而导致混凝土的宏观力学性质表现得像具有强化性质的延性材料那样。 本节介绍的塑性损伤模型并不能有效模拟混凝土在高围压作用下的力学行为。而只能模拟混凝土和其它脆性材料在与中等围压条件（围压通常小于单轴抗压强度的四分之一或五分之一）下不可逆损伤有关的一些特性。这些特性在宏观上表现如下： ?单拉和单压强度不同，单压强度是单拉强度的10倍甚至更多； ?受拉软化，而受压在软化前存在强化； ?在循环荷载(压)下存在刚度恢复； ?率敏感性，尤其是强度随应变率增加而有较大的提高。 概论 混凝土非粘性塑性损伤模型的基本要点介绍如下： 应变率分解 对率无关的模型附加假定应变率是可以如下分解的： 是总应变率，是应变率的弹性部分，是应变率的塑性部分。 应力应变关系 应力应变关系为下列弹性标量损伤关系： 其中是材料的初始（无损）刚度，是有损刚度，是刚度退化变量其值在0（无损）到1（完全失效）之间变化，与失效机制（开裂和压碎）相关的损伤导致了弹性刚度的退化。在标量损伤理论框架内，刚度退化是各向同性的，它可由单个标量d来描述。按照传统连续介质力学观点，有效应力可定义如下：

混凝土塑性损伤模型1

混凝土和其它准脆性材料的塑性损伤模型 这部分介绍的是ABAQUS提供分析混凝土和其它准脆性材料的混凝土塑性损伤模型。ABAQUS 材料库中也包括分析混凝的其它模型如基于弥散裂纹方法的土本构模型。他们分别是在ABAQUS/Standard “An inelastic constitutive model for concrete,” Section 4.5.1, 中的弥散裂纹模型和在ABAQUS/Explicit, “A cracking model for concrete and other brittle materials,” Section 4.5.3中的脆性开裂模型。 混凝土塑性损伤模型主要是用来为分析混凝土结构在循环和动力荷载作用下的提供一个普遍分析模型。该模型也适用于其它准脆性材料如岩石、砂浆和陶瓷的分析；本节将以混凝土的力学行为来演示本模型的一些特点。在较低的围压下混凝土表现出脆性性质，主要的失效机制是拉力作用下的开裂失效和压力作用下的压碎。当围压足够大能够阻止裂纹开裂时脆性就不太明显了。这种情况下混凝土失效主要表现为微孔洞结构的聚集和坍塌，从而导致混凝土的宏观力学性质表现得像具有强化性质的延性材料那样。 本节介绍的塑性损伤模型并不能有效模拟混凝土在高围压作用下的力学行为。而只能模拟混凝土和其它脆性材料在与中等围压条件（围压通常小于单轴抗压强度的四分之一或五分之一）下不可逆损伤有关的一些特性。这些特性在宏观上表现如下： ?单拉和单压强度不同，单压强度是单拉强度的10倍甚至更多； ?受拉软化，而受压在软化前存在强化； ?在循环荷载(压)下存在刚度恢复； ?率敏感性，尤其是强度随应变率增加而有较大的提高。 概论 混凝土非粘性塑性损伤模型的基本要点介绍如下： 应变率分解 对率无关的模型附加假定应变率是可以如下分解的： 是总应变率，是应变率的弹性部分，是应变率的塑性部分。 应力应变关系 应力应变关系为下列弹性标量损伤关系： 其中是材料的初始（无损）刚度，是有损刚度，是刚度退化变量其值在0（无损）到1（完全失效）之间变化，与失效机制（开裂和压碎）相关的损伤导致了弹性刚度的退化。在标量损伤理论框架内，刚度退化是各向同性的，它可由单个标量d来描述。按照传统连续介质力学观点，有效应力可定义如下：

地震反应谱分析实例

结构地震反应谱分析实例 在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0

!进行模态求解 ANTYPE,MODAL MODOPT,LANB,30 SOLVE FINISH !进行谱分析 /SOLU ANTYPE,SPECTR SPOPT,SPRS,30,YES SVTYP,2 !加速度反应谱 SED,1,1 !X与Y向 FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167 FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000 SV,0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716 SV,0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.4853,0.4853,0.4853,0.4853 SV,0.05,0.2588,0.2167 SOLVE FINISH !进行模态求解(模态扩展) /SOLU ANTYPE,MODAL EXPASS,ON MXPAND,30,,,YES,0.005 SOLVE FINISH !进行谱分析(合并模态) /SOLU ANTYPE,SPECTR SRSS,0.15,disp SOLVE FINISH /POST1 SET,LIST !结果1 /INP,,mcom

ABAQUS_混凝土损伤塑性模型_损伤因子

混凝土损伤因子的定义 BY lizhenxian27 1 损伤因子的定义 损伤理论最早是1958年Kachanov提出来用于研究金属徐变的。所谓损伤，是指在各种加载条件下，材料内凝聚力的进展性减弱，并导致体积单元破坏的现象，是受载材料由于微缺陷(微裂纹和微孔洞)的产生和发展而引起的逐步劣化。损伤一般被作为一种“劣化因素”而结合到弹性、塑性和粘塑性介质中去。 由于损伤的发展和材料结构的某种不可逆变化，因而不同的学者采用了不同的损伤定义。一般来说，按使用的基准可将损伤分为: (1) 微观基准量 1,空隙的数目、长度、面积、体积; 2空隙的形状、排列、由取向所决定的有效面积。 (2) 宏观基准量 1、弹性常数、屈服应力、拉伸强度、延伸率。 2、密度、电阻、超声波波速、声发射。 对于第一类基准量，不能直接与宏观力学量建立物性关系，所以用它来定义损伤变量的时候，需要对它做出一定的宏观尺度下的统计处理(如平均、求和等)。 对于第二类基准量，一般总是采用那些对损伤过程比较敏感，在实验室里易于测量的量，作为损伤变量的依据。 由于微裂纹和微孔洞的存在，微缺陷所导致的微应力集中以及缺陷的相互作用，有效承

载面积由

A 减小为A ’。如假定这些微裂纹和微孔洞在空间各个方向均匀分布，A ’与法向无关，这时可定义各向同性损伤变量D 为 D= ( A- A ’ )/ A 事实上，微缺陷的取向、分布及演化与受载方向密切相关，因此材料损伤实际上是各向异性的。为描述损伤的各向异性，可采用张量形式来定义。损伤表征了材损伤是一个非负的因子，同时由于这一力学性能的不可逆性，必然有 0dD dt ≥ 2有效应力 定义Cauchy 有效应力张量'σ ''//(1)A A D σσσ==- 一般情况下，存在于物体内的损伤(微裂纹、空洞)是有方向性的。当损伤变量与受力面法向相关时，是为各向异性损伤;当损伤变量与法向无关时，为各向异性损伤。这时的损伤变量是一标量。 3等效性假设 损伤演化方程推导一般使用两种等效性假设，一种是应变等效性假设，另一种是能量等效性假设。采用能量等效性假设可以避免采用应变等效假设而使得各向异性损伤模型中的有效弹性矩阵不对称的问题.以下对两种假设进行简要的介绍。 (1) 应变等效性假设 1971年 Lematire 提出，损伤单元在应力σ作用下的应变响应与无损单元在定义的有效应力'σ作用下的应变响应相同。在外力作用下受损材料的本构关系可采用无损时的形式，只要

(仅供参考)Abaqus混凝土损伤塑性模型的参数标定

Abaqus 混凝土损伤塑性模型的参数标定 1. 塑性参数（Plasticity ） 1) 剪胀角（Dilation Angle ） = 30° 2) 流动势偏移量（Eccentricity ） 3) 双轴受压与单轴受压极限强度比 = 1.16 4) 不变量应力比 = 0.667 5) 粘滞系数（Visosity Parameter ） = 0.0005 2. 受压本构关系 应力-Yield Stress ：第一行应输入本构模型刚进入非弹性段非弹性应变为0时所对应的应力。 非弹性应变-Inelastic Strain （受拉时为开裂应变-Cracking Strain ）：根据应力按混凝土本构模型得出对应的应变值，并通过 , 和 ，得出非弹性应变。 3. 受压损伤因子（Damage Parameter ）计算 根据《Abaqus Analysis User's Manual (6.10)》 - 20.6.3 “Concrete damaged plasticity ”中公式： 假设非弹性应变 in c ε中塑性应变 pl c ε所占的比例为c β，通过转换可得损伤因子c d 的计算公式： () () 0 011in c c in c c c c E E d βεσβε-=+- 根据《ABAQUS 混凝土损伤塑性模型参数验证》规定，混凝土受压时c β的取值范围为0.35 ~ 0.7。

4. 受拉损伤因子（Damage Parameter ）计算 受拉损伤因子的计算与受压损伤因子的计算方法基本相同，只需将对应受压变量更换为受拉即可： () () 0011in t t in t t t t E E d βεσβε-=+- 而根据参考文献混凝土受拉时t β的取值范围为0.5 ~ 0.95。 5. 损伤恢复因子 受拉损伤恢复因子（Tension Recovery ）：缺省值0t w =。 受压损伤恢复因子（Compression Recovery ）：缺省值1c w =。

结构地震反应谱分析实例

在多位朋友的大力帮助下，经过半个多月的努力，鄙人终于对结构地震反应谱分析有了一定的了解，现将其求解步骤整理出来，以便各位参阅，同时，尚有一些问题，欢迎各位讨论！ 为叙述方便，举一简单实例： 在侧水压与顶部集中力作用下的柱子的地震反应谱分析，谱值为加速度反应谱，考虑X 与Y向地震效应作用。已知地震影响系数a与周期T的关系： a(T)= 0.4853*(0.4444+2.2222*T) 0<T<=0.04 秒 0.4853*(0.10/T)^(-0.686) 0.04<T<=0.1 秒 0.4853 0.1<T<=1.2 秒 0.4853*(1.2/T)^1.5 1.2<T<=4 秒 以下是命令流程序 ---------------------------------------------------------------------------------------------------- /filname,SPEC,1 /PREP7 !定义单元类型及材料特性 ET,1,45 MP,EX,1,2.8E10 MP,DENS,1,2.4E3 MP,NUXY,1,0.18 !建立模型 BLOCK,0,1,0,1,0,5 !网格剖分 ESIZE,0.5 VMESH,all /VIEW,,-0.3,-1,1 EPLOT FINISH /SOLU !施加底部约束 ASEL,,LOC,Z,0 DA,ALL,ALL ALLSEL !施加自重荷载 ACEL,0,0,10 !进行模态求解

ANTYPE,MODAL MODOPT,LANB,30 SOLVE FINISH !进行谱分析 /SOLU ANTYPE,SPECTR SPOPT,SPRS,30,YES SVTYP,2 !加速度反应谱 SED,1,1 !X与Y向 FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167 FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000 SV,0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716 SV,0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.4853,0.4853,0.4853,0.4853 SV,0.05,0.2588,0.2167 SOLVE FINISH !进行模态求解(模态扩展) /SOLU ANTYPE,MODAL EXPASS,ON MXPAND,30,,,YES,0.005 SOLVE FINISH !进行谱分析(合并模态) /SOLU ANTYPE,SPECTR SRSS,0.15,disp SOLVE FINISH /POST1 SET,LIST !结果1 /INP,,mcom lcwrite,11

混凝土损伤模型

塑性应力－应变关系 Compress Stress Inelastic strain Damage Inelastic strain 24019000000 292080000.00040.12990.0004 317090000.00080.24290.0008 323580000.00120.34120.0012 317680000.00160.42670.0016 303790000.0020.50120.002 285070000.00240.5660.0024 219070000.00360.7140.0036 148970000.0050.82430.005 29530000.010.96910.01 Tension recovery factor = 0.0 Tension Stress Cracking strain Damage Cracking strain 1.78E+06000 1.46E+060.00010.30.0001 1.11E+060.00030.550.0003 9.60E+050.00040.70.0004

8.00E+050.00050.80.0005 5.36E+050.00080.90.0008 3.59E+050.0010.930.001 1.61E+050.0020.950.002 7.30E+040.0030.970.003 4.00E+040.0050.990.005 Compression recovery factor = 1.0 **************************************************************************************** ********************************************************* 全应力-应变关系 Compress Stress Strain Tension Stress Strain 0000 0.000126480000.00001264800 0.000252960000.00002529600 0.000379440000.00003794400 0.0004105920000.000041059200 0.0005132400000.000051324000 0.0006158880000.000061588800 0.0007185360000.0000671780000 0.0008211840000.0001550231457000 0.000907062240190000.0003420321113000 0.001503021292080000.000436254960000 0.00199747317090000.000530211800000 0.002421979323580000.000820242536000 0.002799698317680000.001013557359000 0.003147243303790000.00200608161000 0.003476548285070000.00300275773000 0.004427304219070000.00500151140000 0.00556257614897000 0.0101115182953000

混凝土塑性损伤模型表格解读by自习菌

ABAQUS塑性损伤模型计算表格解读by自习菌（wx公众号） 受压本构： fc,r：砼单轴抗压强度标准值，可根据需要取多种值，此处取fck轴心抗压强度标准值 fck：C30,20.1MPa;C35,23.4;C40,26.8;C45,29.6;C50,32.4;C55,35.5;C60,38.5 εcr：与单轴抗压强度fc,r相应的峰值压应变，规范附录公式 αc：单轴受压应力-应变曲线下降段的形状参数，规范附录公式 εcu：应力-应变曲线下降段0.5 fc,r对应的压应变 εcu/εcr：规范附录公式 可适当修正抗压强度代表值fcr，峰值压应变εcr，以及曲线形状参数αc，砼规C.2.4附录。 Ec：弹性模量，只是辅助计算的一个临时取值。C30,3e4MPa;C40,3.25e4;C50,3.45e4 ρc：规范公式 n：规范公式 x：穷举数列，按规范公式与ε、εcr相关 dc：单轴受压损伤演化参数，以x=1为界限，规范为分段公式 ε：由x计算出，规范公式 σ：规范公式 σ修正：在应力-应变曲线上选定弹性阶段与弹塑性阶段的分界点，按Susoo88取0.4 fc,r， 或取1/3~1/2 fc,r，可见这也是一个可调整的值。通过这个选定的点的应力应变，计算弹性阶段的斜率，即E0弹性模量，这个弹性模量就是所采用本构的弹性模量，用E0和ε再重新反算弹性阶段的σ，即得“σ修正”。 对于C30砼，fc,r 取fck=20.1MPa，0.4*20.1=8.04MPa，在表格中插入一行，定义一个ε值，使σ无限逼近8.04（此时尚需重新定义表格这一行x列公式，使之由ε列导出）。根据这个应力应变值，求出E0，再由E0修正弹性阶段的应力值（即插入行之上的部分）。 【Susoo88：受压曲线与受拉曲线弹性临界点不一样，会产生两个弹性模量，需要在输入时选较大值，不然在后面导算等效塑性应变时会出错…】 σtrue，εtrue：之前得到的应力应变是“名义”应力应变，需要在此转换成真实应力应变。 表格中的红色部分为abaqus中的输入数据。

损伤模型为基础的钢筋混凝土构件

Damage Model Based Reinforced-Concrete Element 损伤模型为基础的钢筋混凝土构件 Chee Kiong Soh1; Yu Liu2; Yue Xing Dong3; and Xin Zheng Lu4 Abstract: In this paper, the damage in concrete around a reinforcement bar is examined, and three types of damage in concrete and at the concrete-rebar interface are defined. Based on these definitions, a three-dimensional damaged reinforced-concrete ~DRC! element is proposed. The element consists of a 10-node concrete element, a two-node rebar element, and a four-node concrete-rebar interface element. Experiments were carried out to obtain the parameters for the bond damage evolution equation. The proposed DRC element is implemented as a user-defined element of MSC MARC2000, a commercial finite-element analysis software package. Two numerical examples were examined, and the results show that the DRC element can simulate the bond deterioration in reinforced-concrete structural members with sufficient rationality, accuracy, and efficiency. DOI:10.1061/~ASCE!0899-1561~2003!15:4~371! CE Database subject headings: Damage; Models; Concrete, reinforced. 摘要： 在本论文中，讨论了加强钢筋周围的混凝土的破坏情况，将混凝土钢筋破坏界面的界定为3种类型。根据这些定义，提出了高强混凝土的三维破坏原理。这个原理包含三个元素，分别是一个10节点的混凝土构件，两节点筋单元，和一个四节点的混凝土钢筋接口。进行了试验，以取得弯矩破坏公式的有关参数。数据控制记录被建议作为MSC MARC2000的一个用户定义元素，MSC MARC2000是一个商用有限元分析软件包。两个数值实例研究结果表明，数据控制记录元素可以合理精确有效的模拟高强混凝土黏结失效。 DOI: 10.1061/~ASCE!0899-1561~2003!15:4~371! CE数据库主题词：破坏;模型;混凝土，加固。 Introduction

简支梁的地震响应分析

简支梁的地震响应分析 /PREP7 !进入前处理模块 /TITLE, EX 8.4(3) by Zeng P, Lei L P, Fang G ET,1,BEAM3 !设定1号单元 L=240 $A=273.9726 $H=14 $I=1000/3 !设定几何参数 R,1,273.9726,(1000/3),14 !设定1号实常数(梁单元) MP,EX,1,3E7 $MP,PRXY,1,0.3 $MP,DENS,1,73E-5 !设定弹性模量, 泊松比, 密度 K,1,0,0 $K,2,L,0 !生成两个关键点 L,1,2 !由关键点生成线 ESIZE,,8 !设定单元网格划分的分段数 LMESH,1 !对1号线划分单元网格 NSEL,S,LOC,X,0 !选择位置x=0的节点 D,ALL,UY !对所选择的节点施加位移约束UY=0 NSEL,S,LOC,X,L !选择位置x=L的节点 D,ALL,UX,,,,,UY !对所选择的节点施加位移约束UX=UY=0 NSEL,ALL !选择所有节点 FINISH !结束前处理模块 /SOLU !进入求解模块 ANTYPE,MODAL !设定模态分析方式 MODOPT,REDUC,,,,3 !设置缩减算法，提取3阶模态 MXPAND,1,,,YES ! 设定模态扩展的阶数为1，并计算单元及支反力结果 M,ALL,UY !对所有节点定义主自由度UY OUTPR,BASIC,1 !设置输出结果的方式 SOLVE !进行求解 *GET,F1,MODE,1,FREQ !提取第一阶模态频率，赋给F1 FINISH !结束 /SOLU !进入求解模块 ANTYPE,SPECTR !设定谱分析方式 SPOPT,SPRS !设定单点激励谱分析 SED,,1, !设定单点激励的方向为Y轴 SVTYP,3 !指定单点响应谱类型为地震位移谱 FREQ,.1,10 !设定频率数据表格的频率点 SV,,.44,.44 !设定频率数据表格的对应于频率点的激励值SOLVE !进行求解 *GET,F1_COEF,MODE,1,MCOEF !提取模态1的谱分析结果的模态系数FINISH !结束求解 /POST1 !进入一般性后处理模块 SET,1,1,F1_COEF !调出第1阶模态的结果，并乘以模态系数PRNSOL,DOF !打印节点结果 PRESOL,ELEM !打印单元结果 PRRSOL,F !打印支反力结果

混凝土损伤理论的分析研究

SHANGHAI UNIVERSITY 结构非线性分析课程论文 UNDERGRADUATE PROJECT (THESIS) 题 目:钢筋混凝土结构有限元分析及其断裂损伤理 论应用 学 院 土木工程系 专 业 建筑与土木工程 学 号 xxxxxxxx 学生姓名 xxx 指导教师 xx 日 期 2017.12.24

上海大学2017～2018学年冬季学期研究生课程考试 小论文 课程名称：结构非线性分析课程编号：18Z147004 论文题目:钢筋混凝土结构有限元分析及其断裂损伤理论应用 研究生姓名: xxx 学号: xxxxxxxx 论文评语: 成绩: 任课教师: xx 评阅日期:

目录 一混凝土损伤理论的研究背景 (1) 二国内外对混凝土损伤理论的研究现状 (2) 1）国外混凝土损伤理论研究现状 (2) 2）国内混凝土研究现状 (2) 三混凝土损伤理论研究中的问题和研究方法 (3) 1）试验条件相差较大时混凝土的本构关系将发生变化 (3) 2）复杂的多轴应力状态下的损伤理论 (3) 3）试验难度大 (3) 4）研究方法 (3) 四钢筋混凝土非线性损伤理论及有限元法 (4) 1）混凝土非线性本构模型 (4) 2）规范中的混凝土损伤理论 (5) ①混凝土单轴受压时的本构模型及dc的选取 (5) ②混凝土单轴受拉时的损伤理论 (6) 2）ABAQUS算例 (6) ①混凝土塑形损伤模型 (6) ②数值分析 (7) 五研究成果与创新 (8) 1)当今国际的研究成果 (8) 2)理论研究的新进展 (8) 3)在有限元中的应用 (8) 六研究混凝土损伤理论的意义和结论 (9) 1)社会意义 (9) 2)经济效益 (9) 3)结论 (9) 七展望 (9) 八建议 (10)

ANSYS地震响应分析讨论

地震响应分析 1模态组合就是根据模态分析中的几阶振型（也可以少于这几阶，看你要求的精度）进行组合（类似于结构最不利组合），从而求出地震响应的最大值。 2组合各振型反应的最大值，求得结构地震响应的最大值。 这个问题在论坛上已经有很多人问过，也有各种各样的回答，但是至今没有令人满意的解答。我自己试过很多种方法，加上论坛上其他人提到的方法，大致归类如下： 1.先做静力恒载工况分析，打开预应力pstres开关；然后转到时程分析。 结果：恒载对后面的时程计算不起作用，时程计算依然从0开始。 2.直接在antype，trans中考虑恒载：先把timint，off加acel,,9.81,打开应力刚化，sstif，on，lswrite，1，然后timint，on开始时程计算。 结果：恒载9.81起作用了，但结果是错的，它被积分了。 3.不用什么前处理，直接把9.81加在地震波上acel，9.81+ac（i）。 结果，同2，9.81带入了积分，这个9.81相当于阶跃荷载，而不是产生恒载。 4.ansys帮助中施加初始加速度的方法（篇幅限制请自己看帮助）。 结果，同2、3，9.81还是带进时间积分。 5.这种是我受到别人的启发，通过结构受ramp荷载的特点施加的，可以近似的解决问题。 即1）求出结构的自振一阶频率w 2）令tr＝1/w 3) 定义ramp荷载为从0到tr加到9.81，然后在整个时间积分中保持不变 4）antype,trans中分几个荷载步将荷载从0加到9.81 5) 在随后的荷载步中acel，，9.81+ac（i） 这种做法虽然也是将9.81++加到地震波中，但是因为满足TR的要求，所以这个动力效应被削弱到了静力效应，它作用在结构上就像静载一样。对于单自由度结构理论上跟静载是完全一样的，但是多自由度会子静力效应上下很小的范围内波动，所以可以认为相当于静载的作用，这样我们就可以达到考虑恒载的目的了。 第5种是我至今为止考虑恒载的做法，我也很想知道还有没有更简单精确的方法，或者在前4种方法中就有只是我使用不正确，希望大家能一起来讨论，彻底解决这个问题。谢谢！ 地震反应怎么考虑重力 SOLU ANTYPE, TRANS TRNOPT,FULL TIMINT,OFF !*先关闭时间积分效应 TIME,1E-8 !*设一个极短的积分时间 acel,,9.8 NSUBST,2 !有时候子步数要增大 KBC,1 LSWR,1 !*把这个写入第一步 TIMINT,ON !*然后再时间积分效应开关，以后就正常写载荷步了 这种方法应该是对的,ANSYS帮助文件中也有提到, 可是,有一个问题:由于是阶跃荷载,就会产生动力效应,整个结构的变形大于实际的情况吧?这样与实际结构在重力下受到的变形就不一样了!

地震反应分析：动力方法

地震反应分析：动力方法Structural Response Analysis: Dynamic Methods 教师：李爽副教授 lleshuang@https://www.doczj.com/doc/048246305.html, 2015年4月10日 1

本章导读 ?多维动力分析输入的一般处理方法 ?多次动力分析结果的一般处理方法 ?增量动力分析法(Incremental Dynamic Analysis Method，IDA) ?云图分析方法(Cloud Analysis Method)?结构地震模拟振动台试验基本步骤 2

多维动力分析输入的一般处理方法?当结构采用三维空间模型等需要双向 （两个水平方向）或三向（两个水平一 个竖向）地震动输入时，其加速度峰值 可按1（水平1）：0.85（水平2）：0.65 （竖向）的比例调整 ?具体如何操作？ 3

4 多维动力分析输入的一般处理方法 （2）初步选择若干条地震动，将所选择地震动进行反应谱分析，并与设计反应谱绘制在一起 （3）计算结构振型参与质量达到XX %（如50%~90%）对应各周期点处的地震动谱值（或0.2T 1~1.5T 1）。检查各周期处的包络值与设计反应谱值相差是否不超过20%。如不满足，则回到第二步重新选择地震动 （4）将各地震动在主要周期点处各方向上的值，按1（水平1）：0.85（水平2）：0.65（竖向）加权求和，按该求和值从小到大的顺序输入地震动（仅针对振动台试验，数值 计算不用分先后顺序，因为后者没有损伤）（1）根据研究对象所在场地类型和设防烈度确定地震设计反应谱（加速度反应谱）

多次动力分析结果的一般处理方法 ?《规范》规定 特别不规则的建筑、甲类建筑和下表所列高度范围的高层建筑，应采用时程分析法进行多遇地震下的补充计算；当取三组加速度时程曲线输入时，计算结果宜取时程法的包络值和振型分解反应谱法的较大值。当取七组及七组以上的时程曲线时，计算结果可取时程法的平均值和振型分解反应谱法的较大值 5

混凝土塑性损伤模型 -ABAQUS

4.5.2 混凝土塑性损伤模型ABAQUS ABAQUS 材料库中也包括分析混凝的其它模型如基于弥散裂纹方法的土本构模型。他们分别是在ABAQUS/Standard “An inelastic constitutive model for concrete,” Section 4.5.1, 中的弥散裂纹模型和在ABAQUS/Explicit, “A cracking model for concrete and other brittle materials,” Section 4.5.3中的脆性开裂模型。 混凝土塑性损伤模型主要是用来为分析混凝土结构在循环和动力荷载作用下的提供一个普遍分析模型。该模型也适用于其它准脆性材料如岩石、砂浆和陶瓷的分析；本节将以混凝土的力学行为来演示本模型的一些特点。在较低的围压下混凝土表现出脆性性质，主要的失效机制是拉力作用下的开裂失效和压力作用下的压碎。当围压足够大能够阻止裂纹开裂时脆性就不太明显了。这种情况下混凝土失效主要表现为微孔洞结构的聚集和坍塌，从而导致混凝土的宏观力学性质表现得像具有强化性质的延性材料那样。 本节介绍的塑性损伤模型并不能有效模拟混凝土在高围压作用下的力学行为。而只能模拟混凝土和其它脆性材料在与中等围压条件（围压通常小于单轴抗压强度的四分之一或五分之一）下不可逆损伤有关的一些特性。这些特性在宏观上表现如下： ?单拉和单压强度不同，单压强度是单拉强度的10倍甚至更多； ?受拉软化，而受压在软化前存在强化； ?在循环荷载(压)下存在刚度恢复； ?率敏感性，尤其是强度随应变率增加而有较大的提高。 概论 混凝土非粘性塑性损伤模型的基本要点介绍如下： 应变率分解 对率无关的模型附加假定应变率是可以如下分解的： 是总应变率，是应变率的弹性部分，是应变率的塑性部分。 应力应变关系 应力应变关系为下列弹性标量损伤关系： 其中是材料的初始（无损）刚度，是有损刚度，是刚度退化变量其值在0（无损）到1（完全失效）之间变化，与失效机制（开裂和压碎）相关的损伤导致了弹性刚度的退化。在标量损伤理论框架内，刚度退化是各向同性的，它可由单个标量d来描述。按照传统连续介质力学观点，有效应力可定义如下： Cauchy应力通过标量退化变量（d）转化为有效应力

时程分析时地震波的选取及地震波的反应谱化

时程分析时地震波的选取及地震波的反应谱化 摘要：目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应 谱法。时程分析法作为补充计算方法，在不规则、重要或较高建筑中采用。进行 时程分析时，首先面临正确选择输入的地震加速度时程曲线的问题。时程曲线的 选择是否满足规范的要求，则需要首先将时程曲线进行单自由度反应计算，得到 其反应谱曲线，并按规范要求和规范反应谱进行对比和取舍。本文通过介绍常用 的数值计算方法及计算步骤，实现将地震加速度时程曲线计算转化成反应谱曲线，从而为特定工程在时程分析时地震波的选取提供帮助。 关键词：时程分析，地震波，反应谱，动力计算 1 地震反应分析方法的发展过程 结构的地震反应取决于地震动和结构特性。因此，地震反应分析的水平也是随着人们对 这两个方面认识的深入而提高的。结构地震反应分析的发展可以分为静力法、反应谱法、动 力分析法这三个阶段。在动力分析法阶段中又可分为弹性和非弹性（或非线性）两个阶段。[1] 目前，在我国和其他许多国家的抗震设计规范中，广泛采用反应谱法确定地震作用，其 中以加速度反应谱应用得最多。反应谱是指：单自由度弹性体系在给定的地震作用下，某个 最大反应量（如加速度、速度、位移等）与体系自振周期的关系曲线。反应谱理论是指：结 构物可以简化为多自由度体系，多自由度体系的地震反应可以按振型分解为多个单自由度体 系反应的组合，每个单自由度体系的最大反应可以从反应谱求得。其优点是物理概念清晰， 计算方法较为简单，参数易于确定。 反应谱理论包括如下三个基本假定：1、结构物的地震反应是弹性的，可以采用叠加原理 来进行振型组合；2、现有反应谱假定结构的所有支座处地震动完全相同；3、结构物最不利 的地震反应为其最大地震反应，而与其他动力反应参数，如最大值附近的次数、概率、持时 等无关。[1] 时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由 于此法是对运动方程直接求解，又称直接动力分析法。可直接计算地震期间结构的位移、速 度和加速度时程反应，从而描述结构在强地震作用下弹性和非弹性阶段的内力变化，以及结 构构件逐步开裂、屈服、破坏甚至倒塌全过程。 根据我国《建筑抗震设计规范》(GB5011-2010)(以下简称《抗规》)第5.1.2-3条要求，特 别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑，应采用时程分析法进行多 遇地震下的补充计算。此外《高层建筑混凝土结构技术规程》(JGJ3-2010) (以下简称《高规》)第4.3.4条也有相关要求。 2 时程分析时地震波的选取要求 在进行时程分析时，首先面临地震波选取的问题。所选的地震波需要符合场地条件、设 防类别、震中距远近等因素。《抗规》对于地震波的选取主要有以下几点要求： 1、当取三组加速度时程曲线输入时，计算结果宜取时程法的包络值和振型分解反应谱法 的较大值；当取七组及七组以上的时程曲线时，计算结果可取时程法的平均值和振型分解反 应谱法的较大值（其中实际强震记录的数量不应少于总数的2/3）。 2、弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计 算结果的65%，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计 算结果的80%。 3、多组时程曲线的平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数 曲线在统计意义上相符。根据规范条文说明，所谓“统计意义上相符”指的是，多组时程波的 平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在对应于结构主 要振型的周期点上相差不大于20%。但计算结果也不能太大，每条地震波输入计算不大于135%，平均不大于120%。 4、时程曲线要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间均要符合规

地震响应的反应谱法与时程分析比较

发电厂房墙体地震响应的反应谱法与时程分析比较 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.1 9 0.261 2.5 0.313 0.25 0.047 与RG1.60标准谱对应的两条人工波见文件rg160x.txt与rg160y.txt

2数值分析框图思路与理论简介 2.1理论简介 该问题主要牵涉到结构动力分析当中的时程分析和谱分析。时程分析是用于确定承受任意随时间变化荷载的结构动力响应的一种方法。谱分析是模态分析的扩展，是用模态分析结果与已知的谱联系起来计算模型的位移和应力的分析技术。 2.2 分析框架： 时程分析：在X和Z两个水平方向地震波作用下，提取底部跨中单宽上的剪力、弯矩值和顶部水平位移，并求出最大响应。 谱分析：先做模态分析，再求谱解，由于X和Z两个方向的单点谱激励，因此需进行两次谱分析，分别记入不同的工况最后组合进行后处理得出结够顶部水平位移、底部单宽上剪力和弯矩的最大响应。 3有限元模型与荷载说明 3.1 有限元模型 考虑结构的几何特性建立有限元模型，首先建立平面几何模型，并将模型进行合理的切割，采用plane42单元，使用映射划分网格的方法生产平面单元（XOY平面）。然后，采用solid45单元，设置拖拉方向的单元尺寸并清楚初始平面单元plane42，将平面单元进行拖拉，最后生成发电厂墙体的有限元立体几何模型。单元总数为6060个，总节点数为8174个，有限元模型如图2所示： 图2 发电厂墙体有限元模型 3.2 荷载说明

结构地震反应分析

结构地震反应分析 摘要：结构地震反应分析方法有很多，单自由度体系可以采用duhamel积分法，多自由体系可以采用振型分解法，和直接积分法。在工程实践中，根据建筑物的结构体系，抗震设防烈度，选择合适的方法，计算结构的动力特性和响应。本文对一个7层框架结构进行抗震计算，采用不同的计算方法计算结构动力特性和响应。 关键词：duhamel积分法多自由度体系振型分解法直接积分法 Structural seismic response analysis FeiJianWei Civil and traffic institute structural engineering 200820104470 Abstract: There are many methods for Structural seismic response analysis, single-degree-of-freedom system using duhamel integral method, more free system can use strikeout decomposition method, and the direct integral method. In engineering practice, according to the building of the structure types, the seismic fortification intensity, select the appropriate method to calculate the dynamic characteristics, and response. Article choose a 7 layers framework for earthquake-resistant calculation, using different calculation method to calculate the dynamic characteristics and response. Keywords: duhamel integral method ；multi-freedom system ；vibration mode decomposition method ；direct integral method 1 前言 建筑结构抗震设计首先要计算结构的地震作用，然后再求出结构和构件的地震作用效应。结构的地震作用效应就是指地震作用在结构中所产生的内力和变形，主要有弯矩、剪力、袖向力和位移等，最后将地震作用效应与其他荷载效应进行组合，并验算结构和构件的抗震承载力及变形，以满足“小震不坏，中震可修，大震不倒”的抗震设计要求。 结构的地震反应是指地震引起的结构振动，它包括地震在结构中引起的速度、加速度、位移和内力等。结构的地震反应分析属于结构动力学的范畴，比结构的静力分析要复杂得多。因为结构的地震反应不仅与地震作用的大小及其随时间的变化特性有关，而且还取决于结构本身的动力特性，即结构的自振周期和阻尼等。然而，地震时地面的运动是一种很难确定的随机过程，运动极不规则，而建筑结构又是一个由各种不同构件组成的空间体系，其动力特性也十分复杂。因此，地震引起的结构振动实际上是一种很复杂的空间振动。这样，在进行建筑结构的地震反应分析时，为了便于计算，常需做出一系列简化的假定[1]。 1.1 结构抗震理论的发展 近百年来，经过各国学者的共同努力，结构抗震理论的研究取得了长足的发展。结构抗震理论的发展可以划分为静力理论、反应谱理论和动力理论三个发展阶段。 1.1.1 静力理论 水平静力抗震理论创始于意大利，发展于日本，1900年日本学者大森房吉提出震度法的概念。该理论认为：结构物所受到的地震作用，可以简化为作用于结构的等效水平静力F，其大小等于结构重力荷载G乘以地震系数k，即： F =α G / g = kG(1.1)