Predicting strains, stresses and swelling in nuclear power plant components exposed to irradiation directly from the observed or computed defect and dislocation microstructure is a fundamental problem of fusion power plant design that has so far eluded a practical solution. We develop a model, free from parameters not accessible to direct evaluation or observation, that is able to provide estimates for irradiation-induced stresses and strains on a macroscopic scale, using information about the distribution of radiation defects produced by high-energy neutrons in the microstructure of materials. The model exploits the fact that elasticity equations involve no characteristic spatial scale, and hence admit a mathematical treatment that is an extension to that developed for the evaluation of elastic fields of defects on the nanoscale. In the analysis given below we use, as input, the radiation defect structure data derived from ab initio density functional calculations and large-scale molecular dynamics simulations of high-energy collision cascades. We show that strains, stresses and swelling can be evaluated using either integral equations, where the source function is given by the density of relaxation volumes of defects, or they can be computed from heterogeneous partial differential equations for the components of the stress tensor, where the density of body forces is proportional to the gradient of the density of relaxation volumes of defects. We perform a case study where strains and stresses are evaluated analytically and exactly, and develop a general finite element method implementation of the method, applicable to a broad range of predictive simulations of strains and stresses induced by irradiation in materials and components of any geometry in fission or fusion nuclear power plants.