Ok, then yes, I agree you can pick ##a(0) \neq 0## in your solution, and, as far as I can tell, that also makes ##\dot{a}##, ##\ddot{a}##, and ##\rho## finite at ##t = 0## (basically because you have substituted ##t + B## for ##t##, so all of the values at ##t = 0## are proportional to some power of ##B## instead of diverging).
However, this model is obviously extensible to negative values of ##t##, and when you reach ##t = - B##, your model has ##a = 0## and ##\dot{a}##, ##\ddot{a}##, and ##\rho## all infinite. So your model is not a different model from the standard one, it's just a shift of the ##t## coordinate by ##B## (strictly speaking there is a rescaling of ##t## as well). Considering the patch ##t \ge 0## in this model is simply equivalent to only considering the patch ##t \ge B## in the standard Einstein-de Sitter model. This is not a model in which the singularity theorems are violated; it's just a model in which you have artificially restricted attention to a particular patch.