Monotonicity-constrained Bézier surrogates for thermoelastic multiscale optimization of heterogeneous lattices