I sent this email couple of weeks back to both <a href="mailto:gmsh@geuz.org">gmsh@geuz.org</a> and <a href="mailto:gmsh-announce@geuz.org">gmsh-announce@geuz.org</a> but didn't get any reply forwarding the message again hopefully someone would get back to me on this one.<br>
<br>my email address is <a href="mailto:vaibhavd85@gmail.com">vaibhavd85@gmail.com</a><br><br><div class="gmail_quote"><br><br>Hi,<div><br></div><div>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</div>


<div><br></div>
<div>1-2</div><div>2-3</div><div>3-4</div><div>4-1</div><div><br></div><div>where 1 2 3 4 are the coordinates of geometrical points </div><div>x y coordinates for nodes are</div><div>1, 0 0</div><div>2, 10 0</div><div>3, 10 1</div>


<div>4,0 1</div><div><br></div><div>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.</div>


<div><br></div><div>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.</div>



<div><br></div><div>I am using following geometry file. <br><br>/********************************************************************* <br> *<br> *  Gmsh cantilever beam<br> *<br> *********************************************************************/<br>

<br>lc = 0.5;<br><br>Point(1) = {0, 0, 0, lc};<br>Point(2) = {10, 0,  0, lc} ;<br>Point(3) = {10, 1, 0, lc} ;<br>Point(4) = {0,  1, 0, lc} ;<br><br>Line(1) = {1,2} ;<br>Line(2) = {3,2} ;<br>Line(3) = {3,4} ;<br>Line(4) = {4,1} ;<br>

<br>Line Loop(5) = {4,1,-2,3} ;<br><br>Plane Surface(6) = {5} ;<br><br>MyLineThree = 3;<br>Physical Line (MyLineThree) = {3} ;<br><br>MyLineFour = 4;<br>Physical Line (MyLineFour) = {4} ;<br><br>Physical Surface("Mesh") = {6} ;<br>

<br>the boundary conditions I want to apply are on line 4 is Dirichlet boundary condition, which I was able to figure out,<br><br>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?<br>

<br>Please let me know how I can go about it. Thanks a lot.</div><div><br></div><div>-- <br>Regards,<br>Vaibhav Deshpande<br>

</div>
</div><br><br clear="all"><br>-- <br>Regards,<br>Vaibhav Deshpande<br>