الحساب الهرمي في الهندسة المدنية | 20 تطبيقًا واقعيًا مع أمثلة عددية

الملخص (Abstract)

تقدم هذه الصفحة نسخة Paper أكاديمية لتطبيق الحساب الهرمي في الهندسة المدنية. كثير من الإخفاقات الهندسية ليست نتيجة “زيادة كبيرة في القيم”، بل نتيجة تسارع في التدهور النسبي. يصف \(D_0^1\) الزيادات المحلية، بينما يصف \(D_1^1\) التدهور النسبي والمؤشرات المعيارية مثل الاستغلال (utilization)، أما \(D_2^1\) فيمثل طبقة إنذار مبكر عندما يتغير معدل التغير النسبي بسرعة (نقطة تحول). نقدم 20 حالة واقعية مع أمثلة عددية ومعادلات هندسية إضافية ونقاش موسع لدعم اتخاذ القرار في الرصد الإنشائي والجيوتكنيك والمتانة والنمذجة العددية وإدارة الأصول.

Keywords: الحساب الهرمي؛ الرصد الإنشائي؛ المتانة؛ الهبوط؛ استقرار المنحدرات؛ التعب؛ الاستغلال؛ المخاطر غير الخطية؛ تقارب FEM؛ اتخاذ القرار.

فهرس المحتوى

1) المقدمة
2) المنهجية: \(D_0^1\), \(D_1^1\), \(D_2^1\)
3) جدول تفسير الرتب
4) التطبيقات والأمثلة (1–20)
5) دليل عملي لاختيار الرتبة
6) نقاش عام + قيود + عمل مستقبلي
المراجع

1) المقدمة

تتغير الأنظمة المدنية عبر مقاييس شديدة الاتساع: من نمو شروخ بالمليمتر إلى بنى تحتية واسعة بالكيلومترات، ومن اهتزازات سريعة إلى تدهور تدريجي يمتد لعقود. لذا فإن الاعتماد على التغير المطلق وحده قد يخفي بنية التدهور الحقيقي. في الواقع، قرارات التصميم والصيانة تعتمد كثيرًا على: (1) مؤشرات معيارية (Condition Index, Utilization, Damage Index)، (2) معدلات تدهور نسبية، (3) تسارع هذه المعدلات قرب حدود الفشل.

ينظم الحساب الهرمي هذه الفكرة ضمن سلم للتغير: مطلق (\(D_0^1\))، نسبي (\(D_1^1\))، وبنيوي/لوغاريتمي (\(D_2^1\))، ما يقدم لغة موحدة لقراءة بيانات الرصد وتحولات المخاطر بين الأنظمة.

2) المنهجية: المشتقات الهرمية

نعتمد ثلاث طبقات (درجة 1):

\[ D_0^1 f(x)=\frac{df}{dx}, \qquad D_1^1 f(x)=\frac{d\ln f}{d\ln x}, \qquad D_2^1 f(x)=\frac{d\ln(\ln f)}{d\ln(\ln x)}. \]

تفسير هندسي:
\(D_0^1\): زيادات محلية (mm, kN, points).
\(D_1^1\): تدهور نسبي/مؤشرات معيارية قابلة للمقارنة.
\(D_2^1\): تسارع في التغير النسبي (إنذار مبكر لنقطة تحول).

ملاحظة المجال: الصيغة الصارمة لـ \(D_2^1\) تتطلب مجالًا مناسبًا. عمليًا نستخدم تحويلات لمؤشرات موجبة (Health/Risk indices, utilization) بحيث يبقى المعنى: “تغير معدل التغير النسبي”.

3) كيف نقرأ الأمثلة؟

الرتبة	الصيغة التقريبية	المعنى الهندسي	القراءة العددية
0	\(D_0^1 y \approx \Delta y\)	تغير مطلق	زيادة مباشرة
1	\(D_1^1 y \approx \Delta y/y\)	تغير نسبي	% تدهور / استغلال / مقارنة عادلة
2	\(D_2^1 y\) (تصوري)	تسارع النسَب	نقطة تحول / قرب الفشل

4) التطبيقات والأمثلة (1–20)

طريقة القراءة: اقرأ \(D_0^1\) ثم \(D_1^1\). عندما يتسارع التغير النسبي نفسه، تظهر أهمية \(D_2^1\).

الحالة 1 — نمو طول الشرخ في الخرسانة

المجال: متانة الخرسانة الهدف: إنذار مبكر الرتبة: \(D_1^1\)

ينمو طول الشرخ من 20 مم إلى 22 مم خلال دورة رصد. الزيادة المطلقة تبدو صغيرة، لكن الزيادة النسبية (10%) تحمل دلالة أكبر في ميكانيكا الضرر.

\[ D_0^1:\;\Delta L = 2\text{ mm},\qquad D_1^1:\;\frac{2}{20}=10\%. \]

معادلة إضافية (مؤشر شدة الإجهاد):

\[ K_I = Y\,\sigma \sqrt{\pi a}. \]

Discussion: نمو الشرخ ليس خطيًا مع \(a\). لذلك فإن الزيادة النسبية تعكس اقتراب النظام من نطاق تصبح فيه \(K_I\) أكبر بكثير. هذا ما يجعل \(D_1^1\) أكثر تعبيرًا من الزيادة المطلقة وحدها.

الحالة 2 — تسارع نمو الشروخ (إشارة نقطة تحول)

المجال: الرصد الإنشائي الهدف: اكتشاف تحول الرتبة: \(D_2^1\)

إذا ارتفع معدل النمو النسبي للشرخ من 5% إلى 10% لكل دورة، فإن “معدل التغير النسبي” تضاعف.

\[ \frac{10\%}{5\%}=2 \Rightarrow \text{سلوك رتبة ثانية (تسارع خطر)}. \]

معادلة إضافية (مراقبة تغير المعدل):

\[ r=\frac{\Delta L}{L}\quad \Rightarrow \quad \text{مؤشر رتبة ثانية يراقب }\Delta r. \]

Discussion: كثير من الفشل يحدث عند بداية التسارع، وليس عند بلوغ الضرر قيمة كبيرة. لذلك تمثل الرتبة الثانية طبقة إنذار مبكر تحدد بدء تحول النظام نحو اللاانعكاسية.

الحالة 3 — زيادة زحف الخرسانة (Creep)

المجال: تشوهات طويلة الأمد الهدف: الخدمة الرتبة: \(D_1^1\)

ترتفع انفعالات الزحف من 400 με إلى 460 με خلال نافذة زمنية ثابتة.

\[ D_0^1:\;60\mu\varepsilon,\qquad D_1^1:\;\frac{60}{400}=15\%. \]

معادلة إضافية (صياغة معامل الزحف):

\[ \varepsilon(t)=\frac{\sigma}{E}\left(1+\varphi(t)\right). \]

Discussion: التغير النسبي في الزحف يسمح بالمقارنة بين عناصر مختلفة الإجهاد. إذا زاد المعدل النسبي لاحقًا، فذلك يشير إلى انتقال سلوكي، وهنا تظهر أهمية الرتبة الثانية لتفسير التسارع.

الحالة 4 — انكماش الجفاف وتأثيره على الشروخ

المجال: المتانة الهدف: التحكم بالتشققات الرتبة: \(D_1^1\)

يرتفع انكماش الجفاف من 250 με إلى 300 με.

\[ D_0^1:\;50\mu\varepsilon,\qquad D_1^1:\;\frac{50}{250}=20\%. \]

معادلة إضافية (إجهاد الانكماش المقيد):

\[ \sigma_{sh} \approx E\,\varepsilon_{sh}\,R. \]

Discussion: خطر التشقق يعتمد على درجة التقييد \(R\). لذا فإن القراءة النسبية تعكس زيادة الخطر بصورة مباشرة. إذا أصبح النمو النسبي نفسه يتسارع مع الزمن، فذلك يدل على انتقال نحو نطاق تشققات حاد.

الحالة 5 — تآكل التسليح وتأثيره على مساحة المقطع

المجال: متانة الهدف: فقدان مقاومة الرتبة: \(D_1^1\)

ينخفض قطر السيخ من 16.0 مم إلى 15.2 مم بسبب التآكل.

\[ D_0^1:\;-0.8\text{ mm},\qquad D_1^1:\;\frac{-0.8}{16}=-5\%. \]

معادلة إضافية (حساسية المساحة):

\[ A=\frac{\pi d^2}{4}\Rightarrow \frac{\Delta A}{A}\approx 2\frac{\Delta d}{d}\approx -10\%. \]

Discussion: نقص القطر بنسبة 5% يؤدي إلى نقص المساحة بنسبة تقارب 10% (وهذا مهم جدًا للقدرة). هذا المثال يوضح لماذا التغير النسبي أكثر فاعلية من الأرقام المطلقة في الحكم على الخطر.

الحالة 6 — هبوط الأساس (مطلق vs نسبي مع حدود)

المجال: جيوتكنيك الهدف: خدمة/حدود الرتبة: \(D_1^1\)

يرتفع الهبوط من 12 مم إلى 15 مم. الحد المسموح 25 مم. القرار يعتمد على الاقتراب من الحد.

\[ D_0^1:\;3\text{ mm},\qquad D_1^1:\;\frac{3}{12}=25\%,\qquad \eta=\frac{15}{25}=0.60. \]

معادلة إضافية (مؤشر الاستغلال):

\[ \eta(t)=\frac{S(t)}{S_{\max}},\qquad \eta\to 1 \Rightarrow \text{فشل خدمة}. \]

Discussion: مؤشر \(\eta\) يعبر مباشرة عن قرب المنشأ من حد عدم القبول. إذا بدأت \(\eta\) بالازدياد بوتيرة نسبية متسارعة، فهذه علامة على تحول في سلوك التربة (رتبة ثانية).

الحالة 7 — هبوط تفاضلي عبر بحر

المجال: جيوتكنيك/إنشاء الهدف: التشوه الرتبة: مؤشرات معيارية

هبوط دعامة يسار 8 مم ويمين 14 مم والبحر 10 م. المهم هو التشوه النسبي.

\[ \Delta S = 6\text{ mm},\qquad \beta=\frac{\Delta S}{L}=\frac{6}{10}=0.6\text{ mm/m}. \]

معادلة إضافية (مؤشر التشوه الزاوي):

\[ \beta=\frac{\Delta S}{L}. \]

Discussion: التشوه التفاضلي سبب أساسي للشقوق في العناصر الهشة. مراقبة \(\beta(t)\) أكثر فاعلية من مراقبة الهبوط المطلق لكل نقطة.

الحالة 8 — تقدم الانضغاط (Consolidation) وتسارع غير متوقع

المجال: Consolidation الهدف: تحقق النموذج الرتبة: \(D_1^1\) + \(D_2^1\)

تقدم الهبوط من 40% إلى 55% من القيمة النهائية. التغير النسبي 37.5%.

\[ D_0^1:\;+15\text{ نقطة},\qquad D_1^1:\;\frac{55-40}{40}=37.5\%. \]

معادلة إضافية (Terzaghi):

\[ T_v=\frac{c_v t}{H_d^2},\qquad U \approx f(T_v). \]

Discussion: تسارع التقدم النسبي قد يعني تغيرًا في التصريف أو طبقات التربة أو \(c_v\). إذا أصبح معدل التقدم النسبي يتغير بسرعة، فهذا تحول رتبة ثانية يستدعي مراجعة الفرضيات.

الحالة 9 — استقرار منحدر: تدهور عامل الأمان

المجال: منحدرات الهدف: أمان حدّي الرتبة: \(D_1^1\)

ينخفض FoS من 1.35 إلى 1.20 بعد أمطار. الانخفاض النسبي ≈ 11.1%.

\[ D_0^1:\;-0.15,\qquad D_1^1:\;\frac{-0.15}{1.35}=-11.1\%. \]

معادلة إضافية (Mohr–Coulomb):

\[ \tau_f = c' + \sigma' \tan\phi'. \]

Discussion: الأمطار تقلل \(\sigma'\) فتضعف مقاومة القص. إذا أصبحت الانخفاضات النسبية تتكرر وتزداد، فهذا تسارع رتبة ثانية نحو عدم الاستقرار.

الحالة 10 — هامش تحمل الأساس ومؤشر الاستغلال

المجال: Foundations الهدف: utilization الرتبة: \(D_1^1\)

الحمل المطلوب يرتفع من 900 إلى 980 kN، والقدرة 1200 kN.

\[ D_0^1:\;+80kN,\qquad \eta_1=0.75,\quad \eta_2=0.817, \quad D_1^1(\eta)=\frac{0.817-0.75}{0.75}=8.9\%. \]

معادلة إضافية (شرط حدّي):

\[ Q_d \le Q_u/\gamma\quad \Leftrightarrow\quad \eta \le \eta_{\max}. \]

Discussion: الاستخدام \(\eta\) هو مؤشر معياري. إذا بدأ \(\eta\) بالارتفاع بمعدل نسبي متسارع، فهذه إشارة رتبة ثانية لاقتراب التحول نحو خطر.

الحالة 11 — اختبار ركيزة: فقدان الصلابة

المجال: Piles الهدف: صلابة الرتبة: \(D_1^1\)

عند نفس الحمل 600 kN كان الهبوط 6 مم وأصبح 8 مم.

\[ D_0^1:\;+2mm,\qquad D_1^1:\;\frac{2}{6}=33.3\%. \]

معادلة إضافية (صلابة مقطعية):

\[ k_s=\frac{Q}{s}. \]

Discussion: الصلابة تقل من 100 إلى 75 kN/mm (نقص 25%). التغير النسبي أوضح من الزيادة المطلقة في الهبوط عند اتخاذ قرار تدخل.

الحالة 12 — تراجع مؤشر حالة سطح الجسر

المجال: Bridges الهدف: إدارة أصول الرتبة: \(D_1^1\)

ينخفض مؤشر الحالة من 82 إلى 74 خلال سنة.

\[ D_0^1:\;-8\text{ نقطة},\qquad D_1^1:\;\frac{-8}{82}=-9.76\%. \]

معادلة إضافية (نموذج اضمحلال):

\[ CI(t)=CI_0 e^{-kt}. \]

Discussion: إذا زادت قيمة \(k\) مع الزمن فهذا تسارع رتبة ثانية. وهنا يصبح التدخل ضرورياً قبل السقوط تحت حد سياسة الصيانة.

الحالة 13 — تراكم ضرر التعب (Miner)

المجال: Fatigue الهدف: قرب العتبة الرتبة: \(D_1^1\)

يرتفع مؤشر الضرر من 0.30 إلى 0.36.

\[ D_0^1:\;+0.06,\qquad D_1^1:\;\frac{0.06}{0.30}=20\%. \]

معادلة إضافية (Miner):

\[ D=\sum_i \frac{n_i}{N_i},\qquad D\to1 \Rightarrow \text{فشل}. \]

Discussion: قراءة \(D_1^1\) تعبر عن سرعة الاقتراب من الفشل. إذا زادت نسبة الزيادة مع زيادة \(D\)، فهذا يشير لتسارع رتبة ثانية.

الحالة 14 — زيادة الترخيم الخدمي

المجال: Structures الهدف: Serviceability الرتبة: \(D_1^1\)

الترخيم من 18 إلى 21 مم.

\[ D_0^1:\;+3mm,\qquad D_1^1:\;\frac{3}{18}=16.7\%. \]

معادلة إضافية (Scaling):

\[ \delta \propto \frac{L^4}{EI}. \]

Discussion: بسبب اعتماد الترخيم على \(L^4\)، فإن المقارنة المطلقة غير عادلة عبر بحور مختلفة. التغير النسبي يكشف تدهور الصلابة ويتيح مقارنة علمية.

الحالة 15 — تضخيم ديناميكي: زيادة الاهتزاز

المجال: Dynamics الهدف: سلامة تشغيلية الرتبة: \(D_1^1\)

RMS يزيد من 0.08g إلى 0.11g.

\[ D_0^1:\;+0.03g,\qquad D_1^1:\;\frac{0.03}{0.08}=37.5\%. \]

معادلة إضافية (SDOF):

\[ X(\omega)=\frac{F_0/k}{\sqrt{(1-r^2)^2+(2\zeta r)^2}}. \]

Discussion: ارتفاع الاستجابة قد يعني انخفاض الصلابة أو تغير التخميد. إذا تسارعت هذه الزيادة النسبية، فهذه علامة رتبة ثانية على تقدم الضرر.

الحالة 16 — ارتفاع نسبة الانجراف الزلزالي

المجال: Seismic الهدف: Performance الرتبة: \(D_1^1\)

الانجراف من 0.9% إلى 1.2%.

\[ D_0^1:\;+0.3\%\text{ نقاط},\qquad D_1^1:\;\frac{1.2-0.9}{0.9}=33.3\%. \]

معادلة إضافية:

\[ \theta=\frac{\Delta}{h}. \]

Discussion: الانجراف مؤشر معياري بطبيعته. تسارعه النسبي عبر أحداث متعددة يمثل انتقالًا رتبيًا نحو خطر الانهيار.

الحالة 17 — مؤشر خطر التسييل (موجب)

المجال: Liquefaction الهدف: Risk index الرتبة: \(D_2^1\) ممكن

مؤشر خطر \(R\) من 1.10 إلى 1.40.

\[ D_0^1:\;+0.30,\qquad D_1^1:\;\frac{0.30}{1.10}=27.3\%. \]

معادلة إضافية:

\[ FS_{liq}=\frac{CRR}{CSR}. \]

Discussion: بما أن \(R>1\) فإن المجال مناسب للرتبة الثانية. تسارع الزيادة النسبية عبر تحديثات الخطر يحدد نقطة تحول منهجية في المخاطر.

الحالة 18 — ازدياد النحر (Scour) عند ركائز الجسر

المجال: Hydraulics الهدف: استقرار الركيزة الرتبة: \(D_1^1\)

النحر من 0.80 إلى 0.92 م.

\[ D_0^1:\;+0.12m,\qquad D_1^1:\;\frac{0.12}{0.80}=15\%. \]

معادلة إضافية (Scaling):

\[ d_s \propto \left(\frac{V^2}{gD}\right)^m. \]

Discussion: بسبب لاخطية الاعتماد على السرعة \(V\)، فإن الزيادات النسبية هي الأقوى في التنبؤ. إذا تسارعت النسبة موسميًا فهذا انتقال رتبي يستدعي إنذارًا مبكرًا.

الحالة 19 — رتبة تقارب FEM (استخراج الأس)

المجال: FEM الهدف: تحقق التقارب الرتبة: \(D_1^1\)

\(h\) من 0.20 إلى 0.10، والخطأ من 0.040 إلى 0.010.

\[ p=\frac{\ln(E_2/E_1)}{\ln(h_2/h_1)} =\frac{\ln(0.25)}{\ln(0.5)}=2. \]

معادلة إضافية:

\[ E(h)\approx C h^p. \]

Discussion: هذا مثال واضح أن الرتبة الأولى تستخرج مباشرة أس التقارب. أي أن الحساب الهرمي يظهر طبيعيًا في النمذجة العددية دون إضافات.

الحالة 20 — تضخم تكلفة المشروع وتسارع المخاطر

المجال: إدارة مشاريع الهدف: السيطرة الرتبة: \(D_1^1\) + \(D_2^1\)

التكلفة من 10M إلى 12M: زيادة نسبية 20%.

\[ D_0^1:\;+2M,\qquad D_1^1:\;\frac{2}{10}=20\%. \]

معادلة إضافية (Earned Value):

\[ CPI=\frac{EV}{AC},\qquad CPI<1 \Rightarrow \text{تجاوز}. \]

Discussion: الفشل غالبًا يحدث بسبب التسارع: 20% ثم 35% ثم 50%. الرتبة الثانية هي الأداة الطبيعية لرصد “تزايد معدل التجاوز” مبكرًا.

5) دليل عملي: متى أستخدم أي رتبة؟

السؤال الهندسي	الرتبة	السبب
كم تغيرت القيمة بالمطلق؟	\(D_0^1\)	تقارير مباشرة
ما أهمية التغير نسبةً للحالة؟	\(D_1^1\)	مقارنة عادلة عبر الأحجام
هل التغير النسبي نفسه يتسارع؟	\(D_2^1\)	تحول/إنذار مبكر
سلوك متعدد الأطوار طويل الأمد	\(D_n^1\)	تحليل هرمي أعلى

الخلاصة الهندسية: كثير من الانهيارات ليست “قيمة كبيرة”، بل “تسارع في التدهور النسبي”.

6) نقاش عام + قيود + عمل مستقبلي

6.1 نقاش عام

يتكرر نمط مشترك في كل الحالات: القرارات تعتمد على الاقتراب من حدود، وعلى المؤشرات النسبية المعيارية. لذا تظهر الرتبة الأولى كأداة مركزية. أما الرتبة الثانية فتضيف طبقة إنذار مبكر ترصد تغير سلوك النظام قبل الانهيار.

6.2 قيود

صيغة \(D_2^1\) الصارمة تتطلب مجالًا مناسبًا. عمليًا نستخدم مؤشرات موجبة أو تحويلات. كما أن البيانات قد تكون ضجيجية، ولذلك يلزم ترشيح واحتمالية لاستخراج التسارع بثقة.

6.3 عمل مستقبلي

يتضمن العمل المستقبلي: (1) تعريف تحويلات معيارية موجبة لمؤشرات الصحة، (2) دمج الرتب ضمن نماذج Bayesian SHM، (3) اختبار إنذار الرتبة الثانية على بيانات جسور حقيقية، (4) توسيع التحليل إلى رتب أعلى للأنظمة متعددة الأطوار.

المراجع (References)

المراجع أدناه مصادر معيارية للسياق الهندسي. إطار الحساب الهرمي هو جزء من منهج المؤلف.

\[ \textbf{[1]}\;\text{Terzaghi, K. (1943). Theoretical Soil Mechanics. Wiley.} \] \[ \textbf{[2]}\;\text{Fellenius, B. H. (2001). Basics of Foundation Design.} \] \[ \textbf{[3]}\;\text{Eurocode 2: Design of Concrete Structures.} \] \[ \textbf{[4]}\;\text{Miner, M. A. (1945). Cumulative Damage in Fatigue. Journal of Applied Mechanics.} \] \[ \textbf{[5]}\;\text{Bathe, K. J. (1996). Finite Element Procedures. Prentice Hall.} \] \[ \textbf{[6]}\;\text{GOSSA AHMED (2025). Hierarchical Calculus Framework (Official Website).} \]