ABAQUS中多个inp文件的Matlab语言生成及批处理

54yiwei

ABAQUS中inp文件的非编程性质，导致了变量参数值修改的费时费工，因此，采用matlab文件读与写的形式，来生成大批量的修改参数后的inp文件，最后，采用批处理提交计算。
以下是matlab实现文本读与写的代码：
%   
%    读入ABAQUS的inp文件LB_Q110.inp，修改其中参数，并
%    生成修改参数后的适用于批处理所用文件
%    54yiwei,2010-04-20,CS,CN
%   
%
clear all
%   调用尺寸参数函数parameter()
[x]=parameter();
nfile=length(x(:,1));
%    指定要生成的文件名
Sbl='LB_Q110-';
Se='.inp';
nn=1:nfile;
for i=1:nfile
    %
    diameter=x(:,1);
    thickness=x(:,2);
    fid_inl=fopen('LB_Q110.inp','r');
    %
    if(-1==fid_inl),error(['can''t open fileB_Q110.inp']),end
    %
    Sm=num2str(nn(i));
    S_lb=strcat(Sbl,Sm);          %生成LB_Q110-XX文件名，其中XX为数字1到nfile
    Soutl=strcat(Sbl,Sm,Se);   %生成LB_Q110-XX.inp文件名
    %
    fid_outl=fopen(Soutl,'w');  %生成输出文件LB_Q110-XX.inp
    %
    %读取LB_Q110.inp中的第一行，
    nfile=1;
    Str1=fgetl(fid_inl);
    %写入到LB_Q110-XX.inp
    fprintf(fid_outl,'%s\r\n',Strl);
    %读取LB_Q110-XX.inp文件每一行，修改参数
    while ~feof(fid_inl)          %判断是否为文件尾
        %
        nline=nline+1;
        Sn=fgetl(fid_inl);         %读取第nline
        %判断修改第2、3和第4行
        if nline==2
            datestr=strcat('54yiwei,',date,',CS,CN');
            fprintf(fid_outl,'%s\r\n',datestr);
        elseif nline==3
            fprintf(fid_outl,'diameter=%5.2f\r\n',diameter);
        elseif nline==4
            fprintf(fid_outl,'thickness=%5.2f\r\n',thickness);
        else
            fprintf(fid_outl,'%s\r\n',Sn);
        end
        if(nline>300),break;end
    end
    fclose(fid_inl);
    fclose(fid_outl);
end
关于ABQUS批处理提交计算，SimWe仿真论坛有帖子介绍：
http://forum.simwe.com/viewthread.php?tid=775346&highlight=+zengguo823
inp文件批处理工具Abaqus-bat最新版1.3：
http://forum.simwe.com/viewthread.php?tid=894456&extra
后处理中结果文件.odb的自动读取问题，涉及python脚本语言，本人仍在研究中，有没有做过类似工作的朋友，能否探讨探讨。

54yiwei

后处理中结果文件.odb的自动提出问题已经解决，感谢帖子“Riks法的python后处理程序”所给的启发：
http://forum.simwe.com/viewthrea ... B%BA%F3%B4%A6%C0%ED

ckyyjs2004

感谢，学习

kensington28

程序中nline怎么和fgetl联系起来？
也就是说怎么改变读取的行数， nline并不是fgetl的参数啊

mskywalker

先学了ansys，感觉非常不适应“ABAQUS中inp文件的非编程性质”。这个应该会相当有用

kensington28

[url=http://forum.simwe.com/redirect.php?goto=findpost&pid=2106653&p
tid=929759]4#[/url] kensington28

  原文中：

%读取LB_Q110.inp中的第一行，1 w9 @; G4 u# X* e) I1 [# N
    nfile=1;

nfile应该为nline吧

很不错的资料，谢谢楼主了~~

dindindindindin

很不错的资料，谢谢楼主了~~