How to export FEMesh?

Question

Fabian Günther 2021-8-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1438724-how-to-export-femesh

编辑： Wan Ji 2021-8-23

采纳的回答： Wan Ji

在 MATLAB Online 中打开

Dear Matlab Community,

how can I export a FEMesh created in Matlab?

Alternatively, it would also help me to just regularise the imported stl-file and then export it again as a stl-file.

Here is my code so far:

meshSize=0.05;
meshGrowthRate=1.1;
meshOrder='quadratic';
% import any unordered stl-file
model = createpde(1);
importGeometry(model,'geometry.stl');
% generation of a regular FE mesh
mesh=generateMesh(model,'Hmax',meshSize,'Hmin',meshSize,'Hgrad',meshGrowthRate,'GeometricOrder',meshOrder);
pdeplot3D(model)

I am very grateful for any help.

Best regards,

FG

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Wan Ji 2021-8-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1438724-how-to-export-femesh#answer_772574

在 MATLAB Online 中打开

The mesh is now the field of model

you can get mesh by

Mesh = model.Mesh; 
% nodes are obtained by
Nodes = Mesh.Nodes; % d-by-n matrix where d is the number of dimension
% elements are obtained by
Elements = Mesh.Elements; % s-by-m matrix where s is the number of nodes of one element

7 个评论
显示 5更早的评论隐藏 5更早的评论

Wan Ji 2021-8-23

在 MATLAB Online 中打开

Use the following function to write to an abaqus inp file

function toinp(Mesh, fileName, ElementType)
Mesh.Elements = Mesh.Elements';
Mesh.Nodes = Mesh.Nodes';
    if(~strcmpi(fileName(end-3:end),'.inp'))
        fileName = [fileName,'.inp'];
    end
    while(exist(fileName,'file'))
        s = input('File name already exists, are you sure to overwrite it (1 for yes and 0 for no)?');
        switch s
            case {'yes', 1}
                break;
            otherwise
                fileName = input('Please input a new file name:');
        end
    end
    
    fid = fopen(fileName,'wt');
    fprintf(fid,'*HEADING\n');
    fprintf(fid,'%s\n**\n',fileName);
    fprintf(fid,'*Node\n');
    for i = 1:1:size(Mesh.Nodes,1)
        ndim = size(Mesh.Nodes,2);
        fprintf(fid, ['%d',repmat(',%f', 1, ndim),'\n'], i, Mesh.Nodes(i,:));
    end
    switch lower(ElementType)
        case{'c2d3','c2d4','c2d6','c2d8'}
            eType = ElementType;
            eType(2:3) = 'ps';
        otherwise
            eType = ElementType;
    end
    fprintf(fid,'*Element, Type=%s\n', upper(eType));
    for i = 1:1:size(Mesh.Elements,1)
        nenode = size(Mesh.Elements,2);
        fprintf(fid,['%d',repmat(',%d', 1, nenode),'\n'], i, Mesh.Elements(i,:));
    end
    fclose(fid);
end

For example

toinp(Mesh, 'myinp.inp', 'c3d10')

Wan Ji 2021-8-23

编辑：Wan Ji 2021-8-23

在 MATLAB Online 中打开

Or just have a little modification on the function

function toinp(MeshIn, fileName, ElementType)
Mesh.Elements = MeshIn.Elements';
Mesh.Nodes = MeshIn.Nodes';
if(~strncmpi(fileName(end:-1:1),'pni.',4))
    fileName = [fileName,'.inp'];
end
while(exist(fileName,'file'))
    s = input('File name already exists, are you sure to overwrite it (1 for yes and 0 for no)?');
    switch s
        case {'yes', 1}
            break;
        otherwise
            fileName = input('Please input a new file name:','s');
    end
end
if(~strncmpi(fileName(end:-1:1),'pni.',4))
    fprintf('File name not with suffix ''.inp'', program will add it\n')
    fileName = [fileName,'.inp'];
end
fid = fopen(fileName,'wt');
fprintf(fid,'*HEADING\n');
fprintf(fid,'%s\n**\n',fileName);
fprintf(fid,'*Node\n');
for i = 1:1:size(Mesh.Nodes,1)
    ndim = size(Mesh.Nodes,2);
    fprintf(fid, ['%d',repmat(',%f', 1, ndim),'\n'], i, Mesh.Nodes(i,:));
end
switch lower(ElementType)
    case{'c2d3','c2d4','c2d6','c2d8'}
        eType = ElementType;
        eType(2:3) = 'ps';
    otherwise
        eType = ElementType;
end
fprintf(fid,'*Element, Type=%s\n', upper(eType));
for i = 1:1:size(Mesh.Elements,1)
    nenode = size(Mesh.Elements,2);
    fprintf(fid,['%d',repmat(',%d', 1, nenode),'\n'], i, Mesh.Elements(i,:));
end
fclose(fid);
fprintf('Output Inp File Successfully!\n')
end

Then use

Mesh = model.Mesh;
toinp(Mesh, 'myinp.inp', 'c3d10')

It is OK now

Fabian Günther 2021-8-23

Thanks alot man!

请先登录，再进行评论。

How to export FEMesh?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to export FEMesh?

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论