abaqus讀取數(shù)據(jù)的案例
Abaqus利用region讀取結(jié)果數(shù)據(jù)
Abaqus利用region讀取結(jié)果數(shù)據(jù)
長(zhǎng)安CAE
Abaqus計(jì)算完成后，為輸出某部分模型的結(jié)果數(shù)據(jù)，需要用到region，并且配套使用的是getSubset函數(shù)。
示例一創(chuàng)建一個(gè)變量center，將一個(gè)叫做PUNCH的節(jié)點(diǎn)set賦值給這個(gè)變量，通過(guò)整個(gè)模型的位移場(chǎng)displacement得到該節(jié)點(diǎn)set的輸出數(shù)據(jù)，代碼如下：
center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']
centerDisplacement = displacement.getSubset(region=center)
centerValues = centerDisplacement.values
for v in centerValues:
print v.nodeLabel, v.data
其中的displacement是已經(jīng)創(chuàng)建好的位移場(chǎng)數(shù)據(jù)，centerDisplacement變量得到該節(jié)點(diǎn)集的所有數(shù)據(jù)，centerValues變量得到該節(jié)點(diǎn)集的數(shù)值數(shù)據(jù)，通過(guò)循環(huán)輸出節(jié)點(diǎn)號(hào)和當(dāng)前數(shù)據(jù)。
展開(kāi) Abaqus利用region讀取結(jié)果數(shù)據(jù)
單元質(zhì)心點(diǎn)
轉(zhuǎn)自公眾號(hào)——ABAQUS大世界
旨在分享，若侵即刪.
ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
問(wèn)題2
如何從ODB結(jié)果文件中讀取計(jì)算結(jié)果數(shù)據(jù)
第二個(gè)問(wèn)題是大家經(jīng)常要面對(duì)的，常規(guī)的數(shù)據(jù)提取采用ABAQUS自帶的后處理工具已經(jīng)能滿足需求了，但是如果我們需要進(jìn)行一些寫(xiě)出文本，數(shù)據(jù)篩選，數(shù)據(jù)的二次計(jì)算處理。。。，我們就需要用程序語(yǔ)言的幫助了。
針對(duì)今天的問(wèn)題——
提取結(jié)果并寫(xiě)出到文本
，如果用自帶的工具，可以導(dǎo)出abaqus.rpt文件，不過(guò)這種方式并不能滿足我們較為苛刻的要求，用python程序會(huì)的自由度會(huì)更大一些。
另外需要說(shuō)明的是，如果我們僅僅只是看看某單元或者節(jié)點(diǎn)的各種計(jì)算結(jié)果，或者簡(jiǎn)要的數(shù)據(jù)處理，臨時(shí)用一用，可以采用如下方式，不必費(fèi)心寫(xiě)復(fù)雜的程序，不值當(dāng)浪費(fèi)時(shí)間，
青春易逝，頭發(fā)寶貴
。當(dāng)然了，通過(guò)阿信這里給出的基本程序改一改，能省幾根頭發(fā)。
限于時(shí)間關(guān)系，這里只列出基本的應(yīng)力提取程序，寫(xiě)入文本比較簡(jiǎn)單，隨意百度即可實(shí)現(xiàn)，不再給出。
展開(kāi) VB讀取txt的數(shù)據(jù)
VB讀取txt的數(shù)據(jù)
1 概述
VB通常結(jié)合其他的計(jì)算工具進(jìn)行相關(guān)的二次開(kāi)發(fā)，在調(diào)用軟件進(jìn)行計(jì)算時(shí)產(chǎn)生的結(jié)果文件很多都是txt格式的，或者有時(shí)候?yàn)榱朔奖爿敵鰰?huì)在計(jì)算軟件里面有意將想要的結(jié)果輸出為txt格式的文件進(jìn)行存儲(chǔ)。
在VB里除了需要調(diào)用計(jì)算，還需要對(duì)txt的結(jié)果進(jìn)行讀取并顯示，以及做相關(guān)的計(jì)算。VB在讀取txt數(shù)據(jù)時(shí)用到的幾個(gè)主要函數(shù)是LineInput和InputString。本次以一次實(shí)例，通過(guò)VB讀取并顯示txt文件的某些數(shù)據(jù)。
需要讀取的數(shù)據(jù)如圖1紅色方框標(biāo)示。該txt文件一共100多行，需要的數(shù)據(jù)在第93行和第96行。
圖1 數(shù)據(jù)位置
2 方法
在VB里面拖曳生成一個(gè)簡(jiǎn)單的界面，如圖2，用于顯示數(shù)據(jù)，增加一個(gè)按鈕，執(zhí)行讀取和顯示操作。雙擊按鈕，進(jìn)入代碼編輯界面，首先定義變量。如下所示：
圖2 VB界面
Dim m As Integer, n As Integer, i As Integer, s(1000) As String, str(1000) As String, ss(1000) As String
m = 1
n = 1
Dim input1 As Integer
Dim output1 As Integer
input1 = 1
接著讀取txt文件，讀取時(shí)先獲得VB界面得到的文件位置，采用FileOpen函數(shù)，代碼如下：
FileOpen(1, 文件位置.Text & "\linearStressResults.txt", OpenMode.Input)
得到txt文件的行數(shù)，VB讀取txt文件是逐行逐行的讀，一般輸出的結(jié)果格式都是固定的，具體位置也是固定的，因此可以采用位置控制的方法找到需要的數(shù)據(jù)。
展開(kāi) 
fluent讀取Maxwell磁場(chǎng)數(shù)據(jù)的方法 ￥198
<p>fluent讀取Maxwell磁場(chǎng)數(shù)據(jù)的方法</p><p> 在計(jì)算磁流體的時(shí)候需要考慮磁場(chǎng)的分布，考慮流體的分布，那么fluent中的MHD模塊能夠很好的耦合兩者之間的效果，默認(rèn)的mhd是輸入磁場(chǎng)定值來(lái)確定空間的磁場(chǎng)分布，那么有時(shí)候是變換的磁場(chǎng)，那么就需要其他軟件來(lái)完成磁場(chǎng)的計(jì)算，那么Maxwell軟件作為ANSYS的磁場(chǎng)計(jì)算軟件，越來(lái)越多的得到了應(yīng)用。當(dāng)然其磁場(chǎng)的計(jì)算是比較方便的，fluent的流體計(jì)算也是比較方便的，那么兩者之間如何傳遞磁場(chǎng)數(shù)據(jù)呢？</p><p> 先看一下傳遞的效果，下面為一個(gè)線圈在空間的磁場(chǎng)分布如圖所示</p><div contenteditable="false" width="100%">
<img src="https://img.jishulink.com/202204/imgs/d6c521ffed7848ba8ae2fe02260181a5.png" title="Picture1.png" alt="Picture1.png" style="max-width: 760px; width: 371px; height: 200px;" width="371" height="200" data-mobile-src="https://img.jishulink.com/202204/imgs/d6c521ffed7848ba8ae2fe02260181a5.png?image_process=/format,webp/resize,w_371" data-pc-src="https://img.jishulink.com/202204/imgs/d6c521ffed7848ba8ae2fe02260181a5.png?
展開(kāi) 讀取csv文件中指定的行數(shù)據(jù)-轉(zhuǎn)載
######創(chuàng)建csv文件，并按行讀取數(shù)據(jù)#######
import numpy as np
import pandas as pd
import csv
# 1.創(chuàng)建csv文件
with open('123.csv','w',newline='') as file:
csv_writer=csv.writer(file)
# 2.讀取csv指定行的數(shù)據(jù)
with open('C:\\Users\\86178\\Desktop\\python\\pythonProject\\ethylene_methane.csv', 'r') as f:
reader = csv.reader(f)
index = 0
for row in reader:
index = index + 1
# 3.要求
if index %30000 ==0:
# 4.寫(xiě)入滿足要求的行數(shù)據(jù)
csv_writer.writerow(row)
######添加表頭#######
# 1.讀取csv文件
dg = pd.read_csv('C:\\Users\\86178\\Desktop\\python\\pythonProject\\123.csv',header=None)
# 2.添加表頭
dg = dg.rename(columns = {0:'Num', 1:'Time (seconds)', 2:'Methane conc (ppm)',3:'Ethylene conc(ppm)',4:'Sensor1
展開(kāi) Python讀取Excel表格的數(shù)據(jù)并進(jìn)行冒泡排序
工程計(jì)算的時(shí)候很多數(shù)據(jù)是存放在Excel表格中的，我們可能需要通過(guò)Python去讀取Excel表格中的數(shù)據(jù)，并作相關(guān)的計(jì)算，如何按照想要的方式讀取這些數(shù)據(jù)就很重要，這里整理一下通過(guò)Python的pandans去獲取表格某一列或者某一行的數(shù)據(jù)的方式。
首先，給出如下的一個(gè)表格，文件的位置在C:\Users\workstation\Desktop\234.xlsx，接下來(lái)就要按照
然后再Python中導(dǎo)入pandans
import pandas as pd
通過(guò)如下函數(shù)讀取表格
df = pd.read_excel(r"C:\Users\workstation\Desktop\234.xlsx",sheet_name=0)，其中sheet_name=0表示讀取sheet1，這樣整個(gè)表格就讀取進(jìn)來(lái)了。
年齡 身高
0 2 4
1 5 7
2 7 10
如何獲取年齡這一列，只需要如下的代碼df['年齡'].values就可獲取這一列除開(kāi)表頭的所有數(shù)據(jù)，同樣的身高也可以通過(guò)df['身高'].values獲取.
0 2
1 5
2 7
Name: 年齡, dtype: int64
獲取某一行要也是類似的df.values[0],就可以獲取相應(yīng)行的數(shù)值
[2 4]
這樣就可以很輕松的導(dǎo)入數(shù)據(jù)并進(jìn)行處理了。
展開(kāi) 國(guó)產(chǎn)Ｓ７２００配合485讀取數(shù)據(jù)時(shí)出現(xiàn)的一些問(wèn)題
1、請(qǐng)問(wèn)一臺(tái)PLC是否可以通過(guò)一個(gè)RS485接口接工控機(jī)發(fā)送來(lái)的三路數(shù)據(jù)信號(hào)？工控機(jī)上裝了三個(gè)數(shù)據(jù)采集終端。如果可以的話請(qǐng)問(wèn)怎樣區(qū)分三路信號(hào)？
答： PLC可以通過(guò)485接工控機(jī)發(fā)的出來(lái)的數(shù)據(jù)來(lái)讀取數(shù)據(jù)采集終端的數(shù)據(jù)，通過(guò)自由口通信的RCV命令接收工控機(jī)發(fā)出來(lái)的數(shù)據(jù)。
2、為什么RS485不能完全替代RS232呢？我們的觸摸屏和PLC通訊或PC與PLC通訊都是要用RS232為什么不能都用RS485呢？資料上說(shuō)485的優(yōu)點(diǎn)比232多多啦！為什么不能完全替代呢？
答：因?yàn)槭紫?32是計(jì)算機(jī)的標(biāo)準(zhǔn)信號(hào)，是一個(gè)業(yè)界標(biāo)準(zhǔn)，存在范圍廣泛；其次232可以實(shí)現(xiàn)全雙工通信，在發(fā)送數(shù)據(jù)的同時(shí)可以接收，可以用對(duì)講機(jī)和電話來(lái)理解，485是半雙工通信，相當(dāng)于對(duì)講機(jī)，一個(gè)時(shí)間只能一個(gè)講話，232是電話機(jī)方式，雙方都可以同時(shí)說(shuō)話。
3、我現(xiàn)在用的筆記本電腦裝有MCGS通用軟件，還有一臺(tái)西門(mén)子200PLC，數(shù)據(jù)線有USB轉(zhuǎn)232，還有232轉(zhuǎn)485.另外還有一個(gè)PPI編程電纜，我用自由口監(jiān)控200時(shí)為什么數(shù)據(jù)監(jiān)測(cè)不到呢，讀寫(xiě)都不行。想請(qǐng)問(wèn)高手我應(yīng)該如何設(shè)置呢。我想監(jiān)視和控制輸出點(diǎn)
答：485是基于232串口通信標(biāo)準(zhǔn)衍生出來(lái)的一種遠(yuǎn)距離，半雙工通信的硬件標(biāo)準(zhǔn)，相對(duì)于232標(biāo)準(zhǔn)來(lái)說(shuō)，485采用差分信號(hào)進(jìn)行收發(fā)，可以得到更好的干擾抑制效果，所以通信距離更遠(yuǎn)，理想狀態(tài)可以達(dá)到1.2KM的通信距離，另外485適用于多設(shè)備組成總線使用;由于其是半雙工的工作模式，所以在發(fā)送時(shí)不可接收,接收時(shí)不可發(fā)送,而232信號(hào)屬于全雙工,接收的同時(shí)也可以發(fā)送。自由口做數(shù)據(jù)監(jiān)視，這個(gè)描述表達(dá)本身是缺少很多條件或者問(wèn)題不清晰的， 存在問(wèn)題：a、使用編程軟件監(jiān)測(cè)還是使用MCGS監(jiān)測(cè)？
展開(kāi) 電子產(chǎn)品可靠性仿真 (ECAD數(shù)據(jù)直接讀取、跌落、熱、密封、斷裂等)
培訓(xùn)時(shí)間：
2016年6月7日
14:00 - 15:00
電路的集成規(guī)模越來(lái)越大，I/O數(shù)越來(lái)越多，使得PCB互連密度不斷加大，隨之帶來(lái)許多PCB可靠性問(wèn)題。ANSYS 17.0版本集成針對(duì)PCB的Trace Mapping強(qiáng)大功能，可以快速?gòu)腅ACD中直接導(dǎo)入PCB熱物參數(shù)，從而能在Mechanical中進(jìn)行準(zhǔn)確的PCB板熱力、疲勞、隨機(jī)振動(dòng)、跌落等可靠性問(wèn)題的仿真。本此網(wǎng)絡(luò)培訓(xùn)將介紹Trace Mapping功能，并演示ANSYS解決PCB板可靠性問(wèn)題的案例。
報(bào)名方式
PC端報(bào)名：
在瀏覽器中輸入
http://www.ansys.com/zh-cn/About-ANSYS/Events
在選擇您需要參加的網(wǎng)絡(luò)培訓(xùn)即可
微信端一鍵報(bào)名：
微信已綁定微信的用戶一鍵報(bào)名：
打開(kāi)ANSYS公眾號(hào)，點(diǎn)擊下面的菜單：
“最新活動(dòng)“點(diǎn)擊“活動(dòng)報(bào)名”，選擇活動(dòng)參加報(bào)名即可。
未綁定微信用戶的報(bào)名方式：
1).關(guān)注ANSYS官方微信
2).點(diǎn)擊進(jìn)入到ANSYS微信，點(diǎn)擊“咨詢反饋”-“注冊(cè)綁定”
3).點(diǎn)擊”最新活動(dòng)“-“網(wǎng)絡(luò)培訓(xùn)”，選擇活動(dòng)參加報(bào)名即可。
展開(kāi) Abaqus二次開(kāi)發(fā)讀取變形后的節(jié)點(diǎn)坐標(biāo)并輸出到txt文件中。讀取模態(tài)頻率到txt 文件中。 ￥10
有時(shí)候在abaqus中，我們需要知道某一個(gè)集合中的節(jié)點(diǎn)變形后的坐標(biāo)。以此為輸入來(lái)進(jìn)行一些研究。這里我用一個(gè)自編函數(shù)，將節(jié)點(diǎn)集合變形后的坐標(biāo)寫(xiě)入到一個(gè)txt文件中。格式化的寫(xiě)入文件，方便用此坐標(biāo)來(lái)進(jìn)行計(jì)算研究。
同樣有時(shí)候，需要將計(jì)算的模態(tài)頻率值提取出來(lái)。同樣用一個(gè)函數(shù)將模態(tài)頻率提取出來(lái)，放進(jìn)txt文件中，方便后續(xù)研究。
讀取的節(jié)點(diǎn)結(jié)果如下圖所示：
Abaqus二次開(kāi)發(fā)之Python腳本讀取場(chǎng)輸出位移
本次腳本將要指定數(shù)據(jù)庫(kù)文件，并輸出指定位置（set）的位移信息。
編程思路：
寫(xiě)注釋行，盡量詳細(xì)解釋腳本語(yǔ)句，增強(qiáng)可移植性；
導(dǎo)入相應(yīng)模塊；
本次實(shí)例要輸出場(chǎng)輸出位移信息，必須訪問(wèn)分析步、幀；
創(chuàng)建對(duì)象時(shí)，分步創(chuàng)建多個(gè)變量表示對(duì)象（推薦）：
創(chuàng)建變量表示odb對(duì)象，打開(kāi)輸出數(shù)據(jù)庫(kù)；
創(chuàng)建變量表示第一個(gè)分析步；
創(chuàng)建變量表示第一個(gè)分析步的最后一幀；
創(chuàng)建變量表示節(jié)點(diǎn)集；（本次案例節(jié)點(diǎn)集只有一個(gè)節(jié)點(diǎn)）
創(chuàng)建變量表示訪問(wèn)第一個(gè)分析步最后一幀的位移U.
操作腳本如下：
結(jié)果如下：
Node label: 1000
Displacement in X direction: -7.891572e-34
Displacement in Y direction: -76.45553
Displacement: 76.4555282593
注：本文中涉及的模型文件和完整腳本文件可在Up的公眾號(hào)：易木木響叮當(dāng)，中回復(fù)“讀取場(chǎng)輸出數(shù)據(jù)”自動(dòng)獲取。
展開(kāi) 
ANSYS與Abaqus球坐標(biāo)系下的結(jié)果讀取
ANSYS與Abaqus球坐標(biāo)系下的結(jié)果讀取
1 概述
采用ANSYS和Abaqus軟件計(jì)算的結(jié)果通常默認(rèn)的結(jié)果是在總體笛卡爾坐標(biāo)系下產(chǎn)生的結(jié)果，這對(duì)于應(yīng)力或者應(yīng)變等分量的分析有時(shí)候不方便，比如對(duì)于一個(gè)圓筒體，比較關(guān)心其徑向應(yīng)力和環(huán)向應(yīng)力，而這個(gè)結(jié)果直接讀取使不可能的，需要一定的轉(zhuǎn)換。
這就是結(jié)果坐標(biāo)系轉(zhuǎn)換。
在軟件里，應(yīng)力分量表示為sx,xy,xz（ANSYS），s11,s22,s22（Abaqus），當(dāng)其轉(zhuǎn)換到柱坐標(biāo)或者球坐標(biāo)時(shí)，對(duì)應(yīng)的應(yīng)力分量就發(fā)生變化，sx和s11均表示徑向應(yīng)力。
2 ANSYS
建立一個(gè)球體模型，如圖1，加載求解，得到其總體坐標(biāo)系下的sx應(yīng)力分量。
圖1
在后處理器中，將結(jié)果坐標(biāo)系轉(zhuǎn)換為球坐標(biāo)系，采用的命令為：RSYS。查詢ANSYS幫助文檔，如圖2：
圖2 RSYS
0,1,2分別代表笛卡爾坐標(biāo)系，柱坐標(biāo)系，球坐標(biāo)系。
輸入命令：RSYS，2
顯式結(jié)果sx為圖3，此時(shí)的sx應(yīng)力分量為徑向應(yīng)力。
圖3
3 Abaqus
建立模型加載求解，得到s11應(yīng)力分量如圖4.
圖4
轉(zhuǎn)換結(jié)果坐標(biāo)系，Visualization模塊下選擇 Tools--Create Coordinate Aystem，按指定方法建立局部坐標(biāo)系，然后選擇Result-Option，選擇Transformation標(biāo)簽，User-specified，就可以看到新建立的坐標(biāo)系，選擇新建的坐標(biāo)系即可完成坐標(biāo)轉(zhuǎn)換。
如圖5，圖6
圖5
建立球坐標(biāo)系的時(shí)候根據(jù)Abaqus窗口下方的提示進(jìn)行操作。
圖6
最終轉(zhuǎn)換為徑向應(yīng)力的顯式結(jié)果，如圖7
圖7
展開(kāi) Abaqus利用腳本讀取結(jié)果并創(chuàng)建新的輸出場(chǎng)
Abaqus利用腳本讀取結(jié)果并創(chuàng)建新的輸出場(chǎng)
長(zhǎng)安CAE
Abaqus的后處理有時(shí)候采用Python腳本會(huì)比較方便，此次實(shí)例通過(guò)腳本語(yǔ)言讀取已經(jīng)存在的Odb結(jié)果文件，并且通過(guò)程序?qū)崿F(xiàn)不同載荷步下的結(jié)果場(chǎng)數(shù)據(jù)運(yùn)算，以得到新的場(chǎng)數(shù)據(jù)，并顯示其云圖。
讀取ABAQUS結(jié)果文件中的單元應(yīng)變能并輸出至excel文件的處理腳本 ￥1
腳本內(nèi)容如題。
abaqus cae數(shù)據(jù)文件損壞如何修復(fù)? abaqus data
將jnl文件復(fù)制一份，重命名后綴為.py的文件，打開(kāi)abaqus，file-->run script, 運(yùn)行一遍該py文件。但要注意你原始導(dǎo)入的幾何或孤立網(wǎng)格文件還在原來(lái)的文件夾位置下。
