接觸變形分析的案例
WorkBench 接觸分析變形教程<簡單易學>
效果動畫.rar
以下是圖片
以下是動畫
效果動畫.rar
WorkBench 接觸變形分析教程.pdf
密封圈接觸變形仿真 ￥500
密封圈的密封性能取決于密封圈和接觸構件之間的接觸壓力，當密封圈周圍的液體壓力差超過接觸所提供的抵抗力時，發生泄漏，密封圈失效。本案例仿真了密封圈接觸變形及變形回復過程，模擬結果如圖所示：</p><p><img src="https://img.jishulink.com/upload/202110/97c5e68e339e4619bcba887dc372e416.gif" alt="Untitled1.gif"></p><p>感興趣的朋友可以下載模型，也可以加我，歡迎交流</p><p><br></p>
展開 線接觸彈性變形計算
內容：給定一個圓柱線接觸曲線h=x^2/2,接觸區壓力分布為拋物線，pi=ph*sqrt(1-xi^2),計算其彈性變形。
主程序：
clear
clc
N=200;
X1=1.4;
X0=-4.0;
DX=(X1-X0)/(N-1.0);
X(1:N)=0;
H0(1:N)=0;
H(1:N)=0;
P(1:N)=0;
for I=1:N
X(I)=-4.0+(I-1)*DX;
H0(I)=0.5*X(I)^2;
H(I)=H0(I);
if X(I)>=-1 && X(I)<=1
P(I)=sqrt(1-X(I)^2);
end
end
global AK
AK=SUBAK(N);
V=elastic_deformation(N,DX,P,0);
for I=1:N
H(I)=H(I)+V(I);
end
figure(1) %% 圖片
plot(X,V,'-','LineWidth',1.5,'Color',[0 1 1])
hold on
plot(X,P,'-.'
展開 點接觸彈性變形數值計算
contact_deformation.m
main.m
S.m
SUBAK.m
untitled.jpg
最近又做了個案例，是關于點接觸彈性變形數值計算，具體描述參考《彈性流體動壓潤滑數值計算方法》P23-29，黃平著。以下為一個算例的matlab實現方法(原書為FORTRAN語言)，與大家交流，請大家多提意見。
05-橢圓接觸彈性變形Fortran和MATLAB程序 ￥19.89
05-橢圓接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接
03-線接觸彈性變形Fortran和MATLAB程序 ￥19.89
<p><span style="color: rgb(24, 25, 28);">03-線接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接</span></p><div contenteditable="false" width="100%" class="ql-align-center"><jsk id="C_Play608ae0c3c70271f0be504531958d0102" videoid="608ae0c3c70271f0be504531958d0102" duration="16秒"><img src="https://img.jishulink.com/static/web/youku-case.png"></jsk></div><p class="ql-align-center"><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202511/attachment/14991c8788664884bc0328fd5f81e840.png" style="display: inline-block;" data-regular="true">
<img src="https://img.jishulink.com/202511/attachment/14991c8788664884bc0328fd5f81e840.png" data-mobile-src="https://img.jishulink.com/202511/attachment/14991c8788664884bc0328fd5f81e840
展開 04-點接觸彈性變形Fortran和MATLAB程序 ￥19.89
04-點接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接
06-通過多重積分法求解線接觸彈性變形Fortran和MATLAB程序 ￥19.89
06-通過多重積分法求解線接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接
07-通過多重積分法求解點接觸彈性變形Fortran和MATLAB程序 ￥19.89
07-通過多重積分法求解點接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接
CEL與Lagrange模型在大變形分析時的適用性CEL與Lagrange模型在大變形分析時的適用性
對同一個模型來講，通常，拉格朗日建模方式計算更加準確，計算效率更高，因為所有的幾何體都采用拉格朗日單元類型，而CEL建模方式的計算更加耗時，且產生的文件更大，一個直接的原因是流體或大變形幾何體是歐拉體模型，采用歐拉單元建模，而歐拉單元的數量要明顯多于相應的拉格朗日模型的單元數量。
但是，如果模型要經歷極大變形，那么這兩種建模方式的優劣就要好好評價一下了。在大變形分析中，拉格朗日模型容易發生網格畸變，網格畸變區的計算結果準確性將會大打折扣，產生不可信的結果甚至計算中斷得不到結果；而CEL模型在犧牲一定的幾何模型精度和結果準確性的前提下，計算會非常穩定，網格不會發生畸變，相較于拉格朗日的網格畸變區反而會得到更加合理的計算結果。所以，在選擇建模分析方式時，尤其是大變形分析，兩種方法孰優孰劣，需要結合一定的經驗和以往案例，選擇折中處理或者兩種都用以綜合衡量。
本篇案例是一個鉚接案例，如下面的示意圖所示。 ? 具體的模型長下面這樣：左邊是中央截面圖，右面是實物圖，上下兩部分是沖模，張揚帶孔圓盤是固定模板，上下兩部分沖模同時施力以使鉚釘達到最終的變形。 ? 這個過程很明顯是一個極限大變形過程，我們可能關心這個過程中的三個問題：
1、 鉚釘在成型過程中的變形是否適當？
2、 成型后，鉚釘是否有足夠的力量保持材料的連接？
3、 成型過程工具的壓力是否足夠？
那么這三個關心的問題我們可以考察分析鉚釘的變形位移、成型后的等效塑性變形和成型過程中的沖模受力等變量，去評估我們關心的問題從而做出一些結論或改進。 本案例不再進行step by step的演示，各位小伙伴可以自行練習。下面來具體看一下分析模型和相關結果。 ?
左邊是拉格朗日建模，右邊是CEL建模。兩種建模方式中，接觸全部采用無摩擦通用接觸。
展開 【接觸分析】詳述ABAQUS接觸分析（1/2）
在整個分析過程中，都將保持這些關系，絕不會改變主面部分與從面節點的相互作用關系。如果在模型中包括了幾何非線性，小滑動算法將考慮主面的任何轉動和變形，并更新接觸力傳遞的路徑。如果在模型中沒有考慮幾何非線性，則忽略主面的任何轉動和變形，載荷的路徑保持不變。
有限滑動接觸公式要求ABAQUS/Standard經常地確定與從面的每個節點發生接觸的主面區域。這是一個相當復雜的計算，尤其是當兩個接觸物體都是變形體時。在這種模擬中的結構可以是二維的或者三維的。ABAQUS/Standard也可以模擬一個變形體的有限滑動自接觸問題。
在變形體與剛性表面之間接觸的有限滑動公式。不像兩個變形體之間接觸的有限滑動公式那么復雜。主面是剛性面的有限滑動模擬，可以應用在二維和三維的模型上。
4.4 單元選擇
接觸算法的關鍵在于確定作用于從面節點上的力。在ABAQUS/Standard中為接觸分析選擇單元時，一般來說，最好在那些將會構成從面的模型部分使用 一階單元。在模擬中，二階單元有時可能會出現問題，原因是這些單元從常值壓力計算等效節點載荷的方式。為避免問題的發生，在應用于定義從面的任何二階、三維實體或楔型體的單元中，ABAQUS/Standard自動地增加了一個 中面節點。對于常值壓力作用下采用中面節點的二階單元的面上的等效節點力具有相同的符號，盡管這些節點力的量值仍有很大的差異。
對于施加的壓力，一階單元的等效節點力的符號和量值總是保持一致性，因此，由給定節點力的分布所表示的接觸狀態不是模棱兩可的。
如果幾何形狀是復雜的并需要利用自動網格生成器，應該使用在ABAQUS/Standard中的 修正的二階四面體單元（C3D10M），該單元是為了應用在復雜的接觸模擬問題中而設計的，規則的二階四面體單元（C3D10）在其角點處的接觸力為零，導致很差的接觸預測值。
展開 
【接觸分析】詳述ABAQUS接觸分析（2/2）
束縛約束的公式約束了平移自由度和可選擇的轉動自由度，當應用束縛接觸于結構單元時，必須保證任何沒有約束的轉動不會帶來問題。
7.4 初始過盈接觸
為了消除任何初始的過盈，ABAQUS/Explicit將
自動調整在接觸表面上未變形的節點坐標。當應用平衡主從算法時，兩個表面均被調整；當應用單純主從算法時，僅調整從屬表面。為了消除過盈接觸，與調整表面相關的位移不會對在分析的第1個分析步中定義的接觸引起任何的初始應變或應力。
當存在
矛盾的約束時，重新定位節點可能不會完全地解決初始過盈。這種情況下，當采用接觸對算法時，在分析剛開始階段可能會導致網格的嚴重扭曲。通用接觸算法存儲了任何無法消除的初始侵徹，將其作為偏置量以避免過大的初始加速度。
在隨后的分析步中，為了消除初始過盈而進行的任何節點調整都將引起應變，并常常引起網格的嚴重扭曲，因為整體的節點調整發生在一個單一、非常短暫的增量步內。當采用
動力學接觸方法時，這個問題更為明顯，即便是一個非常小的過盈量，都可能引起極大的加速度。通常，在分析步2以及后續的分析步中，所定義的任何新的接觸表面都不能有過盈，這一點非常重要。
八、ABAQUS/Standard和ABAQUS/Explicit的比較
在ABAQUS/Standard和ABAQUS/Explicit中的力學接觸算法具有本質區別，體現在如何定義接觸條件，主要區別如下：
· ABAQUS/Standard在施加接觸約束時應用
嚴格的主從權重，約束從屬表面的節點不能侵入主控表面，而主控表面上的節點原則上可以侵入從屬表面。ABAQUS/Explicit包括這個公式，但是典型地它默認應用
平衡主從權重。
展開 【免費】workbench中橡膠壓縮變形分析-自適應網格+大變形
workbench中橡膠壓縮變形分析
橡膠壓縮是密封圈當中經常遇到的一種現象，但是仿真分析對于橡膠壓縮有很難收斂的現象，本實例通過兩個簡單模擬（公眾號: CAE_ANSYS）
方法說明橡膠壓縮的過程和方法，通過本實例可以了解到以下知識
1.自適應網格的應用方法
2.橡膠材料參數的設置
3.非線性接觸的設置
5.模型建立方法
模型
建立分析模型，如圖所示，本實例以一個簡單模型為例
2.初始網格隨意劃分
3設置自適應網格，采用mesh方法，計算過程自動加密網格，需要注意的是，必須打開大變形，單元必須去除中間節點
4邊界條件采用，向下強制位移的方法
5提取結果（公眾號: CAE_ANSYS）
可以看到網格發生了重新劃分，網格由三角形劃分成了4變形
另一個模型是模擬橡膠壓縮的過程，上下兩個剛性體擠壓中間的橡膠，結果顯示中間橡膠發生大變形
本實例需要注意的是，橡膠材料的設置，不需要設置彈性模量，還有就是接觸的設置，需要選擇相應的線體為接觸面，最好將模型分割，最后獲取相應的結果。
以下模型為兩個模型的計算原始文件，供大家免費參考
（公眾號: CAE_ANSYS）
供大家免費參考，版本為ansys17.2
rubber.zip
展開 #批量cohesive+顯示分析+通用接觸分析沖擊時未接觸就有力了？？
<p>當分析接觸問題時遇到了奇怪的問題，子彈沒有與基體接觸時基體上就有接觸力了，特別是cohesive的力比較大，而基體的力比較小。</p><p>嘗試分析原因：去掉通用接觸后就不會出現這個問題了</p><p>但是，做沖擊分析我們必須設置通用接觸，現在知道是通用接觸導致出現子彈沒接觸基體就出現了應力的問題，但是嘗試修改接觸屬性等參數，目前還是沒有解決這個問題，如果有遇到這類問題的，找到解決辦法的，歡迎一起討論，附上cae文件（本源文件來自星辰北極星，只為尋找解決通用接觸的方法）</p><p><br></p><p>暫時解決辦法：</p><p>解決方法1，建立兩個分析步，沖頭達到基體前建立一個分析步，不設置通用接觸，沖頭接觸基體后建立一個分析步，設置通用接觸，這必須計算好第一步的運行時間，其實第一步沒什么實際意義，只是模擬了飛行過程而已。</p><p><br></p><p>解決方法2，建立一個分析步，建模時候直接讓沖頭與基體接觸，只分析接觸以后的響應過程。
展開 Abaqus帶螺紋螺栓接觸應力分析淺析 Abaqus帶螺紋螺栓接觸應力分析淺析
在實際情況下，很多結構都采用螺栓連接的方式，如何考慮螺栓連接、對連接螺栓的分析計算是個難點。目前的常規做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實應力，圖1為某結構連接螺栓簡化的beam梁應力云圖，沒有接觸應力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網格建模）、計算量大（接觸收斂困難），如圖為某結構帶螺紋螺栓和連接件模型（圖2）和計算得出的結果（圖3）：
圖3 計算結果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結果？
運用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設置跟實際螺紋形狀有關聯的參數，如牙角、螺距、螺栓小徑等，就可以模擬真實的連接螺栓接觸狀況。既可以得到足夠精確的分析結果，又節省了時間專注進行其他的分析設置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結構直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結構直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 