abaqus球面坐標的案例
批量提取Abaqus的節(jié)點坐標(初始坐標、指定Step下的變形量、變形后節(jié)點坐標) ￥40
如果需要提取變形后的節(jié)點坐標的話，我們還要繼續(xù)進行。</p><h3>3.4提取Set組中初始坐標</h3><p>方法一：提取初始節(jié)點坐標可以通過Abaqus后處理查詢結果，在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ù)處理，將對應值求和即可求得變形后的節(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ù)的表面力（keyword 曲面加載,圓柱坐標,面力）
例如下圖所示，受Y方向某拉力作用，各點應力狀態(tài)為：
在圓孔中心位置建立圓柱坐標系，該應力狀態(tài)在圓柱坐標系下的公式為：
在這種情況下反推物理量，需要對曲面施加基于圓柱坐標系的面力。
案例如下：在圓弧面基于圓柱坐標系施加等效于單向應力狀態(tài)的面力。
加載前先建立圓柱坐標系（注意R軸方向為0度位置，T軸方向為角度增大方向，示意圖見文后的加載圖）
具體設置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設置徑向力，按以下參數(shù)設置：
Distribution：應力分配，點擊后面的f(x)創(chuàng)建一個基于圓柱坐標系的表達式，Local system 要選擇圓柱坐標系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點擊選擇圖標后，依次選擇（0，0，0） （-1，0，0） ，坐標選擇建立的圓柱坐標系。
注：面力方向矢量是基于所選坐標系，（-1，0，0）就是沿圓柱坐標系下的R軸反向。
Magnitude：選擇應力大小為1。
然后在創(chuàng)建一個Load,設置切向力，如下圖所示，也是基于圓柱坐標系。
再創(chuàng)建一個Load，在整體坐標系下對兩側的平面施加Y方向的面力，大小為1，同時對后面的面施加全約束。
最后加載形式為下圖所示：
求解結果如下圖：
大部分位置應力在0.99~1.01之間，為單向應力狀態(tài)，加載方式正確。
本問題的關鍵是面力的方向問題，在選擇面力的方向矢量時，是基于所選坐標系。對于圓柱坐標系，切向力矢量為（0，-1，0）時，即力的方向只沿著theta的反方向。
展開 ABAQUS 輸出節(jié)點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節(jié)點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節(jié)點坐標
展開 abaqus怎么提取極坐標一下的節(jié)點坐標？
如題
ABAQUS二次開發(fā)-參考點坐標自動識別，與更新坐標 ￥80
ABAQUS二次開發(fā)-參考點坐標自動識別，與更新坐標
hypermesh關聯(lián)abaqus-坐標系關聯(lián)邊界 ￥10
有時候使用hm去設置坐標系，都不太清楚邊界是否關聯(lián)上相應的坐標系，只有打開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是單獨將對應的點作為一個Nset，然后調用*Transform關鍵字搞定的。
展開 ABAQUS輸出單元積分點坐標
方法
在ABAQUS CAE的場輸出中選擇的坐標點是節(jié)點的坐標，而節(jié)點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點的坐標也可以在節(jié)點場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結果中創(chuàng)建場輸出時會附帶著一份XYZ坐標，這個應該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數(shù)據(jù)精度的問題。
展開 ABAQUS處理前的節(jié)點坐標提取方法
abaqus在工作后，會自動生成.inp文件存在工作文件夾，通過如下操作找到.inp所在文件夾
2.找到.inp文件 ，建議用word打開，數(shù)據(jù)太多，記事本打不開，這里的數(shù)據(jù)點都是變形前的節(jié)點坐標。變形前后的節(jié)點編號相對應
ANSYS與Abaqus球坐標系下的結果讀取
ANSYS與Abaqus球坐標系下的結果讀取
1 概述
采用ANSYS和Abaqus軟件計算的結果通常默認的結果是在總體笛卡爾坐標系下產(chǎn)生的結果，這對于應力或者應變等分量的分析有時候不方便，比如對于一個圓筒體，比較關心其徑向應力和環(huán)向應力，而這個結果直接讀取使不可能的，需要一定的轉換。
這就是結果坐標系轉換。
在軟件里，應力分量表示為sx,xy,xz（ANSYS），s11,s22,s22（Abaqus），當其轉換到柱坐標或者球坐標時，對應的應力分量就發(fā)生變化，sx和s11均表示徑向應力。
2 ANSYS
建立一個球體模型，如圖1，加載求解，得到其總體坐標系下的sx應力分量。
圖1
在后處理器中，將結果坐標系轉換為球坐標系，采用的命令為：RSYS。查詢ANSYS幫助文檔，如圖2：
圖2 RSYS
0,1,2分別代表笛卡爾坐標系，柱坐標系，球坐標系。
輸入命令：RSYS，2
顯式結果sx為圖3，此時的sx應力分量為徑向應力。
圖3
3 Abaqus
建立模型加載求解，得到s11應力分量如圖4.
圖4
轉換結果坐標系，Visualization模塊下選擇 Tools--Create Coordinate Aystem，按指定方法建立局部坐標系，然后選擇Result-Option，選擇Transformation標簽，User-specified，就可以看到新建立的坐標系，選擇新建的坐標系即可完成坐標轉換。
如圖5，圖6
圖5
建立球坐標系的時候根據(jù)Abaqus窗口下方的提示進行操作。
圖6
最終轉換為徑向應力的顯式結果，如圖7
圖7
展開 abaqus怎么獲取單元積分點上的坐標呢
如題 為什么我查詢的時候坐標都是0
Abaqus 利用FindAt函數(shù)根據(jù)坐標查找點，線，面
--------------------------------------------------------------------------------
-----------------------------------查找線---------------------------------------
示例：
#加載
#一次選擇一條邊進行加載
a = mdb.models['Model-1'].rootAssembly
s1 = a.instances['Part-1-1'].edges
#這個點的坐標只需要在這條線上即可，這個坐標位置處不一定得有關鍵點存在
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)
#####一次選擇兩條邊進行加載
#這個點的坐標只需要在這條線上即可，這個坐標位置處不一定得有關鍵點存在
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ā)腳本
Abaqus-部件的坐標系變換(平移或者旋轉) ￥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>的【部件】的坐標系（平移、旋轉或二者的組合）變換，圖1為平移的示例。
展開 [軟件使用]abaqus殼單元局部坐標系，你學會了嗎？
在使用abaqus進行有限元分析的工作中，確定殼單元局部坐標系是一項重要的工作，其原因之一在于在abaqus中，殼單元的位移輸出基于整體坐標系，應力應變輸出基于局部坐標系，因此如果不能準確地確定殼單元的局部坐標系，在后處理查看計算結果時可能會無法準確理解計算結果。
通常情況下，殼單元的局部坐標系如下圖所示，其包含平面內的1，2軸和平面法線的n軸（3軸）。顯然，n軸由殼單元所在平面確定，但是其有兩種選擇，即由“殼內指向殼外”和由“殼外指向殼內”。
那么在abaqus中，殼單元的局部坐標系依據(jù)以下規(guī)則定義：
（1）對于一個3節(jié)點/4節(jié)點殼單元，按照右手定則，拇指指向即為n軸方向。
殼單元節(jié)點順序為1-2-4-3時的n軸方向。
（2）確定好n軸之后，接下來的1軸和2軸按照以下規(guī)則確定：
將整體坐標系的X軸投影到殼單元上，投影方向即為1軸。再按照右手定則，1-2-n軸形成右手坐標系，即右手拇指指向n軸時，其余4指的旋轉方向從1軸轉向2軸，具體圖解如下：右側為整體坐標系，左手為局部坐標系。
按照上述規(guī)則必然會存在一種特殊情況，即整體1軸與殼單元垂直，則此時整體1軸投影到殼單元上會是一個點，無法確定局部1軸方向，在這種情況下，abaqus采用整體3軸投影到殼單元上作為局部1軸方向。
以上就是殼單元局部坐標系的確定過程，下面以一個例子，來表明殼單元局部坐標系確定的具體作用。
以如圖所示外壓圓環(huán)為例：
計算完成后，后處理S11應力分布如下：
S22分布：
很明顯，應力云圖不符合常規(guī)理解。均勻外壓圓環(huán)的應力分布應當是相對均勻的，而不會出現(xiàn)在“某一格”的單元應力分布明顯不同于其他單元。
展開 Abaqus疑難雜癥——局部坐標系的那些事兒
本篇文章將詳細講解Abaqus/CAE中局部坐標系的一些故事，內容來源于本人平時學習軟件時的心得和官方在線手冊以及曹金鳳老師、石亦平博士編寫的《ABAQUS有限元分析常見問題解答》，分為基礎小白篇（面向初學者）和高手進階篇（面向中級Abaqus仿真師）。
01
基礎小白篇（界面“點點點”）
1.在Abaqus/CAE中建立部件，然后點擊Tools——Datum;
2.點擊CSYS——Three points(以最常用的三點建立坐標系為例)
3.直角坐標系：先指明原點，按Enter鍵確定，然后在圖中選中自定義的X軸點的坐標（或自己直接輸入X軸點的坐標）按Enter鍵確定，最后在圖中選中自定義的X軸點的坐標（或自己直接輸入X軸點的坐標）按Enter鍵確定。
4.柱坐標系以及球坐標系：先指明原點，按Enter鍵確定，然后在圖中選中徑向坐標軸R上的點（或者自己輸入徑向坐標軸R上的點），最后在圖中選中周向坐標軸Theta上的點（或者自己輸入周向坐標軸Theta上的點）。
展開 