tensorflow打印模型结构_钢结构模型3D打印与有限元网格的融合方法
作者:魏鲁双 刘尚蔚 王 颖 魏 群
华北水利水电大学钢结构与工程研究院
中国科学院大学人工智能学院
摘 要
薄壁结构的3D打印STL文件是单侧外表面三角面网络,而钢结构模型体系的3D打印数据图形是由STL三角形沿着结构杆件、连接件各表面围成的空间三角面混合包络的集成,这与薄壁三角面网络存在较大差异。STL文件的复杂程度、单元数目都大为增加,也增加了STL三角面与有限元网格融合一体的难度。基于所建立数字图形模型,很容易得到STL的打印数据。为将STL三角面转换为满足有限元计算网格单元,需要预处理、清除、归并奇异点。通过建立结构整体节点编号和坐标体系,并依据STL三角面和法线向量,形成简化的六面体单元,进而形成整体几何矩阵,为生成有限元网格提供了有效方法,并开发了依据转换流程的计算机应用程序。通过计算实例,表明该方法为3D打印数据与有限元计算单元网格融合提供了一条可用的捷径。
由于计算机图形学软硬件的迅速发展,其处理问题的速度、规模和方式都有了极大的改善和提高,各种复杂结构都可以迅速准确地建立三维数字模型。3D打印的图形处理软件(例如CURA)可以在已建模型的各个表面形成STL格式的三角形面片包络体及相应的数据格式。可以发现,将结构离散化进行数值分析,与3D打印可属同一类的离散方式。相关研究人员曾对3D打印薄壁结构的打印STL三角面与三维有限元网格进行了研究,在基于单纯表面三角形STL网格上,提出了实用有效的薄壁结构3D打印与有限元计算统一网格的方法。对于具有一定厚度、切割打孔加工形状、焊接栓接等连接方式的钢结构体系,若采用有限元计算,其实体单元网格划分则比较困难,通常采用杆单元或简化连接节点的种类形式进行计算,单元处理的工作量相当繁琐巨大。与常规的STL转化有限元网格的思路和方法不同,3D打印结构包络面的STL数据体系和逐层增高的各平面切割数据,组合成实体的离散化计算单元,仅在数据层面上进行调整,避免了单独划分单元的复杂过程,保证了各部分满足有限单元的协调和精度要求。本文以钢结构体系的数字图形模型为例,利用融合方法进行研究,所得到的处理流程和计算结果可推广到更多的结构体系之中。
1 钢结构模型3D打印STL格式与有限元网格格式转换
1.1 3D打印STL格式
任一STL格式包括了三角形三个角点的坐标xi,yi,zi及外法线矢量ni。对于一个具有m个单元的结构而言,其STL文件将包括m段这样的循环表达式,格式如下:
(1,2,…,m个单元)
solid Object01
facet normal-7.078 227e-002 8.323 698e-001
5.496 821e-001
outer loop
vertex 3.325 273e+002 4.636 763e+002
1.153 948e+002
vertex 3.474 805e+002 2.071 548e+002
5.057 643e+002
vertex 3.803 828e+002 2.096 953e+002
5.061 541e+002
endloop
endfacet
……………………………………………………
……………………………………………………
facet normal 9.969 845e-001 7.670 629e-002
1.174 874e-002
outer loop
vertex 1.825 676e+003 2.013 525e+003
9.046 692e+001
vertex 1.872 824e+003 1.396 057e+003
1.208 904e+002
vertex 1.873 184e+003 1.396 084e+003
9.021 741e+001
endloop
endfacet
endsolid Object01
其中,facet normal为三角形的外法线方向单位矢量,关键词“outer loop”和“end loop”之间包含了三角形的三个顶点x,y,z坐标。
图1为STL格式的角点与法线示意,图2为有限元简化的六面体单元及角点示意。由图1、图2可知,对于m个三角面的结构则形成如下矩阵。
图1 STL格式的角点与法线
图2 有限元简化的六面体单元及角点
整体单元表达为:
对此矩阵给定容差值进行近似点合并、奇异点消除,形成该结构的n个整体节点列阵,其中n、m符合欧拉关系数。
在整体节点数表示中,m个三角形单元数列表示为9个x,y,z坐标转换为总节点[P]中相应的3个节点编号。
n个整体节点编号为:
m个单元表达为:
式中:a,b,c,d,e,f,i,j,k,n1,n2,n3为中元素,形成3D打印STL外表面三角网的矩阵形式。
1.2 有限元格式
对于钢结构模型要进行有限元力学计算分析,即要先将钢桁架3D打印STL外表面三角网的数据格式转换为用有限元方法表示的数据格式,3D打印STL表面三角网同m个单元数的表达。
1.3 三维有限元六面体单元的形成
由外表面三角生成三维六面体单元,在STL格式下,八结点六面体单元实体简化为六角点五面体单元,结构体由一系列空间三角形平面表示。产生三维有限元单元需要对外表面三角形网格进行归类调整,消除奇异点,如三角形角度太小或太大,两结点距离太小等。通过自主开发程序对其进行处理后,形成符合计算要求的三维空间三角形网格,以满足计算模式之需要。图3为描述于上述三步骤的处理流程;图4为处理过的钢结构外表面空间三角形网格。
图3 外表面三角网处理流程
图4 处理后的3D打印外表面三角形网格
1.4 形成三维有限元网格
本文仅针对杆件横断面为空心薄壁形状的钢桁架为例,将其外表面三角形的每一个结点沿着该结点连接处所有单元的外法线矢量合成的负方向拉伸出新的结点,形成内表面三角形,如图5所示。
图5 形成三角形内表面示意
以四个三角形(在钢桁架角点处会有更多的三角形)外表面为例,研究其生成内表面点的方法。四个三角形的外法线(归一化)方向分别为其中,余类推,它们的合成矢量为:
(1)
式中:分别为x、y、z方向的单位矢量。
矢量的模为:
(2)
将矢量归一化为:
(3)
沿着的负方向产生新的内表面结点,设外表面与内表面的距离为D,在负方向的截距为DL。D与DL的夹角α近似用与的夹角代替。
和为归一化矢量,有:
(4)
(5)
(6)
产生的新内表面结点坐标为(在的负方向取长度DL):
x′=x-DL×ax (7a)
y′=y-DL×ay (7b)
z′=z-DL×az (7c)
每一个上表面结点都会按上述方法在下表面产生新的结点,并归入已有的总节点标号系列之中,根据上表面单元信息,可得到体单元信息。
2 3D打印与计算实例
与薄壁结构从STL文件转化为具有一定厚度的结构体数值计算网格不同,对钢结构体系数据建模产生的STL格式按照上节所述步骤方法进行处理,得到3D打印STL与有限元计算的统一网格。
2.1 钢桁架3D打印
笋溪河特大桥的三维钢桁架某一跨间的实体结构如图6所示。
图6 钢桁架实体结构
该钢桁架是一对称结构,将其输出为STL格式的文件,经过专用的切片软件“CURA”,将模型在高度方向切片。图7为某一高度切片图形。
图7 钢桁架某一高度切片
钢桁架模型尺寸为400 mm×254 mm×86 mm,每隔0.20 mm切一层,共430层,形成3D打印机打印文件“钢桁架.gcode”,并通过3D打印机打印。
2.2 钢桁架三维有限元建模及奇异点修正
单跨钢桁架杆件较少,STL文件长度约1 500行,通过“ZHSTL程序”将钢桁架的STL文件转换为有限元数据文件格式,形成具有一定厚度的三角形单元。结点数1 918个,单元数2 264个。相对于艺术花瓶的结点数和单元数少得多,由于桁架的杆件主要由平面组成,平面的数量较少,从而形成的三角形网格也很少。但同时也带来了不少问题,例如有许多网格的三角形非常细长,产生很尖或很钝的角,也有一定扭曲的或奇异的三角形单元,需要经过一定的处理,才指清除奇异点。
将有限元数据导入ADINA大型结构分析软件中,形成计算模型,如图8所示。
图8 钢桁架有限元模型
钢桁架的打印材料为聚乳酸(PLA),其力学参数如表1所示。
表1 PLA材料的力学参数
边界条件:钢桁架底部四个角固定,即约束底部四个角x=0、y=0、z=0。
计算工况:1)钢桁架在自重荷载下变形及应力状况;2)钢桁架在自重荷载、顶部承受总合力为10 N的荷载(平均分布于钢桁架顶部)下变形及应力状况。
2.3 钢桁架三维有限元计算结果
利用3D打印STL转换为有限元六面体网格进行了各工况计算。
图9、图10分别为钢桁架自重荷载作用下x向位移云图和最大主应力(σ1)云图。图11为钢桁架在自重荷载、顶部承受总合力为10 N荷载(平均分布于钢桁架顶部)下y向位移云图和最小主应力(σ3)云图。由图9~图11可知,计算结果合理,满足了工程需要。
图9 桁架自重x向位移云图 m
图10 桁架自重最大主应力(σ1)云图 MPa
图11 桁架加载y向位移云图 m
图12 桁架加载最小主应力(σ3)云图 MPa
3 结束语
在作者团队提出的数字图形介质理论当中,已经扩充并更新了数据采集的一些新技术,例如三维激光扫描方法、全景摄影方法、对视几何拍摄方法、无人机倾斜式摄影方法、室内全景摄影方法,它们都可以很容易得到结构的点云数据,据此得到的逆向数字模型可自动转化为STL或OBJ格式的文件,也就是包络整个工程的Polygon Soap,通过单体法分离,各个结构件可以进行三维单体打印。
本文基于STL格式,通过融合转化为数值计算网格,可以迅速进行数值计算,在测算与快速评估工程应力、应变、变形分布规律的过程中,提出了高效率的图形网格转换方法,具有较强的实用性与普遍性。随着研究的继续深入,可以预计该种方法一定会有更多的应用和提高。
来源:魏鲁双,刘尚蔚,王颖,等. 钢结构模型3D打印与有限元网格的融合方法[J]. 钢结构, 2019, 34(12): 88-92.
DOI: 10.13206/j.gjg201912017
点击
阅读原文
下载全文
推荐阅读:
文章精选 | 周绪红院士关于冷弯薄壁型钢(CFS)结构房屋关键技术及其应用的研究进展
装配式钢结构梁柱连接节点研究进展
论文推荐 | 3D打印技术在工程建设中的应用及前景
该篇文章刊载于《钢结构》2019年第12期,如需购买请点击以下链接:
tensorflow打印模型结构_钢结构模型3D打印与有限元网格的融合方法相关推荐
- pla3d打印材料密度_求解:3D打印 PLA 材料究竟会不会被水分解?
几天前,Reddit上一个ID为PuffThePed的网友在Reddit上发帖称,2012年,他用白色PLA线材3D打印了一个小玩意儿,为了给它上色,就将其泡在了一瓶水中,但是很快就把这件事忘得一干二 ...
- 华为手机楷体字体下载_手机查看3D打印模型图纸,国产软件SView来了
南极熊导读:3D打印模型,往往都需要使用电脑上的专业3D软件来打开和处理.现在,国内拥有三维CAD内核技术的山东华天软件,已经开发出手机版3D打印模型图纸查看软件.可有效地释放3D打印的处理空间场景, ...
- 点坐标(点云)表示的模型转为.stl文件用于3D打印
点坐标(点云)表示的模型转为.stl文件用于3D打印) 步骤 这里的模型主要针对优化结束后获得的二维模型.由于模型的不规则,最终得到的优化结果是模型上点的坐标及这个点上的密度值.这个方法可以拓展到任何 ...
- python 3d打印_基于Python的结构拓扑优化与3D打印试验研究
收稿日期:2017 -05 -22 修回日期:2017 -06 -01 第 35 卷 第 8 期 计 算 机 仿 真 2018 年 8 月 文章编号:1006 -9348( 2018) 08 -017 ...
- 不需要支撑结构的3D打印技术:多轴机器人3D打印系统
为期五天的世界上最大的计算机图形和交互技术会议-SIGGRAPH 2018将于下周在温哥华启动.来自世界各地的许多研究团队在会上展示他们的最新项目成果,包括3D形状.增强现实.拓扑优化和3D扫描.今年 ...
- tensorflow打印模型结构_五行代码用图提升模型表现,TensorFlow开源NSL神经结构学习框架...
选自Medium 机器之心编译 参与:一鸣.李亚洲 今日,谷歌 TensorFlow 宣布推出神经结构学习(NSL)开源框架,它使用神经图学习方法来训练带有图和结构化数据的神经网络. 据谷歌 Tens ...
- 安卓开发3d模型展示源码_玩具模型样品3D打印,让玩具开发无限可能
3D打印技术的进步,正在悄无声息的影响着各个行业与领域.对各大玩具爱好者来说,这是一个激发创造力以及将自己想象力变为现实的好东西.据统计,越来越多的玩具厂家用3D打印机来制作玩具样品,以便开发出更好的 ...
- java识别3D打印图_如何识别3D打印FDM机器结构i3、并联臂、UM、MB
导读:FDM打印机牌子辣么多,看得人眼花缭乱,但其实万变不离其宗,按照结构来分也就那么几种,i 3.并联臂.UM.MB等类型.下面3D海豚带大家简单分析下四种结构打印机的优缺点,选选看哪种机器适合你. ...
- 决策树 结构_决策树模型是什么?5步了解好看模型图
在复杂的决策情况中,企业完成一个决策后,后面可能面临n种可能状态的发生,而决策树模型是基础的数学定律,它通过已知信息,通过逻辑推理,将问题中的策略,概率,风险,收益用类似于树状的形式呈现.决策出各个可 ...
最新文章
- ERP实施完了,为什么还要做MES?
- iOS 9应用开发教程之iOS 9新特性
- 【转】解决wine中文乱码的问题
- 华为OJ: 公共字符串计算
- 目前银行利息高还是余额宝利息高?如何选择?
- ios点击有300毫秒延迟,输入框必须重压或长按才能获取焦点唤起软键盘
- 【ARTS】01_12_左耳听风-20190128~20190203
- 三步彻底清除电脑的mysql
- 吴恩达《机器学习》第十七章:大规模机器学习
- Paper pass使用方法总结,毕业论文查重攻略
- suse 11添加阿里源
- 手机wifi已连接但无法访问互联网_我们的手机WiFi出现“已连接但无法上网”时咋办?...
- FFmpeg —— ffplay源码 - 制作桌面动态壁纸
- ICPC North America Qualifier 2017 B.Bumped! (分层图 + spfa)
- mysql怎么设置不要密码_mysql无密码怎么设置密码
- 表格的下拉放大 ----------王朋
- WebApi编程(一)-DOM
- android 支持最低版本是多少g,手机需要多少G的运行内存,才真正够用?
- Unix 时间戳 (Unix Timestamp) 与 Windows 时间转换工具
- Python光的干涉仿真
热门文章
- 【Python】汉诺塔问题
- 【Python】画一个心形
- C#LeetCode刷题之#700-二叉搜索树中的搜索(Search in a Binary Search Tree)
- 通用职责分配软件原则之2-创造者原则
- React Router教程–如何使用代码示例渲染,重定向,切换,链接等
- redux reducer_如何通过约定创建Redux reducer
- vimb java_vim 中更好的編輯 java 文件.
- Python自动化办公知识点整理汇总
- 教你用Python 每日定时推送睡前小故事给你__的人
- 程序员如何 10 分钟用 Python 画出蒙娜丽莎?