Distmesh, mesh refinement

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Distmesh, mesh refinement

Mishal0488
Hi Guys

Does anyone have any experience with the use of Distmesh? I have attached a
link below.

http://persson.berkeley.edu/distmesh/

Below is my code for a mesh. I am not sure to refine the mesh to attain
better results for my finite element code. The 0.01 which is in bold allows
the mesh to get finer when compared to a value of 0.1, however anything
smaller than 0.01 does not create a better mesh.

If anyone has any experience in this package please help.

P1 = [0 0];
P2 = [0 60];
P3 = [30 30]; #load point
P4 = [60 0];
pv = [P1(1,1) P1(1,2);P2(1,1) P2(1,2);P3(1,1) P3(1,2);P4(1,1)
P4(1,2);P1(1,1) P1(1,2)];
fd = { 'l_dpolygon', [], pv };
fh = @(p) (min(p(:,3)));#ones(size(p,1),1);
[p,T] = distmesh( fd, fh, *0.01*, [-2,-2;2,2],pv);
patch( 'vertices', p, 'faces', T, 'facecolor', [.9, .9, .9] )



--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html


Reply | Threaded
Open this post in threaded view
|

Re: Distmesh, mesh refinement

Andreas Stahel
Mishal0488 wrote

> Hi Guys
>
> Does anyone have any experience with the use of Distmesh? I have attached
> a
> link below.
>
> http://persson.berkeley.edu/distmesh/
>
> Below is my code for a mesh. I am not sure to refine the mesh to attain
> better results for my finite element code. The 0.01 which is in bold
> allows
> the mesh to get finer when compared to a value of 0.1, however anything
> smaller than 0.01 does not create a better mesh.
>
> If anyone has any experience in this package please help.
>
> P1 = [0 0];
> P2 = [0 60];
> P3 = [30 30]; #load point
> P4 = [60 0];
> pv = [P1(1,1) P1(1,2);P2(1,1) P2(1,2);P3(1,1) P3(1,2);P4(1,1)
> P4(1,2);P1(1,1) P1(1,2)];
> fd = { 'l_dpolygon', [], pv };
> fh = @(p) (min(p(:,3)));#ones(size(p,1),1);
> [p,T] = distmesh( fd, fh, *0.01*, [-2,-2;2,2],pv);
> patch( 'vertices', p, 'faces', T, 'facecolor', [.9, .9, .9] )
>
>
>
> --
> Sent from:
> https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html

Dear Mishal

With distmesh I have no experience, but will have a look soon.
For my FEM code I used triangle, which woked without major problems.
If you like give FEMoctave a try. Install by
pkg install https://github.com/AndreasStahel/FEMoctave/archive/v2.0.3.tar.gz
or use github at https://github.com/AndreasStahel/FEMoctave

I hope it helps

Andreas






--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html