Sono richiesti 8 punti Gauss per gli elementi finiti esaedrici del secondo ordine?


10

È possibile ottenere la precisione del secondo ordine per gli elementi finiti esaedrici con meno di 8 punti Gauss senza introdurre modalità non fisiche? Un singolo punto Gauss centrale introduce una modalità di taglio non fisico e la disposizione simmetrica standard di 8 punti Gauss è costosa rispetto alle discretizzazioni tetraedriche.

Modifica : qualcuno ha chiesto equazioni. Le equazioni che mi interessano sono l'elasticità non lineare, dinamica o quasistatica. Le equazioni quasistatiche sono

P(ϕ)=0

dove , e è una prima funzione di stress iperelastica di Piola-Kirchoff. Un semplice esempio è comprimibile neo-hookean, dove ϕ:ΩR3ΩR3P:R3×3R3×3

P(F)=μ(FFT)+λFTlogdetF

Cosa stai simulando esattamente?
Dan

Elasticità lineare al momento, ma la domanda riguarda l'elasticità non lineare in generale.
Geoffrey Irving il

1
Probabilmente dovresti includere le equazioni che ti interessano, poiché la definizione di "non fisico" dipende da esse. O almeno definire con precisione lo spazio delle funzioni "fisiche".
David Ketcheson,

Equazioni aggiunte.
Geoffrey Irving il

1
Con dPhi / dx, intendi il gradiente?
Wolfgang Bangerth,

Risposte:



0

È relativamente ovvio che in generale non è possibile cavarsela con meno punti di quadratura per cella rispetto ai gradi di libertà. Nel caso di elementi trilineari su un esaedro 3d, ci sono 8 gradi di libertà (uno per vertice), quindi anche il numero minimo di punti di quadratura sarebbe otto.

che non è invertibile e di conseguenza completamente inutile. Il motivo è che una formula di quadratura a un punto non può distinguere tra tutte le funzioni lineari (parte dello spazio di prova) che hanno lo stesso valore nel punto di quadratura; in altre parole, per la regola del punto medio, la funzione di forma 'x' è la stessa della funzione '0' è la stessa della funzione '-x'. In altre parole, mentre lo spazio di prova ha dimensione 2 con integrali esatti, per la regola del punto medio lo spazio ha dimensione 1, anche se ci sono due gradi di libertà - questa è la definizione di uno spazio che non è insolente.) per la regola del punto medio, la funzione forma 'x' è uguale alla funzione '0' è uguale alla funzione '-x'. In altre parole, mentre lo spazio di prova ha dimensione 2 con integrali esatti, per la regola del punto medio lo spazio ha dimensione 1, anche se ci sono due gradi di libertà - questa è la definizione di uno spazio che non è insolente.) per la regola del punto medio, la funzione forma 'x' è uguale alla funzione '0' è uguale alla funzione '-x'. In altre parole, mentre lo spazio di prova ha dimensione 2 con integrali esatti, per la regola del punto medio lo spazio ha dimensione 1, anche se ci sono due gradi di libertà - questa è la definizione di uno spazio che non è insolente.)


Penso che la domanda di Geoff sia più sottile. Per spazi di elementi finiti continui su tetraedri in domini ben formati (ad es. Senza elementi isolati), è possibile cavarsela con quadrature a punto singolo che sono chiaramente sottointegrazione. La domanda è se lì è anche possibile sottointegrare in qualche modo con elementi hexahedral. Non conosco la risposta, ma non sono sicuro di quanto sia grande dato che i punti di quadratura non richiedono un movimento di memoria aggiuntivo. Dopo aver vettorializzato la valutazione residua degli elementi finiti, è comune che sia associato alla memoria, quindi potresti stare meglio usando i flop.
Jed Brown,

Un buon punto sul movimento della memoria.
Geoffrey Irving il

1
Espandere sul punto di Jed: la ragione per cui l'argomento "ovvio" sopra è falso è che ogni punto di quadratura vede una matrice . Per il tetraedro, che copre tutti i movimenti dei vertici esclusa la traduzione uniforme, che non influenza l'energia o le forze, quindi un punto di quadratura è sufficiente per la precisione del primo ordine. 3×3
Geoffrey Irving il

1
Piuttosto scomodo che i commenti non possano includere nuove righe.
Geoffrey Irving il

@JedBrown: buon punto. Il gradiente delle funzioni lineari sulle tet sono costanti e quindi è sufficiente un singolo punto di quadratura, seguendo l'argomento che ho fatto per la matrice di massa (la matrice di rigidità è la matrice di massa per i gradienti :-). D'altra parte, i gradienti delle funzioni trilineari sull'esahedra sono funzioni quadratiche (anisotropiche), quindi è necessario sicuramente più di un solo punto di quadratura per direzione delle coordinate.
Wolfgang Bangerth,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.