structured quadrilateral mesh in PDE Toolbox

26 ビュー (過去 30 日間)
Natalia
Natalia 2014 年 6 月 5 日
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 件のコメント
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 件のコメント
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 件)

Community Treasure Hunt

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

Start Hunting!

Translated by