This work analyzes the relationship between the shear relaxation modulus of entangled, linear and flexible homopolymer blends and its molecular weight distribution (MWD) when a fraction of the sample contains chains with molecular weight M lower than the effective critical molecular weight between entanglements Mc(eff). This effective critical parameter is defined in terms of the critical molecular weight between entanglements M-c of the bulk polymer that forms the physical network and the effective mass fraction Wc(eff) of the unentangled chains. In the terminal zone of the linear viscoelastic response, the double reptation mixing rule for blended entangled chains and a modified law for the relaxation time of chains in a polydisperse matrix are considered, where the effect of chains with M < Mc(eff) is included. Although chain reptation with contour length fluctuations and tube constraint release are still the relevant mechanisms of chain relaxation in the terminal zone when the polydispersity is high, it is found that the presence of a fraction of molecules with M < Mc(eff) modifies substantially the tube constrain release mode of chain relaxation. In this sense, a modified relaxation law for polymer chains in a polydisperse entangled melt that includes the effect of the MWD of unentangled chains is proposed. This law is validated with rheometric data of linear viscoelasticity for well-characterized polydimethylsiloxane (PDMS) blends and their MWD obtained from size exclusion chromatography. The short time response of PDMS, which involves the glassy modes of relaxation, is modeled by considering Rouse diffusion between entanglement points of chains with M > Mc(eff). This mechanism is independent from the MWD. The unentangled chains with M < Mc(eff) occluded in the polymer network also follow Rouse modes of relaxation although they exhibit dependence on the MWD.