I have a model where I apply an acceleration to a point mass and I’m trying to observe the reaction force values in solid models. I used the following code for the model. However, I’m not sure if the mass element is properly receiving the acceleration load. When I change the mass within the *Mass element, the results vary, but I still think the reaction force is coming out too low. Is there an additional card I should be using in my code? Should I incorporate the gravity constant? Here’s the part of the code I wrote.

I read it from Rf node output. For an example case, reaction forces differs 10N with respect to analytical solution. I expect 50N but i got 40N. The error is about 25 percent which is very high?
Additionally, is it the true definition of applying acceleration on two direction?
*DLOAD
mass,GRAV,9800,-1.5,0.,-4.

Should i define it like in two rows for applying acceleration on both x and z direction?
mass,GRAV,9800,-1.5,0.,0.
mass,GRAV,9800,0.,0.,-4.
I try to apply 4g acceleration on z axis and at the same time apply -1.g acceleration on x axis.

Are we sure about the accuracy of Prepomax in calculating the correct vector? How does Prepomax perform this calculation? Would it be incorrect if I define acceleration in the x-direction as 1.5g and in the z-direction as 4g like this?
*DLOAD
mass,GRAV,13700, 1,0.,0.
mass,GRAV,39200,0.,0.,1.