A direct differentiationmethod is presented for the shape design sensitivity analysis of axisymmeetric elastic solids. Based on the exisymmetric boundary integralequaiton formulation, a new boundary ntegral equatio for sensitivity analysis is derived by taking meterial derivative to the same integral identity that was used in the adjoint variable melthod. Numerical implementation is performed to show the applicaiton of the theoretical formulation. For a simple example with analytic solution, the sensitivities by present method are compared with analytic sensitivities. As an application to the shape optimization, an optimal shape of a gas turbine disc toinimize the weight under stress constraints is found by incorporating the sensitivity analysis algorithm in an optimizatio program.