General Mechanical

General Mechanical

PPATH function includes the origin in every path defined

    • LucasInacioDP

      Hi all,

      I have meshed a general block volume and I'm trying to use PATH along with ESEL to select elements in the neightbourhood (of radius "R") of a given element "i". 

      I use *DO to insert the X Y Z coords of a spherical path around " i " into PPATH,

      then I issue ESEL,S,path,"MYPATH" to select elements, the expected result is a set of elements containing those around "i".

      the problem is that, although the array of coordinates XYZ does not contain the value (0,0,0), the selected set of elements includes a "tail" of elements that goes from the origin to the element "i" as if the path defined by PPATH wasn't only a sphere of points but the sphere plus the origin, and that behavior holds for any element.

      Code used:


      *AFUN, DEG

      !Dimensões do bloco

      L = 0.1 ! Lado [m]

      x1 = 0

      y1 = 0

      z1 = 0

      x2 = L

      y2 = L

      z2 = L

       BLOCK, X1, X2, Y1, Y2, Z1, Z2

      E _0= 2E11

      Poiss = 0.3

      ET, 1,SOLID185  

      MP, EX, 1, E_0

      MP, PRXY, 1,Poiss


      n_div = 30


      VMESH, 1,

      i = 15879 ! "i-th" element

      loc_ei_x = CENTRX(i)

      loc_ei_y = CENTRY(i)

      loc_ei_z = CENTRZ(i)

      PATH, Viz, 70,30, 50

      R_viz = (L/n_div)*(3**(1/2)) !Setting Neightbourhood radius


      !Filling array of spherical coords phi and theta

      ang_phi(1:21) = 0

      *DO, i, 2,21,

      ang_phi(i) = ang_phi(i-1) + 18


      ang_teta(1) = 0,120,240

      !Creating arrays to store values of the x y z coords inserted into path for verification
      *DIM,x_path, ARRAY, 21, 3

      *DIM,y_path, ARRAY, 21, 3

      *DIM,z_path, ARRAY, 21, 3

       !Inserting x y z into PPATH

      k = 1

      *DO, i, 1,21,

       *DO, j, 1,3,

      x_path(i,j) =  loc_ei_x + R_viz*cos(ang_phi(i))*sin(ang_teta(j))

      y_path(i,j) = loc_ei_y  + R_viz*sin(ang_phi(i))

      z_path(i,j) = loc_ei_z  + R_viz*cos(ang_phi(i))*cos(ang_teta(j))

      PPATH, k, , x_path(i,j), y_path(i,j), z_path(i,j) 

      k = k+1



       ESEL, S, path, Viz,



Viewing 0 reply threads
  • You must be logged in to reply to this topic.