Berstmembran prismatischer Zellen: Berstdruck und Reststegdicke

Die Membran ist kreisförmig, anfangs eben und spannungsfrei, mit Radius $a$, voller Dicke $t$ und festgehaltener radialer Verschiebung am Rand. In Anlehnung an Fichter, NASA TP-3658 (1997), pp. 2–7 (attached pages 4–9), werden $x=r/a$, $q=pa/(Et)$ und $N_r/(Et)=q^{2/3}b(x)/4$ definiert. Die Hencky-Gleichung lautet $b^2(b''+3b'/x)=-8$. Ihre reguläre Reihe mit geraden Potenzen wird durch Koeffizientenfaltung bis $x^{48}$ erzeugt. Eine Nullstellenbestimmung mit dem Brent-Verfahren erzwingt $b'(1)+(1-nu)b(1)=0$. Die Membrankraft in Umfangsrichtung folgt aus $N_theta=d(rN_r)/dr$; die Integration von $W'=-2q^{1/3}x/b$ liefert die Durchbiegung.

An der umlaufenden Laserkerbe gilt $sigma_{rm score}=K_n N_r(r_s)/t_r$. Gleichsetzen mit der eingegebenen Bruchfestigkeit ergibt $p_b=(Et/a)[S_u t_r/(K_n Et C_r)]^{3/2}$ mit $C_r=b(r_s/a)/4$. Die inverse Berechnung löst dasselbe Kriterium nach der Reststegdicke auf und gibt die abgetragene Tiefe $t-t_r$ aus. Dieses Nettoquerschnittsmodell für eine schmale Kerbnut ist eine ingenieurmäßige Erweiterung und keine veröffentlichte Korrelation für das Bersten von Batteriezellen. Der effektive Kerbfaktor muss kalibriert werden.

Für unabhängige logarithmisch normalverteilte Grundgesamtheiten von Reststegdicke, Bruchfestigkeit und Kerbfaktor werden die Varianzen der logarithmierten Größen $v=ln(1+CV^2)$ verwendet. Dann gelten $s_{ln p}=1.5sqrt{v_t+v_S+v_K}$ und $p_{50}=p_bexp[-0.75(v_t+v_S-v_K)]$. Das 5. und 95. Perzentil ergeben sich aus $p_{50}exp(mp1.64485s_{ln p})$.

Während des Impulses beträgt die Gasbildungsrate $dot m_g=dot m_{pk}sin^2(pi t/t_g)$. Ihr analytisches Integral bestimmt den Zeitpunkt des erstmaligen Berstens. Der freie Öffnungsquerschnitt wächst gemäß $A(t)=A_osin^2[pi,mathrm{clamp}((t-t_b)/t_o,0,1)/2]$. Verwendet wird die isentrope Massenstromdichtegleichung von NASA Glenn. Für kritische Strömung wird das Verhältnis von Gegen- zu Vordruck nach unten auf $(2/(kappa+1))^{kappa/(kappa-1)}$ begrenzt; unterhalb des Gegendrucks ist der Gasaustritt null. Dieselbe Massenbilanz $dP/dt=RT(dot m_g-dot m_o)/V$ wird mit dem impliziten BDF2-Verfahren bei variabler Schrittweite und vier internen Schritten je Ausgabeintervall integriert. Die Phase mit geschlossener Zelle wird analytisch ausgewertet. Bei abrupten Schrittweitenänderungen oder negativen extrapolierten Verlaufswerten wird mit dem impliziten Euler-Verfahren neu gestartet. Für den Überdruck $y=P-P_b$ löst jeder implizite Schritt $alpha y+beta F(P_b+y)=B$. Der Zweig für kritische Strömung ist linear. Im unterkritischen Zweig werden $H=sqrt{2kappa/[(kappa-1)RT]}$, $f=F/(P_bH)$ und $v=(1+sqrt{1+4f^2})/2$ gesetzt; dann gilt $y=P_b[v^{kappa/(kappa-1)}-1]$. Die Newton-Iteration mit der nichtnegativen Massenstromdichte als Variable vermeidet die unendliche Druckableitung bei Überdruck null. Eine numerisch stabile Auswertung mit log1p/expm1 erhält die Genauigkeit nahe dem Umgebungsdruck. Dies ist eine implizite Integration der vollständigen Strömungsgleichung und kein quasistationäres Ersatzmodell.

In die Strömungsberechnung gehen Absolutdrücke ein; angezeigt werden Überdrücke. Die 1601 Ausgabezeitpunkte sind zu Beginn des Öffnens und der Entspannungsphase quadratisch verdichtet.

Homogene, isotrope, elastische Folie; keine Vorspannung, Biegesteifigkeit oder anfängliche Wölbung. Eine schmale umlaufende Laserkerbe verändert die globale Membransteifigkeit nicht. Radialzug öffnet die Laserkerbe; örtliche Plastizität und wärmebedingte Schädigung werden ausschließlich durch die experimentell bestimmte Bruchfestigkeit und den effektiven Kerbfaktor erfasst. Nach dem Aufreißen bleibt eine angebundene Folienklappe zurück; der endgültige freie Öffnungsanteil ist eine unabhängige Eingabegröße.

Reststegdicke, Bruchfestigkeit und Kerbfaktor sind unabhängig logarithmisch normalverteilt und werden durch arithmetische Mittelwerte und Variationskoeffizienten beschrieben. Es erfolgt keine erneute Ziehung von Dickenwerten unter Einhaltung einer harten Obergrenze bei der vollen Foliendicke. Andere Abmessungen, der Elastizitätsmodul und die Gaseingaben weisen in der instationären Berechnung keine Streuung auf.

Die Zelle ist ein starres, vollständig durchmischtes Idealgasreservoir mit der eingegebenen konstanten Temperatur und anfangs auf Gegendruck. Gaszusammensetzung, Ausflussbeiwert und Gegendruck sind konstant. Die Gastemperatur kann von der Folientemperatur abweichen. Die Voreinstellungen für Gasbildung, Öffnungszeit, freien Öffnungsanteil, Kerbfaktor und Streuung sind beispielhafte Szenariowerte und keine Messwerte einer kommerziellen Zelle.

Verwenden Sie dieses Modell für eine kreisförmige lasergekerbte Folie im Deckel einer prismatischen Zelle. Eine ovale oder rechteckige Öffnung wird damit nicht auf einen gleichwertigen Kreis zurückgeführt. Hencky beschreibt eine gleichmäßige Querbelastung bei großer Durchbiegung und kleiner Dehnung; die radiale Komponente einer mit der verformten Oberfläche mitgehenden Drucklast wird vernachlässigt. Ein Lastparameter über 0,01, ein Verhältnis von voller Foliendicke zu Radius über 0,05, ein Verhältnis von Durchbiegung zu Foliendicke unter 3 oder eine Spannung der ungeschwächten Folie oberhalb der Streckgrenze lösen eine Meldung zum Gültigkeitsbereich aus. Diese Kriterien werden sowohl beim Sollberstdruck als auch beim berechneten Berstdruck geprüft. Ergebniskarten zur Bemessung auf den Sollberstdruck werden ausgeblendet, wenn dieser außerhalb des Gültigkeitsbereichs liegt. Der Grenzwert für die Durchbiegung ist ein praktisches Prüfkriterium und kein exakter Übergang. Plastische Verformung der ungeschwächten Folie erfordert eine elastisch-plastische Platten- oder Membrananalyse. Das örtliche Kerbkriterium kann duktiles Aufreißen, Laserschädigung oder temperaturbedingten Festigkeitsabfall ohne repräsentative Versuche nicht vorhersagen.

Die voreingestellten Festigkeitswerte für 3003-H14 sind beispielhafte Blechwerte bei Raumtemperatur innerhalb des von Smiths angegebenen Bereichs (revision 11 February 2025, p. 2); sie sind keine zertifizierten Folieneigenschaften. Bestimmen Sie die Festigkeit bei der tatsächlichen Folientemperatur unabhängig von der Temperatur des heißen Gases.

Der instationäre Druckverlauf verwendet den Nennberstdruck und nicht das obere Streuungsperzentil. Er setzt den vorgegebenen Gasimpuls und eine konstante Reservoirtemperatur voraus und berücksichtigt weder die Aufheizkinetik noch Flüssigkeitsmitriss, Partikel, zunehmende Verstopfung, Verbrennung, Zellverformung oder dynamische Druckänderungen im Entgasungskanal des Batteriepacks. Der ausgewiesene erforderliche Öffnungsquerschnitt gilt für die maximale Gasbildungsrate bei Sollberstdruck und vollständig geöffneter Entgasungsöffnung; die Drucküberschwingung beim Öffnen wird separat mit der festgelegten Druckgrenze des Gehäuses verglichen. Die Bewertung prüft ausschließlich diese beiden Gaskriterien.

Das implizite BDF2-Verfahren verwendet vier interne Schritte je Ausgabeintervall und eine relative Newton-Toleranz für die Massenstromdichte von 10⁻¹¹. Der angezeigte Höchstdruck ist das Maximum aus 1601 Abtastwerten, die um den Öffnungsbeginn und die frühe Entspannungsphase verdichtet sind. Eine Halbierung der Zeitschritte wurde mit den Voreinstellungen sowie für Fälle mit kleinem Gasvolumen und großem Öffnungsquerschnitt geprüft; die Druckverläufe änderten sich dabei um weniger als 0,003 % des größeren Werts aus Berstdruck und Höchstdruck. Sehr kleine Restüberdrücke gegen Ende der Entspannungsphase können größere relative Fehler aufweisen. Unabhängige Vergleiche mit SciPy Radau sind als Nachweise numerischer Regressionstests beigefügt, nicht als experimentelle Validierung. Das Kerbtiefendiagramm lässt Berechnungswerte außerhalb des angegebenen Gültigkeitsbereichs der elastischen Membrantheorie aus. Das Streuungsintervall erfasst weder systematische Modellabweichungen noch Änderungen des Bruchmechanismus. Die Druck- und Geometrievariationen umfassen jeweils 161 Punkte. Bei der Wiedergabe wird zwischen den berechneten Abtastwerten mit pchip interpoliert; der Klappenwinkel und die Länge der Druckpfeile sind schematisch und keine Vorhersage der Bewegung beim Bruch. Sichern Sie die Einstellungen für die Serienfertigung durch Berstversuche über den Bereich der Prozesstoleranzen und Temperaturen sowie durch repräsentative Gasströmungsversuche unter Fehlbeanspruchung ab.