# تحنيب ذاتي

يمكن للعمود أن يتعرض للتحنيب بسبب ثقله الذاتي دون أي قوىً مباشرة أخرى تؤثر عليه، في وضع انهيار (فشل ميكانيكي) يسمى التحنيب الذاتي. في المسائل التقليدية لتحنيب الأعمدة، يهمل الوزن الذاتي عادةً بما أنه يفترض كونه صغيرًا بالمقارنة مع الأحمال المحورية. ولكن، عندما لا يكون هذا الافتراض صالحًا، من المهم أخذ التحنيب الذاتي بالحسبان.
أول من درس التحنيب المرن لعمود «ثقيل» (أي عمود يتعرض للتحنيب تحت تأثير ثقله الذاتي) كان غرينهيل في عام 1881. وجد أن عمودًا شاقوليًّا حرًّا، بكثافة $\rho$، وعامل يونغ مقداره $E$، ومساحة مقطعه العرضي $A$، يتعرض للتحنيب تحت ثقله الذاتي إذا تجاوز ارتفاعه قيمة حرجة معينة:
$l_{\text{max}}\approx \left(7.8373\,{\frac {EI}{\rho gA}}\right)^{\frac {1}{3}}$
حيث $g$ تسارع الجاذبية، $I$ عزم المساحة الثاني للمقطع العرضي للجائز.
اقترح غرينهيل مثالًا مثيرًا عن استخدام المعادلة في ورقته البحثية. قدر الارتفاع الأعظمي لشجرة صنوبر، ووجد أنها لا يمكن أن تنمو لطول يزيد عن 90 قدمًا. يشكل هذا الطول الارتفاع الأعظمي للشجر على الكرة الأرضية إذا افترضنا أن الأشجار موشورية وبإهمال الأفرع والأغصان.

## الاشتقاق الرياضي

نفترض عمودًا متجانسًا مثبتًا باتجاه شاقولي من نقطته السفلى، ويرتفع بمقدار $l$، وهو الحد الذي يصبح عنده التوضع الشاقولي غير مستقر ويبدأ الانحناء. هناك قوة حجمية $q$ لواحدة الطول $q=\rho gA$، حيث $A$ مساحة المقطع العرضي للعمود، و$g$ التسارع الناتج عن الجاذبية، و$\rho$ الكثافة الكتلية (الكتلة الحجمية).
ينحني العمود بمقدار بسيط تحت ثقله الذاتي، لذا يصف المنحني $w(x)$ انحراف الجائز في الاتجاه $y$ عند موضع ما $x$. بالنظر إلى أي نقطة على العمود، يمكننا كتابة علاقة موازنة العزوم:
$M=-\int _{0}^{x}q(y-w)\mathrm {d} x$
حيث يمثل الجزء اليميني من العلاقة عزم وزن BP حول P.
وفق نظرية أويلر-برنولي للجوائز:
$M=-EI{\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}$
حيث $E$ معامل يونغ لمرونة المادة، و$I$ عزم العطالة.
لذا فإن المعادلة التفاضلية للخط المركزي من BP هي:
$EI{\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}=q\int _{0}^{x}(y-w)\mathrm {d} x$
بالاشتقاق بالنسبة للإحداثي x نحصل على:
${\begin{aligned}EI{\mathrm {d} ^{3}w \over \mathrm {d} x^{3}}&=q\int _{0}^{x}\left(-{\mathrm {d} w \over \mathrm {d} x}\right)\mathrm {d} x\\&=-qx{\mathrm {d} w \over \mathrm {d} x}\end{aligned}}$
نستنتج أن المعادلة الحاكمة هي المعادلة التفاضلية الخطية من المرتبة الثالثة بمعامل متغير. تعتمد طريقة حل المسألة على استخدام متغيرات جديدة $n,z,k$ و $r$:
$k^{2}={4 \over 9}{q \over EI},\quad r^{2}=x^{3},\quad {\sqrt {x}}z={\mathrm {d} w \over \mathrm {d} x},\quad n^{2}={1 \over 3^{2}}$
ثم تتحول العلاقة إلى علاقة بيسل:
$r^{2}{\mathrm {d} ^{2}z \over \mathrm {d} r^{2}}+r{\mathrm {d} z \over \mathrm {d} r}+\left(k^{2}r^{2}-n^{2}\right)z=0$
حل المعادلة المحولة هو:
$z=AJ_{\frac {1}{3}}(kr)+BJ_{-{\frac {1}{3}}}(kr)$
حيث $J_{n}$ تابع بيسل من النوع الأول. ثم يكون حل العلاقة الأصلية:
${\mathrm {d} w \over \mathrm {d} x}={\sqrt {x}}\left(AJ_{\frac {1}{3}}\left(kx^{\frac {3}{2}}\right)+BJ_{-{\frac {1}{3}},1}\left(kx^{\frac {3}{2}}\right)\right)$
الآن، نستخدم الشروط الحدية:
* العزم معدوم عند $x=0\rightarrow {\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}(x=0)=0\rightarrow A=0$
* التثبيت عند $x=l\rightarrow w(l)=0\rightarrow J_{-{\frac {1}{3}}}\left(kl^{\frac {3}{2}}\right)=0$
من الشرط الحدي الثاني، نجد أن الطول الحرج الذي يحصل عنده تحنيب العمود الشاقولي تحت ثقله الذاتي هو:
$l_{\text{max}}=\left({\frac {j_{{\frac {1}{3}},1}}{k}}\right)^{\frac {2}{3}}=\left({\frac {9{j_{{\frac {1}{3}},1}}^{2}}{4}}{\frac {EI}{q}}\right)^{\frac {1}{3}}$
باستخدام التقريب $j_{{\frac {1}{3}},1}\approx 1.86635$، الصفر الأول من تابع بيسل من النوع الأول من الرتبة $-{\frac {1}{3}}$، يمكن تقريب $l_{\text{max}}$ إلى:
$l_{\text{max}}\approx \left(7.8373\,{\frac {EI}{\rho gA}}\right)^{\frac {1}{3}}$

## خطأ أويلر

درس أويلر العمود تحت تأثير ثقله الذاتي في ثلاث أوراق شهيرة (1778a، 1778b، 1778c). للمزيد عن القصة المدهشة والمفيدة لهذه الورقات الثلاث، يُنصح القارئ بالرجوع إلى الكتاب الذي ألفه بانوفكو وغوبانوفا (1965). استنتج أويلر في ورقته الأولى (1778a) أن العمود المستند إلى وزنه الذاتي ببساطة لن يفقد استقراره أبدًا. وصف أويلر في ورقته الثانية عن الموضوع (1778b) نتيجته السابقة بأنها متناقضة ومثيرة للريبة. وجد أويلر في الورقة الثالثة (1778c) أنه أخطأ في فهم المسألة وأن استنتاج «حمل التحنيب اللانهائي» ثبت أنه خاطئ. ولكن، ولسوء الحظ، فقد ارتكب خطًا عدديًّا وبدل أول قيمة ذاتية، حسب واحدةً أخرى. اشتق دينيك الحلول الصحيحة عام 1912، بعد 132 سنة، بالإضافة إلى ويلرز (1941)، إنغلهاردت (1954) وفريتش-فاي (1966). قدم آيزنبرغر حلولًا عدديةً بدقة عشوائية (1991). بعد 222 سنة، في عام 2000، عاد إليشاكوف إلى هذه المسألة الشهيرة واشتق حلول الصيغة المغلقة لأول مرة، لمسائل التحنيب الذاتي، عن طريق اللجوء إلى طريقة شبه المقلوب. عالج عدة مسائل أخرى في دراسته تلك.