Calculation of "Broad-sense heritability" with various GRMs

[HONG-genetics]

I am trying to get dominance and epistasis effects with the custom relation matrix.
I had 5 genomic relationship matrice, A (addictive), D (dominance), AxA, AxD, DxD
and run the below test.

ALL_MODEL <- Multitrait(y=y,ETA=list(list(K=A,   model="RKHS",saveEffects=T),
                                     list(K=D,   model="RKHS",saveEffects=T),
                                     list(K=AxA, model="RKHS",saveEffects=T),
                                     list(K=AxD, model="RKHS",saveEffects=T),
                                     list(K=DxD, model="RKHS",saveEffects=T)
                                    ), verbose=T, nIter=10000,burnIn=5000,saveAt='AD_')

With this model, how can I calculate "Broad-sense heritability"?