Simulation of head-neck dynamic response in −gx and +gy