下载文档原格式
truegrid网格模型导入autodyn计算方法Truegrid和ICEM转化如何导入truegrid的zon文件格式到autodyn中？1.建立truegrid文件，其代码如下：autodynblock 1 10;1 5;-1;0 10;0 5;0; %%建立一个长10，高5的矩形，其中长度方向10个网格，高度上5个网格mergewrite建立完毕的图形见下图；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，然后放入到计算求解程序中，其程序代码为：clcformat long emodeldata=importdata('info.geo')%读入ICEM的geo格式文件KK=modeldata.dataa1=KK(1,1)%模型网格起始ia2=KK(1,2)%模型网格终止ia3=KK(1,3)%模型网格起始ja4=KK(1,4)%模型网格终止jA=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:rfor j=1:cfprintf(fid,' %8.8e',C(i,j));%输入模型的数据，格式为zon格式中的数据endfprintf(fid,'\r\n');endfprintf(fid,'%s\r\n',s2);%endfclose(fid);修改之前的geo文件转变为了zon文件4.导入到autodyn中导入过后的图形。
三峡大学土木水电学院肖泽云
TrueGrid是美国XYZ Scientific Applications公司推出的著名、专业通用的网格划分前处理软件，支持大部分有限元分析(FEA)及计算流体动力学(CFD)软件。
它采用命令流的形式来完成整个建模过程，可以支持外部输入的IGES数据，也可以在TrueGrid中通过Block或Cylinder命令来创建基本块体，然后使用TrueGrid强大的投影功能完成各种复杂的建模。
本文主要参考TrueGrid用户手册等资料，结合作者的学习经验来具体介绍TrueGrid中常见的建模命令。
全文共分十六个部分，从最基础的启动TrueGrid开始介绍，文中凡属于命令部分都用红色字体表示。
由于作者水平有限，文中错误在所难免，敬请各位朋友批评指正！另外文中部分内容参考TrueGrid用户手册等资料，版权归原作者所有！
作者联系地址：湖北省宜昌市三峡大学土木水电学院
邮编：443002
Email: xwebsite@
QQ : 289700062
如果为右尖括号（>），则表示输入一串字符，如下图所示：。
TrueGrid和LSDYNA动力学数值计算详解1.引言T r ue Gr id是一款流行的网格生成软件，而L S-DY NA是一个常用的动力学仿真软件。
本文将详细解释T ru eG ri d和L S-DY NA的动力学数值计算方法，并介绍如何在两款软件中进行配合，以实现精确的模拟结果。
2. Tr ueGrid概述T r ue Gr id是一个用于生成各种网格的强大工具。
在进行动力学仿真前，首先需要准备一个合适的网格模型。
T r ue Gr id提供了丰富的功能和选项，可以帮助用户轻松地生成所需的网格。
该软件支持多种网格类型，包括结构化网格和非结构化网格，这使得用户能够根据模拟需求创建适应性更强的网格。
3. LS DYNA概述L S-D YN A是一个用于解决动力学问题的高性能仿真软件。
它在汽车碰撞、爆炸、结构材料测试等领域得到广泛应用。
L S-DY NA基于有限元分析方法，可以有效地模拟各种材料的力学响应和动态行为。
该软件提供了多种物理材料模型和碰撞算法，并且可以自定义应力-应变关系以满足实际需要。
4. Tr ueGrid和LSDYN A的集成通过将T ru eG ri d生成的网格导入LS-D YN A，可以实现对复杂动力学问题的精确建模和仿真。
以下是集成T rue G ri d和LS-D YN A的步骤：1.使用Tr ue Gr id生成所需的网格模型，并进行必要的后处理。
2.导出网格文件，通常为常见的文件格式如CG NS或N as tra n。
3.在L S-DY NA中导入T ru eG ri d生成的网格文件，并定义材料属性和边界条件。
4.配置仿真参数，如时间步长、停止准则等。
5.运行仿真，观察并分析结果。
5.动力学数值计算方法L S-D YN A使用显式时间积分方法来求解动力学问题。
该方法通过离散化时间和空间来近似求解动态行为。
L S-D Y NA采用C e nt ra lD if fe re nce M et ho d对时间进行离散化，并使用有限元法对空间进行离散化。
AutODana数值仿真计算AutODana数值仿真计算对于一些大型的模型，如爆炸冲击波对建筑物的破坏等。
由于爆炸过程中网格的尺寸效应较为明显，只有很小的网格采用较好的模拟爆炸初期的冲击波，（网格大的话，冲击波的压力峰值会变得很低），但是由于计算机的性能有限，所以我们可以通过先建立1维或者2维的模型，然后通过将1维或者2维的结果导入到3维模型中，方便进行结果查看。
如通过先进行2维爆炸的计算，生成map文件，然后用提取到3d模型中计算即可。
AutODana中一维结果映射到三维还有一些情况如我们需要先用欧拉域计算一些爆炸结果，如金属射流或者EFP等，等计算完成后需要提取结果对靶板进行侵彻，这个时候也可以通过结果映射，将前期计算的结果导入到模型中，完成计算。
网格填充及映射在计算爆炸问题的时候，经常有一些异形炸药的结构，其网格比较难以划分，在autodyn和LS-dyna中都提供又基于一种几何映射的网格划分方式。
对于dyna可以通过initial_volume_fractioan_geometry进行几何的映射。
在autodyn中可以通过在欧拉域通过fill by part的操作进行模型网格的转化。
LSDYNA中的DEM方法对于一些沙土问题、材料混合搅拌等，其具有离散性质，其颗粒不均匀性，在仿真模拟中较难建模。
可以使用sph方法进行建模，其计算时间会较长单元与单元之间没有强度，但是有阻尼和摩擦，这种情况下使用DEM（离散元）方法是最合适的。
Autodyn中SPH和拉格朗日耦合SPH能够较好的模拟混凝土等一些脆性材料，但是其计算时间较长，所以我们会在模拟一些问题的时候在需要考察的关键部位采用sph粒子，在其他区域采用FEM的网格即可，通过对SPH粒子和FEM网格之间定义好接触和绑定，即可完成SPH和FEM网格的耦合，即减少了计算时间，又能够得到计算所需要的效果。
烟幕初始云团最大半径数值模拟陈浩;高欣宝;李天鹏;张开创;杨洋【摘要】基于圆柱体发烟装置,采用Truegrid仿真软件建立带有\"V\"型刻槽的圆柱壳体网格模型,并用Autodyn软件建立发烟装置仿真模型,对烟幕初始云团爆炸分散过程爆轰压力变化规律进行数值模拟.借助于MATLAB提供的分段线性插值函数,将压力作为已知量引入爆炸分散理论模型,采用欧拉法对模型进行数值计算.开展了野外测试实验,采用摄像法和图像处理技术得到试验数据.结果表明,爆炸分散第一阶段发生时间为0～8.9×10-6 s,烟幕初始云团最大半径增大为原来的3～4倍;在爆炸分散第二阶段,烟幕初始云团最大半径增大到约3 m.基于Truegrid与Autodyn混合仿真,并与爆炸分散理论模型相结合的方法,相比仅利用数值积分求解爆炸分散模型的传统理论方法,烟幕初始云团最大半径的计算误差降低了3％～8％,且烟幕初始云团最大半径比单一理论模型法更接近实验结果.采用欧拉法计算模型的收敛性优于四阶龙格-库塔法.【期刊名称】《含能材料》【年(卷),期】2018(026)010【总页数】8页(P820-827)【关键词】烟幕初始云团;爆炸分散;Truegrid模拟;Autodyn模拟;运动规律;欧拉法【作者】陈浩;高欣宝;李天鹏;张开创;杨洋【作者单位】陆军工程大学石家庄校区,河北石家庄 050003;陆军工程大学石家庄校区,河北石家庄 050003;陆军工程大学石家庄校区,河北石家庄 050003;陆军工程大学石家庄校区,河北石家庄 050003;陆军北京军代局驻763厂军代室,山西太原 030000【正文语种】中文【中图分类】O383;TJ551 引言随着精确制导武器的广泛使用，发展与之对抗的干扰理论及相关技术，提高部队战场生存能力，受到世界各国军队的重视［1］。
烟幕干扰具有成本低、效费比高等优点，是干扰领域发展的重要方向。
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中
导入过后的图形。
