Hello Wei,
If lsdyna i=job.k works, why don't you continue using this command?
Which version of LS-DYNA are you using? Does the model work with the regular version of LS-DYNA (without user subroutines). If so, then the problem is likely the user subroutine.
Reno.
You are navigating away from the AIS Discovery experience