FLAC3D 6.0/7.0 代码转换tecplot文件

已于 2024-06-20 09:42:22 修改
阅读量148 收藏
点赞数 1
文章标签： matlab 3d 数学建模
于 2024-06-14 18:48:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wzn1129/article/details/139688250
版权

在flac3d 6.0/7.0里面打开sav文件，并运行以下代码即可在相应文件夹下生成tecplot文件，再根据往期教程在tecplot里面处理云图及等值线。

fish automatic-create on
;set fish sa off
program echo off
program	log off	
def ini_Flac2tec
    global IO_READ  = 0;
    global IO_WRITE = 1;To overwrite the exist file
    global IO_ASCII = 1;To write the file in ASCII format
    global N_RECORD = 10;MODIFY to set the figure count of each line
	global i=0;To control the loop
	global j=0;To mark the new indices of Nodes
    array my_Point(8)
    array buf(1)
	global ScaleFactor = 0;MODIFY to scale the transform
    global tec_file = 'FLAC3DMesh_deformation.dat';To define the write file name
	command;create an output range with the name writeRange
		;MODIFY the group to define the range to deal with,Default the whole file,DON'T USE group XX not
		model range create 'writeRange';group v_fgc any group v_jy any group v_rtjc any
	endcommand
end;End of ini_Flac2tec
@ini_Flac2tec
;;
;;Write Tecplot File Head
def write_FileHead
	local p_gp=gp.head
	i=0
	loop while p_gp # null
		if range.isin('writeRange',p_gp) = 1  then
			i=i+1
		endif
		p_gp = gp.next(p_gp);next grid point
	endloop
	local p_z=zone.head
	j=0
	loop while p_z # null
		if range.isin('writeRange',p_z) = 1  then
		if zone.model(p_z) # 'NULL' then
			j=j+1
		endif
		endif
		p_z = z.next(p_z);next zone
	endloop
    buf(1) =          ' \n'
    buf(1) = buf(1) + 'VARIABLES = "X" "Y" "Z" "DISP" "XDISP" "YDISP" "ZDISP" "SXX" "SYY" "SZZ" "Sig1" "Sig2" "Sig3"\n'
    buf(1) = buf(1) + 'ZONE T="FLAC3D Data" \n'
    buf(1) = buf(1) + 'N='+string(i)+','+'E='+string(j)+','+'ZONETYPE=FEBrick '
	buf(1) = buf(1) + 'DATAPACKING=block\n'+'VARLOCATION=([8-13]=CELLCENTERED)\n'
    buf(1) = buf(1) + 'DT=(SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE)'
    status = file.write(buf,1)
    ii=io.out('File Head has been written!')
end;End of write_FileHead
;;
;;Write Nodes' Coordinates X Y Z
def write_node
    ;Write x position  
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1  then
					buf(1) = buf(1)+string(gp.pos.x(p_gp)+gp.disp.x(p_gp)*ScaleFactor)+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;
	;Write y position  
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1  then
					buf(1) = buf(1)+string(gp.pos.y(p_gp)+gp.disp.y(p_gp)*ScaleFactor)+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	; 
	;Write z position  
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1  then
					buf(1) = buf(1)+string(gp.pos.z(p_gp)+gp.disp.z(p_gp)*ScaleFactor)+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
ii=io.out('Node coordinations have been written!')
end;End of write_node
;;
;;Write the data on Node
def write_nodeData
	;Write total displacement
	p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1 then
					buf(1) = buf(1)+string(math.sqrt(gp.disp.x(p_gp)^2+gp.disp.y(p_gp)^2+gp.disp.z(p_gp)^2))+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write x-displacement
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1 then
					buf(1) = buf(1)+string(gp.disp.x(p_gp))+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;
	;Write y-displcaement 
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1  then
					buf(1) = buf(1)+string(gp.disp.y(p_gp))+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;  
	;write z-displacement
    p_gp = gp.head
    loop while p_gp # null
		buf(1)=''
		n=1
		loop while n <= N_RECORD
			if p_gp # null then
				if range.isin('writeRange',p_gp) = 1  then
					buf(1) = buf(1)+string(gp.disp.z(p_gp))+ '  '
					n=n+1
				endif
				p_gp = gp.next(p_gp)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
ii=io.out('Displacement data have been written!')
end;End of write_nodeData
;;
;;Write the data in Elements
def write_zoneData
	;Write SXX
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.xx(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = zone.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write SYY
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.yy(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = zone.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write SZZ
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.zz(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = z.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write Sig1 Max Princinple
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.min(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = zone.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write Sig2 Mid Princinple
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.int(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = zone.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
	;Write Sig3 Min Princinple
    p_z = zone.head
    loop while p_z # null
		buf(1)=''
		n=1
		loop while n<=N_RECORD
			if p_z # null then
				if range.isin('writeRange',p_z) = 1 then
				if zone.model(p_z) # 'NULL' then
					buf(1) = buf(1)+string(zone.stress.max(p_z))+ '  '
					n=n+1
				endif
				endif
				p_z = zone.next(p_z)
			else
				n=N_RECORD+1
			endif
		endloop
		status = file.write(buf,1)
    endloop
ii=io.out('Element stress data have been written!')
end;End of write_zoneData
;;
;;Zone Connectivity List
;Calculate the maximum index of the array my_Grid()
def my_array
	p_gp=gp.head
	loop while p_gp # null
		ski=gp.id(p_gp)
		p_gp=gp.next(p_gp)
	endloop
end;End of my_array
@my_array
;
def write_zone
	;create an array to store the grid point data
	array my_Grid(ski) 
	;initialize the grid point array
	p_gp=gp.head
	i=1
	loop while p_gp # null
		if range.isin('writeRange',p_gp) = 1  then
			my_Grid(gp.id(p_gp))=i
			i=i+1
		endif
		p_gp=gp.next(p_gp)
	endloop
	;Write Zone Connectivity List
	p_z = zone.head
	loop while p_z # null
	if range.isin('writeRange',p_z) = 1  then
	if zone.model(p_z) # 'NULL' then
		if zone.gp.num(p_z) = 8 then;A
		;1 2 5 3 4 7 8 6
			my_Point(1)=string(my_Grid(gp.id(z.gp(p_z, 1))))
			my_Point(2)=string(my_Grid(gp.id(z.gp(p_z, 2))))
			my_Point(3)=string(my_Grid(gp.id(z.gp(p_z, 3))))
			my_Point(4)=string(my_Grid(gp.id(z.gp(p_z, 4))))
			my_Point(5)=string(my_Grid(gp.id(z.gp(p_z, 5))))
			my_Point(6)=string(my_Grid(gp.id(z.gp(p_z, 6))))
			my_Point(7)=string(my_Grid(gp.id(z.gp(p_z, 7))))
			my_Point(8)=string(my_Grid(gp.id(z.gp(p_z, 8))))
			;
			buf(1) = my_Point(1) + '  ' + my_Point(2)+ '  ' + my_Point(5)+ '  ' + my_Point(3)
			buf(1) = buf(1)+ '  ' +my_Point(4) + '  ' + my_Point(7)+ '  ' + my_Point(8)+ '  ' + my_Point(6)
		else 
		if zone.gp.num(p_z) = 4 then;B
		;1 2 3 3 4 4 4 4
			my_Point(1)=string(my_Grid(gp.id(z.gp(p_z, 1))))
			my_Point(2)=string(my_Grid(gp.id(z.gp(p_z, 2))))
			my_Point(3)=string(my_Grid(gp.id(z.gp(p_z, 3))))
			my_Point(4)=string(my_Grid(gp.id(z.gp(p_z, 4))))
			;
            	buf(1) = my_Point(1) + '  ' + my_Point(2)+ '  ' + my_Point(3)+ '  ' + my_Point(3)
	    		buf(1) = buf(1)+ '  ' +my_Point(4) + '  ' + my_Point(4)+ '  ' + my_Point(4)+ '  ' + my_Point(4)
		else  
		if zone.gp.num(p_z) = 6 then;C
		;1 2 5 3 4 4 6 6 or 1 2 4 4 3 5 6 6
			my_Point(1)=string(my_Grid(gp.id(z.gp(p_z, 1))))
			my_Point(2)=string(my_Grid(gp.id(z.gp(p_z, 2))))
			my_Point(3)=string(my_Grid(gp.id(z.gp(p_z, 3))))
			my_Point(4)=string(my_Grid(gp.id(z.gp(p_z, 4))))
			my_Point(5)=string(my_Grid(gp.id(z.gp(p_z, 5))))
			my_Point(6)=string(my_Grid(gp.id(z.gp(p_z, 6))))
			;
            	buf(1) = my_Point(1) + '  ' + my_Point(2)+ '  ' + my_Point(4)+ '  ' + my_Point(4)
	    		buf(1) = buf(1)+ '  ' +my_Point(3) + '  ' + my_Point(5)+ '  ' + my_Point(6)+ '  ' + my_Point(6)    
		else
		if zone.gp.num(p_z) = 5 then;D
		;1 2 5 3 4 4 4 4
			my_Point(1)=string(my_Grid(gp.id(z.gp(p_z, 1))))
			my_Point(2)=string(my_Grid(gp.id(z.gp(p_z, 2))))
			my_Point(3)=string(my_Grid(gp.id(z.gp(p_z, 3))))
			my_Point(4)=string(my_Grid(gp.id(z.gp(p_z, 4))))
			my_Point(5)=string(my_Grid(gp.id(z.gp(p_z, 5))))
			;
            	buf(1) = my_Point(1) + '  ' + my_Point(2)+ '  ' + my_Point(5)+ '  ' + my_Point(3)
	    		buf(1) = buf(1)+ '  ' +my_Point(4) + '  ' + my_Point(4)+ '  ' + my_Point(4)+ '  ' + my_Point(4)  
		endif;A
		endif;B
		endif;C
		endif;D
        	status = file.write(buf,1)
	endif
	endif;if inrange('writeRange',p_z) = 1 then
		p_z = zone.next(p_z);next zone
	endloop;loop while p_z # null
ii=io.out('Zone Connectivity List has been written!')
end;End of write_zone
;;
;;Main Function
def Flac2tec	
    status = file.close
	status = file.open(tec_file,IO_WRITE,IO_ASCII)
	if status = 0 then
	    write_FileHead
	    write_node
		write_nodeData
		write_zoneData
	    write_zone
	    status = file.close
	    ii = io.out('Data File for Tecplot has been created!  ' + tec_file)
	else
	    ii = io.out('Open File Error! Status = ' + string(status))			
	endif
end;End of Flac2tec	
@Flac2tec
program system("preplot FLAC3DMesh_deformation.dat")
program echo on

 

 

YH数字地球
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
FLAC3D2Tecplot501_flac3D_tecplot_
10-01
将FLAC3D计算文件的应力和位移通过导入tecplot绘制等高线。
ProCAD.2D/3D.Designer.v2008.1 1CD(最大程度地提高管道设计效率)\
software2017的博客
08-17 2124
Lumerical FDTD Solutions 2.2最佳微尺度光學模擬軟體\ MPCCI3.0.5_带加密锁.运行稳定版\ SolidCAM2008_R12_SP1（简繁体）\ ssdd 2006_USB锁.运行稳定版\ SSDD2006及STAAD CHINA 2006中文版组合(新)\ 神机妙算工程资料管理软件2008完整版，各功能好用_USB加密锁\ 珠宝设计JewelCAD
最新CAX/EDA/CFD/GIS/光学/化工/液压软件资源网
asasq112的博客
08-24 3794
 最新CAX/EDA/CFD/GIS/光学/化工/液压软件资源网 阳光软件园 所有软件资料都随时更新,急需软件可以去看看,基本上能找到你想要的! http://zhangqg.51.net e-mail: even2one@vip.sina.com;even2one@sohu.com 将以上任意链接连接起来输入IE 窗口即可进入网站 下面是一部分软件，更多软件在我们的软件列
多年收集的一些稀有软件4
热门推荐
weixin_33728708的博客
10-06 2万+
QQ:365543212　 Geovariances产品: Geovariances Isatis 2013 1CD(地质统计学软件) Geovariances.Isatis.2016.Win64 1CD Geovariances Minestis 2016 v2.0.0Win64 1CD 　 软脑公司产品: 3D-Sigma for WinALL-ISO 1CD（岩土体三维应力分...
AutoForm Master v4.04-ISO 1CD(专业板金优化分析软件，好用！)/
tingcad8的专栏
11-08 1万+
AutoForm Master v4.04-ISO 1CD(专业板金优化分析软件，好用！)/AVL_FIRE_V8.31/BENTLEY.AUTOPIPE.V08.05.01.11/BENTLEY.AUTOPLANT.STRUCTURAL.V08.06.00.11/CAESAR II 4.50-ISO/civil inroads suite 2004/ETop进销存管理系统 (V3.0网络版/Of
FLAC3D6.0转tecplot（含孔隙水压力）.zip
08-17
Tecplot系列软件是由美国Tecplot公司推出的功能强大的数据分析和可视化处理软件。它包含数值模拟和CFD结果可视化软件Tecplot 360，工程绘图软件Tecplot Focus，以及油藏数值模拟可视化分析软件Tecplot RS。
tecplot提取FLAC3D6.0和FLAC3D7.0应力应变数据代码
04-08
tecplot提取FLAC3D6.0和FLAC3D7.0应力应变数据代码
Flac3D60-to-Tecplot-e.rar
06-22
这是购买的FLac3D6.0视频的一个开发插件，flac 到tecplot的，包括以下两步：1,：把此接口dat文件放到项目文件prj一个路径里（也就是和结果文件f3sav一个路径），若是直接call的话，需要在Flac3D中先打开了f3sav...
tecplot初学者 - 副本_flac3D_tecplot_
10-04
1. **数据导入**：首先，你需要将flac3D的输出文件（通常是.FL3D格式）导入到tecplot环境中。这通常可以通过tecplot的“文件”菜单中的“导入”选项完成，选择合适的文件并指定数据格式。 2. **数据布局**：flac3D...
Transformer-LSTM分类 | Matlab实现Transformer-LSTM多特征分类预测/故障诊断
关注并私信文章链接，获取对应文章源码和数据，机器学习之心的博客。
07-10 919
Transformer-LSTM分类 | Matlab实现Transformer-LSTM多特征分类预测/故障诊断
MATLAB——运算符
MelyLenient的博客
07-10 420
关系运算符可以进行数值与数值、数值与数组、数组与数组进行比较。但数值与矩阵进行比较时，以数值为标量，将数组中每个元素与标量进行比较。MATLAB进行运算处理的顺序取决于每个运算符的优先级。在每个优先级中，运算符具有相同的优先级，并从左到右进行处理。MATLAB中算数运算符包括加、减、乘、除、点乘、点除等。除表中三项逻辑运算符外，MATLAB提供了逻辑函数。MATLAB中的运算优先级如下：由高到低排序。)，右除(./)，左除(.)，矩阵的乘法(
负荷预测 | Matlab基于Transformer-LSTM多变量时间序列多步预测
关注并私信文章链接，获取对应文章源码和数据，机器学习之心的博客。
07-14 521
负荷预测 | Matlab基于Transformer-LSTM多变量时间序列多步预测
Matlab中如何添加OptiluX？
weixin_47065524的博客
07-10 285
在Matlab中添加OptiluX
4.10、matlab生成脉冲序列：pulstran()函数
XU157303764的博客
07-10 1114
matlab生成脉冲序列：pulstran()函数
基于Python/MATLAB长时间序列遥感数据处理及在全球变化、植被物候提取、植被变绿与生态系统固碳分析、生物量估算与趋势分析应用
WangYan2022的博客
07-11 1016
定量评估植被时空动态变化是制定生态系统可持续发展目标和衡量生态系统固碳潜力的重要前提，卫星遥感数据衍生的生态参量产品为研究长时间序列全球及区域植被时空变化提供了重要数据源。
Matlab 判断直线上一点
最新发布
dayuhaitang1的博客
07-14 69
Matlab 判断直线上一点
1.27、基于径向基神经网络的曲线拟合(matlab)
XU157303764的博客
07-10 1090
径向基神经网络（Radial Basis Function Neural Network, RBFNN）是一种人工神经网络，它通过在输入空间中使用径向基函数来进行数据处理和模式识别。RBFNN通常由输入层、隐含层和输出层组成，其中隐含层使用径向基函数来处理输入数据。RBFNN的训练过程通常包括两个阶段：中心点确定和权重确定。在中心点确定阶段，通常使用聚类算法（如k-means算法）来确定隐含层中径向基函数的中心点；在权重确定阶段，通过最小化损失函数（通常是均方误差）来确定隐含层到输出层的权重。
flac3d7.0转tecplot
05-15
FLAC3D7.0是三维有限元分析软件，可以模拟地质构造和土木工程中的各种问题。而TECPLOT是一款用于可视化科学数据的软件，可以显示三维图形和动画。 要将FLAC3D7.0的结果转换为TECPLOT格式，需要进行以下步骤： 1.在FLAC3D7.0中运行模拟，并将结果输出为ASCII格式的数据文件，一般为*.dat文件。 2.使用TECPLOT软件打开*.dat文件，并选择需要显示的数据变量和场景。可以根据需要设置轴标签、网格线等参数，以及添加文本或注释。 3.保存TECPLOT图像，并可以导出为各种数据格式，如PDF、JPG或PNG等格式。也可以将TECPLOT图形导入到其他软件中进行后续处理。 总之，将FLAC3D7.0的模拟结果转换为TECPLOT格式，可以更直观地呈现设定的场景，并便于进一步分析和处理数据。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值