How do I simulate the stress distribution that result in shrink fit of a compond cylinder.

Hi,
I model a compound cylinder, the inner cylinder has initial values before the shrink fit; inner radius, ri and outer radius (rc + delta) where delta is the interface displacement due to shrink fit. The outer cylinder has inner radius, rc and outer radius, ro.
I tried to simulate the problem using prescribed displacement using an initial value of delta on the compound cylinder boundary interface(where I place the value of u0=delta*cos(theta) and v0=delta*sin(theta) and the outer cylinder boundary. The result I obtained does not correspond with what I got from the analytic computation of the problem when I assume a value for delta.
Please, does anyone know the right approach that I can use to simulate the stresses in compound cylinder that can results from delta and the different material properties of the shrink fitted cylinders and diameters of the inner and outer cylinders when the value delta is unknown and is also dependent on the different materials of the compound cylinders.
I have attached my model. Any suggestion is greatly appreciated. Thanks

Daniel.
.