999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

彈性非均勻介質應力場動網格算法改進

2015-04-14 18:20:35顧驍
中國科技縱橫 2015年4期

顧驍

【摘 要】 彈性非均勻介質的彈性應力場模擬是工程力學領域的重要應用。本文以相場微彈性理論為基礎,改進了Chen的彈性非均勻體系的動網格法,提升了30%左右的計算效率。

【關鍵詞】 彈性非均勻 相場微彈性理論 動網格法

1 導論

非均勻介質的彈性應力場研究一直是工程力學領域的一個重要問題。絕大多數的工程材料系統都是彈性非均勻的,即使單相多晶中每一個晶粒都是完美晶體,整個材料也是彈性非均勻的,因為從整體坐標系看,依據不同位相放置的單個晶粒的彈性模量張量也是不同的。而多相系統——比如多層膜、功能梯度材料等常見的人造復合工程材料——不僅是彈性非均勻,更是晶體結構非均勻的。存在晶體缺陷的系統也可以看作一類重要的彈性非均勻系統。人們希望建立一個存在外加應力的三維分析模型來研究上述彈性非均勻系統。但是其分析方法一直受到數學理論的限制[1,2]。

相場微彈性理論是一種彈性非均勻場模擬計算的相場分析方法[3]。它可以計算合金條幅分解[4]、位錯附近的溶質偏析[5]、馬氏體相變[6]等多種動力學過程的彈性應力場演化。相場方法善于模擬復雜結構的晶粒、缺陷等,但往往需要很大的計算量。非均勻計算網絡可以很大程度地降低計算成本,通過網格非均勻化,將計算資源集中在應力變化較大的區域,可以在保證計算結果精度的條件下很好地節約計算資源。Chen等人給出了相場模型計算單元非均勻化算法的修正方程,并將這一方法稱為動網格法[7]。本文在他的算法基礎上,建立了改進的非均勻介質相場微彈性理論的動網格算法。

2 非均勻介質相場微彈性理論的動網格算法及其改進

在計算模型中,計算坐標和物理坐標之間存在映射,當時,計算單元是均勻的。彈性介質的應力應變關系、幾何方程和應力平衡方程都基于物理坐標,將這三個方程結合得到:

(1)

其中,是體系平均彈性模量,是局部彈性模量,是位移矢量,是應變張量,是無應力應變張量,代表由于相變、晶體缺陷等因素引起的塑性應變。設,取弗羅貝尼烏斯范數作為比例因子,保證對于任意矩陣都有,將公式(1)改寫為一個迭代方程:

(2)

其中下標NStep和N+1 Step的變量分別代表第N和N+1迭代步的值。對公式(2)進行傅里葉變換,得到:

(3)

所以,修正的迭代算法通過在倒空間迭代位移矢量得到位移的解,再計算系統應變場。這一迭代算法已經實現了計算網格的非均勻化,為了進一步實現網格在每一次迭代后的重新劃分,首先需要確定局部區域分配的計算格點數目的權重。結合應力場的判定公式,可以以應力場梯度的平方作為權重。其次,需要確定如何將舊網格的數據導入新網格,三次樣條插值等插值方法是很好的選擇。通過合理的計算網格邊長調整方案,非均勻化計算網格的算法可以在應變結果偏差小于1%的條件下將計算規??s小為原來的10到40倍,這一倍數與模型具體的計算內容有關。

在此基礎上,我們得到了改進的非均勻計算單云修正方程。注意到公式(3)的迭代過程中需要計算應變張量。而應變張量需要通過位移矢量求偏導得到。該求偏導過程增加數值求偏導或者額外的傅里葉正逆變換步驟。如果用位移對于物理坐標的梯度作為迭代變量,將公式(3)改寫為:

(4)

那么應變張量可以通過迭代變量的加減運算得到,所以這種替代可以節約計算成本。

3 改進效果分析

改進算法存在計算速度上的優勢,也存在數據存儲空間上的劣勢。在計算速度方面,可以進行定性的比較,假設計算單元一共有N個計算單元,通過傅里葉正逆變換來計算位移對于物理坐標的梯度,從而計算應變張量,并忽略數組傅里葉變換的規模效應。注意到傅里葉正逆變換是迭代的主要計算成本。如果N個標量進行一次傅里葉正變換或逆變換需要t的時間,那么利用公式(3)進行一次迭代的過程中,通過位移計算應變需要,公式(3)中的傅里葉變換需要,總用時為。利用公式(4)進行迭代的時間是,可見新的迭代算法對原算法的加速比為。表1是改進算法和原算法分別計算Cu/Ni納米多層膜中單根位錯產生的應力場的迭代過程的實際計算時間,可見實測的加速比在1.2到1.4之間。所以,替代算法具有加速優勢。在計算存儲空間方面,替代算法中的每個計算單元需要始終儲存位移偏導這個九維數組,而原算法只在部分時間會達到這一存儲量,如果適當調整循環運算順序并犧牲部分計算效率,可以減低存儲量,這對于計算規模相對較大的相場微彈性模型是有意義的。

4 結語

本文在Chen的非均勻彈性場動網格算法的基礎上,重新選擇了迭代變量,并修改了迭代算法,使計算效率提升了30%左右。這一算法改進對于以相場微彈性理論為基礎的晶體缺陷和相變過程的數值模擬具有重要意義,特別是對需要計算規模較大的研究復雜構型的相場模型,可以大大減少計算時間。結合高性能計算和計算網絡非均勻化本身的效果,可以獲得理想的加速比。

參考文獻

[1]Demir I,Hirth J P,Zbib H M.The extended stress field around a cylindrical crack using the theory of dislocation pile-ups[J].International journal of engineering science,1992,30(7):829-845.

[2]Leo P H, Lowengrub J S, Jou H J. A diffuse interface model for microstructural evolution in elastically stressed solids[J].Acta materialia,1998,46(6): 2113-2130.

[3]Wang Y U,Jin Y M,Khachaturyan A G.Phase field microelasticity theory and modeling of elastically and structurally inhomogeneous solid[J].Journal of Applied Physics,2002,92(3):1351-1360.

[4]Wang Y,Banerjee D,Su C C,et al.Field kinetic model and computer simulation of precipitation of L12 ordered intermetallics from fcc solid solution[J].Acta materialia, 1998,46(9):2983-3001.

[5]Hu S Y,Chen L Q.Solute segregation and coherent nucleation and growth near a dislocation—a phase-field model integrating defect and phase microstructures[J].Acta materialia,2001,49(3):463-472.

[6]Jin Y M,Artemev A,Khachaturyan A G. Three-dimensional phase field model of low-symmetry martensitic transformation in polycrystal:simulation of ζ′2 martensite in AuCd alloys[J].Acta materialia,2001,49(12):2309-2320.

[7]Feng W M,Yu P,Hu S Y,et al.A Fourier spectral moving mesh method for the Cahn Hilliard equation with elasticity[J].Commun Comput Phys,2009,5:582-599.

主站蜘蛛池模板: 国产一级特黄aa级特黄裸毛片| 欧美成人午夜视频免看| 欧美成人免费午夜全| 久久人人爽人人爽人人片aV东京热| 国产精品区网红主播在线观看| 亚洲狼网站狼狼鲁亚洲下载| 精品视频免费在线| 国产黄网站在线观看| 精品国产电影久久九九| 91色在线观看| 青青草91视频| 国产成人精品视频一区视频二区| 伊人精品成人久久综合| 一级毛片基地| av天堂最新版在线| 欧美日韩国产成人在线观看| 亚洲日韩AV无码精品| 免费一级毛片在线播放傲雪网| 日本影院一区| 亚洲Av激情网五月天| 国产乱码精品一区二区三区中文| 国产成人亚洲无码淙合青草| 青青青草国产| 尤物亚洲最大AV无码网站| 久996视频精品免费观看| 国产素人在线| 日韩欧美在线观看| 九九热免费在线视频| 国产精品美女免费视频大全 | 欧美日韩午夜视频在线观看 | 国产成人1024精品下载| 91九色最新地址| 99久久99这里只有免费的精品| 91在线无码精品秘九色APP| 国产凹凸视频在线观看| 久久精品aⅴ无码中文字幕| 四虎免费视频网站| 四虎永久在线视频| 毛片久久久| 美女国内精品自产拍在线播放| 99视频精品在线观看| 谁有在线观看日韩亚洲最新视频| 欧美成人日韩| 免费在线播放毛片| 久久夜夜视频| 欧美亚洲欧美| 国产欧美日韩精品第二区| 国产一级做美女做受视频| 欧美日韩专区| 在线观看视频一区二区| 99久久国产自偷自偷免费一区| 日韩精品毛片人妻AV不卡| 欧美色亚洲| 特级精品毛片免费观看| 亚洲系列无码专区偷窥无码| 成人免费视频一区| 国产成人综合久久精品下载| 亚洲日韩精品综合在线一区二区 | 亚洲无码精彩视频在线观看| 亚洲精品高清视频| 女人18一级毛片免费观看| 精品精品国产高清A毛片| 99视频在线免费| 国产91小视频在线观看| 99精品国产自在现线观看| 日韩无码白| 国产91丝袜在线观看| 伊人久久青草青青综合| 亚洲人成日本在线观看| 欧美日韩成人| 亚洲视频色图| 国产97色在线| 全免费a级毛片免费看不卡| 无码一区二区三区视频在线播放| 成人在线观看一区| 欧洲一区二区三区无码| 欧美一区福利| 一本久道久综合久久鬼色| 91成人精品视频| 综合成人国产| 成年女人a毛片免费视频| 欧美三级视频网站|