ansys彈性計算的案例
通過ansys利用均勻化理論計算復合材料等效性能--等效彈性模量，剪切模量等
***********計算NIX,NIY,應變矩陣元素
EPX=0
EPY=0
*DO,K,1,4 !***********計算插值點應變
EPX=NIX(K)*V(K)+EPX
EPY=NIY(K)*U(K)+EPY
*ENDDO !***********計算插值點應變
PX=PX+EPX
PY=PY+EPY
*ENDDO !***********計算單元應變PX,PY
*GET,DA,ELEM,M,AREA !提取單元的面積
K=EPSI/(1+NU)/2/AREA
E1212H=K*(DA-PY*KA-PX*KA
)+E1212H
*GET,M,ELEM,M,NXTH
*ENDDO
*enddo
allsel,
finish
G12=E1212H
展開 超彈性與亞彈性，顯式與隱式的HCP多晶滑移+孿晶（主導孿晶重定向（PTR））計算效率比較
參考文獻：《Influence of texture distribution in magnesium welds on their non-uniform mechanical behavior: A CPFEM study》
主導孿晶重定向（PTR）方案作為目前處理HCP晶格結構的多晶材料孿晶模擬中最常使用的方案被廣泛討論，然而晶體取向旋轉過程可能會造成模擬的收斂性問題，選擇一個相對穩定的本構框以及迭代變量對模擬計算效率的提升是有意義的。文章中公式以及其對應的參數總結如下：
這里使用文章的模型和參數對超彈性和亞彈性PTR方案進行比較。
二維（200個晶粒X方向壓縮20%）
以下各個圖中左圖為超彈性結果，右圖為亞彈性結果：
應力分布云圖
應變分布云圖：
孿晶分布云圖：
這里使用文章的模型和參數對顯示和隱式PTR方案進行比較
二維（200個晶粒Y方向剪切變形20%）
以下各個圖中左圖為顯示結果，右圖為隱式結果：
應力分布云圖
應變分布云圖：
孿晶分布云圖
拉壓非對稱與織構演化方面超彈性與亞彈性保持一致：
初始極圖：
RD拉伸20%：
RD壓縮20%：
應力應變曲線
模擬的結果建議，使用PTR方案，超彈性建議使用PK2應力和當前強度為迭代變量，并使用雙重迭代方案，亞彈性建議使用柯西應力為迭代變量，兩者在模擬過程中，計算效率相差較小，無論是局部晶粒的應力應變響應，整體的流動應力，以及變形后的織構結果幾乎保持一致。同時涉及到接觸，碰撞問題，修改為顯式對于收斂性的提升是必要的。
展開 基于ABAQUS超彈性材料橡膠襯套的剛度計算 附基于Abaqus的橡膠和粘彈性建模下載
橡膠材料作為一種具有可逆形變的高彈性、高分子聚合物材料，基于其在彈性特性方面所具有的超彈性與粘彈性一直被廣泛應用于各個工程領域的減振制品中。對于一些結構簡單的橡膠制品，我們可以基于一些理論推導或工程經驗算法在設計初期來獲取其靜剛度特性。但由于橡膠具有非線性粘彈性與超彈性，這種理論計算結果往往與試驗存在一定誤差，并且這種誤差在一般情況下是不可以忽略不計的，其具有一定的工業應用價值。
為減小誤差或實現零誤差的前期預測，我們引入了有限元仿真分析技術，其可以通過控制模型參數與網格質量實現較小誤差的預測計算。其價值也在各個行業實際的生產中得到了很好的驗證。本文基于減振襯套簡單講訴一下基于ABAQUS軟件的橡膠制品靜剛度仿真分析過程。
仿真分析過程可分為三個大過程：前處理、求解計算和后處理。本文基于ABAQUS軟件設定的分析步驟，不再重點區分分析的三個過程，將操作過程拆分為：部件、屬性、裝配、分析步與輸出設置、相互作用、網格、加載、作業提交與監管以及計算結果的可視化處理九個模塊，下面講訴橡膠襯套靜剛度仿真分析過程。
一、部件
由于本文主旨是為介紹橡膠剛度仿真的過程，所以選用了結構較為簡單的橡膠襯套為例，直接借助ABAQUS軟件的部件模塊常見如圖1所示的幾何模型。
圖1、幾何模型結構圖
二、屬性
為了使仿真結果更接近與實驗值或真實值，除了需要一個適合的仿真求解器和一個高質量的網格文件，更需要選擇一個合適的橡膠本構模型，在ABAQUS軟件中內置了許多相對成熟的橡膠本構模型（如圖2所示），我們可以通過指定相關的系數來實現本構模型的定義,當然我們還可以直接提交我們的試驗數據，交由ABAQUS軟件進行擬合，得出相對精準的參數。
展開 交流-ANSYS橡膠材料超彈性本構模型和粘彈性性能仿真和試驗
交流-ANSYS橡膠材料超彈性本構模型和粘彈性性能仿真和試驗
最近在搞橡膠這個方向，單軸拉伸試驗和動態DMA，研究橡膠次本構模型
有研究橡膠超彈性。粘彈性性能的朋友可以聯系，互相交流學習、答疑。
Q254958758
線接觸彈性變形計算
最近做了個案例，具體描述參考《彈性流體動壓潤滑數值計算方法》，黃平著。以下為一個算例的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,'-.'
展開 等效彈性模量計算公式
誰知道不同材料的板粘接在一起后，等效彈性模量計算公式！！比如方艙壁板它由兩層1.5mm的鋁板夾49mm厚的泡沫板粘接在一起，這種板的等效彈性模量計算公式？？那位高手知道請告訴我！！謝謝！！
點接觸彈性變形數值計算
contact_deformation.m
main.m
S.m
SUBAK.m
untitled.jpg
最近又做了個案例，是關于點接觸彈性變形數值計算，具體描述參考《彈性流體動壓潤滑數值計算方法》P23-29，黃平著。以下為一個算例的matlab實現方法(原書為FORTRAN語言)，與大家交流，請大家多提意見。
采用VASP如何計算晶體的彈性常數Cij
最后，有相關需求歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
粘彈性邊界施加及等效節點力計算施加結果 ￥10
粘彈性邊界施加及等效節點力計算施加程序 VBEA；下圖為模擬結果，文字為軟件部分說明書，付費內容僅為軟件完整使用說明書，不提供軟件）
Table of Contents
1. Brief Introduction.
2. Installation.
3. Input Files Preparation.
3.1 Earthquake motion file.
3.2 Node number file.
4. Run the program..
5. Interpretation of output files.
1. Brief Introduction
Viscous-Spring Boundary and Equivalent Node Load Application is a program that automatically applies viscoelastic boundary conditions and inputs seismic wave with equivalent nodal load. It has a graphical user interface that is easy to use and understand, and data visualizations through the program are made possible thanks to the advanced plotting capability of MATLAB.
展開 如何在ANSYS中擬合橡膠材料曲線？ 附Ansys橡膠材料的粘彈性本構模型下載
下載地址：Ansys橡膠材料的粘彈性本構模型
ANSYS知識普及系列15——粘彈性人工邊界在ANSYS中的實現
本人準備出一個ANSYS知識普及系列，將有用的網上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術鄰ANSYS專家
業務咨詢網址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網上；
2、如侵犯知識產權，請聯系ANSYS專家本人或者技術鄰，我將第一時間刪除。
小技巧：加本人關注，可以及時觀看本人發布的技術貼
從半空間無限域取一4X2的矩形平面結構，頂部中間一定范圍內受隨時間變化的均布荷載，荷載如下
p(t)=t
當0< DIV>
p(t)=2-t
當1<=t<=2時
p(t)=0
當t>2時
材料彈性模量E=2.5，泊松比0.25，密度1
網格尺寸0.1X0.1，在網格邊界上所有結點加法向和切向combin14號單元用以模擬粘彈性人工邊界（有關理論可參考劉晶波老師的相關文章）。combine14單元的兩個結點，其中一個與實體單元相連，另一個結點固定。網格圖如圖1所示
時程分析的時間步長為0.02秒，共計算16秒。計算得到四個控制點位移時程圖如圖2所示，控制點坐標A(0,2)、B(0,1)、C(0,0)、D(2,2).
計算所用命令流如下：
/PREP7
L=4
!水平長度
H=2
!豎起深度
E=2.5
!彈性模量
density=1
!密度
nu=0.25
!泊松比
dxyz=0.1
!
展開 
如何在ANSYS中擬合橡膠材料曲線？ 附Ansys橡膠材料的粘彈性本構模型下載
下載地址：Ansys橡膠材料的粘彈性本構模型
計算方法 | 淺析橡膠超彈體與粘彈性仿真(超彈篇)
付穌昇
仿真 xiu專欄作者
本文主要依據個人之前學習和工作積累進行橡膠類超彈體材料本構模型在CAE仿真計算的技術簡要整理，由于非此類科學技術計算專業工程人員，唯恐出錯以致誤導，誠懇大家辨別學習，但還是愿意以此種方式進行編寫分享以能給那些曾經如我一樣沒人指導和參與培訓學習的朋友。
很早之前就曾想過把CAE仿真中對于超彈體和粘彈性計算的準備和部分流程方法進行一個分類整理并分享出來。但自身僅做為機械工程普通從業人員，又不是專業的高分子材料仿真從事人員，底氣始終是少了一點，即使曾經很長一段時間自學過相關理論和做過不少相關計算模型，所幸主導過百萬費用級橡膠密封計算項目，但始終怕誤導剛入門這方面仿真的朋友，遲遲不敢下筆。
幾天前，仿真 xiu平臺一句“如果說是真心的分享技術，那么讀者自然會有自己的理解和判斷，無須擔心”，打消了我的顧慮。所以以這種心態利用閑散時間對超彈體和粘彈性計算的CAE仿真初步工作做一個整理，若是對部分讀者有所幫助，那就是莫大的欣慰。本文分為兩個篇幅，第一個篇幅進行橡膠類超彈體本構仿真計算的內容簡述，另一個篇幅簡述粘彈性仿真計算的準備工作。
展開 有關polyflow計算粘彈性本構方程不收斂情況簡要分析
粘彈性本構方程是研究聚合物的流動性質，polyflow提供了積分型和微分型本構方程，對于初學者在運用這兩種方程的時候經常會遇到一些收斂的問題(滿足網格質量要求情況下)，下面我們簡單來分析一下這其中的原因，為了更好的說明這種現象，研究從KBKZ積分粘彈性方程來說明。
polyflow中KBKZ粘彈性方程
K-BKZ模型能夠很好地描述粘彈性流體剪切變稀，拉伸黏度，以及彈性方面的第一第二法向應力差，其方程中附加應力張量可分為兩個部分：T1黏彈部分，T2純黏部分
其中m(t-s)是記憶函數，反映材料的時間依賴性；i指的是第i個松弛模量，H是阻尼函數，θ是控制法向應力差比值的一個標量
在polyflow中需要定義時間松弛譜，我們定義6個松弛時間對分別如下
物理模型（全長尺寸大概200mm左右）
邊界條件
入口速度100mm/s(紅色)
計算結果
是不是很蛋疼…………………………？是的。
簡要分析：t流動≈200/100=2s,也就是說聚合物在該區域中的流動時間最多為2s（按照壁面無滑移來說的話壁面上的聚合物速度為0），那么對于松弛時間譜上1.999和2.999這兩個時間的話，polyflow到底有沒有參與計算呢？有點懷疑。因此把松弛時間譜的個數降為4個的情況繼續算.
驚奇的發現，在去掉了2個貌似不合理的時間松弛譜之后，計算收斂了。有點讓人費解，為了研究的方便，我們取兩端的壓力降來研究。當然了這過程中涉及到時間松弛譜個數的選擇。
那么我們的懷疑的對象該不該指向這個polyflow處理時間松弛譜上呢？
展開 《軌道車輛橡膠彈性元件設計計算方法>中的fesafe的讀后感
軌道車輛橡膠彈性元件設計計算方法 model.zip
剛買了《軌道車輛橡膠彈性元件設計計算方法》這本書的，
這本書的第六章里面有關于fesafe分析膠芯壽命的講解。
讀了這本書，關于fesafe的講解，還不如網上的教程講的詳細
哎188頁，花了我差不多五十元錢
附件為這本書的光盤中的模型為第4章“橡膠彈性元件設計技術例解”中的案例結構模型，具體模型文件與書中的案例對應關系如下：
4.1 橡膠堆旁承設計: bearing.x_t
4.2 橡膠球鉸設計：bush.dxf
4.3 軸箱橡膠彈簧設計：axel-spring.dxf
4.4 抗側滾扭桿設計：torsing-bar.x_t
4.5 空氣彈簧設計：airspring-model.dxf
提供了橡膠堆旁承設計中的橡膠堆、壓板工裝等結構的三維文件
展開 