עקרון ד'אלמבר
קבלת משוואת אוילר לגרנז'

נבדיל בין כוחות ממשיים  {F_{i}}^{(a)} לכוחות אילוצים  f_{i}:
\sum _{i}({F_{i}}^{(a)}-{\dot {p_{i}}})\delta {r_{i}}+\sum _{i}f_{i}\delta {r_{i}}=0
כאשר האיבר השני אשר עוסק בעבודה וירטואלית אשר מתבצעת על ידי כוחות אילוצים, חייב להתאפס. לדוגמה חלקיק המוגבל לנוע על משטח- כוחות האילוצים פועלים בהכרח בניצב למישור, בעוד ההעתקים הווירטואליים חייבים להיות משיקים לו. הטיפול שלנו לא מתייחס למערכות בהן פועלים כוחות כגון חיכוך לאורך המשטח. למרות זאת במקרים מסוימים, כמו גלגול ללא החלקה, התנאי של אי קיום עבודה עדיין מתקיים.
אנו רוצים עכשיו לאפס בנפרד את המקדים של כל ההעתקות הווירטואליות, אך קיום האילוצים במערכת כופה תלות בין הקוארדינטות השונות. לכן נרצה לעבור לקואורדינטות מוכללות, בלתי תלויות אחת בשנייה.
נתאר את מעבר הקוארדינטות בעזרת מערכת משוואות:
r_{i}=r_{i}(q_{1},q_{2},...,q_{n},t)
באופן כללי מספר הקוארדינטות המוכללות שונה ממספר הקוארדינטות הממשיות (כל אילוץ מוריד את מספר הקוארדינטות המוכללות הבלתי תלויות).

ההעתקה הווירטואלית  \delta {r_{i}} קשורה להעתקה הווירטואלית  \delta {q_{i}} על ידי:
\delta {r_{i}}=\sum _{j}{\frac {\partial r_{i}}{\partial q_{j}}}\delta {q_{j}}
על ידי קשר זה נוכל להגדיר את הכוחות המוכללים  Q_{j}
\sum _{i}F_{i}\delta {r_{i}}=\sum _{i,j}F_{i}{\frac {\partial r_{i}}{\partial q_{j}}}\delta {q_{j}}=\sum _{j}Q_{j}\delta {q_{j}}
Q_{j}=\sum _{i}F_{i}{\frac {\partial r_{i}}{\partial q_{j}}}
שים לב כי היחידות הפיזיקליות של הכוחות המוכללים אינם חייבים להיות של כוח. המגבלה היחידה עליהם היא שהמכפלה שלהם עם הקוארדינטה המוכללת תהיה בעלת יחידות של עבודה.
עוד איבר במשוואה עוסק בעבודה הנעשית על ידי הכוחות המדומים. באופן דומה לכוחות המוכללים:
\sum _{i}{\dot {p}}_{i}\delta r_{i}=\sum _{i,j}m{\ddot {r_{i}}}{\frac {\partial r_{i}}{\partial q_{j}}}\delta {q_{j}}=\sum _{i,j}{\frac {d}{dt}}(m{\dot {r_{i}}}{\frac {\partial r_{i}}{\partial q_{j}}})\delta {q_{j}}-\sum _{i,j}m{\dot {r_{i}}}{\frac {d}{dt}}({\frac {\partial r_{i}}{\partial q_{j}}})\delta {q_{j}}
נוכל לבטא את המהירויות:
v_{i}={\frac {dr_{i}}{dt}}=\sum _{j}{\frac {\partial r_{i}}{\partial q_{j}}}{\dot {q_{j}}}+{\frac {\partial r_{i}}{\partial t}}
נזכור כי  r_{i} אינם תלויים באופן מפורש ב  {\dot {q}}_{j} ולכן:
{\frac {\partial v_{i}}{\partial {\dot {q_{j}}}}}={\frac {\partial r_{i}}{\partial q_{j}}}
בנוסף לכך:
{\frac {d}{dt}}{\frac {\partial r_{i}}{\partial q_{j}}}=\sum _{k}{\frac {\partial ^{2}r_{i}}{\partial q_{j}\partial q_{k}}}{\dot {q}}_{k}+{\frac {\partial ^{2}r_{i}}{\partial q_{j}\partial t}}={\frac {\partial }{\partial q_{j}}}{\frac {dr_{i}}{dt}}
נציב ביטוים אלו באיבר התנע ונפתח:
\sum _{i}{\dot {p}}_{i}\delta r_{i}=\sum _{i,j}{\frac {d}{dt}}(m{\dot {r_{i}}}{\frac {\partial r_{i}}{\partial q_{j}}})\delta {q_{j}}-\sum _{i,j}m{\dot {r_{i}}}{\frac {d}{dt}}({\frac {\partial r_{i}}{\partial q_{j}}})\delta {q_{j}}=\sum _{i,j}{\frac {d}{dt}}(mv_{i}{\frac {\partial v_{i}}{\partial {\dot {q_{j}}}}})\delta {q_{j}}-\sum _{i,j}mv_{i}{\frac {\partial v_{i}}{\partial q_{j}}}\delta {q_{j}}=\sum _{i,j}{\frac {d}{dt}}({\frac {1}{2}}m{\frac {\partial {v_{i}}^{2}}{\partial {\dot {q_{j}}}}})\delta {q_{j}}-\sum _{i,j}{\frac {1}{2}}m{\frac {\partial {v_{i}}^{2}}{\partial q_{j}}}\delta {q_{j}}=
={\frac {d}{dt}}({\frac {\partial T}{\partial {\dot {q_{j}}}}})\delta {q_{j}}-{\frac {\partial T}{\partial q_{j}}}\delta {q_{j}}
כאשר T היא האנרגיה הקינטית מבוטאת בקואורדינטות המוכללות.
אם הכוח במערכת הוא משמר, ניתן להגדיר אנרגיה פוטנציאלית סקלרית אשר מקיימת:
F_{i}=-{\frac {\partial V}{\partial r_{i}}}
לכן בקואורדינטות מוכללות איבר הכוחות נראה כך:
\sum _{i,j}-{\frac {\partial V}{\partial r_{i}}}{\frac {\partial r_{i}}{\partial q_{j}}}\delta {q_{j}}=\sum _{j}-{\frac {\partial V}{\partial q_{j}}}\delta {q_{j}}
והכוח המוכלל ניתן על ידי:
Q_{i}=-{\frac {\partial V}{\partial q_{j}}}
נסכום את האיברים:
\sum _{j}({\frac {\partial (T-V)}{\partial q_{j}}}-{\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q_{j}}}}})\delta {q_{j}}=0
כאשר  \delta {q_{j}} הם העתקים וירטואליים שרירותיים. לכן מתקיים לכל איבר בסכום:
{\frac {\partial (T-V)}{\partial q_{j}}}-{\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q_{j}}}}}=0
ברוב הבעיות האנרגיה הפוטנציאלית אינה תלויה ב {\dot {q}}_{i}.
לכן על ידי הוספת איבר שמתאפס נוכל לבטא מחדש את המשוואה כ:
{\frac {\partial (T-V)}{\partial q_{j}}}-{\frac {d}{dt}}{\frac {\partial (T-V)}{\partial {\dot {q_{j}}}}}=0
נגדיר את פונקציית הלגרנז'יאן כ
L=T-V
ובכך קיבלנו את משוואת אוילר לגרנז':
{\frac {\partial L}{\partial q_{j}}}-{\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {q_{j}}}}}=0