abaqus輸入數據的案例
ABABQUS關于塑性材料輸入數據的轉換及輸入的理解
由于我們從實驗得到的數居成為工程應力或應變（這里說的是常得到的應力-應變）也稱作名義應力-名義應變，在ABAQUS中定義時需要按照公式轉化為真實應力和真實應變，轉化之后，所謂的塑性應變等于真是應變-彈性應變，其中彈性應變值等于真實應力除以彈性模量，然后將所得的數據輸入即可。此時屈服應力一欄仍輸入名義應力，應變一欄需要輸入轉化后的塑塑性應變。
我自己編的ansys輸入文件轉到marc輸入數據文件的APDL程序 ***
將這個材料編號為1,6~25格空白,表示參數的描述為預設值.26~30格在MARC中不使用,預設值為零.31~35格輸入0表示不進行水泥破裂(concrete cracking)的運算.
*VWRITE,YOUNG,POISSON
(E10.4,F10.4) !第四小Card定義楊氏模量，泊松比，1~10格表示楊氏模量,11~20格表示波松比( poisson ratio),21~30格表示密度,31~40格表示熱膨脹系數,41~50格表示Equivalent tensile yield stress,51~60格表示for ORNL yield criteria.
*VWRITE,chrval(Nelem)
('1 TO ',A5) !第五小card,輸入和這個材料有關的單元
*VWRITE
('geometry') !定義幾何數據
*VWRITE
() !空白行,不描述
*VWRITE,t,0,0,0,0,0,0
(F10.5,6F10.5) !第1~10給定材料的厚度,11格以后的參數請參閱Volumn C
ESEL,S,REAL,,1
*GET,ElR1,ELEM,0,NUM,MIN !所選單元的最小編號
*GET,NelemR1,ELEM,0,COUNT !得到所選單元的總單元數
ROWR=NelemR1/10 !以每行10個數據計算得到的行數（實數）
ROWI=NINT(ROWR) !
展開 混凝土CDP模型，直接輸入ABAQUS中，建議方便 ￥20
本excle簡捷易懂，只需在excle表中更改彈模以及軸心抗壓強度自動生成數據，表中列出了公式以及只需要輸入ABAQUS中的數據，十分容易上手
混凝土CDP模型.xlsx
ANSYS非線性分析MISO模型數據輸入的問題
首先來認識一下MISO，它的全名叫做多線性等向強化模型。
所謂“等向強化”，可以用鋼筋的冷拉變形硬化來類比，即達到屈服后繼續加載，出現塑形階段后，卸載，重新加載時應力屈服強度會有所提高，并且是一個方向屈服強度提高的同時，其他方向屈服強度同步提高。
MISO可以使用多線性來表示使用Von Mises屈服準則的等向強化的應力－應變曲線，它適用于比例加載的情況和大應變分析。
但是，應用這個模型有兩點是應當注意的：
1、曲線的第一個點必須與材料彈性模量相對應；
2、不允許有大于彈性模量或小于零的斜率段。
所有的關于MISO模型的報錯，也就是基于上述兩點原因，尤其是第二點。
fc=14.3
ft=1.43
tb,concr,1
tbdata,,0.5,0.95,ft,-1
tb,miso,1,,11
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.001,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
tbpt,,0.0033,fc*0.85
在上面的應力應變曲線中，最后一段是個下降段——但是MISO明明是不能有下降段的。。。
在ANSYS10.0及以前版本中，即便有下降段也可以繼續計算，但ANSYS12.0以后版本遇到下降段就無法計算了。這是因為老版本軟件只是把這個錯誤忽略掉，實際上并未解決，新版本軟件則老老實實地通知了用戶而已。
如何解決這個問題呢？
用上面的實例來說，就將最后的*0.85去掉即可，即把曲線的下降段換做水平段。
以上材料定義的案例，來自王新敏老師著
展開 
ZEMAX | OpticStudio 中如何輸入 Grid Sag 面型中的數據
輸入網格數據
網格點的所有數據必須在 OpticStudio 之外進行制表和計算，其數據格式必須存儲在特定格式的文件中，并在表面屬性中的導入（Import ）菜單欄中加載這個文件才能正確導入網格矢高面。OpticStudio 讀取的文件數據格式為：
其中第一行包含七個參數，其中 nx 和 ny（整數）為每個方向上網格點的數量且 nx 和 ny 必須大于等于5；unitflag 為數據的單位，0表示毫米，1表示厘米，2表示英寸，3表示米；delx 和 dely（浮點數）為網格點陣在 x 和 y 方向上的間距；xdec 和 ydec（浮點數）為坐標偏移量（可選項）。
剩下的 nx*ny 行的數據每行都包括四個浮點數和一個整數（可選項）。其中四個浮點數分別為 Z 矢高；矢高在 X 方向的導數；矢高在 Y 方向的導數；矢高的交叉導數 d/dxdy。第五個可選數據 nodata 用來表示之前的數據是否有效，0或空行表示有效的數據，其他整數表示無效數據。當nodata 為其他時，整行的數據將被替換為0，并且計算只考慮基礎面的矢高值。如果輸入的數據不是數字，則系統默認將該數據替換為零。
第一行的數據對應網格中左上角的點（沿 -Z 軸方向看去），即 -X 和 +Y 點。隨后的數據將從左至右定義這一行的點。定義完一行以后，將順次定義下一行（ny-1）的點直到所有的點都被定義。數據必須存儲為 ASCII 格式并且以 .DAT 為擴展名（如果用于非序列模式，則擴展名為 .GRD ）。
各項導數值用于雙三次樣條曲線插值算法中，用作不同矢高數據點之間的平滑處理。線性插值算法不會使用這些參數。
展開 如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio
經過《如何以數據的形式定義Zemax OpticStudio中的網格矢高表面》以及《如何以數據的形式定義Zemax OpticStudio中的網格相位表面》兩篇文章的鋪墊，今天讓我們來學習：
如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio
引言
本文示范了如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio。
閱讀本篇文章之前，請先參閱《如何以數據的方式定義網格矢高表面》以及《如何以數據的方式定義網格相位表面》。
本文中使用到的文件請從以下鏈接中下載：
鏈接: https://pan.baidu.com/s/12BjXyAlNeu6Ig4sLt9_TpQ 提取碼: ec47
正文
首先，ZYGO的測量數據是可以用Zemax OpticStudio的文件格式輸出的。輸出文件的后綴名為.zxgrd，我們可將其轉換為.DAT并將轉換后的文件，放置到 \Document\Zemax/Objects\Grid Files\ 路徑下。
下面是數據文件開頭部分的截圖，可以看到數據文件變成適用于網格矢高 (Grid Sag) 表面的輸入格式了。
X, Y 數據點的總數為723個，因為數據間距為0.01344，可得知這份數據的直徑為9.72 mm。第一行最后的標志 (flag) 被設定為0，這代表單位是 mm。
不過，如果數據測量的是穿透波前的話，這個數據就必須以相位數據的形式來輸入。
展開 如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio
如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio
如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio
引言
本文示范了如何將ZYGO干涉儀的測量數據輸入Zemax OpticStudio。 閱讀本篇文章之前，請先參閱 《如何以數據的方式定義網格矢高表面》 以及 《如何以數據的方式定義網格相位表面》。 本文中使用到的文件請從以下鏈接中下載：
鏈接: https://pan.baidu.com/s/12BjXyAlNeu6Ig4sLt9_TpQ
提取碼: ec47
正文
首先，ZYGO的測量數據是可以用Zemax OpticStudio的文件格式輸出的。輸出文件的后綴名為.zxgrd，我們可將其轉換為.DAT并將轉換后的文件，放置到 \Document\Zemax/Objects\Grid Files\ 路徑下。
下面是數據文件開頭部分的截圖，可以看到數據文件變成適用于網格矢高 (Grid Sag) 表面的輸入格式了。
X, Y 數據點的總數為723個，因為數據間距為0.01344，可得知這份數據的直徑為9.72 mm。第一行最后的標志 (flag) 被設定為0，這代表單位是 mm。
不過，如果數據測量的是穿透波前的話，這個數據就必須以相位數據的形式來輸入。
展開 python小程序，提取路徑下所有excel表中數據,輸入界面是基于TK語言編寫的
用python識別路徑下所有excel文件，并檢索每個excel文件中的sheet表，并匹配搜索關鍵詞，如果sheet表中含有關鍵詞則提取表中的markov數據到數組中，同時點擊數據輸出可將數據輸出到名稱為output_data_liu的excel表中。
ADAMS/Car 平順性評價指標計算及后處理數據處理方法-隨機輸入 ￥10
根據國標GB/T 4970-2009 汽車平順性試驗方法中規定了隨機輸入行駛指標的計算方法，下面我們探討如何利用ADAMS/car進行隨機輸入行駛指標計算。
首先，小編對國標GB/T 4970歸納總結，加權加速度均方根值是按震動方向并根據人體對震動頻率的敏感程度而進行加權計算的，是人體震動的評價指標。
單軸向加權加速度均方根值計算：
式中：
總加權均方根值計算：
式中：
利用總加速度均方根值進行平順性評價：
其次,我們介紹一下平順性后處理流程：
最后，我們舉例說明平順性后處理數據處理方法（軟件版本Adams 2013）。
使用軟件自帶Vehicle_full_4post_PAC2002.asy,獲得仿真結果文件命名為test。
啟動ADAMS/Postprocessor：
插入Wd（HZ_XY）、Wc（HZ_Z）文件。
分別繪制Wd（HZ_XY）、Wc（HZ_Z），曲線待用。
單擊（1）
繪制整車質心位置（classis_acceleration）縱向加速度曲線。
單擊(1)；
部分單詞翻譯
longitudinal：縱向
lateral：橫向
vertical：垂向
縱向加速度曲線縱坐標單位是g，因此需要換算單位。
單擊Math (1)；
在（2）處輸入*9.8；
單擊Apply（3）。
繪制縱向加速度自功率譜密度函數曲線。
單擊Plot—FFT。
按照下圖輸入參數，并單擊Apply。
生成加速度自功率譜密度函數曲線。
為了方便觀察，創建一個新page，并將自功率譜密度曲線復制（Ctrl+C ）+粘貼（Ctrl+V）到新page。
展開 『下載』如何將拉伸試驗所得的材料數據輸入到材料庫中
如何將拉伸試驗所得的材料數據輸入到材料庫中，建立流動應力曲線
70基于matlab的BP神經網絡多輸入單輸出數據結果預測，輸出結果包括均方根誤差，決定系數。 ￥25.9
基于matlab的BP神經網絡多輸入單輸出數據結果預測，輸出結果包括均方根誤差，決定系數。數據可更換自己的，程序已調通，可直接運行。
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
以下內容也可在我發布的文檔中免費下載
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'
展開 