搜索
提供常用的P50%10（50年超越概率10%），一般的工程设计地震常用这个,时间增量0.02秒。*Amplitude,name=HAMPX0.02,0.014,0.04,0.014,0.06,0.014,0.08,0.0650.1,0.014,0.12,0.016,0.14,0.219,0.16,0.130.18,0.082,0.2,0.3,0.22,0.583,0.24,0.1290.26,-0.263,0.28,-0.948,0.3,-0.105,0.32,-0.5240.34,-0.952,0.36,0.088,0.38,0.843,0.4,1.1520.42,1.716,0.44,2.523,0.46,0.07,0.48,-1.690.5,-0.708,0.52,-1.42,0.54,-1.807,0.56,-1.0910.58,-1.674,0.6,-2.547,0.62,-1.639,0.64,-2.5140.66,-5.463,0.68,-5.08,0.7,-5.128,0.72,-6.9550.74,-7.118,0.76,-5.805,0.78,-3.695,0.8,-1.8710.82,3.558,0.84,6.373,0.86,4.406,0.88,5.7690.9,10.47,0.92,11.534,0.94,10.337,0.96,12.440.98,6.454,1.,8.596,1.02,5.458,1.04,6.4031.06,0.05,1.08,1.007,1.1,-5.859,1.12,-9.4481.14,-6.851,1.16,-8.897,1.18,12.645,1.2,16.451.22,13.529,1.24,12.146,1.26,16.093,1.28,10.121.3,-9.287,1.32,24.022,1.34,22.118,1.36,21.6571.38,17.831,1.4,-1.39,1.42,10.005,1.44,8.1741.46,4.502,1.48,-2.972,1.5,-7.108,1.52,-8.635为方便大家使用，已经将其转化为标准的ABAQUS输入格式，数据文件是加速度，加速度单位是cm，请在加界中按0.01缩放！在INP中加入以下字段：;P,zC:b)V!Ce2y*AMPLITUDE,NAME=HAMPx,INPUT=X.inp*AMPLITUDE,NAME=VAMPy,INPUT=Y.inp$r.X3~d8j'L/x6c5c*AMPLITUDE,NAME=VAMP,INPUT=Z.inp----------------------------------------------------------------------------------------------8C5K(s)f-[#s6x$L1VA+k){8Q,e/]X*Boundary,op=NEW,amplitude=HAMPx,type=ACCELERATION)l0|9F$K$f2C“定义的约束集合名”,1,1,0.01（红字）场地土层反应计算采用的输入地震波是以地震危险性分析结果得到的基岩加速度峰值和基岩加速度反应谱基岩地震相关反应谱作为目标谱，用人工数值模拟方法合成得到的，并以此作为场地地震反应计算的输入地震波。;y!|4E(R6~3p+^8C按照本章6.2.1的技术思路，利用计算机自动产生的不同随机相位谱，按50年超越概率10%、5%和100年超越概率2%三个概率水准合成基岩加速度时程，其中每一个概率水准合成三条时程，分别对应于三个不同的随机相位，时程采样步长为0.02秒，目标谱与合成时程的反应谱(计算谱)之间的相对误差均小于5%。目标谱与计算谱的比较、拟合精度、基岩地震波图见图。g8J-`@-t,N(O$P#g2M;v$k9}7L!s9K离散塔身，对进口段使用三维实体单元整体建模，模型示意图见插图X-7，主要参数：顺水流方向（X轴）长30.0m，上游边界23.4m，下游边界23.4m，塔体总高86.0m，模型节点总数11426个，C3D8R（为获得时间与精度的平衡而选用）单元总数8912个。)g'A1)计算荷载和计算工况①主要考虑的荷载有：自重荷载、正常蓄水位水压力、地震加速度。②计算工况与分析时序相同，共三个工况：工程完建工况(自重荷载施加)→正常蓄水工况(静水压力施加)→正常蓄水+地震工况(动水压力施加和地震加速度施加)。2)地震波输入7C1j(R!u+D%I*\)z计算地震波输入采用地震安评报告所提供的闸址(基岩)场地谱人工合成的地震波，概率水平为50年超越概率10%，地震波时程见插图X-X。;D,k-b1T#W#w6a!l计算时，在模拟岸塔式进口施工完成及蓄水过程的基础上，假定运行期某一时刻发生地震，同时输入顺水流向、垂直水流向和竖直向地震曲线，本次计算将No.1、No.2分别作为顺水流向和垂直水流向输入，No.3作为竖向输入。并将竖直向峰值加速度调整为水平向的2/3倍，即顺水流向、垂直水流向、竖向的峰值加速度分别为100cm/s2、100cm/s2、66.7cm/s2。1}8D:P;l%w+U4)计算成果及分析0r0j.s1\9N5O.U①位移成果：3X1a.h+I3i*T9v蓄水期闸顶(相对于闸基岩)顺水流方向的水平位移为-1.1mm，蓄水+地震工况下最大水平位移达15.0mm，见插图X-8。蓄水期闸顶沉降2.9mm，蓄水+地震工况下沉4.8mm，增加沉降约1.9mm，见插图X-9。:c3w#A3u4j*D(q,~(n5Q#F3p/e4K)K*o0T（一个完整的位移时程如图所示，从小到大，再趋近于0（永久大变形除外））②应力成果：泄洪放空洞进口竣工期：拉应力主要出现在胸墙底部，最大值约0.80MPa(见插图X-10)，压应力主要分布于进口两侧边墙与底板底部，最大约2.68MPa；蓄水期：拉应力极值出现在工作闸门槽两侧边墙处，约1.07MPa(见插图X-11)，压应力主要分布于底板底部，极值出现底板末端靠山侧，约为3.71MPa；蓄水期+地震：拉应力极值出现在底板前端靠进口处，约1.28MPa(见插图X-12)，压应力极值仍然出现在底板末端靠山侧，约5.89MPa。#\9S.iP$L*s4f;T0q#V0_'R2U+K\'t;{A+v$M7R%q&v-X需要说明的是，在地震时程步中水平加速度未计入Y方向（实际意义不大）。%q!A$J0P;e)eFGC-30S4.rar(221.68KB,下载次数:118,售价:20仿真币)2011-1-1810:15:25上传下载次数:118售价:20仿真币[记录][购买]该INP中已包含马尔康加速度，时间步等可自行调整。完整计算时间约为6个小时（电脑E5200，4GB）。2u2X6T-V/}-x&p粗看了一下inp文件，觉得有以下几个方面可以予以考虑：（1）模型中未见阻尼设置选项，是否忽略了阻尼；4Q/h2p'@)T1y*t+~(O（2）最后的dynamic分析是否应该采用自动步长比较好，即增加haftol选项:j'S)~6Fa6I*Step,name=Dynamic,nlgeom=YES,inc=30008JIA*u1[/s*Dynamic0.1,30.,0.0003Suggestedinitialtimeincrement.Forimplicitintegration,thissametimeincrementwillbeusedthroughoutthestepunlesscontactimpactsorreleasesoccurortheautomatictimeincrementationschemeisused.时程步长为0.02sdynamic的参数如下是否好些：0.02,33,1e-3,0.02（3）dynamic分析中是否应该释放掉u1,u2,u3,这一点很奇怪，呵呵，不知道不释放的话算出来是什么结果（相对的？绝对的？）7R9Y4h-V/h7I!L0_a问题2关于dynamic步中约束的释放，ABAQUS自动给完成了，这里我做过一些试验可以有这个结论：为解决“震飞”现象，将模型以INP形式导入以便直观查看，导入后可以发现ABAQUS6.10自动修改在动力分析步上，自动释放了前一步的约束，从计算结果文件可以看出，动态分析步即使加约束与不加约束，所得到的结果完全一致。因此可以总结为：对于三维空间模型，选择acceleration/angularacceleration类型输入地震动时，采用CAE模型接交作业时，有同方向（同作用点）的位移约束将报错；采用INP提交时不管有无约束均能计算，其成果完全一致，可以认为程序自动忽略加速度方向的位移约束。故建议在动力分析步解除地震动方向的约束，以避免不必要的麻烦。关于地震时程分析模型位移过大问题（模型飞走），不知有没有解决了的。老问题再谈。我是一名设计院工作者，面对的是实际工程的地震时程分析之前也做过不少闸室之类的分析。!M%u+^%u0w.?7k2G!]$R当地震输入数据是加速度时，地震动的输入的前提是要释放该方向的约束，如不释放上一个step该边界的约束，该动力分析步会提示过约束的错误，导致计算中断：从导入INP后ABAQUS6.10自动修改后的文件可以看出在动力分析步上，ABAQUS自动释放了以前的约束：7f-b~:]1j%K-a+j*AMPLITUDE,NAME=HAMPx,INPUT=H.inp0G)}JW7v?*AMPLITUDE,NAME=VAMP,INPUT=V.inp(数据格式一排四组：时间,加速度时间,加速度时间,加速度时间,加速度)——————————————————————————————————**STEP:Dynamic8[E$s'L+\6XP5v*Step,name=Dynamic,nlgeom=YES,inc=3000*Dynamic3\4`7O1Q7b)Y(Q0.02,30.,0.0003'I/jA#L7p3[0\8I(?'V*Boundary,op=NEW,amplitude=HAMP,type=ACCELERATION3x;TF.OA!C,U*p2K$W!L*W_PickedSet94,1,1,0.01（这个0.01是因为我的地震数据单位是cm，将基调整为m）**_PickedSet94,2,2,0.01/}!S#q5z7f)]6H!o（为实验，专门去掉了2方向的acceleration，因此下图中为一条直线）$B0V:E0Y9@:~.a'B*Boundary,op=NEW,amplitude=VAMP,type=ACCELERATION#T9G6z9_#D-t2f(y0__PickedSet94,3,3,0.01*Boundary,op=NEW（从这里可以看出，ABAQUS自动释放了以前的约束）这里问题出现了，当地震时程很长时，模型的位移会随增量步的增加变得越来越大，从下图可以看出一个电站建在基础上的取水口直接跑了12m，这与事实是完全不相符的（即使震毁也应该是原地破坏）。一般做法是用特征点（自定）减去基础点的位移值来获得结构在地震动作用下的变形，如图2。7y4j*j:J5lQ(F)n7f;s引用shanhuimin923版主的原话：&i/c.V:R0?7U$\“对于三维空间模型，选择acceleration/angularacceleration类型输入地震动时，建议尽量解除地震动方向的约束，尽管有时不解除约束计算也有可能进行，但更多的时候程序会报错。无论是否解除该方向的约束，只要能进行计算，结构都是正确的，对结果没有影响。”我现在很想弄出在原地震动的结果，即一个点无需操做XY数据（减操做），直接输出波形文件，这样感觉更合情合理。3M3F5q;_4A8F3D*q#U-U)@!o3p8[4g下图1为特征点历时（合计值）的变位，可以看出飞了12米。下图2为特征点减底部的历时变位（看出很合理）9_4L&