@dokken Maybe you can help me out. I have created simple square geometry in Gmsh and converted the .msh file into xdmf. When i import it in fenics using your code I’m getting different results compare to the simple command "UnitSquareMesh(30,30). I’m trying tonsolve the simple linear elastic pro…

here you go to generate the square mesh: //+ Point(1) = {0, 0, 0, 1.0}; //+ Point(2) = {30, 0, 0, 1.0}; //+ Point(3) = {30, 30, 0, 1.0}; //+ Point(4) = {0, 30, 0, 1.0}; //+ Line(1) = {1, 2}; //+ Line(2) = {2, 3}; //+ Line(3) = {3, 4}; //+ Line(4) = {4, 1}; //+ Curve Loop(1) = {3, 4, 1, 2}; //+ Pl…

Please edit your posts and make sure They use 3x encapsulation surrounding codes or files, i.e. ``` //+ Point(….) ``` and ```python from dolfin import * #…. ```

here you go to generate the square mesh: //+ Point(1) = {0, 0, 0, 1.0}; //+ Point(2) = {30, 0, 0, 1.0}; //+ Point(3) = {30, 30, 0, 1.0}; //+ Point(4) = {0, 30, 0, 1.0}; //+ Line(1) = {1, 2}; //+ Line(2) = {2, 3}; //+ Line(3) = {3, 4}; //+ Line(4) = {4, 1}; //+ Curve Loop(1) = {3, 4, 1, 2}; //+ Pl…

[image] taytay33: top = CompiledSubDomain('on_boundary && near(x[1], 1, tol)', tol=1E-5) It is quite obvious as you have hardcoded the top boundary condition, which assumes that your mesh is a unit square, while what you actually have in your mesh is a [0,30]\times[0,30] mesh. Thus using t…

@dokken thanks for your kind help. I really appreciate that. A little modification. I have this image domain in which i have 2 different materials the circles in domain have different material. how can I give different values of E to different domains in fenics? The applied boundary conditions are…

One question do i have to change this xdmf file? import matplotlib.pyplot as plt from dolfin import * import meshio mesh_from_file = meshio.read("mesh.msh") def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) cell_data = mesh.get_cell_data("gmsh:physic…

[image] taytay33: cell_data = mesh.get_cell_data("gmsh:physical", cell_type) points = mesh.points[:, :2] if prune_z else mesh.points out_mesh = meshio.Mesh(points=points, cells={cell_type: cells}, cell_data={ "name_to_read": [cell_data]}) As long as you hav…

Let’s suppose the image given above is meshed using GMSH. The circles in images are extracted using OpenCv Contours detection. And I add physical tags to only boundary than whats the procedure to detect subdomains. In simple I use Gmsh for meshing the Image which has different domains how will fenic…

You can loop over the read in facets, find the connecting cells on each side and mark them. However, I would strongly advice you to marker the surfaces/subdomains (cell volumes) in Gmsh/your mesh generator, as it is alot more expensive to do recursive searches through neighbouring cells if you do n…

How to import xdmf mesh generated by gmsh into fenics

General

dokken December 29, 2022, 7:50am 12

See How to extract and use physical labels from GMSH surfaces on FEniCS (2D mesh) - #5 by dokken
and
How to define different materials / importet 3D geometry (gmsh) - #2 by dokken

Topic		Replies	Views
Assigning different physical groups to different color in gmsh and mark them as different material in fenics General	34	3512	June 6, 2023
Want to import already defined physical groups of Microstructure in Gmsh to FEnics mesh mesh	40	1554	January 13, 2024
Need help converting GMSH to FEniCS installation	61	7006	September 28, 2022
How to define different materials / importet 3D geometry (gmsh)	14	5407	November 4, 2020
Extracting mesh regions for gmsh mesh mesh	24	5391	December 5, 2023

How to import xdmf mesh generated by gmsh into fenics

Related topics