«CAE» είναι το βήμα όπου ρωτάς το μοντέλο κάτι φυσικό: θα αντέξει, πόσο θα ζεσταθεί, πώς θα ρεύσει, με ποια συχνότητα θα ταλαντωθεί. Οι εξισώσεις που απαντούν είναι συνεχείς και δεν λύνονται αναλυτικά πάνω σε πραγματικό σχήμα — άρα το σχήμα διακριτοποιείται: γίνεται πεπερασμένα στοιχεία με κόμβους, και η λύση γράφεται ως άθροισμα συναρτήσεων σχήματος πάνω τους. “CAE” is the step where you ask the model something physical: will it hold, how hot will it get, how will it flow, at what frequency will it vibrate. The equations that answer are continuous and cannot be solved analytically over a real shape — so the shape is discretised: it becomes finite elements with nodes, and the solution is written as a sum of shape functions over them.
Ανάμεσα στο μοντέλο του CAD και στον λύτη μεσολαβεί, λοιπόν, μια μετάφραση. Και αυτή η μετάφραση δεν είναι ουδέτερη: παράγει ένα σφάλμα που υπάρχει πριν γραφτεί η πρώτη εξίσωση ισορροπίας. Αυτό είναι το θέμα της σελίδας. Between the CAD model and the solver, then, sits a translation. And that translation is not neutral: it produces an error that exists before the first equilibrium equation is written. That is what this page is about.
Τάση, μετατόπιση, θερμοκρασία, ροή, ιδιοσυχνότητα. Διατυπώνεται πάνω σε ένα σώμα: μια περιοχή του χώρου με σύνορο, υλικό και φορτία. Το σύνορο είναι γεωμετρία — και έρχεται από το CAD. Stress, displacement, temperature, flow, natural frequency. It is posed over a body: a region of space with a boundary, a material and loads. The boundary is geometry — and it comes from CAD.
Το σώμα κόβεται σε στοιχεία και η άγνωστη συνάρτηση αντικαθίσταται από έναν πεπερασμένο συνδυασμό συναρτήσεων σχήματος. Στα κλασικά στοιχεία αυτές είναι πολυώνυμα ανά στοιχείο, συνεχή μόνο κατά C⁰ στα σύνορά τους. The body is cut into elements and the unknown function is replaced by a finite combination of shape functions. In classical elements these are per-element polynomials, continuous only to C⁰ across element boundaries.
Το ένα είναι της λύσης: το πολυώνυμο δεν είναι η αληθινή απόκριση. Το άλλο είναι της γεωμετρίας: το πλέγμα δεν είναι το σχήμα. Το πρώτο το παρακολουθούν όλοι· το δεύτερο σχεδόν κανείς, επειδή δεν το σχεδιάζει κανείς. One belongs to the solution: the polynomial is not the true response. The other belongs to the geometry: the mesh is not the shape. Everyone watches the first; almost nobody watches the second, because nobody draws it.
Το διάγραμμα από κάτω σχεδιάζει το δεύτερο. Δεν είναι σφάλμα του λύτη, δεν διορθώνεται με καλύτερο αλγόριθμο, και δεν εξαφανίζεται με πύκνωση — μόνο μικραίνει, με ρυθμό που μετριέται. The diagram below draws the second one. It is not the solver’s error, it is not fixed by a better algorithm, and it does not vanish under refinement — it only shrinks, at a rate that can be measured.
Η γεωμετρίαThe geometry · η ακριβής καμπύλη (μπλε) και το πλέγμα που πάει στον λύτη (κόκκινο)the exact curve (blue) and the mesh handed to the solver (red)
Το σφάλμαThe error · κατά μήκος του τόξου, με κορυφήalong the span, peaking at 1.90 mm · ο άξονας αυτοκλιμακώνεταιthe axis auto-scales
το πλέγμα · και το σφάλμα τουthe mesh · and its error η ακριβής γεωμετρίαthe exact geometry
Το σφάλμα αυτό δεν είναι σφάλμα επίλυσης. Υπάρχει πριν αρχίσει η επίλυση, το φτιάχνει η πλεγματοποίηση, και δεν το βλέπει σχεδόν ποτέ κανείς επειδή δεν σχεδιάζεται. Σύρε τα στοιχεία: πέφτει σαν h², δηλαδή διπλασιάζοντάς τα το σφάλμα τετραπλασιάζεται προς τα κάτω — αλλά δεν φτάνει ποτέ στο μηδέν, όσα στοιχεία κι αν βάλεις. Ο τρίτος μετρητής είναι μηδέν εξαρχής: ισογεωμετρικά δεν υπάρχει δεύτερη γεωμετρία να αποκλίνει, γιατί η ανάλυση τρέχει πάνω στην ίδια βάση που έγραψε το σχήμα. This error is not a solver error. It exists before the solve begins, meshing creates it, and almost nobody ever sees it because it is never drawn. Drag the elements: it falls as h², so doubling them divides it by four — but it never reaches zero, however many you add. The third counter is zero from the start: isogeometrically there is no second geometry to drift, because the analysis runs on the same basis that authored the shape.
Το γιατί πέφτει σαν h² είναι γεωμετρία, όχι φυσική. Ένα ευθύγραμμο στοιχείο είναι μια χορδή του πραγματικού τόξου: ακουμπά τους δικούς του κόμβους και αποκλίνει στο ενδιάμεσο, και το βέλος της χορδής μεγαλώνει με το τετράγωνο του μήκους της. Μισό μήκος, ένα τέταρτο σφάλμα — ακριβώς αυτό μετράει ο κάτω καμβάς, και ακριβώς αυτό δείχνει ο μετρητής. The reason it falls as h² is geometry, not physics. A straight element is a chord of the true arc: it touches its own nodes and departs in between, and the sagitta of a chord grows with the square of its length. Half the length, a quarter of the error — which is exactly what the lower canvas measures and the counter reports.
Στοιχεία ανώτερης τάξης με καμπύλα σύνορα το μικραίνουν — δεν το μηδενίζουν: το σύνορο παραμένει άλλη καμπύλη από εκείνη που έγραψε το CAD. Και το χειρότερο σημείο είναι το πιο ακριβό: εκεί όπου η καμπυλότητα είναι έντονη, δηλαδή σε φιλέτα και εγκοπές, είναι και εκεί που κορυφώνεται η τάση. Η ισογεωμετρική ανάλυση προτείνει την άλλη λύση: να μη γίνει καθόλου η μετάφραση, τρέχοντας την ανάλυση πάνω στην ίδια βάση που έγραψε το σχήμα — οπότε το σφάλμα γεωμετρίας δεν μικραίνει, απλώς δεν γεννιέται. Higher-order elements with curved edges shrink it — they do not zero it: the boundary remains a different curve from the one CAD authored. And the worst place is the most expensive one: where curvature is strong, that is at fillets and notches, is also where stress peaks. Isogeometric analysis proposes the other answer: do not perform the translation at all, running the analysis on the same basis that authored the shape — so the geometric error is not reduced, it is simply never created.
Τρεις περιπτώσεις όπου το σφάλμα της πλεγματοποίησης δεν είναι ακαδημαϊκή λεπτομέρεια αλλά αλλάζει την απάντηση ή το κόστος της: Three cases where the meshing error is not an academic detail but changes the answer, or what the answer costs:
Η συγκέντρωση τάσεων εξαρτάται από την καμπυλότητα του συνόρου. Ένα πλέγμα αντικαθιστά τη λεία καμπύλη με επίπεδες όψεις που συναντιούνται υπό γωνία: η κάθετος πηδά από όψη σε όψη, και η τάση που διαβάζεις στην ακμή ακολουθεί το πλέγμα και όχι το κομμάτι. Πυκνώνοντας συγκλίνει — αλλά πρέπει να ξέρεις ότι δεν είχε συγκλίνει ήδη. Stress concentration depends on the curvature of the boundary. A mesh replaces the smooth curve with flat facets meeting at angles: the normal jumps from facet to facet, and the stress you read at the edge follows the mesh rather than the part. Refining converges — but you have to know it had not already converged.
Κάθε φορά που αλλάζει το σχέδιο, το πλέγμα ξαναχτίζεται από την αρχή — και συχνά με το χέρι. Η πρόταση των Hughes, Cottrell & Bazilevs (2005) ήταν να καταργηθεί αυτό το βήμα: η ανάλυση να τρέχει πάνω στην ακριβή γεωμετρία του CAD, και η εκλέπτυνση για την ανάλυση να γίνεται χωρίς να αλλάζει η γεωμετρία ούτε η παραμετροποίησή της. Every time the design changes, the mesh is rebuilt from scratch — often by hand. The proposal of Hughes, Cottrell & Bazilevs (2005) was to abolish that step: run the analysis on the exact CAD geometry, and refine the basis for analysis without changing the geometry or its parameterisation.
Η θεωρία Kirchhoff–Love για λεπτά κελύφη ζητά πεδίο μετατοπίσεων C¹, ενώ τα κλασικά στοιχεία δίνουν C⁰. Οι Cirak, Ortiz & Schröder (2000) έδειξαν ότι η βάση της υποδιαίρεσης Loop είναι εγγενώς C¹ και άρα κατάλληλη — η ίδια αναπαράσταση που γράφει το σχήμα, να λύνει και το πρόβλημα. Kirchhoff–Love thin-shell theory asks for a C¹ displacement field, while classical elements deliver C⁰. Cirak, Ortiz & Schröder (2000) showed that the Loop subdivision basis is inherently C¹ and therefore suitable — the same representation that authors the shape solving the problem too.
| ΣφάλμαError | Πού γεννιέται — και τι το μικραίνειWhere it is born — and what reduces it |
|---|---|
| γεωμετρίαςgeometric | Στην πλεγματοποίηση, πριν από κάθε φυσική. Πυκνότερο πλέγμα ή καμπύλα στοιχεία — ή καμία μετάφραση, οπότε μηδέν.At meshing, before any physics. A denser mesh or curved elements — or no translation at all, hence zero. |
| διακριτοποίησηςdiscretisation | Στην επιλογή των συναρτήσεων σχήματος. Μικρότερα στοιχεία (h) ή ανώτερος βαθμός (p).In the choice of shape functions. Smaller elements (h) or higher degree (p). |
| μοντέλουmodelling | Στις παραδοχές: υλικό, φορτία, συνοριακές συνθήκες. Δεν το αγγίζει καμία πύκνωση.In the assumptions: material, loads, boundary conditions. No refinement touches it. |
| ανοχήςtolerance | Κληρονομημένο: το μοντέλο έφτασε ήδη «στεγανό μέσα σε ανοχή» — CAD.Inherited: the model already arrived “watertight within tolerance” — CAD. |
Και η ισογεωμετρική δεν είναι δωρεάν. Το ότι το σφάλμα γεωμετρίας μηδενίζεται δεν σημαίνει ότι μηδενίζονται τα υπόλοιπα: το σφάλμα διακριτοποίησης παραμένει, τα μητρώα είναι πυκνότερα επειδή ο φορέας κάθε συνάρτησης πιάνει περισσότερα στοιχεία, και ο λύτης πρέπει να δέχεται τη βάση του CAD αντί για τυπικά στοιχεία Lagrange — υποστήριξη που στα εμπορικά πακέτα είναι ακόμη περιορισμένη. Το επιχείρημα δεν είναι «λύθηκε»· είναι ότι ένα από τα τέσσερα σφάλματα του πίνακα δεν χρειάζεται καν να γεννηθεί. And isogeometric analysis is not free either. Zeroing the geometric error does not zero the rest: the discretisation error remains, the matrices are denser because each function’s support spans more elements, and the solver has to accept the CAD basis instead of standard Lagrange elements — support that is still limited in commercial packages. The claim is not “solved”; it is that one of the four errors in the table need never be born.
Πείτε μας τι φτιάχνετε και θα σας πούμε αν η υποδιαίρεση είναι το σωστό εργαλείο — και αν όχι, ποιο είναι.Tell us what you are building and we will tell you whether subdivision is the right tool — and if not, what is.