鄧奕銘，王彥濤，惠衛(wèi)華

(1.北京機電工程研究所，北京 100074；2.鄭州機電工程研究所，鄭州 450015；3.西北工業(yè)大學 燃燒、熱結(jié)構(gòu)與內(nèi)流場重點實驗室，西安 710072)

0 引言

固體推進劑作為固體火箭發(fā)動機的重要構(gòu)成部件之一，其力學性能將直接關系到發(fā)動機能否正常工作[1]。固體火箭發(fā)動機裝藥在溫度沖擊載荷下容易出現(xiàn)裂紋和脫粘，使發(fā)動機結(jié)構(gòu)完整性遭到破壞，從而影響發(fā)動機的安全性和使用可靠性[2]。但是長期以來，溫度沖擊下的裝藥力學響應仿真分析過程復雜，三維建模與有限元計算軟件前處理之間的數(shù)據(jù)傳遞過程費時費力、效率低下，給發(fā)動機裝藥設計工作帶來了極大的不便[3]。而且溫度沖擊實驗是每個型號的發(fā)動機在列裝前必須進行的，如果能找到一種快速有效的方法來分析溫度沖擊裝藥結(jié)構(gòu)響應，在發(fā)動機設計完成前驗證裝藥的溫度沖擊性能，可有效提高發(fā)動機設計效率，避免時間浪費。在藥柱結(jié)構(gòu)完整性方面國內(nèi)外學者有很多研究，H C Yildirim和S Ozupek以及Heller等[4-6]研究了低溫貯存溫度下固體推進劑的應力應變和累積損傷；鄧康清等[7]通過熱-機耦合仿真分析，分別模擬了固體火箭發(fā)動機在降至常溫的慢速固化降溫過程和降至低溫的快速降溫過程。在參數(shù)化建模方面，李磊等[8]利用MSC.Patran軟件的二次開發(fā)工具PCL，實現(xiàn)了圓管傘盤型藥柱、星型藥柱、車輪型藥柱的二維參數(shù)化建模分析。

上述固體火箭發(fā)動機參數(shù)化建模研究中對推進劑藥柱進行參數(shù)化建模時可能會過多地簡化藥柱模型或忽略藥柱的一些微小幾何特征，無法全面地反映藥柱的幾何形狀，而對裝藥結(jié)構(gòu)響應的研究中對不均勻溫度場對裝藥結(jié)構(gòu)響應的影響也分析甚少。因此，研究溫度沖擊下的裝藥結(jié)構(gòu)響應快速分析方法具有重要意義。其能簡化繁瑣的發(fā)動機裝藥設計仿真流程，提高設計、分析效率，降低試驗、計算成本，而且能夠?qū)崿F(xiàn)特定工況下的分析計算功能。

1 裝藥參數(shù)化建模

1.1 裝藥幾何分析及參數(shù)化

固體火箭發(fā)動機從初樣設計到最終定型是一個不斷重復的過程，要歷經(jīng)多輪設計，但每次修改僅是小范圍的尺寸修改，如殼體厚度，藥柱肉厚等，如果每次修改都要求設計人員重新開始，將極大地增加設計人員的負擔，而固體發(fā)動機的參數(shù)化建模能很好地解決這一問題。

參數(shù)化建模[9-10]有基于特征、尺寸約束、關系驅(qū)動等特點，因此在參數(shù)化建模的過程中要注意以下要點：

(1)參數(shù)化幾何模型要力求精確還原實際物體基本特征，如星型藥柱要準確反映星根倒角、星尖弧度等基礎特征；

(2)尺寸鏈完備且閉合，要有完備的尺寸鏈約束且不能過約束，過約束會導致參數(shù)傳遞錯誤，造成尺寸鏈無法閉合，從而引起參數(shù)無法修改的情況；

(3)選擇適當?shù)年P系驅(qū)動，盡量選取關鍵參數(shù)舍棄次要參數(shù)。

參數(shù)化設計起始于草圖，進而形成實體模型。利用草圖進行相關參數(shù)化尺寸驅(qū)動，然后對草圖上的對象進行位置關系約束。需要注意的是，建模時有且只能有一個主特征，其他的特征都利用它的基準點等進行定位，從而確定位置關系[11-13]。

藥柱通常是由藥柱基體和藥柱芯模通過布爾差構(gòu)造的。在對藥柱型面進行修改時，往往藥柱外輪廓不動而對內(nèi)型面進行微小改動，這樣就只需修改藥柱芯模，再運用一次布爾差即可，可以避免藥柱外型面和絕熱層之間的聯(lián)動而造成的錯誤。

1.1.1 藥柱基體

藥柱基體即為實心藥柱。它的主要參數(shù)包括長度L、橢球封頭的長軸L3、短軸L2、藥柱外半徑R等參數(shù)。典型的星型藥柱外形輪廓如圖1所示。藥柱基體是典型的回轉(zhuǎn)體，在得到藥柱外輪廓后，將此外輪廓沿藥柱中軸旋轉(zhuǎn)360°后可得到三維藥柱基體。

圖1 藥柱基體外輪廓尺寸

1.1.2 藥柱芯模

星孔藥柱的芯?？梢酝ㄟ^星形草圖拉伸得到，星孔參數(shù)較為固定，如表1所示。星孔的幾何參數(shù)化輪廓及得到的實體芯模如圖2所示。

1.1.3 藥柱裝配

在完成藥柱基體和芯模的參數(shù)化建模后，將芯模端面和藥柱基體端面進行匹配定位，最后采用布爾差運算完成星型藥柱的建模。圖3為運用該方法得到的星型藥柱。

表1 星型內(nèi)腔參數(shù)化方案

圖2 二維星孔建模參數(shù)及星孔芯模

圖3 星型藥柱實體模型

1.1.4 裝配體建模

在藥柱、殼體、絕熱層模型參數(shù)化完成后，為方便模型的修改，建立裝配體參數(shù)，然后將各個參數(shù)關聯(lián)到零件的參數(shù)上。這樣就可以實現(xiàn)只需修改裝配參數(shù)，零件參數(shù)就可自動修改，進而整個模型(包括裝配體中所有零件)實現(xiàn)參數(shù)聯(lián)動，即實現(xiàn)重生所有零件，并且還能保證所有零件配合依然正確。

經(jīng)過參數(shù)與關系的設置，即可最終完成整個模型的參數(shù)聯(lián)動，完整的裝配體模型如圖4所示。

圖4 星型藥參數(shù)化完成后的裝配體

1.2 裝藥建模自動化實現(xiàn)

在完成裝藥的參數(shù)化建模后，運用Pro/E二次開發(fā)來實現(xiàn)裝藥模型的修改、再生、保存等，二次開發(fā)涉及到調(diào)用函數(shù)對應關系如圖5所示?？梢钥闯觯诙伍_發(fā)中比較重要功能有加載參數(shù)、修改參數(shù)、模型再生、導出模型等幾個關鍵步驟，下面對著幾個功能進行介紹[14-15]。

圖5 Pro/E與J-Link函數(shù)對應使用流程

(1)加載參數(shù)

設計模板涉及多個裝藥參數(shù)模型，每個裝藥模型又涉及多個參數(shù)，在單個設計模板中Pro/E隨機返回參數(shù)，存在隨機性，這樣程序不夠穩(wěn)健。為加強穩(wěn)健性，采用“參數(shù)名-取值”的結(jié)構(gòu)，返回列表中的數(shù)據(jù)也采用這種結(jié)構(gòu)，可有效聯(lián)動圖形用戶界面，方便用戶操作，實現(xiàn)模型參數(shù)界面修改

(2)修改參數(shù)

在參數(shù)加載完成后，同一模板各參數(shù)不可能一直保持不變，否則就失去了參數(shù)化建模的意義，如果想對特定尺寸進行修改，這時就要用到Pro/E二次開發(fā)中的接口Parameter.SetValue()函數(shù)，該函數(shù)可將參數(shù)修改值進行傳遞。

(3)執(zhí)行再生

在數(shù)據(jù)參數(shù)修改完成后，執(zhí)行Pro/E的模型再生功能，讀入修改后尺寸數(shù)據(jù)，進行重新造型。再生步驟利用Pro/E提供的Solid. Regenerate()函數(shù)來完成，該函數(shù)可對Solid類型的模型進行操作。

(4)模型導出

Pro/E生成的部件和裝配體的默認格式為*.prt和*.asm，但本文對裝藥溫度沖擊的仿真計算時需要模型往往是*.x_t或*.step格式文件，這就需要模型導出時進行文件格式的轉(zhuǎn)化，借助J-Link中的Model.Export()函數(shù)，即可將零件或裝配體轉(zhuǎn)化為所需文件類型。

1.3 溫度沖擊過程

根據(jù)GJB 5021—2001有關溫度沖擊實驗的規(guī)定，溫度沖擊循環(huán)次數(shù)一般規(guī)定為3次，如圖6所示。發(fā)動機溫度沖擊保溫時間與保溫厚度之間的關系如表2所示。保溫時間指藥柱內(nèi)腔溫度達到保溫溫度后，還需要繼續(xù)保溫的時間。

圖6 溫度沖擊實驗流程

本次仿真所用星形裝藥發(fā)動機的裝藥肉厚為60 mm，故根據(jù)溫沖實驗保溫要求，選擇保溫時間為t≥24 h。

表2 保溫時間與厚度關系

本文采用第三類熱傳導邊界條件進行溫度沖擊仿真，由流體橫掠圓管表面平均表面?zhèn)鳠嵯禂?shù)的關聯(lián)式、相關物性參數(shù)和發(fā)動機幾何參數(shù)可得：調(diào)溫期間對流換熱系數(shù)h=8.0 W/(m2·K)；高溫沖擊期間對流換熱系數(shù)h=13.0 W/(m2·K)；低溫沖擊期間對流換熱系數(shù)h=9.0 W/(m2·K)。

對應圖6實驗流程設定溫度沖擊裝藥的仿真流程如下，裝藥模型初始溫度為20 ℃，然后對發(fā)動機殼體外表面施加8.0 W/(m2·K)的對流邊界條件，持續(xù)48 h；緊接著高溫沖擊對發(fā)動機殼體外表面施加13 W/(m2·K)的對流邊界條件，也保持48 h；最后進行低溫沖擊時對殼體外表面施加9.0 W/(m2·K)的對流邊界條件，同樣保持48 h，第一次高低溫仿真沖擊完成，然后按照此流程重復3次即可，這時溫度沖擊裝藥結(jié)構(gòu)完整性仿真分析完成。

2 溫度沖擊裝藥結(jié)構(gòu)響應計算模板

在進行典型裝藥模板的設計時，為保證模板的通用性和精準性，要建立一套標準的分析流程，用通用的方法完成網(wǎng)格的剖分和離散、邊界條件的設置、載荷工況的加載等，在保證精確度的情況下，盡可能的減少人工干預和工作量，保持模板的合理性。

典型的溫度沖擊裝藥模板的建立流程如下：

(1)對目標裝藥的溫度沖擊過程進行詳細分析，確定所建立的有限元模型的理想假設情況；(2)確定有限元模型的幾何參數(shù)、材料參數(shù)、邊界條件、載荷工況等，根據(jù)裝藥特征確定重點關注部位，如星型藥柱的星尖位置；(3)確定裝藥的網(wǎng)格離散方法及通用建模方法；(4)在通用建模方法的基礎上，提取*.rpy日志文件，并在提取的*.rpy日志文件的基礎上添加相關命令流，實現(xiàn)參數(shù)化的定義；(5)編寫GUI，使Pro/E的參數(shù)化模型和Python腳本得以聯(lián)動，實現(xiàn)數(shù)據(jù)的內(nèi)部傳遞及計算參數(shù)的靈活輸入。

2.1 控制方程和物性參數(shù)

推進劑材料的熱粘彈性本構(gòu)方程為[15]：

(1)

式中t為具體某時刻；τ為時間變量；T0為參考溫度；εT(τ)為熱應變，可通過熱膨脹系數(shù)α0和溫度變化來表示；ζ(t)為時間折減因子；ζ(t)=t/αT，αT為溫度-時間移位因子。

其他裝藥具體材料參數(shù)如表3所示。

該藥柱為典型的粘彈性材料，采用大應變模型來描述，用Prony級數(shù)的形式來表征其松弛模量，如式(2)所示。

(2)

其前6階參數(shù)如表4所示。

藥柱材料的時間-溫度等效因子的WLF方程為

(3)

式中C1，C2為材料常數(shù)，由材料自身特性決定，C1=17，C2=52；T0為基準溫度，T0=20 ℃。

表3 裝藥材料參數(shù)

表4 藥柱粘彈性材料參數(shù)

確定完裝藥各部件的材料參數(shù)后，以藥柱材料的傳導率為例，定義材料參數(shù)代碼如下：

mdb.models[′Model-1′].materials[′grain′].Conductivity(table=daorexishu_yao)

對實體部件賦予不同材料屬性：

p=mdb.models[′Model1′].parts[′GRAIN_XING′]

c=p.cells

region=p.Set(cells=cells, name=′Set-1′)

c1=mdb.models[′Model-1′].parts[′grain′] .cells

p.SectionAssignment(region=region, sectionName=′grain′,offset=0.0,

offsetType=MIDDLE_SURFACE, offsetField=′′,

thicknessAssignment=FROM_SECTION)

2.2 網(wǎng)格劃分與邊界條件

對于簡單的幾何體，可以采用掃掠法。但是設計藥柱時，為內(nèi)彈道的平穩(wěn)性通常會有含有翼槽等幾何結(jié)構(gòu)，就不能直接生成全六面體網(wǎng)格，需要做適當?shù)钠史郑謩e對剖分后的每一部分幾何元素采用掃掠法生成六面體網(wǎng)格。

如翼柱形裝藥的網(wǎng)格劃分的難點集中在翼槽的位置，這里幾何形狀較為復雜，有眾多線段、倒角、圓弧等幾何，需要將這一段單獨剖分出來，與中間的圓柱段分離。如圖7所示，利用Face1可將翼槽單獨剖分出來，之后利用Face3沿著軸向剖分翼槽。需要說明的是，F(xiàn)ace3與翼槽的拉伸面平行，且距離翼槽中心的面距離為Thick/2-3，Thick為翼寬，如圖8所示。

圖7 翼柱形藥柱(1/2n)的幾何剖分

圖8 Face3示意圖

其中Face1和Face2可以通過創(chuàng)建模型時由Pro/E輸出面的x點坐標確定，F(xiàn)ace3由偏移翼槽剖面得到。三個面的創(chuàng)建相關代碼如下：

Face1=mdb.models[′Model-1′].parts[′grain′].DatumPlaneByPrincipalPlane(offset=offset1,principalPlane=YZPLANE)

Face2=mdb.models[′Model-1′].parts[′grain′].DatumPlaneByPrincipalPlane(offset=offset2,principalPlane=YZPLANE)

Face3=mdb.models[′Model-1′].parts[′grain′].DatumPlaneByOffset(flip=SIDE2,offset=(Wing_Thick*0.5-3),plane=plane1)

Offset1和Offset2由Pro/E輸出面的x坐標來定位，plane1為翼槽的剖面。

通過上述方法的剖分，翼槽的兩個幾何元素就可采用掃掠法生成六面體網(wǎng)格，掃掠路徑見圖9。

圖9 翼槽掃掠示意圖

本文采用的通星型裝藥網(wǎng)格劃分較為容易，從前端到后端即是掃掠路徑，藥柱的網(wǎng)格單元選擇雜交公式的熱力耦合C3D8HT單元。生成的網(wǎng)格模型如圖10所示。

圖10 藥柱翼槽有限元網(wǎng)格

其實現(xiàn)命令如下：

p=mdb.models[′Model-1′].parts[′GRAIN_XING′]

p.seedPart(size=seedSize,deviationFactor=0.1,minSizeFactor=0.1)

elemType1 = mesh.ElemType(elemCode=C3D8HT,elemLibrary=STANDARD)

c = p.cells

pickedRegions =(c,)

p.setElementType(regions=pickedRegions,elemTypes=(elemType1)

生成的星型藥柱網(wǎng)格如圖11所示。

圖11 星型裝藥網(wǎng)格

由于裝藥載荷和幾何均具有對稱性，計算通常只選取1/2n模型，為限制剛性位移，通常還需設置固支邊界條件，為識別這些表面，可采用RGB編碼來標識。其實現(xiàn)如下：

p1.SetFromColor("shell_duichen_set",(255,255,0))

p2.SetFromColor("jrc_duichen_set",(255,255,0))

p3.SetFromColor("grain_duichen_set",(255,255,0))

a=mdb.models[′Model-1′].rootAssembly

p1_in=a.instances[′XING_SHELL′]

p2_in=a.instances[′XING_JRC′]

p3_in=a.instances[′GRAIN_XING′]

a.SetByBoolean("duichen_set",(p1_in.sets[′shell_duichen_set′],

p2_in.sets[′jrc_duichen_set′],p3_in.sets[′grain_duichen_set′]))

由以上語句可知，對不同的表面進行顏色標識時，不同顏色的表面建立不同的集合，然后進行標識儲存。

3 快速計算平臺實現(xiàn)

3.1 界面設計思想

固體火箭發(fā)動機裝藥結(jié)構(gòu)完整性快速計算分析軟件包括了模型生成、材料參數(shù)輸入、仿真計算及后處理等模塊，如圖12所示。

圖12 溫度沖擊裝藥仿真分析平臺

(1)圖形用戶界面：用于完成用戶與程序間的數(shù)據(jù)傳遞，顯示仿真結(jié)果。

(2)模型生成：用于輸入模型幾何參數(shù)，生成發(fā)動機模型。

(3)材料參數(shù)輸入：用于輸入裝藥各部件材料參數(shù)。

(4)仿真及后處理：用于提交inp文件給ABAQUS以及統(tǒng)計結(jié)構(gòu)完整性仿真結(jié)果，包括應力、應變、位移。

軟件采用模塊化的方式建立了界面，軟件分為4個模塊，分別負責三維視圖顯示、連接Pro/E、處理Python腳本文件與界面可視化。采用模塊化設計的優(yōu)勢不言而喻，模塊之間相互獨立，一個模塊內(nèi)部發(fā)生變化時，其他模塊幾乎不受影響。

3.2 GUI界面調(diào)用

利用JAVA GUI接口建立軟件界面，實現(xiàn)軟件可視化。這部分主要負責連接Pro/E與ABAQUS，并且提供與用戶的交互功能。其工作流程大致如下：在后臺連接Pro/E，打開模型模板，讀取模板參數(shù)至界面供用戶修改，在用戶修改完畢后，重新生成模型，并轉(zhuǎn)換為ABAQUS支持的格式；讀取相應計算模塊的ABAQUS腳本參數(shù)，用戶修改后生成新的腳本；后臺執(zhí)行腳本，分析開始。

軟件共有3個大的功能界面。每個大的功能界面底下有5個獨立的小的功能，分別為模型幾何參數(shù)修改、材料參數(shù)設置、網(wǎng)格與計算參數(shù)設置、結(jié)果云圖查看、安全系數(shù)查看。

下文以材料參數(shù)設置和計算及后處理為例進行展示。

3.3 材料參數(shù)

設置材料參數(shù)是溫度沖擊仿真分析的關鍵步驟，精確完備的材料參數(shù)是完成高精度仿真分析的基礎，如圖13所示，該界面實現(xiàn)了藥柱、絕熱層、殼體的密度、比熱、熱傳導率等材料參數(shù)的設置，其中絕熱層與藥柱的粘彈特性用Prony級數(shù)與時溫等效因子來表征，簡明清晰，可以實現(xiàn)材料參數(shù)的快速設定和修改。

圖13 材料參數(shù)設置界面

3.4 計算及后處理

圖14為求解計算界面，該界面完成溫度沖擊裝藥仿真的計算及后處理部分，在計算完成后可以對溫度沖擊分析的固化降溫階段、高溫沖擊及低溫沖擊階段的溫度和Mises應力及應變，無需進行繁瑣的轉(zhuǎn)換，用戶可以自行選擇所關注部件進行2D或者3D結(jié)果顯示。

圖14 計算與后處理

3.5 模型優(yōu)化及收斂性分析

在建立通用建模方法時，應根據(jù)裝藥的形狀和尺寸，以從底向上的方法，從低維圖元向高維圖元建立藥柱幾何模型。保證藥柱模型中圖元的數(shù)量和相對位置的固定，保證建模過程中圖元的生成順序固定，確保圖元之間具有正確的拓撲約束關系。同時還需要對建模過程中用到的關鍵參數(shù)進行標記，固定其編號，避免因圖元編號混亂而導致建模意外終止。為后續(xù)通過圖元編號進行實體剖分和網(wǎng)格加密等步驟提供基礎。

裝藥模型優(yōu)化方法[16]框架中主要包含優(yōu)化算法更新參數(shù)、模型重生及導出、溫度沖擊裝藥結(jié)構(gòu)完整性計算。優(yōu)化算法和模型重生及導出模塊之間通過預留的參數(shù)配置文件作為數(shù)據(jù)交換的接口；模型重生與溫度沖擊裝藥計算模塊之間采用顏色標識的方法傳遞初始信息，事先要給表面標記上特定的顏色，裝藥結(jié)構(gòu)完整性計算運用模型中的顏色標識進行識別。

計算結(jié)果的收斂性主要受模型構(gòu)造復雜度、網(wǎng)格劃分精度、分析步類型、工況、邊界條件設置五個方面，提高仿真分析的收斂性主要從以下幾個方面入手：

(1)在進行幾何建模時，可忽略藥柱外表面的凸臺等對結(jié)構(gòu)完整性影響不大的細小幾何特征，將藥柱外表面看成光滑、連續(xù)的曲面，這樣以來就無需考慮和凸臺有關的幾何模型參數(shù)，簡化了建模過程。

(2)劃分成精度更高的全六面體網(wǎng)格；采用映射法或掃掠法等高質(zhì)量網(wǎng)格生成算法進行網(wǎng)格劃分，對于復雜造型的藥柱，需先分割成簡單的子實體；對重點關注區(qū)域進行適當?shù)木W(wǎng)格加密操作。

(3)依據(jù)實際工作狀態(tài)選擇恰當?shù)姆治霾筋愋?，設置適當大小的分析步長，根據(jù)理論預估情況確定幾何非線性的開關，可有效提高仿真分析的收斂性。

(4)定義工況時，確保工況參數(shù)的合理正確性，如溫沖過程要準確設定高低溫沖擊時間、裝藥轉(zhuǎn)移時間。初始溫度、對流換熱系數(shù)等，特別注意的是要保證各參數(shù)單位的一致性。

(5)邊界條件設置既要合理又要簡化，在進行裝藥溫度沖擊分析時，裝藥的邊界條件是不變的，即在裝藥柱外表面上施加固支條件，約束外表面節(jié)點前三個方向的自由度，在藥柱對稱面上施加對稱約束條件，約束對稱面節(jié)點第二個方向的自由度。另外需要根據(jù)具體工況添加必要的邊界條件，如溫度沖擊時需要添加溫度邊界條件。

通過以上五點可有效解決溫度沖擊裝藥仿真分析的收斂性問題。

3.6 快速性分析

在保證裝藥模型、材料參數(shù)、邊界條件、網(wǎng)格單元等一致的情況下，在已驗證手動計算結(jié)果相對正確的條件下，以該特定星型裝藥為例進行溫度沖擊裝藥結(jié)構(gòu)完整性的快速性對比分析，各階段用時對比結(jié)果如表5所示。

表5 原始手段與快速計算對比

在計算結(jié)果基本正確的前提下，從表5總的時間與手動計算的比對中可以看到，在計算精度一致的情況下，計算效率提升了54.03%。有效提高了溫度沖擊裝藥結(jié)構(gòu)完整性的仿真分析效率。

4 結(jié)論

(1)通過分析裝藥參數(shù)化建模的核心思想和實現(xiàn)途徑，基于J-Link的Pro/E二次開發(fā)技術，實現(xiàn)了藥柱的修改、再生、保存、導出等一系列功能，建立起設計變量與藥柱模型尺寸約束之間的映射關系，從而實現(xiàn)變量驅(qū)動藥柱更改的功能，與計算模板做到了有效的數(shù)據(jù)傳遞。

(2)本章利用ABAQUS二次開發(fā)技術，與參數(shù)化建模模塊形成了接口。建立了幾種典型藥型的結(jié)構(gòu)完整性計算模板，通過修改模板內(nèi)的參數(shù)，即可自動進行模型導入、網(wǎng)格離散、載荷及邊界條件添加、有限元計算等功能。

(3)完成了溫度沖擊裝藥結(jié)構(gòu)完整性的快速計算平臺設計，整合了Pro/E的參數(shù)化建模、溫度沖擊計算模板，采用圖形用戶界面實現(xiàn)簡明高效的溫度沖擊數(shù)值仿真，在保證計算精度情況下，最終實現(xiàn)溫度沖擊裝藥結(jié)構(gòu)完整性的快速計算，能極大提高工作效率。