對于聚合物合金及其復合材料，其力學性能在很大程度上取決于界面（如黏合、結構形態和應力分布等）。一方面，界面起到了傳遞顆粒與基體之間應力的作用；另一方面，界面也決定了復合材料在載荷作用下的斷裂方式。應力分析是有限元分析法最常見的應用之一。
本工作擬應用有限元方法對POM/EVA/CaCO3復合材料拉伸載荷下界面應力分布進行二維數值模擬，考察復合材料單元中相鄰EVA粒子和CaCO3粒子之間的相互作用對界面處應力及其分布的影響。
一、基本假設
（1）EVA粒子與CaCO3粒子在POM基體中的分布是均勻的，且為簡立方分布；
（2）EVA粒子與CaCO3的粒徑都很小，可近似視為球形，且兩個粒子的直徑相同；
（3）材料在彈性范圍內變形，應力與應變關系遵循胡克定律。
二、物理模型
現考察v(EVA)和v(CaCO3)均為10%的POM/EVA/CaCO3復合材料在拉伸載荷（F）作用下的情形。根據上述假定，取中央僅含有一個球形EVA粒子和CaCO3粒子的正方形單元體，上下兩端受簡單拉伸作用，其物理模型如圖5-11所示，其中，左邊為EVA粒子，右邊為CaCO3粒子。
 
 
 
 
 
三、界面正應力分析
圖5-17中橫坐標原點設在粒子赤道的左端點（L=0μm處），L為從相應粒子的最左第五章拉伸過程中界面應力的有限元模擬端開始算起逆時針方向的弧長，而L=2.20μm處代表粒子的極點。從圖4-17看出：EVA粒子界面正應力分布存在不對稱性。在坐標原點附近，X軸方向的界面拉伸應力不斷增大，且在L=0.73μm處拉伸應力達到值，并在附近形成應力集中。而后，拉伸應力緩慢減小，直到L=1.07μm，拉伸應力開始迅速減小，并在L=1.28μm臨界處應力為零。
然后，應力開始轉為負值（受壓），壓縮應力隨之迅速增加，且在L=2.20μm極點處達到值，在極點位置附近形成應力集中。在粒子的右半側，從極點開始，壓縮應力迅速減小，在L=3.01μm臨界處，應力開始轉為正值，受拉伸應力，在L=3.54μm位置處達到值，即拉伸應力值。之后拉伸應力逐漸減小。
這說明在粒子的右側，由于受到CaCO3粒子的作用，使其在X軸方向的拉伸應力變大，同時，赤道處的應力集中區域有所擴大。從圖5-17還看出：在EVA粒子赤道的坐標原點處的Y軸方向應力，然后朝粒子極點方向迅速減小，并于L=1.62μm位置達到極小值，然后緩慢增加，到極點處不再增加。從EVA粒子極點開始，應力緩慢減小，并與L=2.56μm處達到極小值，之后應力開始以較快速度平穩增加，直到L=3.84μm處出現轉折，應力迅速增加，在L=4.40μm處到達極大值。這是因為與CaCO3粒子相鄰的一側，由于CaCO3粒子對其的影響，使其在Y軸方向的拉伸應力變小，同時，極點處的應力集中區域有所收縮。
 


 
 
從圖5-18看出：CaCO3粒子X軸方向界面應力分布情況同樣存在明顯的不對稱性。在坐標原點附近，其應力值，為正值。之后，應力迅速減小，并在臨界點L=1.02μm處應力值由正變負，壓縮應力增大，直到L=1.06μm，增大速率趨于平緩，且在L=1.35μm處達到極大值。而后，在L=1.46μm處，壓縮應力轉為拉伸應力，并迅速增大，直到在L=2.20μm極點處達到。在粒子右側，從極點開始，應力以較小的速率減小，在臨界點L=3.24μm處應力值由正變負，壓縮應力增大，在L=3.74μm處到達值。
而后，壓縮應力緩慢減小，在L=3.74μm，壓縮應力轉為拉伸應力，并緩慢增加，直到赤道終點達到應力值。這表明EVA粒子對CaCO3粒子界面應力及分布有明顯的影響。從圖5-18還看出：Y軸方向的分布情況也存在明顯的不對稱性。在赤道坐標原點處，拉伸應力先增大，在L=1.35μm處達到值，而后逐漸減小，直到極點處。但是，在粒子的右側，拉伸應力始終保持減小的趨勢。這表明由于兩粒子之間的相互作用，使與EVA粒子相鄰的赤道區形成了應力集中的現象，且CaCO3粒子的影響更大。



 
四、界面剪切應力分析
從圖5-19看出：在L=0.37μm位置，EVA粒子正方向的剪切應力，并在附近區域形成應力集中；在L=3.82μm位置，負方向的剪切應力，并在附近區域形成應力集中。這表明，EVA粒子靠近CaCO3粒子的一側，由于受到CaCO3粒子的牽制作用，其剪切應力集中區域向CaCO3粒子方向偏移了。
從圖5-19還看出：CaCO3粒子剪切應力正負方向與EVA粒子剛好相反。CaCO3粒子在L=0.37μm處，負方向的剪切應力，并在附近區域形成應力集中；L=3.45μm處，正方向的剪切應力，并在附近區域形成應力集中。CaCO3粒子的剪切應力分布不對稱性更加顯著，這表明CaCO3粒子的界面剪切應力分布受兩粒子之間的相互作用的影響更大。同時，由于在拉伸載荷下基體樹脂與剛性粒子和彈性粒子之間界面的剪切應力正負方向相反，正、負剪切應力正好融合在一起，使剪切應力在赤道中間形成明顯的應力集中。
 

 
 
POM/EVA/CaCO3復合材料承受拉伸載荷時，產生相應的拉伸形變。在此情形下，EVA粒子和CaCO3粒子在垂直于拉伸方向（即X軸方向）受到樹脂基體的擠壓作用，產生橫向和縱向的形變。由于EVA的彈性模量遠小于POM樹脂，故EVA粒子的極點區受到樹脂基體的擠壓，應力為負值；而赤道區則受到樹脂基體的提拉，故應力為正值，因此，EVA粒子及周邊的基體在拉伸載荷下發生了較大的形變。與此相反，CaCO3的彈性模量遠大于POM樹脂，故粒子的赤道區受到樹脂基體的擠壓，而極點區受到樹脂基體的提拉；但同時受相鄰EVA粒子的相互作用較大，與EVA粒子相鄰的赤道區也受到提拉作用，故CaCO3粒子的極點區及與EVA粒子相鄰的赤道區的界面應力均為正值。
此外，由于EVA粒子與CaCO3粒子間存在相互作用效應，且CaCO3粒子的彈性模量遠大于EVA粒子，故兩粒子之間于赤道區處，尤其是靠近CaCO3粒子的區域，其正應力集中區域明顯擴大。
5.4.4結論
a)在POM/EVA/CaCO3復合材料的拉伸過程中，EVA粒子界面受到的剪切作用區域都偏離了極點區和赤道區。
b)EVA粒子與CaCO3粒子間界面存在明顯的正應力集中，尤其是在CaCO3粒子靠近EVA粒子的周圍，應力集中面積大；由于粒子之間的相互牽制作用，剪切應力區域都偏離了45°，且呈現不對稱性。剪切應力和壓縮應力都發生在EVA粒子界面上。此外，兩粒子之間的相互作用對CaCO3粒子界面應力及分布的影響更大。
應用ANSYS軟件對POM/CaCO3復合材料、POM/EVA/CaCO3復合材料在拉伸過程中界面應力及其分布進行了數值模擬。結果表明：在對POM/CaCO3復合材料在拉伸載荷作用下的界面應力進行數值模擬，發現相鄰碳酸鈣粒子之間存在的相互作用效應，對界面應力及其分布有明顯的影響。剪切應力區域都偏離了45°，相鄰粒子之間的赤道區的界面應力較小。正應力發生在粒子的極點處。在對POM/EVA/CaCO3復合材料在拉伸載荷作用下的界面應力進行數值模擬，發現POM、EVA粒子界面受到的剪切作用區域都偏離了極點區和赤道區。
同時，EVA粒子與CaCO3粒子間界面存在明顯的正應力集中，尤其是在CaCO3粒子靠近EVA粒子的周圍，應力集中面積大；由于粒子之間的相互牽制作用，剪切應力區域都偏離了45°，且呈現不對稱性。剪切應力和壓縮應力都發生在EVA粒子界面上。此外，兩粒子之間的相互作用對CaCO3粒子界面應力及分布的影響更大。http://www.hbard.cn