
發布
注冊
/
登錄ansys中積分點
關注創建者:王靖雯 創建時間:2023-03-07
ansys中積分點的視頻教程
有限元中單元積分點與節點應力相互轉換(二維和三維)
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取. 但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
免費 2小時16分鐘 3603播放
查看
ansys中積分點的實例教程
節點和積分點是有限單元法(FEM)的兩個基本概念,初涉有限元計算的同志往往在這點上產生混淆,假設導師面試的時候,問單元應力是什么,若回答不慎,將貽笑大方,得不償失。本文試圖以簡略易懂的說法來闡述節點和積分點的區別。
1.節點位移是有限元法的基本未知量。節點構筑了問題域的幾何離散化形狀,節點是形函數的零點,通常形函數是以節點為依據進行假設的。形函數決定了單元內部各點運動的位移模式(常用帕斯卡三角形來選擇單元位移模式),這樣就形成了數學上所說的插值。
有限元法的原理就是將問題域分割成N多小單元,在每個單元內采用簡單的函數來近似表達單元的真實位移,將各單元再連接起來,就可以近似描述整個問題域的運動。因此,有限元法從根本上就是精確的,而不是準確的。
2.積分點是單元進行數值積分的已知量。有限元法中一般采用高斯積分,但是積分方法不限于高斯積分,如果有人用了Irons積分或者Hammer積分,請不要驚訝。在形成單元剛度矩陣和進行節點應力磨平的時候,需要高斯積分。
以等參單元為例,其剛度矩陣
,這個就需要數值積分來快速計算,高斯點坐標及權系數如表4.2[王勖成]所示。
老師授課時一般對常應力單元進行推導,而常應力單元只有一個積分點,被積函數是常數,因此體現不出高斯積分來。很多老師對高斯積分在單元剛度矩陣的應用不予細述,導致部分同學對單元積分點認識不足。
3.單元應力指的是高斯積分點的應力,而非節點上的應力。有了位移模式,再通過虛功原理得到單元剛度矩陣,然后聚合總剛,求解平衡方程,就會把基本未知量——節點位移求出來了。通過節點位移得到單元應變結果,利用物理方程求得單元應力結果。
在等參元中,單元中n+1階(n=p-m)高斯積分點上的應變或應力近似解比其它部位具有較高的精度,因此我們稱(n+1)階高斯積分點是等參元中的最佳應力點。
展開 可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為:
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為:
展開 經常有小伙伴問獲取積分點坐標的方法,今天給大家介紹三種獲取積分點坐標的方式,希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標,例如我們熟知的UMAT子程序中包含COORD參數,即為積分點坐標。順帶一提的是,當打開了幾何非線性時,該積分點是當前構形下的坐標,如果未打開幾何非線性則為初始坐標。
2通過history output輸出積分點坐標
Abaqus可以直接在歷程變量history output中輸出積分點坐標。直接在history output中勾選COORD選項,但是這里需要注意的是,Domain中的Set集合如果是node set,這里輸出來的是節點坐標,當這里是element set的時候,輸出來的才是積分點坐標。
3通過等參單元映射函數計算
等參元中,為了方便計算,把整體坐標映射到自然坐標,然后在自然坐標下進行高斯積分。如果知道了自然坐標下的高斯積分點,通過映射函數反算,便能得到整體坐標下的高斯積分點坐標。以四邊形等參單元為例,其以自然呢坐標為變量的插值形函數如下
坐標變換采取同樣的插值函數(叫做等參的原因),整體坐標和自然坐標的關系式如下,如果知道自然坐標下的高斯積分點,直接通過此公式計算其在整體坐標下的坐標。
展開 按照正常的理解,毫無.疑問,abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看,似乎也是如此,每個單元存在8個積分點。
然而,如果自己動手跑一遍程序,就會發現事實遠非如此,采用全積分計算得到的結果與abaqus 存在差異,原因何在?
事實賞,abaqus C3D8 采用的選擇積分方式(selective intergation schema),即對于偏應變,采用8個積分,對于球應變,采用中心點積分。這樣計算得到的結果才能與abaqus 完全對標,亦可從abaqus 幫助文檔得到答案。
展開 在下半節課中,詳細地分析了扭曲單元與有限元精度之間的關系。我們常聽到單元網格質量要劃好,不然精度會不行,甚至會求解失敗,但這是為什么呢?我們透過表面來看有限元方法的本質,用簡潔易懂的數學推導來展現誤差與單元形狀之間的關系。
本課從實際問題出發,帶著問題去講解有限元中的高斯點與數值積分。一開始拋出了以下3個關鍵問題:
1.對于一個任意函數怎么去得到它的積分?
2.數值積分的本質是什么?為什么簡單地取幾個點就可得到積分值?此種方法的立足點在哪?
3.很多資料上都說“有限元求解精度嚴重依賴于網格質量,過度扭曲的單元會導致結果不收斂或者精度極度惡化”,這只是為什么呢?扭曲單元到底影響的是有限元方法中的哪一步?
圍繞這3個問題,本課分別講了一下三個內容:
1. 數值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
希望有興趣的同學多多支持下,你們的支持是我更新的動力
展開 
ansys中積分點的相關專題、標簽、搜索
ansys中積分點的最新內容
附件下載
聯系工作人員獲取附件
概要
OpticStudio中的坐標間斷是非常靈活的。坐標間斷可用于傾斜或偏心任何光學表面,或光學表面組,圍繞任何軸點,而不干擾光學系統的其余部分。本文將利用坐標間斷來重新定義順序系統的光軸。
簡介
坐標間斷是一個非常通用的工具,可以用來傾斜或偏心一個或多個光學表面。它是非常有用的,能夠選擇光學表面將圍繞什么點旋轉或偏心,我們將在這篇文章中展示如何指定該點
最近突然遇到一個有意思的問題,一時不知道如何操作,想著Ansys 應該比較容易實現,但是用了很長時間才找到一種方案(lll¬ω¬)。不知道大家是如何操作的。
已知:X坐標系和Y坐標系,和A點 相對Y坐標系的位置。查看A點相對X坐標系的位置,A點可以不是幾何點或網格節點。
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為:
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
注:由于技術鄰排版風格有限,故部分內容顯示不全,感興趣的小伙伴可點擊原文進行閱覽:
有限元計算過程中積分點應力如何外插至節點處?【公式推導篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是:有限元計算過程中,單元積分點應力如何外推至節點
寫在前文
盡管減隔震技術與有限元結合取得了眾多成果,但仍面臨諸多挑戰,如材料非線性、模型不確定性等等。減隔震設計除了常規的宏觀結構設計采用SAP2000、Etabs、Midas、SSG、Paco-SAP 或 YJK\PKPM等。
【JY】各類有限元軟件計算功能賞析與探討
我們需要更清楚減隔震元件的破壞模式,對減隔震元件進行破壞分析,除了對減隔震元件在正常工況下的性能進行評估
按照正常的理解,毫無.疑問,abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看,似乎也是如此,每個單元存在8個積分點。
然而,如果自己動手跑一遍程序,就會發現事實遠非如此,采用全積分計算得到的結果與abaqus 存在差異,原因何在?
事實賞,abaqus C3D8 采用的選擇積分方式(selective intergation schema),即對于偏應變,采用
在下半節課中,詳細地分析了扭曲單元與有限元精度之間的關系。我們常聽到單元網格質量要劃好,不然精度會不行,甚至會求解失敗,但這是為什么呢?我們透過表面來看有限元方法的本質,用簡潔易懂的數學推導來展現誤差與單元形狀之間的關系。
本課從實際問題出發,帶著問題去講解有限元中的高斯點與數值積分。一開始拋出了以下3個關鍵問題:
1.對于一個任意函數怎么去得到它的積分?
2.數值積分的本質是什么?為什么簡單地取幾個點就可得到積分值
本課從實際問題出發,帶著問題去講解有限元中的高斯點與數值積分。一開始拋出了以下3個關鍵問題:
1.對于一個任意函數怎么去得到它的積分?
2.數值積分的本質是什么?為什么簡單地取幾個點就可得到積分值?此種方法的立足點在哪?
3.很多資料上都說“有限元求解精度嚴重依賴于網格質量,過度扭曲的單元會導致結果不收斂或者精度極度惡化”,這只是為什么呢?扭曲單元到底影響的是有限元方法中的哪一步
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正
經常有小伙伴問獲取積分點坐標的方法,今天給大家介紹三種獲取積分點坐標的方式,希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標,例如我們熟知的UMAT子程序中包含COORD參數,即為積分點坐標。順帶一提的是,當打開了幾何非線性時,該積分點是當前構形下的坐標,如果未打開幾何非線性則為初始坐標。