螺旋桨流体仿真：SolidWorks参数化建模到推力曲线

简介：面向机械工程、航空航天与船舶推进方向的学习者和工程师，这份资源围绕螺旋桨流体仿真展开，结合SolidWorks环境下的CFD建模、网格划分与结果分析，适合用于课程设计、课题预研及设计验证参考。压缩包共105个文件，约92.06MB，包含sldprt/sldasm三维模型与装配文件、fld/res/plot/vort等仿真结果文件，以及txt日志、参数与配置数据，可据此复现计算流程并对照流场输出。已有772人学习下载，说明其在螺旋桨性能分析场景中具有一定参考价值。读者可从中了解推力、扭矩、阻力等指标的仿真组织方式，查看速度、压力分布等结果表达，并借鉴模型设置、边界条件与后处理思路，为优化叶片形状和工况分析提供线索。

1. 螺旋桨流体仿真：一份能从建模一路走到推力曲线的完整流程

手上有张桨叶型线图，要回答的问题很具体：直径 254 mm、转速 3000 rpm、来流 5 m/s 时，这桨能出多大推力、耗多少扭矩、效率落在哪个区间。没有风洞，也不能真去加工一付桨试。多数人第一反应是上通用 CFD，从画贴体网格开始折腾两周，最后卡在网格质量上。这套流程走的是另一条路：直接在 SolidWorks 里用全局变量驱动桨叶几何，再交给 SolidWorks Flow Simulation 做外流场求解，把建模、网格、求解、后处理串成一条可复现的链路。它适合两类人：一是课程设计或毕设要做螺旋桨性能分析的，二是需要在方案评审里拿出推力曲线、又不想把周期耗在网格工具链上的从业者。

2. 需求拆解与技术选型：三类交付指标与建模链路的职责划分

在动手建模之前，最重要的是把“做个螺旋桨仿真”翻译成几个能写进验收表的指标。需求阶段模糊，后面网格怎么划、算几个工况点、结果怎么呈现，全都会跟着飘。这一章把需求拆开，再解释为什么这条链路把几何交给 SolidWorks、把流场交给 Flow Simulation，以及稳态与瞬态之间该怎么选。

2.1 把“仿真一下螺旋桨”拆成三类可验收的指标

第一类是性能指标，也是最终要写进报告的核心：在 3 到 5 个进速比 J 上给出推力 T、扭矩 Q 和效率 η，并且要形成 J-Kt、J-Kq、J-η 三条曲线，而不是只算一个点。单点结果没法判断趋势，评审时一句“换个转速会怎样”就能问住。

第二类是流场指标，用来解释性能数据为什么长这样：桨叶吸力面的低压区分布、叶梢附近的涡结构位置、桨盘前后的轴向诱导速度。这部分不需要算得极准，但必须能自圆其说，比如叶梢低压区如果跑到弦长 70% 就明显不合理。

第三类是校核指标，用于证明结果可信：网格无关性是否收敛、同一工况下不同湍流模型差多少、与叶素-动量理论（BEMT）估算值的偏差范围。通常的做法是至少准备一套独立估算作为参照，哪怕精度不高，也能兜住“你的数是怎么来的”这个问题。

提示：三类指标要在开算前定好，尤其是工况点数量。算完再补工况，网格和求解设置往往要重来一遍。

2.2 几何、网格、求解三段的工具分工

SolidWorks 在这条链路里承担两件事：一是参数化几何，二是仿真的宿主环境。参数化是后面的关键——直径、弦长分布、扭转角、螺距比全部用全局变量和方程式驱动，改一个数，桨叶形状跟着变，仿真配置不用重搭。这一点在需要扫多个螺距比做对比时省下的时间非常可观。

Flow Simulation 的技术路线和通用 CFD 不一样：它用的是求解域内的正交六面体网格加部分单元处理，固体壁面不需要贴体网格，靠“最小壁厚”“最小缝隙”这类参数控制壁面附近的分辨率。好处是省掉了画贴体网格这一步，对课程设计这种周期紧、人手少的场景非常友好；代价是壁面附近的边界层分辨率不如专门的贴体网格精细，对分离强烈的工况要更谨慎地看结果。

这也是选型时最容易被误解的地方。不是 Flow Simulation “更简单所以精度差”，而是它把复杂度从网格工具转移到了物理设置上：网格级别、局部加密、湍流模型、旋转区域尺寸，这几项定错了，结果照样跑偏。翻车的往往是这一步，不是软件本身。

2.3 稳态 MRF 还是瞬态滑移网格：先想清楚要回答什么问题

旋转部件的处理有两条路。稳态的多重参考系（MRF）把旋转区域内的流动在旋转坐标系下求解，区域外仍是静止坐标系，交界面上做通量传递。它给出的是一段时间内的平均流场，算得快、收敛稳，适合出性能曲线。瞬态滑移网格则让旋转区域真正转起来，能捕捉叶梢涡的脱落和推力脉动，代价是计算时间翻几倍到十几倍。

课程设计和方案汇报的场景，我一般直接选稳态 MRF。理由是：要交付的是时均推力和效率，脉动量不是重点；而且瞬态算例对时间步长、界面插值的要求更高，一旦设置不当，收敛曲线会很乱，很容易让人怀疑是模型错了。选稳态的同时要在报告里写清楚局限——MRF 给的是时均值，不能用来分析叶片通过频率下的激励特性。这句话写上去，比硬着头皮算瞬态更让人信服。

如果确实需要瞬态结果，建议先用 MRF 把稳态场跑收敛，再以此为初始场切到瞬态，能少走很多弯路。

3. 螺旋桨几何建模：从翼型坐标到能计算的实体

几何是这条流程的地基。桨叶不是简单拉伸体，截面沿展向既有弦长收缩，又有扭转角变化，某些设计还带后掠和梢部修形。手工在草图里一个个画截面，改一次参数就得重画一遍，这是最常见的时间黑洞。更稳妥的做法是先用脚本生成截面坐标，再在 SolidWorks 里放样。

3.1 用全局变量驱动直径、弦长与扭转角

先把设计变量定成一张表，全部写进 SolidWorks 的全局变量里，后面建模、仿真、后处理都引用同一套值，避免“模型里是 254 mm、后处理里写成 250 mm”这种低级但致命的错误。

截面坐标用一个 NACA 四位数翼型生成器算，再按展向站位施加弦长缩放和扭转角旋转。扭转角按等螺距假设取 β = arctan(P / (2πr))，这是螺旋桨几何里最常用的一阶近似，够用且便于解释。

import numpy as np # ---- 与 SolidWorks 全局变量一一对应的设计参数 ---- D = 0.254 # 桨叶直径 m R_hub = 0.022 # 轮毂半径 m N_STATION = 9 # 展向截面数 PITCH_RATIO = 0.8 # 螺距比 P/D CHORD_ROOT = 0.18 # 根部弦长 / D CHORD_TIP = 0.09 # 梢部弦长 / D def naca4(m, p, t, n=61): """生成 NACA 四位数翼型闭合点列，前缘余弦加密。 m: 最大弯度, p: 最大弯度位置, t: 最大厚度, n: 单侧点数""" beta = np.linspace(0, np.pi, n) x = 0.5 * (1 - np.cos(beta)) # 前缘加密，尾缘不至于稀疏 yt = 5 * t * (0.2969 * np.sqrt(x) - 0.1260 * x - 0.3516 * x**2 + 0.2843 * x**3 - 0.1015 * x**4) yc = np.where(x < p, m / p**2 * (2 * p * x - x**2), m / (1 - p)**2 * ((1 - 2 * p) + 2 * p * x - x**2)) dyc = np.where(x < p, 2 * m / p**2 * (p - x), 2 * m / (1 - p)**2 * (p - x)) theta = np.arctan(dyc) xu, yu = x - yt * np.sin(theta), yc + yt * np.cos(theta) xl, yl = x + yt * np.sin(theta), yc - yt * np.cos(theta) # 从尾缘沿上表面到前缘，再沿下表面回到尾缘，保证闭合 return (np.concatenate([xu[::-1], xl[1:]]), np.concatenate([yu[::-1], yl[1:]])) def stations(): """按展向站位输出弦长与扭转角，曲柄为线性过渡，可按需换分布""" r = np.linspace(R_hub, D / 2, N_STATION) rn = (r - R_hub) / (D / 2 - R_hub) # 归一化 0~1 chord = D * (CHORD_ROOT * (1 - rn) + CHORD_TIP * rn) pitch = PITCH_RATIO * D beta = np.degrees(np.arctan(pitch / (2 * np.pi * r))) # 等螺距扭转角 return r, chord, beta if __name__ == "__main__": r, chord, beta = stations() xa, ya = naca4(0.04, 0.4, 0.12) # 4% 弯度、12% 厚度 for i, (ri, ci, bi) in enumerate(zip(r, chord, beta)): # 弦向按弦长缩放，绕 Z 轴按扭转角旋转，展向直接放到 y = ri xr = (xa - 0.3) * ci * np.cos(np.radians(bi)) - ya * ci * np.sin(np.radians(bi)) zr = (xa - 0.3) * ci * np.sin(np.radians(bi)) + ya * ci * np.cos(np.radians(bi)) np.savetxt(f"station_{i}.csv", np.column_stack([xr, np.full_like(xr, ri), zr]), delimiter=",", header="x,y,z", comments="")

代码分两段：naca4负责生成单个截面的二维型线，stations负责给出每个截面所在的半径、弦长和扭转角。弦长用根部到梢部的线性过渡，实际设计里也可能用椭圆或指数分布，改chord那一行即可。(xa - 0.3)这一项是把弦向坐标平移到以 30% 弦长为变距轴，这个位置接近多数螺旋桨的变距中心，放样出来的桨叶不会因为扭转而整体偏移。输出的 CSV 三列是 x、y、z，y 就是展向半径。

在 SolidWorks 里用“曲线通过 XYZ 点”逐条导入这些 CSV，得到 9 条截面曲线，再用放样凸台把它们连成实体，同时把前缘和后缘各连一条引导线，保证放样过程中截面不会扭曲。

参数上还有两个地方值得调：一是n_station，9 个截面对于光滑放样够用，但如果扭转角变化剧烈，建议加到 12 到 15 个，否则曲面会出现肉眼可见的折痕；二是翼型本身的弯度和厚度，naca4(0.04, 0.4, 0.12)是一组常见的低速桨叶参数，改这两个值会直接影响升阻特性和吸力面低压区的形状。

3.2 截面放样与根部过渡的处理顺序

放样顺序有讲究。正确的做法是先做桨叶主体，从根部截面一路放样到梢部，再单独做一个轮毂圆柱，最后用圆角或小放样把桨叶根部和轮毂连起来。很多人图省事把根部截面直接埋进轮毂里做布尔并集，几何上没问题，但会在交界处留下一个尖锐的内凹角。

这个内凹角在流场里就是一把刀——网格在这里会被强行压缩，局部速度容易出现异常峰值，后处理看到的桨叶根部压力云图会有一块莫名其妙的深色区域，排查半天才发现是几何尖角导致的。通常的做法是在根部加 R2 到 R3 的过渡圆角，让流动平滑过渡。

倒角的取舍也值得说一下。桨叶梢部不建议做大圆角，梢部形状对叶梢涡影响明显，圆角会把尖梢本来清晰的涡结构抹平，和后面对比实验数据时容易产生系统性偏差。梢部做 0.5 mm 左右的去锐边处理就够了。

放样完成后要检查曲面的连续性。用“曲率梳”看一下截面连接处，如果出现明显的折点，说明引导线控制不够，需要补一条或者调整截面的对应点顺序。截面点顺序不一致是放样扭曲的最常见原因，生成 CSV 时保证每条曲线的起点都在前缘同一侧，能避免这个问题。

3.3 进仿真前的模型体检清单

几何建完不能直接扔进仿真。下面这几项每次都要过一遍，漏一项后面都可能变成几小时的排查：

• 实体数量：整个模型必须是一个实体。曲面、开放壳体、多个不相连的实体都会让流体域识别失败。
• 最小壁厚：用“检查”工具量一下桨叶最薄处。如果小于网格最小壁厚设定，壁面附近的网格会被压扁甚至穿透。
• 干涉检查：桨叶和轮毂之间不能有体积重叠，重叠区在划分流体域时会产生内部面，直接导致计算域不连通。
• 坐标系：原点放在桨毂中心，桨轴与 Z 轴重合，旋转方向按右手定则。仿真里的旋转轴直接引用这个坐标系，偏一点就是偏心旋转。
• 单位：文档单位设成毫米还是米，后面所有仿真设置都跟着走。这里出过一次问题，模型用毫米建，仿真里默认米，推出来的力差了六个数量级，排查了很久才想起来看单位。

4. 流体域与前处理：计算域尺寸、旋转区域与网格的定档

前处理这几步最考验经验，因为它们没有唯一解，只有“够用”和“不够用”。计算域开小了，尾流没充分发展就撞到出口边界，推力会被边界反射干扰；旋转区域贴太紧，桨尖和静止区域的交界面会出问题。这一章给出我常用的取值区间和判断方法。

4.1 计算域：外流场的三向尺寸怎么定

外流场的计算域以桨盘为基准。轴向方向上，入口放在桨盘前方 5D 处，出口放在后方 10D 处，这是螺旋桨数值模拟里比较通行的比例。入口选 5D 是因为来流在前方基本均匀，不需要更长的缓冲区；出口要 10D，是因为尾流在后方的收缩和旋涡耗散需要足够距离，出口太近会让压力边界条件“看见”桨盘，反过来影响推力。

径向取 5D，也就是计算域是一个边长 5D 到 10D 的长方体。无人机的螺旋桨类应用里，桨盘附近的诱导流场影响范围大致在这个尺度内。

注意：计算域一旦缩小，收敛曲线会变得很难看，而且是那种“看起来收敛了但数值不对”的假象。宁可算久一点，别在域尺寸上省。

4.2 旋转区域与静止区域的分界面

旋转区域是一个包住桨叶的圆柱，直径取 1.1D，厚度取 0.35D。直径取 1.1D 是为了让桨尖和旋转区域的静止交界面之间留出约 0.05D 的间隙——如果交界面正好落在桨尖上，叶梢的高速度梯度会被交界面的插值抹平，叶梢涡直接消失。

厚度 0.35D 是从桨叶前缘最前点到后缘最后点再各留出一点余量。太薄会导致桨叶部分区域落在旋转区外，那部分就不转了，结果会明显偏低；太厚则旋转区域的体积变大，交界面上要插值的网格变多，算得慢，还容易引入不必要的数值耗散。

设置时在 Flow Simulation 里选择“局部旋转区域”，选中包住桨叶的圆柱实体（如果建模时没画，可以先在仿真里手动创建），转速填 3000 rpm，方向按右手定则确认一遍。

4.3 全局网格级别与桨叶表面局部加密

Flow Simulation 的网格是正交六面体，基础分辨率由“全局网格级别”控制，通常取 3 级起步，配合局部网格对关键区域加密。全局级别 3 意味着基础网格在最短方向上有 8 个单元，对于 254 mm 的桨叶太粗，所以必须做局部加密。

局部网格分两块：一是旋转区域整体，级别提到 5 到 6；二是桨叶表面，用“局部网格”功能按面选，控制“最小壁厚”和“最小缝隙”两个参数。最小壁厚建议取桨叶最薄处的 1/2，最小缝隙取 1 mm 左右，让壁面附近至少有 3 到 4 层单元。

初始化之后要做一次“网格预览”，重点看两处：桨尖附近有没有出现被压扁的狭长单元，根部圆角处有没有单元穿透壁面。这两处出问题，后面结果一定异常。

4.4 边界条件、初始条件与目标设定

边界条件的默认值大部分可以直接用，但有两处必须改：

• 入口：从默认的“压力开口”改成“速度入口”，速度填 5 m/s，方向沿 +Z。用速度入口是因为桨的进速是已知量，压力入口会变成求解的一部分，增加不确定性。
• 出口：保持“压力开口”，静压设为环境压力 0 Pa（相对值）。
• 外部壁面：设为“理想壁面”，模拟无界外流。如果计算域尺寸足够大，这里用滑移壁面和理想壁面差别不大。

初始条件里要把“旋转”选项打开，转速填与旋转区域一致的 3000 rpm，否则求解器可能把它当成静止区域。

目标（Goals）设三个就够：桨叶表面的轴向力（推力）、绕旋转轴的扭矩、以及一个体积目标用于监控质量守恒。目标的收敛是判断收敛状态最直接的依据。

5. 求解收敛、后处理与避坑排查：六类翻车现场的处理路径

算得对不对，一半看几何和网格，一半看你会不会读收敛曲线和后处理数据。这一章先把收敛判据和结果提取讲清楚，再把项目里最容易出问题的几种情况摊开说。

5.1 收敛判据：残差、目标值与行程

Flow Simulation 里判断收敛不能只看残差降到多少。正交网格加部分单元的处理方式，残差曲线天然会比贴体网格更“毛躁”，残差降到一个量级就不太动了是正常现象。

真正的判据是两个：一是监控的目标值（推力、扭矩）在最后 20% 的迭代里波动小于 1%；二是“行程”达到计算域轴向长度的 3 倍以上，也就是流体从入口到出口走完三个来回。行程不够意味着流场还没充分发展，此时目标值即使看起来平了也是假的。

提示：如果目标值在最后阶段仍然缓慢单调下降，不要当作收敛。这通常说明旋转区域或网格还不够，再多迭代也收不住。

5.2 推力、扭矩、效率的提取口径与换算

Flow Simulation 的“表面目标”直接给出的是 X、Y、Z 三个方向的力和力矩，需要手动换算成无量纲系数。这里最容易出单位错误，所以我把换算写成一个脚本，每次算完直接把导出的 CSV 喂进去。

import numpy as np # ---- 与模型、工况严格一致的物理量 ---- RHO = 1.225 # 空气密度 kg/m3，按工况高度查表 RPM = 3000 # 转速 N = RPM / 60 # 转每秒 rev/s，注意不是 rad/s D = 0.254 # 直径 m V = 5.0 # 来流速度 m/s # 从 Flow Simulation 导出的表面目标，取稳定段最后 50 个迭代点 T = np.loadtxt("goal_force_z.csv", delimiter=",", skiprows=1, usecols=1) # N Q = np.loadtxt("goal_torque_z.csv", delimiter=",", skiprows=1, usecols=1) # N·m T, Q = T[-50:].mean(), Q[-50:].mean() # 无量纲化：Kt = T/(ρn²D⁴), Kq = Q/(ρn²D⁵), J = V/(nD) Kt = T / (RHO * N**2 * D**4) Kq = Q / (RHO * N**2 * D**5) J = V / (N * D) eta = J * Kt / (2 * np.pi * Kq) # 效率，量纲一 print(f"J={J:.3f} Kt={Kt:.4f} Kq={Kq:.5f} eta={eta:.3f}") print(f"T={T:.3f} N Q={Q:.4f} N·m")

几个参数必须核对：转速要用转每秒（N = RPM/60），不是角速度。往 Kt 公式里代角速度是初学者最容易犯的错，结果会差 (2π)² 倍。直径用米，且是完整直径而非半径。密度按实际工况取，常温常压取 1.225，如果做的是高空工况要换。

稳定段取最后 50 个点做平均，不要只取最后一个点。最后一个点可能正好落在振荡的波峰上，均值更稳。

量级参考：直径 254 mm、转速 3000 rpm、进速 5 m/s 这个工况下，推力大约在 1 到 2 N，扭矩在 0.05 N·m 左右，效率在 0.4 到 0.5 之间。如果你算出来推力是 140 N 或者 0.014 N，先别怀疑物理模型，回去查单位。

5.3 网格无关性：三套网格怎么收敛到可信区间

网格无关性不是“越密越好”，而是证明结果对网格不再敏感。做法是固定其他所有设置，只改全局网格级别和局部加密等级，跑三套：粗（全局 3、局部 4）、中（全局 4、局部 5）、细（全局 5、局部 6）。

判据是相邻两套之间的推力差异小于 3%，扭矩差异小于 5%。如果粗到中差 15%、中到细差 2%，说明中等网格已经够用，报告里就写中等网格的结果，并附上这张对比表。这张表比任何文字说明都有说服力。

5.4 翻车点一：推力随迭代持续振荡，怎么调都不收敛

现象是目标值曲线在某个区间里上下跳，幅度超过 5%，迭代到几千步也不见收窄。

原因通常有两个。一是旋转区域直径太接近桨尖，交界面的插值把叶梢的高梯度抹来抹去，每次迭代结果都不一样。二是局部网格的最小缝隙设得太大，壁面附近单元太粗，分离点和再附点位置在网格尺度上跳变。

解决顺序是：先把旋转区域直径从 1.05D 放宽到 1.1D 到 1.15D，重算；如果还振荡，把桨叶表面的最小缝隙减半，重做一次网格预览确认壁面附近至少有 3 层单元。这两步做完还没收敛，再去看是不是出口太近导致压力边界反射，把出口从 8D 拉到 12D 试试。

5.5 翻车点二：推力数量级明显不对

现象是算出来的推力比预期大或小几十倍甚至上千倍，但流场云图看起来“挺正常”。

原因几乎都出在单位上。常见的三种：模型用毫米建、仿真里按米算，长度差 1000 倍，力差 10⁶ 倍；转速填了 3000 当成转每秒，实际应该是 50；直径填了半径 127 mm。

解决办法是建一张量纲核对表，把长度、速度、转速、密度、力、力矩六项的单位和取值列出来，每次开算前对一遍。这一步花两分钟，能省下几小时。

5.6 翻车点三：桨叶表面出现不明原因的异常压力区

现象是压力云图在桨叶根部或某个截面上出现一块与周围明显不连续的深色或亮色区域，形状规整得不自然。

原因是几何内部有问题。最常见的是放样时截面点顺序不一致，导致曲面自交；其次是桨叶和轮毂布尔并集时留下内部面，流体域识别时把这块当成了壁面。还有一种隐蔽情况是模型里有极小的一条边或一个面，肉眼看不见，但网格生成时会在这里强行加密。

解决靠三个检查：曲面曲率梳看有没有折点，实体检查看实体数量是否唯一，“检查”工具里的最小边和最小面。发现异常就回去改几何，别指望靠加密网格绕过去。

5.7 翻车点四：改一个参数就要全部重来

现象是导师说“把螺距比从 0.8 改到 1.0 再看看”，结果几何要重画、边界条件要重设、目标要重选，一下午没了。

原因是最初建模时所有尺寸都是硬编码，没有用全局变量和方程式驱动。SolidWorks 的方程式功能可以做这件事：把直径、弦长、扭转角写成全局变量，草图尺寸和放样引导线都引用变量名。改一个值，整个模型自动重建。

仿真端也一样，把工况参数（转速、来流速度）单独建一个配置，用派生配置管理多个工况点。这样扫 5 个进速比不需要手动改 5 次，配置切过去就行。

6. 验收与答辩：怎么证明这套仿真结果站得住

结果算出来只是开始，真正决定评价的是你怎么呈现它。评审最常问的三个问题，其实都有标准答案的结构：网格无关吗、湍流模型为什么选这个、结果可信度多少。提前把答案组织成表格和曲线，比临场解释有效得多。

第一张表是网格无关性对比表，就是 5.3 那张。三套网格、推力扭矩、相对偏差，一眼就能看出结果已经收敛。这张表放在报告前部，评审翻到就知道你做过校核。

第二张表是多工况性能表，把 3 到 5 个进速比下的 J、Kt、Kq、η 列出来，配上三条曲线。曲线的形状比单点数值更能说明问题：Kt 随 J 增大而单调下降、η 在中段出现峰值，这两条趋势符合螺旋桨的一般规律，如果曲线形状反常，说明某个工况点算错了。

第三张表是与独立估算的对比。用叶素-动量理论把同一付桨的推力和效率估一遍，和仿真结果并排放。偏差在 10% 到 15% 以内是正常的，因为 BEMT 本身忽略了很多三维效应；如果偏差到 30% 以上，要么 BEMT 的升阻数据不对，要么仿真设置有问题，两种都要查。

答辩时被问到湍流模型，不要背定义。回答的思路是：外流场以附着流为主、分离区集中在根部，用 k-ε 的稳健性好、收敛快；如果关心叶梢涡的精细结构，会切到 k-ω SST。然后补一句“我用两种模型各算了一个工况点，推力差异在 4% 以内”，这句话比任何理论解释都有分量。

最后还有一个习惯值得养成：每次交付前，把原始几何、仿真文件、导出数据、后处理脚本一起归档，文件名带上日期和工况编号。有次改了一个全局变量忘了存副本，结果第二天想复现前一天的曲线，怎么都跑不出同样的数，只能重算一遍。从那以后我每次改参数前都强制走一遍“另存为 + 记录变量值”的流程，哪怕只是改一个弦长。这套流程连同建模参数表、检查清单和后处理脚本，我整理进了一份归档里，需要照着做的可以直接拿去改。

希望帮到你。

本文还有配套的精品资源，点击获取