在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)
在ABAQUS中对框架结构施加地震波（对初学者普及，同时向大虾们求教）
初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。但是我只学了石老师《实例》的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。
在找完很多帖子后，建模并计算，基本算是成功的。所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。
不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。
但同时还有许多自己不知其所以然的，想要请教各位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答）。
在网上查了些方法：
module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界
---done----出现对话框edit bondary condition对话框，在amplitude 里选择你所定义的时间和加速度。点ok就完工了。
这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。
首先在CAE里建好模型，定义两个分析步。
第一个分析步是加自重，采用线性加载的方式。
(a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous （瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。
(b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。
#1
liudaomuxi ng
助理工程师
精华0
积分16
帖子12
水位15
技术分0 第二个分析步就是加地震波。
输入地震波有两种方法：
1
波
据
%%%%%%%%%%%%%%%%%%%%%%
*End Assembly
*Amplitude, name=Amp-1
0.005, -7.5e-08, 0.01, -3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-07 0.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-05
0.045, 3.54e-05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-05………………%%%%%%%%%%%%%%%%%%%%%%
2
种是直接把地震波放在
逗号的（
这种方法是将地震波另外放置在一个
件得跟
%%%%%%%%%%%%%%%%%%%%%%
*End Assembly
*Amplitude, name=Amp-1, definition=SMOOTH STEP
0., 0., 0.5, 0.5, 1., 1.
**
*Amplitude, name=Amp-2, input=DIZHENBO--E-W.txt
**
** MATERIALS %%%%%%%%%%%%%%%%%%%%%%
其中
%%%%%%%%%%%%%%%%%%%%%%
0,0.0650966,0.02,0.0485266,0.04,-0.00118357,0.06,-0.094686, 0.08,-0.102971,0.1,-0.0514855,0.12,-0.0153865,0.14,0.0224879, 0.16,0.0432005,0.18,0.0538527,0.2,0.0420169,0.22,-0.00532609, 0.24,-0.0426087,0.26,-0.0467512,0.28,-0.0272222,0.3,-0.00769324,………………%%%%%%%%%%%%%%%%%%%%%%
#2
liudaomuxi
ng
助理工程师
精华 0
积分 16
帖子 12
水位 15
技术分 0 对于边界条件，在初始步时将柱底约束住： %%%%%%%%%%%%%%%%%%%%%% ** BOUNDARY CONDITIONS ** ** Name: BC-1 Type: Acceleration/Angular acceleration *Boundary, type=ACCELERATION _PickedSet1399, 1, 1 _PickedSet1399, 2, 2 _PickedSet1399, 3, 3 %%%%%%%%%%%%%%%%%%%%%% 但关键在于这里边界条件的类型是Displacement 还是如上所示的Acceleration 我不能确定（第二步中的边界得是Acceleration ，
第一步好像是不可以用Displacement 的，不知道对不对）。 在第二个分析步中，将地震波以加速度的形式加在柱底上： %%%%%%%%%%%%%%%%%%%%%%
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Acceleration/Angular acceleration
*Boundary, type=ACCELERATION
_PickedSet1399, 2, 2
_PickedSet1399, 3, 3
*Boundary, amplitude=Amp-2, type=ACCELERATION
_PickedSet1399, 1, 1, 1
%%%%%%%%%%%%%%%%%%%%%%
其中最后的那个“1”参数是指将地震波数据乘以此系数得到最终地震波的值。这里我有点不明白的是蓝色字体部分，这里不加这些语句可以吗？因为在初始步中已经定义过了，但我不知道这时是不是得重新把2、3方向的给再次定义下。（我知道有些东西我自己改下再运算一遍比较下结果也许可以知道，但一来不太明白的地方太多，初学也不知道改的对不对，二来我觉得这样可能知道结果，却不知道原因，所以发上来希望有大侠告知）
Inp 文件添加到这里也就可以放入Command 里进行计算了。 结果如附件所示：
20倍变形图
位移曲线（在ABAQUS 中将顶点跟底点位移曲线相减得到相对位移）
我想介绍的就是这些，希望对于需要用到加载地震波的各位能有所帮助。
同时，这里面除了我希望有疑问的地方能得到大侠们解答之外，如有其他错误或不足之处，请大家指出。谢谢！O(∩_∩)O 哈哈~
在 ansys 中如何 施加 地震波
三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据： -0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据： 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件，并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行，同上 ACCEXYZ(0,3)=3 !08行，同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1 acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度 SOLVE
Abaqus中应力应变的理解
在ABAQUS 中对应力的部分理解 1、三维空间中任一点应力有6个分量yz xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点的任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用的面为主截面，其上的正应力为主应力，主截面的法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max. Principal 、Mid. Principal 、Min. Principal ，这三个量在任何坐标系统下都是不变量。 可利用最大主应力判断一些情况：比如混凝土的开裂，若最大主应力（拉应力）大于混凝土的抗拉强度，则认为混凝土开裂，同时通过显示最大主应力的法线方向，可以大致表示出裂缝的开裂方向等。 利用最小主应力，可以查看实体中残余压应力的大小等。 3、弹塑性材料的屈服准则 3.1、Mises 屈服准则 22 13232 2 212)()()(S σ σσσσ σσ=-+-+- 其中s σ为材料的初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力的定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就是我们常见的 )(31z y x p σσ σ ++= 。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑 性变形形状的变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义的不同而变化）S11，S22，S33，S12，S13，S23 3.2、Trasca 屈服准则 主应力间的最大差值=2k
时程分析中地震波输入位置的讨论
时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[1]。目前有限元软件可以实现结构的时程分析，但是在不同的软件中，其实现方式不同，主要区别在地震波的输入位置不同。本文通过有限元软件ABAQUS采用不同的地震波输入位置对同一结构进行时程分析分析，对比结构相同位置的时程位移曲线，结果表明结构在采用不同地震波输入位置的时程分析中，结构的地震响应基本一致。 关键词：时程分析、有限元软件、钢筋混凝土剪力墙 Abstract: The time history analysis method to analyze the available structure through direct power to the relationship between the corresponding changes over time, truly reflect the structure of earthquake corresponding to the whole process of change over time, is an important method of seismic analysis [1]. Finite element software can be time-history analysis of the structure, but in different software in different ways, the main difference between the different positions in the seismic wave input. In this paper the finite element software ABAQUS using different seismic wave input location on the same structure, process analysis analysis, contrast structure the same location of when the process displacement curve, the results show that the structure using different seismic waves enter the position time history analysis, the seismic response basically the same. Keywords: time history analysis, finite element software, reinforced concrete shear walls 一、引言 在时程分析等动力学问题中，地震力以加速度形式从基础固定处输入。由于结构的刚度不是无限大，在结构上的加速度反应与基础输入的加速度并不相同。在很多时候，结构的加速度比基础输入的加速度更大，即对输入的加速度有一个动力放大效应。在单自由度弹性体系中，体系最大绝对加速度与地面运动最大加速度的比值，即称为动力系数[2] （1） 动力系数与结构的动力学特性和输入的地震波的频率特性有关。它与地震系数k的乘积即为单自由度体系的地震影响系数。 因此，从原理上讲，时程分析是将地震波的加速度时程曲线作用到结构的基础约束处，得到上部结构的各种地震反应。但是在不同的软件中，其实现方
ABAQUS定义真实应力和真实应变
ABAQUS 中定义真实应力和真实应变 在ABAQUS 中必须用真实应力和真实应变定义塑性.ABAQUS 需要这些值并对应地在输入文件中解释这些数据。 然而，大多数实验数据常常是用名义应力和名义应变值给出的。这时，必须应用公式将塑性材料的名义应力（变）转为真实应力（变）。 考虑塑性变形的不可压缩性，真实应力与名义应力间的关系为： 00l A lA =， 当前面积与原始面积的关系为： 00l A A l = 将A 的定义代入到真实应力的定义式中，得到： 00 ()nom F F l l A A l l σσ=== 其中0 l l 也可以写为1nom ε+。 这样就给出了真实应力和名义应力、名义应变之间的关系： (1)nom nom σσε=+ 真实应变和名义应变间的关系很少用到，名义应变推导如下： 0001nom l l l l l ε-= =- 上式各加1，然后求自然对数，就得到了二者的关系： ln(1)nom εε=+ ABAQUS 中的*PLASTIC 选项定义了大部分金属的后屈服特性。ABAQUS 用连接给定数据点的一系列直线来逼近材料光滑的应力－应变曲线。可以用任意多的数据点来逼近实际的材料性质;所以，有可能非常逼真地模拟材料的真实性质。在*PLASTIC 选项中的数据将材料的真实屈服应力定义为真实塑性应变的函数。选项的第一个数据定义材料的初始屈服应力，因此，塑性应变值应该为零。 在用来定义塑性性能的材料实验数据中，提供的应变不仅包含材料的塑性应变，而是包括材料的总体应变。所以必须将总体应变分解为弹性和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值，从总体应变中减去弹性应变，就得到了塑性应变，其关系为： /pl t el t E ε εεεσ=-=- 其中pl ε是真实塑性应变，t ε是总体真实应变，el ε是真实弹性应变。
abaqus中显示动力学分析步骤
准静态分析——ABAQUS/Explicit 准静态过程（guasi-static process） 在过程进行的每一瞬间，系统都接近于平衡状态，以致在任意选取的短时间dt 内，状态参量在整个系统的各部分都有确定的值，整个过程可以看成是由一系列极接近平衡的状态所构成，这种过程称为准静态过程。无限缓慢地压缩和无限缓慢地膨胀过程可近似看作为准静态过程。准静态过程是一种理想过程，实际上是办不到的。 准静态原为一个热力学概念，在这里引用主要是指模型在加载的过程中任意时刻所经历的中间状态都可近似地视为静力状态，因此当加载过程进行得无限缓慢时，在各个时刻模型所处的状态就可近似地看作是静态，该过程便是准静态过程。准静态啮合过程仿真主要考虑的是弧齿锥齿轮副在加载时的接触状态，以及齿面和齿根的应力变化规律，其前提是不考虑齿轮副惯性的影响。 ABAQUS/Explicit准静态分析 显式求解方法是一种真正的动态求解过程，它的最初发展是为了模拟高速冲击问题，在这类问题的求解中惯性发挥了主导性作用。当求解动力平衡的状态时，非平衡力以应力波的形式在相邻的单元之间传播。由于最小稳定时间增量一般地是非常小的值，所以大多少问题需要大量的时间增量步。 在求解准静态问题上，显式求解方法已经证明是有价值的，另外ABAQUS/Explicit在求解某些类型的静态问题方面比ABAQUS/Standard更容易。在求解复杂的接触问题时，显式过程相对于隐式过程的一个优势是更加容易。此外，当模型很大时，显式过程比隐式过程需要较少的系统资源。 将显式动态过程应用于准静态问题需要一些特殊的考虑。根据定义，由于一个静态求解是一个长时间的求解过程，所以在其固有的时间尺度上分析模拟常常在计算上是不切合实际的，它将需要大量的小的时间增量。因此，为了获得较经济的解答，必须采取一些方式来加速问题的模拟。但是带来的问题是随着问题的加速，静态平衡的状态卷入了动态平衡的状态，在这里惯性力成为更加起主导作用的力。目标是在保持惯性力的影响不显著的前提下用最短的时间进行模拟。
abaqus如何施加地震波
施加地震波： 1 *amplitude,name=amp,input=seismicdata.dat 输入地震波 2 *boundary,type=acceleration,amplitude=amp施加荷载 方法：module选load，在tools-----amplitude-----creat默认的continue在Edit A mplitude里面输入时间和加速度，点OK。点creat boundary condition，涌现对 话框creat boundary condition，选择acceleration/angular acceleration,continu e---选择要施加的边界---done----涌现对话框edit bondary condition对话框，在 amplitude里选择你所定义的时间和加速度。点ok就完工了。 在网上查了些方法： module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude 里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框，在amplitude里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型，定义两个分析步。 第一个分析步是加自重，采用线性加载的方式。 (a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。 (b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。 第二个分析步就是加地震波。 输入地震波有两种方法： 1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波,时间,地震波,时间, 地震波,时间, 地震波…………每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）。%%%%%%%%%%%%%%%%%%%%%% *End Assembly *Amplitude, name=Amp-1 0.005, -7.5e-08, 0.01, -3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-07 0.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-05 0.045, 3.54e-05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-05 ………………
ABAQUS建模如何施加预应力(残余应力)
ABAQUS建模如何施加预应力 本文参考了百度文库中的文章： https://www.360docs.net/doc/924209029.html,/link?url=dt_VLOGCUf8hUo7A9THhyv7BuSHry71EbLVtBtkWpoiYtkm Lxbfk0Io63jsygs6vWbFU7x22HHFv8pIGgPMYkv1lyFXWbgPJqvCodSioUqa 关键字格式： “*initial conditions, type=stress, input=bb.dat” 上面的关键字，即绿色部分，全部插于*STEP语句之前（如下图），两语句之间不能有空格。施加预应力场只是initial conditions关键字的一个应用，详见abaqus6.8帮助文档，《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。 实例：点焊所产生的焊点中存在着残余应力，本文就是教大家如何完成焊点残余应力的模拟。原理说明：先在模型上施加一个任意载荷（记为状态1），可得出此载荷作用下模型上的等效节点载荷，然后通过keywords让这个等效节点载荷作用于模型上，使它与之前施加在模型上的载荷相平衡，便得到了一个位移为0的初始状态（记为状态2），该状态下，模型中含有的应力场与状态1相同，只是模型不发生变形而已。 具体操作步骤： 1、建立有限元模型，部件类型为轴对称
2、设置材料常数（自己任意设） 3、分析步，设置两个分析步 4、设置任意一个自己需要的载荷，此载荷即为与初始应力对应的载荷。让该载荷从分析步2开始作用。分析步1空着。原因不明。 （让载荷从step2开始，得到的分析结果图显示很光滑。但若是让载荷从step1开始并延续到step2，或者从没有step2 的step1开始，得到的分析图都不是光滑，原因未知。） 设置边界条件，从状态Initial开始。
ABAQUS时程分析实例
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.
ABAQUS瑞利阻尼
关于ABAQUS中的质量比例阻尼 总结论： ABAQUS中的质量比例阻尼是和绝对速度有关的，即质量比例阻尼产生的阻尼力由绝对速度引起。以阻尼系数表达的阻尼，产生的阻尼力由相对速度引起。 Abaqus Analysis User's Manual—Mass proportional damping：The factor introduces damping forces caused by the absolute velocities of the model and so simulates the idea of the model moving through a viscous “ether” (a permeating, still fluid, so that any motion of any point in the model causes damping).帮助手册也说明了质量比例阻尼是和绝对速度有关。 问题： 1、应用直接积分法进行时程分析，地震波一般以边界条件 的形式加到支座处，结构阻尼只能使用Rayleigh阻尼， 而这时产生的阻尼力是绝对速度产生的，而运动方程中 的阻尼项产生的阻尼力是与相对速度有关。 2、SAP2000中施加地震波，支座处相对位移为0，绝对位 移不为0，其相对位移相对哪一点来说的？
算例： 单自由度体系，如图（1），质量m=，k=1N/m，阻尼比ξ=，对应的阻尼系数c=，若应用直接积分法进行时程分析，结构的阻尼需要转换成Rayleigh阻尼，使用如下公式： 22 n n n βω α ξ ω =+ 如果只使用质量比例阻尼（结构只有一阶振型），即 2 n n α ξ ω =，容易得出α=。. 图（1） 情况（1）： 在ABAQUS中用spring单元模拟竖向的直杆，水平刚度k=1N/m，采用Rayleigh阻尼，通过*mass,alpha=（质量比例阻尼）施加，地震波需用Elcentrol波，以边界条件的形式加在支座处(竖向杆下端)。 为了作对比，在SAP2000中的结构阻尼在分析工况中以质量比例阻尼的形式施加。MATLAB中变成使用NewMark-beta方法。
最新Abaqus中应力应变的理解
在ABAQUS 中对应力的部分理解 1、三维空间中任一点应力有6个分量y z xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点的任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用的面为主截面，其上的正应力为主应力，主截面的法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max. Principal 、Mid. Principal 、Min. Principal ，这三个量在任何坐标系统下都是不变量。 可利用最大主应力判断一些情况：比如混凝土的开裂，若最大主应力（拉应力）大于混凝土的抗拉强度，则认为混凝土开裂，同时通过显示最大主应力的法线方向，可以大致表示出裂缝的开裂方向等。 利用最小主应力，可以查看实体中残余压应力的大小等。 3、弹塑性材料的屈服准则 3.1、Mises 屈服准则 22132322212)()()(S σσσσσσσ=-+-+- 其中s σ为材料的初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力的定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就是我们常见的)(3 1z y x p σσσ++=。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑性 变形形状的变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义的不同而变化）S11，S22，S33，S12，S13，S23 3.2、Trasca 屈服准则 主应力间的最大差值=2k
ABAQUS施加预应力
ABAQUS施加预应力 参考百度文库： https://www.360docs.net/doc/924209029.html,/link?url=dt_VLOGCUf8hUo7A9THhyv7BuSHry71EbLVtBtkWpoiYtkm Lxbfk0Io63jsygs6vWbFU7x22HHFv8pIGgPMYkv1lyFXWbgPJqvCodSioUqa 关键字格式： “*initial conditions, type=stress, input=bb.dat” 上面的关键字，即绿色部分，全部插于*STEP语句之前（如下图），两语句之间不能有空格。施加预应力场只是initial conditions关键字的一个应用，详见abaqus6.8帮助文档，《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。 实例：点焊所产生的焊点中存在着残余应力，本文就是教大家如何完成焊点残余应力的模拟。原理说明：先在模型上施加一个任意载荷（记为状态1），可得出此载荷作用下模型上的等效节点载荷，然后通过keywords让这个等效节点载荷作用于模型上，使它与之前施加在模型上的载荷相平衡，便得到了一个位移为0的初始状态（记为状态2），该状态下，模型中含有的应力场与状态1相同，只是模型不发生变形而已。 具体操作步骤： 1、建立有限元模型，部件类型为轴对称
2、设置材料常数（自己任意设） 3、分析步，设置两个分析步 4、设置任意一个自己需要的载荷，此载荷即为与初始应力对应的载荷。让该载荷从分析步2开始作用。分析步1空着。原因不明。 （让载荷从step2开始，得到的分析结果图显示很光滑。但若是让载荷从step1开始并延续到step2，或者从没有step2 的step1开始，得到的分析图都不是光滑，原因未知。） 设置边界条件，从状态Initial开始。
[土木] 在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)
[土木]在ABAQUS中对框架结构施加地震波（对初学者普及，同时向大虾们求教） 初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。但是我只学了石老师《实例》的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。 在找完很多帖子后，建模并计算，基本算是成功的。所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。 不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。 但同时还有许多自己不知其所以然的，想要请教各位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答）。 在网上查了些方法： module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框，在amplitude里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型，定义两个分析步。 第一个分析步是加自重，采用线性加载的方式。 (a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。 (b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。 以上方法也是得自论坛，我感觉说的有道理，就用了，不知道具体对不对。 第二个分析步就是加地震波。 输入地震波有两种方法： 1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波,时间,地震波,时间, 地震波,时间, 地震波…………每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）。 %%%%%%%%%%%%%%%%%%%%%%
abaqus后处理中各应力解释个人收集修订版
a b a q u s后处理中各应力 解释个人收集修订版 IBMT standardization office【IBMT5AB-IBMT08-IBMT2C-ZZT18】
ABAQUS中的壳单元S33代表的是壳单元法线方向应力，S11 S22 代表壳单元面内的应力。因为壳单元的使用范围是“沿厚度方向应力为0”，也即沿着法相方向应力为0，且满足几何条件才能使用壳单元，所以所有壳单元的仿真结果应力查看到的S33应力均为0。 S11 S22 S33 实体单元是代表X Y Z三个方向应力，但壳单元不是，另外壳单元只有S12，没有S13，S23。 LE----真应变（或对数应变） LEij---真应变 ... 应变分量； PE---塑性应变分量； PEEQ---等效塑性应变 ABAQUS Field Output Stresses S stress components and invariants 应力分量和变量 SVAVG volume-averaged stress components and invariants (Eulerian only) MISESMAX 最大 Mises 应力 TSHR transverse shear stress(for thick shells)横向剪切应力 CTSHR transverse shear stress in stacked continuum shells 连续堆垛壳横向剪切应力 TRIAX stress triaxiality 应力三轴度 VS stress in the elastic-viscous network 弹粘性网格应力 PS stress in the plastic-viscous
Abaqus中应力应变的理解
在ABAQUS 中对应力得部分理解 1、三维空间中任一点应力有6个分量y z xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点得任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用得面为主截面，其上得正应力为主应力，主截面得法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max 、 Principal 、Mid 、 Principal 、Min 、 Principal ，这三个量在任何坐标系统下都就是不变量。 可利用最大主应力判断一些情况：比如混凝土得开裂，若最大主应力（拉应力）大于混凝土得抗拉强度，则认为混凝土开裂，同时通过显示最大主应力得法线方向，可以大致表示出裂缝得开裂方向等。 利用最小主应力，可以查瞧实体中残余压应力得大小等。 3、弹塑性材料得屈服准则 3、1、Mises 屈服准则 22132322212)()()(S σσσσσσσ=-+-+- 其中s σ为材料得初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力得定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就就是我们常见得)(31z y x p σσσ++=。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑 性变形形状得变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义得不同而变化）S11，S22，S33，S12，S13，S23 3、2、Trasca 屈服准则 主应力间得最大差值=2k 若明确了321σσσ≥≥，则有k =-)(2 131σσ，若不明确就需要分别两两求差值，瞧哪个最大。 ABAQUS 中得Trasca 等效应力就就是“主应力间得最大差值” 3、3 ABAQUS 中得Pressure----等效压应力 即为上面提到得p ：, 也就就是我们常见得)(3 1z y x p σσσ++=。
abaqus问题集锦
:我老板要求我用的是concrete damage plasticity model 1。在做非线性分析的时候，怎么知道结构破坏了，如何得知什么时候获取了结构的极限承载力 2。如何确定concrete compression damage和concrete tension damage？二者似乎都无法通过试验得到。如果说可以忽略compression damage，那么由于tension stiffening的存在（下降段相对而言变长了），tension damage的定义就显得很重要了。在ABAQUS关于重力坝的例子中，tension stiffening和tension damage的数据，似乎有一些关系。 3。定义concrete tension stiffening的时候，我看一些论文上面关于tesion stiffening说明的时候，说如果tensiong stiffening 定义得太小，收敛会很难。但是tension stiffening的大小在用*concrete tension stiffening定义如何体先呢，在前两个数据是remaining direct stress after cracking和direct cracking strain. 4。tension stiffenig 与配筋率是否有关系呢？如果有的话，关系是怎么样的呢？谢谢大家了。 A:我认为damage plasticity model主要用于混凝土受循环往复荷载的情况，concrete compression damage和concrete tension damage也只在这时需要定义，tension stiffening好像就是定义混凝土单轴受拉本构关系，与tension damage应该没有关系。与配筋率好像也没有关系。damage plasticity model既可以模拟循环拟静荷载的情况，也可计算结构在动力荷载（地震、风以及冲击荷载的作用），通过concrete compression damage和concrete tension damage考虑混凝土的拉压异性，tension stiffening则可考虑拉伸屈服后的软还阶段本构行为，同样是混凝土损伤全过程的一部分，它与配筋率及网格划分密度的相互关系很大程度上影响到计算的收敛性。A:损伤主要针对循环加载问题，如果你不是循环加载的话可以不考虑损伤,我做过比较 A:楼上的大侠好象问题考虑的不是很全面吧.如果在单向加载时同时存在拉压,就应该考虑.而其中的恢复系数是循环加载下考虑损伤恢复A:我现在也在搞concrete damage plasticity model,希望能和各位大虾交流交流A:。在做非线性分析的时候，怎么知道结构破坏了，如何得知什么时候获取了结构的极限承载力? 这个问题牵涉到结构在损伤到多少时即认为是损坏.也就是通常说的怎么把结构算倒.不知道我理解的是否正确.如果是这个问题就牵涉到损伤系数的取法.因为是整个结构,就须要把所有的损伤系数加权平均.加权的取法也很有讲究,有人认为达到0.7即认为倒塌,可以参考王光远的文章.现有有限元程序还不能完成,是个急待解决的难题. 当然是有关系的了concrete compression damage和concrete tension damage分别是受压损伤系数和受拉损伤系数函数,在单项受力时是直接乘入混凝土刚度的.在往复荷载下要乘损伤恢复系数，在三轴应力下还要乘加权因子。考虑了该参数后，开裂应变就不等于塑性应变。开裂应变是考虑弹性模量折减的非弹性应变。因此此时塑性应变不等于非弹性应变。这是关键啊A:谢谢各位的回复, 帮助很大. 现在riks 去求解的时候,可以获取极限何在,结果与实验结果差别不大,但是在极限荷载下的位移与实验结果差别很大,无论用shell还是solid 来模拟钢筋混凝土的话, 位移误差都很大. 请问大家有没有碰到着个问题呢? abaqus问答精华 Q:预拉钢筋怎样施加预应力，请各位指点~~~~ Q:我在文档里看到要在inp文件定义一个rebar，但是rebar只能用于shell, membrane, and solid elements 。我现在想做的是一个预应力拉索，不是镶嵌在shell, membrane, and solid 这些单元里的，而是独立的一根拉锁。拉索单元打算用truss,但是怎样在truss上使用rebar啊？请高手指点还有个问题，我看到别人的inp文件，如下：*rebar,element=continuum,material=rebar2,name=ubar top1,1.005e-4,0.15,0.0,0.5,1 第二行第一
Abaqus中应力应变的理解
1、三维空间中任一点应力有6个分量y z xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点的任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用的面为主截面，其上的正应力为主应力，主截面的法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max. Principal 、Mid. Principal 、Min. Principal ，这三个量在任何坐标系统下都是不变量。 可利用最大主应力判断一些情况：比如混凝土的开裂，若最大主应力（拉应力）大于混凝土的抗拉强度，则认为混凝土开裂，同时通过显示最大主应力的法线方向，可以大致表示出裂缝的开裂方向等。 利用最小主应力，可以查看实体中残余压应力的大小等。 3、弹塑性材料的屈服准则 、Mises 屈服准则 22132322212)()()(S σσσσσσσ=-+-+- 其中s σ为材料的初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力的定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力，I 为单位矩阵，p 为等效压 应力（定义如下）：, 也就是我们常见的)(3 1z y x p σσσ++=。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑性变形形状的变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义的不同而变化）S11，S22，S33，S12，S13，S23 、Trasca 屈服准则 主应力间的最大差值=2k 若明确了321σσσ≥≥，则有k =-)(2 131σσ，若不明确就需要分别两两求差值，看哪个最大。 ABAQUS 中的Trasca 等效应力就是“主应力间的最大差值” ABAQUS 中的Pressure----等效压应力
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 计算对象 操作过程为： 1 打开ABAQUS/CAE，点击create model database。 2 进入Part模块，点击create part，命名为column，3D、deformation、wire。continue
3 Create lines，在 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.
5 Create section，name：Section-1，category：beam，type：beam, Continue
create profile, name: Profile-1, shape:I, 按图1尺寸输入界面尺寸，ok。 在profile name选择I，material name 选择steel。Ok
6 Assign section，选择全部，done，弹出的对话框选择section：Section-1，ok。 7 Assign beam orientation，选择全部，默认值确定。 8 View-->>part display options，在弹出的对话框里勾选，render beam profiles，以可视化梁截面形状。
在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)
在ABAQUS中对框架结构施加地震波（对初学者普及，同时向大虾们求教） 初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。但是我只学了石老师《实例》的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。 在找完很多帖子后，建模并计算，基本算是成功的。所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。 不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。 但同时还有许多自己不知其所以然的，想要请教各位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答）。 在网上查了些方法： module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界 ---done----出现对话框edit bondary condition对话框，在amplitude 里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型，定义两个分析步。 第一个分析步是加自重，采用线性加载的方式。 (a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous （瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。 (b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。