【ABAQUS 二次开发笔记】使用keyword 、python和matlab一起处理Odb数据

news2024/9/22 1:36:57

我的主页：

技术邻：小铭的ABAQUS学习的技术邻主页
博客园 : HF_SO4的主页
哔哩哔哩：小铭的ABAQUS学习的个人空间
csdn：qgm1702

博客园文章链接：

【ABAQUS 二次开发笔记】使用keyword 、python和matlab一起处理Odb数据 - hayden_william - 博客园

用conversion shell element （S4R单元）建模层合板，有6层ply,每个lamina（ply）有3个 integration point,共计18个integration point。我想得到集合SET-Middle-elem中所有integration point的E S TSHR13 TSHR23的output。

提取出结果后，我还需要根据剪切模量计算出13 23方向的shear strain。得到10个vars的output后，利用他们计算出整个laminate的弹性应变能(将每个单元的每个ply的弹性应变能U计算出来进行累加).计算很麻烦。

方法一:python访问odb文件

通过python 脚本访问odb文件，然后多个循环进行嵌套提取出output结果。这个办法可行，能够一下子访问到每个element的每个point.

import odbA
ccessE=session.odbs['30d-50hz.odb'].steps["Step-1"].frames[10].fieldOutputs['E'].values[10]
print(E)

执行上述代码：

>>> import odbAccess
>>> E=session.odbs['30d-50hz.odb'].steps["Step-1"].frames[10].fieldOutputs['E']
>>>>E.values[10]
session.openOdb(r'G:/SIMULIA/workspace/Tanslate_repetition/DMAsimulia/30d-50hz.odb').steps['Step-1'].frames[10].fieldOutputs['E'].values[10]
>>>> print(E.values[10])
({    'baseElementType': 'S4R', 
      'conjugateData': None, 
      'conjugateDataDouble': 'unknown',     
      'data': array([0.0, -0.0, 0.0, 0.0], 'f'), 
      'dataDouble': 'unknown',     
      'elementLabel': 11,     
      'face': None,     
      'instance': 'OdbInstance object',    
       'integrationPoint': 1,   
       'inv3': 0.0,     
       'localCoordSystem': ((0.866025388240814, 0.5, 0.0), (-0.5, 0.866025388240814, 0.0), (0.0, 0.0, 0.999999940395355)),     
       'localCoordSystemDouble': 'unknown',     
       'magnitude': None,     
       'maxInPlanePrincipal': 0.0,     
       'maxPrincipal': 0.0,     
       'midPrincipal': 0.0,     
       'minInPlanePrincipal': 0.0,     
       'minPrincipal': 0.0,     
       'mises': 0.0,     
       'nodeLabel': None,     
       'outOfPlanePrincipal': 0.0,     
       'position': INTEGRATION_POINT,     
       'precision': SINGLE_PRECISION,     
       'press': -0.0,     
       'sectionPoint': 'SectionPoint object',     
       'tresca': 0.0,     
       'type': TENSOR_3D_PLANAR})
>>> len(E.values)
10800

可以看出，在用odb的话就会循环非常多次才能算完，并且不容易定位到单个element。我感觉用python的话会code会非常混乱，逻辑不容易理清，因此放弃。

方法二:edit keywords 输出到dat文件，然后用MATLAB编程计算

【Shell Element】一次获得多个积分点的输出到dat - hayden_william - 博客园

通过关键字*EL PRINT可以将多个变量在指定积分点上的结果进行输出。这样输出有个好处就是，输出的结果ABAQUS会整理为类似 table的样式。

这样的格式我可以用python进行处理，写入CSV文件。

【ABAQUS 二次开发笔记】读入TXT分析结果&输出csv文件 - hayden_william - 博客园

但是使用*EL PRINT有一个问题比较麻烦：*EL PRINT关键字语句的 first data line 一次最多16个section integration point ,多于16个integration point 则需要重复使用*EL PRINT来输出结果。而我有18个integration point 如果分两次输出的话。那输出的就是两个 table .得到结果我还需要整合成一个matrix。不方便。

但是幸好，composite laminate 的每个ply之间是应力应变满足连续性条件。我可以不输出几个重合的integration point：

Aindex=[3 5 7 9 11]

后续需要在matlab中写一个function，恢复这五个integration point 的output。

function B = expanMat(A,Aindex)
%将13个积分点数据展开为18个积分点的矩阵
    B=[];
    sizeAindex=size(Aindex);
    if isempty(Aindex)
        B=A;
        disp('A矩阵未改变')
    else
        %如果only移动一行就不循环了
        if sizeAindex(1)+sizeAindex(2)==2
            B=cat(1,A(1:Aindex(1),:),A(Aindex(1),:),A(Aindex(1)+1:end,:));
        else
            %把index向量升序排序
            Aindex=sort(Aindex);
            %得到A的行数
            rows=size(A);
            rows=rows(1);
            %loop A matrix
            for i=1:rows
                temp=A(i,:);
                % if the row operated is the member of Aindex(a array containts indexs of all rows needed to copy)
                if ismember(i,Aindex)
                    B=cat(1,B,[temp;temp]);
                else
                    B=cat(1,B,temp);
                end           
            end
        end
    end
end