Constraining energy-momentum-squared gravity by binary pulsar observations