ANSYS與ABAQUS比較之實(shí)例4---圓壓頭與平板的接觸分析2

2016-09-16 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

本文是ANSYS與ABAQUS比較之系列篇,本文是第四篇,關(guān)注的是在接觸分析方面二者的異同。

由于分析比較復(fù)雜,該比較分為兩篇來(lái)說(shuō)明。上篇1是使用ABAQUS進(jìn)行求解的過(guò)程,本篇2則是用ANSYS求解的過(guò)程,比較的結(jié)果將在本篇2中給出。

--------------------------------------------------------------------------------------

1.問(wèn)題描述

一半徑為10mm的圓形薄片,壓在一90mm*30mm的矩形板頂邊中間。在圓片上施加一個(gè)6KN的集中力,使得圓片壓在矩形板上,現(xiàn)在要求分析模型的受力狀態(tài)。

已知:矩形板材料為鋼材,彈性模量為210GPA,泊松比為0.3;圓形薄片相當(dāng)剛硬,可以看做是剛體;在圓片和矩形板之間沒有摩擦。

(該問(wèn)題來(lái)自于張建華,丁磊的《ABAQUS基礎(chǔ)入門與案例精選》,電子工業(yè)出版社,2012.6)

---------------------------------------------------------------------------------------

2.分析思路

(1)本問(wèn)題是一個(gè)靜力學(xué)問(wèn)題,而且屬于一個(gè)平面應(yīng)力問(wèn)題,對(duì)于板使用平面應(yīng)力單元,對(duì)于圓形薄片則使用剛體。

(2)該問(wèn)題還是一個(gè)接觸分析問(wèn)題,它是非線性問(wèn)題的一種。對(duì)于接觸設(shè)置為無(wú)摩擦的接觸。

(3)由于結(jié)構(gòu)關(guān)于Y軸對(duì)稱,為了提高計(jì)算效率,可以進(jìn)一步只取右半邊來(lái)分析。

3.步驟

(1)打開ANSYS WORKBENCH,創(chuàng)建靜力學(xué)分析模塊,并設(shè)置為2D分析類型。

(2)創(chuàng)建材料屬性,設(shè)置彈性模量為2.1e11Pa,泊松比為0.3.

(3)創(chuàng)建幾何模型

(4)賦予材料屬性,長(zhǎng)方形為上面創(chuàng)建的材料,圓形為剛性體,設(shè)置厚度為0.001m,與上篇一致。

(5)設(shè)置對(duì)稱約束,對(duì)稱的法線方向?yàn)閄軸

(6)創(chuàng)建遠(yuǎn)程點(diǎn),創(chuàng)建遠(yuǎn)程點(diǎn)是為下文施加載荷所用。遠(yuǎn)程點(diǎn)關(guān)聯(lián)的幾何體選擇半圓面,位置坐標(biāo)如下圖所示。

(7)設(shè)置接觸,接觸類型為無(wú)摩擦接觸,保持默認(rèn)的接觸算法等設(shè)置。

(8)劃分網(wǎng)格,設(shè)置面網(wǎng)格尺寸為0.0015m,與上篇博文里的網(wǎng)格尺寸一致

(9)施加約束,對(duì)長(zhǎng)方形下面的線進(jìn)行固定約束。

(10)施加載荷,施加圓形剛性體豎直向下的集中載荷6KN,這時(shí)我們發(fā)現(xiàn),載荷對(duì)象無(wú)法選擇剛性體,這就需要借助創(chuàng)建的創(chuàng)建遠(yuǎn)程點(diǎn),載荷類型選擇遠(yuǎn)程力,設(shè)置如下圖。

(11)求解設(shè)置,打開大變形,其余保持默認(rèn)設(shè)置。

(12)進(jìn)行求解.

(13)求解完成后查看后處理,查看變形和應(yīng)力,這時(shí)我們發(fā)現(xiàn),雖然可以求解完成,但是半圓形剛性體的位移為2038.6m,產(chǎn)生了剛性位移,也就是說(shuō),半圓和長(zhǎng)方形并沒有接觸上。接著將需要修改接觸設(shè)置,先修改接觸表面類型,設(shè)置為Adjust to touch,也就是允許調(diào)整模型的間隙到剛剛接觸,然后再次進(jìn)行計(jì)算。

(14)查看變形,變形最大值為0.17534mm

(15)查看應(yīng)力,應(yīng)力最大值為902.81MPa

4.Abaqus與ANSYS的結(jié)果分析及討論

(1)從收斂上看,二者都能快速收斂;

(2)從結(jié)果上看,由上篇在ABAQUS中計(jì)算得出的應(yīng)力為606.5MPa,而在ANSYS中得出的應(yīng)力為902.81MPa,ANSYS計(jì)算的米塞斯應(yīng)力是ABAQUS計(jì)算的米塞斯應(yīng)力的1.5倍,這相對(duì)誤差是比較大的,至于哪個(gè)軟件計(jì)算的值更接近真實(shí)應(yīng)力,此時(shí)無(wú)法比較,因?yàn)椴牧弦呀?jīng)屈服,除非設(shè)置材料的塑性行為,且與實(shí)驗(yàn)做比較才知道誰(shuí)更接近真實(shí)值,有興趣的朋友可以對(duì)比下。

開放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才