simulate_chains function

Simulate transmission chains