删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于CUDA的超声二维声场EFIT仿真*

本站小编 Free考研考试/2021-12-25

闂傚倷娴囬褏鈧稈鏅犻、娆撳冀椤撶偟鐛ラ梺鍦劋椤ㄥ懐澹曟繝姘厵闁告挆鍛闂佹娊鏀遍崹鍫曞Φ閸曨垰绠涢柛鎾茬劍閸嬔囨⒑绾懏顥夐柣鈽嗘珋闂傚倸鍊搁崐鐑芥倿閿旈敮鍋撶粭娑樻噽閻瑩鏌熼悜姗嗘畷闁稿孩顨嗛妵鍕棘閸喒鎸冮梺鍛婎殕瀹€鎼佸蓟濞戙垹绠绘俊銈傚亾闁硅櫕鍔欓幃锛勨偓锝庡枟閸婂灚顨ラ悙鑼虎闁告梹宀搁弻锝呂旈埀顒勫疮閺夋垹鏆﹂柟杈剧畱缁犲鎮归崶顏勭毢妞わ负鍔岄—鍐Χ閸涱垳顔囧┑鈽嗗亝椤ㄥ﹪骞冮垾鏂ユ闁靛繆妾ч幏娲⒑缁嬭法绠伴柛姘儑濞嗐垽寮婚妷锔惧幈闂佺粯姊绘灙闁崇粯娲熼弻鐔碱敍濞嗘垵濡介悗鍨緲鐎氼喗绂掗敃鍌氱畾鐟滃酣鎮楅銏♀拻濞达綀娅i妴濠囨煕閹惧绠栫紒妤冨枎铻栭柛娑卞幗鏉堝牓姊洪棃娑辨Т闁哄懏绮撻崺娑㈠箳濡や胶鍘卞┑鐐村灥瀹曨剟鐛Ο鑲╃<闁绘瑥鎳愰幃鑲╃磼鏉堛劌绗уù鐙呯畵閹墽浠﹂悾灞藉箺闂傚倷娴囬鏍垂閺夋5娲Ω閳轰胶鍔﹀銈嗗笒閿曪箓宕氶弶娆剧唵鐟滃秴霉妞嬪骸鍨濋柛顐ゅ枎缁剁偤鎮楅敐搴′簻鐎殿喖娼″娲箰鎼达絿鐣甸梺纭呮珪閸旀鍒掓繝姘煑濠㈣泛鐬奸惁鍫ユ⒑闁偛鑻晶鎾煙椤斻劌娲ら柋鍥煃閸ㄦ稒娅嗛柣蹇撳暣閺岋絾鎯旈敍鍕殯闂佺ǹ閰f禍鍫曠嵁婵犲啯鍎熼柕濠忕畱娴滄姊虹粙鎸庢拱闁糕晛鍟嵄闁割偁鍎查悡娆戠磽娴i潧鐏╅柣蹇曞Т闇夋繝濠傜墢閻g儤鎱ㄦ繝鍐┿仢鐎规洏鍔嶇换婵嬪磼濠婂懏鍣┑鐘殿暯濡插懘宕戦崨顖涘床闁割偁鍎冲畵渚€鏌熼幍顔碱暭闁哄懏鎮傞弻銊╂偆閸屾稑顏�
547闂傚倸鍊搁崐椋庣矆娴i潻鑰块梺顒€绉查埀顒€鍊圭粋鎺斺偓锝庝簽閿涙盯姊洪悷鏉库挃缂侇噮鍨堕崺娑㈠箳濡や胶鍘遍梺鍝勬处椤ㄥ棗鈻嶉崨瀛樼厽闊浄绲奸柇顖炴煛瀹€瀣埌閾绘牠鏌嶈閸撶喖寮绘繝鍥ㄦ櫜濠㈣泛锕﹂悿鍥⒑鐟欏嫬绀冩い鏇嗗懐鐭嗛柛鎰ㄦ杺娴滄粓鐓崶銊﹀鞍妞ゃ儲绮撻弻锝夊箻鐎靛憡鍒涘┑顔硷攻濡炶棄鐣峰Δ鍛闁兼祴鏅涢崵鎺楁⒒娴e憡鎲搁柛锝冨劦瀹曟垿宕熼娑樹患闂佺粯鍨兼慨銈夊疾閹间焦鐓ラ柣鏇炲€圭€氾拷1130缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎹囬弻锝夊閵忊晜姣岄梺绋款儐閹瑰洤鐣疯ぐ鎺濇晝闁挎繂娲﹂濠氭⒒娓氣偓閳ь剛鍋涢懟顖涙櫠閸欏浜滄い鎰╁焺濡叉椽鏌涢悩璇у伐妞ゆ挸鍚嬪鍕節閸愵厾鍙戦梻鍌欑窔閳ь剛鍋涢懟顖涙櫠閹绢喗鐓涢悘鐐登规晶鑼偓鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊勭矒閿濈偞鎯旈埦鈧弨浠嬫煟閹邦垰鐨哄褎鐩弻娑㈠Ω閵壯傝檸闂佷紮绲块崗姗€寮幘缁樺亹闁肩⒈鍓﹀Σ浼存煟閻斿摜鐭婄紒缁樺笧閸掓帒鈻庨幘宕囧€為梺鍐叉惈閸熶即鏁嶅⿰鍕瘈闁靛骏绲剧涵楣冩煥閺囶亪妾柡鍛劦濮婄粯鎷呴崨濠傛殘闁煎灕鍥ㄧ厱濠电姴鍟版晶杈╃磽閸屾稒宕岄柟绋匡攻缁旂喖鍩¢崒娑辨閻庤娲︽禍婵嬪箯閸涱垱鍠嗛柛鏇ㄥ幗琚欓梻鍌氬€风粈浣革耿闁秴鍌ㄧ憸鏃堝箖濞差亜惟闁宠桨鑳堕鍥⒑閸撴彃浜濇繛鍙夌墵閹偤宕归鐘辩盎闂佺懓顕崑娑㈩敋濠婂懐纾煎ù锝呮惈椤eジ鏌曢崶褍顏い銏℃礋婵偓闁宠桨绀佹竟澶愭⒒娴g懓顕滅紒瀣浮瀹曟繂鈻庨幘璺虹ウ闁诲函缍嗛崳顕€寮鍡欑瘈濠电姴鍊规刊鍏间繆閺屻儲鏁辩紒缁樼箞閹粙妫冨☉妤佸媰闂備焦鎮堕崝宀€绱炴繝鍌ゅ殨妞ゆ劑鍊楅惌娆愪繆椤愩倖鏆╅柛搴涘€楅幑銏犫攽鐎n亞鍊為梺闈浨归崕鏌ヮ敇濞差亝鈷戦柛婵嗗濡叉悂鏌eΔ浣虹煉鐎规洘鍨块獮鎺懳旈埀顒勫触瑜版帗鐓涢柛鎰╁妿婢ф盯鏌i幘宕囩闁哄本鐩崺鍕礃閳哄喚妲烽梻浣呵圭换鎰版儔閼测晜顫曢柟鐑橆殢閺佸﹪鏌涜箛鎿冩Ц濞存粓绠栭幃娲箳瀹ュ棛銈板銈庡亜椤︾敻鐛崱娑樻閹煎瓨鎸婚~宥夋⒑閸︻厼鍔嬮柛銊ㄦ珪缁旂喖寮撮悢铏诡啎闁哄鐗嗘晶浠嬪箖婵傚憡鐓涢柛婊€绀佹禍婊堝础闁秵鐓曟い鎰Т閸旀粓鏌i幘瀛樼闁哄瞼鍠栭幃婊兾熺拠鏌ョ€洪梻浣呵归鍥ㄧ箾閳ь剟鏌$仦鐣屝ユい褌绶氶弻娑滅疀閺冨倶鈧帗绻涢崱鎰仼妞ゎ偅绻勯幑鍕洪鍜冪船婵犲痉鏉库偓褏寰婃禒瀣柈妞ゆ牜鍋涚粻鐘虫叏濡顣抽柛瀣崌閻涱噣宕归鐓庮潛闂備礁鎽滈崰鎾寸箾閳ь剛鈧娲橀崹鍧楃嵁濡皷鍋撳☉娅亪顢撻幘缁樷拺缂備焦锚閻忥箓鏌ㄥ鑸电厓鐟滄粓宕滃☉銏犵;闁绘梻鍘ч悞鍨亜閹烘垵鏋ゆ繛鍏煎姍閺岀喖顢欓懖鈺佺厽閻庤娲樺ú鐔笺€佸☉銏″€烽柤纰卞墮婵附淇婇悙顏勨偓鏍垂婵傜ǹ纾垮┑鐘宠壘缁€鍌炴倶閻愭澘瀚庡ù婊勭矒閺岀喖骞嗚閹界娀鏌涙繝鍐ㄥ闁哄瞼鍠栭、娆撴嚃閳轰胶鍘介柣搴ゎ潐濞茬喐绂嶉崼鏇犲祦闁搞儺鍓欐儫闂侀潧顦崐鏇⑺夊顑芥斀闁绘劘鍩栬ぐ褏绱掗懠顒€浜剧紒鍌氱Ч閹崇偤濡疯濞村嫰姊洪幐搴㈢5闁稿鎹囧Λ浣瑰緞閹邦厾鍘遍棅顐㈡处濞叉牜鏁崼鏇熺厵闁稿繐鍚嬮崐鎰版煛鐏炵晫啸妞ぱ傜窔閺屾稖绠涢弮鍌楁闂傚洤顦甸弻娑㈠Ψ椤旂厧顫╃紒鐐劤閵堟悂寮婚弴鐔虹瘈闊洦娲滈弳鐘差渻閵堝棙绀夊瀛樻倐楠炲牓濡搁妷搴e枔缁瑩宕归纰辨綍闂傚倷鑳舵灙妞ゆ垵妫濋獮鎰節濮橆剛顔嗛梺鍛婁緱閸ㄩ亶宕伴崱娑欑厱闁哄洢鍔屾晶浼存煛閸℃ê鍝烘慨濠勭帛閹峰懘宕崟顐$帛闁诲孩顔栭崰妤呭磿婵傜ǹ桅闁圭増婢樼粈鍐┿亜韫囨挻顥犲璺哄娣囧﹪濡惰箛鏇炲煂闂佸摜鍣ラ崹璺虹暦閹达附鍋愮紓浣贯缚閸橀亶姊洪弬銉︽珔闁哥噥鍋呴幈銊╁焵椤掑嫭鈷戠紒瀣儥閸庢劙鏌熺粙娆剧吋妤犵偛绻樺畷銊р偓娑櫭禒鎯ь渻閵堝棛澧柤鐟板⒔缁骞嬮敂瑙f嫽婵炶揪绲介幉锟犲箚閸儲鐓曞┑鐘插閸︻厼寮查梻渚€娼х换鍫ュ磹閺囥垺鍊块柛顭戝亖娴滄粓鏌熺€电ǹ浠滄い鏇熺矌缁辨帗鎷呯憴鍕嚒濡炪値鍙€濞夋洟骞夐幘顔肩妞ゆ巻鍋撶痪鐐▕閹鈻撻崹顔界亾闂佽桨绀侀…鐑藉Υ娴g硶妲堟俊顖涚矌閸犲酣鎮鹃埄鍐跨矗濞达絽澹婂Λ婊勭節閻㈤潧浠╅柟娲讳簽缁辩偤鍩€椤掍降浜滄い鎰╁焺濡偓闂佽鍣换婵嬪春閳ь剚銇勯幒鎴濐仾闁抽攱甯¢弻娑氫沪閹规劕顥濋梺閫炲苯鍘哥紒顔界懇閵嗕礁鈻庨幇顔剧槇闂佸憡娲﹂崜锕€岣块悢鍏尖拺闁告挻褰冩禍婵囩箾閸欏澧辩紒顔垮吹缁辨帒螣闂€鎰泿闂備礁婀遍崑鎾翅缚濞嗘拲澶婎潩閼哥數鍘遍柣搴秵閸嬪懐浜告导瀛樼厵鐎瑰嫮澧楅崵鍥┾偓瑙勬礈閸忔﹢銆佸Ο琛℃敠闁诡垎鍌氼棜濠电姷鏁告慨鏉懨洪敃鍌氱9闁割煈鍋嗙粻楣冩煙鐎涙ḿ绠橀柡瀣暟缁辨帡鍩€椤掑倵鍋撻敐搴℃灍闁绘挸鍟伴幉绋库堪閸繄顦у┑鐐村灦濮樸劑鎯岄崱妞曞綊鏁愰崼鐔粹偓鍐煟閹烘埊韬柡宀€鍠庨埢鎾诲垂椤旂晫浜愰梻浣呵归鍡涘箰閹间礁鐓″璺哄閸嬫捇宕烽鐐愩儲銇勯敂鍨祮婵﹥妞介弻鍛存倷閼艰泛顏梺鍛娒幉锛勬崲濞戙垹绾ч柟瀵稿仜閺嬬姴顪冮妶鍐ㄧ仾闁挎洏鍨归悾鐑筋敃閿曗偓鍞悷婊冪灱缁厽寰勬繛鐐杸闁圭儤濞婂畷鎰板箻缂佹ê鈧潡鏌ㄩ弮鈧畷妯绘叏閾忣偅鍙忔俊顖氱仢閻撴劙鏌i幘宕囩闁哄本鐩崺鍕礃閳哄喚妲舵俊鐐€х拋锝嗕繆閸ヮ剙鐒垫い鎺嗗亾婵犫偓鏉堛劎浠氭俊鐐€ら崢濂稿床閺屻儲鍋╅柣鎴eГ閺呮煡鏌涢妷顖炴闁告洖鍟村铏圭矙閹稿孩鎷卞銈冨妼閹冲繒绮嬪澶婄畾妞ゎ兘鈧磭绉洪柡浣瑰姍瀹曘劑顢欓崗鍏肩暭闂傚倷绀侀幉鈥趁洪悢铏逛笉闁哄稁鍘奸拑鐔兼煥濠靛棭妲归柛濠勫厴閺屾稑鈻庤箛锝嗏枔濠碘槅鍋呴崹鍨潖濞差亝鐒婚柣鎰蔼鐎氫即鏌涘Ο缁樺€愰柡宀嬬秮楠炴帡鎮欓悽鍨闁诲孩顔栭崳顕€宕滈悢椋庢殾闁圭儤鍩堝ḿ鈺呮煥濠靛棙顥犻柛娆忓暞缁绘繂鈻撻崹顔界亾闂佺娅曢幐鍝ュ弲闂佺粯枪椤曆呭婵犳碍鐓欓柟顖嗗懏鎲兼繝娈垮灡閹告娊寮诲☉妯锋婵鐗婇弫楣冩⒑闂堚晝绋婚柟顔煎€垮濠氭晲閸℃ê鍔呴梺闈涚箳婵挳寮稿▎鎾寸厽闁绘ê鍟挎慨澶愭煕閻樺磭澧电€规洘妞介崺鈧い鎺嶉檷娴滄粓鏌熺€电ǹ浠滄い鏇熺矋閵囧嫰鏁冮崒銈嗩棖缂備浇椴搁幐鎼侇敇婵傜ǹ妞藉ù锝嚽规竟搴ㄦ⒒娴d警鏀版繛鍛礋閹囨偐鐠囪尙鐤囬梺缁樕戝鍧楀极閸℃稒鐓曢柟閭﹀枛娴滈箖鏌﹂幋婵愭Ш缂佽鲸鎹囧畷鎺戔枎閹存繂顬夐梻浣告啞閸旀洟鈥﹂悜鐣屽祦闊洦绋掗弲鎼佹煥閻曞倹瀚�28缂傚倸鍊搁崐鎼佸磹妞嬪孩顐介柨鐔哄Т绾捐顭块懜闈涘Е闁轰礁顑囬幉鎼佸籍閸稈鍋撴担鑲濇棃宕ㄩ闂寸盎闂備焦鍎崇换鎰耿闁秵鍋傞悗锝庡枟閳锋垿鎮峰▎蹇擃仾闁稿孩顨婇弻娑㈠Ω閵壯嶇礊缂備緡鍠栭悧濠勬崲閸涱厸鏋庣紒顖涚缁绘繈鎮介棃娑楁勃濠电偛鍚嬮悷鈺佺暦閵壯€鍋撻敐搴℃灍闁搞倕鑻灃闁挎繂鎳庨弳娆戠磼閻橀潧鏋涢柡灞熷棛鐤€闁挎繂鎳嶇花鍏肩節濞堝灝鐏¢柟鍛婂▕閻涱噣寮介‖銉ラ叄椤㈡鍩€椤掑嫬鐒垫い鎺嶇劍缁€鈧梺宕囩帛閹瑰洭寮幘缁樻櫢闁跨噦鎷�128.00闂傚倸鍊搁崐鐑芥嚄閼哥數浠氱紓鍌欒兌缁垶宕濆▎蹇曟殾闁哄绨卞Σ鍫熸叏濮楀棗寮炬繛鑲╁枛濮婅櫣鎹勯妸銉︾彚闂佺懓鍤栭幏锟�
在超声工业检测中,一般采用数值仿真技术对超声换能器在不同材料中的声场进行模拟计算[1]。其中,弹性动力学有限积分算法(Elastodynamic Finite Integration Thecnology,EFIT)采用交叉网格模型,对均匀和非均匀介质中声场的模拟具有较好的计算稳定性[2]。但随着仿真对象几何尺寸的增大,计算单元数量迅速增长,导致数值计算耗时过长,时间成本大幅增加,严重限制了EFIT在实际工程中的应用。
通用图形处理器(Graphic Processing Unit, GPU)能够提供高速的并行计算能力,而计算设备统一构架(Compute Unified Device Architecture, CUDA)使基于NVIDIA公司GPU产品的并行计算能够方便地推广到一般工程应用。NVIDIA公司GPU和CUDA相结合,为数值计算效率的提升提供了新的实现手段,在成像算法、物理场模拟等方面得到了广泛的应用[3-6]
本文基于CUDA将并行计算技术应用于超声声场的EFIT,对钢材料中的声场分布进行仿真计算,测试基于GPU并行计算和传统CPU计算的效率差异。
1 弹性动力学有限积分算法 1.1 算法的积分形式 EFIT由Fellinger和Langenberg[7]在1989年提出。该算法基于牛顿动量定理和广义胡克定律给出了弹性动力学运动方程的积分形式,并对方程组进行离散化,基于数值计算的方式在时间和空间上进行迭代求解。通过研究表明,EFIT不但能够模拟计算各向同性固体介质中弹性波的传播和散射情况,而且对于各向异性固体材料也具有较好的算法收敛性[8-10]
假设一个封闭的固体单元,体积以V表示,封闭表面的面积以S表示, 则根据牛顿动量定理可得
(1)

根据动量密度矢量定义:
(2)

将式(2)代入式(1)可得
(3)

根据变形定律可得[2, 9]
(4)

式中:R为质点的空间位置;t为时间;n为曲面外法向向量;p(R, t)为质点的动量密度矢量;σ(R, t)为应力张量;ε(R, t)为应变张量;v(R, t)为质点速度矢量;ρ(R)为密度。
根据广义胡克定律可知,在弹性假设前提下,应力与应变之间为单值线性关系,从而建立了式(3)与式(4)中应力与应变之间的联系。
(5)

式中:s为顺度系数矩阵。
1.2 积分形式的离散化 在获得弹性动力学运动方程积分表达式的基础上,为进行数值计算需要对相关物理量进行离散化。
假设整个仿真计算区域由紧密排布的相邻物质单元构成,每个物质单元均为直角坐标系下边长相等的立方体单元网格,同一个物质单元网格内各相关物理量是恒定的。按照EFIT,在整个计算区域内需要将离散应力点和离散速度点交叉配置。在一个网格单元中,3个正应力分量处于网格单元的中心,3个切应力分量分别处于3个边的中点,3个速度分量则分别处于物质单元的3个面的中心位置,从而构成速度分量与应力分量相互交错的网格放置方式(见图 1)。
图 1 单个计算单元内速度与应力的分布 Fig. 1 Velocity and stress distribution in single calculation unit
图选项




对于各向同性材料,弹性常数矩阵中独立常数个数减少为2个,各向同性材料弹性常数矩阵形式如下:

式中:λμ为拉梅常数的2个参数。
因此,当在二维条件下计算各向同性材料中超声波的传播和散射时,弹性动力学运动方程离散后可以简化为

式中:n为迭代步数;MxMz分别表示xy方向计算单元的步进量。
在此基础上,在时间序列k上进行迭代计算,从而计算得到各时刻超声在固体介质中的传播和散射情况。

1.3 吸收边界 在实际工程应用中,有时检测对象体积较大,如果对完整对象建模计算,将需要巨大的网格单元数量,由于硬件资源的限制,实际很难实现,因此,在实际工程中经常只对有限区域的声场进行计算。
但当以有限区域模拟较大区域时,由于存在计算区域边界,将引入额外的计算边界反射。为了避免模型中由于计算区域的限制形成边界反射回波,本文采用完全匹配层(Perfect Matched Layer,PML)吸收边界条件[11-14],在计算区域边界处设置了一定数量的吸收层,达到较好的边界吸收效果。图 2为未添加吸收边界时的声场分布情况,图 3为添加吸收边界后的声场分布情况。
图 2 未设置吸收边界的声场分布 Fig. 2 Sound field distribution without absorption boundary
图选项




图 3 设置吸收边界的声场分布 Fig. 3 Sound field distribution with absorption boundary
图选项




2 基于CUDA的EFIT模型优化 2.1 GPU与CUDA 与传统CPU流水线型的数据处理方式不同,为了适应大规模并行计算的需要,GPU采用并行多线程的内部构架,有效提升并行计算任务的执行效率。
2016年,NVIDIA公司发布了第一款基于CUDA的通用图形处理器(GPGPU),为GPU广泛应用于通用计算领域奠定了基础。
在提供高速并行计算能力的同时,GPU也存在自身的技术限制。为了达到高速的并行计算能力,GPU减少了逻辑计算单元的数量,因此逻辑判断能力相对于CPU存在差距;同时,只有相同的计算操作才能并行化,不同外部条件导致的算法变更可能会由GPU自行转为串行处理,过多的逻辑判断将极大地降低计算效率。在极端情况下,不合理的程序设计会导致GPU计算速度远低于CPU。
因此,在使用GPU进行并行计算时,首先需要对计算模型进行分析,判断是否具备并行计算的条件; 然后对程序逻辑进行优化,合理分配GPU与CPU之间的工作内容,从而基于CUDA实现高质量的并行计算[15]
2.2 EFIT模型的并行实现 根据EFIT,当前速度场每个单元的变化量由前一时间节点应力场相关单元的数值计算得出,而当前应力场变化量则由计算后的速度场相关单元的数值计算得出。可见,在EFIT的计算模型中,待计算场量仅与已获得的计算结果相关,而待计算场量各个单元之间在当前计算步骤中没有关联。因此,并行计算方式可以简单、高效地实现EFIT,而无需考虑计算过程中不同计算结果间的同步问题。
基于以上分析,建立EFIT并行计算模型,算法实现流程如图 4所示。由于GPU主要进行并行计算,GPU中存储单元的分配仍通过CPU进行控制。在EFIT计算模型初始化时,分别为vxvzσxxσzzσxz在GPU内存中分配二维的存储空间。在计算过程中,为实现仿真计算中间变量的灵活保存,通过CPU控制GPU分别完成速度场和应力场的分步计算。
图 4 并行计算流程 Fig. 4 Flowchart of parallel calculation
图选项




在程序的具体实现中,分别针对vxvzσxxσzzσxz设计各自的核函数。对于空间各单元,首先将二维矩阵转换为一维向量,然后通过idx=blockIdx.x*blockDim.x+threadIdx.x方式对所有空间计算单元进行遍历;对于时间迭代则通过CPU控制程序循环实现。
2.3 纹理内存 在传统图像处理GPU中,纹理内存是专门针对内存访问模式中存在大量空间局部性的图形应用程序所设计的,对线程读取相邻近的内存位置进行了优化。实现方式是通过将全局内存绑定为纹理内存,对其读写操作通过专门的纹理缓存(texture cache)进行。而将纹理内存应用于通用计算,则可以利用纹理缓存针对二维空间的局部性访问优化,实现二维矩阵邻域访问速度的有效提升。
在EFIT计算模型中,每个单元的值都是通过另一个二维矩阵中相邻单元的值计算获得,如图 5所示。由于地址非连续,在采用一般内存进行计算时,这些地址没有被缓存,因此每次访问时需要重新进行寻址,导致内存访问速度较低。而采用纹理内存时,相邻单元的地址将被整体缓存,极大降低访存时间,从而使程序性能得到有效提升[16]
图 5 纹理内存线程寻址示意图 Fig. 5 Schematic diagram of texture memory thread addressing
图选项




本文计算模型中,将速度场与应力场相关矩阵均与纹理内存进行绑定,从而实现模型计算速度的优化。由于GPU中纹理内存容量有限,在设计实际算法时需要考虑纹理内存容量限制,当计算区域超出纹理内存容量时,需要进行分块计算。
2.4 吸收边界优化 为对有限区域进行模拟计算,在EFIT模型边界设置一定数量的吸收层,形成吸收区域。与模型中一般的计算单元不同,每个吸收层需要采用包含衰减因子的波动方程,通过多个吸收层的衰减近似达到边界声场吸收的效果。
引入衰减系数,导致吸收区域波动方程形式变化。模型中可以按照单元位置判断是否处于吸收区域,从而决定是否增加衰减系数。但这种实现方式可能导致GPU在同一线程束(wrap)内的分支语句执行不同的指令路径,此时条件判断的不同指令路径线程只能被顺序执行,程序的并行性被破坏,导致计算速度降低[17]
因此,在基于CPU的模型初始化过程中,增加了吸收系数矩阵的定义,对于一般计算区域,吸收系数预置为1.0,而对于吸收区域单元则设置为小于1.0的衰减系数。在基于GPU进行计算时,只需调用相应的衰减系数矩阵元素即可,从而保证模型计算的并行性。
2.5 数据传输优化 GPU具有高速的并行计算能力,但CPU与GPU之间的数据传输却受限于传输的带宽,成为GPU并行计算的瓶颈。因此在进行模型并行计算过程中,需要合理设计CPU和GPU之间的数据传输,使数据传输量与传输带宽相匹配。
对于EFIT计算模型,由于计算过程为迭代计算,当前计算参数为前一时间节点的计算结果,因此在算法实现时,只需在程序初始化时完成相关参数和参数矩阵的初始化,在并行计算过程中,无需外部数据传入GPU,避免数据传输带来的速度降低。
如果每次都将计算得到声场分布结果输出显示,数据由GPU输出到CPU的时间将成为主要的时间消耗,同时在实际应用中过于密集的声场成像也没有意义。因此,在算法实现中,应根据实际需要合理设置计算结果的输出频率,从而平衡实际应用需求和数据传输效率。
3 仿真结果与效率分析 为了定量测量并行计算对于EFIT的加速效果,本文建立二维点源激励声场模型。模型中,检测对象为钢材料,弹性模量为207 GPa,泊松比为0.289,材料密度为7.9×103 kg/m3,纵波声速设置为5 900 m/s。计算区域选取单元数为512×512,单个单元边长设置为纵波波长的1/10,时间步长为换能器周期的1/20。激励源选用中心频率为100 kHz的点源,振动5个周期,振动方向为x方向,激励源施加于计算区域的中心点。选取计算区域边缘64层单元作为吸收边界。各单元系数采用单精度浮点数进行存储。
在并行计算时,选取GPU型号为NIVIDIA NVS 5400M,CUDA核心数为96。CUDA运算在Visual Studio 2010环境下基于C语言编程,计算结果通过MATLAB界面进行显示。
进行CPU和GPU计算结果对比之前,首先确定GPU计算时选用的核函数参数。分别设置GPU线程块内线程数为16×16和32×32时进行对比,选取最优的GPU计算参数。在进行200次迭代计算时,线程格内线程数为16×16的条件下,并行计算时间为1 060.826 ms,而采用线程格内线程数为32×32的条件下,并行计算时间为1 088.902 ms。可见,两者耗时基本相同,线程格内线程数为16×16时计算耗时略短,因此,后续计算中均采用线程格内线程数16×16的参数设置。
分别通过CPU和GPU运行程序进行计算,选取第450次迭代和第650次迭代时主应力的计算结果进行观察,显示图像如图 6所示。
图 6 GPU与CPU计算结果对比 Fig. 6 Comparison of GPU and CPU computing results
图选项




通过对比可以看出,基于CPU和GPU计算结果完全一致。在此基础上,分别统计不同迭代次数时计算耗时。
在统计计算时间时,无论是CPU程序或GPU程序,仅统计计算本身耗时,程序开始运行时的初始化和计算结果导出、显示时间等部分不计入统计时间。统计结果如图 7图 8所示。
图 7 GPU与CPU计算耗时对比 Fig. 7 Comparison of computing time consumption of GPU and CPU
图选项




图 8 GPU与CPU平均每秒迭代次数对比 Fig. 8 Comparison of average iteration times in one second of GPU and CPU
图选项




通过图 7图 8统计结果可以看出,采用GPU并行计算对EFIT模型进行迭代计算时,相对于CPU计算速度有明显的提升。本文模型中,对于不同的迭代次数,计算速度提升在23~31倍之间。
4 结论 本文基于CUDA编程,实现了EFIT模型的并行计算,并在此基础上对于EFIT模型的GPU和CPU计算速度进行了对比,可以得到如下结论:
1) EFIT模型中,当前单元的计算只与对应位置的已知计算结果有关,具有良好的并行计算条件。
2) 基于CUDA框架能够完成EFIT模型的并行计算程序开发,计算结果与CPU程序计算结果一致。
3) EFIT模型采用并行计算方法进行计算,能够有效提升计算速度,而随着GPU性能的快速提升,并行计算方法对于复杂声场仿真应用具有广阔的应用前景。
由于本文中模型主要用于测试,需要获取部分中间变量,如在实际应用中主要关注最终计算结果,EFIT模型的并行计算仍可进行进一步优化。

参考文献
[1] 徐娜, 李洋, 周正干, 等. FDTD方法的改进及在超声波声场计算中的应用[J]. 北京航空航天大学学报, 2013, 39(1): 78-82.
XU N, LI Y, ZHOU Z G, et al. Improvement of finite difference time domain method and its application to calculation of ultrasonic sound fields[J]. Journal of Beijing University of Aeronautics and Astronautics, 2013, 39(1): 78-82. (in Chinese)
[2] FELLINGER P, MARKLEIN R, LANGENBERG K J, et al. Numerical modeling of elastic wave propagation and scattering with EFIT-Elastodynamic finite integration technique[J]. Wave Motion, 1995, 21: 47-66. DOI:10.1016/0165-2125(94)00040-C
[3] 张霞, 何兴无. CUDA平台下的超声弹性成像并行处理算法[J]. 计算机与数字工程, 2012, 40(9): 113-116.
ZHANG X, HE X W. A Parallel algorithm of ultrasound strainimaging based on CUDA[J]. Computer & Digital Engineering, 2012, 40(9): 113-116. DOI:10.3969/j.issn.1672-9722.2012.09.038 (in Chinese)
[4] 贾春刚, 郭立新, 刘伟. 基于GPU的并行FDTD方法在二维粗糙面散射中的应用[J]. 电波科学学报, 2016, 31(4): 683-687.
JIA C G, GUO L X, LIU W. GPU-based FDTD method for analysis of electromagnetic scattering from a 2D rough surface[J]. Chinese Journal of Radio Science, 2016, 31(4): 683-687. (in Chinese)
[5] 付小波, 马中高, 余嘉顺, 等. 基于多图形处理单元加速的各向异性弹性波正演模拟[J]. 科学技术与工程, 2018, 18(11): 16-22.
FU X B, MA Z G, YU J S, et al. Anisotropic elastic wave forward modeling based on multiple graphics processing unit[J]. Science Technology and Engineering, 2018, 18(11): 16-22. DOI:10.3969/j.issn.1671-1815.2018.11.002 (in Chinese)
[6] 杨尚琴. 地震正演数值模拟仿真计算的并行优化设计方法[J]. 地球物理学进展, 2017, 32(3): 1290-1296.
YANG S Q. Parallel optimization design method for seismic forward modeling numerical simulation calculation[J]. Progress in Geophysics, 2017, 32(3): 1290-1296. (in Chinese)
[7] FELLINGER F, LANGENBERG K J.Numerical techniques for elastic wave propagation and scattering[C]//Proceedings of the IUTAM Symposium on Elastic Wave Propagation and Ultrasonic Evaluation, 1990: 81-86.
[8] SCHUBERT F. Numerical time-domain modeling of linear and nonlinear ultrasonic wave propagation using finite integration techniques-Theory and applications[J]. Ultrasonics, 2004, 42(1-9): 221-229. DOI:10.1016/j.ultras.2004.01.013
[9] 丁辉. 计算超声学——声场分析及应用[M]. 北京: 科学出版社, 2010: 33-36.
DING H. Computational ultrasonics-Analysis and application of ultrasonic fiels[M]. Beijing: Science Press, 2010: 33-36. (in Chinese)
[10] 余涛.超声波在混凝土中传播的数值模拟[D].长沙: 中南大学, 2013: 4-9.
YU T.Numerical simulation of ultrasonic wave propagation in concrete[D].Changsha: Central South University, 2013: 4-9(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2426440
[11] BERENGER J. A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 1994, 114(2): 185-200. DOI:10.1006/jcph.1994.1159
[12] 廉西猛, 单联瑜, 隋志强, 等. 地震正演数值模拟完全匹配层吸收边界条件研究综述[J]. 地球物理学进展, 2015, 30(4): 1725-1733.
LIAN X M, SHAN L Y, SUI Z Q, et al. An overview of research on perfectly matched layers absorbing boundary condition of seismic forward numerical simulation[J]. Progress in Geophysics, 2015, 30(4): 1725-1733. (in Chinese)
[13] 刘洋. 波动方程时空域有限差分数值解及吸收边界条件研究进展[J]. 石油地球物理勘探, 2014, 49(1): 35-46.
LIU Y. The review of finite difference numerical solution for wave equation in time domain and obsorption boundary conditions[J]. Oil Geophysical Prospecting, 2014, 49(1): 35-46. (in Chinese)
[14] 秦臻, 任培罡, 姚姚, 等. 弹性波正演模拟中PML吸收边界条件的改进[J]. 地球科学——中国地质大学学报, 2009, 34(4): 658-664.
QIN Z, REN P G, YAO Y, et al. Improvement of PML absorbing boundary conditions in elastic wave forward modeling[J]. Earth Science-Journal of China University of Geosciences, 2009, 34(4): 658-664. (in Chinese)
[15] 卢风顺, 宋君强, 银福康, 等. CPU/GPU协同并行计算研究综述[J]. 计算机科学, 2011, 38(3): 5-9.
LU F S, SONG J Q, YIN F K, et al. Survey of CPU/GPU synergetic parallel computing[J]. Computer Science, 2011, 38(3): 5-9. DOI:10.3969/j.issn.1002-137X.2011.03.002 (in Chinese)
[16] SANDERS J, KANDROT E.GPU高性能编程CUDA实战[M].聂雪军, 等, 译.北京: 机械工业出版社, 2011: 84-100.
SANDERS J, KANDROT E.CUDA by example——An introduction to general-purpose GPU programming[M].NIE X J, et al., translated.Beijing: China Machine Press, 2011: 84-100(in Chinese).
[17] 方民权, 张卫民, 方建滨, 等. GPU编程与优化:大众高性能计算[M]. 北京: 清华大学出版社, 2016: 273-276.
FANG M Q, ZHANG W M, FANG J B, et al. GPU programming and code optimization:High performance computing for the masses[M]. Beijing: Tsinghua University Press, 2016: 273-276. (in Chinese)


闂傚倸鍊峰ù鍥х暦閻㈢ǹ绐楅柟閭﹀劒濞差亜绠i柨鏃囨閳ь剙娼¢弻銊╁籍閸屾稒鐝梺缁樼箚閸╂牠骞堥妸銉富閻犲洩寮撴竟鏇㈡⒒娴d警鏀版繛鍛礋楠炴垿宕惰閺嗭附鎱ㄥ璇蹭壕濡炪們鍨哄Λ鍐ㄧ暦閸楃倣鏃堝礃閵娿儳顔囬梻鍌氬€搁崐椋庣矆娓氣偓瀹曘儳鈧綆浜堕悢鍡樻叏濡炶浜鹃悗瑙勬处閸ㄥ爼銆侀弴銏狀潊闁靛繆鍓濋鏇㈡⒒娓氣偓濞佳嚶ㄩ埀顒€鈹戦垾铏枠闁糕斂鍨藉顕€宕奸悢鍝勫箥闂備胶顢婇~澶愬礉閺囥垺鍎嶆繛宸簼閻撶喖鏌i弮鍫缂佸妞介弻鐔碱敊閼测晝楔閻庤娲橀崕濂杆囬弻銉︾厓闂佸灝顑嗛埛鎰磼缂佹ḿ绠為柛鈹惧亾濡炪倖甯掔€氼剟鎮炲ú顏呯厱闁归偊鍨辩涵鍫曟煟椤撶噥娈旈柍瑙勫灴閸┿儵宕卞Δ鍐у寲缂備胶铏庨崢濂告偉閻撳海鏆﹂柟杈剧畱缁犲鎮楅棃娑欐喐闁硅櫕绻堝娲焻閻愯尪瀚板褎鎸抽弻鏇㈠幢濡も偓閳ь剙娼¢獮濠囨倷閸濆嫀銊╂煥閺囨浜鹃梺鍛婃煥閹虫ê顫忓ú顏勬嵍妞ゆ挾鍋樻竟鏇熺節閳封偓閸愵喖寮伴悗瑙勬礃閸ㄥ潡鐛Ο灏栧亾濞戞顏堫敁閹惧墎纾藉ù锝呮贡閳藉銇勯幋婵囶棤缂佽京鍋炵换婵嬪礃瑜忕粻姘渻閵堝棗濮傞柛濠冩礋閵嗗倹绺介崨濠勫幈闂佸疇妗ㄧ欢銈夊礉閻旈晲绻嗛柛娆忣槸婵秶鈧娲栫紞濠囥€侀弴銏℃櫜闁告洦鍓欓弨顓炩攽閻樺灚鏆╅柛瀣☉铻炴繛鍡樻尭缁犳壆绱掔€n亞姘ㄩ柡瀣Ч閺屻劑鎮㈤崫鍕戯綁鏌涚€e墎鎮奸柟鑼归濂稿炊閿濆懍澹曞┑顔结缚閸樠囩嵁濮椻偓閺屻倕煤椤忓啯鍠氬銈冨灪閻熲晠骞冨▎鎿冩晢濞达絿枪婵酣姊婚崒娆戭槮闁规祴鍓濈粭鐔肺旈崨顓炲亶闂佺粯鎸哥€垫帡宕甸弴銏$厱鐟滃酣銆冮崱娑欏亗闁靛濡囩粻楣冩煙鐎电ǹ鍓辨繛鍫熒戦妵鍕晜婵傚憡顎嶇紓浣虹帛缁诲牓宕洪埀顒併亜閹烘垵顏撮柡浣稿€块弻娑㈠焺閸愵亝鍣梺鎸庣☉缁绘ê顫忓ú顏勭闁肩⒈鍓欑敮銉╂⒑閸濄儱校闁绘绮撳顐︻敋閳ь剟鐛惔銊﹀殟闁靛鍎伴崠鏍⒒娴e摜绉洪柛瀣躬瀹曟粌鈽夊杈╃厠闂佺粯鍨剁湁缂佽妫濋弻鐔虹磼濡櫣鐟ㄥ銈呮禋閸嬪﹪寮婚悢鍓叉Ч閹肩补瀵敐澶嬪€甸柣銏ゆ涧鐢爼鏌嶇拠鏌ュ弰妤犵偛顑夐幃顏堝川椤栨氨鍝�
2婵犵數濮烽弫鎼佸磻閻愬搫鍨傞柛顐f礀缁犳壆绱掔€n偓绱╂繛宸簻鎯熼梺鍐叉惈椤戝洨绮欒箛娑欌拺闁革富鍘奸崝瀣亜閵娿儲顥㈢€规洜鏁婚崺鈧い鎺戝閳锋垿鏌涘☉姗堝伐濠殿噯绠戦湁婵犲﹤鎳庢禒杈┾偓瑙勬礃濡炰粙寮幘缁樺亹鐎规洖娲ら獮姗€姊绘担鍛婃儓妞わ富鍨堕幃妯衡攽鐎n亝杈堝┑鐐叉閸ㄧ喖宕戦幘鑸靛枂闁告洦鍓涢敍姗€姊虹粙鎸庢崳闁轰礁顭烽悰顕€宕橀埡鍐炬祫闁诲函缍嗛崑鎺懳涢崘銊㈡斀闁绘劖娼欓悘锔剧棯閺夎法效妤犵偛顦垫俊鍫曞幢濞嗘埈鍟庨梻浣烘嚀椤曨參宕戦悙鏍稿寮介鐔蜂化婵炴挻鑹鹃敃锕傚箖閸忛棿绻嗛柛娆忣槸婵洭鏌嶇拠鏌ュ弰妤犵偛娲畷婊勬媴閾忕懓骞€婵犵數濮烽。钘壩i崨鏉戝瀭闂傚牊绋堥弸宥夋煥濠靛棙濯兼繛灏栨櫊閺屾洘绻涢悙顒佺彆闂佺粯鎸哥换姗€寮诲☉銏犵労闁告劏鏅濋ˇ銊х磽娴g柉鍏屽褎顨婃俊鐢稿礋椤栵絾鏅濋梺闈涚箚閺呮盯宕滈柆宥嗏拺闁圭ǹ娴烽埊鏇㈡煕閺冣偓閻熲晠鎮伴閿亾閿濆簼绨撮柡瀣叄閹﹢鎮欓弶鎴犱喊婵犮垼顫夊ú妯兼崲濞戞瑦缍囬柛鎾楀嫬浠归梺鍝勵儛娴滅偤鍩€椤掍緡鍟忛柛鐘崇墵閹ê鈹戠€n亞鍘撮梺纭呮彧缁犳垹澹曢崗鑲╃闁硅揪缍侀崫娲煟閿濆洦娅曠紒杈ㄦ尰閹峰懏顦版惔銈囩崶闂備胶枪椤戝懘鏁冮妶澶樻晪闁挎繂顦介弫鍡涙煃瑜滈崜鐔煎春閵忊剝鍎熼柍閿亾闁哄閰i弻鐔衡偓鐢登归灞句繆閸欏灏︽慨濠勭帛閹峰懘宕ㄦ繝鍌涙畼闂備浇宕甸崰鍡涘磿閻㈢ǹ鏋侀柛鎰靛枛鍥撮梺鍛婄矆缁€渚€寮查悩宕囩瘈闁汇垽娼ф牎濠殿喗菧閸斿矂鎮鹃崹顐ょ懝闁逞屽墴楠炲啫螖閸涱噮妫冨┑鐐村灦閻熴儵寮抽崼銉︹拺缂侇垱娲樺▍鍛存煕婵犲啰鎽冮柣蹇擃儔濮婅櫣绱掑鍡樼暥闂佺粯顨堥崑銈呯暦濡ゅ拋鏁嗛柛灞剧矌閿涙粓鏌℃径濠勫闁告柨娴风划濠氬箻缂佹ḿ鍘撻梺鑺ッˇ钘壩熼埀顒勬⒑鐠団€虫灆缂侇喗鐟ラ悾鐤亹閹烘繃鏅濋梺闈涚墕缁绘帡鎯岄敓锟�547闂傚倸鍊搁崐椋庣矆娴i潻鑰块梺顒€绉查埀顒€鍊圭粋鎺斺偓锝庝簽閿涙盯姊洪悷鏉库挃缂侇噮鍨堕崺娑㈠箳濡や胶鍘遍梺鍝勬处椤ㄥ棗鈻嶉崨瀛樼厽闊浄绲奸柇顖炴煛瀹€瀣埌閾绘牠鎮楅敐搴′簻妞ゅ骏鎷�4婵犵數濮烽弫鎼佸磻閻愬搫鍨傞柛顐f礀缁犳壆绱掔€n偓绱╂繛宸簼閺呮煡鏌涢妷銏℃珖妞わ富鍨跺娲偡闁箑娈堕梺绋款儑閸犳牠宕洪姀銈呯睄闁逞屽墴婵$敻宕熼鍓ф澑闂佽鍎抽顓⑺囬柆宥嗏拺缂佸顑欓崕鎰版煙閻熺増鎼愰柣锝呭槻椤粓鍩€椤掑嫨鈧線寮崼婵嗚€垮┑掳鍊曢崯顐︾嵁閹扮増鈷掗柛灞剧懅椤︼箓鏌涘顒夊剰妞ゎ厼鐏濋~婊堝焵椤掆偓閻g兘顢涢悜鍡樻櫇闂侀潧绻堥崹鍝勨枔妤e啯鈷戦梻鍫熶緱濡狙冣攽閳ヨ櫕鍠橀柛鈹垮灲瀵噣宕奸悢鍝勫箥闂備胶顢婇~澶愬礉閺囥垺鍎嶆繛宸簼閻撶喖鏌i弮鍫熸暠閻犳劧绱曠槐鎺撴綇閵娿儳鐟查悗鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊呭仦缁傛帡鎮℃惔妯绘杸闂佺粯鍔樺▔娑氭閿曞倹鐓曟俊銈呭閻濐亜菐閸パ嶅姛闁逞屽墯缁嬫帟鎽繝娈垮灡閹告娊骞冨畡鎵虫瀻婵炲棙鍨甸崺灞剧箾鐎涙ḿ鐭掔紒鐘崇墵瀵鈽夐姀鐘电杸闂佺ǹ绻愰幗婊堝极閺嶎厽鈷戠紒顖涙礃濞呮梻绱掔紒妯肩疄鐎殿喛顕ч埥澶娾堪閸涱垱婢戦梻浣瑰缁诲倿骞婃惔顭掔稏闁冲搫鎳忛埛鎴︽煕濞戞﹫鍔熼柟铏礈缁辨帗娼忛妸锔绢槹濡ょ姷鍋涚换姗€骞冮埡鍐╁珰闁肩⒈鍓﹂崯瀣⒒娴e憡鍟炲〒姘殜瀹曞綊骞庨崜鍨喘閸╋繝宕ㄩ瑙勫闂佽崵鍋炵粙鍫ュ焵椤掆偓閸樻牗绔熼弴銏♀拻濞达絽鎲$拹锟犲几椤忓棛纾奸柕濞垮妼娴滃湱绱掗鍛箺鐎垫澘瀚伴獮鍥敇閻樻彃绠婚梻鍌欑閹碱偆鈧凹鍓涢幑銏ゅ箳閺冨洤小闂佸湱枪缁ㄧ儤绂嶅⿰鍫熺厸闁搞儺鐓侀鍫熷€堕柤纰卞厴閸嬫挸鈻撻崹顔界彯闂佺ǹ顑呴敃銈夘敋閿濆洦宕夐悶娑掑墲閻庡姊虹拠鈥崇€婚柛蹇庡嫎閸婃繂顫忕紒妯诲闁荤喖鍋婇崵瀣磽娴e壊鍎愰柛銊ㄥ劵濡喎顪冮妶鍡樺蔼闁搞劌缍婇幃鐐哄垂椤愮姳绨婚梺鍦劋閸╁﹪寮ㄦ繝姘€垫慨妯煎亾鐎氾拷40缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎸鹃埀顒冾潐濞叉牕煤閵娧呬笉闁哄啫鐗婇悡娆撴煙椤栧棗鑻▓鍫曟⒑瀹曞洨甯涙慨濠傜秺楠炲牓濡搁妷顔藉缓闂侀€炲苯澧版繛鎴犳暬楠炴牗鎷呴崨濠勨偓顒勬煟鎼搭垳绉靛ù婊冪埣閹垽宕卞☉娆忎化闂佹悶鍎荤徊娲磻閹捐绀傞柛娑卞弾濡粎绱撻崒姘偓宄懊归崶銊d粓闁归棿鐒﹂崑锟犳煃閸濆嫭鍣归柦鍐枔閳ь剙鍘滈崑鎾绘煕閺囥劌浜炴い鎾存そ濮婃椽骞愭惔锝囩暤濠电偠灏欐繛鈧€规洘鍨块獮妯肩磼濡鍔掗梺鑽ゅ枑閻熴儳鈧凹鍓熷畷銏c亹閹烘挴鎷洪梺鍛婄箓鐎氼厼顔忓┑瀣厱閹兼番鍨归悘鈺備繆閸欏濮囨顏冨嵆瀹曞ジ鎮㈤崫鍕闂傚倷鑳剁涵鍫曞礈濠靛枹娲冀椤愩儱小缂備緡鍋勭€殿剟姊婚崒姘偓椋庢濮橆兗缂氱憸宥堢亱闂佸搫鍟崐濠氭儗閸℃褰掓晲閸偄娈欓梺鑽ゅ枑鐎氬牓寮崼婵嗙獩濡炪倖妫侀~澶屸偓鍨墵濮婄粯鎷呴崨濠傛殘婵炴挻纰嶉〃濠傜暦閵忋倖瀵犲璺烘閻庢椽鎮楅崗澶婁壕闂佸憡娲﹂崜娑㈠储闁秵鈷戦柛婵嗗閺嗙偤鏌熺粙鍨挃濠㈣娲熼獮鎰償濞戞鐩庨梻渚€娼ф蹇曟閺団偓鈧倿鎳犻鍌滐紲闂佸搫鍟崐鎼佸几濞戞瑣浜滈柕蹇婂墲缁€瀣煙椤旇娅婃い銏℃礋閿濈偤顢橀悜鍡橆棥濠电姷鏁搁崑鐘诲箵椤忓棛绀婇柍褜鍓氶妵鍕敃閵忊晜鈻堥梺璇″櫙缁绘繈宕洪埀顒併亜閹烘垵顏柍閿嬪浮閺屾稓浠﹂幑鎰棟闂侀€炲苯鍘哥紒顔界懇閵嗕礁鈻庨幇顔剧槇闂佸憡娲﹂崜锕€岣块悢鍏尖拺闁告挻褰冩禍婵囩箾閸欏澧辩紒顔垮吹缁辨帒螣闂€鎰泿闂備浇顫夊畷妯衡枖濞戙埄鏁佺€光偓閸曨剛鍘告繛杈剧到婢瑰﹪宕曡箛鏂讳簻妞ゆ挴鍓濈涵鍫曟煙妞嬪骸鈻堥柛銊╃畺瀹曟宕ㄩ娑樼樆闂傚倸鍊风欢姘跺焵椤掑倸浠滈柤娲诲灦瀹曘垽骞栨担鍦幘闂佸憡鍔樼亸娆撳春閿濆應鏀介柨娑樺閺嗩剟鏌熼鐣屾噰鐎殿喖鐖奸獮瀣敇閻愭惌鍟屾繝鐢靛У椤旀牠宕板Δ鍛櫇闁冲搫鎳庣粈鍌涚箾閹寸偟顣叉い顐f礋閺屻劌鈹戦崱妯轰痪閻熸粎澧楃敮妤呭疾閺屻儲鐓曢柍鈺佸暟閹冲懘鏌i幘鍐测偓鎼佲€旈崘顔嘉ч柛鎰╁妿娴犲墽绱掗悙顒佺凡缂佸澧庨崚鎺楀煛閸涱喖浜滅紒鐐妞存悂寮插┑瀣拺闂傚牊绋撴晶鏇熺箾鐠囇呯暤妤犵偛妫濋弫鎰緞鐎Q勫闂備礁婀辨灙婵炲鍏橀崺銉﹀緞鐎c劋绨婚梺鎸庢椤曆冾嚕椤曗偓閺屾盯鍩為幆褌澹曞┑锛勫亼閸婃牜鏁幒妤佹櫇闁靛/鈧崑鎾愁潩閻愵剙顏�28缂傚倸鍊搁崐鎼佸磹妞嬪孩顐介柨鐔哄Т绾捐顭块懜闈涘Е闁轰礁顑囬幉鎼佸籍閸垹绁﹂梺鍛婂姦閸犳牜绮绘繝姘厱闁规崘灏欑粣鏃堟煃閻熸壆绠茬紒缁樼箞婵偓闁挎繂妫涢妴鎰斿Δ濠佺凹闁圭ǹ鍟块悾宄扳攽鐎n亜绐涢柣搴㈢⊕宀e潡宕㈤柆宥嗏拺闁告繂瀚弳濠囨煕鐎n偅灏电紒杈ㄥ笧閳ь剨缍嗛崑鍛暦瀹€鈧埀顒侇問閸n噣宕戞繝鍥х畺濞寸姴顑呴崹鍌涖亜閹扳晛鐏╂鐐村灴濮婄粯鎷呴崨濠冨創濠电偠顕滅粻鎴︼綖濠靛惟闁冲搫鍊告禒顓㈡⒑鐎圭姵銆冮悹浣瑰絻鍗遍柛顐犲劜閻撴瑩鏌i幇闈涘缂傚秵鍨块弻鐔煎礂閸忕厧鈧劙鏌$仦鐣屝ユい褌绶氶弻娑㈠箻閸楃偛顫囧Δ鐘靛仜缁绘﹢寮幘缁樻櫢闁跨噦鎷�1130缂傚倸鍊搁崐鎼佸磹妞嬪海鐭嗗〒姘e亾閽樻繃銇勯弽銊х煂闁活厽鎹囬弻娑㈠箻閼碱剦妲梺鎼炲妽缁诲牓寮婚妸鈺傚亜闁告繂瀚呴姀銏㈢<闁逞屽墴瀹曟帡鎮欑€电ǹ骞堟繝鐢靛仦閸ㄥ爼鏁冮锕€缁╃紓浣贯缚缁犻箖鏌涢锝囩畼闁绘帗鎮傞弻锛勪沪缁嬪灝鈷夐悗鍨緲鐎氼噣鍩€椤掑﹦绉靛ù婊勭矒閿濈偤宕堕浣叉嫼闂備緡鍋嗛崑娑㈡嚐椤栨稒娅犳い鏇楀亾闁哄本绋掔换婵嬪礋椤掆偓濞堝矂鏌ч懡銈呬槐闁哄本娲熷畷鐓庘攽閸パ屾П濠电姵顔栭崹浼村Χ閹间礁钃熼柨婵嗩槸缁狅綁鏌h箛鏃€銇熷ù婊勭箘缁顓兼径濠囧敹闂佸搫娲ㄩ崰鎾诲储閹间焦鈷戦柛娑橈攻缁€瀣箾娴e啿娲﹂崐鍫曟煥濠靛棙顥犵紒鈾€鍋撻梻渚€鈧偛鑻晶鎾煕閳轰礁顏€规洘岣块幑鍕Ω瑜忛崢顒€鈹戦悩娈挎毌婵℃彃鎳樺畷褰掓濞戞碍娈鹃梻鍌氱墛娓氭宕崨顔剧瘈闂傚牊绋撴晶銏㈢棯閹规劖顥夐棁澶愭煥濠靛棙顕氱憸鐗堝笒濮规煡鏌¢崘銊у闁绘挸鍟伴幉绋库堪閸繄顦у┑鐐村灦绾板秹顢曢懞銉﹀弿婵☆垵娉曢崼顏堟煃瑜滈崜娆撳磹閸ф绠犻柣鎰惈鍞梺鎸庢磵閸嬫捇鏌涢弬鎸庡殗婵﹨娅i幉鎾礋椤愩垹袘闂備礁鎽滈崰搴♂缚閿熺姴绠栭柨鐔哄Т閸楁娊鏌i幇顒傂i柣鈺婂灦楠炲啳顦圭€规洖銈搁、妤呭焵椤掑嫬鐭楅柡鍥╁Х绾捐棄霉閿濆嫮鐭欓柛婵婃缁辨帞鎷犻幓鎺撴鐎光偓閿濆牊纭堕柟椋庡█瀹曪綁宕掑⿰鍐╃亖缂備浇椴哥敮鎺曠亽闂佺粯鍨靛ú銊╂偩鏉堛劊浜滄い鎾墲绾爼鏌熼悷鏉款伃闁圭厧缍婂畷婊堝箛閸撲胶锛熺紓浣介哺鐢繝宕规ィ鍐ㄧ鐎瑰壊鍠曠槐鏃堟煟鎼淬埄鍟忛柛锝庡櫍瀹曟垿宕熼姘卞幒闁瑰吋鐣崝宀€绮诲☉娆嶄簻闁规崘娉涘暩闂佸疇顕ч惌鍌氼潖缂佹ɑ濯撮柣鐔煎亰閸ゅ绱撴担鍓插剱闁搞劌澧庣紓鎾寸鐎n亞鐫勯梺绋挎湰缁苯鐣烽妷鈺傜厽閹兼番鍨婚埊鏇炵暆閿濆懏鍋ラ柛鈹垮灲瀵噣宕奸悢鍝勫箺婵犳鍠楅敃鈺呭礈濮樿缍栨繝鏇炲暞缁绘繈鎮介棃娑楁勃濠电偛鍚嬮悷鈺佺暦閵壯€鍋撻敐搴℃灍闁搞倕鑻灃闁挎繂鎳庨弳娆戠磼閻橀潧鏋涢柡宀€鍠栭獮鍡涙偋閸偅顥夐梻鍌欐祰閿熴儵宕愬┑瀣摕婵炴垯鍨瑰敮闂佺懓鐡ㄧ换鍌炲汲閵徛颁簻闁哄倹顑欏Σ鎼佹煃鐟欏嫬鐏撮柟顔规櫅闇夐悗锝庡亾缁辨岸姊绘笟鈧ḿ褔鎮ч崱娑樼疇閹兼番鍔岄崙鐘崇箾閸℃ɑ灏伴柍閿嬪灴閹綊宕堕敐鍌氫壕闁惧浚鍋嗘禍鏍⒒娴e憡鍟為柡灞诲妽椤ㄣ儵骞栨担鍝ュ幒闂佽宕橀褏绮婚敐澶嬬叆闁哄啫鍊婚幗鍌涙叏婵犲倹鍋ユ慨濠冩そ閺屽懘鎮欓懠璺侯伃婵犫拃灞芥珝鐎殿噮鍋婂畷濂稿Ψ閿旀儳骞堥梺璇茬箳閸嬫稓鏁崫鍕垫僵閻犻缚娅i敍娑㈡⒑閻熸澘鈷旂紒顕呭灦瀵煡骞栨担鍦弳闂佺粯娲栭崐鍦偓姘炬嫹
相关话题/计算 程序 优化 工程 空间

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 三维点阵结构等效热分析与优化方法*
    三维点阵结构,亦被称为“类桁架结构”[1],因具有轻质、高比刚度、高比强度的机械性能和主动散热/阻热、减震、吸能等功能特性,在航空、航天等领域受到广泛的关注。而三维点阵结构在研制与应用过程中,涉及到的关键问题包括制造工艺与性能分析方法等,对此,国内外****展开了一系列的研究。由于传统工艺的限制,点 ...
    本站小编 Free考研考试 2021-12-25
  • 高速开关阀的复合PWM控制策略分析与优化*
    电液伺服阀凭借其高频响、高精度等优势在航空航天、机器人及精密驱动等领域得到了广泛的应用[1]。但由于受到自身结构特性的约束,且对油液清洁度要求较高,电液伺服阀在高温、高压及强振动的环境下极易产生零漂或卡死等问题[2],长期工作可靠性难以保证。高速开关阀(High-SpeedOn/OffValve,H ...
    本站小编 Free考研考试 2021-12-25
  • 基于能量观点的混合层流优化设计*
    为了降低燃油消耗,缓解能源危机和环境问题,航空领域提出了“绿色航空”的概念,包含了提高发动机效率、完善飞机系统的一体化设计以及改善飞机气动特性等方案[1-2]。层流控制技术通过尽可能地维系飞机在巡航状态表面的层流来降低燃油消耗。对于A320客机,在机翼、尾翼以及短舱表面维系40%的层流,全机可以收益 ...
    本站小编 Free考研考试 2021-12-25
  • 城市中心区非机动车系统设计优化与探索*
    在城市发展中,积极推进非机动车、公共交通等绿色交通系统的建设、倡导绿色交通出行是可持续发展的城市交通理念的重要内容。随着共享单车在北京、上海、广州、深圳等大城市的逐渐推广,非机动车又重新走进了人们的视野,为了解决目前大部分城市非机动车道存在的机非混行、机动车随意占道停车、自行车与电力驱动的非机动车( ...
    本站小编 Free考研考试 2021-12-25
  • 复杂转子系统支点动载荷模型及其优化设计*
    随着航空燃气轮机向高转速、高负荷的方向发展,转子系统工作于多阶临界转速之上,工作中可能经过或靠近弯曲临界,导致转子不可避免地产生弯曲变形[1-2]。此时,转子受离心载荷及轮盘陀螺力矩的影响,使得支点动载荷不再是一个只由不平衡量单一因素引起的“较小的稳态交变力”。尤其是对于广泛使用的带有中介轴承的双转 ...
    本站小编 Free考研考试 2021-12-25
  • 基于矩阵分解的空间系绳系统不完全反馈控制*
    随着人类航天事业的不断发展,越来越多的航天器被不断发射到太空,因此近地空间已逐渐变成十分拥挤的场所。与此同时,由于复杂的太空扰动导致航天器轨道高度发生不同程度的变化,进而引起不同航天器之间的碰撞,产生大量空间碎片。因此,安全、高效地捕获空间碎片对于航天任务的安全完成具有重要意义。利用空间系绳系统(S ...
    本站小编 Free考研考试 2021-12-25
  • 高空太阳能无人机三维航迹优化*
    近年来,能源危机和环境污染已经成为人类社会面临的主要挑战。太阳能无人机依靠铺设在机翼表面的太阳能电池将太阳辐射能转化为电能,通过电机驱动螺旋桨产生前进动力,具有可持续、无污染的特点。高空太阳能无人机能够在平流层飞行数天乃至数月,可以广泛应用于遥感、侦察、通信等诸多领域,成为近年来各国研究的热点[1- ...
    本站小编 Free考研考试 2021-12-25
  • 一种面向模块化可重构机翼的分步补偿优化方法*
    随着航空结构设计与制造水平的不断提高,现代飞机尤其是军用飞机的任务剖面越来越复杂,模块化飞机设计作为一种可同时兼顾多任务需求和经济性约束的飞机设计思想而得到了国内外****的广泛关注[1-4]。模块化飞机是指将飞机主要结构划分成相对独立的模块进行设计、制造和使用的设计理念,基于模块化的飞机设计思想产 ...
    本站小编 Free考研考试 2021-12-25
  • 基于智能优化算法和有限元法的多线圈均匀磁场优化设计*
    磁场模拟装置在地磁导航等航空航天工程中具有重要的应用价值[1-3]。理想的磁场模拟装置应能够在足够大的空间内产生高均匀度的磁场,因此,磁场均匀性是衡量磁场模拟装置性能的关键技术指标之一。近年来,随着磁场模拟装置的应用领域越来越广泛,其磁场分布的均匀性得到了越来越多科研工作者的重视[4-6]。文献[4 ...
    本站小编 Free考研考试 2021-12-25
  • 质子交换炉温度场均匀性分析与优化*
    铌酸锂(LiNbO3)晶体具有优异的电光和声光等特性,被广泛用于制备各种光波导器件[1-2]。LiNbO3光波导是高精度光纤陀螺系统的核心器件,LiNbO3光波导的大批量生产是实现光纤陀螺产业化的关键之一。质子交换炉是用质子交换法制备LiNbO3光波导的主要设备,为了提高生产效率,往往将批量LiNb ...
    本站小编 Free考研考试 2021-12-25
闂傚倸鍊搁崐椋庣矆娴h櫣绀婂┑鐘插€寸紓姘辨喐閺冨牄鈧線寮介鐐茶€垮┑锛勫仧缁垶寮悩缁樷拺闂侇偆鍋涢懟顖涙櫠閹绢喗鐓熼柟鍨暙娴滄壆鈧娲栨晶搴ㄥ箲閸曨剚濮滈柡澶嬪閻庢娊姊婚崒娆戠獢闁逞屽墰閸嬫盯鎳熼娑欐珷濞寸厧鐡ㄩ悡鏇㈡煟濡崵鍙€闁告瑥瀚埀顒冾潐濞插繘宕归懞銉ょ箚闁割偅娲栭悙濠囨煏婵炲灝鍔村ù鍏兼礋濮婃椽鎳¢妶鍛€鹃梺鑽ゅ枂閸庢娊鍩€椤掍礁鍤柛娆忓暙椤曪綁骞庨挊澶愬敹闂侀潧顧€婵″洭宕㈤柆宥嗏拺鐟滅増甯掓禍浼存煕閹惧鎳囬柕鍡楀暙閳诲酣骞嬮悩纰夌床闂備礁鎲¢悷锕傛晪閻庤鎸稿Λ娑㈠焵椤掑喚娼愭繛鎻掔箻瀹曟繈骞嬮敃鈧弸渚€鏌熼崜褏甯涢柡鍛倐閺屻劑鎮ら崒娑橆伓闂傚倸鍊搁崐鐑芥倿閿旈敮鍋撶粭娑樻噽閻瑩鏌熸潏楣冩闁稿孩鏌ㄩ埞鎴﹀磼濮橆厼鏆堥梺绋款儑閸犳劗鎹㈠☉銏犵婵炲棗绻掓禒鑲╃磼缂併垹骞愰柛瀣崌濮婅櫣鎷犻弻銉偓妤佺節閳ь剚娼忛妸锕€寮块柣搴ㄦ涧閹芥粍绋夊鍡愪簻闁哄稁鍋勬禒锕傛煟閹惧瓨绀冪紒缁樼洴瀹曞崬螣閸濆嫷娼曞┑鐘媰鐏炶棄顫紓浣虹帛缁诲牓宕洪埀顒併亜閹烘垵顏╃紒鐘劜閵囧嫰寮埀顒勫磿閸愯尙鏆﹂柕澶堝劗閺€浠嬫煟閹邦剙绾фい銉у仱閹粙顢涘⿰鍐ф婵犵鈧磭鍩fい銏℃礋閺佸倿鎮剧仦钘夌闂傚倷鑳舵灙闁哄牜鍓涚划娆撳箻鐠囪尙鐤囬梺绯曞墲閻燂箓宕戦弽銊х闁糕剝蓱鐏忎即鏌i幘瀛樼濞e洤锕、娑樷枎閹烘繂濡抽梻浣呵圭€涒晠宕归崷顓燁潟闁规崘顕х壕鍏兼叏濡搫鎮戝Δ鏃堟⒒娓氣偓閳ь剛鍋涢懟顖涙櫠鐎涙ḿ绠惧ù锝呭暱閸氭ê鈽夊Ο閿嬵潔闂侀潧绻嗛埀顒€鍘栭崙鑺ョ節閻㈤潧孝闁挎洏鍊濋幃褎绻濋崶褏鏌у銈嗗笒鐎氼參鎮¢妷鈺傜厽闁哄洨鍋涢埀顒€婀遍埀顒佺啲閹凤拷