The density crossover scaling of various thermodynamic properties of solutions and melts of self-avoiding and highly flexible polymer chains without chain intersections confined to strictly two dimensions is investigated by means of molecular dynamics and Monte Carlo simulations of a standard coarse-grained bead-spring model. In the semidilute regime we confirm over an order of magnitude of the monomer density rho the expected power-law scaling for the interaction energy between different chains e_intersimrho^(21/8), the total pressure Psimrho^3 and the dimensionless compressibility gT=lim(q->0)(S(q)sim1/rho^2). Various elastic contributions associated to the affine and non-affine response to an infinitesimal strain are analyzed as functions of density and sampling time. We show how the size xi(rho) of the semidilute blob may be determined experimentally from the total monomer structure factor S(q) characterizing the compressibility of the solution at a given wavevector q. We comment briefly on finite persistence length effects.