Thin plate finite element models of the neotectonic deformation of the Australian plate have been calculated in order to estimate the stress and strain rate within the plate, specifically concentrating on the Australian continent. The model includes plate-bounding faults, an anelastic brittle-ductile layered rheology and the option of laterally varying elevation and heat flow. The results of the models are compared to (1) the velocity of geodetic benchmarks on the Australian plate, (2) the spreading rate of the mid-oceanic ridges along the Australian plate's margins, (3) the direction of the maximum horizontal principal stress, (4) the stress regime within the plate, and (5) the crustal thickness estimated from the depth to the base of Mohorovicic discontinuity's transition zone. A variety of models are tested with a wide range of input parameters. The model with the smallest misfit with observations predicts that the strain rate for most of the Australian continent is approximately 10-17 s-1. This model has a slightly lower strain rate in the central Australia and is higher off the northern coast of Australia than for the rest of the continent. Strain rates of this magnitude would be difficult to observe from geodetic or geologic data for most parts of Australia but would be enough to generate much of the seismicity that has been observed over the last century.