FEM网格玩法：四边形网格一键变为三角形网格？

在有限元分析中，网格单元是个非常重要的概念，贯穿始终。

今天给大家分享一个有趣的网格小工具，可以将四边形网格一键转换为三角形网格。该工具是由Siva Srinivas Kolukula博士编写，感兴趣的访问他的主页：https://sites.google.com/site/kolukulasivasrinivas/my-cv。

该工具是基于Matlab语言编写的，用户需提供网格节点坐标信息、单元节点构成信息，即可实现网格类型的转变，完整代码获取方式在文末。

接下来，Show time~

主函数

% QUA4 FEM mesh data
coor = load('coordinates.dat') ;
nodes = load('nodes.dat') ;
PlotMesh(coor,nodes) ;
% conver QUA4 mesh to TRI3 mesh 
tri = Qua4toTri3(nodes) ;
PlotMesh(coor,tri)

输入节点坐标信息：coordinates.dat

输入单元节点信息：nodes.dat

运行Qua4toTri3，即可实现网格类型的转变！

Qua4toTri3函数

function nodes2 = Qua4toTri3(nodes1)
if ~any(size(nodes1)==4)
    error('Input the nodal connectivity data for Qua4 elements')
end
nel1 = size(nodes1,1) ;          % total number of QUA4 elements 
%% Initilaize TRI3 data 
nel2 = 2*nel1 ;
nodes2 = zeros(nel2,3) ;
%
nodes2(1:2:end,:) = nodes1(:,1:3) ;
nodes2(2:2:end,:) = nodes1(:,[1 3 4]) ;