ABAQUS坐標(biāo)的案例
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
展開 批量提取Abaqus的節(jié)點(diǎn)坐標(biāo)(初始坐標(biāo)、指定Step下的變形量、變形后節(jié)點(diǎn)坐標(biāo)) ￥40
如果需要提取變形后的節(jié)點(diǎn)坐標(biāo)的話，我們還要繼續(xù)進(jìn)行。</p><h3>3.4提取Set組中初始坐標(biāo)</h3><p>方法一：提取初始節(jié)點(diǎn)坐標(biāo)可以通過(guò)Abaqus后處理查詢結(jié)果，在CSDN上可以查詢到，此處不再贅述。</p><p>方法二：通過(guò)Python腳本，按節(jié)點(diǎn)編號(hào)提取初始節(jié)點(diǎn)坐標(biāo)。</p><p><span style="color: rgb(18, 18, 18);">代碼為附件：2 </span>Python腳本-初始節(jié)點(diǎn)提取</p><p><img src="https://img.jishulink.com/202405/attachment/90a3bc50ff224f1887815b8d0943ed64.jpg"></p><h2>4.數(shù)據(jù)處理</h2><p>通過(guò)數(shù)據(jù)處理，將對(duì)應(yīng)值求和即可求得變形后的節(jié)點(diǎn)坐標(biāo)。以編號(hào)1號(hào)節(jié)點(diǎn)為例，節(jié)點(diǎn)初始坐標(biāo)（X=19.42577，Y=2.96978），變形量（U1=0.54588，U2=0.48286），可求得編號(hào)1號(hào)節(jié)點(diǎn)變形后的坐標(biāo)（X+U1，Y+U2），即為（19.97165，3.45264）。</p><p>驗(yàn)證坐標(biāo)如下圖所示，使用Hyperview后處理，可以看出提取節(jié)點(diǎn)坐標(biāo)與Python腳本后處理的節(jié)點(diǎn)坐標(biāo)一致。綜上所述可以看出該腳本可滿足需求。
展開 利用python找到abaqus模型中離坐標(biāo)點(diǎn)最近的線條邊——代碼示例 ￥5
利用python找到abaqus模型中離坐標(biāo)點(diǎn)最近的線條邊——代碼示例，親測(cè)代碼有效。
在ABAQUS中基于圓柱坐標(biāo)系設(shè)置關(guān)于坐標(biāo)函數(shù)的表面力（keyword 曲面加載,圓柱坐標(biāo),面力）
例如下圖所示，受Y方向某拉力作用，各點(diǎn)應(yīng)力狀態(tài)為：
在圓孔中心位置建立圓柱坐標(biāo)系，該應(yīng)力狀態(tài)在圓柱坐標(biāo)系下的公式為：
在這種情況下反推物理量，需要對(duì)曲面施加基于圓柱坐標(biāo)系的面力。
案例如下：在圓弧面基于圓柱坐標(biāo)系施加等效于單向應(yīng)力狀態(tài)的面力。
加載前先建立圓柱坐標(biāo)系（注意R軸方向?yàn)?度位置，T軸方向?yàn)榻嵌仍龃蠓较颍疽鈭D見文后的加載圖）
具體設(shè)置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設(shè)置徑向力，按以下參數(shù)設(shè)置：
Distribution：應(yīng)力分配，點(diǎn)擊后面的f(x)創(chuàng)建一個(gè)基于圓柱坐標(biāo)系的表達(dá)式，Local system 要選擇圓柱坐標(biāo)系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點(diǎn)擊選擇圖標(biāo)后，依次選擇（0，0，0） （-1，0，0） ，坐標(biāo)選擇建立的圓柱坐標(biāo)系。
注：面力方向矢量是基于所選坐標(biāo)系，（-1，0，0）就是沿圓柱坐標(biāo)系下的R軸反向。
Magnitude：選擇應(yīng)力大小為1。
然后在創(chuàng)建一個(gè)Load,設(shè)置切向力，如下圖所示，也是基于圓柱坐標(biāo)系。
再創(chuàng)建一個(gè)Load，在整體坐標(biāo)系下對(duì)兩側(cè)的平面施加Y方向的面力，大小為1，同時(shí)對(duì)后面的面施加全約束。
最后加載形式為下圖所示：
求解結(jié)果如下圖：
大部分位置應(yīng)力在0.99~1.01之間，為單向應(yīng)力狀態(tài)，加載方式正確。
本問(wèn)題的關(guān)鍵是面力的方向問(wèn)題，在選擇面力的方向矢量時(shí)，是基于所選坐標(biāo)系。對(duì)于圓柱坐標(biāo)系，切向力矢量為（0，-1，0）時(shí)，即力的方向只沿著theta的反方向。
展開 
ABAQUS二次開發(fā)-參考點(diǎn)坐標(biāo)自動(dòng)識(shí)別，與更新坐標(biāo) ￥80
ABAQUS二次開發(fā)-參考點(diǎn)坐標(biāo)自動(dòng)識(shí)別，與更新坐標(biāo)
ABAQUS 輸出節(jié)點(diǎn)坐標(biāo)和積分點(diǎn)坐標(biāo)
總結(jié)inp中添加關(guān)鍵字
輸出單元的積分點(diǎn)坐標(biāo)：*EL FILE
COORD
輸出節(jié)點(diǎn)坐標(biāo)：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個(gè)COORD生成的結(jié)果文件是一樣的，都是節(jié)點(diǎn)坐標(biāo)
展開 abaqus怎么提取極坐標(biāo)一下的節(jié)點(diǎn)坐標(biāo)？
如題
hypermesh關(guān)聯(lián)abaqus-坐標(biāo)系關(guān)聯(lián)邊界 ￥10
有時(shí)候使用hm去設(shè)置坐標(biāo)系，都不太清楚邊界是否關(guān)聯(lián)上相應(yīng)的坐標(biāo)系，只有打開abaqus查看才發(fā)現(xiàn)有點(diǎn)bug，重新校核下。
現(xiàn)下看下abaqus默認(rèn)的
*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是單獨(dú)將對(duì)應(yīng)的點(diǎn)作為一個(gè)Nset，然后調(diào)用*Transform關(guān)鍵字搞定的。
展開 [軟件使用]abaqus殼單元局部坐標(biāo)系，你學(xué)會(huì)了嗎？
在使用abaqus進(jìn)行有限元分析的工作中，確定殼單元局部坐標(biāo)系是一項(xiàng)重要的工作，其原因之一在于在abaqus中，殼單元的位移輸出基于整體坐標(biāo)系，應(yīng)力應(yīng)變輸出基于局部坐標(biāo)系，因此如果不能準(zhǔn)確地確定殼單元的局部坐標(biāo)系，在后處理查看計(jì)算結(jié)果時(shí)可能會(huì)無(wú)法準(zhǔn)確理解計(jì)算結(jié)果。
通常情況下，殼單元的局部坐標(biāo)系如下圖所示，其包含平面內(nèi)的1，2軸和平面法線的n軸（3軸）。顯然，n軸由殼單元所在平面確定，但是其有兩種選擇，即由“殼內(nèi)指向殼外”和由“殼外指向殼內(nèi)”。
那么在abaqus中，殼單元的局部坐標(biāo)系依據(jù)以下規(guī)則定義：
（1）對(duì)于一個(gè)3節(jié)點(diǎn)/4節(jié)點(diǎn)殼單元，按照右手定則，拇指指向即為n軸方向。
殼單元節(jié)點(diǎn)順序?yàn)?-2-4-3時(shí)的n軸方向。
（2）確定好n軸之后，接下來(lái)的1軸和2軸按照以下規(guī)則確定：
將整體坐標(biāo)系的X軸投影到殼單元上，投影方向即為1軸。再按照右手定則，1-2-n軸形成右手坐標(biāo)系，即右手拇指指向n軸時(shí)，其余4指的旋轉(zhuǎn)方向從1軸轉(zhuǎn)向2軸，具體圖解如下：右側(cè)為整體坐標(biāo)系，左手為局部坐標(biāo)系。
按照上述規(guī)則必然會(huì)存在一種特殊情況，即整體1軸與殼單元垂直，則此時(shí)整體1軸投影到殼單元上會(huì)是一個(gè)點(diǎn)，無(wú)法確定局部1軸方向，在這種情況下，abaqus采用整體3軸投影到殼單元上作為局部1軸方向。
以上就是殼單元局部坐標(biāo)系的確定過(guò)程，下面以一個(gè)例子，來(lái)表明殼單元局部坐標(biāo)系確定的具體作用。
以如圖所示外壓圓環(huán)為例：
計(jì)算完成后，后處理S11應(yīng)力分布如下：
S22分布：
很明顯，應(yīng)力云圖不符合常規(guī)理解。均勻外壓圓環(huán)的應(yīng)力分布應(yīng)當(dāng)是相對(duì)均勻的，而不會(huì)出現(xiàn)在“某一格”的單元應(yīng)力分布明顯不同于其他單元。
展開 ABAQUS輸出單元積分點(diǎn)坐標(biāo)
方法
在ABAQUS CAE的場(chǎng)輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來(lái)的，單元積分點(diǎn)的信息相對(duì)真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無(wú)法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場(chǎng)輸出下面添加COORD（這和CAE中場(chǎng)輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場(chǎng)輸出時(shí)會(huì)附帶著一份XYZ坐標(biāo)，這個(gè)應(yīng)該也可以當(dāng)做單元的坐標(biāo)，，但是我比較過(guò)這個(gè)附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo)，有時(shí)候有點(diǎn)差別，可能是數(shù)據(jù)精度的問(wèn)題。
展開 abaqus二次開發(fā)：批量更新硬點(diǎn)（RP點(diǎn)）坐標(biāo)及名稱 ￥198
abaqus二次開發(fā)：批量更新硬點(diǎn)（RP點(diǎn)）坐標(biāo)及名稱：
實(shí)現(xiàn)的功能：
1、自動(dòng)識(shí)別<a href="/major/<a href="/major/abaqus">abaqus/CAE模型中存在的所有硬點(diǎn)名稱、坐標(biāo)XYZ，硬點(diǎn)包括“子從點(diǎn)”或”父點(diǎn)”
2、父點(diǎn)坐標(biāo)可重新自定義修改坐標(biāo)XYZ，及新命名，也可以選擇點(diǎn)擊copy是否創(chuàng)建復(fù)制的子從點(diǎn)
3、子從點(diǎn)不可更改坐標(biāo)，僅可以重命名，也可以選擇點(diǎn)擊copy是否創(chuàng)建復(fù)制的子從點(diǎn)
4、附件：加密文件
5:、需要源碼文件的，可私聊961765699
展開 
abaqus怎么獲取單元積分點(diǎn)上的坐標(biāo)呢
如題 為什么我查詢的時(shí)候坐標(biāo)都是0
ABAQUS處理前的節(jié)點(diǎn)坐標(biāo)提取方法
abaqus在工作后，會(huì)自動(dòng)生成.inp文件存在工作文件夾，通過(guò)如下操作找到.inp所在文件夾
2.找到.inp文件 ，建議用word打開，數(shù)據(jù)太多，記事本打不開，這里的數(shù)據(jù)點(diǎn)都是變形前的節(jié)點(diǎn)坐標(biāo)。變形前后的節(jié)點(diǎn)編號(hào)相對(duì)應(yīng)
abaqus二次開發(fā)：前處理批量更改RP點(diǎn)坐標(biāo)及名稱（源碼+注釋） ￥168
abaqus二次開發(fā)：前處理批量更改RP點(diǎn)坐標(biāo)及名稱（源碼+注釋）
購(gòu)買后，請(qǐng)加V：wzd_1021_
更多插件內(nèi)容，請(qǐng)私聊
同時(shí)提供定制hypermesh/hyperview/abaqus前后處理插件開發(fā)。
Abaqus 利用FindAt函數(shù)根據(jù)坐標(biāo)查找點(diǎn)，線，面
--------------------------------------------------------------------------------
-----------------------------------查找線---------------------------------------
示例：
#加載
#一次選擇一條邊進(jìn)行加載
a = mdb.models['Model-1'].rootAssembly
s1 = a.instances['Part-1-1'].edges
#這個(gè)點(diǎn)的坐標(biāo)只需要在這條線上即可，這個(gè)坐標(biāo)位置處不一定得有關(guān)鍵點(diǎ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)
#####一次選擇兩條邊進(jìn)行加載
#這個(gè)點(diǎn)的坐標(biāo)只需要在這條線上即可，這個(gè)坐標(biāo)位置處不一定得有關(guān)鍵點(diǎ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坐標(biāo)的相關(guān)專題、標(biāo)簽、搜索
ABAQUS坐標(biāo)導(dǎo)入abaqus 坐標(biāo)abaqus 導(dǎo)出坐標(biāo)abaqus查看坐標(biāo)abaqus創(chuàng)建坐標(biāo)abaqus坐標(biāo)網(wǎng)格 Abaqus abaqus積分點(diǎn)坐標(biāo)abaqus子程序節(jié)點(diǎn)坐標(biāo)abaqus讀取節(jié)點(diǎn)坐標(biāo)并與實(shí)際坐標(biāo)abaqus讀取節(jié)點(diǎn)坐標(biāo)并與實(shí)際坐標(biāo)矯正abaqus將坐標(biāo)轉(zhuǎn)換為柱坐標(biāo)abaqus直角坐標(biāo)系與柱坐標(biāo)轉(zhuǎn)換abaqus提取極坐標(biāo)下的節(jié)點(diǎn)坐標(biāo)
