THMMY.gr

Tο παρόν topic αφορά την εργασία του μαθήματος κατά το χειμερινό εξάμηνο 2016-2017. Stay on topic!

Δε ξερω αν εχει ξεκινησει να ασχολειται κανεις, αλλα εχω μια απορια.. πως μετατρεπουμε μια διαφορικη εξισωση σε εξισωση διαφορων?

http://prnt.sc/d9moso (http://prnt.sc/d9moso)

Τσέκαρε ΨΕΣ εκεί έχει την θεωρία.
Γενικά αυτό που εφαρμόζω στα γρήγορα είναι dy/dt =x έστω δίνει y[n]-y[n-1]=x[n] αλλά δεν είναι και 100% σωστό απλά για να κάνω μια γρήγορη ανάλυση στο μυαλό μου.

επισης πως βρισκουμε τις παραμετρους που ζηταει  για το pid
?

Θα βρεις τη γραφικη παρασταση που απεικονιζει η συναρτηση μεταφορας και απο αυτη θα υπολογισεις τις απαραιτητες παραμετρους μεσω καπιου αλγοριθμου (π.χ. nichols ziegler  https://books.google.gr/books?id=2WQP5JGaJOgC&pg=PA818&lpg=PA818&dq=nichols+ziegler+dead+time+0.63&source=bl&ots=P5mOVhu8CL&sig=Ara0cLDeMUE2YlO1MjUVsF2Ufyo&hl=el&sa=X&ved=0ahUKEwj28aX_vbrQAhXHBBoKHdozCM8Q6AEIKjAC#v=onepage&q=nichols%20ziegler%20dead%20time%200.63&f=false (https://books.google.gr/books?id=2WQP5JGaJOgC&pg=PA818&lpg=PA818&dq=nichols+ziegler+dead+time+0.63&source=bl&ots=P5mOVhu8CL&sig=Ara0cLDeMUE2YlO1MjUVsF2Ufyo&hl=el&sa=X&ved=0ahUKEwj28aX_vbrQAhXHBBoKHdozCM8Q6AEIKjAC#v=onepage&q=nichols%20ziegler%20dead%20time%200.63&f=false))

Παιδιά ξέρει κανείς τι παίζει με την δεύτερη εργασία προγραμματισμός σε περιβάλλον 4DIAC που ανέβηκε στο ethmmy?

Έχει κανείς το email του μεταπτυχιακού που κάνει τα εργαστήρια για να στείλω mail?

Ρε παιδια να ρωτησω εγω κατι, λεει στο 6) στην εκφώνηση οτι με την νεα τιμη που θα υπολογίσουμε για την τροφοδοσια μεσω του PID( ας παραβλεψω οτι εχω καψει το κεφαλι μου στο πως θα φτασω εδω και 1 μερα εκει αλλα αστο αυτο, γιατι τα ΨΕΣ δεν πολυβοηθανε και ... :/ ) να παμε στην σχεση (1) και να βρουμε μια νεα τιμη για την θερμοκρασια.Η σχεση (1) υποτιθεται ειναι οταν δεν λειτουργει ο PID και για τιμη 15kW(μεγιστη τιμη). Έστω οτι εγω βρήκα τιμη 10kW , τι ακριβως θα το κανω αυτο το νουμερο στην σχέση (1)?

Ρώτησα τον Χασαπη και μου είχε απαντήσει αυτο :
"
Εφέτος δεν θα υπάρξει απαλλακτική εργασία, αλλά το 4DIAC μπορείτε να το χρησιμοποιήσετε γα να κατανοήσετε καλύτερα πως γίνεται η συγγραφή προγραμμάτων σύμφωνα με το IEC61499. Στη γραπτή εξέταση θα υπάρχουν ερωτήσεις από προγραμμ που θα έχει συνταχθεί και αποσφαλματωθεί με τη βοήθεια του 4DIAC.
"
Οποτε ουσιαστικα αν θες να κανεις την εργασία ειναι πιο πολυ για σενα γιατι θα σε βοηθησει στις εξετασεις

Αρα η προφορικη εξέταση είναι για την εργασια που θα κάνουμε στο codesys (αυτη με τον χημικο αντιδραστηρα)

Η διορία για την εργασία ποιά είναι;

Πως φαίνεται η επιμελής φοιτήτρια, παραμονή πρωτοχρονιάς κανοντας εργασία ΒΠ..

Μάλλον το Φλεβάρη που είναι και η εξέταση της εργασίας αλλά δεν είναι 100% σίγουρο, πρέπει να ρωτησουμρ το Χασάπη..  

χαχαχαχα σίίίγουραααα

Τελικα καταλαβες τι ακριβως θελει γιατι εγω απο την απαντηση που ανεβασε στο ετημμυ δεν...

Οχι δεν μπορω να πω οτι καταλαβα στα σιγουρα αλλα αν το κανουμε χειροκινητα και με τον Yokogawa για PID, ουσιαστικα πρεπει να κανεις manual tuning των τιμών Heater για ( t-1 ) και (t-2)καθώς και του Temperature(t-1) καθώς και το χ(κ),χ(κ-1) το οποιο δεν ειναι καθολου εύκολο(δοκιμαζα πολλές ώρες ακομα ακρη δεν βρηκα). Πήγα να χρησιμοποιήσω τώρα τον ετοιμο PID του Codesys αλλα παλι το output του ετοιμου block ( δηλαδη την νεα τιμή της ισχύς του κινητήρα) πρέπει να την χρησιμοποιήσεις μέσα απο μια σχέση για να βρεις καινουργια τιμή θερμοκρασίας. Δοκίμασα εξίσωση διαφορών και την συνάρτηση μεταφοράς που είπε ο χασάπης στο ethmmy παλι δεν τρέχει σωστά.ΕΠΙΣΗΣ ΞΕΡΕΙ ΚΑΝΕΙΣ ΠΩΣ ΧΡΗΣΙΜΟΠΟΙΕΙΣ ΣΑΝ ΜΕΤΑΒΛΗΤΗ ΤΟΝ ΧΡΟΝΟ ?δηλαδη το Τ(και δεν εννοω να βάλεις πχ το main task να τρέχει ανα second και ενα έχεις ενα μετρητή που αυξάνει).
Βοηθήστε :)

μια απ'τα ίδια... Οσον αφορα τη μεταβλητη στο χρονο στο τελευταιο εργαστηριο απλα πολλαπλασιαζει το τρεχων k με τον χρονο δειγματοληψιας.

κοιτα έχω μια εντύπωση ότι αν διαβασουμε τις απαντήσεις απο χρονιές 2011-2012 και 2012-2013 που απο ότι φαινεται είναι η ίδια εργασία ίσως βγάλουμε μια άκρη παραπάνω, αλλα για τον χρόνο βρήκα πατέντα και βάζω TIME_TO_REAL(Step#.t)/1000 όπου # το Step που βρίσκεσαι και το /1000 γιατι το γυρναει σε ms!!!

Που μπορουμε να βρουμε αυτες τις απαντήσεις ;

Δεν θέλει καμιά αναφορά ή κάτι τέτοιο, ε; Μόνο τα αρχεία του project.

Σελίδα 9 και 11 στο Protipa_last.ppt λέει πως πας από την transfer function σε εξίσωση διαφορών.
Μετά Σελίδα 21 μέθοδος Ziegler Nichols ή σελίδες 22, 23, 25, 26 και 29 (PID) μέθοδος Persson και Astrom για υπολογισμό κερδών ελεγκτή.

_{Δεν βάζω και το χέρι μου στη φωτιά αλλά λογικά αυτά θέλει,}

Παιδιά τον αναδευτήρα πάνω λέει να τον ανοίξουμε 5sec (στην αρχική περιγραφή) και κάτω 5min (στο βήμα 9).
Τυπογραφικό είναι φαντάζομαι.
Ρώτησε κανείς ποιο είναι το σωστό;

Επίσης τις μεταβλητές που χρειαζόμαστε συχνά μπορούμε να τις βάλουμε Global;
Γενικά μας απαγορεύει το πρότυπο να έχουμε Global variables;

Προχώρησε κανένας με τον PID? Δώστε κανα tip ρε παιδιάααααα  :( :-\

Δεν το έχω υλοποιήσει ακόμα για να ξέρω ότι δουλεύει αλλά η G(s) που έχει στην διαφάνεια είναι η TEMP(s)/HEATER(s). Από εκεί παίρνεις τα α, β, γ για να υπολογίσεις την τιμή που πήρε η θερμοκρασία στην διάρκεια του κύκλου που πέρασε. Επίσης αυτήν την τιμή την αφαιρείς από 100 για να βρεις το σφάλμα. Επίσης χρειάζεται να βρεις τα κέρδη του ελεγκτή με έναν από τους 2 τρόπους (Ziegler Nichols ή Persson Astrom) ώστε να πολλαπλασιάσεις το παραπάνω σφάλμα (όρος P), το ολοκλήρωμα του σφάλματος (όρος I) και την παράγωγο του σφάλματος (όρος D) για να υπολογίσεις την έξοδο του ελεγκτή u(k), δηλαδή την ισχύ που θα δώσεις στον θερμαντή. Το ολοκλήρωμα φαντάζομαι μπορείς να το υπολογίζεις ως (προηγούμενη_τιμή + e(k)*ΔΤ) ενώ την παράγωγο ως (e(k)-e(k-1))/ΔΤ.

edit: Δεν δουλεύει :-\

Λοιπόν, νομίζω κάτι πάει λάθος με την εξίσωση που μας έδωσε. Βάζοντας την εξίσωση G στη matlab και τρέχοντας lsim για είσοδο heater = 15 έχει μέγιστη τιμή 84. Δηλαδή ο θερμαντής δεν έχει την απαραίτητη μέγιστη ισχύ για να κρατάει τον αντιδραστήρα στους 100 ^οC σύμφωνα με την G.

Για τον PID πρέπει να χουμε 2 συναρτήσεις, 1 που να υπολογίζει την θερμοκρασία με βάση την τάση που έχει βάλει ο PID και μία για τον PID που ν αλλάζει την καινούρια τάση ή με μία συνάρτηση του PID θα γίνονται όλα αυτά?

Εξέτααη εργασιών.
13 Φεβ 2017 3:16 μμ
Χασάπης

Από την τέταρτη και πέμπτη ομάδα διαγράφηκαν οι υπεράριθμοι φοιτητές. Παρακαλούνται οι φοιτητές αυτοί και όσοι δεν έχουν δηλώσει ακόμα ομάδα, να εγγραφούν σε μια από τις δυο νέες ομάδες που δημιουργήθηκαν(δείτε σχετικά στην ενότητα "Εργαστήριο".)

Παιδια εχω μια απορια.

Στις 23 και 24 Φερβουαριου η σχολη ειναι κλειστη λογω αργιας του πανεπιστημιου.

Εμεις πως σκατα θα εξεταστουμε στην εργασια?  :D

Στο ethmmy τις έχει, στο συχνές ερωτήσεις

Παίζει  κανείς να ανεβάσει πως έχει στήσει την βασική process? Για να πάρουμε μια ιδέα του πως να αρχίσουμε!

έχει κρατήσει κανείς τις ασκήσεις που κάναμε στο εργαστήριο 5; Αν ναι μπορεί να την ανεβάσει;  ::)

Γιατί στο ethmmy δεν τις έχουν ανεβάσει, βασικά τις έχουν ανεβάσει αλλά είναι αυτή που συμπληρώναμε πάνω! Έπρεπε να την αποθήκευα!   :-\

νομιζω ειναι ετσι

Ευχαριστω!! :)

Εχει καταλαβει κανεις πως θα κανουμε τον pid?Απλα θα υπολογιζει ποσα kW πρεπει να δωσει στον heater για να φτασει τους 100C η θα εχουμε και διαταραχη και θα πρεπει να προσεγγιζει παλι τους 100C???

+1 αν καταλαβε καποιος να βοηθησει λιγοο

Αρχικά η Συνάρτηση Μεταφοράς που έδωσε ο Χασάπης είναι από ισχύ σε θερμοκρασία αλλά όχι 15 kW απλά 15 W. Στην συνέχεια εισάγοντας υγρό λειτουργεί ο pid  με αρχική τιμή στα 15W (η προσέγγιση 3 όρων της ΣΜ βγαίνει στα 87C μπορείς να δώσεις παραπάνω ισχύ για να έχεις 100C σαν πρώτη τιμή) και μετά αλλάζει την τιμή της ισχύς βάσει του νόμου ελέγχου που έχεις φτιάξει. Μπορείς να πάρεις του εργαστηρίου τέσσερα, εάν θυμάμαι, καλά τον ελεγκτή και να του δώσεις τα νέα κέρδη που θα υπολογίσεις και την διακριτή προσέγγιση της θερμοκρασίας ή την συνεχή ΣΜ με κάποιες μετατροπές. Τέλος δεν βάζεις όρια στις τιμές εξόδου του ελεγκτή (Ισχύς που δίνεις στην ΣΜ) καθώς το δοκίμασα και στην matlab και είναι ασταθές το σύστημα με PID και όρια πχ [0,15] W.

η εξεταση της εργασιας ξερει κανεις τι φαση ειναι

Στιγμιαία πετάει σε κανέναν αρνητική θερμοκρασία ο PID???

Αυτό ίσα που φαίνεται και μετά δουλεύει κανονικά.


edit: Αλλάζοντας το T πηγαίνει καλύτερα πάντως. Αλλά μόνο με δοκιμές το βρήκα αυτό.

παιδια στην εξίσωση δαφορών οι συντελεστές πόσο σας βγήκαν?γιατί σε μένα βγήκαν α=0.2396
β=-5.8392 και γ=0.685,οπότε για αρχική ισχύ 15Watt μου βγαίνει θερμοκρασία -87 βαθμοι περίπου.κάπου είδα ότι την βγάζει +87.κοίταξα αλλά δεν φαίνεται να έχω κάνει λάθη στα πρόσημα με βάση τις σημειώσεις του.ξέρει κανείς τι παίζει εδώ?

παιζει να μπορει καποιος να μου πει πως υλοποιουμε τον pid βημα-βημα και βασικα πως να υπολογισω τις απαραίτητες εξισωσεις? δεν το εχω νοιώσει καθόλου.. :/ 

και εγω αυτα βρισκω αλλα μετα δεν ξερω να προσχωρήσω.
η εξισωση αυτη μας δινει την θερμοκρασία συναρτήσει της ισχύς, σωστα? δλδ  Temperature(t) = 0.24*Temperature(t-1) - 5,84*Heater(t-1) + 0,64*Heater(t-2)
βρηκα και τα Ti, Td, Kp, K χρησιμοποιώντας Astrom..  
μετα απο εδω τι κανω?

γιατί παίρνεις 5.84 θετικό?αφού αρνητικό το δεν το βρίσκεις και εσύ?νομίζω πρέπει σε κάθε εκτέλεση του προγράμματος σου(δείγμα χρόνου ας πούμε τα 20 ms) να εκτελείς τον αλγόριθμο Pid να βρίσκεις τη τιμή της ισχυός,μετά την αντίστοιχη τιμή της θερμοκρασίας και απλά να αφαιρείς από τη θερμοκρασία που βρίσκεις μία σταθερή τιμή πχ. -3.(αυτό το λέει στο ethmmy) στις συχνές ερωτήσεις και το κάνει για να δείξει πως έχουμε και επιδραση εξωτερικών διαταραχών.Αυτό όλο το κάνεις μέχρι το υγρό να φτάσει στο σημείο του LLS1.αυτή είναι η γενική ιδέα.το έχεις?

αρνητικο ειναι το διορθωσα!
ναι νομιζω καταλαβα.. το πολυ πολυ αν δεν μου βγει θα ξαναποσταρω  :D
thaaanks!

σας έτρεξε με αρνητικό β και θετικό γ?

Ρε παιδιά στο visualization έχω κολλήσει.
Δεν μπορώ να δημιουργήσω Image pool για να βάλω δικές μου εικόνες.
Αν κάνω create global image pool δεν μου ανοίγει τίποτα.

Κάτω αριστερά θα έχει Devices και POUs. Άνοιξε την καρτέλα POU. Εναλλακτικά πήγαινε View -> POU. Αν έκανες create global pool θα είναι εκεί.

Για τον PID πήρατε τις τιμές από τον πίνακα τριών όρων για Μ = 2?

Μπορει καποιος να εξηγησει τι κανουμε για τον PID σε σχέση με το παράδειγμα του τελευταίου εργαστηρίου?

Για να υλοποιησω τον pid στο κωδικα κανω το εξης.
Εισαγω μια μειωση θερμοκρασιας πχ 3 βαθμους λογω της ροης του υγρου Α. Αρα δημιουργω ενα σφαλμα. Αμεσως μετα υπολογιζω συμφωνα με τις εξισωσης απο το Yokogava την ισχυ που θα παρει ο θερμαντηρας. Στην προκειμενη περιπτωση η ισχυς ειναι η μεταβλητη p(kT)??
Τελος απο την εξισωση διαφορων υπολογιζω την νεα θερμοκρασια??? Ειναι σωστο ολο αυτοή δεν εχω καταλαβει τπτ?
Επισης το Τ στις εξισωσεις yokogava ποσο το οριζουμε και γιατι?
Στο φυλλαδιο στο βημα 6 γιατι λεει να υπολογισουμε την τιμη της θερμοκρασιας απο την εξισωση 1? Αν υπολογιστει απο εκει ή εξισωση διαφορων που φτιαξαμε που χρησιμευει?

Καταφερατε να βαλετε trace? Γιατι εμενα μου βγαζει σφαλμα "The device does not support the trace component in the runtime system"

χρειαζεται να βαλουμε Trace? αν ναι, που?

Καταφέρα να κάνω trace κανονικά. Ποιο device έχεις? Νομίζω το "CODESYS Control Win V3" διατίθεται μόνο στο students edition και αυτό έχουμε και στο εργαστήρι

Για τον PID το αξιοποίησα στην επαλήθευση του ελεγκτή... σκοπεύω να το έχω στο Project στην εξέταση για να δείξω ότι το σύστημα είναι κομπλέ....

δεν ειναι ομως απαραίτητο ε? επειδη στην εκφωνηση δεν ειδα να λεει καπου να το κανουμε

Στο γ)7) λέει βασικά "αρχιστε να δείχνετε πως μεταβάλλεται η θερμοκρασία ανά μονάδα χρόνου". Υποθέτει κανείς οτι μάλλον θέλει διάγραμμα εδω πέρα ...

Αλλά όχι δεν το λέει πουθενά ξεκάθαρα. Τσέκαρε αριστερά στην οθόνη σου, που είναι το tree με τα POUs. Η μορφή του θα είναι:

<YourProjectName>
   -Device (CODESYS Control Win V3)
      -Plc Logic
         -Application
            -κλπ κλπ

Αν το Device ΔΕΝ ειναι το παραπάνω, το σφάλμα σου είναι το αυτό

Μπορεί να πει κάποιος πως ήταν η εξέταση σήμερα? Και αν πρέπει να κάνουμε και trace για τον PID...

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

   
ΠΡΟΣΟΧΗ_ΑΝΑΚΟΙΝΩΣΗ
20 Φεβ 2017 11:48 πμ
Χασάπης

Η  εξέταση  της  εργαστηριακής  άσκησης  γίνεται  στο  εργαστήριο  Αρχιτεκτονικής  Υπολογιστικών
Συστημάτων  που  βρίσκεται  στο  ισόγειο  του  κτιρίου  Ε  της  Πολυτεχνικής(απέναντι  από  τη  νησίδα
Βεργίνα).  Οι  φοιτητές  να  έχουν  μαζί  τους  στυλό  διαρκείας  για  να  απαντήσουν  γραπτώς  στις
ερωτήσεις  που  θα  τεθούν.

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

Super SOS παιδια...

Αυτο με το στυλο με αγχωσε...

Μας έδωσε ένα φυλλάδιο να συμπληρώσουμε με ποιες μεταβλητές γίνεται το start και ανοίγουν οι βαλβίδες, μετά έπρεπε να εξηγήσεις τι γίνεται στο μιμικό διάγραμμα (ποιες μεταβλητές γίνονται true, false αντίστοιχα στον κώδικα σου κλπ) και στο τέλος είχε να εξηγήσεις και πώς λειτουργεί ο PID σου. Τίποτε άλλο. Του το έδινες στο τέλος και ανέβαζες την εργασία σου στο ethmmy.

Παιδια, στον PID  τις αρχικες τιμες των Bk,B(k-1),B(k-2), xk,x(k-1),x(k-2) pk, p(k-1),p(k-2), πως τις επιλέγουμε; Γιατι οταν τις βαζω 0 δε συγκλίνει. Παίζει μηπως ρολο που στην εξίσωση διαφορών βρήκα αρνητικό
 β(= -5.839);

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

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

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

• Έφτιαξα το μοντέλο στο 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ο εργαστηριο δεν μπορω να καταλαβω πως προκύπτουν οι αρχικες τιμες  που εχει και αρα τι πρεπει να αρχικοποιήσουμε εμεις ..

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

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

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

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

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

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

δηλαδη βαζοντας τις παραμετρους που βρηκαμε απο τη μεθοδο 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.

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

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

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

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

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

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

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

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

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

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

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

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

Διαταραχή? Ποια διαταραχή?

Τη σταθερη διαταραχη που αναφερει και στις συχνες ερωτησεις, τουλαχιστον αυτο καταλαβα...

Όταν η Blinker στον PID δηλώνεται τύπου BLINK βγάζει error ότι τοBLINΚ είναι αγνωστος τύπος. Εχει κανείς καμιά ιδέα πώς διορθώνεται αυτό;

έχεις την βιβλιοθήκη ; util ;

Το εργαστήριο και την εργασία μπορούμε να τα δώσουμε και τον Σεπτέμβρη?

Εχετε καταλαβει για ποιον λογο οι περισσοτεροι δεν εχουμε καν βαθμο στην εργασια?

Του στειλα μαιλ. περιμενουμε....

+1

+1
πραγματικα ελπιζω να ειναι απο λαθος...
Ας δουμε τι θα πει στον sunrise και πραττουμε αναλογως...

πάντως το κατέβασε από το ethmmy το βαθμολόγιο

Μαλλον ειχε γινει καποιο λαθος...

Μα να μην παρω ουτε μηδεν?  :D

Τώρα που θα το ξανανεβάσει, απλά θα προστεθεί η εργασία έτσι; Μην κοπώ έτσι άδικα από κει που είχα περάσει ;D

Θελω να πιστευω πως οι βαθμοι γραπτου οι ιδιοι θα μεινουν... Μονο θα προστεθουν οι εργασιες. Δεν ειαι προετοιμασμενος για ψυχρολουσια..  :P

Να φαμε και μαγειρεμα προς τα κατω δηλαδη?  Δεν τραβηξα και printscreen το πρωτο βαθμολογιο ρε γμτ  :(

Προσωπικά όπως έκανα το βαθμό μου βγαίνει ο ίδιος με το 50% να πιάνει η εξέταση 25% τα εργαστήρια και 25% η εργασία(που δεν είχα βαθμό).  Οπότε λογικά απλά προστίθεται. Επίσης αυτό που το 5.3 πχ μου το κάνουν 5 αλλά το 5.7 5.5 αντί για 6 ποτέ δεν το κατάλαβα. Μου έχουν κάνει και το 6,8 6,5 στο παρελθόν :D

Ναι κι εγω το υπολογισα ειναι 50% βαθμος γραπτου+25% μο εργαστηριων+ 0 απ ' την εργασια που δεν ειχε καν βαθμο. Οσο για τις στρογγυλοποιησεις μην το ψαχνεις καν θα ελεγα... Μονο θα συγχιστεις  :P

Γνωρίζει κανεις τελικά?