
發布
注冊
/
登錄abaqus輸入角度的案例
基于ABAQUS的python不同噴丸角度的隨機彈丸噴丸插件
Python實現可變噴丸角度的隨機彈丸分布
噴丸是一種典型的強化手段,采用高速彈丸反復撞擊材料表面,使表層發生塑性變形,從而引入顯著殘余壓應力,增大零件表面硬度,細化晶粒,有利于提高材料的抗疲勞性能。
噴丸覆蓋率和噴丸角度作為噴丸工藝中最重要的工藝參數之一,直接影響噴丸工藝的加工質量。覆蓋率C表示為被噴零部件表面上的彈痕面積與零部件總面積的比值,可通過阿夫拉米公式計算:
單彈丸模擬結果
r為單個彈丸撞擊后留下的彈坑半徑;N 為彈丸個數。
ABAQUS軟件作為一款強大的非線性軟件,在沖擊領域具有重要應用。本帖主要通過ABAQUS的python二次開發實現機械噴丸強化過程,并根據單彈丸噴丸后的結果結合阿夫拉米公式近似得到一定覆蓋率下的彈丸個數。
噴丸速度與噴丸強度有關,實際實驗中以ALMEN試片的弧高值表征噴丸強度,在模擬中以100%覆蓋率下的噴丸速度表征噴丸強度,下面是ALMEN試片噴丸后的變形特征。
ALMEN試片噴丸后的位移變化
下圖為通用機械噴丸插件,可以根據覆蓋率輸入彈丸參數:彈丸個數、彈丸大小,彈丸本插件采用的是可變形體,需要材料參數,噴丸角度以90度最佳,小角度噴丸在相同噴丸強度下容易出現損傷;靶材參數:靶材大小與靶材的材料參數。其次還有分析步時間以及網格大小,分析步時間其實可以通過噴丸速度與最高彈丸位置點計算,本插件仍然可以優化,彈丸采用解析剛體減少計算量。
機械噴丸插件
下面是一些插件的案例。
展開 abaqus中旋轉角度?
abaqus中我想給個長方體兩端面固定,然后給長方體除了兩端面添加旋轉角度,這個長方體會發生變形,最后想看變形后的力?該如何做?
abaqus中旋轉角度?
abaqus中我想給長方體兩端面施加固定約束,除了長方體兩端面添加旋轉角度,看長方體發生變形,最后看變形中的力?該如何設置?
從編程角度闡述有限元最佳入門方法:以Abaqus 和OpenSees 為例
而我們常見的商業有限元軟件Abaqus和ansys 的求解模塊均由Fortran 語言編寫的,這也是為什么,我們經??吹胶芏鄬W者基于自己課題研究的需要,通過Fortran語言編寫了很多材料的umat 或者vumat的子程序,實際上就是對abaqus 的材料庫進行了擴充。但是Fortran 語言相比C在調用計算機計算資源這方面,卻遠遠不及,這也是為什么,我們發現相同的自由度規模的分析,Opensees 的速度較Abaqus 快很多。更不用說,Opensees 中擁有最為強大的材料庫和單元庫,每年都會有很多來自全球的學者開發相當的代碼來擴充OS 的求解模塊。
而對于有限元軟件的前后模塊,則對計算性能要求不高,所以通常用其他的編程語言進行GUI的設計,比如我們常用的abaqus 就是采用python 語言進行前后處理器的設計的。前處理模塊,主要功能就是實現和使用者的交互,你在窗口上的點擊交互,會激發它自動生成某些代碼,比如你在用鼠標操作abaqus的每一步,都會在對應的abaqus.py文件中輸出相應的代碼如圖1:
圖1 Abaqus.py中對應GUI 操作代碼
最終abaqus 會自動把這些由GUI出發產生的代碼轉換為通用格式ASCII .inp 文件,由這些命令流去觸發abaqus 求解模塊,使其采用相應的材料,單元,作用,算法,荷載,來形成有限元模型,并對其求解。
上述分析,也就說明,我們在建立abaqus 模型時,可以采用三種方式建立模型,一是GUI操作,二是基于python 寫模型的script, 三是直接寫.inp 文件。 但是對于上述過程,對于初學者而言,或者大多數已經使用abaqus 的用戶,如果不做二次開發的話,想必對后兩者的方法都相當陌生。
展開 
abaqus用戶材料參數輸入
abaqus安裝目錄下有一個文件夾,其中有一行代碼是控制用戶材料模塊中參數輸入的問題。
有誰知道這行代碼在哪里嗎?
基于python編程操作ABAQUS輸入文件生成PD3D單元顆粒
在空間中生成剛性顆粒(注意是剛性顆粒)有下列幾種方法:
1.修改關鍵字,構建粒子生成器模型生成隨機分布剛性顆粒
2.使用python語言直接在ABAQUS中生成顆粒,并進行剛體綁定,使其成為剛性顆粒,或者直接生成解析剛體或離散剛體。
方法1生成顆粒的隨機性較好,操作簡單。方法2直接在ABAQUS界面生成顆粒,當所需顆粒數量以萬為計量單位時,在前處理界面時就會卡死,對顯卡要求極高。因此,在僅考慮到這些弊端情況下,就已經使研究人員頭皮發麻,無從下手。
在一些特定應用場合下,比如所需顆粒數量數以萬計,我們只能采用方法1生成顆粒,但我們不僅僅是需要顆粒,還需將這些顆粒與其它模型進行耦合求解計算,這個時候粒子生成器就會有局限性。此外,考慮到顆粒在空間中排布的多樣性,比如最典型的高斯分布,那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。
我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒,其可操作性更強,我們可以不采用粒子生成器內部定義的隨機算法生成顆粒,用戶可以根據需求自定義顆粒分布算法,以契合實際工況。此外,可省去粒子生成顆粒的分析步,直接進行工況建模求解計算。
本貼只是個人興趣,只提供思路,不提供源碼,用戶需了解ABAQUS的inp文件的書寫規則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信,提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒,數量為1000。具體生成結果如下圖所示。
展開 c30~c45塑性損傷本構ABAQUS輸入 ¥1
快速查詢混凝土彈模和拉壓強度標準值.xlsx
以下內容也可在我發布的文檔中免費下載
混凝土CDP模型,直接輸入ABAQUS中,建議方便 ¥20
本excle簡捷易懂,只需在excle表中更改彈模以及軸心抗壓強度自動生成數據,表中列出了公式以及只需要輸入ABAQUS中的數據,十分容易上手
混凝土CDP模型.xlsx
abaqus Python二次開發之 交互輸入和提示框
#單輸入框 提示框 from abaqus import getInput from math import sqrt number = float(getInput('Enter a number:')) print sqrt(number) #多輸入提示框 from abaqus import getInputs fields = (('Width:','10'), ('Length:', '20'), ('Height:', '30')) length, width, height = getInputs(fields=fields, label='Specify block dimensions:', dialogTitle='Create Block', ) print length, width, height #警告提示框 from abaqus import getWarningReply, YES, NO reply = getWarningReply(message='Okay to continue?', buttons=(YES,NO)) if reply == YES: print 'YES clicked' elif reply == NO: print 'NO clicked'
展開 abaqus鋼管混凝土的本構關系輸入
各位大佬,請問abaqus里面關于鋼管混凝土的核心混凝土的本構關系怎么輸入?是按照一般的分成彈性和塑性分別輸入嗎,但是彈性模量不知道怎么取。本構關系是根據韓林海的鋼管混凝土結構計算的。
hypermesh與abaqus聯合仿真:模型導入(單component多次輸入)
一、hypermesh模型導出
單個component(hypermesh)對應part(abaqus),利用inp文件做中間文件。
在操作前注意:
將需要的網格放入hypermesh中的一個component中,比如本例中只需要三維網格,但有兩個component,因此將對應的三維網格放入對應的component中,刪除多余的網格和幾何信息。
Abaqus的命名規則有特殊要求,因此在文件導出hypermesh前需要檢查component的命名,注意命名不能有小數點、開頭必須是英文字符、不能是abaqus關鍵詞。
做好設置后按如下步驟操作:
二、abaqus導入
按component順序依次導入模型,成為abaqus中的part。
File—import-model 對應的inp文件。
依次導入后,導入幾次就產生多少個model,然后按以下操作將模型進行規整。完成操作后所有parts就匯集到了最終的一個model中
展開 
ABAQUS混凝土損傷塑性模型損傷因子對本構關系影響 附c40~c45混凝土損傷因子ABAQUS輸入
下載地址:c40~c45混凝土損傷因子ABAQUS輸入
ABAQUS實現一致激勵和多點激勵輸入的結構動力彈塑性時程分析
在7度0.15g區在罕遇地震作用下,采用位移輸入模式,采用南北向的EL-centro波,峰值加速度取值為310cm/s2,分別采用一致激勵輸入和多點激勵輸入方法,進行動力彈塑性時程分析。對時程曲線的時間步長縮短一倍,即采用時間間隔為0.01s,整體時間縮短一倍,由53.48s縮短為26.74s。
加速度時程曲線
位移時程曲線
結構模型
第600步是應力云圖
頂層邊、角節點的相對柱底的X向位移
D1初始輸入端(C1組);D2結構中部(C3組);D3結構中部(C4組);
D4最后輸入端(C6組);S1一致激勵輸入角點
展開 如何使用Abaqus輸入隨時間變化的材料屬性,是否需要編寫用戶程序?
以用Field Variable+Amplitude實現,具體看
作者:謝杏子
鏈接:https://www.zhihu.com/question/51392853/answer/126127926
來源:知乎
在Abaqus中超出定義范圍的插值都是常數。比如time<86400, FV1=0; time>2.42e+05, FV1=2. 所有插值都是同理。