橢圓接觸彈性變形仿真的案例
05-橢圓接觸彈性變形Fortran和MATLAB程序 ￥19.89
05-橢圓接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網盤鏈接
線接觸彈性變形計算
最近做了個案例，具體描述參考《彈性流體動壓潤滑數值計算方法》，黃平著。以下為一個算例的matlab實現方法，與大家交流。
內容：給定一個圓柱線接觸曲線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語言)，與大家交流，請大家多提意見。
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程序，程序請見下文附件及百度網盤鏈接
密封圈接觸變形仿真 ￥500
密封圈的密封性能取決于密封圈和接觸構件之間的接觸壓力，當密封圈周圍的液體壓力差超過接觸所提供的抵抗力時，發生泄漏，密封圈失效。本案例仿真了密封圈接觸變形及變形回復過程，模擬結果如圖所示：</p><p><img src="https://img.jishulink.com/upload/202110/97c5e68e339e4619bcba887dc372e416.gif" alt="Untitled1.gif"></p><p>感興趣的朋友可以下載模型，也可以加我，歡迎交流</p><p><br></p>
展開 基于粘彈性本構模型的熱固性樹脂基復合材料固化變形數值仿真模型
背景介紹
熱固性樹脂基復合材料在制件成型過程中會產生殘余應力，引起固化變形，從而增加裝配和制造的難度，因此，合理預測預制件固化過程中的殘余應力的發展具有重要意義。
早期的研究主要集中于彈性理論來研究復材的固化成型，現今，越來越多的文獻考慮了樹脂的固化放熱以及材料的各向異性等因素的影響，發展了基于粘彈性模型的數值仿真計算方法，證明了粘彈性的結果固化變形量小于線彈性的結果，且樹脂含量越高的復材，其粘彈性效果越明顯。
RTM成型工藝示意圖
二。粘彈性模型在Abaqus中的實現
本文作者在參考文獻【1】的基礎上，使用廣義Maxwell粘彈性本構模型，聯合編寫了HETVAL、USDFLD、DISP、UMAT及UEXPAN子程序，在abaqus軟件平臺中實現了復材固化成型的仿真模擬，其基本編程思路如下圖所示：
其中，最關鍵的粘彈性本構公式為：
參考上述公式和子程序的編寫流程，可以完成上述模型。最后得到仿真Mises應力云圖和S33云圖如下：
得到的S33關于時間的曲線趨勢如下所示：
該曲線結果和文獻有出入，但是榮的文獻中關于底數的取值有錯誤，亦即下列公式的底數應以e為底數，而不是10
【1】
基于黏彈性本構模型的熱固性樹脂基復合材料固化變形數值仿真模型.pdf
最后，歡迎大家關注“320科技工作室”微信公眾號，有相關需求可以添加管理員聯系方式~
展開 