ansys 節點應力的案例
ansys中的節點應力
我想知道ansys中的節點應力是如何得到的？因為理論上講應力應該是針對微元體來講的，單純的節點是不存在應力的，那么ansys中結果所提供的節點應力是怎樣得到的？與單元表所顯示的應力往往存在較大差別，那實際進行強度分析的時候應該以哪個為準呢？
【原創】能夠生成ANSYS中節點應力釋放所需文件的小程序
有節點限制，請注意！
根據上一步的計算結果（nodecal）數據文件生成ansys應力釋放所需要的節點應力文件，
可以按不同比例生成應力文件.歡迎大家給出意見。
使用方法：將結果數據文件，命名為exam.dat，具體格式如例子。
運行node_force.exe，即可。
生成的nodeforce.dat就是ansys所需格式的文件，用input讀入即可。
可以大大提高對各節點進行應力釋放的效率！為平面應變的隧道開挖而設計！
New Folder.rar
ANSYS如何提取某一節點的應力時程 ￥100
那么如何提取某一個節點的von Mises stress呢？
首先明確ANSYS的節點附加在單元上，可以通過選擇單元上節點的方法提取節點應力。
1 確定節點所在單元，顯示節點編號。
例單元號8560，節點號8678。
2 進入TimeHist Postpro, 定義變量。
3變量顯示。
付費內容為相關命令流。
通過Abaqus python腳本批量獲取節點的應力 ￥25
背景
有限單元法計算單元積分點的應力應變值，而對于節點的應力應變值是通過外插得到的，Abaqus中云圖顯示的就是經過插值和平均后的節點的值。通過工具欄的Query-Probe values可以查看單元或節點的應力應變等結果。
對于自動化的后處理場景，通常需要自動批量地獲取單元/節點的結果，通常都需要通過python腳本來實現。通過類似odb.steps['Step-1'].frames[-1].fieldOutputs['S']的場輸出可以比較方便地直接獲得單元的積分點應力，但沒有直接的API可以獲取節點的應力應變等結果。
如果需要獲取部件表面節點應力，可以通過創建路徑+XYData的方式實現，但想要獲得最大節點應力，則該方式不便實現。
2. 通過python腳本獲取節點應力結果
本文通過fieldOutput.getSub()函數獲取所有單元的節點結果，并對每一節點關聯的多個單元的節點值進行平均后得到節點的結果。以下以某個簡單的odb結果進行驗證。
（1）批量獲得節點的mises應力值
（2）批量獲得節點的X方向正應力值
（3）批量獲得節點的最大主應力值
（4）獲取節點的最大mises應力及編號
3. 獲取節點應變等結果
只需將腳本程序中的應力場改為應變成E等即可，此處不再演示。
以下為本文的python腳本代碼（代碼中作了必要的簡單注釋）。
展開 
基于節點位移的應力強度因子外推法
今天木木給大家分享的是基于節點位移求解應力強度因子，相比于上一期出的基于單元應力求解應力強度因子得出的結果更加接近解析解。這一期包括以下內容：（1）簡要講述INP文件（2）運用最小二乘法進行線性擬合（3）對裂尖數據進行特殊處理。
ANSYS中單元解、節點解以及節點單元解的概念解析
理論上，任何結構任何位置處的應力應變應該都是連續的，而上面所說的單元應力應變解并不連續，因而就出現了另外一個解，我個人稱之為節點單元解，它是單元解在公共節點上應力應變值的平均值，通過平均化就使得公共節點上的應力應變值變得唯一，但這樣會帶來另外一個問題，就是節點單元解和節點有關，也即是和單元數目有關。在某些情況下，可能會由于網格劃分的影響，導致畸變較大。
總結起來，三個解的概念如下：
節點解：節點位移解，原始解，最為精確的解；
單元解：單元的應力應變，派生解，通過節點解推導得到；
節點單元解：節點的應力應變，派生解的平均化顯示。
祝好
ANSYS結構院
2017.12.25
展開 ANSYS中單元解、節點解以及節點單元解該怎么理解
總結起來，三個解的概念如下：
節點解：節點位移解，原始解，最為精確的解；
單元解：單元的應力應變，派生解，通過節點解推導得到；
節點單元解：節點的應力應變，派生解的平均化顯示。
來源：ANSYS學習與應用
單元積分點應力如何外插至節點上 | 數值實現篇
如果我們還想看一下細節方面的，以1號單元的節點應力s11為例：
自研程序與Abaqus的結果也是一致的，在提取Abaqus單元節點應力時，應該將應力平滑選項取消勾選，即：
單元積分點應力外插matlab函數
function [StressElem,StressNode] = QuadNodeStress(node, element, prop, U, averageType,elemType,guassType)
% 通過節點位移計算節點應力，正應力：Sxx、Syy、Sxy、VonMises
% 增加節點應力均勻化標識：averageType，==1時，采用繞節點直接平均，==2時采用繞節點面積加權平均
E = prop(1);
NU = prop(2);
ID = prop(4);
[numberNodes, ~] = size(node);
[numberElements, ~] = size(element);
StressElem = zeros(numberElements, 3); % 只計算出正應力Sxx、Syy、Sxy即可
StressNode = zeros(numberNodes, 4);
WeightSum = zeros(numberNodes, 1); % 用于加權平均的權重總和
% 根據平面應力/應變狀態ID選擇應力-應變矩陣
if ID == 1
D = (E/(1-NU^2)) * [1, NU, 0; NU, 1, 0; 0, 0, (1-NU)/2];
elseif ID == 2
D = (E/(1+NU)/(1-2*NU)) * [1-NU, NU, 0;
展開 ABAQUS Fortran基于提供的代碼 P53，完善 8 節點單元程序，并增加對應力結果的處理。
1.改寫輸入數據格式，使之能適應任意幾何（可利用節點坐標輸 入節點，利用單元-節點關系輸入單元）； 2. 計算節點應力，給出并實現至少一種應力處理方案，提供處理 前后的應力結果（可用表格和云圖表示），可與其它軟件對比； 3.提交總結報告（包括方法/方案描述、帶詳細注釋的代碼、程序框圖、算例描述、結果比較分析等）、可編譯源代碼、可執行文件、 數據文件、結果文件
program p53
!-----------------------------------------------------------------------------
! program 5.3 plane strain of an elastic solid using uniform
! 8-node quadrilateral elements numbered in the x direction
!-----------------------------------------------------------------------------
use new_library ; use geometry_lib ; implicit none
integer::nels,nxe,neq,nband,nn,nr,nip,nodof=2,nod=8,nst=3,ndof,loaded_nodes,&
i,k,iel,ndim=2
real::aa,bb,e,v,det ; character(len=15) :: element = 'quadrilateral'
!
展開 LMS Virtual.Lab Motion_教程53之如何尋找最大/最小應力/應變節點
之前發過幾個關于剛柔耦合的帖子，今天介紹一個剛柔耦合后處理的使用方法，它的功能是用來找出應力/應變/位移的最大值/最小值，同時還能找出出現該值的節點的編號等信息，可能會對有些人有幫助。
使用的模型還是大家都熟知的雷達的模型，本帖以尋找最大應力值以及出現的節點為例。
剛柔耦合分析完成之后，Export Modal Participation factors。
新窗口中打開柔性化的部件，插入一個Linear Superposition Feature。
在Linear Superposition Set中選擇Load-FE Assignment。
創建一個Stress Image。
再創建一個Time-Series Extrema，并且選擇剛才創建的Stress Image。
計算剛剛創建的Time-Series Extrema。
Report Time-Series Extrema，確定最大值出現的時間點，在Stress Image中選擇該時間點。
顯示Image Extrema，就可以指出最大應力點的位置了。
更多下載資料請關注百度網盤LMS_VL_Motion，Moiton交流群：324201728
展開 ABAQUS-如何求結構的節點位移單元應力分量和支反力
這兩個輸出的是節點位移與支反力。
圖15
圖16
在位置
下拉框中選擇積分點，在新窗口選擇S下拉菜單中的S11。這輸出的是單元應力分量。
圖17
點擊窗口的設置按鈕，更改名稱為link.rpt。
圖18
在保存目錄中用記事本打開
link
.
rpt
，即可得到相應的結果。
圖19
微信公眾號：CAE小花生
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
注：由于技術鄰排版風格有限，故部分內容顯示不全，感興趣的小伙伴可點擊原文進行閱覽：
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是：有限元計算過程中，單元積分點應力如何外推至節點？
有關積分點與節點的概念可點擊跳轉閱讀歷史推文：有限元基本概念-【節點和積分點】，現科普一下Q4單元、Q8單元、Q9單元的形函數和高斯積分方案。
Q4單元
Q8/9單元
應力外插
核心理念：坐標系的轉換。
假設是母單元的自然坐標系，是由高斯積分點控制的坐標系（術語可能不專業），假設高斯積分方案為。坐標系轉換關系：
單元內任一點的應力，由4個高斯積分點應力進行插值時，可表示為
其中，是基于高斯積分點的形函數，第一個積分點的坐標在母單元坐標系下為（-1，-1），根據上述的坐標系轉換的方式，在高斯積分點的坐標系下，第一個單元節點在高斯積分點坐標系下坐標為，將此坐標值代入第一個形函數，得，相同的道理，可推導至四個節點在4個形函數下的外插矩陣：
對于Q8、Q9單元，依然可采用高斯積分方案（減縮積分）。
展開 有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 ansys導入節點坐標數據 附80多種ANSYS常用材料的參數文件下載
有時候，再用ansys做一些復雜的模型分析時候（如：桁架，拱形架，繩網等），因為其模型數量很多，模型空間位置相對復雜，采用apdl語言實現可能比較繁瑣或者會遇到調試方面的不便。所以，我們可以用數據處理功能更為強大的matlab或者c++進行編程，將節點坐標直接導入到ansys中進行分析。
matlab可用如下格式導出節點坐標：
接下來，采用apdl語言定義存放數據的數組：（如下圖）注意：（3F5.2要和matlab的fprintf中%5.2f對應）
將存放數組的.txt文件與坐標.txt放在工作目錄下：
在菜單中選擇file——read to file——選擇“wang.txt”，程序自動搜索到存放在nn.txt的坐標數據。
接下來，我們就可以在數組文件中看到導入的數據了：
下載地址：80多種ANSYS常用材料的參數文件
展開 ANSYS使用APDL語言提取節點編號及對應坐標 ￥10
首先選取好你想選取的節點
NSEL,S,…………………..
然后使用*vget讀取節點編號及相應坐標
*Get,nnod,NODE,0,COUNT
*vget,nl,node,,nlist !得到表面節點編號
*vget,locx,node,,loc,x
…………………….
*DIM,locx1,array,nnod,1 !定義一個數組，其為nnod行1列
………………………….
要注意，這里面得到的nl是從小到大排列的，只包含一部分節點，而我們得到的locx卻是所有節點的坐標，所以我們還需要定義一個locx1，再用一個循環把你想選擇的節點編號和其坐標一一對應起來。具體的關系從下面的圖可以看出。
*DO, j,1,nnod,1
locx1(j)=locx(nl(j)) !節點對應坐標
…………………………….
*ENDDO
這時我們就已經得到了想選取的節點坐標及對應編號，此時我們需要運行一個Output.mac文件，把得到的數組輸出。
Output.mac 中包含的內容
!----------------------------------!
*cfopen,node_number.dat, ! Generate Ist File
*vwrite,nl(1)
(1F6.0)
*cfclos
*cfopen,node_locx.dat,
*vwrite,locx1(1)
(1E15.6)
*cfclos
………………….剩下的按照同樣格式寫
!----------------------------------!
最后得到的txt文件的內容分別如下：
展開 