This paper presents the three-dimensional stress analysis of orthotropic thick plates using the three-dimensional spline strip method based on the theory of elasticity. The orthotropic plates are made of Aragonite crystal and sitka spruce. To demonstrate the convergence and accuracy of the present method, several examples are solved, and results are compared with those obtained by other exact and numerical methods based on the theory of elasticity. Good convergence and accuracy are obtained. The effects of thickness/width ratio, aspect ratio and boundary conditions on normal stress distributions of Aragonite crystal plates and sitka spruce plates are investigated. Moreover, the difference of weak orthotropic and strong orthotropic properties given to the characteristics of stress distributions are also shown.