abaqus柱狀坐標的案例
ABAQUS柱狀晶模型基于泰森多邊形建模
建立柱狀晶幾何模型進行有限元分析有助于深入理解材料的微觀結(jié)構(gòu)與宏觀性能之間的關(guān)系，為材料設(shè)計、制造工藝優(yōu)化及失效預(yù)測提供了強有力的工具。本案例介紹采用AutoCAD基于泰森多邊形算法生成柱狀晶三維幾何部件，并導(dǎo)入Abaqus有限元軟件內(nèi)建立包含晶粒及晶界在內(nèi)的柱狀晶模型。
采用CAD Voronoi V2.1插件建立泰森多邊形草圖，插件內(nèi)的孔壁厚參數(shù)即為需要建立的晶界厚度。
只保留紅色圖層內(nèi)容，對其他圖層進行隱藏，并通過面域REGION命令建立面域。
在三維建模中，通過拉伸并設(shè)置高度，將面域拉伸為三維實體模型。
通過并集，將所有晶粒合并為一個部件，并將整個模型在原位置復(fù)制一份。
手動建立一個長方體部件，尺寸與原模型相同，使長方體與原模型重合。
通過差集，將長方體與復(fù)制的模型進行差集，以建立晶界部件。建立完成后將整個模型導(dǎo)出為.iges格式文件備用。
在Abaqus內(nèi)將建立好的模型文件以部件的形式導(dǎo)入。
可對晶粒及晶界進行材料截面指派，也可對每個晶粒分別進行材料指派。
將模型單獨進行裝配，也可配合其他自定義部件裝配為整體。
對柱狀晶模型進行網(wǎng)格劃分。
展開 ABAQUS三維晶體結(jié)構(gòu)柱狀晶等軸晶建模
通過ABAQUS三維晶體塑性有限元建模，深入揭示柱狀晶微觀結(jié)構(gòu)（如晶粒尺寸、取向）與力學(xué)性能的關(guān)聯(lián)，為鑄造、焊接工藝優(yōu)化提供關(guān)鍵理論依據(jù)，顯著提升材料可靠性與使用壽命。本案例介紹在ABAQUS內(nèi)建立三維晶體結(jié)構(gòu)有限元模型。
柱狀晶體模型采用CAD Voronoi V2.1插件建模，首先建立二維Voronoi模型，并在CAD內(nèi)通過拉伸命令形成三維柱狀晶體。
等軸晶模型采用CAD Voronoi 3D插件直接生成。
CAD晶體結(jié)構(gòu)建模完成后，導(dǎo)出為iges文件，并以部件的形式導(dǎo)入到ABAQUS內(nèi)。
在ABAQUS內(nèi)建立裝配即可完成三維晶體結(jié)構(gòu)的建模。
展開 批量提取Abaqus的節(jié)點坐標(初始坐標、指定Step下的變形量、變形后節(jié)點坐標) ￥40
如果需要提取變形后的節(jié)點坐標的話，我們還要繼續(xù)進行。</p><h3>3.4提取Set組中初始坐標</h3><p>方法一：提取初始節(jié)點坐標可以通過Abaqus后處理查詢結(jié)果，在CSDN上可以查詢到，此處不再贅述。</p><p>方法二：通過Python腳本，按節(jié)點編號提取初始節(jié)點坐標。</p><p><span style="color: rgb(18, 18, 18);">代碼為附件：2 </span>Python腳本-初始節(jié)點提取</p><p><img src="https://img.jishulink.com/202405/attachment/90a3bc50ff224f1887815b8d0943ed64.jpg"></p><h2>4.數(shù)據(jù)處理</h2><p>通過數(shù)據(jù)處理，將對應(yīng)值求和即可求得變形后的節(jié)點坐標。以編號1號節(jié)點為例，節(jié)點初始坐標（X=19.42577，Y=2.96978），變形量（U1=0.54588，U2=0.48286），可求得編號1號節(jié)點變形后的坐標（X+U1，Y+U2），即為（19.97165，3.45264）。</p><p>驗證坐標如下圖所示，使用Hyperview后處理，可以看出提取節(jié)點坐標與Python腳本后處理的節(jié)點坐標一致。綜上所述可以看出該腳本可滿足需求。
展開 在ABAQUS中基于圓柱坐標系設(shè)置關(guān)于坐標函數(shù)的表面力（keyword 曲面加載,圓柱坐標,面力）
例如下圖所示，受Y方向某拉力作用，各點應(yīng)力狀態(tài)為：
在圓孔中心位置建立圓柱坐標系，該應(yīng)力狀態(tài)在圓柱坐標系下的公式為：
在這種情況下反推物理量，需要對曲面施加基于圓柱坐標系的面力。
案例如下：在圓弧面基于圓柱坐標系施加等效于單向應(yīng)力狀態(tài)的面力。
加載前先建立圓柱坐標系（注意R軸方向為0度位置，T軸方向為角度增大方向，示意圖見文后的加載圖）
具體設(shè)置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設(shè)置徑向力，按以下參數(shù)設(shè)置：
Distribution：應(yīng)力分配，點擊后面的f(x)創(chuàng)建一個基于圓柱坐標系的表達式，Local system 要選擇圓柱坐標系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點擊選擇圖標后，依次選擇（0，0，0） （-1，0，0） ，坐標選擇建立的圓柱坐標系。
注：面力方向矢量是基于所選坐標系，（-1，0，0）就是沿圓柱坐標系下的R軸反向。
Magnitude：選擇應(yīng)力大小為1。
然后在創(chuàng)建一個Load,設(shè)置切向力，如下圖所示，也是基于圓柱坐標系。
再創(chuàng)建一個Load，在整體坐標系下對兩側(cè)的平面施加Y方向的面力，大小為1，同時對后面的面施加全約束。
最后加載形式為下圖所示：
求解結(jié)果如下圖：
大部分位置應(yīng)力在0.99~1.01之間，為單向應(yīng)力狀態(tài)，加載方式正確。
本問題的關(guān)鍵是面力的方向問題，在選擇面力的方向矢量時，是基于所選坐標系。對于圓柱坐標系，切向力矢量為（0，-1，0）時，即力的方向只沿著theta的反方向。
展開 
ABAQUS 輸出節(jié)點坐標和積分點坐標
總結(jié)inp中添加關(guān)鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節(jié)點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結(jié)果文件是一樣的，都是節(jié)點坐標
展開 abaqus怎么提取極坐標一下的節(jié)點坐標？
如題
ABAQUS二次開發(fā)-參考點坐標自動識別，與更新坐標 ￥80
ABAQUS二次開發(fā)-參考點坐標自動識別，與更新坐標
ABAQUS輸出單元積分點坐標
方法
在ABAQUS CAE的場輸出中選擇的坐標點是節(jié)點的坐標，而節(jié)點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點的坐標也可以在節(jié)點場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場輸出時會附帶著一份XYZ坐標，這個應(yīng)該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數(shù)據(jù)精度的問題。
展開 hypermesh關(guān)聯(lián)abaqus-坐標系關(guān)聯(lián)邊界 ￥10
有時候使用hm去設(shè)置坐標系，都不太清楚邊界是否關(guān)聯(lián)上相應(yīng)的坐標系，只有打開abaqus查看才發(fā)現(xiàn)有點bug，重新校核下。
現(xiàn)下看下abaqus默認的
*Nset, nset=_T-PART-1-1-WW, internal
_M18,
_M19,
_M20,
_M21,
_M22,
_M23,
_M24,
*Transform, nset=_T-PART-1-1-WW
0.650207662680776, -0.749978756954605, 0.121498393026525, -0.152359230545774, 0.027956435087132, 0.987929705295229
嗯，abaqus是單獨將對應(yīng)的點作為一個Nset，然后調(diào)用*Transform關(guān)鍵字搞定的。
展開 ABAQUS處理前的節(jié)點坐標提取方法
abaqus在工作后，會自動生成.inp文件存在工作文件夾，通過如下操作找到.inp所在文件夾
2.找到.inp文件 ，建議用word打開，數(shù)據(jù)太多，記事本打不開，這里的數(shù)據(jù)點都是變形前的節(jié)點坐標。變形前后的節(jié)點編號相對應(yīng)
abaqus怎么獲取單元積分點上的坐標呢
如題 為什么我查詢的時候坐標都是0
ANSYS與Abaqus球坐標系下的結(jié)果讀取
ANSYS與Abaqus球坐標系下的結(jié)果讀取
1 概述
采用ANSYS和Abaqus軟件計算的結(jié)果通常默認的結(jié)果是在總體笛卡爾坐標系下產(chǎn)生的結(jié)果，這對于應(yīng)力或者應(yīng)變等分量的分析有時候不方便，比如對于一個圓筒體，比較關(guān)心其徑向應(yīng)力和環(huán)向應(yīng)力，而這個結(jié)果直接讀取使不可能的，需要一定的轉(zhuǎn)換。
這就是結(jié)果坐標系轉(zhuǎn)換。
在軟件里，應(yīng)力分量表示為sx,xy,xz（ANSYS），s11,s22,s22（Abaqus），當其轉(zhuǎn)換到柱坐標或者球坐標時，對應(yīng)的應(yīng)力分量就發(fā)生變化，sx和s11均表示徑向應(yīng)力。
2 ANSYS
建立一個球體模型，如圖1，加載求解，得到其總體坐標系下的sx應(yīng)力分量。
圖1
在后處理器中，將結(jié)果坐標系轉(zhuǎn)換為球坐標系，采用的命令為：RSYS。查詢ANSYS幫助文檔，如圖2：
圖2 RSYS
0,1,2分別代表笛卡爾坐標系，柱坐標系，球坐標系。
輸入命令：RSYS，2
顯式結(jié)果sx為圖3，此時的sx應(yīng)力分量為徑向應(yīng)力。
圖3
3 Abaqus
建立模型加載求解，得到s11應(yīng)力分量如圖4.
圖4
轉(zhuǎn)換結(jié)果坐標系，Visualization模塊下選擇 Tools--Create Coordinate Aystem，按指定方法建立局部坐標系，然后選擇Result-Option，選擇Transformation標簽，User-specified，就可以看到新建立的坐標系，選擇新建的坐標系即可完成坐標轉(zhuǎn)換。
如圖5，圖6
圖5
建立球坐標系的時候根據(jù)Abaqus窗口下方的提示進行操作。
圖6
最終轉(zhuǎn)換為徑向應(yīng)力的顯式結(jié)果，如圖7
圖7
展開 Abaqus-部件的坐標系變換(平移或者旋轉(zhuǎn)) ￥3.99
<div contenteditable="false" width="100%">
<p><a href="https://www.yqgqt.org.cn/major/abaqus" style="text-decoration: none;"></a><a href="https://www.yqgqt.org.cn/major/abaqus" style="text-decoration: none;"></a><a href="https://www.yqgqt.org.cn/major/abaqus" style="text-decoration: none;"></a><a href="https://www.yqgqt.org.cn/major/abaqus" style="text-decoration: none;"></a><a href="https://www.yqgqt.org.cn/major/abaqus" class="jsk-anchor">Abaqus</a>的【部件】的坐標系（平移、旋轉(zhuǎn)或二者的組合）變換，圖1為平移的示例。
展開 Abaqus 利用FindAt函數(shù)根據(jù)坐標查找點，線，面
--------------------------------------------------------------------------------
-----------------------------------查找線---------------------------------------
示例：
#加載
#一次選擇一條邊進行加載
a = mdb.models['Model-1'].rootAssembly
s1 = a.instances['Part-1-1'].edges
#這個點的坐標只需要在這條線上即可，這個坐標位置處不一定得有關(guān)鍵點存在
side1Edges1 =s1.findAt(((20.0,5.0,0.0),))
region = a.Surface(side1Edges=side1Edges1, name='Surf-1')
mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-1',
region=region, distributionType=UNIFORM, field='', magnitude=-pp,
amplitude=UNSET)
#####一次選擇兩條邊進行加載
#這個點的坐標只需要在這條線上即可，這個坐標位置處不一定得有關(guān)鍵點存在
side1Edges2 =s1.findAt(((10.0,10.0,0.0),),((-10.0,10.0,0),))
region2 = a.Surface(side1Edges=side1Edges2, name='Surf-1')
mdb.models['Model-1'].Pressure(name='Load-1', createStepName
展開 ABAQUS批量提取部件節(jié)點集節(jié)點編號及坐標腳本 ￥30
基于python的ABAQUS批量提取部件節(jié)點集節(jié)點編號及坐標二次開發(fā)腳本
