You could check out PLINK2's simulation function - note that this doesn't simulate the LD patterns that you would find in real data, however - it is worth reading the original PLINK documentation on this as well

We created a set of programs called gargammel that can simulate the reads once you have the fasta file of the genome:

gargammel on github

here is the paper

For simulating demographies, I have been using msprime.