ABAQUS ODB文件楼层剪力和位移文本输出PYTHON代码弹塑性理论yjk 转 ABAQUS 使用体会20160909

Abaqus的Python后处理详解: /view/5cb923fc1a37f111f0855b1f.html利用Python对Abaqus进行后处理结果输出: /s/blog_6465f2ed0102x0kn.htmlhttp:/tool.lu/pyc 这个网站在线工具点按钮“选择文件”直接可免费把PYC文件转为PY文件。BAT 文件可加：call abaqus cae noGUI=singleAba_YJK.pyc这个文件安装YJK1.8后，在目录C:ProgramDatayjkSoftYJKS1.8ABAQUS下，YJK转ABAQUS目录下也有此Python文件，这个文件需要读取下图内容，最后一排为切割面的Z坐标。这个批处理能产生文件 BodyCut.dat 格式如下：*Free Body Cuts*Step label = AccelStep = 0Frame = 0Number of active free bodies = 1Name = View-CutResultant force = +8.984e-02 +4.688e-02 -6.986e+06Resultant moment = 0 0 +6.836e-02Summation point = +4.500e+00 +4.500e+00 +1.000e-01*Free Body Cuts*Step label = AccelStep = 0Frame = 1Number of active free bodies = 1Name = View-CutResultant force = +1.751e+01（VX每层剪力） +1.172e-02（VY每层剪力） -6.490e+06Resultant moment = 0 +1.800e+01（My每层总体弯矩） +3.438e-01Summation point = +4.500e+00 +4.500e+00 +1.000e-01（切割面Z坐标）根据这个文件，可以处理出楼层剪力图和弯矩图，第二步总体重量G=6490KN,比静力分析的G=6986KN,小了500kn, 纤维模型积分不精确吧，没加竖向地震加速度，但FZ一直在变，理论上FZ应该保持6986KN不变，弹塑性分析也就是骗骗专家的。由于ABAQUS显格式动力分析，不能用 *NODE print，节点位移只能从ODB文件读取根据网上资料编了一个ABAQUS 脚本文件,可以运行，也可以BAT文件加如下：call abaqus job=testaba double int cpus=1call abaqus job=testaba_xpl oldjob=testaba double int cpus=1call abaqus cae noGUI=singleAba_YJK.pyccall abaqus cae noGUI=test_NODE_DISP_OUTPUTPause在D:TEMP 产生一个TEST1.TXT文件 格式如下：121 frames length 304 0 total nodes number and & frame No. 1 -0.00000 0.00000 0F N0 2 0.00000 0.00000 0F N1 3 0.00000 -0.00000 0F N2 4 -0.00000 -0.00000 0F N3 1000001 -0.00001 -0.00001 0F N4 1000002 0.00001 -0.00001 0F N5 1000003 -0.00001 0.00001 0F N6 1000004 0.00001 0.00001 0F N7 1000005 -0.00003 -0.00001 0F N8第一个数为ABAQUS模型节点编号，第2个为X方向侧移，第3个为Y方向侧移。根据这个文件，可以处理变成层间位移比的数据，也可画顶层位移时程曲线，需扣除支座的位移。注意把上图文件改成只读，YJK转ABAQUS好像自动改成ABAQUS6141目录，导致ABAQUS 6101 动态链接库定位错而不能计算。注意上图用户材料每排只能8个数，YJK转的第9个数放入下一排。可以用ABAQUS 导入INP文件，第9个数为0,上图的ALPHA阻尼比 = （结构总质量/模型构件总自重）*混凝土大震阻尼比0.08*4*PI/第一振型周期T1，因为CLOAD附加的楼面装修荷载+50%活 的阻尼比 靠放大材料的阻尼比总体等效。上图不能输入BETA阻尼，ABAQUS不允许定义，定义了会出错警告，问了一个华工博士，说会使时间步长小，计算太慢，但SAUSAGE的界面是有BETA阻尼定义， 剪力墙结构不宜考虑BETA阻尼。_Material.inp文件中墙壳元材料的36为膨胀角，陆新征书P125126也有解释。YJK这个正确的，可以理解为：墙的最大轴压比/1.2/1.4*20100000为初始应力，陆新征的书建议初始应力为1/31/2 * fc。0.5/1.2/1.4=0.3, 初始为0.3fc是合理的，非弹性应变值 = 应变值 - 应力/Ec。红框内数据正确，非弹性应变值 = 应变值 - 应力/Ec。上图板壳元可输入BETA阻尼，无剪力墙，但输入的结果是总体地震基底剪力反而增加了，如果有剪力墙的话，或许总体会小。这种材料就是陆新征书P121 说的concrete damaged plasticity模型,36为膨胀角，陆的书p126建议为3035，6030000为混凝土初始应力0.3fck，对标准值相当于墙的轴压比/1.2的重力荷载代表值系数/1.4，fck/fc=1.4把上图的粘滞系数从0.06改大到0.08，基底剪力只小了1%。SDV6 是梁柱混凝土的受压损伤 SDV7 是梁柱混凝土的受拉损伤DAMAGEC是墙壳单元的受压损伤。YJK 转的 ABAQUS 模型 第一层叫ADD1，第2层叫ADD2，模拟施工非常清楚，不需要用*elcopy，和PKPM或YJK模拟施工3方法一样。上图YJK转ABA把梁上荷载简化为节点集中力了，所以算周期的模型梁单元可以不细分，但计算内力的模型梁单元跨度不能大于2米，否则内力计算不够精确。弹性时程分析，YJK转的模型没有BETA阻尼，加上BETA阻尼总体基底剪力减小13%。ABAQUS周期比YJK刚，梁都没考虑放大，弹性模型也没考虑钢筋，即使把楼板减薄，总质量不变，周期还是比YJK小20%，总质量两个模型是一样的，有点奇怪。只有4个柱子，X和Y刚度一样，ABAQUS为X平动，YJK为45度平动， 可能0，90度和45度惯性矩基本相等，所以YJK的解是45度。动力显格式求解，好像梁纤维单元不能用*REBAR加钢筋， YJK通过两个节点之间建一个混凝土单元和BOX钢截面单元两个单元模拟，格式如下：*ELEMENT, TYPE=B32, ELSET=C-JX-5 B32为2次插值的铁木辛科梁单元23,1000007,1000040,1 单元号，节点1，节点2，节点325,1000008,1000041,2*ELEMENT, TYPE=B32, ELSET=C-XX-724,1000007,1000040,126,1000008,1000041,2BOX 截面数据例子墙暗柱简单转为B31梁，和墙不太协调，目前没和墙象陆新征的书P133 写的通过ABAQUS的EMBED方法添加于墙体内。ABAQUS可IMPORT INP文件 查看截面属性，从上图可知0.000542为厚度，暗柱简化为箱型钢梁了。这个名义拉应力 = ( -恒-0.5活+1.0|中震|）/(B*H+(Es-Ec)/Ec*AS) 2*Ftk 广东省2没问题，外省有的专家要求1.5节点位移时程曲线越来越大，主要是支座有节点位移，如果对支座的相对位移，应该和下图类似：304 1118 total nodes number and & frame No. 1 0.31655 0.00000 1118F N0 支座有节点位移 2 0.31655 -0.00000 1118F N1 3 0.31655 0.00000 1118F N2