[Bιομηχανική Πληροφορική] Εργασία 2016/2017

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

[Admins]

Για αλλαγή του public name σας, επικοινωνήστε με έναν από τους Admins.

[tniton]

Ναι, με α=0,2397, β=-5,8393 και γ=0,6851. Σωστο δεν ειναι; Εχουμε λαθος λογική;

[Goemon]

+1
Επίσης μπορεί κάποιος να εξηγήσει πως κάνει τον pid γιατί έμενα απειρίζονται οι όροι ??

• Έφτιαξα το μοντέλο στο matlab γιατί παίρνει ΠΟΛΥ λιγότερη ώρα από τον να τον ρυθμίσεις μέσω trace στο CODSYS. Χρησιμοποίησε την tf αν θυμάσαι τίποτα από ΣΑΕ...
• Δες αυτές τις οδηγίες: https://en.wikipedia.org/wiki/Ziegler%E2%80%93Nichols_method και https://www.youtube.com/watch?v=nvAQHSe-Ax4 (το βιντεο έχει μερικά λαθάκια, απλά αξιοποίησε το για να καταλάβεις ΤΙ ειναι το Ku και το Tu)
• Με χρήση της step ΚΑΙ ΤΙΠΟΤΑ πιο πολύπλοκο βρες τα Ku και Tu στο matlab (βάσει των οδηγιών πιο πάνω)
• Βάλε το PID στο CODESYS. Τοε΄φτιαξες μονο σου ή πήρες το έτοιμο? Πες για να σε πω αν θες
• Βάλε τις παραμέτρους που έχεις στο PID του CODESYS: ΠΡΟΣΟΧΗ: Το PID στο codesys δεν είναι το κλασσικό με τα Kp, Kd, Ki αλλά στην μορφή Κ, Τd, Ti όπως τα έχει στην βικιπέδια.

[sotiristsar]

• Έφτιαξα το μοντέλο στο matlab γιατί παίρνει ΠΟΛΥ λιγότερη ώρα από τον να τον ρυθμίσεις μέσω trace στο CODSYS. Χρησιμοποίησε την tf αν θυμάσαι τίποτα από ΣΑΕ...
• Δες αυτές τις οδηγίες: https://en.wikipedia.org/wiki/Ziegler%E2%80%93Nichols_method και https://www.youtube.com/watch?v=nvAQHSe-Ax4 (το βιντεο έχει μερικά λαθάκια, απλά αξιοποίησε το για να καταλάβεις ΤΙ ειναι το Ku και το Tu)
• Με χρήση της step ΚΑΙ ΤΙΠΟΤΑ πιο πολύπλοκο βρες τα Ku και Tu στο matlab (βάσει των οδηγιών πιο πάνω)
• Βάλε το PID στο CODESYS. Τοε΄φτιαξες μονο σου ή πήρες το έτοιμο? Πες για να σε πω αν θες
• Βάλε τις παραμέτρους που έχεις στο PID του CODESYS: ΠΡΟΣΟΧΗ: Το PID στο codesys δεν είναι το κλασσικό με τα Kp, Kd, Ki αλλά στην μορφή Κ, Τd, Ti όπως τα έχει στην βικιπέδια.

Αν το εχουμε φτιαξει μονοι μας ??

Δες τον κωδικα που καναμε στο 5ο εργαστηριο, λεει εκει πως τα οριζεις.

Για να βρειτε την τελικη τιμη της θερμοκρασιας μετα τον ελεγχο pid πηρατε την εξισωση διαφορων?

Απο το 5ο εργαστηριο δεν μπορω να καταλαβω πως προκύπτουν οι αρχικες τιμες  που εχει και αρα τι πρεπει να αρχικοποιήσουμε εμεις ..

[Goemon]

Αν το εχουμε φτιαξει μονοι μας ??

Δες ποια είναι ημορφή με τα K,Ti,Td, δες ποια έιναι η δικιά σου μορφή, κάνε την άλγεβρα και βγάλε τα κέρδη του δικού σου ελεγκτή συναρτήση των υπολογισμενων Κ,Τi,Td

[ori0ngel]

Εχω βρει την εξισωση διαφορων την βαζω στο matlab, το 100 που ψαχνουμε ειναι στο πεδιο τιμων της... ολα καλα μεχρι εδω.

Αυτο που δεν μου βγαινει ακομα ειναι οι οροι του PID. Εχουμε απο την εξισωση (2) της εκφωνησης οτι Kss=5.6, L=0.5 και Τ=0.7.

Εστω οτι το υλοποιουμε με ziegler nichols. απο την διαφανεια 21 παιρνουμε τις παραμετρους του pid. οταν ομως τις βαζω στις εξισωσεις yokogava απειριζονται ολα. (η θερμοκρασια προς το -απειρο και η θερμοτητα στο +απειρο).

Ακομα και με την συναρτηση μεταφορας βαζοντας την στο ματλαβ βρισκοντας τελειως διαφορετικες παραμετρους η ισχυς του θερμαντηρα αυξανοταν και παλι.

Χρησιμοποιω την υλοποιηση που ειχαμε κανει στο 5ο εργαστηριο με την διαφορα οτι για τον τελικο υπολογισμο της θερμοκρασιας χρησιμοποιω την εξισωση διαφορων και στην εναρξη καθε βηματος του ελεγχου pid εισαγω ενα απλο σφαλμα temp=temp-1 ωστε να προσομοιωσω την πτωση της θερμοκρασιας.
ειναι σωστο αυτο το σκεπτικο, αρα εχω προβλημα μονο στην ευρεση των παραμετρων ή ειναι ολο λαθος?

[nikos1]

• Έφτιαξα το μοντέλο στο matlab γιατί παίρνει ΠΟΛΥ λιγότερη ώρα από τον να τον ρυθμίσεις μέσω trace στο CODSYS. Χρησιμοποίησε την tf αν θυμάσαι τίποτα από ΣΑΕ...
• Δες αυτές τις οδηγίες: https://en.wikipedia.org/wiki/Ziegler%E2%80%93Nichols_method και https://www.youtube.com/watch?v=nvAQHSe-Ax4 (το βιντεο έχει μερικά λαθάκια, απλά αξιοποίησε το για να καταλάβεις ΤΙ ειναι το Ku και το Tu)
• Με χρήση της step ΚΑΙ ΤΙΠΟΤΑ πιο πολύπλοκο βρες τα Ku και Tu στο matlab (βάσει των οδηγιών πιο πάνω)
• Βάλε το PID στο CODESYS. Τοε΄φτιαξες μονο σου ή πήρες το έτοιμο? Πες για να σε πω αν θες
• Βάλε τις παραμέτρους που έχεις στο PID του CODESYS: ΠΡΟΣΟΧΗ: Το PID στο codesys δεν είναι το κλασσικό με τα Kp, Kd, Ki αλλά στην μορφή Κ, Τd, Ti όπως τα έχει στην βικιπέδια.

δηλαδη βαζοντας τις παραμετρους που βρηκαμε απο τη μεθοδο Zieglers - Nichols στο ετοιμο block PID του CODESYS θα λειτουργει? Εμενα μου εβγαζε οτι να 'ναι...εκτος αν εκανα καποια ανοησια...

[Goemon]

δηλαδη βαζοντας τις παραμετρους που βρηκαμε απο τη μεθοδο Zieglers - Nichols στο ετοιμο block PID του CODESYS θα λειτουργει? Εμενα μου εβγαζε οτι να 'ναι...εκτος αν εκανα καποια ανοησια...

Εμένα δεν είχε θέμα πάντως...  τις παραμέτρους του PID τις μετέφερες σωστά στο Function Block?

Η εξίσωση διαφορών είναι σωστή?
Αν έβγαλες την εξίσωση διαφορών προσσεγίζοντας την διαφορική με τον τύπο ydot = (y[n+1]-y[n])/Ts, ο τύπος δεν θα έπρεπε να απέχει πολύ από αυτόν που βγαίνει με τον μετασχηματισμό Tustin (s = 2/Ts*(1+z)/(1-z)), οπότε αναμένω no-overshoot εκδοχή του Zigler Nic να δουλέψει μια χαρά πάλι (δεν το δοκίμασα αλλά έβγαλα την εξ. διαφορών και είδα οτι μοιάζουν τρομερά)

Αν σου πηγαίνει μεταβίας μέχρι τους 85C, αύξησε το Y_MAX στα 20. Eγώ το έβαλα 15000 γιατί θεώρησα οτι μερικά Watt/sec αρκούν να αντισταθμίσουν τα 10m^3/sec Pου θα μπαίνουν, οπότε θεώρησα οτι η συνάρτηση μεταφοράς είναι σε Celsius/Watt.

[tniton]

Παιδια, δυο βασικα πραγματα που βρηκα σημερα, εστω και τελευταια στιγμη για καποιους, για να μην τρωτε τσαμπα χρονο. Η εξισωση στη διαφανεια 11 εχει το προσημο του β λαθος, το σωστο ειναι -β(για να βγει θετικος ο ορος). Επισης το χρονο δειγματοληψιας τον εβαλα Τ/5=0,7/5=0,14 ( Το είχε πει ο Χασαπης στις συχνες ερωτησεις μιας παλιας εργασιας(2011-12) νομιζω στο ethmmy.)

[Goemon]

Παιδια, δυο βασικα πραγματα που βρηκα σημερα, εστω και τελευταια στιγμη για καποιους, για να μην τρωτε τσαμπα χρονο. Η εξισωση στη διαφανεια 11 εχει το προσημο του β λαθος, το σωστο ειναι -β(για να βγει θετικος ο ορος). Επισης το χρονο δειγματοληψιας τον εβαλα Τ/5=0,7/5=0,14 ( Το είχε πει ο Χασαπης στις συχνες ερωτησεις μιας παλιας εργασιας(2011-12) νομιζω στο ethmmy.)

Τόσο λίγο? Εγώ το έβαλα 50ms καθώς η όλη διαδικασία δεν είναι ιδιαίτερα ραγδιαία... θα αλλάξω τα νούμερά μου ίσως για να με σε καλή μεριά

[Goemon]

Εμένα δεν είχε θέμα πάντως...  τις παραμέτρους του PID τις μετέφερες σωστά στο Function Block?

Η εξίσωση διαφορών είναι σωστή?
Αν έβγαλες την εξίσωση διαφορών προσσεγίζοντας την διαφορική με τον τύπο ydot = (y[n+1]-y[n])/Ts, ο τύπος δεν θα έπρεπε να απέχει πολύ από αυτόν που βγαίνει με τον μετασχηματισμό Tustin (s = 2/Ts*(1+z)/(1-z)), οπότε αναμένω no-overshoot εκδοχή του Zigler Nic να δουλέψει μια χαρά πάλι (δεν το δοκίμασα αλλά έβγαλα την εξ. διαφορών και είδα οτι μοιάζουν τρομερά)

Αν σου πηγαίνει μεταβίας μέχρι τους 85C, αύξησε το Y_MAX στα 20. Eγώ το έβαλα 15000 γιατί θεώρησα οτι μερικά Watt/sec αρκούν να αντισταθμίσουν τα 10m^3/sec Pου θα μπαίνουν, οπότε θεώρησα οτι η συνάρτηση μεταφοράς είναι σε Celsius/Watt.

Sυμπλήρωση για το PID (σύμφωνα με την δικιά μου προσέγγιση με matlab και ZN):
όχι απλώς το PID δούλεψε, αλλά ο χρονος αποκατάστασης βγήκε περίπου 4 sec όπως προέβλεψε και το matlab...

[daydreamer]

Τα κουμπιά STOP και EMPTY που με το που πατηθει το πρωτο πρεπει η διαδικασια να διακοπει σε οποιοδηποτε σταδιο κι αν βρισκεται και μετα και μονο μετα αν πατηθει το EMPTY να αδειασει το δοχειο πώς τα υλοποιούμε;

[princess_of_the_dawn]

Όταν βγαίνω από ένα στεπ δεν κρατιούνται οι τιμές τις οποίες όρισα μες στο στεπ. Έχει κανείς ιδέα γιατί συμβαίνει αυτό;

[Goemon]

Όταν βγαίνω από ένα στεπ δεν κρατιούνται οι τιμές τις οποίες όρισα μες στο στεπ. Έχει κανείς ιδέα γιατί συμβαίνει αυτό;

Τώρα θα σε γελάσω. Παρατήρησα 2 περιπτώσεις:

1) Όταν το action του step είναι απευθείας μεταβλητή, τότε αυτή "ενεργοποιείται" μόνο όσο ισχύει το στεπ. Αυτό είναι εκ κατασκευής έτσι. Αξιοποίησέ το

2) Όταν το action είναι κάτι σε function ή άλλο step program ή κάτι σύνθετο, παρατήρησα οτι δεν επανέρχονται οι μεταβλητές.

Αν είσαι μέσα σε μία από τις παραπάνω περιπτώσεις, δεν νομίζω πως έχει κάποιο λάθος... το λέω με 87.45% σιγουριά

[driverbulba]

Αν βαλεις σκετη μεταβλητη κι απο διπλα S0 αντι για Ν, η μεταβλητή σετάρεται true και δεν αλλαζει βγαινοντας απο το step. Για να την κανεις παλι false, σε αλλο step, βαζεις R0 διπλα. Το Ν κανει ο,τι κανει μονο οσο το step ειναι ενεργο.

[tniton]

Πληροφοριακα για τους επομενους οι ερωτησεις του εργαστηριου ειναι του στυλ,
-Ποιες μεταβλητες χρησιμοποιησατε για να προσομοιωσετε τα start, stop , τη σταθμη του υγρου κτλ στο visualization;
-Την εξισωση διαφορων που χρησιμοποιησαμε;
-Να τρεξουμε το visualization ως ενα σημειο, οπου θα παταμε το stop και θα καταγραψουμε τις καταστασεις των μεταβλητων(βαλβιδων θερμοκρασιας μιξερ) μεσω του visualization.
-Πως λειτουργει ο PID(αυτο νομιζω ειναι και το σος που ρωταει ολους) (+πως προσομοιωσαμε τη διαταραχη);

Αυτα περιπου, ο Χασαπης ειπε οτι θα τρεξει τα προγραμματα μας και θα δει αν επαληθευονται αυτα που γραψαμε...