[土木] 在ABAQUS中对框架结构施加地震波(对初学者普及,同时向大虾们求教)

[土木]在ABAQUS中对框架结构施加地震波（对初学者普及，同时向大虾们求教）初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。但是我只学了石老师《实例》的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。在找完很多帖子后，建模并计算，基本算是成功的。所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。但同时还有许多自己不知其所以然的，想要请教各位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答）。在网上查了些方法：module选load，在tools-----amplitude-----creat默认的continue在EditAmplitude里面输入时间和加速度，点OK。点creatboundarycondition，出现对话框creatboundarycondition，选择acceleration/angularacceleration,continue---选择要施加的边界---done----出现对话框editbondarycondition对话框，在amplitude里选择你所定义的时间和加速度。点ok就完工了。这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。首先在CAE里建好模型，定义两个分析步。第一个分析步是加自重，采用线性加载的方式。(a)加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。(b)新建加载方式：创建一个新的Amplitude，Type=smoothtpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smoothtpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。以上方法也是得自论坛，我感觉说的有道理，就用了，不知道具体对不对。第二个分析步就是加地震波。输入地震波有两种方法：1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波,时间,地震波,时间,地震波,时间,地震波…………每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）。%%%%%%%%%%%%%%%%%%%%%%*EndAssembly*Amplitude,name=Amp-10.005,-7.5e-08,0.01,-3.55e-07,0.015,-7.03e-07,0.02,-4.53e-070.025,1.82e-06,0.03,7.01e-06,0.035,1.5e-05,0.04,2.49e-050.045,3.54e-05,0.05,4.5e-05,0.055,5.2e-05,0.06,5.5e-05………………%%%%%%%%%%%%%%%%%%%%%%2、在如下位置加入下面加黑的字体部分。区别就在于上面那一种是直接把地震波放在inp文件里，而且每一行的最后面是没有逗号的（我也不知道这是不是必须的，帮助文件里没找着）。而这种方法是将地震波另外放置在一个txt文件里，当然，这个文件得跟inp文件在同一个文件夹里。%%%%%%%%%%%%%%%%%%%%%%*EndAssembly*Amplitude,name=Amp-1,definition=SMOOTHSTEP0.,0.,0.5,0.5,1.,1.***Amplitude,name=Amp-2,input=DIZHENBO--E-W.txt****MATERIALS%%%%%%%%%%%%%%%%%%%%%%其中DIZHENBO--E-W.txt是地震波文件。%%%%%%%%%%%%%%%%%%%%%%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,………………%%%%%%%%%%%%%%%%%%%%%%对于边界条件，在初始步时将柱底约束住：%%%%%%%%%%%%%%%%%%%%%%**BOUNDARYCONDITIONS****Name:BC-1Type:Acceleration/Angularacceleration*Boundary,type=ACCELERATION_PickedSet1399,1,1_PickedSet1399,2,2_PickedSet1399,3,3%%%%%%%%%%%%%%%%%%%%%%但关键在于这里边界条件的类型是Displacement还是如上所示的Acceleration我不能确定（第二步中的边界得是Acceleration，第一步好像是不可以用Displacement的，不知道对不对）。在第二个分析步中，将地震波以加速度的形式加在柱底上：%%%%%%%%%%%%%%%%%%%%%%**BOUNDARYCONDITIONS****Name:BC-1Type:Acceleration/Angularacceleration*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哈哈~原创：ABAQUS地震时程分析的一些经验总结做了挺长时间地震，自己做了个总结，省得忘记，也希望能给大家提供参考，减少摸索的时间。Abaqus地震分析的总结——时步、单元尺寸、滤波、等效非线性、无限元1、自由场地震反应经典的自由场地震反应用shake91或proshake等进行分析，在分析完可以导出各层土的等效线性参数，包括阻尼（粘滞阻尼系数）和剪切模量，用剪切模量可以计算弹性模量，shake中假定泊松比为常数，对地震反应没影响。其实FLAC中有自带自由场边界，计算地震很方便。如以下是shake91中自带的例子输出得到的参数ITERATIONNUMBER8VALUESINTIMEDOMAINNOTYPEDEPTHUNIFRM.----DAMPING--------SHEARMODULUS-----G/Go(FT)STRAINNEWUSEDERRORNEWUSEDERRORRATIO------------------------------------------------------------122.5.00077.007.007.03851.53851.5.0.992227.5.00295.014.014.03020.03020.0.0.9603215.0.00634.023.023.02803.82803.8.0.8924225.0.00976.028.028.02985.82985.8.0.8525135.0.01099.030.030.03621.73621.6.0.9336145.0.01403.035.035.03540.53540.4.0.9127155.0.01362.034.034.04296.04296.0.0.9158165.0.01566.037.037.04239.84239.8.0.9039275.0.01356.034.034.05402.75402.7.0.79210285.0.01505.037.037.05266.05266.0.0.77211295.0.01336.034.034.06288.26288.2.0.795122105.0.01413.035.035.06203.46203.4.0.784132115.0.01233.032.032.07357.27357.2.0.810142125.0.01281.033.033.07290.87290.8.0.803152135.0.01115.030.030.08570.48570.4.0.829162145.0.00865.026.026.011292.611292.6.0.863Shake中的outcrop指出露基岩，baserock指土层底部的基岩，因此不考虑波的衰减的情况下，在outcrop处用加速度计测得的地震加速度幅值为baserock处的2倍。如果在abaqus中的土层与基岩的交界面处输入地震波，在跟shake91对比时，要是用baserock输入。在abaqus中建立的土层数要和shake91中的土层数相等，使用同样的参数，注意是弹模，泊松比，密度和瑞利阻尼。注意shake91中必须是8列或4列数据，如果使用两列数据，则无法读入加速度。Shake91是频域有限元法，对于加速度幅值是先放大，再滤波。自带的例子中滤掉大于25Hz的波。2、阻尼，时程分析中一般使用瑞利阻尼瑞利阻尼的一般公式M、K为质量和刚度矩阵，α、β分别为质量系数和刚度系数。式中，ζi和ωi分别为结构的第i阶模态的阻尼比和圆频率，关于两个频率的取值有很多取法，详见详见Parametricstudyonseismicgroundresponsebyfiniteelementmodelling详见Analysisforseismicresponseofdrystoragefacilityforspentfuel3、滤波使用seismosignal进行滤波，把高频波滤掉，shake2000手册中提到可以把大于15Hz的波都滤掉，而shake91中则把大于25Hz的滤掉，只对高频滤波，低频不滤波，因为低频也就是长周期的波对地震响应起控制作用。故在把地震波用于有限元分析前，先滤掉高频波。4、关于单元尺寸、最小增量步时间4.1单元尺寸有限单元长度必须小于最大频率对应波长的1/10-1/8，因此滤掉高频波后可以增大单元尺寸详见Parametricstudyonseismicgroundresponsebyfiniteelementmodeling4.2时步选择常选择0.005，但精确的选择如下：详见Soil–structureinteractioneffectsonse