# Newton's law of universal gravitation **URL:** \ **Category:** Uncategorized\ **Created:** [August 14, 2023, 8:35am UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768 "2023-08-14T08:35:13Z")\ **Posts on this page:** 6\ **Page:** 1
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [August 14, 2023, 8:35am UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/1 "2023-08-14T08:35:13Z")
Hello everybody. Is there any way to apply the gravitational force on a body in its most general expression? ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/9/97cbdef715241536d3dbe29e2cf73b49270ccade.png) I would impose m1. m2 is my model which has a finite dimension. Regards ---
**Author:** ![vicmw](https://avatars.discourse-cdn.com/v4/letter/v/6de8d8/32.png) [@vicmw](https://calculix.discourse.group/u/vicmw)\ **Post date:** [August 15, 2023, 8:06am UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/2 "2023-08-15T08:06:48Z")
There’s the NEWTON option for \*DLOAD. I always wondered who would use that but unsurprisingly it’s you 😉 It looks like you could do orbital mechanics with this and dynamic. ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [August 15, 2023, 9:23am UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/3 "2023-08-15T09:23:41Z")
Oh!! How did I miss that? I looked up NEWTON in the manual but only saw references to the Newton-Raphson numerical method and Newton’s second law. It really didn’t make sense that the gravitational constant would appear under \*PHYSICAL CONSTANTS and then not be used. This is good news because the Coulomb force have exactly the same shape. There is an analogy to the electric field just converting the different constants. Thank you very much Victor !! ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/126cd3a7a5679512f02be8520236e14399d6214f.png) ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [August 15, 2023, 1:16pm UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/4 "2023-08-15T13:16:44Z")
Hi Victor, I think your help deserves some promo footage for Mecway 😉 and Calculix. This simple test runs in 2:30 min. It’s linear dynamics with arbitrary initial velocities in the three components to check convergence. Enough to understand why so few planets survived in the solar system. Only those areas at the right distance, density and initial velocities are able to find an stable orbit one around each other. The impacts could probably be simulated moving to Nonlinear dynamics + contact. I have no doubt that any student or fan of celestial dynamics will be happy about this possibility in CalculiX. NOTE: Some distance and time units scaling are needed to avoid rounding errors. Keep in mind we are managing very large dimensions with gravitational constant O(1Ee-11) , millions of Km and days/years time periods ![Newton-Gravitational-force](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/c/c7d1e738ab84f63b93e0162e7ec9e6bde3a9eb43.gif) ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [August 15, 2023, 4:31pm UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/5 "2023-08-15T16:31:18Z")
It is a bit absurd to answer myself, but I do not want to raise false expectations. Applying the analogy to an electrostatic calculation has not been so simple. First because the charge is not uniformly distributed over the surface as density is in the volume and second because, even if only as a first approximation, Calculix does not admit neither densities nor a negative gravitational constant to be able to model the electrostatic repulsion of charges of the same sign. Anyway, the gravitational part is not bad at all. ---
**Author:** ![vicmw](https://avatars.discourse-cdn.com/v4/letter/v/6de8d8/32.png) [@vicmw](https://calculix.discourse.group/u/vicmw)\ **Post date:** [September 6, 2023, 10:56am UTC](https://calculix.discourse.group/t/newtons-law-of-universal-gravitation/1768/6 "2023-09-06T10:56:50Z")
Beautiful animation @Disla !