Adding a Mesh Grid and Colour to Structural Cylinder

Hello

I am trying to create 3 cylinders inside each other with different mechanical properties for each cyliner.

I would ideally like to have a cuboidal mesh grid placed over the whole thing so that each cyclinder is subdivided into cubes.

I'd also like the colour of the different cylinders to be different.

How would I need to edit this code to do this:

structuralModel = createpde('structural','static-solid');
gm = multicylinder([0.01,0.04,0.05],0.05);
surface(x,y,z, 'FaceColor','texturemap',...
'EdgeColor','none','Cdata', rand(10))
structuralModel.Geometry = gm;
pdegplot(structuralModel,'CellLabels','on','FaceAlpha',0.4)
structuralProperties(structuralModel,'Cell',1,'YoungsModulus',17E9, ...
                                              'PoissonsRatio',0.3)
structuralProperties(structuralModel,'Cell',2,'YoungsModulus',0.0045, ...
                                              'PoissonsRatio',0.3)
structuralProperties(structuralModel,'Cell',3,'YoungsModulus',0.0045, ...
                                              'PoissonsRatio',0.47)

Thank you in advance!