專利名稱：納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種納米材料力學(xué)性能測試技術(shù)，特別是涉及一種納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法。
背景技術(shù)：
隨著微機(jī)電和納機(jī)電系統(tǒng)的發(fā)展，零件和元器件趨于微型化，甚至達(dá)到納米數(shù)量 級，這要求傳統(tǒng)的加工工業(yè)由大尺度、粗獷式的工作方式發(fā)展到微納米尺度的精工細(xì)作。隨 著工程尺度的減小，加工的對象進(jìn)入到納米尺寸——產(chǎn)生了一個(gè)科學(xué)與工程交叉的新興領(lǐng) 域——納米工程。它包括納米切削、納米潤滑、納米磨削以及納米材料之間作用等。在納米 工程中，隨著工藝水平的提高，操作空間將進(jìn)一步縮小，直至進(jìn)入納米尺度。這時(shí)，我們操 作的對象不再是宏觀的塊體材料，而是具有納米尺寸的納米點(diǎn)、納米線或納米帶。在納米尺 度，由于化學(xué)效應(yīng)、小尺寸效應(yīng)、量子效應(yīng)和表面效應(yīng)的影響，材料的力學(xué)特性、缺陷、彈性 模量、載荷特性和失效機(jī)理等都將發(fā)生質(zhì)的變化，因此宏觀的機(jī)械定律不再適用，而基于量 子理論的原子、分子理論將逐漸占據(jù)主導(dǎo)地位。當(dāng)物質(zhì)的尺寸小于lOOnm，材料的總原子數(shù) 將小于1.0X IO7個(gè)原子或分子。在這個(gè)由原子和分子組成的幾何體中，化學(xué)鍵的非周期性 將對材料的性質(zhì)起主要作用。例如，在傳統(tǒng)的機(jī)械加工過程中，表面的變質(zhì)層一般有幾十納 米，如果我們加工只有IOOnm的零件，幾十納米的變質(zhì)層就會(huì)對器件的性質(zhì)帶來致命的影 響。因此，了解和把握納米尺度加工過程中表面層的化學(xué)性質(zhì)對納米零部件的加工和處理 至關(guān)重要，而建立在連續(xù)介質(zhì)基礎(chǔ)上的傳統(tǒng)加工理論和研究方法將不再適用，因此必須從 原子和分子水平上建立新的模型和新的研究方法——計(jì)算機(jī)仿真。仿真計(jì)算在物理學(xué)、化學(xué)、生物學(xué)和材料科學(xué)等許多科學(xué)領(lǐng)域中得到廣泛應(yīng)用，被 認(rèn)為是上世紀(jì)以來除理論分析和實(shí)驗(yàn)觀察之外的第三種科學(xué)研究手段，稱之為“計(jì)算機(jī)實(shí) 驗(yàn)”手段。目前，基于微觀理論的計(jì)算機(jī)仿真主要以分子動(dòng)力學(xué)和分子蒙特卡羅模擬為主。 但只有當(dāng)計(jì)算規(guī)模達(dá)到一定的水平，仿真計(jì)算才能真正起到溝通宏觀特性與微觀結(jié)構(gòu)的作 用，對于許多在理論分析和試驗(yàn)觀察上都難以了解的現(xiàn)象才可以做出一定的微觀解釋。任何物理現(xiàn)象的宏觀特性都是通過構(gòu)成該系統(tǒng)的大量原子、分子或電子等統(tǒng)稱為 粒子的運(yùn)動(dòng)狀態(tài)所決定的。而真實(shí)的物理系統(tǒng)所包含的微觀粒子數(shù)量往往高達(dá)IO23數(shù)量 級，另外由于原子振動(dòng)模擬的時(shí)間步長必須在飛秒(fs)級，模擬幾皮秒(PS)的現(xiàn)實(shí)時(shí)間就 需上千步，所以由于普通計(jì)算機(jī)CPU的性能和內(nèi)存空間的限制，在能接受的模擬時(shí)間范圍 內(nèi)，模擬規(guī)模只能達(dá)到數(shù)萬個(gè)原子，遠(yuǎn)遠(yuǎn)低于實(shí)驗(yàn)規(guī)模，因而分子動(dòng)力學(xué)模擬是典型的計(jì)算 瓶頸問題。因此，需要改進(jìn)算法，從而提高單機(jī)的計(jì)算能力，通過并行計(jì)算，大大提高分子動(dòng) 力學(xué)的仿真規(guī)模。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，通過 計(jì)算機(jī)仿真方法中分子動(dòng)力學(xué)的大規(guī)模優(yōu)化方法，實(shí)現(xiàn)納米工程中納米級切削、磨削、研磨等加工方式以及納米級壓痕、拉伸等納米材料力學(xué)性能測試。本發(fā)明的目的是通過以下技術(shù)方案實(shí)現(xiàn)的納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，該方法包括:MD串行法的系統(tǒng)劃 分和數(shù)據(jù)存儲(chǔ)方法的步驟(1)根據(jù)勢能函數(shù)的截?cái)喟霃絩。ut，首先將系統(tǒng)劃分成邊長1大于等于截?cái)喟霃?r-的η個(gè)三維立方體小網(wǎng)格，建立三維立方體小網(wǎng)格數(shù)組，將位于三維立方體小網(wǎng)格中的 原子信息存儲(chǔ)在鏈表數(shù)組中；(2)根據(jù)系統(tǒng)計(jì)算過程中工件原子位置變化的幅度和與刀具或壓頭原子的相互作 用的大小，確定工件與刀具或壓頭的作用距離R ；(3)建立系統(tǒng)中所有原子的相互作用列表，首先確定與系統(tǒng)原子i相互作用的原 子j，即在包括該原子i在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)網(wǎng)格中、二維是9個(gè)網(wǎng)格，確定 與該原子i距離小于等于原子截?cái)喟霃絩。ut+Ar的所有原子j，將j原子以鏈表形式存儲(chǔ)在 i原子相互作用列表中；(4)在作用距離R內(nèi)建立工件原子與刀具或壓頭原子間的計(jì)算列表，即在作用距 離R內(nèi)，在包括該原子i在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)網(wǎng)格中、二維是9個(gè)網(wǎng)格，確定 與該原子i距離小于等于原子截?cái)喟霃絩。ut+Ar的原子j，將j原子存儲(chǔ)在相應(yīng)的i原子鏈 表中；(5)根據(jù)仿真過程中工件原子位置變化和ΔΓ大小，確定計(jì)算列表的更新頻率，根 據(jù)刀具或壓頭的位置和相互作用距離R，更新刀具或壓頭與工件的作用范圍；(6)在系統(tǒng)計(jì)算過程中，計(jì)算每個(gè)原子與其他原子之間的作用，僅需計(jì)算原子與其 相應(yīng)列表中原子之間的作用，而與系統(tǒng)中其他原子之間的作用無需判定。所述的納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，其所述的MD串行法的系 統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)方法，根據(jù)仿真過程中工件原子位置變化的頻率，確定原子的相互作用 列表更新頻率；根據(jù)刀具或壓頭的位置，以刀具或壓頭為中心，更新刀具或壓頭與工件的作 用范圍R。所述的納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，其所述的MD串行法的系 統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)方法，對MD并行算法節(jié)點(diǎn)機(jī)計(jì)算任務(wù)的劃分和數(shù)據(jù)存儲(chǔ)的具體步驟如 下(1)根據(jù)工件計(jì)算過程中工件的變形情況，確定變形最小的方向；(2)在單機(jī)MD串行算法的系統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)的基礎(chǔ)上，將系統(tǒng)劃分成η個(gè)三維 立方體小網(wǎng)格；為了并行計(jì)算過程中負(fù)載平衡，沿著變形最小的方向?qū)⑾到y(tǒng)分成N等分，這 N個(gè)部分被稱作計(jì)算子區(qū)域，每個(gè)子區(qū)域?qū)?yīng)一臺(tái)節(jié)點(diǎn)機(jī)、共有N臺(tái)節(jié)點(diǎn)機(jī)；每個(gè)子區(qū)域計(jì) 算任務(wù)以三維立方體小網(wǎng)格為單位，由m >= 1 ( Σ NXm = η)個(gè)網(wǎng)格組成，根據(jù)節(jié)點(diǎn)機(jī)的 性能，每臺(tái)處理機(jī)的網(wǎng)格數(shù)m可不相等，只要每臺(tái)節(jié)點(diǎn)機(jī)每次計(jì)算任務(wù)時(shí)間大致相等；(3)系統(tǒng)中原子位置坐標(biāo)、速度 和所受的力等信息分機(jī)保存，各節(jié)點(diǎn)機(jī)之間只傳遞 鄰近原子位置坐標(biāo)等計(jì)算必需的信息，數(shù)據(jù)處理所需的切削力、原子徑向分布、能量和表面 能等重要信息也采取分機(jī)保存，最后用一個(gè)合并程序?qū)⒏鳈C(jī)數(shù)據(jù)合并，用可視化程序動(dòng)態(tài) 顯示切削過程。


圖1為本發(fā)明為減少計(jì)算時(shí)間，建立系統(tǒng)中原子相互作用原子表的二維示意圖；圖2為本發(fā)明并行計(jì)算過程中節(jié)點(diǎn)機(jī)計(jì)算任務(wù)的劃分示意圖；
圖3為本發(fā)明表示并行計(jì)算過程中節(jié)點(diǎn)機(jī)之間數(shù)據(jù)和信息相互傳遞的示意圖。
具體實(shí)施例方式下面參照附圖對本發(fā)明進(jìn)行詳細(xì)說明。從理論上來說，在計(jì)算原子間相互作用力時(shí)，每個(gè)原子與系統(tǒng)中其余所有原子的 相互作用都需要進(jìn)行計(jì)算；但是由于原子間的作用力隨距離增加而迅速減小，因此，沒有必 要計(jì)算每個(gè)原子與系統(tǒng)中其余所有原子的相互作用。本發(fā)明提出一種減少不必要計(jì)算量的 串行算法，同時(shí)通過優(yōu)化數(shù)據(jù)存儲(chǔ)，節(jié)約內(nèi)存，提高串行MD計(jì)算的效率。通過并行計(jì)算，從 而達(dá)到大大提高M(jìn)D仿真規(guī)模。本發(fā)明根據(jù)原子的物理位置將系統(tǒng)中所有原子劃分成三維 小網(wǎng)格，根據(jù)系統(tǒng)原子在仿真過程中的運(yùn)動(dòng)情況，確定原子間相互作用的范圍。在原子所在 網(wǎng)格和周圍的網(wǎng)格中，根據(jù)相互作用范圍，建立系統(tǒng)中每個(gè)原子的相互作用原子表。當(dāng)系統(tǒng) 中涉及刀具或壓頭時(shí)，根據(jù)刀具或壓頭對工件的作用，確定兩者之間的作用范圍。在兩者相 互作用范圍內(nèi)，建立工件原子或刀具、壓頭原子的相互作用原子表。根據(jù)原子運(yùn)動(dòng)情況，確 定更新原子相互作用原子表頻率。根據(jù)刀具或壓頭的位置，更新刀具或壓頭與工件的作用 范圍。根據(jù)工件仿真過程中工件的變形情況，確定工件變形最小的方向；為了并行計(jì)算過程 中負(fù)載平衡，沿著變形最小的方向?qū)⑾到y(tǒng)分成N等分，這N個(gè)部分被稱作計(jì)算子區(qū)域，每個(gè) 子區(qū)域?qū)?yīng)一臺(tái)節(jié)點(diǎn)機(jī)(共有N臺(tái)節(jié)點(diǎn)機(jī))，根據(jù)節(jié)點(diǎn)機(jī)的性能，每臺(tái)處理機(jī)的計(jì)算任務(wù)可 不相等，只要每臺(tái)節(jié)點(diǎn)機(jī)每次計(jì)算任務(wù)時(shí)間大致相等就可。本發(fā)明圖1是為減少計(jì)算時(shí)間，建立系統(tǒng)中原子相互作用原子表的二維示意圖。 表示根據(jù)系統(tǒng)選取的勢函數(shù)截?cái)喟霃絩。ut將系統(tǒng)劃分成邊長為1大于等于截?cái)喟霃絩。ut的 小網(wǎng)格，在原子i所在網(wǎng)格1和周圍的網(wǎng)格2 9共計(jì)9個(gè)網(wǎng)格(三維是27個(gè)網(wǎng)格)中， 建立原子i的相互作用原子列表。本發(fā)明圖2是指并行計(jì)算過程中節(jié)點(diǎn)機(jī)計(jì)算任務(wù)的劃分示意圖。在仿真過程中， 原子密度不斷發(fā)生變化，在Z軸方向工件變形很小，基本可以忽略，因此為了計(jì)算過程中負(fù) 載平衡，首先將整個(gè)工件劃分成η個(gè)三維正方體小網(wǎng)格(網(wǎng)格邊長等于晶格常數(shù))，沿Z軸 方向分成N等分，如圖中虛長方體所示，這N個(gè)部分被稱作計(jì)算子區(qū)域，每個(gè)子區(qū)域?qū)?yīng)一 臺(tái)節(jié)點(diǎn)機(jī)(共有N臺(tái)節(jié)點(diǎn)機(jī))，每個(gè)子區(qū)域由m> = KNXm = η)個(gè)網(wǎng)格組成，根據(jù)節(jié)點(diǎn)機(jī) 的性能，每臺(tái)處理機(jī)的網(wǎng)格數(shù)m可不相等，只要每臺(tái)節(jié)點(diǎn)機(jī)每次計(jì)算任務(wù)時(shí)間大致相等。本發(fā)明圖3表示并行計(jì)算過程中節(jié)點(diǎn)機(jī)之間數(shù)據(jù)和信息相互傳遞的示意圖。判別 新時(shí)刻每臺(tái)節(jié)點(diǎn)機(jī)上的原子位置，目的在于更新每個(gè)網(wǎng)格的原子鏈，同時(shí)還要判斷越出邊 界的原子點(diǎn)。在更新原子坐標(biāo)的同時(shí)，將邊長為1的boundary邊界區(qū)域的原子信息打包 (分上下boundary邊界區(qū)域打兩個(gè)包)，如圖a所示，同時(shí)將從P節(jié)點(diǎn)機(jī)運(yùn)動(dòng)到P+1機(jī)的原 子信息和P節(jié)點(diǎn)機(jī)中靠近P+1機(jī)boundary網(wǎng)格行中原子信息一起打包發(fā)送到P+1機(jī)，同時(shí) 將運(yùn)動(dòng)到P-I機(jī)的原子信息和P節(jié)點(diǎn)機(jī)中靠近P-I機(jī)boundary網(wǎng)格行中原子信息一起打 包發(fā)送到P-I機(jī)。為了防止將剛發(fā)送出去的原子信息被發(fā)送回原節(jié)點(diǎn)機(jī)，將從P節(jié)點(diǎn)機(jī)上 打包、收拆包裹的順序需按照圖b所示的順序。
1.首先根據(jù)系統(tǒng)材料選取合適的勢函數(shù)，再根據(jù)勢函數(shù)，選取合適的截?cái)喟霃絩cut.2.根據(jù)截?cái)喟霃絩。ut，將系統(tǒng)劃分成邊長為1大于等于截?cái)喟霃絩。ut的三維小網(wǎng) 格，建立三維網(wǎng)格的三維數(shù)組，將位于三維立方體小網(wǎng)格中原子信息以鏈表結(jié)構(gòu)存儲(chǔ)在網(wǎng) 格三維數(shù)組中，以便于原子的更新和快速查找；3.根據(jù)系統(tǒng)計(jì)算過程中工件原子位置變化的幅度和與刀具或壓頭原子的相互作 用的大小，確定工件與刀具或壓頭的作用距離R ；4.建立工件所有原子i的計(jì)算列表，首先確定工件原子間需要計(jì)算的原子，即在 包括該原子在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)(二維是9個(gè)網(wǎng)格)網(wǎng)格中，確定與該原子 距離小于等于原子截?cái)喟霃絩。ut+Ar的原子以鏈表結(jié)構(gòu)形式存儲(chǔ)在計(jì)算列表中；5.在作用距離R內(nèi)建立工件原子與刀具或壓頭原子間的計(jì)算列表，即在作用距離 R內(nèi)，在包括該原子在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)(二維是9個(gè)網(wǎng)格)網(wǎng)格中，確定與 該原子距離小于等于原子截?cái)喟霃絩。ut+Ar的原子以鏈表結(jié)構(gòu)形式存儲(chǔ)在計(jì)算列表中；6.根據(jù)仿真過程中工件原子位置變化和ΔΓ大小，確定計(jì)算列表的更新頻率。根 據(jù)刀具或壓頭的位置和相互作用距離R，更新刀具或壓頭與工件的作用范圍；7.根據(jù)工件計(jì)算過程中工件的變形情況，確定變形最小的方向；為 了并行計(jì)算過 程中負(fù)載平衡，將已被劃分成η個(gè)三維立方體小網(wǎng)格的系統(tǒng)沿著變形最小的方向再劃分成 N等分，這N個(gè)部分被稱作計(jì)算子區(qū)域，每個(gè)子區(qū)域?qū)?yīng)一臺(tái)節(jié)點(diǎn)機(jī)(共有N臺(tái)節(jié)點(diǎn)機(jī))，每 個(gè)子區(qū)域由m >= 1 ( Σ NXm = η)個(gè)網(wǎng)格組成，根據(jù)節(jié)點(diǎn)機(jī)的性能，每臺(tái)處理機(jī)的網(wǎng)格數(shù) m可不相等，只要每臺(tái)節(jié)點(diǎn)機(jī)每次計(jì)算任務(wù)時(shí)間大致相等。
權(quán)利要求
納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，其特征在于該方法包括MD串行法的系統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)方法的步驟(1)根據(jù)勢能函數(shù)的截?cái)喟霃絩cut，首先將系統(tǒng)劃分成邊長1大于等于截?cái)喟霃絩cut的n個(gè)三維立方體小網(wǎng)格，建立三維立方體小網(wǎng)格數(shù)組，將位于三維立方體小網(wǎng)格中的原子信息存儲(chǔ)在鏈表數(shù)組中；(2)根據(jù)系統(tǒng)計(jì)算過程中工件原子位置變化的幅度和與刀具或壓頭原子的相互作用的大小，確定工件與刀具或壓頭的作用距離R；(3)建立系統(tǒng)中所有原子的相互作用列表，首先確定與系統(tǒng)原子i相互作用的原子j，即在包括該原子i在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)網(wǎng)格中、二維是9個(gè)網(wǎng)格，確定與該原子i距離小于等于原子截?cái)喟霃絩cut+Δr的所有原子j，將j原子以鏈表形式存儲(chǔ)在i原子相互作用列表中；(4)在作用距離R內(nèi)建立工件原子與刀具或壓頭原子間的計(jì)算列表，即在作用距離R內(nèi)，在包括該原子i在內(nèi)的周圍立方體小網(wǎng)格共計(jì)27個(gè)網(wǎng)格中、二維是9個(gè)網(wǎng)格，確定與該原子i距離小于等于原子截?cái)喟霃絩cut+Δr的原子j，將j原子存儲(chǔ)在相應(yīng)的i原子鏈表中；(5)根據(jù)仿真過程中工件原子位置變化和Δr大小，確定計(jì)算列表的更新頻率，根據(jù)刀具或壓頭的位置和相互作用距離R，更新刀具或壓頭與工件的作用范圍；(6)在系統(tǒng)計(jì)算過程中，計(jì)算每個(gè)原子與其他原子之間的作用，僅需計(jì)算原子與其相應(yīng)列表中原子之間的作用，而與系統(tǒng)中其他原子之間的作用無需判定。
2.根據(jù)權(quán)利要求1所述的納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，其特征在 于，所述的MD串行法的系統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)方法，根據(jù)仿真過程中工件原子位置變化的頻 率，確定原子的相互作用列表更新頻率；根據(jù)刀具或壓頭的位置，以刀具或壓頭為中心，更 新刀具或壓頭與工件的作用范圍R。
3.根據(jù)權(quán)利要求2所述的納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，其特征在 于，所述的MD串行法的系統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)方法，對MD并行算法節(jié)點(diǎn)機(jī)計(jì)算任務(wù)的劃分和 數(shù)據(jù)存儲(chǔ)的具體步驟如下(1)根據(jù)工件計(jì)算過程中工件的變形情況，確定變形最小的方向；(2)在單機(jī)MD串行算法的系統(tǒng)劃分和數(shù)據(jù)存儲(chǔ)的基礎(chǔ)上，將系統(tǒng)劃分成n個(gè)三維立方 體小網(wǎng)格；為了并行計(jì)算過程中負(fù)載平衡，沿著變形最小的方向?qū)⑾到y(tǒng)分成N等分，這N個(gè) 部分被稱作計(jì)算子區(qū)域，每個(gè)子區(qū)域?qū)?yīng)一臺(tái)節(jié)點(diǎn)機(jī)、共有N臺(tái)節(jié)點(diǎn)機(jī)；每個(gè)子區(qū)域計(jì)算任 務(wù)以三維立方體小網(wǎng)格為單位，由m> = 1( E NXm = n)個(gè)網(wǎng)格組成，根據(jù)節(jié)點(diǎn)機(jī)的性能， 每臺(tái)處理機(jī)的網(wǎng)格數(shù)m可不相等，只要每臺(tái)節(jié)點(diǎn)機(jī)每次計(jì)算任務(wù)時(shí)間大致相等；(3)系統(tǒng)中原子位置坐標(biāo)、速度和所受的力等信息分機(jī)保存，各節(jié)點(diǎn)機(jī)之間只傳遞鄰近 原子位置坐標(biāo)等計(jì)算必需的信息，數(shù)據(jù)處理所需的切削力、原子徑向分布、能量和表面能等 重要信息也采取分機(jī)保存，最后用一個(gè)合并程序?qū)⒏鳈C(jī)數(shù)據(jù)合并，用可視化程序動(dòng)態(tài)顯示 切削過程。
全文摘要
納米加工中大規(guī)模分子動(dòng)力學(xué)的并行優(yōu)化方法，涉及一種納米材料力學(xué)性能測試技術(shù)。首先根據(jù)描述系統(tǒng)相互作用勢的截?cái)喟霃?，將系統(tǒng)劃分成三維立方體小網(wǎng)格。除了邊界網(wǎng)格外，在相鄰27個(gè)小網(wǎng)格(二維9個(gè))組成的立方體(二維正方形)上，建立中心網(wǎng)格中原子的臨近原子相互作用表。通過該法，建立系統(tǒng)中所有原子的臨近原子相互作用表。然后，根據(jù)系統(tǒng)中刀具、壓頭(如有刀具和壓頭)與工件的相互作用范圍，在以刀具、壓頭為中心的相互作用范圍內(nèi)，建立該范圍內(nèi)工件原子與刀具、壓頭原子的臨近原子相互作用表。最后，以三維立方體小網(wǎng)格為單位，進(jìn)行并行程序中節(jié)點(diǎn)機(jī)計(jì)算任務(wù)的劃分。來提高納米工程中MD的仿真規(guī)模。
文檔編號G06F17/50GK101859336SQ20101017754
公開日2010年10月13日 申請日期2010年5月20日 優(yōu)先權(quán)日2010年5月20日
發(fā)明者吳玉厚, 唐玉蘭, 孫紅, 張珂 申請人:沈陽建筑大學(xué)