Numerical modeling of weakly fused fiber-optic polarization beamsplitters. Part II: the three-dimensional electromagnetic model