هندسة التحكم 2: تصميم التحكم بالتغذية الراجعة: حزمة المراجعة المكثفة

أوراق القوانين، وأهم خطوات الحل لكل موضوع، و10 مسألة تدريبية، مع الحلول الكاملة في الخلف.

مسائل جديدةالعودة إلى المقرّر

في نافذة الطباعة اختر «حفظ بصيغة PDF» لحفظ نسخة على جهازك.

1القوانين وطرق الحل

الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور

الاستقرار
جميع جذور كثير الحدود المميِّز في النصف الأيسر من المستوى
جدول راوث، عنصر الصف s¹ في كثير حدود تكعيبي
b1 = (a2a1 − a3a0)/a2
جذور النصف الأيمن من المستوى
عدد تغيّرات الإشارة في العمود الأول
كثير الحدود التكعيبي s³ + a₂s² + a₁s + a₀
مستقر حين a > 0 لجميع المعاملات وa₂·a₁ > a₀
K/(s(s + a)(s + b))
Kmax = ab(a + b), ω = √(ab)
المعادلة المميزة
1 + K·G(s) = 0
الخطوط المقاربة
σa = (Σp − Σz)/(n − m), θ = (2q + 1)·180°/(n − m)
نقطتا الانفصال والالتقاء
dK/ds = 0
شرط المقدار
K = 1/|G(s)|
زاوية المغادرة
θd = 180° − Σ∠(from other poles) + Σ∠(from zeros)

الطريقة

  1. اكتب كثير الحدود المميِّز للحلقة المغلقة: مقام 1 + K·G(s) = 0، بعد فكّه.
  2. ضع المعاملات بالتناوب في الصفين الأولين من جدول راوث وابنِ بقية الصفوف.
  3. عُدّ تغيّرات الإشارة في العمود الأول: كل تغيّر قطب في النصف الأيمن من المستوى. ومع وجود K في الجدول اشترط أن يبقى كل عنصر في العمود الأول موجبًا.
  4. للمحل الهندسي للجذور عيّن الأقطاب والأصفار وقطع المحور الحقيقي والخطوط المقاربة من المركز ونقاط الانفصال حيث dK/ds = 0. واستخدم K = 1/|G(s)| لإيجاد الكسب عند نقطة.

خطأ الحالة المستقرة والتحكم PID

نوع النظام
عدد أقطاب G(s) عند s = 0
ثوابت الخطأ
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
الخطوة A
ess = A/(1 + Kp)
المنحدر A·t
ess = A/Kv
القطع المكافئ A·t²/2
ess = A/Ka
المتحكم PID
C(s) = Kp + Ki/s + Kd·s
الهدف من الرتبة الثانية
s² + 2ζωns + ωn²
التجاوز والتوقيت
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

الطريقة

  1. عُدّ أقطاب الحلقة المفتوحة C(s)·G(s) عند s = 0: فذلك هو نوع النظام.
  2. احسب ثابت الخطأ المناسب للدخل: Kp للخطوة، وKv للمنحدر، وKa للقطع المكافئ.
  3. طبّق ess = A/(1 + Kp) أو A/Kv أو A/Ka، أو حُلّها لإيجاد الكسب المطلوب.
  4. لمواصفة التخميد أو السرعة اكتب كثير الحدود المميِّز للحلقة المغلقة وطابقه مع s² + 2ζω_n·s + ω_n².

2مسائل تدريبية

1.الاستجابة الخطويةمتوسطة
t (s)y(t)3?

لنظام دالة النقل G(s) = 66/(s² + 7s + 22) ويتلقى دخلًا خطويًا وحدويًا. أوجد النسبة المئوية للتجاوز %OS في الاستجابة الخطوية.

G(s) = 66/(s² + 7s + 22)
%OS = %
2.الاستجابة الخطويةصعبة
t (s)y(t)1tp​ = 0,25 s20 %

تُظهر الاستجابة الخطوية لنظام من الرتبة الثانية بلا أصفار تجاوزًا قدره 20 %، وتقع الذروة الأولى عند t = 0,25 s. أوجد التردد الطبيعي ωn.

%OS = 20 %tp = 0,25 s
ωn = rad/s
3.مخططات بودي والهوامشمتوسطة
|H| (dB)0,11101001000−80−60−40−20020ω (rad/s)ω = 30

أوجد المقدار الدقيق |G(jω)| لـ G(s) = 616/((s + 8)(s + 11)) عند ω = 30 rad/s، بوحدة dB.

G(s) = 616/((s + 8)(s + 11))ω = 30 rad/s
|G(jω)| = dB
4.مخططات بودي والهوامشصعبة
|H| (dB)0,10,1111010100100−80−60−40−2002040ω (rad/s)φ (°)ω (rad/s)−270−180−900ω180​

حلقة تغذية عكسية وحدوية دالة نقلها المفتوحة L(s) = K/(s(s + 5)(s + 6)). أوجد الكسب K الذي يعطي هامش كسب قدره 6 dB.

L(s) = K/(s(s + 5)(s + 6))GM = 6 dB
K =
5.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذورمتوسطة
R(s)+K1s(s + 4)(s + 10)C(s)−

لحلقة المخطط دالة نقل الحلقة المفتوحة L(s) = K/(s(s + 4)(s + 10)). أوجد أكبر كسب Kmax تكون الحلقة المغلقة عنده مستقرة.

L(s) = K/(s(s + 4)(s + 10))
Kmax =
6.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذورصعبة
ReIm−1−2−3−4−5−6−7−81j1j2

يغادر المحل الهندسي للجذور لـ L(s) = K(s + 7)/(s(s + 2)) المحور الحقيقي بين s = 0 وs = −2، ثم يدور ويلتقي بالمحور الحقيقي مجددًا إلى يسار الصفر. أوجد نقطة الانفصال sb.

L(s) = K(s + 7)/(s(s + 2))
sb =
7.خطأ الحالة المستقرة والتحكم PIDمتوسطة
r(t) = 5t+Ks(s + 5)C(s)−

يجب أن تتبع حلقة تغذية عكسية وحدوية فيها G(s) = K/(s(s + 5)) المنحدر r(t) = 5t بخطأ حالة مستقرة قدره 0,05. أوجد الكسب K.

G(s) = K/(s(s + 5))r(t) = 5tess = 0,05
K =
8.خطأ الحالة المستقرة والتحكم PIDصعبة
R(s)+K = 402s(s + 4)C(s)−

يتحكم متحكم تناسبي K = 40 في المحطة G(s) = 2/(s(s + 4)) ضمن حلقة تغذية عكسية وحدوية. أوجد زمن الاستقرار ts عند 2 % للاستجابة الخطوية باستخدام ts ≈ 4/(ζωₙ).

G(s) = 2/(s(s + 4))K = 40
ts = s
9.المخططات الصندوقيةمتوسطة
R(s)+4(s + 7)(s + 1)(s + 2)G(s)C(s)−

حلقة تغذية راجعة سالبة أحادية لها G(s) = 4(s + 7)/((s + 1)(s + 2)). تُطبَّق خطوة مقدارها 1 عند المرجع. أوجد خطأ الحالة المستقرة ess.

G(s) = 4(s + 7)/((s + 1)(s + 2))R = 1
ess =
10.المخططات الصندوقيةصعبة
R(s)+17s + 5G₁+6s + 3G₂C(s)0,5H₁−0,2H₂−

اختزل المخطط الصندوقي مع G₁ = 17/(s + 5) و G₂ = 6/(s + 3) وتغذية راجعة سالبة داخلية H₁ = 0,5 حول G₂ وتغذية راجعة سالبة خارجية H₂ = 0,2. أوجد كسب التيار المستمر T(0) للنسبة C(s)/R(s).

G₁ = 17/(s + 5)G₂ = 6/(s + 3)H₁ = 0,5H₂ = 0,2
T(0) =

3المحاليل

1.الاستجابة الخطوية

الإجابة: %OS = 2,956 %

الاستجابة للدخل الخطوي: النظام من الرتبة الأولى G(s) = K/(τs + 1) يستجيب لخطوة مقدارها U بالعلاقة y(t) = K·U·(1 − e−t/τ): يبلغ 63,2 % من القيمة النهائية K·U بعد ثابت زمني واحد τ و98 % بعد 4τ. والنظام من الرتبة الثانية G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) له القطبان s = −ζωₙ ± jωd مع ωd = ωₙ√(1 − ζ²). وعندما 0 < ζ < 1 تتجاوز الاستجابة الخطوية القيمة النهائية: النسبة المئوية للتجاوز %OS = 100·e−ζπ/√(1−ζ²)، وزمن الذروة tp = π/ωd، وزمن الاستقرار عند 2 % هو ts ≈ 4/(ζωₙ). قارن المقام مع s² + 2ζωₙs + ωₙ² لقراءة ωₙ وζ. وفي نظام كتلة–نابض–مخمّد ωₙ = √(k/m) وζ = c/(2√(km))؛ وفي دائرة RLC على التسلسل ωₙ = 1/√(LC) وζ = (R/2)·√(C/L).

قارن مع الصيغة القياسية

يجب أن يطابق المقام s² + 2ζωₙs + ωₙ².

ωₙ² = 22 → ωₙ = 4,69 rad/s
2ζωₙ = 7 → ζ = 7/(2·4,69) = 0,7462

النتيجة

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−3,521 = 2,956 %

الإجابة

%OS = 2,956 %
2.الاستجابة الخطوية

الإجابة: ωn = 14,12 rad/s

الاستجابة للدخل الخطوي: النظام من الرتبة الأولى G(s) = K/(τs + 1) يستجيب لخطوة مقدارها U بالعلاقة y(t) = K·U·(1 − e−t/τ): يبلغ 63,2 % من القيمة النهائية K·U بعد ثابت زمني واحد τ و98 % بعد 4τ. والنظام من الرتبة الثانية G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) له القطبان s = −ζωₙ ± jωd مع ωd = ωₙ√(1 − ζ²). وعندما 0 < ζ < 1 تتجاوز الاستجابة الخطوية القيمة النهائية: النسبة المئوية للتجاوز %OS = 100·e−ζπ/√(1−ζ²)، وزمن الذروة tp = π/ωd، وزمن الاستقرار عند 2 % هو ts ≈ 4/(ζωₙ). قارن المقام مع s² + 2ζωₙs + ωₙ² لقراءة ωₙ وζ. وفي نظام كتلة–نابض–مخمّد ωₙ = √(k/m) وζ = c/(2√(km))؛ وفي دائرة RLC على التسلسل ωₙ = 1/√(LC) وζ = (R/2)·√(C/L).

نسبة التخميد

ζ = −ln(0,2)/√(π² + ln²(0,2)) = 1,609/√(9,87 + 2,59) = 0,4559

التردد الطبيعي

ωd = π/tp = π/0,25 = 12,57 rad/s
ωₙ = ωd/√(1 − ζ²) = 12,57/0,89 = 14,12 rad/s

الإجابة

ωn = 14,12 rad/s
3.مخططات بودي والهوامش

الإجابة: |G(jω)| = −4,139 dB

الاستجابة الترددية: عوّض s = jω في دالة النقل. المقدار بالديسيبل هو 20·log₁₀|H(jω)| والطور هو مجموع زوايا عوامل البسط ناقص زوايا عوامل المقام. مرشح RC أو RL من الرتبة الأولى له تردد الركن (−3 dB) عند ωc = 1/(RC) أو R/L، وfc = ωc/(2π)؛ ومرشح تمرير الترددات المنخفضة يعطي |H| = 1/√(1 + (ω/ωc)²) وφ = −atan(ω/ωc). وفي مخطط بودي التقريبي (الخطوط المستقيمة) ينحني كل قطب بالمقدار إلى الأسفل بـ 20 dB/decade عند تردد ركنه، وينحني كل صفر إلى الأعلى بـ 20 dB/decade، ويعطي المكامل 1/s ميلًا قدره −20 dB/decade عند كل الترددات. تردد تقاطع الكسب ωc: |L(jω)| = 1 (0 dB)؛ وهامش الطور PM = 180° + ∠L(jωc). تردد تقاطع الطور ω180: ∠L(jω) = −180°؛ وهامش الكسب GM = −20·log₁₀|L(jω180)| dB. تكون الحلقة المغلقة مستقرة حين يكون الهامشان موجبين.

اضرب المقادير

|G(j30)| = 616/(√(30² + 8²)·√(30² + 11²)) = 616/(31,05·31,95) = 0,6209

النتيجة

20·log₁₀(0,6209) = −4,139 dB

الإجابة

|G(jω)| = −4,139 dB
4.مخططات بودي والهوامش

الإجابة: K = 165,4

الاستجابة الترددية: عوّض s = jω في دالة النقل. المقدار بالديسيبل هو 20·log₁₀|H(jω)| والطور هو مجموع زوايا عوامل البسط ناقص زوايا عوامل المقام. مرشح RC أو RL من الرتبة الأولى له تردد الركن (−3 dB) عند ωc = 1/(RC) أو R/L، وfc = ωc/(2π)؛ ومرشح تمرير الترددات المنخفضة يعطي |H| = 1/√(1 + (ω/ωc)²) وφ = −atan(ω/ωc). وفي مخطط بودي التقريبي (الخطوط المستقيمة) ينحني كل قطب بالمقدار إلى الأسفل بـ 20 dB/decade عند تردد ركنه، وينحني كل صفر إلى الأعلى بـ 20 dB/decade، ويعطي المكامل 1/s ميلًا قدره −20 dB/decade عند كل الترددات. تردد تقاطع الكسب ωc: |L(jω)| = 1 (0 dB)؛ وهامش الطور PM = 180° + ∠L(jωc). تردد تقاطع الطور ω180: ∠L(jω) = −180°؛ وهامش الكسب GM = −20·log₁₀|L(jω180)| dB. تكون الحلقة المغلقة مستقرة حين يكون الهامشان موجبين.

تقاطع الطور

(ω/5)·(ω/6) = 1 → ω180 = √(5·6) = 5,477 rad/s

المقدار عند تقاطع الطور

|L(jω180)| = K/(ω·√(ω² + 25)·√(ω² + 36)) = K/(5·6·(5 + 6)) = K/330

الكسب المطلوب

330/K = 106/20 = 1,995
K = 330/1,995 = 165,4

الإجابة

K = 165,4
5.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور

الإجابة: Kmax = 560

معيار راوث–هرفيتز: يكون النظام الخطي مستقرًا حين يقع كل جذر من جذور كثير حدوده المميِّز في النصف الأيسر من المستوى s. اكتب المعاملات بالتناوب في الصفين الأولين من جدول راوث. يأتي كل عنصر في صف جديد من الصفين اللذين فوقه؛ وبالنسبة إلى كثير حدود تكعيبي a₃s³ + a₂s² + a₁s + a₀ يكون عنصر الصف s¹ هو b₁ = (a₂·a₁ − a₃·a₀)/a₂. عدد تغيّرات الإشارة في العمود الأول يساوي عدد الجذور في النصف الأيمن من المستوى. وحين يدخل الكسب K في المعاملات يجب أن يبقى كل عنصر في العمود الأول موجبًا: الكسب الذي يجعل عنصر الصف s¹ صفرًا هو حدّ الاستقرار، ثم يعطي الصف s² تردد التذبذب المستمر.

المعادلة المميزة

s³ + 14s² + 40s + K = 0

جدول روث

s³: 1, 40
s²: 14, K
s¹: (14·40 − K)/14
s⁰: K

الشروط

يجب أن يكون كل عنصر في العمود الأول موجبًا.

K > 0
14·40 − K > 0 → K < 560

الإجابة

Kmax = 560

عند K = 560 تتذبذب الحلقة المغلقة بالتردد ω = √40 = 6,325 rad/s.

6.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور

الإجابة: sb = −1,084

المحل الهندسي للجذور: حين يتغير الكسب K من 0 إلى ∞ تتحرك أقطاب الحلقة المغلقة، وهي جذور 1 + K·G(s) = 0، من أقطاب الحلقة المفتوحة إلى أصفار الحلقة المفتوحة أو إلى ما لا نهاية. تنتمي نقطة على المحور الحقيقي إلى المحل الهندسي حين يكون عدد الأقطاب والأصفار الحقيقية على يمينها فرديًا. تتبع الفروع n − m المتجهة إلى ما لا نهاية خطوطًا مقاربة بالزوايا (2q + 1)·180°/(n − m)، وتنطلق من المركز σₐ = (Σ الأقطاب − Σ الأصفار)/(n − m). تغادر الفروع المحور الحقيقي أو تنضم إليه حيث dK/ds = 0. يحدد شرط الزاوية ∠G(s) = 180° الشكل، ويعطي شرط المقدار K = 1/|G(s)| الكسب عند أي نقطة من المحل الهندسي.

الكسب على طول المحور الحقيقي

K(s) = −(s² + 2s)/(s + 7)

النقاط الحرجة

اشتق بقاعدة القسمة وساوِ البسط بالصفر.

(2s + 2)(s + 7) − (s² + 2s) = s² + 14s + 14 = 0
s = −7 ± √35 = −1,084, −12,92

نقطة الانفصال

sb = −1,084

الإجابة

sb = −1,084

بين نقطتي الانفصال والالتقاء يكون المحل الهندسي دائرة مركزها الصفر، ونصف قطرها √35 = 5,916.

7.خطأ الحالة المستقرة والتحكم PID

الإجابة: K = 500

خطأ الحالة المستقرة: في حلقة تغذية عكسية وحدوية بدالة نقل الحلقة المفتوحة G(s) يكون الخطأ E(s) = R(s)/(1 + G(s))، وبنظرية القيمة النهائية ess = lims→0 s·E(s). نوع النظام هو عدد أقطاب G(s) عند s = 0. ثوابت الخطأ هي Kp = lim G(s) وKv = lim s·G(s) وKa = lim s²·G(s)، كلها عندما s → 0. تعطي دالة خطوة بحجم A الخطأ ess = A/(1 + Kp)، وتعطي دالة منحدرة A·t الخطأ A/Kv، وتعطي دالة قطع مكافئ A·t²/2 الخطأ A/Ka. تتبع حلقة من النوع 0 دالة الخطوة بخطأ منتهٍ، أما حلقة من النوع 1 فتتبعها تمامًا وتتبع دالة المنحدر بخطأ منتهٍ.

ثابت السرعة المطلوب

Kv = A/ess = 5/0,05 = 100

الكسب

K = 5·Kv = 500

الإجابة

K = 500

الكسب الكبير يجعل هذه الحلقة قليلة التخميد جدًا: ζ = 0,1118. وهذا التعارض هو سبب استخدام المعوِّضات.

8.خطأ الحالة المستقرة والتحكم PID

الإجابة: ts = 2 s

التحكم P وPI وPID: يعمل المتحكم C(s) = Kp + Ki/s + Kd·s على الخطأ. الكسب التناسبي يجعل الحلقة أسرع ويقلل خطأ محطة من النوع 0، لكنه لا يزيله أبدًا. يضيف الحد التكاملي قطبًا عند s = 0: فيرتفع نوع النظام بواحد، وتُتبع الخطوة دون خطأ. ويضيف الحد التفاضلي تخميدًا. لاختيار الكسوب اكتب كثير الحدود المميِّز للحلقة المغلقة وقارنه بـ s² + 2ζωns + ωn². عندئذٍ %OS = 100·e−ζπ/√(1−ζ²)، وtp = π/ωd مع ωd = ωn√(1 − ζ²)، وts ≈ 4/(ζωn).

كثير الحدود المميِّز

s(s + 4) + 40·2 = s² + 4s + 80

ωₙ وζ

ωₙ = √80 = 8,944 rad/s
ζ = 4/(2·8,944) = 0,2236

زمن الاستقرار

ζωₙ = 4/2 = 2
ts ≈ 4/2 = 2 s

الإجابة

ts = 2 s

لا يعتمد زمن الاستقرار هنا على K: فرفع K لا يجعل الاستجابة إلا أسرع تذبذبًا.

9.المخططات الصندوقية

الإجابة: ess = 0,06667

جبر المخططات الصندوقية: الصناديق على التوالي تُضرب (G₁·G₂)، والصناديق على التوازي تُجمع (G₁ + G₂)، والحلقة ذات المسار الأمامي G ومسار التغذية الراجعة H تُختزل إلى T = G/(1 + G·H) للتغذية الراجعة السالبة و G/(1 − G·H) للتغذية الراجعة الموجبة. اختزل الحلقات الداخلية أولًا، ثم اتجه إلى الخارج. أقطاب الحلقة المغلقة هي جذور المعادلة المميّزة 1 + G(s)H(s) = 0. لحلقة مستقرة، دخل خطوة مقداره A يعطي خرج الحالة المستقرة A·T(0)؛ ومع التغذية الراجعة الأحادية يكون خطأ الحالة المستقرة لخطوة A/(1 + Kp) مع Kp = G(0).

ثابت خطأ الموضع

Kₚ = G(0) = 4·7/(1·2) = 14

خطأ الحالة المستقرة

ess = A/(1 + Kₚ) = 1/(1 + 14) = 0,06667

الإجابة

ess = 0,06667
10.المخططات الصندوقية

الإجابة: T(0) = 2,024

جبر المخططات الصندوقية: الصناديق على التوالي تُضرب (G₁·G₂)، والصناديق على التوازي تُجمع (G₁ + G₂)، والحلقة ذات المسار الأمامي G ومسار التغذية الراجعة H تُختزل إلى T = G/(1 + G·H) للتغذية الراجعة السالبة و G/(1 − G·H) للتغذية الراجعة الموجبة. اختزل الحلقات الداخلية أولًا، ثم اتجه إلى الخارج. أقطاب الحلقة المغلقة هي جذور المعادلة المميّزة 1 + G(s)H(s) = 0. لحلقة مستقرة، دخل خطوة مقداره A يعطي خرج الحالة المستقرة A·T(0)؛ ومع التغذية الراجعة الأحادية يكون خطأ الحالة المستقرة لخطوة A/(1 + Kp) مع Kp = G(0).

الحلقة الداخلية

Gᵢ = G₂/(1 + G₂H₁) = 6/(s + 6)

الحلقة الخارجية

T(s) = G₁Gᵢ/(1 + G₁GᵢH₂) = 102/((s + 5)(s + 6) + 20,4)

كسب التيار المستمر

T(0) = 102/(5·6 + 20,4) = 102/50,4 = 2,024

الإجابة

T(0) = 2,024

ExamProblems · examproblems.com/courses/control-engineering-2