truegrid网格模型导入autodyn计算方法

truegrid网格模型导入autodyn计算方法
Truegrid和ICEM转化
如何导入truegrid的zon文件格式到autodyn中？
1.建立truegrid文件，其代码如下：
autodyn
block 1 10;1 5;-1;0 10;0 5;0; %%建立一个长10，高5的矩形，其中长度方向10个网格，高度上5个网格
merge
write
建立完毕的图形见下图；
2.找到truegrid存储文件trugrdo，复制改其格式为trugrdo.zon 文件。

3.修改数据文件
原始文件修改过的文件
其修改的主要是将第三排的数字修改，最后一排列删除。

4.导入到autodyn中
导入结束，添加材料等。

ICEM网格
1.建立ICEM文件。

在output中的output-solver中选择Multiple-info，求解器选择autodyn，
进行输出网格，在响应的文件夹中形成info.geo文件
2.打开geo文件其数据类型为图所示，可以看出他和右图的zon
格式还是有很多不一样的
地方。

Geo文件格式（按照第二列升序排列）zon格式（按照第一列升序排列）
3.修改geo文件，主要是在第一行进行一个enter，然后放入到计算求解程序中，其程序
代码为：
clc
format long e
modeldata=importdata('info.geo')%读入ICEM的geo格式文件
KK=modeldata.data
a1=KK(1,1)%模型网格起始i
a2=KK(1,2)%模型网格终止i
a3=KK(1,3)%模型网格起始j
a4=KK(1,4)%模型网格终止j
A=KK(:,1)%读取模型第一列数据
B=KK(:,2);%读取模型第二列数据
[r,c]=size(KK);
M=A(2:r,:);%读取矩阵的第二行开始的是网格节点x坐标文件
N=B(2:r,:);%读取矩阵的第二行开始的是网格节点y坐标文件
C_1=[M,N]%C_1矩阵式节点对应的坐标（x,y）
C=sortrows(C_1);%按照升序进行排序,与truegrid文件对应；
[r,c]=size(C);%r为矩阵的行数，c为矩阵的列数
s1='IMPLICIT';%要写入的字符串，truegrid开始符号
s2='END';%要写入的字符串，truegrid结束符号
s3='modeldata'%要写入的字符串，模型保存名字
fid=fopen('modeldata.zon','w');
fprintf(fid,'%s\r\n',s1);
fprintf(fid,'%s\r\n',s3);
fprintf(fid,' %g',[a1 a2 a4 a3]);%输出模型的i,j起始和终止
fprintf(fid,'\r\n');
for i=1:r
for j=1:c
fprintf(fid,' %8.8e',C(i,j));%输入模型的数据，格式为zon格式中的数据end
fprintf(fid,'\r\n');
end
fprintf(fid,'%s\r\n',s2);%end
fclose(fid);
修改之前的geo文件转变为了zon文件4.导入到autodyn中
导入过后的图形。