1. 項(xiàng)目概述CRTSⅡ型軌道精細(xì)模型的定位這幾年做高鐵無(wú)砟軌道結(jié)構(gòu)分析繞不開(kāi)的一個(gè)問(wèn)題就是怎么把列車移動(dòng)荷載真實(shí)地加到有限元模型里而不是用靜態(tài)荷載糊弄過(guò)去。用Abaqus建立CRTSⅡ型軌道精細(xì)模型配合DLOAD子程序?qū)崿F(xiàn)列車移動(dòng)荷載的施加是目前工程和科研里比較主流的做法。這套方案既能拿到軌道板、底座板、CA砂漿層內(nèi)部的應(yīng)力應(yīng)變細(xì)節(jié)又能反映移動(dòng)荷載下結(jié)構(gòu)動(dòng)力響應(yīng)的真實(shí)規(guī)律比單純?cè)陟o力分析里加一個(gè)固定力要靠譜得多。這個(gè)項(xiàng)目適合誰(shuí)參考一類是做無(wú)砟軌道結(jié)構(gòu)設(shè)計(jì)驗(yàn)算和疲勞評(píng)估的工程師另一類是研究車輛-軌道耦合動(dòng)力學(xué)、軌道結(jié)構(gòu)損傷演化的研究生和科研人員還有剛接觸Abaqus二次開(kāi)發(fā)、想搞清楚DLOAD子程序怎么用于移動(dòng)荷載的同行。不管你屬于哪一類這套從建模到子程序再到后處理的完整流程都是可以直接落地復(fù)用的。我下面寫的內(nèi)容全部來(lái)自一個(gè)實(shí)際項(xiàng)目的完整復(fù)現(xiàn)過(guò)程涉及幾何尺寸、材料參數(shù)、網(wǎng)格策略、子程序代碼和各類報(bào)錯(cuò)排坑一步步來(lái)。2. 建模方案與結(jié)構(gòu)簡(jiǎn)化2.1 CRTSⅡ型軌道的結(jié)構(gòu)組成與幾何參數(shù)CRTSⅡ型板式無(wú)砟軌道和CRTSⅠ型、CRTSⅢ型最大的區(qū)別在于軌道板是連續(xù)結(jié)構(gòu)沿縱向有預(yù)應(yīng)力筋貫穿中間沒(méi)有斷開(kāi)的板縫。整體從下往上依次是底座板、CA砂漿調(diào)整層、軌道板、扣件系統(tǒng)和鋼軌。這種結(jié)構(gòu)形式在高速鐵路線上很常見(jiàn)尤其是長(zhǎng)軌鋪設(shè)的區(qū)段整體性好、剛度均勻?qū)Ω咚龠\(yùn)行的列車來(lái)說(shuō)平順性優(yōu)勢(shì)很明顯。建模前先把尺寸理清楚。鋼軌采用CHN60軌高176mm底部寬150mm扣件采用WJ-8B型小阻力扣件間距一般在630mm左右軌道板厚度一般取200mm寬度2550mmCA砂漿層厚度約30mm底座板厚度約300mm比軌道板寬一些具體寬度和配筋方式根據(jù)線下基礎(chǔ)類型略有差異。需要說(shuō)明的是不同線路的局部尺寸會(huì)有調(diào)整建模前先拿到設(shè)計(jì)圖確定具體數(shù)值別直接抄本文的數(shù)值否則算出來(lái)跟實(shí)際結(jié)構(gòu)對(duì)不上。如果做整段區(qū)間計(jì)算量會(huì)非常大。一般做法是取軌道結(jié)構(gòu)的一段代表性長(zhǎng)度來(lái)建模比如含4到5塊軌道板范圍的連續(xù)段總長(zhǎng)10到20m然后施加電周期性邊界或者用黏彈性邊界模擬無(wú)限長(zhǎng)的軌道基礎(chǔ)。我這里取了一個(gè)18m長(zhǎng)的軌道段模型包含了完整的底座板、砂漿層、兩塊連續(xù)軌道板以及對(duì)應(yīng)范圍的扣件和鋼軌。2.2 材料參數(shù)與單位體系的選擇Abaqus本身沒(méi)有固定的單位制需要自己保證量綱一致。軌道結(jié)構(gòu)分析常用mm-N-s-tonne這套單位體系長(zhǎng)度用mm力用N時(shí)間用s質(zhì)量用tonne對(duì)應(yīng)的應(yīng)力單位就是MPaN/mm2密度單位是tonne/mm3。這套單位的優(yōu)勢(shì)在于和工程圖紙尺寸直接兼容不需要做數(shù)值換算。材料參數(shù)按常見(jiàn)工程取值來(lái)給。軌道板是C50混凝土彈性模量3.55×10? MPa泊松比0.2密度2.5×10?? tonne/mm3。底座板C30或C40混凝土彈性模量取3.0×10?到3.25×10? MPa。CA砂漿層是彈模相對(duì)低很多的材料一般在100到300 MPa之間泊松比0.2左右這層是整個(gè)軌道結(jié)構(gòu)里比較薄弱的環(huán)節(jié)也是車致?lián)p傷分析的重點(diǎn)關(guān)注對(duì)象。鋼軌用鋼E2.06×10? MPa泊松比0.3密度7.85×10?? tonne/mm3?？奂到y(tǒng)用彈簧阻尼單元模擬豎向剛度一般在20到80 kN/mm之間阻尼20到40 kN·s/m。提示如果你的模型單位體系不同注意把速度單位也換算過(guò)來(lái)。比如mm-s單位下300 km/h要寫成83333.3 mm/s寫錯(cuò)單位是最常見(jiàn)的低級(jí)錯(cuò)誤。2.3 網(wǎng)格劃分與單元類型選擇網(wǎng)格策略是精細(xì)模型的關(guān)鍵取舍點(diǎn)直接決定計(jì)算精度和成本。鋼軌如果只關(guān)心荷載傳遞效果可以用梁?jiǎn)卧肟吹戒撥壗孛鎯?nèi)部接觸應(yīng)力就得用實(shí)體單元。這里我建議用實(shí)體單元建模鋼軌軌頭接觸區(qū)網(wǎng)格加密到2到3mm軌底過(guò)渡區(qū)5到8mm這樣移動(dòng)荷載下鋼軌的彎曲應(yīng)力和接觸區(qū)應(yīng)力都能算得比較準(zhǔn)。軌道板和底座板都是大體積混凝土構(gòu)件用C3D8R六面體縮減積分單元常規(guī)區(qū)域網(wǎng)格尺寸30到50mm扣件區(qū)域局部加密到10mm左右。CA砂漿層厚度只有30mm沿厚度方向至少要?jiǎng)澐?到3層單元否則彎曲響應(yīng)出不來(lái)。網(wǎng)格數(shù)量控制在50萬(wàn)到100萬(wàn)之間配合隱式求解和子程序計(jì)算單次動(dòng)力響應(yīng)分析可在可接受的時(shí)間內(nèi)完成。單元類型選C3D8R時(shí)要注意沙漏控制尤其是在沖擊荷載作用下縮減積分單元容易產(chǎn)生沙漏變形增加單元網(wǎng)格密度或使用增強(qiáng)沙漏控制選項(xiàng)可以緩解。軌道板受彎為主的受力模式下也可以考慮C3D20R二次單元提高彎曲精度但計(jì)算成本會(huì)明顯上升權(quán)衡下來(lái)C3D8R加合理網(wǎng)格密度是更劃算的方案。3. DLOAD子程序的編寫與核心實(shí)現(xiàn)3.1 DLOAD的工作機(jī)制與調(diào)用原理DLOAD是Abaqus標(biāo)準(zhǔn)求解器中用于施加與位置、時(shí)間相關(guān)分布荷載的接口。每個(gè)帶分布荷載的積分點(diǎn)都會(huì)在增量步開(kāi)始時(shí)調(diào)用一次子程序Abaqus把當(dāng)前計(jì)算時(shí)間的積分點(diǎn)坐標(biāo)傳遞給子程序子程序根據(jù)坐標(biāo)和時(shí)間判斷這個(gè)點(diǎn)是否處于荷載作用范圍內(nèi)并返回對(duì)應(yīng)的荷載值F。F的單位是壓強(qiáng)即力除以面積所以最終施加在鋼軌頂面上的荷載是通過(guò)一個(gè)壓力帶的形式來(lái)模擬移動(dòng)輪載的。搞清楚這個(gè)機(jī)制后實(shí)現(xiàn)移動(dòng)荷載的思路就很清晰了每個(gè)時(shí)刻計(jì)算一組輪對(duì)的空間位置然后判斷鋼軌頂面各積分點(diǎn)是否落在輪載作用區(qū)段內(nèi)。落在區(qū)段內(nèi)的點(diǎn)施加相應(yīng)的壓力值區(qū)段外的點(diǎn)返回零。如果采用隱式求解這個(gè)荷載會(huì)跟隨時(shí)間步長(zhǎng)更新從空間上看就是一組沿鋼軌縱向移動(dòng)的壓力帶。一個(gè)需要提前說(shuō)明的地方是DLOAD作用范圍是按當(dāng)前荷載作用區(qū)域來(lái)識(shí)別的每個(gè)積分點(diǎn)只能屬于有荷載或者無(wú)荷載兩種狀態(tài)。如果直接寫一個(gè)if判斷荷載區(qū)段首尾會(huì)出現(xiàn)階躍變化也就是壓力從零瞬間跳到滿值這種突變?cè)陔[式求解里容易造成收斂困難。所以實(shí)際編寫時(shí)要給荷載區(qū)段的邊界設(shè)置過(guò)渡區(qū)域用線性過(guò)渡或平滑過(guò)渡的方式讓壓力值漸變上升和下降計(jì)算穩(wěn)定性會(huì)好很多。3.2 單個(gè)輪對(duì)移動(dòng)荷載的數(shù)學(xué)表達(dá)先解決單輪對(duì)的實(shí)現(xiàn)問(wèn)題。假設(shè)列車沿軌道縱向X方向行駛初始時(shí)刻輪對(duì)位于X0處運(yùn)行速度為V那么在時(shí)間T時(shí)輪對(duì)的位置可以寫成Xwheel X0 V × T需要明確的是DLOAD子程序里TIME(1)指的是分析步的累計(jì)時(shí)間。多分析步情況下如果荷載只在第二個(gè)分析步開(kāi)始施加要注意TIME(1)的基準(zhǔn)點(diǎn)是從當(dāng)前分析步還是整個(gè)分析開(kāi)始計(jì)算可以先用寫入外部文件的方式確認(rèn)一下時(shí)間基準(zhǔn)避免出現(xiàn)荷載位置錯(cuò)位的問(wèn)題。在Abaqus代碼里判斷當(dāng)前積分點(diǎn)是否在荷載作用區(qū)內(nèi)的邏輯如下如果ABS(COORDS(1) - Xwheel)小于等于荷載分布半長(zhǎng)L/2那么這個(gè)積分點(diǎn)位于荷載帶內(nèi)F取為輪載壓強(qiáng)否則F取0。輪載壓強(qiáng)怎么定假設(shè)軸重為14噸單個(gè)輪載為70kN即7×10? N。荷載分布區(qū)長(zhǎng)200mm鋼軌頂面荷載作用寬度按50mm估算壓力帶的承載面積就是200×5010000mm2對(duì)應(yīng)壓強(qiáng)F70000/100007 MPa。這個(gè)壓強(qiáng)值在上述單位制下恰好為7 N/mm2。完整的單輪對(duì)DLOAD子程序如下SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NODE,NOEL,NPT,LAYER, 1 KSPT,COORDS,JLTYP,SNAME) C INCLUDE ABA_PARAM.INC C DIMENSION COORDS(3), TIME(2) CHARACTER*80 SNAME C REAL*8 V, X0, XWHEEL, LZONE, PRESS, DLOADWIDTH C C 參數(shù)定義 V 83333.3D0 ! 列車速度 mm/s對(duì)應(yīng)300km/h X0 1000.0D0 ! 初始輪對(duì)位置 mm LZONE 200.0D0 ! 荷載分布區(qū)長(zhǎng)度 mm DLOADWIDTH 50.0D0 ! 荷載分布區(qū)寬度 mm PRESS 70000.0D0 / (LZONE * DLOADWIDTH) ! 輪載壓強(qiáng) N/mm2 C C 當(dāng)前時(shí)刻輪對(duì)位置 XWHEEL X0 V * TIME(1) C C 判斷積分點(diǎn)是否在荷載作用區(qū)內(nèi) IF (DABS(COORDS(1) - XWHEEL) .LE. LZONE / 2.0D0) THEN F PRESS ELSE F 0.0D0 END IF C RETURN END這里把荷載定義成一個(gè)200mm長(zhǎng)的均布帶實(shí)際車輪鋼軌接觸斑沿縱向也就十幾毫米200mm是一個(gè)等效分布長(zhǎng)度目的是在網(wǎng)格尺寸不小的情況下也能保證荷載帶覆蓋至少一個(gè)單元的積分點(diǎn)避免壓力只落在個(gè)別單元上導(dǎo)致局部應(yīng)力失真。網(wǎng)格越密這個(gè)分布長(zhǎng)度可以越接近真實(shí)接觸斑尺寸。3.3 多個(gè)輪對(duì)與整列車荷載的擴(kuò)展實(shí)現(xiàn)高速列車一個(gè)轉(zhuǎn)向架帶兩輪對(duì)兩輪對(duì)之間軸距通常為2.5m一節(jié)車廂兩端的轉(zhuǎn)向架中心距約17.5m實(shí)際編組車?yán)锩總€(gè)輪對(duì)的絕對(duì)位置都會(huì)隨時(shí)間變化。擴(kuò)展寫法是把所有輪對(duì)初始位置存成數(shù)組每個(gè)輪對(duì)按同樣的速度移動(dòng)然后任意積分點(diǎn)只要落在任何一個(gè)輪對(duì)的作用區(qū)段內(nèi)就施加對(duì)應(yīng)的輪載壓強(qiáng)。具體實(shí)現(xiàn)里要注意一軸兩端輪對(duì)分布在兩根鋼軌上每根鋼軌只承擔(dān)左側(cè)或右側(cè)的輪載所以子程序施加時(shí)按鋼軌位置區(qū)分。如果模型的鋼軌編號(hào)和坐標(biāo)固定可以直接在子程序里判斷COORDS(2)或者COORDS(3)來(lái)區(qū)分是哪根鋼軌。多輪對(duì)擴(kuò)展的推薦寫法是用循環(huán)例如REAL*8 DIST(4) DATA DIST /0.0D0, 2500.0D0, 17500.0D0, 20000.0D0/ C xBase X0 V * TIME(1) F 0.0D0 C DO I 1, 4 XWHEEL X0 V*TIME(1) - DIST(I) IF (DABS(COORDS(1) - XWHEEL) .LE. LZONE/2.0D0) THEN F PRESS GOTO 100 END IF END DO C 100 CONTINUE RETURN END注意DIST數(shù)組存的是相對(duì)首輪對(duì)的偏移量xBase相當(dāng)于首輪對(duì)的當(dāng)前位置后面每個(gè)輪對(duì)的位置在此基礎(chǔ)上減去偏移量。一個(gè)常見(jiàn)的錯(cuò)誤是直接把所有輪對(duì)的絕對(duì)位置寫死這樣車一動(dòng)起來(lái)輪對(duì)之間的間距就不對(duì)了。3.4 子程序的編譯驗(yàn)證與調(diào)試技巧在Abaqus中使用子程序前先確認(rèn)Fortran編譯環(huán)境和Abaqus版本匹配。過(guò)一遍這個(gè)流程安裝Intel Fortran Compiler和Microsoft Visual Studio配置好環(huán)境變量后在命令行執(zhí)行abaqus verify -user_std如果顯示successful則說(shuō)明編譯鏈路是通的。這個(gè)驗(yàn)證步驟不要跳過(guò)否則經(jīng)常在提交任務(wù)時(shí)報(bào)一堆找不到編譯器的錯(cuò)誤浪費(fèi)時(shí)間又查不到根因。調(diào)試DLOAD子程序最直接的辦法是在子程序里把關(guān)鍵變量寫入外部文件比如每調(diào)用一次就記錄當(dāng)前節(jié)點(diǎn)坐標(biāo)、時(shí)間、計(jì)算出的XWHEEL和F值。我在實(shí)際調(diào)試中是把這些信息寫入一個(gè)文本文件然后導(dǎo)入Excel里檢查荷載帶的位置時(shí)序是否與理論值一致。這個(gè)方法雖然笨但往往幾分鐘就能定位到問(wèn)題。另外一個(gè)調(diào)試技巧是先做一個(gè)靜態(tài)驗(yàn)證把速度設(shè)為0讓輪對(duì)固定在一個(gè)位置提交一個(gè)靜力分析步看看鋼軌變形和應(yīng)力分布是否對(duì)稱合理。對(duì)稱性檢查能快速發(fā)現(xiàn)模型坐標(biāo)系錯(cuò)誤、荷載作用位置偏移等問(wèn)題比直接上動(dòng)態(tài)分析好查得多。4. 邊界條件、接觸設(shè)置與求解控制4.1 層間接觸與約束策略CRTSⅡ型軌道層間連接是建模中影響結(jié)果很大的環(huán)節(jié)。鋼軌和扣件之間、扣件和軌道板之間采用彈簧阻尼單元連接一般用Spring2/Dashpot2單元單獨(dú)建立扣件系統(tǒng)替代實(shí)際的扣件部件。這樣做的好處是可以通過(guò)調(diào)整彈簧剛度和阻尼參數(shù)來(lái)模擬不同扣件類型比如WJ-8B和WJ-7型差異就直接改參數(shù)不需要重新建模。軌道板與CA砂漿層、CA砂漿層與底座板之間可以采用綁定約束Tie來(lái)簡(jiǎn)化處理。但如果研究目標(biāo)是軌道板與砂漿層的離縫損傷就必須用帶損傷本構(gòu)的界面單元或者面面接觸這樣才能模擬層間拉應(yīng)力超過(guò)粘結(jié)強(qiáng)度后的脫開(kāi)行為。這個(gè)選擇取決于你的研究目標(biāo)不要盲目追求精細(xì)。用Tie約束時(shí)要注意主面和從面的網(wǎng)格密度協(xié)調(diào)從面網(wǎng)格應(yīng)比主面細(xì)一些或至少相當(dāng)否則約束面上會(huì)出現(xiàn)應(yīng)力集中和偽振蕩。CA砂漿層本身是薄弱層在Tie處理后雖然不會(huì)脫開(kāi)但應(yīng)力結(jié)果相對(duì)均勻適合做整體響應(yīng)分析。4.2 邊界條件的合理截?cái)嘬壍澜Y(jié)構(gòu)的縱向尺度遠(yuǎn)遠(yuǎn)大于建模范圍如果直接把有限長(zhǎng)度的模型兩端約束死會(huì)產(chǎn)生嚴(yán)重的邊界效應(yīng)移動(dòng)荷載接近端部時(shí)結(jié)果失真。正確處理辦法是采用半無(wú)限域近似或者黏彈性邊界。最簡(jiǎn)單的方案是把底座板底面固結(jié)在長(zhǎng)度方向兩端外側(cè)再加一段過(guò)渡底座板并在端面施加彈性地基彈簧來(lái)模擬周圍土體與相鄰結(jié)構(gòu)的約束作用。在動(dòng)力學(xué)計(jì)算中還可以在端部加黏性邊界即通過(guò)阻尼單元模擬能量的逸散避免反射波在模型里來(lái)回彈跳導(dǎo)致結(jié)果振蕩。對(duì)18m長(zhǎng)的模型我給底座板底面全部固結(jié)縱向兩端設(shè)置彈性彈簧彈簧剛度根據(jù)地基系數(shù)和等效面積估算。如果做的是具體線路評(píng)估最好按實(shí)際線下基礎(chǔ)條件來(lái)標(biāo)定這部分參數(shù)。4.3 分析步設(shè)置與求解器參數(shù)隱式分析里移動(dòng)荷載是強(qiáng)非線性輸入分析步參數(shù)設(shè)置直接影響收斂性和計(jì)算效率。建議把分析步設(shè)置成固定增量步長(zhǎng)一般取荷載帶走過(guò)一個(gè)單元長(zhǎng)度所需時(shí)間的1/5到1/10。比如網(wǎng)格尺寸20mm速度83333mm/s走過(guò)一個(gè)單元需要0.00024s增量步取2×10??到5×10??s比較合適。時(shí)間增量步太大荷載跳變劇烈容易不收斂步長(zhǎng)太小計(jì)算時(shí)間成倍增加。我實(shí)際試算下來(lái)的經(jīng)驗(yàn)是先用一個(gè)較粗的網(wǎng)格和較大的增量步跑通全流程確認(rèn)結(jié)果合理后再加密網(wǎng)格并細(xì)化步長(zhǎng)不要一上來(lái)就追求極限精度。阻尼方面需要特別注意。軌道結(jié)構(gòu)的實(shí)際阻尼遠(yuǎn)小于一般建筑結(jié)構(gòu)瑞利阻尼的Alpha和Beta參數(shù)要根據(jù)結(jié)構(gòu)自振頻率來(lái)標(biāo)定不要隨便取默認(rèn)值?？梢韵茸瞿B(tài)分析獲取軌道結(jié)構(gòu)的一階豎向彎曲頻率再用頻率值反推阻尼系數(shù)。質(zhì)量阻尼Alpha對(duì)低頻響應(yīng)影響大剛度阻尼Beta對(duì)高頻振蕩影響大給得太高會(huì)把高頻輪軌動(dòng)力響應(yīng)抹平給得太低又會(huì)出現(xiàn)數(shù)值振蕩需要反復(fù)對(duì)比。5. 常見(jiàn)報(bào)錯(cuò)與排坑實(shí)錄5.1 CPU數(shù)量超過(guò)許可限制的報(bào)錯(cuò)Abaqus 提交并行任務(wù)時(shí)報(bào)錯(cuò)“the number of cpus (20) exceeds the number of cpus available”是很多新手容易卡殼的地方搜索量也一直很高。這個(gè)錯(cuò)誤核心原因是兩種一是求解器分配的CPU數(shù)量超過(guò)了當(dāng)前許可證允許的核數(shù)二是軟件讀取的系統(tǒng)邏輯核數(shù)與實(shí)際可用核數(shù)不符。后者在Windows系統(tǒng)下比較常見(jiàn)比如虛擬機(jī)環(huán)境只分配了部分邏輯核Abaqus卻識(shí)別到了更多。處理方法在Job模塊點(diǎn)擊Edit把Parallelization里的CPU數(shù)量改小一般先設(shè)2或者4跑通流程確認(rèn)沒(méi)問(wèn)題再逐步增加。同時(shí)可以命令行執(zhí)行abaqus informationlicenses查看許可證授權(quán)的核心數(shù)。如果是虛擬機(jī)或遠(yuǎn)程桌面環(huán)境檢查系統(tǒng)CPU親和性設(shè)置是否限制了Abaqus實(shí)際可用的核心數(shù)。5.2 安裝后無(wú)桌面啟動(dòng)文件和許可證不能啟動(dòng)“Abaqus安裝后桌面沒(méi)有啟動(dòng)程序文件”和“Abaqus許可證不能啟動(dòng)”屬于安裝配置階段的高頻問(wèn)題。桌面沒(méi)有快捷方式一般不是安裝失敗而是安裝程序沒(méi)默認(rèn)創(chuàng)建圖標(biāo)。解決辦法是找到安裝目錄下的啟動(dòng)腳本比如CAE的bat文件或者Exec文件夾里的abaqus.bat直接雙擊運(yùn)行或者手動(dòng)創(chuàng)建快捷方式指向該腳本。還有一種情況是CAE啟動(dòng)時(shí)依賴的Python環(huán)境路徑配置錯(cuò)誤檢查環(huán)境變量PYTHONHOME是否被其他軟件改寫。許可證不能啟動(dòng)的原因比較多集中在幾個(gè)方向許可證服務(wù)器服務(wù)沒(méi)起來(lái)、環(huán)境變量LM_LICENSE_FILE和服務(wù)器的端口設(shè)置不對(duì)、防火墻阻斷了Abaqus License Server的通信。先確認(rèn)許可證服務(wù)已在服務(wù)管理器里啟動(dòng)再用命令行執(zhí)行abaqus licensing檢查當(dāng)前許可證狀態(tài)。如果服務(wù)器是遠(yuǎn)程的確認(rèn)客戶端環(huán)境變量里填寫的端口和主機(jī)名與服務(wù)器設(shè)置一致注意端口號(hào)必須和服務(wù)器配置的端口完全匹配。5.3 DLOAD子程序編譯與運(yùn)行期故障DLOAD子程序最常見(jiàn)的編譯錯(cuò)誤是找不到Fortran編譯器。Abaqus版本和Intel編譯器版本之間兼容性要求很強(qiáng)不匹配就會(huì)出現(xiàn)“cannot find ifort”或類似的報(bào)錯(cuò)。先執(zhí)行abaqus verify -user_std驗(yàn)證整套編譯鏈路這是最快定位問(wèn)題的方法。如果驗(yàn)證失敗對(duì)照Abaqus官方兼容性表格重新安裝匹配的編譯器版本。運(yùn)行期還有一個(gè)很隱蔽的問(wèn)題DLOAD子程序里的局部變量沒(méi)有初始化。Fortran中未初始化的局部變量在不同編譯環(huán)境下可能是隨機(jī)值導(dǎo)致F輸出異常。強(qiáng)烈建議子程序入口處把F默認(rèn)為0所有局部變量顯式賦值。這類問(wèn)題排查起來(lái)特別耗時(shí)因?yàn)槟Ｐ途W(wǎng)格、材料參數(shù)都沒(méi)問(wèn)題但荷載就是不對(duì)。5.4 荷載帶階躍導(dǎo)致的計(jì)算不收斂這個(gè)問(wèn)題前面提到過(guò)但值得專門拿出來(lái)說(shuō)。DLOAD子程序用if判斷實(shí)現(xiàn)的荷載帶邊界上是從0直接跳到滿值在隱式求解器中容易造成應(yīng)力波傳播異常和收斂迭代次數(shù)激增。我實(shí)際遇到的案例是同樣的模型和材料參數(shù)加了漸變過(guò)渡的荷載帶后計(jì)算時(shí)間縮短了一半還多而且結(jié)果更平滑。推薦做法是把判斷條件改成按相對(duì)位置計(jì)算過(guò)渡系數(shù)比如離荷載帶中心越遠(yuǎn)荷載值按線性或余弦曲線遞減到0。在子程序里用一個(gè)過(guò)渡半寬定義比如過(guò)渡區(qū)取20mm那么F PRESS × max(0, 1 - |x - xwheel - LZONE/2| / TRANSWIDTH)。注意單獨(dú)處理荷載帶前后兩個(gè)邊界不要寫死對(duì)稱邏輯否則頭尾過(guò)渡不對(duì)稱。6. 結(jié)果解讀與模型驗(yàn)證經(jīng)驗(yàn)6.1 軌道板與鋼軌動(dòng)力響應(yīng)判讀移動(dòng)荷載算完之后第一步是看鋼軌的豎向位移時(shí)程。單輪荷載下鋼軌最大動(dòng)位移一般在1到2mm量級(jí)如果速度提高后位移明顯增大且伴隨高頻振蕩說(shuō)明輪軌動(dòng)力作用增強(qiáng)結(jié)果在物理上說(shuō)得通。軌道板的彎曲應(yīng)力重點(diǎn)關(guān)注板底受拉區(qū)因?yàn)闊o(wú)砟軌道損傷最容易從板底開(kāi)裂開(kāi)始。后處理時(shí)沿軌道板縱向取幾條路徑輸出彎矩應(yīng)力的分布對(duì)比不同時(shí)刻的應(yīng)力峰值位置可以識(shí)別出列車輪載作用下軌道板的受荷循環(huán)特征。這里注意區(qū)分靜載作用和動(dòng)載沖擊作用產(chǎn)生的應(yīng)力增量動(dòng)載沖擊導(dǎo)致的應(yīng)力增幅一般在10%到30%之間如果遠(yuǎn)超這個(gè)范圍要檢查阻尼參數(shù)是否給得過(guò)大。6.2 模型驗(yàn)證的幾條關(guān)鍵指標(biāo)模型做出來(lái)對(duì)不對(duì)不能只看云圖顏色好看得有對(duì)照依據(jù)。幾條可用的驗(yàn)證路徑第一與理論解析解對(duì)比鋼軌在集中力作用下的彈性彎曲位移可以用Winkler地基梁公式估算對(duì)比有限元結(jié)果和理論值的偏差如果超過(guò)10%優(yōu)先檢查扣件剛度和網(wǎng)格密度第二與文獻(xiàn)中類似參數(shù)的CRTSⅡ型軌道實(shí)測(cè)數(shù)據(jù)對(duì)比重點(diǎn)關(guān)注軌道板加速度峰值區(qū)間和鋼軌動(dòng)位移范圍第三做收斂性驗(yàn)證用兩套不同粗細(xì)的網(wǎng)格計(jì)算同一工況如果關(guān)鍵響應(yīng)偏差在5%以內(nèi)說(shuō)明網(wǎng)格密度足夠。6.3 模型擴(kuò)展的方向與建議這套模型框架可以非常方便地?cái)U(kuò)展。想研究鋼軌波磨與輪軌力的關(guān)系可以修改子程序里的輪載表達(dá)式把車輪扁疤或軌道不平順的影響加進(jìn)去。想分析CA砂漿層離縫擴(kuò)展可以把砂漿層單元換成內(nèi)聚力模型給界面一個(gè)損傷起始強(qiáng)度和斷裂能配合DLOAD移動(dòng)荷載反復(fù)掃掠就能模擬疲勞累積損傷的演化過(guò)程。想考慮橋上無(wú)砟軌道則需要在底座板下方增加橋梁梁段和支座的建模計(jì)算量會(huì)再次提升但方法完全一致。我在實(shí)際做這個(gè)項(xiàng)目時(shí)最大的體會(huì)是不要把精力全放在追求模型“多精細(xì)”上而是先明確研究問(wèn)題需要的精度等級(jí)。純粹算整體動(dòng)力響應(yīng)CA砂漿層簡(jiǎn)化成Tie就能得到很好的結(jié)果要研究層間損傷就必須上內(nèi)聚力模型和精細(xì)網(wǎng)格。工具就擺在那里關(guān)鍵是舍得花時(shí)間在子程序調(diào)試和模型驗(yàn)證上這兩個(gè)環(huán)節(jié)做扎實(shí)了后面出結(jié)果和分析都是水到渠成的事。