基于ANSYS/LS-DYNA的便攜式電子產(chǎn)品的跌落研究
2013-06-06 by:廣州有限元分析、培訓(xùn)中心-1CAE.COM 來(lái)源:仿真在線
為降低開(kāi)發(fā)成本,更好的檢驗(yàn)電子產(chǎn)品抗跌落沖擊性能,本文以數(shù)碼相機(jī)為例,利用有限元軟件ANSYS/LS-DYNA對(duì)其外殼跌落過(guò)程進(jìn)行仿真,得到結(jié)構(gòu)的瞬態(tài)動(dòng)力響應(yīng),清晰的展示在不同時(shí)刻的應(yīng)力、變形情況。此方法可以取代傳統(tǒng)的跌落沖擊試驗(yàn),同時(shí)工程師也可以依據(jù)其設(shè)計(jì)出更好的抗跌落、沖擊產(chǎn)品,提前新產(chǎn)品的上市日期。
作者: 李學(xué)藝*劉鳳景*辛紀(jì)光 來(lái)源: e-works
關(guān)鍵字: CAE 跌落仿真 數(shù)碼相機(jī) 有限元分析
0 前 言
隨著電子技術(shù)特別是小型封裝和高密度組裝技術(shù)的迅速發(fā)展,便攜式電子產(chǎn)品得到了極大的普及。消費(fèi)者在使用或攜帶過(guò)程中將產(chǎn)品失手跌落至地面或不慎撞擊到堅(jiān)硬物體上的意外時(shí)有發(fā)生。跌落或撞擊產(chǎn)生的機(jī)械沖擊可能引起產(chǎn)品嚴(yán)重的機(jī)械或功能損壞。而電子產(chǎn)品往往價(jià)格不菲,故給消費(fèi)者帶來(lái)極大的損失。因此,在滿足基本功能的前提下,消費(fèi)者更加青睞有良好抗沖擊性能的產(chǎn)品。生產(chǎn)廠家在電子產(chǎn)品開(kāi)發(fā)階段一般采用跌落試驗(yàn)機(jī)對(duì)其抗沖擊性能進(jìn)行測(cè)試,但這種實(shí)驗(yàn)花費(fèi)時(shí)間較長(zhǎng)且耗資較大。近年來(lái),由于計(jì)算機(jī)硬件和非線性有限元軟件技術(shù)的迅速發(fā)展,企業(yè)應(yīng)用CAE軟件特別是ANSYS在開(kāi)發(fā)階段對(duì)電子產(chǎn)品進(jìn)行跌落仿真分析已成為一種流行的趨勢(shì),不僅可以有效地提高產(chǎn)品的可靠性,而且能降低開(kāi)發(fā)成本,提高產(chǎn)品的市場(chǎng)競(jìng)爭(zhēng)能力。
本文利用ANSYS/LS-DYNA中的DTM模塊對(duì)一數(shù)碼相機(jī)的外殼進(jìn)行了跌落仿真研究,可以迅速準(zhǔn)確地得到其應(yīng)力、應(yīng)變各種參數(shù),進(jìn)而可對(duì)其安全性進(jìn)行分析。
1 跌落分析動(dòng)力學(xué)分析理論
電子產(chǎn)品跌落是指電子產(chǎn)品在極短的時(shí)間內(nèi),受到碰撞動(dòng)態(tài)載荷作用而產(chǎn)生的復(fù)雜非線性動(dòng)態(tài)響應(yīng)過(guò)程。系統(tǒng)除了具有幾何非線性和材料非線性以外,還有接觸界面的非線性,這是接觸問(wèn)題所特有的。因此,對(duì)電子產(chǎn)品跌落碰撞的過(guò)程進(jìn)行仿真,得到碰撞瞬間的應(yīng)力歷程數(shù)值,一般都采用顯式積分的求解算法。在總體坐標(biāo)系下,建立的碰撞運(yùn)動(dòng)方程可以表示為:
式中:[M]為結(jié)構(gòu)的質(zhì)量矩陣;[C]為結(jié)構(gòu)的阻尼矩陣;[K]為結(jié)構(gòu)的剛度矩陣;{a}為加速度向量;{v}為速度向量;{u}為位移向量;
則碰撞方程可以寫(xiě)為:
若采用集中質(zhì)量,質(zhì)量矩陣[M]為對(duì)角矩陣,則各個(gè)自由度的方程應(yīng)是相互獨(dú)立,故結(jié)構(gòu)系統(tǒng)各節(jié)點(diǎn)在第n個(gè)時(shí)間步結(jié)束時(shí)刻tn的加速度向量為:
然后對(duì)時(shí)間積分求得速度,再積分一次就會(huì)獲得位移,這里采用中心差分的顯式格式來(lái)進(jìn)行時(shí)間積分。節(jié)點(diǎn)速度和位移通過(guò)下面兩式計(jì)算:
時(shí)間步和時(shí)間點(diǎn)的定義為:
在整個(gè)時(shí)域范圍內(nèi),可由上述積分遞推公式求得各個(gè)離散時(shí)間點(diǎn)處的位移、速度和加速度。顯式積分不需要進(jìn)行矩陣分解或求逆,無(wú)需求解聯(lián)立方程組,也不存在收斂性問(wèn)題,因此計(jì)算速度快,其穩(wěn)定性準(zhǔn)則能自動(dòng)控制計(jì)算時(shí)間步長(zhǎng)的大小,保證時(shí)間積分的進(jìn)度。應(yīng)用顯示積分算法求解碰撞問(wèn)題時(shí),一個(gè)特別值得注意的問(wèn)題就是時(shí)間步長(zhǎng)不能超過(guò)臨界時(shí)間步長(zhǎng)。保證收斂的臨界時(shí)步需滿足如下公式:
表1 材料性能參數(shù)表 2)另外,取前殼模型上的一個(gè)節(jié)點(diǎn)node256繪制其各個(gè)方向的位移、速度、加速度以及各種變形結(jié)果曲線圖,觀察其每一個(gè)子步各種參數(shù)的變化情況。如圖5、6分別為節(jié)點(diǎn) node256的沿Y方向上的位移、速度、加速度曲線變化圖。圖8為前殼模型在跌落撞擊過(guò)程中總的能量變化曲線。 相關(guān)標(biāo)簽搜索:基于ANSYS/LS-DYNA的便攜式電子產(chǎn)品的跌落研究 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn) Abaqus培訓(xùn) Autoform培訓(xùn) 有限元培訓(xùn) Solidworks培訓(xùn) UG模具培訓(xùn) PROE培訓(xùn) 運(yùn)動(dòng)仿真
其中,
由此得到的最大特征值即為2.1 建立模型
便攜式電子產(chǎn)品制作精巧,通常結(jié)構(gòu)都比較復(fù)雜。ANSYS軟件雖建模能力比較強(qiáng)大,但是對(duì)于建立精確的電子產(chǎn)品模型來(lái)說(shuō)相對(duì)比較困難。因此,我們?cè)诮⒛P蜁r(shí)通常選擇三維造型功能比較強(qiáng)大的軟件,如PROE、UG等。然后利用ANSYS與這些三維軟件的無(wú)縫集成功能,將其導(dǎo)入到ANSYS。如果導(dǎo)入ANSYS 中的模型過(guò)于復(fù)雜,進(jìn)行合適的有限元網(wǎng)格劃分相對(duì)比較困難,即使劃分好了網(wǎng)格,也可能由于劃分不當(dāng)導(dǎo)致無(wú)法求解。因此,在進(jìn)行碰撞仿真之前,有必要對(duì)其幾何模型進(jìn)行合適的簡(jiǎn)化,以便提高運(yùn)算質(zhì)量。由于CAE和CAD存在本質(zhì)區(qū)別,很多加工工藝上需要的過(guò)渡圓角、小孔以及尖銳的曲面過(guò)渡都會(huì)增加有限元建模的難度,影響單元質(zhì)量,對(duì)于用于碰撞計(jì)算的顯式有限元模型,這個(gè)問(wèn)題尤其突出。所以在進(jìn)行劃分網(wǎng)格之前要對(duì)原幾何模型進(jìn)行修模,除去過(guò)渡圓角、小孔以及使相連的曲面盡量化零為整等。
在對(duì)數(shù)碼相機(jī)的外殼建模的過(guò)程中,綜合考慮各種因素,最終導(dǎo)入到ANSYS/LS-DYNA中的模型的前外殼如圖1所示:
2.2 有限元網(wǎng)格劃分
在有限元分析中,模型中網(wǎng)格劃分越細(xì)致,單元階次越高,計(jì)算結(jié)果越精確,但是要耗費(fèi)大量的CPU時(shí)間。因此,選擇什么樣的單元以及如何劃分網(wǎng)格,是有限元建模中一個(gè)極其重要的問(wèn)題。因ANSYS LS-DYNA 程序中包含的顯示動(dòng)力單元類型較少,且均是線性單元,具有一個(gè)線性位移函數(shù),卻支持所有的非線性選項(xiàng),能夠很好的用于大變形和材料失效的問(wèn)題,所以在選擇單元類型時(shí)相對(duì)簡(jiǎn)單。每種單元具有較多的算法,重點(diǎn)是選擇哪種算法以使求解的結(jié)果更接近真實(shí)。
對(duì)于數(shù)碼相機(jī)的殼體來(lái)說(shuō),通常殼壁較薄,為了使離散結(jié)果更接近于真實(shí)化,可采用充分小實(shí)體SOLID164單元,但是這種離散的結(jié)果是直接導(dǎo)致計(jì)算處理時(shí)間成幾何級(jí)數(shù)上升,因此應(yīng)找到平衡點(diǎn),使得離散結(jié)果接近真實(shí)化,而又不浪費(fèi)資源。在對(duì)外殼的處理過(guò)程中,可用曲面代替殼體,采用SHELL163殼單元,設(shè)定厚度值。另外,設(shè)定總體單元尺寸的大小來(lái)控制網(wǎng)格密度,以便生成較為均勻的網(wǎng)格,劃分時(shí)采用四邊形單元進(jìn)行映射網(wǎng)格劃分。定義殼單元厚度T1為1mm,剪切因數(shù)SHRF為5/6,定義積分點(diǎn)數(shù)NIP為3。離散后前殼模型共有1476單元,1664節(jié)點(diǎn)(如圖2所示)。所用到的材料性能參數(shù)如表1。
2.3 跌落參數(shù)設(shè)置及求解
在進(jìn)行跌落分析之前,必須對(duì)跌落測(cè)試相關(guān)的參數(shù)進(jìn)行設(shè)置,包括重力加速度、跌落高度及參考點(diǎn)、測(cè)試物體的方位、計(jì)算時(shí)間以及輸出文件步數(shù)。設(shè)定加速度與高度時(shí)應(yīng)注意單位之間的協(xié)調(diào),而且物體以何種角度進(jìn)行跌落也是分析中需要設(shè)置的一個(gè)重要參數(shù),要依照具體情況而定。為縮短計(jì)算時(shí)間,可以設(shè)置物體跌落至某一高度的初速度。另外,可以通過(guò)修改目標(biāo)面的特性,使模型與不同的目標(biāo)面發(fā)生跌落碰撞,使某一產(chǎn)品的跌落仿真更全面。
根據(jù)國(guó)家標(biāo)準(zhǔn)規(guī)定跌落高度的優(yōu)先選擇值25 mm、100 mm、500 mm、1 000 mm等。出于安全性的考慮,選取跌落高度為1 000 mm,且均為自由落體,與地面剛性墻撞擊。跌落方案為:前殼為豎直Y方向跌落,無(wú)傾斜角,跌落模型如圖3所示。
有限元建模完成以后,可以輸出關(guān)鍵字K文件,必要時(shí)對(duì)其進(jìn)行修改和完善,后遞交DYNA970求解器進(jìn)行求解;同時(shí)可以設(shè)置結(jié)果文件的輸出格式,便于進(jìn)行多種形式的后處理。同時(shí),顯示算法計(jì)算時(shí)間一般較長(zhǎng),但占用的內(nèi)存較小,可以選擇SW系列開(kāi)關(guān)對(duì)求解過(guò)程進(jìn)行監(jiān)控,遇到問(wèn)題時(shí)需要重新啟動(dòng)求解程序以獲得滿意解。
2.4 仿真結(jié)果分析
LSTC公司專門(mén)開(kāi)發(fā)的后處理功能非常強(qiáng)大,能夠進(jìn)行更加全面的結(jié)果數(shù)據(jù)分析。可以以動(dòng)畫(huà)的形式觀察整個(gè)跌落過(guò)程,繪制各種計(jì)算結(jié)果云圖,如變形、應(yīng)力應(yīng)變?cè)茍D,還可以對(duì)結(jié)果數(shù)據(jù)計(jì)算和分析。在LS-PREPOST中對(duì)數(shù)碼相機(jī)外殼跌落仿真結(jié)果進(jìn)行如下簡(jiǎn)單分析:
1)數(shù)碼相機(jī)前殼在出現(xiàn)最大應(yīng)力時(shí)刻的最大等效應(yīng)力等值線分布如圖4所示:
由材料的沖擊強(qiáng)度為97MPa,而仿真計(jì)算的結(jié)果,前殼在Y向跌落時(shí),最大應(yīng)力為75.08MPa,故前殼的抗沖擊強(qiáng)度足夠。
通過(guò)觀察上述曲線可知,它們符合最基本的自由跌落過(guò)程各參數(shù)的變化規(guī)律,因此 模型仿真結(jié)果與實(shí)際相符。
3 結(jié)論
電子行業(yè)是一個(gè)飛速發(fā)展的行業(yè),市場(chǎng)容量極其巨大,便攜式電子產(chǎn)品已經(jīng)滲透到我們生活的各個(gè)角落。隨著社會(huì)的發(fā)展和技術(shù)的進(jìn)步,人們對(duì)電子相關(guān)行業(yè)提出了更高的要求:精確、穩(wěn)定、輕巧、保密、可靠;同時(shí),電子產(chǎn)品行業(yè)又具有產(chǎn)品更新快,研發(fā)周期短的特點(diǎn)。采用跌落仿真技術(shù)對(duì)其在使用過(guò)程中可能碰到的跌落碰撞情況進(jìn)行研究,以量化的形式展示了產(chǎn)品受沖擊情況。通過(guò)分析仿真結(jié)果,工程師可以依據(jù)其設(shè)計(jì)出更好的抗跌落、沖擊產(chǎn)品,最終降低了成本,加快了新產(chǎn)品的上市速度。