The RPIM is one of the meshless methods that by enriching and extending Base Functions in the standard form of MQ-RPIM, can be developed and solving the Fracture Mechanics problems. In this research, with using the Hermite-RBFs and optimizing the Shape Parameters in the RBFs by the Multi-Objective PSO algorithm, the Enriched-MQ-RPIM method has been converted to the OE-MQ-HRPIM method. The appropriate distribution of nodal points in the analysis domain has also been done with the Voronoi computational geometry technique. A cracked structure will be used to discretize the analysis domain and calculation of the displacement and stress fields to estimate the 3D-SIFs in Mode-I under Uni-axial Cyclic Loading. The research findings, indicate that the OE-MQ-HRPIM 's answers are closer to the experimental results compared to the MQ-RPIM, OE-MLPG Penalty and E-PIM methods. Although in this method the size of the matrices, the convergence speed and the number of unknown coefficients from n+m+l to 2n+m+l are increased, but the answers are calculated with higher density nodal point in the Sub-Local domains and fewer nodal points in the domain analysis .Also Combining Second-Order sentences in the Khayyam-Pascal pyramid (Ns=7) to defining PBFs in OE-MQ-HRPIM, will lead to more accurate answers.