The motion of the Earth's layers due to internal pressures is simulated in this research with an efficient mathematical model. The Earth, which revolves around its axis of rotation and is under internal pressure, will change the shape and displacement of the internal layers and tectonic plates. Applied mathematical models are based on a new approach to shell theory involving both two and three-dimensional approaches. It is the first time studying all necessary measures that increase the accuracy of the obtained results. These parameters are essential to perform a completely nonlinear analysis and consider the effects of the Earth’s rotation around its axis. Unlike most modeling of nonlinear partial differential equations in applied mechanics that only considers nonlinear effects in a particular direction, the general nonlinear terms are considered in the present study, which increases the accuracy of the amount of displacement of the Earth's inner layers. Also, the fully nonlinear and dynamic differential equations are solved by a semi-analytical polynomial method which is an innovative and efficient method. Determining the amount of critical pressure at the fault location that will cause phenomena such as earthquakes is one of the useful results that can be obtained from the mathematical modeling in this research.