abaqus 常用 技巧 總結(jié).doc

AbaqusAbaqus 常用技巧總結(jié)常用技巧總結(jié) 本手冊(cè)是由 simwe 等論壇上精華帖以及本人下載的其他資料整理 由于很多資料搜集已經(jīng) 很久 而且時(shí)間有限 都沒有注明原作者 也沒有最資料進(jìn)行分類整理 見諒 建議閱讀方式 目錄目錄 ABAQUSABAQUS 常用技巧總結(jié)常用技巧總結(jié) 1 目錄目錄 2 1 對(duì)對(duì) TIME INCREMENT 的根本理解的根本理解 4 2 ABAQUS 請(qǐng)問請(qǐng)問 MOMENT 的加載的加載 4 3 ABAQUS 計(jì)算時(shí)計(jì)算時(shí) C 盤的臨時(shí)文件太大了 怎么改目錄 盤的臨時(shí)文件太大了 怎么改目錄 5 4 CAE 中如何加預(yù)應(yīng)力中如何加預(yù)應(yīng)力 5 5 HYPERMESH 里面看到里面看到 ABAQUS 分析的結(jié)果分析的結(jié)果 5 6 X Y PLOTS 5 7 把上一次的分析結(jié)果作為下一次分析的初始條件該怎么做把上一次的分析結(jié)果作為下一次分析的初始條件該怎么做 6 8 材料方向與增量步材料方向與增量步 7 9 多個(gè)多個(gè) INP 文件如何實(shí)現(xiàn)批處理文件如何實(shí)現(xiàn)批處理 7 10 關(guān)于關(guān)于 ABAQUS 的任務(wù)管理的任務(wù)管理 8 11 關(guān)于數(shù)據(jù)的輸入輸出關(guān)于數(shù)據(jù)的輸入輸出 9 12 后處理積分后處理積分 10 13 接觸分析激活殺死接觸分析激活殺死 10 14 利用利用 QUEUE 的功能由本地機(jī)器向遠(yuǎn)程的功能由本地機(jī)器向遠(yuǎn)程 UNIX 機(jī)器提交機(jī)器提交 ABAQUS 作業(yè)的方作業(yè)的方 法法 精華精華 10 15 利用命令進(jìn)行計(jì)算時(shí)如何設(shè)置調(diào)用內(nèi)存量利用命令進(jìn)行計(jì)算時(shí)如何設(shè)置調(diào)用內(nèi)存量 13 16 清華大學(xué)清華大學(xué) BBS 的的 ABAQUS 精華精華 13 17 請(qǐng)問怎么實(shí)現(xiàn)雙曲線請(qǐng)問怎么實(shí)現(xiàn)雙曲線 36 18 取消坐標(biāo)系等的顯示取消坐標(biāo)系等的顯示 37 19 如何在計(jì)算中修改材料特性如何在計(jì)算中修改材料特性 37 20 輸出計(jì)算過程中的總質(zhì)量和總剛度矩陣輸出計(jì)算過程中的總質(zhì)量和總剛度矩陣 40 21 先張預(yù)應(yīng)力 先張預(yù)應(yīng)力 40 22 用戶子程序的使用用戶子程序的使用 41 23 怎樣設(shè)定用雙怎樣設(shè)定用雙 CPU 機(jī)器進(jìn)行機(jī)器進(jìn)行 ABAQUS 計(jì)算計(jì)算 41 24 中途停止正在運(yùn)算的中途停止正在運(yùn)算的 JOB 42 25 自適應(yīng)網(wǎng)格技術(shù)自適應(yīng)網(wǎng)格技術(shù) 42 26 ABAQUS 計(jì)算與內(nèi)存計(jì)算與內(nèi)存 42 27 質(zhì)量縮放質(zhì)量縮放 43 28 ABAQUS 多處理器進(jìn)行并行計(jì)算的效果研究多處理器進(jìn)行并行計(jì)算的效果研究 54 29 YAHOO 討論組摘錄討論組摘錄 CONTACT OVERCLOSURE 55 30 原創(chuàng) 無限元建立方法 希望得到加分原創(chuàng) 無限元建立方法 希望得到加分 70 31 分享分享 ABAQUS 使用問答使用問答 76 32 轉(zhuǎn)帖轉(zhuǎn)帖 ABAQUS6 4 導(dǎo)入外來模型的幾點(diǎn)小經(jīng)驗(yàn) 導(dǎo)入外來模型的幾點(diǎn)小經(jīng)驗(yàn) 87 33 ABAQUS 的多圖層繪圖的多圖層繪圖 88 34 子結(jié)構(gòu)子結(jié)構(gòu) 89 35 如何在不同的分析步改變材料的參數(shù)如何在不同的分析步改變材料的參數(shù) 89 36 模型的重啟動(dòng)分析 模型的重啟動(dòng)分析 RESTART 91 37 ABAQUS 的單位心得的單位心得 92 1 對(duì)對(duì) time increment 的根本理解的根本理解 abaqus 的 step 里有 maximum number of increment initial increment minimum increment maximum increment 四個(gè)量許多網(wǎng)友不知怎樣設(shè)置合理 合理設(shè)置是建立在深刻理解基 礎(chǔ)上的 要理解這個(gè)問題 首先需要了解 abaqus 的計(jì)算過程和有限元計(jì)算收斂性問題 abaqus 首先用 initial 值輸入進(jìn)行疊代計(jì)算 如果計(jì)算結(jié)果收斂 則繼續(xù)以這個(gè)值代入計(jì)算下一步 如果不收斂 則自動(dòng)減小時(shí)間步長(zhǎng) time increment 重新計(jì)算直到收斂然后計(jì)算下一步 但是如果時(shí)間步長(zhǎng)減小到最小值 minimum 時(shí)計(jì)算結(jié)果還是不收斂 則 abaqus 將停止 計(jì)算 由此可知 maximum 值和 minimum 值分別是 abaqus 在收斂計(jì)算時(shí)時(shí)間步長(zhǎng)的上下限 同時(shí) total time 求和 time increment number 當(dāng)時(shí)間步長(zhǎng)很小時(shí) 需要計(jì)算的步數(shù) number 相應(yīng)增大 電腦計(jì)算花的時(shí)間也隨之增大 因此 number 一般要設(shè)置較大值 minimum 并不 是越小越好 因?yàn)?1 number 即計(jì)算時(shí)間增大 2 abaqus 計(jì)算精度約在 10 5 當(dāng)時(shí)間步長(zhǎng)小 于這個(gè)值 計(jì)算結(jié)果已經(jīng)沒什么意義了 有限元計(jì)算收斂性與 最小空間步長(zhǎng) 時(shí)間步長(zhǎng) 值有關(guān) 若 minimum 設(shè)為 10 5 還 是不收斂 可適當(dāng)減小空間步長(zhǎng) 即把網(wǎng)格畫細(xì)點(diǎn) 當(dāng)然還有一些其他辦法 如果實(shí)在計(jì) 算不了 也許是模型本身有點(diǎn)問題 或改為顯示 explicit 計(jì)算 總而言之 maximun number 要適當(dāng)設(shè)置較大值 initial 可適當(dāng)改小 如 2 3 量級(jí) minimum 5 量級(jí) 不要修改 maximum 值影響不大 可不改 2 moment 的加載的加載 請(qǐng)教各位大牛 我最近在算一個(gè)題目 一個(gè)大筒體上有三個(gè)接管端面固定 大筒體兩端加載扭距 但是不知道 如何加載 請(qǐng)各位高手能賜教 謝謝 1 將大筒體兩端要施加扭矩的節(jié)點(diǎn)分別定義為兩個(gè) Nset left right 2 分別在大筒體兩端的圓心處定義兩個(gè) reference node rp left rp right 3 用如下命令將兩個(gè)節(jié)點(diǎn)集繞 3 軸旋轉(zhuǎn)的自由度與參考點(diǎn)耦合起來 其他自由度度是否耦合 根據(jù)具體問題而定 KINEMATIC COUPLING REF NODE rp left left 6 6 KINEMATIC COUPLING REF NODE rp right right 6 6 4 在兩個(gè)參考點(diǎn)上施加繞 3 軸旋轉(zhuǎn)的彎矩 友情提醒 reference node 也有自由度 注意相應(yīng)的邊界條件 3 abaqus 計(jì)算時(shí)計(jì)算時(shí) c 盤的臨時(shí)文件太大了 怎么改目錄 盤的臨時(shí)文件太大了 怎么改目錄 臨時(shí)目錄是 Windows 自己定義的 可以在系統(tǒng)環(huán)境變量中修改 4 CAE 中如何加預(yù)應(yīng)力中如何加預(yù)應(yīng)力 具體沒作過 看看 PRESTRESS HOLD 和 INITIAL CONDITIONS TYPE SOLUTION REBAR 這兩個(gè)命令以及 ABAQUS Analysis User s Manual Defining reinforcement Section 2 2 3 Defining rebar as an element property Section 2 2 4 這兩節(jié)吧 5 hypermesh 里面看里面看 abaqus 分析的結(jié)果分析的結(jié)果 你在 abaqus 中計(jì)算完成后 將結(jié)果文件輸出到 fil 2 利用 hyperworks 提供的 hmabaqus exe 在安裝目錄下的 Altair hw7 0 translators 中 3 在控制臺(tái)下運(yùn)行 hmabaqus fil res 執(zhí)行完成后就生成了相應(yīng)的 res 文件 4 在 hyperview 中打開你的模型文件 inp 和結(jié)果文件 res 就可以查看你的結(jié)果了 gook luck 6 X Y Plots STEP pre load Step name pre load nlgeom pre loading Static 0 01 1 1e 05 0 1 LOADS Name pt load Type Concentrated force Cload G5 2 200 E6 Output history frequency 1 node output nset G5 CF2 U2 element output elset G5 E22 S22 monitor node G5 dof 2 7 如何把上一次分析結(jié)果作為下一次分析的初始條件如何把上一次分析結(jié)果作為下一次分析的初始條件 Q 如果想把上一次的分析結(jié)果作為下一次分析的初始條件該怎么做 A 使用 LDREAD 命令 首先需要注意下面兩個(gè)問題 1 每一個(gè) ANSYS 的實(shí)體模型的面或體都要定義對(duì)應(yīng)的單元類型編號(hào)材料屬性 編號(hào)實(shí)常數(shù)編號(hào)單元坐標(biāo)系編號(hào)這些參數(shù)在整個(gè)分析過程中保持不變而這些 編號(hào)對(duì)應(yīng)的屬性在各個(gè)步驟中是不同的 2 網(wǎng)格劃分要滿足所有步驟的要求單 元類型必須兼容 步驟 1 創(chuàng)建實(shí)體模型 2 創(chuàng)建多個(gè)物理環(huán)境 設(shè)定一個(gè)物理環(huán)境中的單元類型材料屬性實(shí)常數(shù)坐標(biāo)系等 將這些參數(shù)的編號(hào)賦給實(shí)體模型的面或體 施加基本物理載荷和邊界條件 設(shè)定求解選項(xiàng) 選擇一個(gè)標(biāo)題使用 PHYSICS WRITE 命令將物理環(huán)境存入文件中 3 清楚當(dāng)前的物理環(huán)境命令是 PHYSICS CLEAR 4 重復(fù)第二步準(zhǔn)備下一個(gè)物理環(huán)境 8 材料方向與增量步材料方向與增量步 材料方向 針對(duì)各向異性材料 如板金材料 復(fù)合材料等 變形體 材料方向定義材料的某一特定方 向如纖維方向 該方向隨著變形體的移動(dòng)而移動(dòng) 旋轉(zhuǎn)而旋轉(zhuǎn) 輸出的場(chǎng)變量值都以材料 方向?yàn)閰⒖?有時(shí)便于數(shù)據(jù)處理 分析計(jì)算結(jié)果等 abaqus 中增量步 abaqus 中把所有載荷按一定的要求分成若干載荷步 step 每一步 step 根據(jù) abaqus 自動(dòng)載荷 增量 分成若干增量 increments 每一增量施加一定的載荷 然后每一增量通過若干迭代 步 iteration 進(jìn)行迭代 當(dāng)系統(tǒng)達(dá)到平衡時(shí) 迭代結(jié)束 完成一個(gè)增量 當(dāng)所有的增量都完 成后 計(jì)算結(jié)束 反之 計(jì)算可能出現(xiàn)發(fā)散 這時(shí) 可以通過采用多鐘方法 如調(diào)整放大 質(zhì)量系數(shù) 單元網(wǎng)格優(yōu)化等 調(diào)整增量大小 使計(jì)算繼續(xù)進(jìn)行 9 多個(gè)多個(gè) inp 文件如何實(shí)現(xiàn)批處理文件如何實(shí)現(xiàn)批處理 我喜歡用 fortran 寫成這樣 nodedeform f90 FUNCTIONS nodedeform Entry point of console application Example of displaying Hello World at execution time PROGRAM nodedeform PURPOSE Entry point for Hello World sample console application program nodedeform USE DFLIB implicit none result SYSTEMQQ abaqus job 3d05ctJ interactive result SYSTEMQQ abaqus job 3d05ctJmodel 2 interactive result SYSTEMQQ abaqus job test 1 interactive result SYSTEMQQ abaqus job test 2 interactive result SYSTEMQQ abaqus job test 3 interactive result SYSTEMQQ abaqus job test 4 interactive result SYSTEMQQ abaqus job test 5 interactive result SYSTEMQQ abaqus job test 6 interactive result SYSTEMQQ abaqus job test 7 interactive result SYSTEMQQ abaqus job test 8 interactive end program nodedeform 10 關(guān)于關(guān)于 ABAQUS 的任務(wù)管理的任務(wù)管理 關(guān)于 ABAQUS 的任務(wù)管理 精華 chibiko 發(fā)帖 26 積分 1 于 2004 07 03 03 56 前些日子的計(jì)算量很大 自己的機(jī)器又不是很快 一個(gè)任務(wù)要算上幾天 一開始 在算題 的時(shí)候由于 CPU 占用很大 基本上做不了別的 捉摸了好幾天 現(xiàn)在終于有點(diǎn)掌握 ABAQUS 的任務(wù)管理方法了 主要有兩點(diǎn) 這里和大家共享 我的系統(tǒng)是 WindowsXP 使用的是 ABAQUS6 4 1 ABAQUS6 4 新增了任務(wù)管理的命令 可以暫停 恢復(fù) 和終止一個(gè)正在背景運(yùn)行的任 務(wù) 方法如下 在命令行輸入并運(yùn)行 任務(wù)暫停 abaqus suspend job job name 任務(wù)恢復(fù) abaqus resume job job name 任務(wù)終止 abaqus terminate job job name 其中任務(wù)暫停 suspend 的時(shí)候 windows 任務(wù)管理欄中仍會(huì)保留 standard explicit 的計(jì)算 線程 只是不再使用 CPU 資源 當(dāng)任務(wù)恢復(fù) resume 的時(shí)候繼續(xù)工作 任務(wù)終止則就像 CAE 中提交的任務(wù)的 KILL 功能類似 直接 cut 掉正在運(yùn)行的任務(wù) 不可 恢復(fù) 其實(shí)也就和在 windows 任務(wù)管理欄中強(qiáng)行終止差不多 但屬于合法操作 2 就是利用 windos 的任務(wù)管理欄了 使用它的 processes 管理功能時(shí) 不僅僅可以強(qiáng)行關(guān) 掉一個(gè)正在運(yùn)行的任務(wù) 還可以通過改變它的優(yōu)先級(jí)別來改變其對(duì) CPU 的占用程度 共有 Realtime High AboveNormal Normal BelowNormal Low 五個(gè)可選等級(jí) windows 默 認(rèn)的等級(jí)是 Normal 此時(shí)所有的任務(wù)都在隨機(jī)搶占 CPU 資源 一般 ABAQUS 在運(yùn)行時(shí)想 要運(yùn)行別的程序就比較困難了 特別是一些大程序 在覺得機(jī)器使用時(shí)有明顯的延遲時(shí) 就可以把 ABAQUS 任務(wù)的優(yōu)先級(jí)別設(shè)的低一些 就可以騰出 CPU 資源給別的級(jí)別高的任 務(wù)了 不用機(jī)器的時(shí)候再把 ABAQUS 任務(wù)的級(jí)別調(diào)上去 這樣就可以?shī)蕵?工作兩不誤 了 You may try this method which was posted at simwe forum but not my original idea Anyway I will feel very happy even if it probably just helps you just a little bit 11 關(guān)于數(shù)據(jù)的輸入輸出關(guān)于數(shù)據(jù)的輸入輸出 1 輸出數(shù)據(jù)到 dat 文件 NODE PRINT NSET nset name FREQ 1 COORD 得到的是變形前的坐標(biāo)還是變形變形后的坐標(biāo) 偶在前面已發(fā)文問過 偶有時(shí)得到的是變形前的 有時(shí)得到是變形后的 一直沒有弄明白 這個(gè)問題 偶想得到的是變形后的坐標(biāo) 2 其實(shí) abaqus 自己就帶有相關(guān)的功能 abaqus job job 1 suspend 可以將計(jì)算掛起 如果你需要重新進(jìn)行運(yùn)算 輸入 abaqus job job 1 resume 3 如何輸出大量節(jié)點(diǎn)的時(shí)間歷程曲線 在環(huán)境文件 abaqus v6 env 中添加一句 max history requests 0 即可 12 后處理積分后處理積分 CAE 自帶此功能 比如對(duì)速度時(shí)程積分得到位移時(shí)程 可以這樣操作 利用 XY Data ODB History Output 將速度時(shí)程存為 V 然后 XY Data Operate on XY Data Operations 看到那個(gè) integrate X 這個(gè)可能還是只是實(shí)現(xiàn)對(duì)時(shí)間的積分 如果對(duì) dx dy 的積分 如何實(shí)現(xiàn)呢 比方說 Q v A 已經(jīng)知道流速的矢量分布 怎么得到流量的分布了 這也是可以的 關(guān)鍵在于你的 x y data 數(shù)據(jù)的兩軸數(shù)據(jù)是什么 其積分就是 y 數(shù)據(jù)值對(duì) x 值的積分 看看 CAE user manual 33 4 4Overview of X Y data operations 13 接觸分析激活殺死接觸分析激活殺死 Usage Use both of the following options MODEL CHANGE TYPE CONTACT PAIR ADD surface 1 surface 2 CONTACT INTERFERENCE surface 1 surface 2 Usage MODEL CHANGE TYPE CONTACT PAIR REMOVE 14 利用利用 queue 的功能由本地機(jī)器向遠(yuǎn)程的功能由本地機(jī)器向遠(yuǎn)程 unix 機(jī)器提交機(jī)器提交 abaqus 作業(yè)的方法作業(yè)的方法 精華精華 利用 queue 的功能由本地機(jī)器向遠(yuǎn)程 unix 機(jī)器提交 abaqus 作業(yè)的方法 假設(shè) 1 遠(yuǎn)程 unix 機(jī)器的安裝目錄為 abaqus 2 遠(yuǎn)程計(jì)算機(jī) IP 地址為 192 168 0 166 主機(jī)名為 ibmlinux 3 本地計(jì)算機(jī) IP 地址為 192 168 0 18 主機(jī)名為 training 下面的討論涉及的內(nèi)容相應(yīng)改變 Step 1 分別設(shè)置本地計(jì)算機(jī)和遠(yuǎn)程計(jì)算機(jī)之間的主機(jī)名和 IP 地址對(duì)應(yīng) 1 修改本地 hosts 文件 windows 上一般位于 c windows system32 drivers etc 下 unix 上位 于 etc 下 加入遠(yuǎn)程計(jì)算機(jī)主機(jī)名和 IP 地址對(duì)應(yīng)行 如 192 168 0 166 ibmlinux 2 修改遠(yuǎn)程 hosts 文件 位于 etc 下 加入本地計(jì)算機(jī)主機(jī)名和 IP 地址對(duì)應(yīng)行 如 192 168 0 18 training Step 2 設(shè)置本地計(jì)算機(jī)對(duì)遠(yuǎn)程計(jì)算機(jī)運(yùn)行 rsh 和 rcp 命令正確 1 兩臺(tái)計(jì)算機(jī)登陸名和密碼一致 舉例而言 若遠(yuǎn)程計(jì)算機(jī)登陸名和密碼為 root root 則 設(shè)置本地計(jì)算機(jī)也用同樣的用戶名和密碼登陸 2 修改遠(yuǎn)程計(jì)算機(jī) rhosts 文件 該文件位于對(duì)應(yīng)用戶名登陸后的主目錄下 如 home root 加入本地計(jì)算機(jī)的 IP 地址使得本地計(jì)算機(jī)可以訪問到遠(yuǎn)程計(jì)算機(jī) 3 在本地計(jì)算機(jī)的提示符下鍵入 rsh ibmlinux ls 測(cè)試 rsh 是否能夠正常使用 如果可以列出遠(yuǎn)程計(jì)算機(jī)主目錄下文件 代表 rsh 和 rcp 工作 正常 Step 3 設(shè)置遠(yuǎn)程計(jì)算機(jī)該用戶的默認(rèn)登陸程序?yàn)?csh 修改 passwd 文件 位于 etc 下 如 root Ty91eFGzybEE2 0 3 usr bin csh Step 4 修改遠(yuǎn)程計(jì)算機(jī)的 abaqus 環(huán)境文件 abaqus v6 env 在環(huán)境文件最后加入 queue name nice queue after prefix b queue prefix q nice queue bin sh c nice S 1 L 2 而在增量理論中 S D E 在有限變形中 其實(shí)應(yīng)該為應(yīng)力的客觀率 似乎不對(duì)吧 大變形下此 D 非彼 D 你看過黃克智的固體本構(gòu)關(guān)系這本書么 如果你從全量理論和增量理論的角度上講 那似乎第一個(gè) Digkl 就不對(duì) 你有第一個(gè)式子么 如果有 求導(dǎo)不久完了 那應(yīng)該就是一樣的 因?yàn)槿坷碚?Sij DijklEkl 滿足張量求和約定 即 Stress D Strain 而在增量理論中 S D E 在有限變形中 其實(shí)應(yīng)該為應(yīng)力的客觀率 是啊 大變形下的 D 與普通意義下的 D 在構(gòu)型上是不一樣的 畢竟 D 大是變形歷史的函數(shù) 而 D 小則不是 我推導(dǎo)一種新的本構(gòu)關(guān)系 Sij Dijkl Ekl 其中 為 Jaumann 率 假設(shè)材料一開始就屈服 即屈服面為 0 想用殼單元 似乎不對(duì)吧 大變形下此 D 非彼 D 你看過黃克智的固體本構(gòu)關(guān)系這本書么 如果你從全量理論和增量理論的角度上講 那似乎第一個(gè) Digkl 就不對(duì) 你有第一個(gè)式子么 如果有 求導(dǎo)不久完了 是啊 大變形下的 D 與普通意義下的 D 在構(gòu)型上是不一樣的 畢竟 D 大是變形歷史的函數(shù) 而 D 小則不是 我推導(dǎo)一種新的本構(gòu)關(guān)系 Sij Dijkl Ekl 其中 為 Jaumann 率 假設(shè)材料一開始就屈服 即屈服面為 0 想用殼單元 唉 別提了 問題就出在 在 abaqus 中 明明寫著可以考慮剪切效應(yīng) 可我打印出剪切力個(gè)數(shù)是 nshr 1 即只有 S12 那我的 S13 S23 就不知怎么計(jì)算 DDSDDE 5 5 無法計(jì)算 因?yàn)?ntens 3 最多只能計(jì)算 DDSDDE 3 3 你編 umat 編進(jìn)去不久行了 他讓用 NDI NSHR NTENS 表示變量 你就用這些表示變量 這樣他就可以任意的計(jì)算了呀 而不在意實(shí)際計(jì)算的變量數(shù) 唉 別提了 問題就出在 在 abaqus 中 明明寫著可以考慮剪切效應(yīng) 可我打印出剪切力個(gè)數(shù)是 nshr 1 即只有 S12 那我的 S13 S23 就不知怎么計(jì)算 DDSDDE 5 5 無法計(jì)算 因?yàn)?ntens 3 最多只能計(jì)算 DDSDDE 3 3 因?yàn)槲业?UMAT 從 abaqus 傳來的變量 當(dāng)我選用殼單元時(shí) ntens 3 ndi 2 nshr 1 這樣的話 STRESS 為 3 個(gè) STRESS 3 不可能計(jì)算 STRESS 4 STRESS 5 同理 STRAIN DDSDDE 也存在同樣的情況 在中厚度板殼元中 MINDLIN 中 DDSDDE 為 5 5 你的 FOR 文件收到 可惜研究的對(duì)象相差太大 看不懂 不過你在文件中定義的變量有的根本就沒有用到 還有有的地方似乎是毫無用處的代碼 比如 SMISES 一段 我覺得你不應(yīng)該太急 第一步應(yīng)該保證代碼的正確性 另外看了你上面的討論 我覺得你應(yīng)該換一種單元試一試 比如用梁?jiǎn)卧?很有可能 NDI 和 NSHR 會(huì)不一樣 不知道你是不是一定要用到殼單元上 如果必須用在殼單元上 可能 S23 S31 就必須自己在 UAMT 中輸出了 具體的做法我也不清楚 還有對(duì) DDSDDE 我的理解也是 dS DDSDDE dE 我的意思是不管他是幾乘幾的矩陣 一般的 D 矩陣就是 6 6 你編的 umat 需要適應(yīng)于任何單元 所以當(dāng)然 6 6 的矩陣應(yīng)該給完全 當(dāng)算體的時(shí)候 如果是殼單元 他只有 ntens 3 那就等于三唄 他用哪個(gè)就是哪個(gè) 反正你的 D6 6 的矩陣放在哪了 如果是三維實(shí)體單元 ntens 6 所以你的 D 矩陣就全用上啦 記住你編的 umat 應(yīng)該適應(yīng)于所有的單元 所以不用考慮哪幾個(gè)參數(shù)的值 而且這個(gè) stress 3 也許是人家只用到了你 D6 6 的第一行 第二行 和的四行哪 sigh 不知道我說明白了沒有 因?yàn)槲业?UMAT 從 abaqus 傳來的變量 當(dāng)我選用殼單元時(shí) ntens 3 ndi 2 nshr 1 這樣的話 STRESS 為 3 個(gè) STRESS 3 不可能計(jì)算 STRESS 4 STRESS 5 同理 STRAIN DDSDDE 也存在同樣的情況 在中厚度板殼元中 MINDLIN 中 DDSDDE 為 5 5 對(duì)于 DDSDDE 其實(shí)就是給了一個(gè)應(yīng)力應(yīng)變的斜率 我感覺 真正計(jì)算的是利用更新的應(yīng)力 DDSDDE 我覺得就是為了確保收斂 真開心有這么多人討論 ABAQUS 特別是 umat 我國(guó)外的師兄說 Umaterial 在國(guó)外應(yīng)用得也很廣 對(duì)于 DDSDDE 其實(shí)就是給了一個(gè)應(yīng)力應(yīng)變的斜率 我感覺 真正計(jì)算的是利用更新的應(yīng)力 DDSDDE 我覺得就是為了確保收斂 真開心有這么多人討論 ABAQUS 特別是 umat 我國(guó)外的師兄說 Umaterial 在國(guó)外應(yīng)用得也很廣 各位大俠 看到我提出的 uamt 及 statev 有這么多人感興趣 實(shí)在很 高興 希望大家能把摸索到的東西介紹給大家 現(xiàn)在我的程序能通了 因沒有屈服面的概念 雖然結(jié)果是錯(cuò)誤的 但總可以靜下心來 不像以前 一運(yùn)行自己的 umat 程序就 standard exe error 并且程序 dat msg 中無 任何錯(cuò)誤信息 log 說 error see msg but msg 中無 現(xiàn)在我明白了 在用戶程序中 WRITE 則寫入 log write 6 則寫入 dat 等 在監(jiān)控程序時(shí) 發(fā)現(xiàn) statev 是每次迭代和增量運(yùn)行它 他就為 0 還是為上一步的值 因?yàn)槲颐看味紴?0 懷疑有誤 還望大俠提醒 dat 等 在監(jiān)控程序時(shí) 發(fā)現(xiàn) statev 是每次迭代和增量運(yùn)行它 他就為 0 應(yīng)該不是 0 如果在每個(gè)增量步的最后你重新更新過他 輸入文件 inp 中 depvar 不要忘了寫 statev 的維數(shù)哦 還是為上一步的值 因?yàn)槲颐看味紴?0 懷疑有誤 還望大俠提醒 你只有一個(gè)單元么 在 saturn 好好工作 的大作中提到 這我以前還不知道 應(yīng)該不是 0 如果在每個(gè)增量步的最后你重新更新過他 輸入文件 inp 中 depvar 不要忘了寫 statev 的維數(shù)哦 在 air1 無影如風(fēng) 勤奮羊羊 課題在即 的大作中提到 可是這個(gè)寫的命令有用么 對(duì)于輸出而言 你只有一個(gè)單元么 已經(jīng)寫了 定義了 statev 的維數(shù) depvar 20 看來我的錯(cuò)誤不在 statev 因?yàn)槲业妮敵鲲@示在第三增量步中 statev 不為 0 但我還是有些不清楚 uamt 程序的執(zhí)行過程 比如 一板一次殼元 s4r 面內(nèi)采用縮減積分 即 1 1 厚度方向?yàn)?5 點(diǎn) simpson 積分 按道理 沒一增量步的迭代步 對(duì)每一積分點(diǎn)都會(huì)運(yùn)行一次 umat 即上問題 每一迭代該運(yùn)行 5 次 那上面的 statev 是每次迭代都會(huì)改變 不是每次運(yùn)行 請(qǐng) air1 beauty life saturn 等各位相助 看過書上例子 如果定義 statev 為 nblock 的數(shù)組 應(yīng)該對(duì)每個(gè)積分點(diǎn)都進(jìn)行迭代計(jì)算 從而更新變量值 寫 stateNew nblock 請(qǐng)問 Mikejwg 變量 statev 怎么輸出 可以象應(yīng)力應(yīng)變那樣進(jìn)行繪圖 動(dòng)畫等后處理嗎 不好意思 我只在 log 文件中輸出 還沒有在 odb 中輸出 想同 element output 等那樣輸出 只能在 UVARM 中定義 但 同 statev 是不同的 老兄你可以試試 發(fā)信人 Mikejwg mikejwg 信區(qū) FEA 標(biāo) 題 用 abaqus 中的 UMAT 的大俠看過來 發(fā)信站 BBS 水木清華站 Tue Dec 25 15 36 27 2001 各位大俠 在加入自己彈塑性本構(gòu)模型時(shí) 如何加入?yún)^(qū)分彈性部分與塑性部分 或 非彈性部分 也就是說 在輸出時(shí) EE LE IE PE 因?yàn)樵?umat 中 我們只發(fā)現(xiàn) dstran stress 除此并沒有發(fā)現(xiàn)其它有關(guān)應(yīng)力應(yīng)變得變量 你所說的區(qū)分彈性部分和塑性部分是指應(yīng)變的彈塑性部分 E EE PE 嗎 還是指加載過程的彈塑性 E EE PE 我覺得想從 umat 的輸入?yún)?shù)里面區(qū)分是不可能的 可不可以你先算出彈性的 然后用 stran 減去彈性的 得到塑性的 我猜的哦 我自己是明白的 但不知怎么讓程序明白 比如 在一定條件下可解出位移增量 abaqus 得到應(yīng)變?cè)隽?dstran stran 而彈性部分 dEEij dSij 2G dEPij dstranij dEEij 但怎么控制程序了 默認(rèn)得情況又是怎樣計(jì)算的 謝謝 我不清楚你的本構(gòu) 我大概覺得你如果是用屈服面的話 你是不是可以用 statev 記錄每一步的 EE EP 然后通過應(yīng)力偏量來計(jì)算 DEE DEP 然后自己更新 EE EP 就不與 stran 扯上關(guān)系 我以前見過用 Umat 算 J2 流動(dòng)理論的時(shí)候 是這么做的 首先我的本構(gòu)不能采用屈服面的概念 一種新型的彈塑性本構(gòu)模型 我也看了那個(gè)采用 j2 流動(dòng)法則的例子 我可以在 umat 中的 statev 重定義 ee pe 但 statev 不能像 EE PE IE 那樣在 odb 中輸出 再說在 umat 中的 ee pe 與 inp 中的 ee pe 是不同的 謝謝 看來你只能試試 uvarm 了 呵呵 發(fā)信人 Mikejwg mikejwg 信區(qū) FEA 標(biāo) 題 Re ABAQUS 中的彈塑性應(yīng)變的計(jì)算問題 發(fā)信站 BBS 水木清華站 Thu Dec 27 14 19 27 2001 用 uvar 勉強(qiáng)成功 標(biāo) 題 在 abaqus 中 DSTRAN STRAN 為工程剪應(yīng)變嗎 Sat Dec 29 09 44 03 2001 同題 這涉及到 ddsdde peeq 的計(jì)算 對(duì) 我也是最近才發(fā)現(xiàn)的 peeq 怎么計(jì)算 在 umat 中計(jì)算 peeq 因用 umat 后 系統(tǒng)缺省 peeq 0 然后再 uvarm 輸出 uvarm 可以輸出到哪里 odb 可以么 另外那個(gè) strav 就是自己定義用于 umat 的那個(gè)數(shù)組里的數(shù)可以輸出到 odb 里么 謝了 你這個(gè)成長(zhǎng)中的大俠 在 umat 中 statev 是不能用在 output 中的 statev 只是作為一個(gè)解的 狀態(tài)變量 說來慚愧 我是在 umat 中定義 peeq peeq 的求解一般在 彈塑性力學(xué)書上有 等效塑性應(yīng)變 用 write 寫入一個(gè)臨時(shí)文件 切記 這個(gè)文件 unit 號(hào)不要與 abaqus 中的重合 因?yàn)樗幸恍┫到y(tǒng)默認(rèn) 的文件號(hào) 然后在 uvarm 中讀取 以 uvarm 輸出 因?yàn)?uvarm 可以以 odb 的形式輸出 支持 output field output history 不知誰還有更好的方法 聽說清華有人 abaqus 用的很熟了 大家過獎(jiǎng)了 我是沒辦法 要出論文 總得學(xué)習(xí) 上午版主還在 qq 中指點(diǎn)了我一番 YOU CAN TRY AS FOLLOWS STEP EXTRAPOLATION NO INC 2000000 STATIC 0 001 500 0 0 00001 0 1 可以 并且在 step 里好像有選項(xiàng)選擇 你是手寫的程序么 新建一個(gè) step 然后定義邊條 在上一 step 定義邊條的節(jié)點(diǎn)上重新定義邊條 而且 好像 是不是有個(gè)選項(xiàng)是 new 這樣就可以取消上一 step 的邊條設(shè)置了 在以后的載荷步里重新定義邊界條件 加參數(shù) new 比如在第二個(gè)載荷步里 STEP BOUNDARY CONDITIONS OP NEW END STEP Command 方式 abaqus job job name user source file 或者 obj file CAE 方式 在最后 submit 的時(shí)候可以選擇 source file 或者 obj file 不太明白你的意思 如果只是要得到所需要的如應(yīng)力 應(yīng)變等 那么在數(shù)據(jù)行請(qǐng)求 輸出就行了啊 而且個(gè)人覺得 abaqus 的后處理功能基本上也夠用了 直接輸出到 odb 文 件 中 用 CAE 或者 Viewer 后處理即可 而 el print 輸出的 dat 文件包含其他的信息 如版本號(hào) 處理內(nèi)容等 若想要直接輸出到其他后處理軟件的話 可以輸出到 result file 文件 fil 不過 這可是二進(jìn)制文件或者 ASCII 文件 還是需要處理的 或者直接用 Hypermesh 等后處理 程序打開做圖 好像是在 tool display create 再選定一個(gè)單元集或者是節(jié)點(diǎn)集 然后操作是 取代 還是其他等等 就是里面的兩個(gè)圓圈相交 還是想減 你試試看看 這個(gè)是要到后處理才行吧 非要計(jì)算以后形成 odb 嗎 前處理里面怎么才可以啊 Python 在 abaqus 中可以運(yùn)行 scrpit 文件 你在 CAE manual 中就 可以將其附錄中的例子存為 py 文件 然后在 CAE 的菜單中 run script 就可以運(yùn)行了 BTW 請(qǐng)問你的是 abaqus 的那個(gè)版本 6 2 or 6 3 6 3 的有盜版了 1 請(qǐng)問如何得到 M C K 矩陣 用什么命令 試一試 element matrix output 2 想要用 Newmark 方法求解 求解問題 在 ABAQUS 中如何實(shí)現(xiàn) 或者 那里能找到類似的例子 看看 standard manuals analysis 的 dynamics analysis 和 explicit manuals 3 ABAQUS 中能實(shí)現(xiàn)這樣的東西嗎 就是我需要平滑 ABAQUS 產(chǎn)生的位移場(chǎng) 還是這個(gè)平滑只能在其他環(huán)境中完成 試試這個(gè)輸出選項(xiàng) El print Position average at Nodes abaqus CAE 好像不能生成節(jié)點(diǎn)集 只能定義某個(gè)實(shí)體 后處理用這個(gè)實(shí)體上的節(jié)點(diǎn) 看到 CAE 的例題都是這樣做的 NODE SET 是 abaqus 的基本功能 也是其方便使用之處 在 CAE 中是支持的 這個(gè)例子我做過 沒有問題 我用 tool set create 居然不能選擇由 CAE 生成的節(jié)點(diǎn) 不是頂點(diǎn) 這是因?yàn)槟銢]有正確操作 Partition Edge 這一步 在 sunhaolan 的大作中提到 abaqus CAE 好像不能生成節(jié)點(diǎn)集 只能定義某個(gè)實(shí)體 后處理用這個(gè)實(shí)體上的節(jié)點(diǎn) ABAQUS divides all of its analysis procedures into two main groups linear perturbation and general ABAQUS makes a very clear distinction between perturbation and general analysis procedures because loading conditions and time are defined differently for the two cases Furthermore the results from each type of procedure should be interpreted differently The response of the model during a general analysis procedure known as a general step may be either nonlinear or linear In a step that uses a perturbation procedure which is called a perturbation step the response can only be linear ABAQUS treats such steps as a linear perturbation about the preloaded predeformed state known as the base state created by any previous general steps therefore its capability for doing linear simulations is rather more general than that of a purely linear analysis program 各位大俠 我在定義 step 的時(shí)候 前三個(gè)用的是 standard 中的相關(guān)命令 在 step 4 的時(shí)候 我想用 explicit 來求解 并且定義了 amplitude 來模擬一個(gè) 位移隨時(shí)間的簡(jiǎn)單變化 命令如下 step name nlgeom inc 1000 好像 explicit 里沒有 inc 這個(gè)參數(shù) dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 剛體參照點(diǎn) 2 2 0 03 end step 有問題嗎 另外 standard 的命令和 explicit 的命令能出現(xiàn)在同一個(gè) inp 文件中嗎 要轉(zhuǎn)換的 ABAQUS Standard and ABAQUS Explicit are separate program modules with different data structures hence the explicit dynamics procedure cannot be used in the same analysis as any of the procedures in ABAQUS Standard However ABAQUS provides a capability to import a deformed mesh and associated material state from ABAQUS Explicit into ABAQUS Standard and vice versa This procedure is described in Transferring results between ABAQUS Explicit and ABAQUS Standard Section 7 6 2 of the ABAQUS Standard User s Manual 在 linyuwho 蒼龍有悔 的大作中提到 各位大俠 我在定義 step 的時(shí)候 前三個(gè)用的是 standard 中的相關(guān)命令 在 step 4 的時(shí)候 我想用 explicit 來求解 并且定義了 amplitude 來模擬一個(gè) 位移隨時(shí)間的簡(jiǎn)單變化 命令如下 step name nlgeom inc 1000 好像 explicit 里沒有 inc 這個(gè)參數(shù) dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 剛體參照點(diǎn) 2 2 0 03 end step feature 菜單下有刪除一項(xiàng) 基準(zhǔn)點(diǎn)是一種 feature 工具欄不有個(gè)像鉛筆一樣的的東東么 我想把計(jì)算出來的應(yīng)力 應(yīng)變 位移等數(shù)據(jù)提取出來 在其他軟件里處理 請(qǐng)問該怎么辦 我在 dat 里沒有找到這些數(shù)據(jù)阿 res 不是文本格式 其他的好像也沒有啊 計(jì)算結(jié)果 沒有文本格式 只能自己把數(shù)據(jù)倒出來 在顯示結(jié)果的模塊中 名字忘了 就是 sketch 前面那個(gè) 在 report 菜單中 可以到處自己需要的文件 文件格式是 rpt 謝謝 但是這個(gè)為什么只有應(yīng)力和應(yīng)變呢 計(jì)算出的位移的數(shù)據(jù)怎么會(huì)沒有呢 在 santong supernova 的大作中提到 計(jì)算結(jié)果 沒有文本格式 只能自己把數(shù)據(jù)倒出來 在顯示結(jié)果的模塊中 名字忘了 就是 sketch 前面那個(gè) 在 report 菜單中 可以到處自己需要的文件 文件格式是 rpt 發(fā)信人 wjytj 賴子 信區(qū) FEA 你要在 STEP 模塊 field output 和 history output 改變默認(rèn) 請(qǐng)求輸出位移 U 或者在 step 中加上 NODE PRINT U 計(jì)算完成后在 DAT 文件中就可以提取節(jié)點(diǎn)位移 請(qǐng)問 field output 和 history output 什么區(qū)別 關(guān)鍵字 node output 和 node print 有什么區(qū)別 謝謝 我的 inp 數(shù)據(jù)文件里有這么幾行 Output field frequency 99999 Node Output U RF Element Output S E EE Output history variable PRESELECT El Print freq 999999 Node Print freq 999999 其中在 history output 中有 node print 關(guān)鍵字 但是沒有指定輸入 U 最后導(dǎo)出的數(shù)據(jù)沒有 U 這一選項(xiàng) 請(qǐng)問是什么原因 在 wjytj 賴子 的大作中提到 你要在 STEP 模塊 field output 和 history output 改變默認(rèn) 請(qǐng)求輸出位移 U 或者在 step 中加上 NODE PRINT U 計(jì)算完成后在 DAT 文件中就可以提取節(jié)點(diǎn)位移 發(fā)信人 Mikejwg mikejwg 信區(qū) FEA field output 和 history output 這在 viewer 模塊中很明顯 field output 是場(chǎng)量輸出 history output 為歷史數(shù)據(jù)輸出 會(huì)記錄場(chǎng)量隨時(shí)間的曲線 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 在 wsun 的大作中提到 請(qǐng)問 field output 和 history output 什么區(qū)別 關(guān)鍵字 node output 和 node print 有什么區(qū)別 謝謝 我的 inp 數(shù)據(jù)文件里有這么幾行 Output field frequency 99999 Node Output U RF Element Output S E EE Output history variable PRESELECT El Print freq 999999 發(fā)信人 wjytj 賴子 信區(qū) FEA 然也 補(bǔ)充一句 node output 是給 CAR 或 VIEWER 做后處理用的 是二進(jìn)制文件 而 node print 則是寫到 DAT 文件中 你可以自己看的文本文件 在 Mikejwg 的大作中提到 field output 和 history output 這在 viewer 模塊中很明顯 field output 是場(chǎng)量輸出 history output 為歷史數(shù)據(jù)輸出 會(huì)記錄場(chǎng)量隨時(shí)間的曲 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 在 wsun 的大作中提到 請(qǐng)問 field output 和 history output 什么區(qū)別 關(guān)鍵字 node output 和 node print 有什么區(qū)別 謝謝 我的 inp 數(shù)據(jù)文件里有這么幾行 Output field frequency 99999 Node Output 發(fā)信人 octopuszy 猴哥 信區(qū) FEA node output 是在 odb 文件中輸出 br node print 是在 dat 文件中輸出 在 Mikejwg mikejwg 的大作中提到 field output 和 history output 這在 viewer 模塊中很明顯 field output 是場(chǎng)量輸出 history output 為歷史數(shù)據(jù)輸出 會(huì)記錄場(chǎng)量隨時(shí)間的曲 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 file print file 文件類型選擇 png 或 tiff 在 Dshmy 的大作中提到 謝謝 一般說可以這樣理解 ddsdde 有點(diǎn)像我們常說的彈塑性矩陣 Dep 在 wjytj 的大作中提到 在 USER MATERIAL 的定義中 對(duì) 的偏倒數(shù) 即 DDSDDE 被稱為 material stiffness matrix 而在 UMAT 中 DDSDDE 被稱為 material s Jacobian matrix 請(qǐng)問 DDSDDE 和材料的切線剛度矩陣的關(guān)系是什么 二者是一個(gè)概念么 第二個(gè) 在 octopuszy 猴哥 的大作中提到 偶在 umat 中調(diào)用求主應(yīng)力函數(shù) CALL SP