Wouldn't you do it the same way? Calculate net gas force vs. time on the piston, calculate net inertia force vs time on the piston, calculate net torque vs time due to rod whipping, translate those into net torque vs. time at the crankshaft, break into half-orders of harmonic content, and apply to your mass-elastic system?