Với phương pháp giải và các bước giải được trình bày trong 3.1.1, luận án xây dựng chương trình tính toán số giải bài toán động lực học ngược cơ hệ pháo có sơ đồ thuật toán chương trình chính được trình bày trên Hình 3.1.
Hình 3.1. Thuật toán chương trình chính
Trong bước giải thứ nhất, sử dụng phương pháp Runge-Kutta, một số tài liệu đã giới thiệu thuật toán số giải hệ PTVP thuật phóng trong xác định quy luật khí thuốc trong lòng nòng và tính lực khí thuốc tác động lên khối lùi cho PPK 37mm-2N, đồng thời cũng đã đưa ra và các kết quả tính toán cho một phát bắn [1], [10], [12], [13], [15], [16]. Trong luận án này, thuật toán số tiếp tục được phát triển vào tính toán lực khí thuốc tác động lên khối lùi cho trường hợp bắn liên thanh. Sơ đồ thuật toán được trình bày trên Hình 3.2.
Quy luật chuyển động của khối lùi thân pháo cũng đã được đưa ra trong một số tài liệu nghiên cứu về máy tự động của PPK 37mm-2N. Trong luận án này, quy luật chuyển động của khối lùi hai thân pháo được tính toán có kể đến ảnh hưởng của chuyển động tầm và hướng pháo. Sơ đồ thuật toán xác định quy luật chuyển động của khối lùi được trình bày trên Hình 3.3(a).
(a) (b)
Hình 3.3. Thuật toán xác định quy luật chuyển động của khối lùi và tính mô-men truyền động
Sau khi xác định được lực khí thuốc tác động lên khối lùi và quy luật chuyển động của khối lùi hai thân pháo. Bước giải tiếp theo là tính toán mô- men truyền động. Sơ đồ thuật toán được trình bày trên hình 3.3(b).