structured quadrilateral mesh in PDE Toolbox

30 次查看（过去 30 天）
Natalia
Natalia 2014-6-5
评论： Alejandro Arrizabalaga 2017-11-29
Hello, I would like to ask if in the PDE Toolbox it is possible to generate a structured quadrilateral mesh instead of the triangular one? If so how can it be implemented in the code? My geometry is simply a rectangular of dimensions 3cm x 0.1mm and the mesh I would like to create is 0.5mm x 0.005mm. Here is a part of my code:
% Geometry description:
pderect([0 0.03 0.0001 0],'R1');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R1')
% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
setappdata(pde_fig,'MesherVersion','preR2013a');
pdetool('initmesh')
pdetool('refine')
Thank you very much in advance, N
  2 个评论
Precise Simulation
Precise Simulation 2017-11-22
编辑：Precise Simulation 2017-11-22
Although an old question, as it seems to resurface again and again, for future reference the FEA Toolbox does support structured meshes with quadrilateral and hexahedral finite elements. They can either be generated with the built-in quadmesh function, generated manually through a set of grid primitives, or imported from an external mesh generator.
Alejandro Arrizabalaga
Alejandro Arrizabalaga 2017-11-29
Yes, but the FEA Toolbox requires an additional license.

请先登录，再进行评论。

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

采纳的回答

Bill Greene
Bill Greene 2014-6-5
There is no support for quadrilateral elements in PDE Toolbox-- either in the meshers or the computational modules. However, if you have a rectangular geometry, you can generate a regular triangular mesh using the poimesh function. This option isn't available in the pdetool GUI so you would need to define your problem using the command line interface.
Bill
  3 个评论
显示 1更早的评论
Bill Greene
Bill Greene 2014-6-5
I don't know of any way to use the poimesh function with the pdetool GUI.
The only way I know to use poimesh is by defining your entire analysis as a MATLAB script using the command line functions that are part of PDE Toolbox. All of the examples on this page: http://www.mathworks.com/help/pde/examples/index.html
take that approach.
You might take a look at the example: http://www.mathworks.com/help/pde/examples/clamped-square-isotropic-plate-with-a-uniform-pressure-load.html
In this example the geometry is a square but you can easily make that a rectangle. That example creates the mesh with this function call:
[p, e, t] = initmesh(g, 'Hmax', hmax);
but you could replace that with
[p, e, t] = poimesh(g, 20);
Bill
PS: You didn't say why you needed a regular grid but if, for example, it is so you can get the results defined on a rectangular grid, you can perform the analysis with the unstructured mesh and then call the tri2grid function to interpolate those results to a grid.
Natalia
Natalia 2014-6-11
Thank you again for all the help Bill. The reason for defining a regular grid was to obtain the results already on such grid without using the unstructured grid and doing the transformation later.
In any case for one rectangle I defined it as:
width = 0.3; height = 0.0001;
gdmTrans = [3 4 0 width width 0 0 0 height height];
g = decsg(gdmTrans', 'S1', ('S1')');
[p,e,t]=poimesh(g,width/5e-3,height/0.005e-3);
figure; pdemesh(p,e,t);
pdeplot(p,e,t);
I have another question though, if I want to have a structure composed of two rectangles, one on top of another with the coordinates [0 width width 0 height height 2*height 2*height], what would be the best way to define the grid for such system as a whole?
Many thanks, Natalia

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Geometry and Mesh 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by