This thesis aims to model the effects of interaction and buckling upon pairs of micro-shells embedded within an elastic medium under far field hydrostatic pressure. This analysis is motivated by the role shell buckling plays in the nonlinear nature of the pressure relative volume curve of elastomers containing micro-shells. Current models of the effective properties of these types of composites assume shells are in a dilute distribution within the host medium, and as such assume shells will buckle at the pressure of the associated isolated embedded shell model. For composites with a high volume fraction of micro-shells, or in poorly mixed composites, the dilute distribution model may provide a first approximation to the effective properties of the composite, however, interaction between shells must be considered to find a more accurate model. We begin the process of modelling the buckling of interacting embedded shells by considering the buckling of an isolated embedded thin spherical shell. For a host medium undergoing far field hydrostatic pressure we demonstrate the parameter ranges in which Jones et al. thin shell buckling theory agrees with the thin shell buckling theory of Fok and Allwright. We then use scalings to increase the range of validity of the thin shell approximation used in the Jones et al. theory to include composites with a high contrast between medium and shell materials. This enables more accurate predictions of buckling pressures of embedded shells under far field axially symmetric pressures to also be found, as is demonstrated for an embedded shell under far field axial compression. We model the linear elastic deformation of pairs of embedded micro-shells using the Boussinesq-Papkovich stress function method, before employing the thin shell linear analysis method developed in previous chapters to calculate the critical buckling pressure and buckling patterns of the pair of embedded shells.