本文探討了輸出濾波在 Abaqus/Explicit 和 Abaqus/Viewer 中的關鍵作用，強調了其在減少雜訊和確保準確解讀模擬資料方面的重要性。它解釋了混疊（aliasing）的挑戰，這是一種因取樣率不足而扭曲資料的現象，並透過執行階段（runtime）與後處理（post-processing）濾波器提供了實用的解決方案。透過掌握這些技術，工程師可以獲得乾淨、可靠的資料，準確反映其模型的物理行為，進而做出明智的設計決策。
主要重點
什麼是輸出濾波？ 輸出濾波可減少模擬資料中不必要的雜訊，突顯相關的物理響應，以提升清晰度與可靠性。
混疊的問題： 當資料取樣率過低，無法擷取高頻訊號時，就會發生混疊，導致誤導性的結果。奈奎斯特頻率（Nyquist frequency）是避免此問題的關鍵基準。
Abaqus 中的濾波解決方案：
- 執行階段濾波器 (Abaqus/Explicit)： 透過在儲存資料前進行濾波來防止混疊。內建的抗混疊濾波器非常有效。
- 後處理濾波器 (Abaqus/Viewer)： 允許在模擬後對資料進行進一步的細化與平滑化。
輸出濾波的最佳實踐：
- 執行頻率分析以識別重要的動態特徵。
- 對輸出取樣率使用 6 倍經驗法則，以確保資料保真度。
- 避免過度濾波，這可能會消除有意義的結構響應。
- 使用執行階段濾波器操作（如 MAX 或 ABSMAX）有效地監控極值。
減輕濾波器假影（Filter Artifacts）： 意識到由濾波引起的相位偏移（time shifts）和端點失真（end distortions），並透過設計模擬來最小化其影響。
穩健的工作流程： 結合頻率分析、執行階段濾波和後處理，以確保為關鍵工程決策提供乾淨、可解讀的資料。
什麼是輸出濾波？
輸出濾波是對模擬結果應用數學技術的過程，以減少不必要的雜訊並突顯系統相關的物理響應。透過過濾模擬輸出，工程師可以提高資料的清晰度與可靠性，使其更容易解讀真實的系統行為並得出準確的結論。
簡介
在有限元素分析中，特別是動態顯式模擬，產生的資料可能龐大且複雜。工程師面臨的一個常見挑戰是區分真實的結構響應與高頻求解雜訊。來自 Abaqus/Explicit 的原始資料可能充滿了振盪，這些振盪雖然在數學上是求解的一部分，卻可能掩蓋模型潛在的物理行為。這就是輸出濾波成為不可或缺工具的地方。正確應用濾波器使分析師能夠衰減這種雜訊，防止嚴重的資料誤讀，並獲得對模型動態效能清晰、可化為行動的理解。
本文提供了對 Abaqus/Explicit 和 Abaqus/Viewer 中輸出濾波的技術原理與實際應用的深入探討。我們將探討混疊的基本問題，詳細說明可用的廣泛濾波選項，並為配置模擬以產生乾淨、可靠且具洞察力的資料提供穩健的指導。
混疊的危機：當資料欺騙人時
就其核心而言，混疊是一種由輸出頻率不足引起的資料損壞形式。想像一下，如果您每秒只看一次，試圖去了解一個快速旋轉的輪子。取決於您看的時間點，輪子可能看起來是靜止的、緩慢旋轉的，甚至是向後旋轉的。您正在以過慢的速率對現實進行取樣，而無法捕捉到真實的運動。在模擬中，同樣的現象也會發生。
奈奎斯特-夏農取樣定理（Nyquist-Shannon sampling theorem）為這個概念提供了數學基礎。
這個定理指出，為了準確地重建訊號，取樣頻率必須至少是該訊號中存在的最高頻率分量的兩倍。這個閾值被稱為奈奎斯特頻率，是資料擷取中的一個關鍵基準。
當您以低於此速率的要求輸出模擬時，您就冒著混疊的風險。您儲存的資料點個別來看是正確的，並且位於真實的高頻曲線上，但當把它們繪製出來時，它們會創造出一個現實中不存在的、具誤導性的新波形。
視覺化混疊：從幻影箱到渦輪機
一個經典、直觀的混疊範例是幻影箱（zoetrope）。這個 19 世紀的裝置從一系列靜態影像中創造出運動的錯覺。在現代模擬環境中，想像一個以特定影格率動畫化的快速旋轉剛體。由於旋轉速率和輸出速率之間的關係，物體可能會看起來發生變形和彎曲，就像是由柔軟材料製成的一樣。這種「影格率效應」是混疊的直接視覺表現，導致對物體剛體運動的完全誤解。
一個更直接的工程範例涉及模擬渦輪風扇，其中一片葉片脫落並撞擊外殼。當動畫影格儲存得太不頻繁時，風扇可能看起來像在向後旋轉。不熟悉混疊的分析師可能會浪費寶貴的時間去調查一個不存在的邊界條件問題，認為規定的旋轉沒有被正確執行。問題不在於物理模擬；而在於觀察它的方式。為了看到真實的逆時針旋轉，必須儲存更多的影格，以更小的旋轉增量來捕捉風扇的位置。
XY 圖表中的混疊：懸臂樑案例研究
混疊對時間歷程圖的影響同樣深遠。考慮一個簡單的懸臂樑，承受突然釋放的壓力負載，導致其自由振動。
- 「真實」訊號： 如果我們在 Abaqus/Explicit 模擬的每一個時間增量（time increment）都記錄樑尖端的垂直速度，我們將獲得一條乾淨、高解析度的曲線，準確地描繪了振動運動。這代表了真實數據。
- 混疊訊號： 現在，如果我們要求相同的輸出變數，但頻率低得多（例如，每隔幾百個增量），結果圖表將截然不同。我們收集的少數資料點落在真實曲線上，但它們錯過了中間所有的波峰和波谷。連接這些點創造了一個新的、較低頻率的訊號，完全誤傳了樑的行為。
在這樣一個簡單的案例中，我們有幸擁有高解析度資料進行比較。在複雜、真實世界的分析中，我們通常只有低頻資料。如果沒有對混疊的了解，工程師可能會錯誤地得出結論，認為結構的振動響應比實際上慢得多，且振幅也較低，從而導致有缺陷的設計決策。
哪些結果最容易受到混疊影響？
在混疊方面，並非所有輸出變數都是平等的。易受影響的程度與變數的「頻率內容」直接相關。涉及積分的操作傾向於平滑訊號，消除高頻分量。相反地，涉及微分的操作傾向於放大它們。
這在結構分析中導致了清晰的易受影響程度光譜：
- 低易受影響度（位移）： 位移是速度隨時間積分的結果，而速度本身又是加速度的積分。這種雙重積分過程自然會過濾掉高頻雜訊。因此，位移圖通常很平滑，較不容易出現混疊問題。
- 中易受影響度（速度）： 作為加速度的第一次積分，速度包含比位移更多的頻率內容，但通常比加速度的雜訊少。它佔據了中間地帶。
- 高易受影響度（加速度和反作用力）： 這些量通常是顯式動態分析中「最原始」的訊號。它們直接反映了通過網格元素傳播的高頻波，因此極其嘈雜且極易受到混疊影響。
準靜態模擬儘管在 Abaqus/Explicit 中執行，但通常具有非常緩慢的負載率，這將動態效應降至最低，並使混疊成為不成問題的問題。然而，對於任何涉及衝擊、碰撞或高頻振動的分析，仔細考慮對加速度和力的輸出要求是至關重要的。
雙管齊下的解決方案：Abaqus/Explicit 與 Abaqus/Viewer 中的濾波
Abaqus 提供了全面的濾波工具套件來對抗雜訊與混疊。這些工具可以在兩個不同的階段部署：在模擬執行期間（執行階段濾波）以及在模擬完成後（後處理）。
1. Abaqus/Explicit 中的執行階段濾波
執行階段濾波器是第一道也是最關鍵的防線。它們在資料被寫入輸出資料庫（ODB）檔案之前對資料進行處理。這種方法從根本上講很重要，因為一旦訊號被降頻取樣且混疊的資料被儲存到 ODB 中，遺失的高頻資訊就無法復原了。
Abaqus/Explicit 提供了幾種工業標準的濾波器類型，您可以使用 *FILTER 關鍵字來定義：
- 巴特沃斯（Butterworth）： 以其在通帶（passband）中極致平坦的頻率響應而聞名，使其成為非常受歡迎的通用濾波器。
- 切比雪夫一型（Chebyshev Type I）： 展現出比巴特沃斯濾波器更陡峭的衰減（roll-off），但代價是在通帶中引入漣波（ripples）。
- 切比雪夫二型（Chebyshev Type II）： 在阻帶（stopband）中而非通帶中有漣波。
內建的抗混疊濾波器：您的預設選擇
對於大多數應用來說，並不需要我們定義自訂濾波器。Abaqus 提供了一個穩健的、預先配置的抗混疊濾波器，這非常有效。當您要求此濾波器時，Abaqus 會自動實施一個二階、低通的巴特沃斯濾波器。其智慧在於它如何設定截止頻率。
截止頻率會被自動定義為您在輸出要求中指定的輸出取樣率的六分之一（1/6）。
讓我們用一個例子來分解它。假設您要求每 0.1667 毫秒輸出一次歷史資料。
- 輸出取樣率： 1 / (0.1667 x 10⁻³ s) = 6,000 Hz (6 kHz)。
- 奈奎斯特頻率： 根據取樣定理，您可能解析出的最高頻率是取樣率的一半，即 3 kHz。原始訊號中任何高於 3 kHz 的頻率內容都有引起混疊的風險。
- 濾波器截止頻率： 抗混疊濾波器將其截止頻率設定為取樣率的 1/6：6,000 Hz / 6 = 1,000 Hz (1 kHz)。
透過將截止頻率設定為 1 kHz，濾波器強烈衰減了所有高於此頻率的訊號，確保在訊號被取樣並寫入 ODB 之前，任何接近有問題的 3 kHz 奈奎斯特頻率的內容都會被消除。這提供了訊號乾淨、無混疊的表示。使用此濾波器的一個顯著好處是，如果您選擇的輸出率過低而導致濾波器無法發揮作用，Abaqus 將發出警告，作為防止混疊的內建安全檢查。
2. Abaqus/Viewer 中的後處理濾波器
如果您已成功將豐富的高頻訊號儲存到 ODB 中（理想情況下是使用執行階段濾波器），您可以在 Abaqus/Viewer 內執行額外的濾波。這透過 Operate on XY data 對話方塊完成。後處理對於以下情況很有用：
- 在不重新執行分析的情況下比較不同的濾波器類型和截止頻率。
- 為了簡報目的應用平滑化。
- 進一步減少在執行階段未完全消除的雜訊。
一個關鍵區別：單次傳遞與雙向濾波
執行階段和後處理濾波器之間存在一個非常重要的技術區別。
- Abaqus/Explicit 執行階段濾波器： 這些是單次傳遞、因果濾波器（causal filters）。隨著模擬的進行，它們只在時間的前進方向上處理資料。這是一種物理上的必然性——求解器無法「預見未來」以知道即將到來的資料。這種單次傳遞操作的一個結果是，在過濾後的輸出中會產生相位偏移（phase shift），或時間延遲。與原始訊號相比，過濾後曲線的波峰將會稍微向右移動（在時間上較晚）。
- Abaqus/Viewer 濾波器： 這些是雙向（或雙重傳遞）濾波器。它們在時間上向前和向後處理整個資料序列。向後傳遞有效地校正了由向前傳遞引入的相位偏移。結果，Viewer 過濾後的訊號在時間上與原始未過濾訊號完美對齊。
雖然執行階段濾波器引起的這種時間偏移通常很小，但它是一種真實的失真，分析師應該意識到這一點，特別是在試圖精確關聯時間上的事件時。
最佳實踐
理論是有價值的，但有效的濾波需要務實的方法。目標是在資料保真度與可管理的檔案大小之間取得平衡。
· 如何選擇您的輸出頻率
理想的輸出頻率應夠高以捕捉結構上顯著的響應，但不應高到產生難以管理的龐大 ODB 檔案。系統化的方法是最好的：
- 執行自然頻率分析： 在執行您的顯式動態模擬之前，在同一個模型上執行 Abaqus/Standard 的自然頻率提取分析（*FREQUENCY）。這將告訴您結構主要振動模態的特徵頻率。
- 識別最高結構顯著頻率： 您的結構可能有數百個模態，但通常只有前幾個在結構上是顯著的。識別您真正關心的工程問題中最高頻率的面相。對於許多結構振動問題，此範圍介於 50 Hz 和 5,000 Hz 之間。這通常遠低於整個模型中的最高頻率，該頻率控制著穩定時間增量（通常在數百千赫茲）。
- 設計您的輸出率： 使用最高的顯著頻率來決定您的輸出率。一個好的經驗法則（rule of thumb）是將濾波器截止頻率設定在等於或略高於此數值。如果您使用內建的抗混疊濾波器，您可以反向推算：
- 期望的截止頻率： 假設為 5,000 Hz。
- 所需的輸出取樣率： 這必須是截止頻率的 6 倍，所以 6 * 5,000 Hz = 30,000 Hz。
- 所需的輸出時間間隔： 1 / 30,000 Hz ≈ 3.33 x 10⁻⁵ 秒 (0.0333 ms)。
使用抗混疊濾波器以此間隔要求歷史輸出，提供了一個穩健的起點。如果磁碟空間允許，增加這個取樣率（例如，增加 2 倍）通常是務實的。資料太多總比不夠好。
過度濾波的危險
正如取樣不足是有問題的一樣，過度濾波也是。如果您將濾波器截止頻率設定得太低，您可能會無意中消除結構響應中具有物理意義的部分。
重新審視嘈雜加速度訊號的例子，如果我們應用一個截止頻率非常低的濾波器，結果曲線可能看起來平滑且乾淨，但它會過濾掉整個結構振動特徵。過度過濾的訊號可能會使峰值振幅急劇降低並產生顯著的時間偏移，使其看起來更像是一個重度阻尼、單自由度系統的響應，而不是複雜的結構振動。這是另一種形式的誤解。
· 進階技術：監控極值
執行階段濾波器一個強大且常被低估的功能是監控極值的能力。透過在您的 *FILTER 或 *OUTPUT 定義中添加 OPERATION 參數，您可以追蹤峰值而無需儲存整個時間歷程。
- OPERATION = MAX： 報告直到該時間點遇到的最大值。
- OPERATION = MIN： 報告最小值。
- OPERATION = ABSMAX： 報告最大絕對值。
這對於場輸出（field output）異常有用。想像一下在碰撞事件中試圖尋找大型模型中的最大主應變。您不需要為數千個時間增量儲存等高線圖，而是可以在單一時間點（步驟結束時）以 OPERATION = MAX 要求應變變數的場輸出。產生的等高線圖將顯示每個元素在整個模擬過程中經歷的峰值。Abaqus 還會將發生此峰值的時間寫入結果檔案，讓您能有效地識別最關鍵的時刻，以便進行進一步的詳細分析。
· 減輕固有的濾波器失真
濾波器很強大但並不完美。除了前面討論的時間偏移外，另一個需要注意的假影（artifact）是端點失真（end distortion）。
濾波演算法需要一定視窗大小的資料點來進行計算。在訊號的最初和最後，這個視窗是不完整的。這可能導致過濾後的輸出在開始和結束時間附近出現失真或「振鈴效應（ringing）」。如果您有一個關鍵感興趣的事件正好發生在模擬步驟的最後，該事件過濾後的表示可能會因為端點失真而受到影響。一個簡單的解決方案是延長模擬時間，使感興趣的事件發生在分析期間的充裕範圍內，遠離端點。
結論
輸出濾波不是可有可無的附加功能；它是顯式動態分析中良好工程實踐的基礎組成部分。透過精通這些工具，您可以從嘈雜、模稜兩可的資料轉向清晰、可解讀的結果。
一個穩健的工作流程是關鍵：
- 了解您的結構： 使用頻率分析來決定模型中顯著的動態特徵。
- 從源頭濾波： 永遠在 Abaqus Explicit 中使用執行階段濾波器，以建立一個豐富、無混疊的 ODB 檔案。內建的抗混疊濾波器因其簡單性與有效性而成為推薦的起點。
- 明智地選擇輸出率： 根據您期望的結構響應設定輸出頻率，為抗混疊濾波器使用 6 倍經驗法則。有疑問時，儲存更多資料。
- 善用後處理： 在 Abaqus/Viewer 中使用雙向濾波器進行額外的平滑化和比較分析，確信您的底層資料是可靠的。
- 注意假影： 了解時間偏移和端點失真的影響，並設計您的模擬以最小化它們對您解讀的影響。
透過採用這種有條理的方法，工程師可以充滿信心地駕馭動態資料的複雜性，確保其模擬結果為關鍵的設計決策提供堅實、準確的基礎。若要了解更多關於輸出濾波的資訊，請在此處查看我們的網路研討會。
在 3DEXPERIENCE 平台上延伸輸出濾波工作流程
在 3DEXPERIENCE 平台內，Abaqus 中的輸出濾波工作流程可作為基於模型、以資料為中心的模擬環境的一部分進行管理。平台並非將濾波視為孤立的求解器或 Viewer 操作，而是允許將其包含在可追溯的模擬生命週期中。
透過 SIMULIA 角色（例如物理模擬或結構模擬工程師）執行的 Abaqus 分析，會將濾波定義（包括 *FILTER 參數、輸出要求和取樣率）保留在模擬物件或關聯的輸入資料中。這確保了：
- 濾波配置是持久且受版本控制的： 對截止頻率、抗混疊設定和輸出間隔的更改，會與模型修訂一起被追蹤。
- 維護模擬來源： 結果可以與用於產生它們的輸入配置關聯起來。
- 集中管理模擬資料： 以高取樣率產生的大型資料集在平台環境內被儲存與治理，提高了可及性並減少了重複。
通常在 Abaqus/Viewer 中執行的後處理操作（包括 XY 資料濾波與訊號處理），在需要時可透過指令碼或流程自動化工作流程來重現。這支援了：
- 跨團隊濾波方法論的標準化。
- 當透過使用者定義的流程實施時，結果提取的自動化，包括過濾後的訊號和極值監控（例如 MAX/ABSMAX 操作）。
- 重複使用模擬流程，提高跨專案的一致性。
透過將輸出濾波納入受治理的資料環境中，3DEXPERIENCE 平台支援對混疊和訊號保真度更一致的處理，同時保持模擬結果如何產生與解讀的可追溯性。
原文轉載來自：Mastering Output Filtering in Abaqus/Explicit & Abaqus/Viewer