【摘要】：在基于物理的仿真領域,近幾年流行的物質(zhì)點法已經(jīng)成為一種越來越重要的仿真方法.物質(zhì)點法已經(jīng)可以仿真雪、物質(zhì)狀態(tài)變化、顆粒狀的沙子以及粘彈性的泡沫等物理現(xiàn)象.現(xiàn)有的各種復雜物理現(xiàn)象的模擬已經(jīng)展示了物質(zhì)點法作為一種通用的求解方法的潛質(zhì),然而對于仿真形變體尤其是自然界普遍存在的非均質(zhì)材質(zhì)還存在一些不足.當同一網(wǎng)格內(nèi)部包含多種材質(zhì)的粒子時,由于背景網(wǎng)格與粒子之間的插值只考慮粒子的幾何位置而不包含物質(zhì)屬性信息,導致不同屬性的粒子在網(wǎng)格內(nèi)部表現(xiàn)為"平均化"的運動趨勢.為了改進傳統(tǒng)物質(zhì)點法模擬多種材質(zhì)的缺陷,該文提出了一種模擬非均質(zhì)彈性材料的仿真方法.首先,在預處理過程中將仿真物體離散表示成粒子,針對物質(zhì)的不同屬性粒子的分布,建立物質(zhì)邊界并確定邊界粒子,仿真過程中動態(tài)更新邊界粒子的位置;其次,引入粒子影響域,區(qū)分網(wǎng)格內(nèi)部具有不同屬性的粒子,增加網(wǎng)格內(nèi)部計算自由度;最后,提出一個判斷準則,根據(jù)物質(zhì)邊界粒子的位置判定仿真粒子是在背景計算網(wǎng)格上求解還是粒子區(qū)域上求解.實驗結果表明,利用改進的物質(zhì)點法能夠有效地模擬非均質(zhì)材質(zhì)彈性體的形變.與有限元法相比,該文方法在模擬非均質(zhì)彈性體時可以獲得相似的效果,并且在模擬拓撲結構變化的問題以及碰撞處理上更有優(yōu)勢.
【圖文】：

子的體積和密度，表示為ρ０ｐ＝∑ｉｍ０ｉω０ｉｐ／ｈ３，Ｖ０ｐ＝ｍｐ／ρ０ｐ（１８）其中，ρ０ｐ、Ｖ０ｐ分別表示粒子ｐ在初始時刻的密度和體積，ｍ０ｉ表示第ｉ個背景計算網(wǎng)格在初始時刻的質(zhì)量，ｈ為網(wǎng)格間距．整個預處理過程包含建立粒子區(qū)域、計算嵌入表面網(wǎng)格與粒子區(qū)域結點的插值權重等操作，消耗時間根據(jù)粒子數(shù)不同有所差異，本文的所有實驗都只需要幾分鐘就可以完成．圖６預處理階段示意圖（２）仿真過程（如圖７所示）．①仿真過程中，每個時間步利用粒子區(qū)域的權重函數(shù)將粒子的所有物理量映射到粒子區(qū)域；②根據(jù)５．１節(jié)提出的算法，查找增強物質(zhì)點，動態(tài)判斷仿真粒子狀態(tài)．對于增強的粒子在粒子區(qū)域上求解，對于普通的粒子，將粒子區(qū)域上的物理量映射到背景網(wǎng)格上，更新背景網(wǎng)格以及粒子區(qū)域上的速度等信息；③根據(jù)網(wǎng)格結點處的外力以及第一Ｐｉｏｌａ－Ｋｉｒｃｈｈｏｆｆ應力Ｐ，判斷粒子的類型分別在背景計算網(wǎng)格或粒子區(qū)域上積分求解，得到新一時刻網(wǎng)格和粒子區(qū)域結點處新的速度，速度更新采用文獻［１２］的半隱式方式計算．力的計算如下式：ｆｉ＝－∑ｐＶ０ｐＰｐ鄀ωｎｉｐ（１９）ｆα＝－∑ｐＶ０ｐＰｐ鄀Ｎαｎｐ（２０）其中，ｆｉ、ｆα分別表示背景計算網(wǎng)格和粒子區(qū)域結點處的力，鄀ωｎｉｐ、鄀Ｎαｎｐ分別表示第ｎ時刻粒子與網(wǎng)格結點ｉ的插值函數(shù)的梯度、第α個粒子區(qū)域結點的插值函

圖７仿真過程中物理量投影和計算示意圖外計算開銷可以接受．本文所有仿真實驗的時間效率詳見實驗部分表１．６實驗結果與分析實驗在Ｌｉｎｕｘ操作系統(tǒng)下，使用Ｃ＋＋語言，利用Ｐｈｙｓｉｋａ①開源的物理仿真框架，Ｐｏｖ－ｒａｙ離線渲染．硬件系統(tǒng)為：單核Ｉｎｔｅｌｉ５ＣＰＵ，２．８ＧＨｚ，ＮｖｉｄｉａＧｅＦｏｒｃｅＧＴ４３０顯卡．所有的實驗粒子密度均為１．０×１０３ｋｇ／ｍ３，除圖１４中雪的本構模型外，所有非均質(zhì)彈性體的泊松系數(shù)均為０．３，所有實驗都是每３００個時間步輸出一次數(shù)據(jù)．６．１有效性實驗圖８是一個分層的立方體的仿真結果，立方體被放置在地面上，下端固定，上端施加向下的壓力．模型被分為兩種屬性，隔層的物質(zhì)具有相同屬性．相對透明的材質(zhì)較軟，楊氏模量Ｅ＝１．０×１０５，其余層材質(zhì)較硬Ｅ＝１．０×１０７．圖８分層的立方體的壓縮實驗左側的上下兩張圖用嵌套在粒子區(qū)域內(nèi)的精細的表面網(wǎng)格繪制，分別表示初始形狀和在力作用下的形變形狀，很好的表現(xiàn)了兩種不同材質(zhì)的壓縮形變．右側兩幅圖，上面的表示初始時刻構建的物質(zhì)邊界的示意圖（用粒子球的形式表示，不同材質(zhì)臨近區(qū)域的粒子為邊界粒子），下方的圖對應其左側效果圖的形變時刻，邊界粒子球所在網(wǎng)格中的粒子為在粒子區(qū)域上求解的增強的粒子，黑色線框表示求解的背景網(wǎng)格．由于這個例子展現(xiàn)的是壓縮形變，，不同的粒子以及粒子邊界都被壓縮到比較緊湊的背景網(wǎng)格內(nèi)．根據(jù)本文的追蹤物質(zhì)邊界方法以及增強的粒子的設置策略，物

【相似文獻】

相關期刊論文 前10條

1 趙靜;唐勇;李勝;汪國平;;基于增強的物質(zhì)點法的非均質(zhì)彈性材料仿真方法研究[J];計算機學報;2018年11期

2 朱明英;姚文勝;邢豫;陳春華;;基于功能點法的軟件工作量評估方法研究與應用探討[J];現(xiàn)代計算機(專業(yè)版);2017年33期

3 駱巖林,汪國昭;控制四點法“形狀參數(shù)”實現(xiàn)地景仿真[J];中國圖象圖形學報;1999年04期

4 王宇新;李曉杰;閆鴻浩;王小紅;;基于物質(zhì)點法的水流運動三維數(shù)值仿真[J];科學技術與工程;2018年01期

5 鄧海衛(wèi);方漪;齊全;;六點法攝像機標定的研究[J];青島大學學報(自然科學版);2009年01期

6 蘇改明;;試論網(wǎng)格技術對數(shù)字圖書館的影響[J];文教資料;2006年17期

7 葛健平;;計算機地理網(wǎng)格技術在現(xiàn)代農(nóng)業(yè)檢測與控制中的應用[J];電腦知識與技術;2016年28期

8 鄒益民;張智雄;曲云鵬;;基于實體網(wǎng)格的語篇表示模型研究[J];情報理論與實踐;2013年06期

9 張愛芳;;網(wǎng)格技術在數(shù)字圖書館中的應用[J];成功(教育);2011年05期

10 趙秀敏;;網(wǎng)格技術在數(shù)字圖書館建設中的應用及存在的問題[J];知識經(jīng)濟;2011年12期

相關會議論文 前10條

1 潘煒;單艷玲;李鐵瑞;馬騰;江存;高博青;;自由曲面網(wǎng)格劃分的研究現(xiàn)狀[A];第十五屆全國現(xiàn)代結構工程學術研討會論文集[C];2015年

2 孫來福;王朋嬌;;網(wǎng)格在網(wǎng)絡教育資源管理中的應用[A];基于信息技術的教學新模式研究[C];2006年

3 劉云生;張傳富;張童;查亞兵;;基于網(wǎng)格的分布式仿真研究[A];'2003系統(tǒng)仿真技術及其應用學術交流會論文集[C];2003年

4 馮統(tǒng)成;;基于網(wǎng)格計算的遠程教育資源庫設計[A];計算機與教育——全國計算機輔助教育學會第十二屆學術年會論文集[C];2005年

5 符尚武;鄔吉明;戴自換;沈隆鈞;;二維拉格朗日網(wǎng)格的一種積分守恒重映方法[A];中國工程物理研究院科技年報（2002）[C];2002年

6 孟敏;計忠平;劉利剛;;基于保特征調(diào)和場的交互式網(wǎng)格分割[A];中國計算機圖形學進展2008--第七屆中國計算機圖形學大會論文集[C];2008年

7 張燕;賈焰;周斌;;一種基于自適應k近鄰聚類的網(wǎng)格資源發(fā)現(xiàn)方法[A];第二十三屆中國數(shù)據(jù)庫學術會議論文集（技術報告篇）[C];2006年

8 馮欽;羅健欣;倪桂強;;一種新型漸進網(wǎng)格在多模型顯示中的應用[A];第五屆全國幾何設計與計算學術會議論文集[C];2011年

9 王志成;吳頌平;;多塊結構網(wǎng)格并行算法研究[A];北京力學會第20屆學術年會論文集[C];2014年

10 廖毅;程志全;黨崗;;一種基于顯著性分析的網(wǎng)格分割算法[A];中國圖學新進展2007——第一屆中國圖學大會暨第十屆華東六省一市工程圖學學術年會論文集[C];2007年

相關重要報紙文章 前10條

1 realrocking;精彩CG廣告之網(wǎng)格特效揭秘[N];中國電腦教育報;2003年

2 本報記者 鄒大斌;云計算莫走網(wǎng)格的老路[N];計算機世界;2008年

3 龔強;歐洲各國的網(wǎng)格項目[N];中國測繪報;2009年

4 祁金華;網(wǎng)格與虛擬化會融合么？[N];網(wǎng)絡世界;2007年

5 甘麗華;中國挺進網(wǎng)格研究領域[N];經(jīng)濟日報;2002年

6 王元放 湯靜波（作者系上海市互聯(lián)網(wǎng)經(jīng)濟咨詢中心研究員）;網(wǎng)格，讓城市生活更美好[N];文匯報;2003年

7 電腦商報記者 彭敏;應用網(wǎng)格：另一種虛擬化？[N];電腦商報;2009年

8 謝文硯;網(wǎng)格是教育信息化的必由之路[N];中國計算機報;2004年

9 雨青;取消網(wǎng)格限制[N];計算機世界;2004年

10 倪洪章;用網(wǎng)格構建企業(yè)IT系統(tǒng)[N];計算機世界;2006年

相關博士學位論文 前10條

1 于曉燕;基于網(wǎng)格的數(shù)字圖書館模式構建研究[D];南開大學;2010年

2 張花艷;基于變分方法的網(wǎng)格處理[D];中國科學技術大學;2015年

3 王坤;乳腺癌普查醫(yī)學網(wǎng)格研究[D];清華大學;2009年

4 計忠平;輕松交互式網(wǎng)格曲面建模方法的研究[D];浙江大學;2008年

5 杜萍;高質(zhì)量動態(tài)四邊網(wǎng)格序列生成算法研究[D];山東大學;2017年

6 施曉晗;高質(zhì)量大規(guī)模網(wǎng)格形變技術[D];浙江大學;2009年

7 沈建國;面向逆向創(chuàng)新設計的網(wǎng)格特征編輯技術研究[D];浙江大學;2009年

8 鐘鋒;基于網(wǎng)格的煤炭企業(yè)信息化關鍵技術研究[D];中國礦業(yè)大學（北京）;2009年

9 錢江;網(wǎng)格分割算法和相關技術研究[D];浙江大學;2008年

10 胡建偉;保持特征的網(wǎng)格形變和編輯算法的研究[D];浙江大學;2009年

相關碩士學位論文 前10條

1 鄧海衛(wèi);攝像機標定的六點法研究[D];青島大學;2009年

2 孟慶敏;網(wǎng)格技術對軍隊信息化建設的影響研究[D];廈門大學;2008年

3 譚富文;保持稀疏特征的網(wǎng)格拼接[D];浙江大學;2012年

4 王建峰;基于網(wǎng)格的圖像內(nèi)容檢索研究[D];山東理工大學;2009年

5 齊杰;網(wǎng)格技術在動畫渲染中的應用與研究[D];湖南師范大學;2008年

6 鄧元凱;網(wǎng)格幾何編輯處理研究[D];中國計量學院;2014年

7 王偉;地學應用中的網(wǎng)格數(shù)據(jù)可視化[D];吉林大學;2013年

8 謝銘瑤;基于網(wǎng)格技術的跨庫檢索技術研究[D];中南大學;2011年

9 吳瑞豐;基于六邊形網(wǎng)格的幾何處理方法[D];吉林大學;2009年

10 張洋;視覺顯著特征引導的網(wǎng)格有意義分割[D];大連理工大學;2010年

本文編號：2586942