This Paper discusses computing principle procedure of indirect adjustment of observation equations by point relaxation method.

Aiming to improve the calculating speed of self calibration bundle adjustment in digital industrial photogrammetry, this paper presented a speedy calculation algorithm based on point relaxation method,which could not only achieve high precision but accelerate the calculation effectively.

The application of dynamic relaxation method will complete this solution without generating the global stiffness matrix,which is difficult for traditional FEM when the topology of mode.

Based on dynamic relaxation of discrete element method(DEM),a control displacement solution for displacement integral is proposed to solve the problem caused by unreasonable time-step interval impacting on the stability of the result.

Based on three dimensional discrete element method which uses dynamic relaxation, a three dimentional method is proposed to analyse underground water and bolts, the corresponding algorithm and formulas are given, and two examples are computed.

The discrete element method based on static relaxation is extended to three dimensional case.

Based on static relaxation fast Lagarangian analysis of continue (S-FLAC ), a new method, generalized-interface-element method, is proposed to solve the problem of rock joints and interfaces.

A new num erical analysis m ethod for geotechnical engineering called static relaxation fast La- garangian analysis of continua(S-FLAC) is proposed.