[Gmsh] Applying Boundary conditions using gmsh
Vaibhav Deshpande
vaibhavd85 at gmail.com
Fri Nov 13 19:07:26 CET 2009
I sent this email couple of weeks back to both gmsh at geuz.org and
gmsh-announce at geuz.org but didn't get any reply forwarding the message again
hopefully someone would get back to me on this one.
my email address is vaibhavd85 at gmail.com
Hi,
I am Vaibhav Deshpande a Graduate student at University of Wisconsin
Madison, In one of my projects I have implemented Quad4 plane stress element
and I am using gmsh as my mesher, the geometry I am using is of a cantilever
beam made of lines
1-2
2-3
3-4
4-1
where 1 2 3 4 are the coordinates of geometrical points
x y coordinates for nodes are
1, 0 0
2, 10 0
3, 10 1
4,0 1
what I want to do now is apply a pressure on line that
joins geometrical points 3 and 4 I was able to find the mesh lines
corresponding to geometrical line 3-4 using the code given below. I need to
find Jacobian of all the elements whose sides have an applied pressure. The
Jacobian is a function of coordinates of that element.
Simply put given a side I want to know which elements are made up of that
side. I can imagine that gmsh has that piece of information in some or the
other way, but I was not able to get that information as an output.
I am using following geometry file.
/*********************************************************************
*
* Gmsh cantilever beam
*
*********************************************************************/
lc = 0.5;
Point(1) = {0, 0, 0, lc};
Point(2) = {10, 0, 0, lc} ;
Point(3) = {10, 1, 0, lc} ;
Point(4) = {0, 1, 0, lc} ;
Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;
Line Loop(5) = {4,1,-2,3} ;
Plane Surface(6) = {5} ;
MyLineThree = 3;
Physical Line (MyLineThree) = {3} ;
MyLineFour = 4;
Physical Line (MyLineFour) = {4} ;
Physical Surface("Mesh") = {6} ;
the boundary conditions I want to apply are on line 4 is Dirichlet boundary
condition, which I was able to figure out,
the second boundary condition is on line 3 I want to apply a pressure on all
the elements on the top that share an edge with line 3 how should I access
this information?
Please let me know how I can go about it. Thanks a lot.
--
Regards,
Vaibhav Deshpande
--
Regards,
Vaibhav Deshpande
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20091113/e667f818/attachment.html>