簡介這份資料是一份ABAQUS隨機(jī)振動分析的實(shí)操講義面向需要完成結(jié)構(gòu)動力學(xué)仿真與產(chǎn)品可靠性評估的工程師、高校師生及有限元初學(xué)者圍繞隨機(jī)激勵(lì)下結(jié)構(gòu)響應(yīng)預(yù)測這一典型工程問題展開。全包僅含1個(gè)PDF文檔體積約2.82MB輕量便于查閱打印內(nèi)容以圖文步驟形式呈現(xiàn)無需額外腳本或工程文件即可跟隨練習(xí)。講義以某電子產(chǎn)品簡化模型為例完整覆蓋IGES實(shí)體導(dǎo)入與shell模型創(chuàng)建、CAE中刪除面與延展表面、修補(bǔ)底座及圓孔、按厚度分割加強(qiáng)區(qū)、材料與截面屬性賦值、結(jié)構(gòu)裝配以及隨機(jī)振動分析設(shè)置等環(huán)節(jié)并給出四底座固支承受隨機(jī)激勵(lì)的載荷邊界條件與結(jié)果解讀思路。文中對導(dǎo)入后密集點(diǎn)影響網(wǎng)格劃分、孔半徑略大于原孔以清除冗余點(diǎn)等細(xì)節(jié)均有提示可幫助讀者規(guī)避常見建模障礙。目前已有686人學(xué)習(xí)下載適合希望系統(tǒng)掌握ABAQUS隨機(jī)振動流程、補(bǔ)齊前處理建模短板的讀者參考。1. 從 IGES 到 shell 幾何ABAQUS隨機(jī)振動分析的前置修模鏈路很多人第一次做隨機(jī)振動卡的不是求解器而是幾何。項(xiàng)目給的是 IGES 實(shí)體交付前卻要轉(zhuǎn)成殼中間得靠 ABAQUS/CAE 一層層刪面、補(bǔ)面、切孔——任何一處幾何有裂縫后面的頻率步就會先崩。這篇以某電子產(chǎn)品在 PSD 基礎(chǔ)激勵(lì)下的響應(yīng)分析為主線把從導(dǎo)入 IGES 到提交random作業(yè)的完整鏈路拆開講。它適合已經(jīng)會點(diǎn) ABAQUS 基本操作、但對隨機(jī)響應(yīng)分析一知半解的工程師也適合拿它當(dāng)中文版操作參照的學(xué)生。所謂隨機(jī)振動本質(zhì)是把結(jié)構(gòu)放到一段功率譜密度描述的基礎(chǔ)激勵(lì)里先求模態(tài)再用模態(tài)疊加算響應(yīng)復(fù)雜度不在公式而在單位、截面、耦合和關(guān)鍵字這四件事上。IGES 文件按實(shí)體的邏輯存面ABAQUS 讀進(jìn)來默認(rèn)會給一批 Solid 體。本例結(jié)構(gòu)是電子產(chǎn)品的殼類件單元類型要按 shell 走所以導(dǎo)入這一關(guān)就得改拓?fù)渲鞑藛蜦ile → Import → Part選文件后在Create Part From IGES File里把Topology切成ShellName 寫random。切完以后你會發(fā)現(xiàn)模型常常比想象的亂——原始 CAD 遺留面、圓孔邊的密集點(diǎn)、肋板和外表面的缺口全在。主菜單Shape → Shell → Remove Face是關(guān)鍵工具點(diǎn)面變紅后中鍵確認(rèn)即可刪除逐塊清到只剩需要的承載表面。刪不是目的補(bǔ)才是。產(chǎn)品外形上常缺 3 個(gè)外表面加 6 個(gè)內(nèi)表面用Shape → Shell → Extrude選中邊緣進(jìn)入 Sketch用Create Lines畫出輪廓線再在Edit Extrusion里把Type設(shè)成Up to Face用Flip調(diào)好箭頭方向最后點(diǎn)底面作為終止面。四角的底座孔更不能馬虎原 IGES 在孔周會留下十幾到幾十個(gè)密集的角點(diǎn)劃分網(wǎng)格時(shí)直接把你逼瘋。做法是用Shape → Cut → Extrude在孔圈上畫一個(gè)半徑略大于原孔的圓——例子里原孔 2.75切到 2.755 就夠把密集點(diǎn)整體切掉。幾何后期還要處理兩類事切薄區(qū)和分割加強(qiáng)區(qū)。厚度不一致的位置要在Tools → PartitionTypeFace、MethodSketch在面上畫線把區(qū)域切出來這樣后面賦不同 Section 才有落腳點(diǎn)。肋板處同樣切一個(gè)矩形加強(qiáng)區(qū)。到這一步能劃分網(wǎng)格和能算準(zhǔn)其實(shí)是兩回事但邊界已經(jīng)畫好了。整個(gè)修模階段最容易翻車的是面與面之間留縫、法向反了、或者 Extrude 方向選錯(cuò)檢查手段是打開Tools → Query → Geometry逐塊點(diǎn)看有沒有 Free edge。2. 單位、截面與耦合約束把電子產(chǎn)品模型裝配成可算的殼單元系統(tǒng)幾何修完接手的是 Property 模塊。這一步看著無聊實(shí)際是隨機(jī)振動最常出事的地方。IGES 給的是毫米量綱材料參數(shù)就必須全程用 mm 制否則頻率算出來會差三個(gè)數(shù)量級隨機(jī)響應(yīng)的量級也跟著一起錯(cuò)。2.1 毫米制材料參數(shù)與四種殼厚截面Material → Create建材料在Mechanical → Elasticity → Elastic填楊氏模量7e7、泊松比0.3再在General → Density填2.7e-6。這里 7e7 是 MPaN/mm22.7e-6 對應(yīng)噸每立方毫米量級的密度配合毫米幾何正好自洽。接著Section → CreateCategory 選Shell依次建出厚度1.6、3、6、8四個(gè)截面。建完別急著Assign先把 Section 名字帶厚度寫清楚比方Section-1_6MM不然后面挑截面時(shí)對著列表發(fā)蒙。部位建議厚度(mm)說明主體外殼1.6最薄最容易在 PSD 下局部共振側(cè)壁3過渡段肋板6承力區(qū)加強(qiáng)底座8接地點(diǎn)剛度大分配用Assign → Section選面群一次性刷完。注意殼單元沒有厚度方向?qū)嶓w厚度只影響彎曲剛度和質(zhì)量如果發(fā)現(xiàn)某塊算出來偏剛先看是不是厚度寫反了部位。2.2 參考點(diǎn)、Coupling 與固支邊界四個(gè)底座要做固支但直接給底座面加密約束容易引起局部應(yīng)力集中慣用做法是打一個(gè)參考點(diǎn)通過 Coupling 把底座面綁到這個(gè)點(diǎn)上。例子里Tools → Reference Point建 RP-1坐標(biāo)默認(rèn)(0,0,0)Constraint → Create選Coupling控制點(diǎn)選 RP-1被控區(qū)域選底座表面上帶網(wǎng)格的那塊面。對話框里記得取消不需要的自由度勾選避免把整個(gè)面自由度全鎖死。約束建完邊界條件用BC → Create直接把Encastre加在 RP-1 上四個(gè)底座通過耦合同時(shí)被固支。# 耦合約束與固支在 INP 里的等價(jià)寫法對照檢查用 *Coupling, Constraint nameConstraint-1, Ref node1, SurfaceSurf-1 *Kinematic 1, 1 1, 2 1, 3 *Boundary 1, ENCASTRE這段代碼檢查的是你 CAE 里耦合和約束到底生成了什么*Coupling后面接*Kinematic是運(yùn)動耦合1,1到1,3表示從節(jié)點(diǎn)跟隨參考點(diǎn)的三個(gè)平動自由度*Boundary里ENCASTRE一次鎖六個(gè)自由度的固定邊界。如果發(fā)現(xiàn)激勵(lì)算完位移量級奇怪先回來看這一對約束是不是把某個(gè)方向多鎖了。3. 頻率步 隨機(jī)響應(yīng)步模態(tài)求解與 PSD 基礎(chǔ)激勵(lì)的參數(shù)配置隨機(jī)響應(yīng)分析在前處理層面永遠(yuǎn)拆成兩步先頻率再隨機(jī)響應(yīng)。第一步求模態(tài)頻率和振型第二步用模態(tài)結(jié)果去疊加 PSD 輸入。少了頻率步求解器會直接報(bào)模態(tài)缺失。3.1 頻率分析步取多少階模態(tài)Step → CreateProcedure type選Linear perturbation子選項(xiàng)選Frequency。進(jìn)入Edit Step在Number of eigenvalues requested的eigen里選Value填20。為什么是 20因?yàn)?PSD 激勵(lì)覆蓋到 2000Hz前 20 階模態(tài)通常已經(jīng)把主要參與質(zhì)量吃掉了——真正確認(rèn)夠不夠得在結(jié)果里看有效模態(tài)質(zhì)量占比。對電子產(chǎn)品這種緊湊結(jié)構(gòu)前幾階往往是整體彎曲高幾階才是局部殼板共振階數(shù)太少會漏掉關(guān)鍵響應(yīng)。3.2 隨機(jī)響應(yīng)步頻率范圍、點(diǎn)數(shù)與偏置第二步Step → CreateProcedure type選Linear perturbation子選項(xiàng)選Random response。在Edit Step的 Basic 頁填參數(shù)取值含義ScaleLinear頻率軸線性Lower Frequency10起始激勵(lì)頻率(Hz)Upper Frequency2000結(jié)束激勵(lì)頻率(Hz)Number of Points20頻響采樣點(diǎn)數(shù)Bias3采樣點(diǎn)向低頻聚集偏置 Bias 是關(guān)鍵。頻率軸取對數(shù)時(shí)10 到 2000 跨度兩個(gè)多數(shù)量級如果點(diǎn)數(shù)均勻鋪開低頻共振區(qū)采樣過疏峰值直接被削平。偏置設(shè) 3 讓采樣點(diǎn)往低頻擠配合 PSD 曲線在 150、350、700、1400Hz 的臺階變化峰值附近才有足夠密度。3.3 阻尼參數(shù)直接模態(tài)阻尼怎么填同一個(gè)對話框切到Damping頁勾using direct damping data在Direct Modal里依次填1、20、0.001。這三個(gè)數(shù)是起始模態(tài)號、結(jié)束模態(tài)號和對應(yīng)的模態(tài)阻尼比含義是第 1 到第 20 階模態(tài)統(tǒng)一用 0.001 的阻尼比。0.001 對電子產(chǎn)品整機(jī)大致落在合理區(qū)間但結(jié)構(gòu)件、裝配剛度不同會差很多做對標(biāo)時(shí)這一步是第一順位要調(diào)的量。填錯(cuò)阻尼最典型的后果是共振峰高度離譜——不是求解器錯(cuò)是你給的能量耗散不對。# 用 Python 直接改 INP 里的模態(tài)阻尼段做掃參 import re with open(random.inp) as f: txt f.read() # 對應(yīng) *MODAL DAMPING, DEFINITIONDIRECT 段 pattern re.compile(r(\*Modal Damping.*?\n)([\d\s\.,\-\n]), re.S) for zeta in [0.001, 0.005, 0.02]: body f1, 20, {zeta}\n txt_new pattern.sub(r\1 body, txt, count1) with open(frandom_zeta_{zeta}.inp, w) as g: g.write(txt_new)這段腳本利用*Modal Damping塊的文本特征對 0.001、0.005、0.02 三檔阻尼比批量生成 INP跑完對比同一節(jié)點(diǎn) U3 的 PSD 響應(yīng)峰。比值一擺出來你就能判斷實(shí)驗(yàn)測得的實(shí)際阻尼大概落在哪一檔比一次性猜參數(shù)靠譜。4. 關(guān)鍵字硬核補(bǔ)丁用 *PSD-DEFINITION 和 *BASE MOTION 打通 CAE 盲區(qū)前處理到這一步會撞上一件讓人頭疼的事ABAQUS/CAE 的圖形界面沒有 PSD 譜定義、基礎(chǔ)運(yùn)動和相關(guān)性的入口。這幾個(gè)關(guān)鍵字必須手工寫進(jìn) INP否則求解器讀不到激勵(lì)只是空轉(zhuǎn)模態(tài)。入口是Model → Edit keywords → Model-1。4.1 PSD 譜定義臺階型基礎(chǔ)加速度譜譜本身是個(gè)臺階函數(shù)本例給到 2000Hz階梯關(guān)鍵點(diǎn)在 10、149、150、349、350、699、700、1399、1400、2000Hz對應(yīng)的G9.81e3是按毫米單位換算的重力加速度。往 INP 的 Model 數(shù)據(jù)段插入*PSD-DEFINITION, NAMEPSD, TYPEBASE, G9.81e3 0.04, 0.0, 10 0.04, 0.0, 149 0.065, 0.0, 150 0.065, 0.0, 349 0.094, 0.0, 350 0.094, 0.0, 699 0.130, 0.0, 700 0.130, 0.0, 1399 0.072, 0.0, 1400 0.072, 0.0, 2000每行三個(gè)數(shù)PSD 值、實(shí)部與虛部虛部為 0代表純實(shí)功率譜、頻率。相鄰行頻率成對給出正是為了讓譜在 149→150 之間做垂直跳變而不是斜線過渡。TYPEBASE表示這是基礎(chǔ)加速度激勵(lì)而不是力或位移譜。4.2 基礎(chǔ)運(yùn)動與相關(guān)性光有譜不夠還要指定在哪個(gè)自由度上施加、用哪個(gè) load case*BASE MOTION, DOF1, LOAD CASE1 *CORRELATION, PSDPSD, TYPEUNCORRELATED 1, 1.DOF1表示基礎(chǔ)激勵(lì)作用在 X 向平動。*CORRELATION里TYPEUNCORRELATED與后面1, 1.配合表示多方向激勵(lì)之間互不相關(guān)、自相關(guān)系數(shù)為 1。如果工況需要三向同時(shí)激勵(lì)這里要擴(kuò)展成互不相關(guān)的多組并對應(yīng)新增 load case。寫完關(guān)鍵字必須回輸入文件確認(rèn)位置——它要落在 Modal 數(shù)據(jù)段之后、History 數(shù)據(jù)段之前插錯(cuò)位置求解器會以unknown keyword報(bào)錯(cuò)別看錯(cuò)在別處去。提示每次改完 Edit keywords 都回 CAE 里再點(diǎn)一次 Write Input讓圖形界面把后續(xù)改動落到 INP 里別只在關(guān)鍵字編輯窗口改而不落盤。5. 網(wǎng)格、輸出與后處理結(jié)果能不能用全看這三步幾何、材料、關(guān)鍵字都齊了隨機(jī)響應(yīng)能不能給出可用結(jié)果最后壓在三件事上網(wǎng)格質(zhì)量、輸出請求、后處理取值。任意一處偷懶報(bào)出來的 PSD 都可以很美但沒意義。5.1 全局種子與 Quad-dominated 控制Mesh 模塊里Seed → Instance選全部三個(gè)部件Global Seeds的Approximate global size填3.5。這個(gè)數(shù)針對毫米量綱大致意味著單元邊長 3.5mm對厚度 1.6mm 的殼來說偏粗但能壓低總自由度。要更細(xì)可以降到 2代價(jià)是頻率步和隨機(jī)響應(yīng)步都明顯變慢。肋板區(qū)另用Mesh → Controls單獨(dú)處理Element Shape選Quad-dominated把四塊肋板選上避免四邊形被劈成四面體導(dǎo)致殼單元退化。劃分完成后先看兩點(diǎn)一是有沒有報(bào)未連接節(jié)點(diǎn)直接搜 INP 里的*Nset數(shù)量或回 CAE 用Tools → Query → Mesh找孤立點(diǎn)二是孔邊圓整度如果圓孔變成多邊形返回 4.3 的切除步驟把切圓半徑再調(diào)大一點(diǎn)。5.2 歷史輸出與場輸出一起抓隨機(jī)振動最有信息量的輸出是歷史輸出本例關(guān)心薄弱環(huán)節(jié)的位移。Tools → Set → Create建 Set-1 圈住關(guān)鍵區(qū)域再Output → History output requests → CreateDomain選Set輸出U1、U2、U3。這里選的是 Step-2才不會把頻率步的無效輸出混進(jìn)來。同時(shí)Output → Field output requests加一份默認(rèn)場輸出方便在云圖上找熱點(diǎn)。*Output, history, setSet-1, frequency1 U1, U2, U3frequency1表示每一步增量都寫一次歷史輸出隨機(jī)響應(yīng)里就是每個(gè)采樣頻率都記。Set-1 若不小心圈到耦合約束面上U1、U2、U3 會被參考點(diǎn)位移污染取結(jié)果前先看節(jié)點(diǎn)編號確認(rèn)圈選范圍。5.3 后處理先看頻率再讀 RMS 與峰值提交random作業(yè)Job Manager 里點(diǎn)Monitor盯住Standard日志重點(diǎn)看兩行模態(tài)求解是否收斂、隨機(jī)響應(yīng)有沒有出現(xiàn)zero mass類警告。算完進(jìn) VisualizationResult → Field Output里點(diǎn)Step/Frame先確認(rèn)低階頻率是否落在合理區(qū)間——電子產(chǎn)品整機(jī)一階通常幾十到幾百赫茲如果一階跑到幾千赫茲八成是單位錯(cuò)了或者截面厚度貼反。接著切到History OutputX 軸選頻率、Y 軸選 U3 的 PSD 曲線找峰值對應(yīng)的頻率與幅值逐級對比 PSD 輸入曲線的臺階。真正拿去做設(shè)計(jì)判斷的是 RMS。對PSD U3曲線在 10 到 2000Hz 積分再開方得到該節(jié)點(diǎn)的位移均方根用它去對標(biāo)產(chǎn)品規(guī)范里的加速度或位移限值。峰值頻率和 RMS 一起看峰值頻率落在激勵(lì)臺階內(nèi)、RMS 明顯偏大的區(qū)域就是下一輪改設(shè)計(jì)要加厚、加筋或者加裝減振的位置。注意修改幾何或截面后關(guān)鍵字塊要重新落到 INP。CAE 會保留你手寫的關(guān)鍵字但一旦切了 Step 或重寫輸入位置就容易跑偏每次重提作業(yè)前掃一眼*PSD-DEFINITION是否還在正確的數(shù)據(jù)段。本文還有配套的精品資源點(diǎn)擊獲取