ansys單元坐標(biāo)提取
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-03-07
ansys單元坐標(biāo)提取的實(shí)例教程
然后使用*vget讀取節(jié)點(diǎn)編號(hào)及相應(yīng)坐標(biāo)
*Get,nnod,NODE,0,COUNT
*vget,nl,node,,nlist !得到表面節(jié)點(diǎn)編號(hào)
*vget,locx,node,,loc,x
…………………….
*DIM,locx1,array,nnod,1 !定義一個(gè)數(shù)組，其為nnod行1列
………………………….
要注意，這里面得到的nl是從小到大排列的，只包含一部分節(jié)點(diǎn)，而我們得到的locx卻是所有節(jié)點(diǎn)的坐標(biāo)，所以我們還需要定義一個(gè)locx1，再用一個(gè)循環(huán)把你想選擇的節(jié)點(diǎn)編號(hào)和其坐標(biāo)一一對(duì)應(yīng)起來(lái)。具體的關(guān)系從下面的圖可以看出。
*DO, j,1,nnod,1
locx1(j)=locx(nl(j)) !節(jié)點(diǎn)對(duì)應(yīng)坐標(biāo)
…………………………….
*ENDDO
這時(shí)我們就已經(jīng)得到了想選取的節(jié)點(diǎn)坐標(biāo)及對(duì)應(yīng)編號(hào)，此時(shí)我們需要運(yùn)行一個(gè)Output.mac文件，把得到的數(shù)組輸出。
Output.mac 中包含的內(nèi)容
!----------------------------------!
*cfopen,node_number.dat, ! Generate Ist File
*vwrite,nl(1)
(1F6.0)
*cfclos
*cfopen,node_locx.dat,
*vwrite,locx1(1)
(1E15.6)
*cfclos
………………….剩下的按照同樣格式寫(xiě)
!----------------------------------!
最后得到的txt文件的內(nèi)容分別如下：
展開(kāi) 在ANSYS中有些數(shù)據(jù)無(wú)法直接訪問(wèn)，需要通過(guò)定義單元表完成單元的結(jié)果的訪問(wèn)。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù)的詳細(xì)過(guò)程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數(shù)據(jù)，這里以Beam188的彎矩為例
2.1 啟動(dòng)ANSYS幫助菜單, 在索引框輸入Beam188然后搜索, 在單元輸出介紹找到彎矩的名稱(chēng)（代號(hào)）。
2.2 回到ANSYS界面，比如要輸出Mz, 則需要在添加SMISC,3 和SMISC,16 ,如圖
3. 輸出數(shù)據(jù)：Main Menu>General Postproc>Element Table> List E T, 選擇前面定義的SMISC,3 和SMISC,16 輸出單元I和J節(jié)點(diǎn)的Mz數(shù)值，如圖
4. 顯示彎矩云圖：Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res, 這里要注意要在LabI 選SMISC,3 LabJ 選SMSCI,16。
輸出彎矩到這就結(jié)束了，小編突然發(fā)現(xiàn)，輸出的彎矩值在每個(gè)單元的I和J處是一樣的（Beam188為2節(jié)點(diǎn)單元），彎矩圖也就成了鋸齒形，于是去問(wèn)了度娘一波，各路盆友給出解決方法，然而并沒(méi)有起作用的，于是乎我又想起來(lái)了“幫助文檔大法”，于是認(rèn)認(rèn)真真將Beam188的幫助文檔閱讀了一遍，功夫不負(fù)有心人，最終。。。
展開(kāi) 采用python語(yǔ)言提取rst 文件結(jié)果提取
問(wèn)題：
在有限元仿真中有時(shí)需要提取某些結(jié)構(gòu)的扭轉(zhuǎn)角度。Ansys workbench的結(jié)果后處理中可以設(shè)定圓柱坐標(biāo)系，然后按圓柱坐標(biāo)讀取Y軸的變形結(jié)果，再進(jìn)行扭轉(zhuǎn)角度的換算。
本文這里將該過(guò)程利用APDL命令進(jìn)行處理，避免一下步驟重復(fù)操作。
? 每次要單獨(dú)記錄變形量，
? 還要測(cè)量關(guān)鍵節(jié)點(diǎn)到坐標(biāo)系原點(diǎn)的距離，
? 將變形量和距離進(jìn)行角度換算（弧度）
? 弧度角轉(zhuǎn)角度
APDL后處理命令功能介紹：
1. 在坐標(biāo)系中創(chuàng)建所需的圓柱坐標(biāo)系，并在屬性ADPL name中進(jìn)行命名：aix （用戶隨意命名）
2. 在Named selection 定義需要查看的區(qū)域，并命名：load（用戶隨意命名）
3. 在后處理中插入command 命令，并將上述坐標(biāo)系和NS的名稱(chēng)修改。
4. 在command的結(jié)果屬性中就會(huì)有最大/最小/平均扭轉(zhuǎn)角度。并且為了方便校核準(zhǔn)確性還提供了沿圓柱坐標(biāo)系Y軸的變形量。
并且，除了界面顯示的結(jié)果外，還會(huì)在WB的結(jié)果文件夾中，顯示named Selection區(qū)域所有節(jié)點(diǎn)的編號(hào)/距離選定坐標(biāo)系的距離/沿坐標(biāo)系Y軸的變形量/換算后的角度值等信息，以便進(jìn)行其它數(shù)據(jù)處理。
展開(kāi) 在ANSYS 中3維坐標(biāo)下的 shell structure 使用2D 平面單元（僅考慮平面內(nèi)的位移）劃分，應(yīng)該使用哪個(gè)單元型號(hào)的單元？
ansys單元坐標(biāo)提取的相關(guān)專(zhuān)題、標(biāo)簽、搜索
ansys單元坐標(biāo)提取的最新內(nèi)容
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計(jì)算幾何特征尺寸
問(wèn)題：
在FKM關(guān)于結(jié)構(gòu)疲勞評(píng)估計(jì)算方法中指出：零部件特征尺寸，影響疲勞結(jié)果評(píng)估。原因是材料的應(yīng)力壽命曲線是由標(biāo)準(zhǔn)試樣進(jìn)行試驗(yàn)測(cè)試獲得的。當(dāng)零部件的特征尺寸與測(cè)試樣件不一致時(shí)，需要考慮零部件的特征尺寸這一因素。（一般而言，當(dāng)零部件的尺寸大于材料標(biāo)準(zhǔn)測(cè)試樣件時(shí)，零部件的表面或內(nèi)部缺陷發(fā)生的概率會(huì)增加
問(wèn)題：
在有限元仿真中有時(shí)需要提取某些結(jié)構(gòu)的扭轉(zhuǎn)角度。Ansys workbench的結(jié)果后處理中可以設(shè)定圓柱坐標(biāo)系，然后按圓柱坐標(biāo)讀取Y軸的變形結(jié)果，再進(jìn)行扭轉(zhuǎn)角度的換算。
本文這里將該過(guò)程利用APDL命令進(jìn)行處理，避免一下步驟重復(fù)操作。
? 每次要單獨(dú)記錄變形量，
? 還要測(cè)量關(guān)鍵節(jié)點(diǎn)到坐標(biāo)系原點(diǎn)的距離，
? 將變形量和距離進(jìn)行角度換算（弧度）
? 弧度角轉(zhuǎn)角度
采用python語(yǔ)言提取rst 文件結(jié)果提取
首先選取好你想選取的節(jié)點(diǎn)
NSEL,S,…………………..
然后使用*vget讀取節(jié)點(diǎn)編號(hào)及相應(yīng)坐標(biāo)
*Get,nnod,NODE,0,COUNT
*vget,nl,node,,nlist !得到表面節(jié)點(diǎn)編號(hào)
*vget,locx,node,,loc,x
…………………….
*DIM,locx1,array,nnod,1 !定義一個(gè)數(shù)組，其為nnod行1列
在ANSYS中有些數(shù)據(jù)無(wú)法直接訪問(wèn)，需要通過(guò)定義單元表完成單元的結(jié)果的訪問(wèn)。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù)的詳細(xì)過(guò)程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數(shù)據(jù)，這里以Beam188
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長(zhǎng)度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號(hào)為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號(hào)為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號(hào)為100的線的長(zhǎng)度
以上對(duì)應(yīng)的GUI操作： Utility Menu>Parameters
在ANSYS 中3維坐標(biāo)下的 shell structure 使用2D 平面單元（僅考慮平面內(nèi)的位移）劃分，應(yīng)該使用哪個(gè)單元型號(hào)的單元？
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib
1、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫(kù)文件
