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

AbaqusAbaqus 常用技巧總結(jié)常用技巧總結(jié) 本手冊是由 simwe 等論壇上精華帖以及本人下載的其他資料整理 由于很多資料搜集已經(jīng) 很久 而且時間有限 都沒有注明原作者 也沒有最資料進(jìn)行分類整理 見諒 建議閱讀方式 目錄目錄 ABAQUSABAQUS 常用技巧總結(jié)常用技巧總結(jié) 1 目錄目錄 2 1 對對 TIME INCREMENT 的根本理解的根本理解 4 2 ABAQUS 請問請問 MOMENT 的加載的加載 4 3 ABAQUS 計算時計算時 C 盤的臨時文件太大了 怎么改目錄 盤的臨時文件太大了 怎么改目錄 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 多個多個 INP 文件如何實現(xiàn)批處理文件如何實現(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)行計算時如何設(shè)置調(diào)用內(nèi)存量利用命令進(jìn)行計算時如何設(shè)置調(diào)用內(nèi)存量 13 16 清華大學(xué)清華大學(xué) BBS 的的 ABAQUS 精華精華 13 17 請問怎么實現(xiàn)雙曲線請問怎么實現(xiàn)雙曲線 36 18 取消坐標(biāo)系等的顯示取消坐標(biāo)系等的顯示 37 19 如何在計算中修改材料特性如何在計算中修改材料特性 37 20 輸出計算過程中的總質(zhì)量和總剛度矩陣輸出計算過程中的總質(zhì)量和總剛度矩陣 40 21 先張預(yù)應(yīng)力 先張預(yù)應(yīng)力 40 22 用戶子程序的使用用戶子程序的使用 41 23 怎樣設(shè)定用雙怎樣設(shè)定用雙 CPU 機(jī)器進(jìn)行機(jī)器進(jìn)行 ABAQUS 計算計算 41 24 中途停止正在運(yùn)算的中途停止正在運(yùn)算的 JOB 42 25 自適應(yīng)網(wǎng)格技術(shù)自適應(yīng)網(wǎng)格技術(shù) 42 26 ABAQUS 計算與內(nèi)存計算與內(nèi)存 42 27 質(zhì)量縮放質(zhì)量縮放 43 28 ABAQUS 多處理器進(jìn)行并行計算的效果研究多處理器進(jìn)行并行計算的效果研究 54 29 YAHOO 討論組摘錄討論組摘錄 CONTACT OVERCLOSURE 55 30 原創(chuàng) 無限元建立方法 希望得到加分原創(chuàng) 無限元建立方法 希望得到加分 70 31 分享分享 ABAQUS 使用問答使用問答 76 32 轉(zhuǎn)帖轉(zhuǎn)帖 ABAQUS6 4 導(dǎo)入外來模型的幾點小經(jīng)驗 導(dǎo)入外來模型的幾點小經(jīng)驗 87 33 ABAQUS 的多圖層繪圖的多圖層繪圖 88 34 子結(jié)構(gòu)子結(jié)構(gòu) 89 35 如何在不同的分析步改變材料的參數(shù)如何在不同的分析步改變材料的參數(shù) 89 36 模型的重啟動分析 模型的重啟動分析 RESTART 91 37 ABAQUS 的單位心得的單位心得 92 1 對對 time increment 的根本理解的根本理解 abaqus 的 step 里有 maximum number of increment initial increment minimum increment maximum increment 四個量許多網(wǎng)友不知怎樣設(shè)置合理 合理設(shè)置是建立在深刻理解基 礎(chǔ)上的 要理解這個問題 首先需要了解 abaqus 的計算過程和有限元計算收斂性問題 abaqus 首先用 initial 值輸入進(jìn)行疊代計算 如果計算結(jié)果收斂 則繼續(xù)以這個值代入計算下一步 如果不收斂 則自動減小時間步長 time increment 重新計算直到收斂然后計算下一步 但是如果時間步長減小到最小值 minimum 時計算結(jié)果還是不收斂 則 abaqus 將停止 計算 由此可知 maximum 值和 minimum 值分別是 abaqus 在收斂計算時時間步長的上下限 同時 total time 求和 time increment number 當(dāng)時間步長很小時 需要計算的步數(shù) number 相應(yīng)增大 電腦計算花的時間也隨之增大 因此 number 一般要設(shè)置較大值 minimum 并不 是越小越好 因為 1 number 即計算時間增大 2 abaqus 計算精度約在 10 5 當(dāng)時間步長小 于這個值 計算結(jié)果已經(jīng)沒什么意義了 有限元計算收斂性與 最小空間步長 時間步長 值有關(guān) 若 minimum 設(shè)為 10 5 還 是不收斂 可適當(dāng)減小空間步長 即把網(wǎng)格畫細(xì)點 當(dāng)然還有一些其他辦法 如果實在計 算不了 也許是模型本身有點問題 或改為顯示 explicit 計算 總而言之 maximun number 要適當(dāng)設(shè)置較大值 initial 可適當(dāng)改小 如 2 3 量級 minimum 5 量級 不要修改 maximum 值影響不大 可不改 2 moment 的加載的加載 請教各位大牛 我最近在算一個題目 一個大筒體上有三個接管端面固定 大筒體兩端加載扭距 但是不知道 如何加載 請各位高手能賜教 謝謝 1 將大筒體兩端要施加扭矩的節(jié)點分別定義為兩個 Nset left right 2 分別在大筒體兩端的圓心處定義兩個 reference node rp left rp right 3 用如下命令將兩個節(jié)點集繞 3 軸旋轉(zhuǎn)的自由度與參考點耦合起來 其他自由度度是否耦合 根據(jù)具體問題而定 KINEMATIC COUPLING REF NODE rp left left 6 6 KINEMATIC COUPLING REF NODE rp right right 6 6 4 在兩個參考點上施加繞 3 軸旋轉(zhuǎn)的彎矩 友情提醒 reference node 也有自由度 注意相應(yīng)的邊界條件 3 abaqus 計算時計算時 c 盤的臨時文件太大了 怎么改目錄 盤的臨時文件太大了 怎么改目錄 臨時目錄是 Windows 自己定義的 可以在系統(tǒng)環(huán)境變量中修改 4 CAE 中如何加預(yù)應(yīng)力中如何加預(yù)應(yīng)力 具體沒作過 看看 PRESTRESS HOLD 和 INITIAL CONDITIONS TYPE SOLUTION REBAR 這兩個命令以及 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 中計算完成后 將結(jié)果文件輸出到 fil 2 利用 hyperworks 提供的 hmabaqus exe 在安裝目錄下的 Altair hw7 0 translators 中 3 在控制臺下運(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 命令 首先需要注意下面兩個問題 1 每一個 ANSYS 的實體模型的面或體都要定義對應(yīng)的單元類型編號材料屬性 編號實常數(shù)編號單元坐標(biāo)系編號這些參數(shù)在整個分析過程中保持不變而這些 編號對應(yīng)的屬性在各個步驟中是不同的 2 網(wǎng)格劃分要滿足所有步驟的要求單 元類型必須兼容 步驟 1 創(chuàng)建實體模型 2 創(chuàng)建多個物理環(huán)境 設(shè)定一個物理環(huán)境中的單元類型材料屬性實常數(shù)坐標(biāo)系等 將這些參數(shù)的編號賦給實體模型的面或體 施加基本物理載荷和邊界條件 設(shè)定求解選項 選擇一個標(biāo)題使用 PHYSICS WRITE 命令將物理環(huán)境存入文件中 3 清楚當(dāng)前的物理環(huán)境命令是 PHYSICS CLEAR 4 重復(fù)第二步準(zhǔn)備下一個物理環(huán)境 8 材料方向與增量步材料方向與增量步 材料方向 針對各向異性材料 如板金材料 復(fù)合材料等 變形體 材料方向定義材料的某一特定方 向如纖維方向 該方向隨著變形體的移動而移動 旋轉(zhuǎn)而旋轉(zhuǎn) 輸出的場變量值都以材料 方向為參考 有時便于數(shù)據(jù)處理 分析計算結(jié)果等 abaqus 中增量步 abaqus 中把所有載荷按一定的要求分成若干載荷步 step 每一步 step 根據(jù) abaqus 自動載荷 增量 分成若干增量 increments 每一增量施加一定的載荷 然后每一增量通過若干迭代 步 iteration 進(jìn)行迭代 當(dāng)系統(tǒng)達(dá)到平衡時 迭代結(jié)束 完成一個增量 當(dāng)所有的增量都完 成后 計算結(jié)束 反之 計算可能出現(xiàn)發(fā)散 這時 可以通過采用多鐘方法 如調(diào)整放大 質(zhì)量系數(shù) 單元網(wǎng)格優(yōu)化等 調(diào)整增量大小 使計算繼續(xù)進(jìn)行 9 多個多個 inp 文件如何實現(xiàn)批處理文件如何實現(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ī)器又不是很快 一個任務(wù)要算上幾天 一開始 在算題 的時候由于 CPU 占用很大 基本上做不了別的 捉摸了好幾天 現(xiàn)在終于有點掌握 ABAQUS 的任務(wù)管理方法了 主要有兩點 這里和大家共享 我的系統(tǒng)是 WindowsXP 使用的是 ABAQUS6 4 1 ABAQUS6 4 新增了任務(wù)管理的命令 可以暫停 恢復(fù) 和終止一個正在背景運(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 的時候 windows 任務(wù)管理欄中仍會保留 standard explicit 的計算 線程 只是不再使用 CPU 資源 當(dāng)任務(wù)恢復(fù) resume 的時候繼續(xù)工作 任務(wù)終止則就像 CAE 中提交的任務(wù)的 KILL 功能類似 直接 cut 掉正在運(yùn)行的任務(wù) 不可 恢復(fù) 其實也就和在 windows 任務(wù)管理欄中強(qiáng)行終止差不多 但屬于合法操作 2 就是利用 windos 的任務(wù)管理欄了 使用它的 processes 管理功能時 不僅僅可以強(qiáng)行關(guān) 掉一個正在運(yùn)行的任務(wù) 還可以通過改變它的優(yōu)先級別來改變其對 CPU 的占用程度 共有 Realtime High AboveNormal Normal BelowNormal Low 五個可選等級 windows 默 認(rèn)的等級是 Normal 此時所有的任務(wù)都在隨機(jī)搶占 CPU 資源 一般 ABAQUS 在運(yùn)行時想 要運(yùn)行別的程序就比較困難了 特別是一些大程序 在覺得機(jī)器使用時有明顯的延遲時 就可以把 ABAQUS 任務(wù)的優(yōu)先級別設(shè)的低一些 就可以騰出 CPU 資源給別的級別高的任 務(wù)了 不用機(jī)器的時候再把 ABAQUS 任務(wù)的級別調(diào)上去 這樣就可以娛樂 工作兩不誤 了 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ā)文問過 偶有時得到的是變形前的 有時得到是變形后的 一直沒有弄明白 這個問題 偶想得到的是變形后的坐標(biāo) 2 其實 abaqus 自己就帶有相關(guān)的功能 abaqus job job 1 suspend 可以將計算掛起 如果你需要重新進(jìn)行運(yùn)算 輸入 abaqus job job 1 resume 3 如何輸出大量節(jié)點的時間歷程曲線 在環(huán)境文件 abaqus v6 env 中添加一句 max history requests 0 即可 12 后處理積分后處理積分 CAE 自帶此功能 比如對速度時程積分得到位移時程 可以這樣操作 利用 XY Data ODB History Output 將速度時程存為 V 然后 XY Data Operate on XY Data Operations 看到那個 integrate X 這個可能還是只是實現(xiàn)對時間的積分 如果對 dx dy 的積分 如何實現(xiàn)呢 比方說 Q v A 已經(jīng)知道流速的矢量分布 怎么得到流量的分布了 這也是可以的 關(guān)鍵在于你的 x y data 數(shù)據(jù)的兩軸數(shù)據(jù)是什么 其積分就是 y 數(shù)據(jù)值對 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ī) IP 地址為 192 168 0 166 主機(jī)名為 ibmlinux 3 本地計算機(jī) IP 地址為 192 168 0 18 主機(jī)名為 training 下面的討論涉及的內(nèi)容相應(yīng)改變 Step 1 分別設(shè)置本地計算機(jī)和遠(yuǎn)程計算機(jī)之間的主機(jī)名和 IP 地址對應(yīng) 1 修改本地 hosts 文件 windows 上一般位于 c windows system32 drivers etc 下 unix 上位 于 etc 下 加入遠(yuǎn)程計算機(jī)主機(jī)名和 IP 地址對應(yīng)行 如 192 168 0 166 ibmlinux 2 修改遠(yuǎn)程 hosts 文件 位于 etc 下 加入本地計算機(jī)主機(jī)名和 IP 地址對應(yīng)行 如 192 168 0 18 training Step 2 設(shè)置本地計算機(jī)對遠(yuǎn)程計算機(jī)運(yùn)行 rsh 和 rcp 命令正確 1 兩臺計算機(jī)登陸名和密碼一致 舉例而言 若遠(yuǎn)程計算機(jī)登陸名和密碼為 root root 則 設(shè)置本地計算機(jī)也用同樣的用戶名和密碼登陸 2 修改遠(yuǎn)程計算機(jī) rhosts 文件 該文件位于對應(yīng)用戶名登陸后的主目錄下 如 home root 加入本地計算機(jī)的 IP 地址使得本地計算機(jī)可以訪問到遠(yuǎn)程計算機(jī) 3 在本地計算機(jī)的提示符下鍵入 rsh ibmlinux ls 測試 rsh 是否能夠正常使用 如果可以列出遠(yuǎn)程計算機(jī)主目錄下文件 代表 rsh 和 rcp 工作 正常 Step 3 設(shè)置遠(yuǎn)程計算機(jī)該用戶的默認(rèn)登陸程序為 csh 修改 passwd 文件 位于 etc 下 如 root Ty91eFGzybEE2 0 3 usr bin csh Step 4 修改遠(yuǎn)程計算機(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 在有限變形中 其實應(yīng)該為應(yīng)力的客觀率 似乎不對吧 大變形下此 D 非彼 D 你看過黃克智的固體本構(gòu)關(guān)系這本書么 如果你從全量理論和增量理論的角度上講 那似乎第一個 Digkl 就不對 你有第一個式子么 如果有 求導(dǎo)不久完了 那應(yīng)該就是一樣的 因為全量理論 Sij DijklEkl 滿足張量求和約定 即 Stress D Strain 而在增量理論中 S D E 在有限變形中 其實應(yīng)該為應(yīng)力的客觀率 是啊 大變形下的 D 與普通意義下的 D 在構(gòu)型上是不一樣的 畢竟 D 大是變形歷史的函數(shù) 而 D 小則不是 我推導(dǎo)一種新的本構(gòu)關(guān)系 Sij Dijkl Ekl 其中 為 Jaumann 率 假設(shè)材料一開始就屈服 即屈服面為 0 想用殼單元 似乎不對吧 大變形下此 D 非彼 D 你看過黃克智的固體本構(gòu)關(guān)系這本書么 如果你從全量理論和增量理論的角度上講 那似乎第一個 Digkl 就不對 你有第一個式子么 如果有 求導(dǎo)不久完了 是啊 大變形下的 D 與普通意義下的 D 在構(gòu)型上是不一樣的 畢竟 D 大是變形歷史的函數(shù) 而 D 小則不是 我推導(dǎo)一種新的本構(gòu)關(guān)系 Sij Dijkl Ekl 其中 為 Jaumann 率 假設(shè)材料一開始就屈服 即屈服面為 0 想用殼單元 唉 別提了 問題就出在 在 abaqus 中 明明寫著可以考慮剪切效應(yīng) 可我打印出剪切力個數(shù)是 nshr 1 即只有 S12 那我的 S13 S23 就不知怎么計算 DDSDDE 5 5 無法計算 因為 ntens 3 最多只能計算 DDSDDE 3 3 你編 umat 編進(jìn)去不久行了 他讓用 NDI NSHR NTENS 表示變量 你就用這些表示變量 這樣他就可以任意的計算了呀 而不在意實際計算的變量數(shù) 唉 別提了 問題就出在 在 abaqus 中 明明寫著可以考慮剪切效應(yīng) 可我打印出剪切力個數(shù)是 nshr 1 即只有 S12 那我的 S13 S23 就不知怎么計算 DDSDDE 5 5 無法計算 因為 ntens 3 最多只能計算 DDSDDE 3 3 因為我的 UMAT 從 abaqus 傳來的變量 當(dāng)我選用殼單元時 ntens 3 ndi 2 nshr 1 這樣的話 STRESS 為 3 個 STRESS 3 不可能計算 STRESS 4 STRESS 5 同理 STRAIN DDSDDE 也存在同樣的情況 在中厚度板殼元中 MINDLIN 中 DDSDDE 為 5 5 你的 FOR 文件收到 可惜研究的對象相差太大 看不懂 不過你在文件中定義的變量有的根本就沒有用到 還有有的地方似乎是毫無用處的代碼 比如 SMISES 一段 我覺得你不應(yīng)該太急 第一步應(yīng)該保證代碼的正確性 另外看了你上面的討論 我覺得你應(yīng)該換一種單元試一試 比如用梁單元 很有可能 NDI 和 NSHR 會不一樣 不知道你是不是一定要用到殼單元上 如果必須用在殼單元上 可能 S23 S31 就必須自己在 UAMT 中輸出了 具體的做法我也不清楚 還有對 DDSDDE 我的理解也是 dS DDSDDE dE 我的意思是不管他是幾乘幾的矩陣 一般的 D 矩陣就是 6 6 你編的 umat 需要適應(yīng)于任何單元 所以當(dāng)然 6 6 的矩陣應(yīng)該給完全 當(dāng)算體的時候 如果是殼單元 他只有 ntens 3 那就等于三唄 他用哪個就是哪個 反正你的 D6 6 的矩陣放在哪了 如果是三維實體單元 ntens 6 所以你的 D 矩陣就全用上啦 記住你編的 umat 應(yīng)該適應(yīng)于所有的單元 所以不用考慮哪幾個參數(shù)的值 而且這個 stress 3 也許是人家只用到了你 D6 6 的第一行 第二行 和的四行哪 sigh 不知道我說明白了沒有 因為我的 UMAT 從 abaqus 傳來的變量 當(dāng)我選用殼單元時 ntens 3 ndi 2 nshr 1 這樣的話 STRESS 為 3 個 STRESS 3 不可能計算 STRESS 4 STRESS 5 同理 STRAIN DDSDDE 也存在同樣的情況 在中厚度板殼元中 MINDLIN 中 DDSDDE 為 5 5 對于 DDSDDE 其實就是給了一個應(yīng)力應(yīng)變的斜率 我感覺 真正計算的是利用更新的應(yīng)力 DDSDDE 我覺得就是為了確保收斂 真開心有這么多人討論 ABAQUS 特別是 umat 我國外的師兄說 Umaterial 在國外應(yīng)用得也很廣 對于 DDSDDE 其實就是給了一個應(yīng)力應(yīng)變的斜率 我感覺 真正計算的是利用更新的應(yīng)力 DDSDDE 我覺得就是為了確保收斂 真開心有這么多人討論 ABAQUS 特別是 umat 我國外的師兄說 Umaterial 在國外應(yīng)用得也很廣 各位大俠 看到我提出的 uamt 及 statev 有這么多人感興趣 實在很 高興 希望大家能把摸索到的東西介紹給大家 現(xiàn)在我的程序能通了 因沒有屈服面的概念 雖然結(jié)果是錯誤的 但總可以靜下心來 不像以前 一運(yùn)行自己的 umat 程序就 standard exe error 并且程序 dat msg 中無 任何錯誤信息 log 說 error see msg but msg 中無 現(xiàn)在我明白了 在用戶程序中 WRITE 則寫入 log write 6 則寫入 dat 等 在監(jiān)控程序時 發(fā)現(xiàn) statev 是每次迭代和增量運(yùn)行它 他就為 0 還是為上一步的值 因為我每次都為 0 懷疑有誤 還望大俠提醒 dat 等 在監(jiān)控程序時 發(fā)現(xiàn) statev 是每次迭代和增量運(yùn)行它 他就為 0 應(yīng)該不是 0 如果在每個增量步的最后你重新更新過他 輸入文件 inp 中 depvar 不要忘了寫 statev 的維數(shù)哦 還是為上一步的值 因為我每次都為 0 懷疑有誤 還望大俠提醒 你只有一個單元么 在 saturn 好好工作 的大作中提到 這我以前還不知道 應(yīng)該不是 0 如果在每個增量步的最后你重新更新過他 輸入文件 inp 中 depvar 不要忘了寫 statev 的維數(shù)哦 在 air1 無影如風(fēng) 勤奮羊羊 課題在即 的大作中提到 可是這個寫的命令有用么 對于輸出而言 你只有一個單元么 已經(jīng)寫了 定義了 statev 的維數(shù) depvar 20 看來我的錯誤不在 statev 因為我的輸出顯示在第三增量步中 statev 不為 0 但我還是有些不清楚 uamt 程序的執(zhí)行過程 比如 一板一次殼元 s4r 面內(nèi)采用縮減積分 即 1 1 厚度方向為 5 點 simpson 積分 按道理 沒一增量步的迭代步 對每一積分點都會運(yùn)行一次 umat 即上問題 每一迭代該運(yùn)行 5 次 那上面的 statev 是每次迭代都會改變 不是每次運(yùn)行 請 air1 beauty life saturn 等各位相助 看過書上例子 如果定義 statev 為 nblock 的數(shù)組 應(yīng)該對每個積分點都進(jìn)行迭代計算 從而更新變量值 寫 stateNew nblock 請問 Mikejwg 變量 statev 怎么輸出 可以象應(yīng)力應(yīng)變那樣進(jìn)行繪圖 動畫等后處理嗎 不好意思 我只在 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)模型時 如何加入?yún)^(qū)分彈性部分與塑性部分 或 非彈性部分 也就是說 在輸出時 EE LE IE PE 因為在 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)變增量 dstran stran 而彈性部分 dEEij dSij 2G dEPij dstranij dEEij 但怎么控制程序了 默認(rèn)得情況又是怎樣計算的 謝謝 我不清楚你的本構(gòu) 我大概覺得你如果是用屈服面的話 你是不是可以用 statev 記錄每一步的 EE EP 然后通過應(yīng)力偏量來計算 DEE DEP 然后自己更新 EE EP 就不與 stran 扯上關(guān)系 我以前見過用 Umat 算 J2 流動理論的時候 是這么做的 首先我的本構(gòu)不能采用屈服面的概念 一種新型的彈塑性本構(gòu)模型 我也看了那個采用 j2 流動法則的例子 我可以在 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)變的計算問題 發(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 的計算 對 我也是最近才發(fā)現(xiàn)的 peeq 怎么計算 在 umat 中計算 peeq 因用 umat 后 系統(tǒng)缺省 peeq 0 然后再 uvarm 輸出 uvarm 可以輸出到哪里 odb 可以么 另外那個 strav 就是自己定義用于 umat 的那個數(shù)組里的數(shù)可以輸出到 odb 里么 謝了 你這個成長中的大俠 在 umat 中 statev 是不能用在 output 中的 statev 只是作為一個解的 狀態(tài)變量 說來慚愧 我是在 umat 中定義 peeq peeq 的求解一般在 彈塑性力學(xué)書上有 等效塑性應(yīng)變 用 write 寫入一個臨時文件 切記 這個文件 unit 號不要與 abaqus 中的重合 因為他有一些系統(tǒng)默認(rèn) 的文件號 然后在 uvarm 中讀取 以 uvarm 輸出 因為 uvarm 可以以 odb 的形式輸出 支持 output field output history 不知誰還有更好的方法 聽說清華有人 abaqus 用的很熟了 大家過獎了 我是沒辦法 要出論文 總得學(xué)習(xí) 上午版主還在 qq 中指點了我一番 YOU CAN TRY AS FOLLOWS STEP EXTRAPOLATION NO INC 2000000 STATIC 0 001 500 0 0 00001 0 1 可以 并且在 step 里好像有選項選擇 你是手寫的程序么 新建一個 step 然后定義邊條 在上一 step 定義邊條的節(jié)點上重新定義邊條 而且 好像 是不是有個選項是 new 這樣就可以取消上一 step 的邊條設(shè)置了 在以后的載荷步里重新定義邊界條件 加參數(shù) new 比如在第二個載荷步里 STEP BOUNDARY CONDITIONS OP NEW END STEP Command 方式 abaqus job job name user source file 或者 obj file CAE 方式 在最后 submit 的時候可以選擇 source file 或者 obj file 不太明白你的意思 如果只是要得到所需要的如應(yīng)力 應(yīng)變等 那么在數(shù)據(jù)行請求 輸出就行了啊 而且個人覺得 abaqus 的后處理功能基本上也夠用了 直接輸出到 odb 文 件 中 用 CAE 或者 Viewer 后處理即可 而 el print 輸出的 dat 文件包含其他的信息 如版本號 處理內(nèi)容等 若想要直接輸出到其他后處理軟件的話 可以輸出到 result file 文件 fil 不過 這可是二進(jìn)制文件或者 ASCII 文件 還是需要處理的 或者直接用 Hypermesh 等后處理 程序打開做圖 好像是在 tool display create 再選定一個單元集或者是節(jié)點集 然后操作是 取代 還是其他等等 就是里面的兩個圓圈相交 還是想減 你試試看看 這個是要到后處理才行吧 非要計算以后形成 odb 嗎 前處理里面怎么才可以啊 Python 在 abaqus 中可以運(yùn)行 scrpit 文件 你在 CAE manual 中就 可以將其附錄中的例子存為 py 文件 然后在 CAE 的菜單中 run script 就可以運(yùn)行了 BTW 請問你的是 abaqus 的那個版本 6 2 or 6 3 6 3 的有盜版了 1 請問如何得到 M C K 矩陣 用什么命令 試一試 element matrix output 2 想要用 Newmark 方法求解 求解問題 在 ABAQUS 中如何實現(xiàn) 或者 那里能找到類似的例子 看看 standard manuals analysis 的 dynamics analysis 和 explicit manuals 3 ABAQUS 中能實現(xiàn)這樣的東西嗎 就是我需要平滑 ABAQUS 產(chǎn)生的位移場 還是這個平滑只能在其他環(huán)境中完成 試試這個輸出選項 El print Position average at Nodes abaqus CAE 好像不能生成節(jié)點集 只能定義某個實體 后處理用這個實體上的節(jié)點 看到 CAE 的例題都是這樣做的 NODE SET 是 abaqus 的基本功能 也是其方便使用之處 在 CAE 中是支持的 這個例子我做過 沒有問題 我用 tool set create 居然不能選擇由 CAE 生成的節(jié)點 不是頂點 這是因為你沒有正確操作 Partition Edge 這一步 在 sunhaolan 的大作中提到 abaqus CAE 好像不能生成節(jié)點集 只能定義某個實體 后處理用這個實體上的節(jié)點 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 的時候 前三個用的是 standard 中的相關(guān)命令 在 step 4 的時候 我想用 explicit 來求解 并且定義了 amplitude 來模擬一個 位移隨時間的簡單變化 命令如下 step name nlgeom inc 1000 好像 explicit 里沒有 inc 這個參數(shù) dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 剛體參照點 2 2 0 03 end step 有問題嗎 另外 standard 的命令和 explicit 的命令能出現(xiàn)在同一個 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 的時候 前三個用的是 standard 中的相關(guān)命令 在 step 4 的時候 我想用 explicit 來求解 并且定義了 amplitude 來模擬一個 位移隨時間的簡單變化 命令如下 step name nlgeom inc 1000 好像 explicit 里沒有 inc 這個參數(shù) dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 剛體參照點 2 2 0 03 end step feature 菜單下有刪除一項 基準(zhǔn)點是一種 feature 工具欄不有個像鉛筆一樣的的東東么 我想把計算出來的應(yīng)力 應(yīng)變 位移等數(shù)據(jù)提取出來 在其他軟件里處理 請問該怎么辦 我在 dat 里沒有找到這些數(shù)據(jù)阿 res 不是文本格式 其他的好像也沒有啊 計算結(jié)果 沒有文本格式 只能自己把數(shù)據(jù)倒出來 在顯示結(jié)果的模塊中 名字忘了 就是 sketch 前面那個 在 report 菜單中 可以到處自己需要的文件 文件格式是 rpt 謝謝 但是這個為什么只有應(yīng)力和應(yīng)變呢 計算出的位移的數(shù)據(jù)怎么會沒有呢 在 santong supernova 的大作中提到 計算結(jié)果 沒有文本格式 只能自己把數(shù)據(jù)倒出來 在顯示結(jié)果的模塊中 名字忘了 就是 sketch 前面那個 在 report 菜單中 可以到處自己需要的文件 文件格式是 rpt 發(fā)信人 wjytj 賴子 信區(qū) FEA 你要在 STEP 模塊 field output 和 history output 改變默認(rèn) 請求輸出位移 U 或者在 step 中加上 NODE PRINT U 計算完成后在 DAT 文件中就可以提取節(jié)點位移 請問 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 這一選項 請問是什么原因 在 wjytj 賴子 的大作中提到 你要在 STEP 模塊 field output 和 history output 改變默認(rèn) 請求輸出位移 U 或者在 step 中加上 NODE PRINT U 計算完成后在 DAT 文件中就可以提取節(jié)點位移 發(fā)信人 Mikejwg mikejwg 信區(qū) FEA field output 和 history output 這在 viewer 模塊中很明顯 field output 是場量輸出 history output 為歷史數(shù)據(jù)輸出 會記錄場量隨時間的曲線 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 在 wsun 的大作中提到 請問 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 是場量輸出 history output 為歷史數(shù)據(jù)輸出 會記錄場量隨時間的曲 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 在 wsun 的大作中提到 請問 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 是場量輸出 history output 為歷史數(shù)據(jù)輸出 會記錄場量隨時間的曲 至于 node output 與 node print 的區(qū)別在于他們寫入不同的后綴文件 比如 odb dat 等 在 help 中有詳細(xì)介紹 file print file 文件類型選擇 png 或 tiff 在 Dshmy 的大作中提到 謝謝 一般說可以這樣理解 ddsdde 有點像我們常說的彈塑性矩陣 Dep 在 wjytj 的大作中提到 在 USER MATERIAL 的定義中 對 的偏倒數(shù) 即 DDSDDE 被稱為 material stiffness matrix 而在 UMAT 中 DDSDDE 被稱為 material s Jacobian matrix 請問 DDSDDE 和材料的切線剛度矩陣的關(guān)系是什么 二者是一個概念么 第二個 在 octopuszy 猴哥 的大作中提到 偶在 umat 中調(diào)用求主應(yīng)力函數(shù) CALL SP