抄録
This paper presents a new algorithm to obtain principal stress trajectory grid. Grid points are calculated as the intersections of three mutually orthogonal isoplethic surfaces(two mutually orthogonal contour lines in 2D domain)of direction functions which are evaluated from the result of FEM. The contour lines cannot be obtained, when the pattern of trajectories is not topologically equivalent to an orthogonal grid. This problem with 2D pattern is solved by partitioning the triangular mesh with a line passing through the singular point where the maximum shear stress reaches the minimum value.