Ισογεωμετρική ανάλυση είναι η πρόταση να λυθεί το πρόβλημα πάνω στις ίδιες συναρτήσεις που περιγράφουν ήδη τη γεωμετρία. Το μοντέλο CAD δεν μεταφράζεται σε πλέγμα πεπερασμένων στοιχείων· γίνεται το ίδιο το μοντέλο ανάλυσης. Διατυπώθηκε έτσι από τους Hughes, Cottrell και Bazilevs το 2005. Isogeometric analysis is the proposal that the problem be solved on the same functions that already describe the geometry. The CAD model is not translated into a finite element mesh; it becomes the analysis model itself. It was stated in this form by Hughes, Cottrell and Bazilevs in 2005.
Στην κλασική διαδρομή υπάρχουν δύο χώροι συναρτήσεων και ανάμεσά τους στέκεται το πλέγμα: ο ένας γράφει το σχήμα, ο άλλος το άγνωστο πεδίο, και η μετάφραση από τον πρώτο στον δεύτερο είναι χειρωνακτική, χρονοβόρα και — το κρίσιμο — φέρνει δικό της σφάλμα. Η ισογεωμετρική δεν βελτιώνει τη μετάφραση. Την καταργεί. The classical route has two function spaces with the mesh standing between them: one writes the shape, the other the unknown field, and the translation from the first to the second is manual, slow and — crucially — brings an error of its own. Isogeometric analysis does not improve that translation. It removes it.
Ένας κύκλος παραμένει κύκλος σε κάθε επίπεδο εκλέπτυνσης — και στο πιο αδρό. Δεν υπάρχει ακολουθία που να τον πλησιάζει, γιατί δεν χρειάστηκε ποτέ να τον προσεγγίσει κανείς. Το σφάλμα γεωμετρίας δεν μικραίνει: δεν υπάρχει. A circle stays a circle at every refinement level — including the coarsest. There is no sequence converging to it, because nobody ever had to approximate it. The geometry error does not shrink: there is none.
Τα όρια των στοιχείων είναι οι εσωτερικοί κόμβοι του διανύσματος κόμβων — γραμμένοι μέσα στο αρχείο CAD πριν καν τεθεί το ερώτημα της ανάλυσης. Δεν υπάρχει βήμα δημιουργίας πλέγματος να αποτύχει, ούτε παραμορφωμένο στοιχείο να διορθωθεί με το χέρι. Element boundaries are the interior knots of the knot vector — written inside the CAD file before the analysis question was even asked. There is no meshing step to fail, and no distorted element to fix by hand.
Ανάμεσα σε δύο κλασικά στοιχεία η λύση είναι συνεχής αλλά η κλίση της όχι. Εδώ η συνέχεια είναι Cp−1: στο κυβικό, δύο συνεχείς παράγωγοι. Γι’ αυτό μια θεωρία κελύφους που απαιτεί C¹ λύνεται χωρίς τεχνάσματα. Between two classical elements the solution is continuous but its slope is not. Here continuity is Cp−1: for cubics, two continuous derivatives. Which is why a shell theory that demands C¹ can be solved without tricks.
Οι δύο μεσαίες στήλες είναι το κέρδος και δεν είναι ανεξάρτητες μεταξύ τους: και οι δύο βγαίνουν από το ότι δεν άλλαξε χώρος συναρτήσεων. Το μηδέν της πρώτης το βλέπει κανείς μετρημένο απέναντι σε πλέγμα στη σελίδα CAE· η δεύτερη είναι ο λόγος που η μέθοδος έφτασε στα λεπτά κελύφη — και εκεί συναντά τις επιφάνειες υποδιαίρεσης, που έχουν την ίδια ιδιότητα από άλλον δρόμο. The two middle figures are the gain, and they are not independent of each other: both follow from the fact that the function space never changed. The zero in the first is shown measured against a mesh on the CAE page; the second is why the method reached thin shells — and there it meets subdivision surfaces, which have the same property by another route.
ΣχήμαShape · x(t) = Σ Nᵢ(t) · Pᵢ ΛύσηSolution · u(t) = Σ Nᵢ(t) · uᵢ
Η βάσηThe basis · Nᵢ(t) — ίδια και για τις δύο· ο διακόπτης δεν την αγγίζειthe same for both; the switch never touches it
Σύρε τις λαβές ή τον άξοναDrag the handles or the axis
συντελεστές · ό,τι γράφειςcoefficients · what you author ό,τι υπολογίζει η βάσηwhat the basis computes
Πάτησε Λύση και κοίτα τον κάτω καμβά: δεν κουνιέται. Άλλαξαν οι συντελεστές, όχι οι συναρτήσεις — και αυτό είναι όλη η ισογεωμετρική ανάλυση σε μία κίνηση. Πάτησε μετά p = 1: οι συναρτήσεις γίνονται τα «καπέλα» του κλασικού γραμμικού πεπερασμένου στοιχείου, η καμπύλη πέφτει πάνω στο πολύγωνο των συντελεστών, και η συνέχεια πέφτει από C² σε C⁰. Ίδια μηχανή, ένα κουμπί διαφορά. Press Solution and watch the lower canvas: it does not move. The coefficients changed, the functions did not — and that is the whole of isogeometric analysis in one gesture. Now press p = 1: the functions become the hat functions of the classical linear finite element, the curve collapses onto the coefficient polygon, and continuity drops from C² to C⁰. Same machine, one button apart.
Ως εδώ η βάση έκανε τη μία της δουλειά, το σχήμα. Η δεύτερη δουλειά δεν είναι καμπύλη — είναι πίνακας. Κάθε άγνωστος του προβλήματος αντιστοιχεί σε μία συνάρτηση, και κάθε κελί του πίνακα είναι το ολοκλήρωμα του γινομένου δύο συναρτήσεων. Ό,τι ισχύει για τους φορείς τους, ισχύει αυτούσιο για τα κελιά. So far the basis has done one of its two jobs, the shape. The second job is not a curve — it is a matrix. Every unknown of the problem corresponds to one function, and every cell of the matrix is the integral of a product of two functions. Whatever holds for their supports holds verbatim for the cells.
Mij = ∫ Ni Nj · σκούρο = μεγάλο στοιχείο· λευκό = αυστηρό μηδέν, όχι μικρόdark = large entry; white = strict zero, not small
μη μηδενικό · οι δύο φορείς τέμνονταιnon-zero · the two supports overlap μηδέν · δεν συναντιούνται ποτέzero · they never meet
Mij = ∫01 Ni(t) · Nj(t) dt
Σύρε τον διακόπτη από τις 6 συναρτήσεις στις 40. Ο πίνακας μεγαλώνει, η ζώνη δεν παχαίνει: μένει 2p+1 κελιά, δηλαδή 7 στο κυβικό, ανεξάρτητα από το πόσο μεγάλο έγινε το πρόβλημα. Τα μη μηδενικά είναι n(2p+1) − p(p+1), οπότε η πληρότητα πέφτει σαν 1/n — από 50% στις 12 συναρτήσεις σε 17% στις 40, και συνεχίζει. Drag the slider from 6 functions to 40. The matrix grows, the band does not thicken: it stays 2p+1 cells, that is 7 for cubics, no matter how large the problem became. The non-zeros number n(2p+1) − p(p+1), so the fill falls like 1/n — from 50% at 12 functions to 17% at 40, and onwards.
Η αραιότητα δεν επιβλήθηκε πουθενά. Είναι ο τοπικός φορέας του πάνω διαγράμματος, γραμμένος ως πίνακας: δύο συναρτήσεις που δεν μοιράζονται ούτε ένα διάστημα κόμβων δίνουν ολοκλήρωμα ακριβώς μηδέν, όχι μικρό. Η ίδια ιδιότητα που κάνει τη σχεδίαση να συγκλίνει, κάνει και το γραμμικό σύστημα λύσιμο. Sparsity was imposed nowhere. It is the local support of the diagram above, written as a matrix: two functions that share not a single knot span integrate to exactly zero, not to something small. The very property that makes design converge is the one that makes the linear system solvable.
Και τα δύο διαγράμματα στέκονται πάνω σε μία παραδοχή που δεν έχει ειπωθεί ρητά: ότι το αρχείο CAD κουβαλάει ήδη αυτή τη βάση. Την κουβαλάει — είναι τα NURBS, με το διάνυσμα κόμβων και τα βάρη τους. Εκεί όμως μπαίνει και η μία πραγματική περιπλοκή της μεθόδου: η ρητή βάση δεν ολοκληρώνεται ακριβώς από κανέναν πολυωνυμικό τετραγωνισμό. Both diagrams rest on an assumption that has not been said out loud: that the CAD file already carries this basis. It does — it is NURBS, with its knot vector and its weights. And that is also where the method’s one real complication enters: a rational basis is integrated exactly by no polynomial quadrature rule.
Όχι κλάδοι — αυτούς τους απαριθμεί η γενική σελίδα. Εδώ είναι τρεις περιπτώσεις όπου η διαφορά δεν είναι ταχύτητα αλλά το αν λύνεται καθόλου το πρόβλημα έτσι όπως γράφτηκε. Not industries — the overview page lists those. Here are three cases where the difference is not speed but whether the problem can be solved at all as it was written.
Η θεωρία Kirchhoff–Love απαιτεί C¹ από τη λύση, και τα κλασικά C⁰ στοιχεία δεν την έχουν — γι’ αυτό χρειάστηκαν μεικτές διατυπώσεις και πρόσθετοι άγνωστοι. Μια βάση με συνέχεια Cp−1 την έχει εξαρχής. Το 2000 οι Cirak, Ortiz και Schröder το έδειξαν με βάση υποδιαίρεσης, πριν καν ονομαστεί η μέθοδος. Kirchhoff–Love theory demands C¹ of the solution, and classical C⁰ elements do not have it — hence mixed formulations and extra unknowns. A basis with Cp−1 continuity has it from the start. In 2000 Cirak, Ortiz and Schröder showed this with a subdivision basis, before the method even had its name.
Σε ροή γύρω από καμπύλο σώμα, ή σε επαφή κυλίνδρου που κυλά, η διακριτοποιημένη επιφάνεια έχει γωνίες — και η λύση τις βλέπει ως γεωμετρικά χαρακτηριστικά, όχι ως σφάλμα. Το αποτέλεσμα είναι περιοδικός κυματισμός στην πίεση επαφής ή στο σημείο αποκόλλησης, που πυκνώνει με το πλέγμα αλλά δεν φεύγει. Με ακριβή γεωμετρία δεν εμφανίζεται καθόλου. In flow around a curved body, or in rolling contact, the discretised surface has corners — and the solution reads them as geometric features, not as error. The result is a periodic ripple in contact pressure or in the separation point, which gets finer with the mesh but does not go away. With exact geometry it never appears.
Οι μεταβλητές σχεδιασμού είναι οι ίδιοι οι συντελεστές που κρατάει το αρχείο CAD. Ο βελτιστοποιητής μετακινεί λαβές — όχι κόμβους πλέγματος — οπότε δεν χρειάζεται ούτε νέο πλέγμα σε κάθε επανάληψη ούτε αντίστροφη μετάφραση του αποτελέσματος πίσω σε γεωμετρία. Το βέλτιστο σχήμα βγαίνει ήδη ως μοντέλο CAD. The design variables are the very coefficients the CAD file holds. The optimiser moves handles — not mesh nodes — so there is neither a re-meshing step per iteration nor a translation of the result back into geometry. The optimal shape comes out as a CAD model already.
| ΒήμαStep | Τι αλλάζει — και πού το βλέπειςWhat changes — and where you see it |
|---|---|
| προετοιμασίαpreparation | Το βήμα δημιουργίας πλέγματος φεύγει: τα στοιχεία είναι οι εσωτερικοί κόμβοι — ο ορισμός.The meshing step disappears: the elements are the interior knots — the definition. |
| γεωμετρίαgeometry | Σφάλμα μηδέν σε κάθε επίπεδο, όχι μικρό σφάλμα — μετρημένο απέναντι σε πλέγμα στο CAE.Zero error at every level, not small error — measured against a mesh under CAE. |
| εκλέπτυνσηrefinement | Εισαγωγή κόμβου: νέοι άγνωστοι, ίδιο ακριβώς σχήμα — ο κανόνας του Böhm.Knot insertion: new unknowns, the exact same shape — Böhm’s rule. |
| σύστημαsystem | Ζώνη 2p+1 που δεν εξαρτάται από το μέγεθος — το δεύτερο διάγραμμα.A band of 2p+1 independent of size — the second diagram. |
| ολοκλήρωσηquadrature | Εδώ χρεώνεται: η ρητή βάση των NURBS δεν είναι πολυώνυμο και θέλει περισσότερα σημεία.Here it costs: the rational NURBS basis is not a polynomial and wants more points. |
| τριμάρισμαtrimming | Το ανοιχτό πρόβλημα: μια κομμένη επιφάνεια δεν έχει πια κανονικό πλέγμα κόμβων — CAD.The open problem: a trimmed surface no longer has a regular knot grid — CAD. |
Γιατί δεν είναι ακόμη η προεπιλογή. Δύο πράγματα κρατάνε τη μέθοδο πίσω, και κανένα από τα δύο δεν είναι η θεωρία. Το πρώτο είναι το τριμάρισμα: τα πραγματικά μοντέλα CAD είναι γεμάτα κομμένες επιφάνειες, και μια κομμένη επιφάνεια δεν έχει πια το κανονικό πλέγμα κόμβων που θέλει η μέθοδος. Το δεύτερο είναι το κόστος ολοκλήρωσης: μια ρητή συνάρτηση δεν ολοκληρώνεται ακριβώς από κανέναν πολυωνυμικό τετραγωνισμό, οπότε ό,τι κερδίζεται σε αγνώστους χρεώνεται εν μέρει σε σημεία Gauss. Το πρώτο είναι ενεργό ερευνητικό πεδίο· το δεύτερο είναι μετρήσιμος συμβιβασμός, και μετριέται. Why it is not the default yet. Two things hold the method back, and neither is the theory. The first is trimming: real CAD models are full of trimmed surfaces, and a trimmed surface no longer has the regular knot grid the method wants. The second is quadrature cost: a rational function is integrated exactly by no polynomial rule, so what is saved in unknowns is partly charged back in Gauss points. The first is an active research field; the second is a measurable trade, and it gets measured.
Πείτε μας τι φτιάχνετε και θα σας πούμε αν η υποδιαίρεση είναι το σωστό εργαλείο — και αν όχι, ποιο είναι.Tell us what you are building and we will tell you whether subdivision is the right tool — and if not, what is.