• 2020年第35卷第3期文章目次
    全 选
    显示方式: |
    • >燃烧、传热、传质
    • 某型航空活塞发动机调贫阈值的实验测定

      2020, 35(3):449-456. DOI: 10.13224/j.cnki.jasp.2020.03.001

      摘要 (976) HTML (0) PDF 2.23 M (735) 评论 (0) 收藏

      摘要:基于某型航空活塞发动机在运行中因过度调贫而导致大量排气门烧蚀的现象,为控制调贫幅度获取便于工程实践操作的安全调贫阈值,通过测量排气门密封面附近气门体的实际运行温度,获取了该型发动机上气门运行温度与排气温度之间的对应关系,解决了以排气温度(EGT) 为参变量间接衡量气门体实际运行温度的工程问题。结果表明:利用该对应关系在确保排气门密封面上保护膜不被破坏的前提下得出该型发动机上安全调贫阈值为排气温度为732℃,经70000飞行小时的实践验证表明此阈值可靠。

    • 加力用扇形喷嘴雾化特性试验

      2020, 35(3):457-470. DOI: 10.13224/j.cnki.jasp.2020.03.002

      摘要 (500) HTML (0) PDF 6.36 M (694) 评论 (0) 收藏

      摘要:根据加力燃烧室内锥凹腔点火与联焰要求,设计了扇形喷嘴并开展相应的雾化试验,研究了供油压差、扇形角度及扇形出口高度等参数对流量特性和雾化特性的影响以及加力环境下横向气流的温度、速度和供油压差对索太尔平均直径(SMD)及穿透深度的影响。采用称质量法测量流量系数,利用马尔文粒度仪和高速摄影仪对下游SMD、雾化角度及穿透深度进行测量。结果表明:①供油压差增大,流量系数先减少,后稳定;②供油压差一定,扇形出口角度越大,流量系数和雾化角度也越大;③扇形出口高度增加,雾化效果变好;④出口位置对雾化特性影响不大;⑤供油压差越大,穿透深度越大,SMD减小;⑥横向气流速度越大、温度越高,穿透深度越浅,油雾场越靠近下游;⑦横向气流温度越高, SMD越小。

    • 航空活塞二冲程发动机的可变排气阀换气过程

      2020, 35(3):471-481. DOI: 10.13224/j.cnki.jasp.2020.03.003

      摘要 (518) HTML (0) PDF 4.50 M (687) 评论 (0) 收藏

      摘要:航空活塞二冲程发动机的换气过程直接影响发动机的性能,采用可变排气阀可以改善发动机的动力性能。通过GT-Power软件搭建发动机一维整机模型,进行可变排气阀的仿真研究,分析了发动机在不同转速、不同排气阀状态下的换气过程及排气阀状态对换气过程的影响。在此基础上,还对定速巡航工况下(6500r/min转速、50%节气门开度)的排气阀状态进行了优化。结果表明:不同的排气阀状态可以改变发动机排气过程压力波动的相位,从而影响发动机的换气品质及燃烧过程;在定速巡航工况下,排气阀为中等开度时(初级阀关闭、次级阀开启),其发动机动力性能达到最优。

    • 平板式预膜喷嘴初次雾化特性试验

      2020, 35(3):482-492. DOI: 10.13224/j.cnki.jasp.2020.03.004

      摘要 (592) HTML (0) PDF 5.31 M (764) 评论 (0) 收藏

      摘要:为深入了解平板式预膜喷嘴初次雾化特性,试验研究不同进口条件对平板式预膜喷嘴初次雾化特性的影响规律。试验采用了背光照明和片状激光照明相结合的高频拍摄手段分别获得液膜俯视和侧视破碎形态,同时运用本征正交分解(POD)法和液膜边缘定位等分析方法进行光学图像结果后处理,获得表征平板式预膜喷嘴初次雾化特性三个物理量:液膜波动频率、破碎距离和横向不稳定波长。试验结果表明:①通过分析侧视和俯视破碎形态,平板式预膜喷嘴液膜破碎形态可分为三类:末端破碎、波浪脱落和表面剥离,进口韦伯数对预膜喷嘴破碎形态的影响占主导地位;②把POD法和液膜边缘定位方法等相结合方法应用到高频非接触光学喷嘴雾化图像的后处理分析中,是一种非常有效的数据后处理方法;③液膜初次雾化特性主要受到进口韦伯数和气液动量比的影响,破碎距离和横向不稳定波长都随进口韦伯数的增加而降低,液膜表面不稳定波动频率随进口韦伯数的增加而增加,所获得的经验关系式与试验数据吻合较好。所获得的进口参数对液膜破碎形态和雾化特性的影响规律为喷嘴后续优化设计提供了依据。

    • 主流攻角对双射流孔气膜冷却特性的影响

      2020, 35(3):493-501. DOI: 10.13224/j.cnki.jasp.2020.03.005

      摘要 (648) HTML (0) PDF 5.07 M (740) 评论 (0) 收藏

      摘要:应用压力敏感漆技术,在平板上测量了不同主流攻角(i=-30°,-20°,-10°,0°,10°,20°,30°)下双射流孔的气膜冷却效率,并利用计算流体动力学(CFD)计算得到的流场对气膜冷却效率的规律进行了解析。所研究的双射流孔结构的孔间无量纲横向距离为0.5,孔间无量纲流向距离为3;射流与主流密度比为1.0,吹风比分别为0.5、1.0、1.5、2.0。结果表明小的主流攻角(i=-10°,0°,10°)下,流场中存在反肾型涡对或挤压作用,气膜层与壁面贴附良好,气膜冷却效率最高;大正值攻角(i=20°,30°)下,虽然气膜覆盖面积大,但反肾型涡对退化,气膜冷却效率下降;大负值攻角(i=-20°,-30°)下,流场中有肾型涡对,且气膜横向覆盖受限,气膜冷却效率最低。

    • 盘腔进气位置对径向预旋系统的影响

      2020, 35(3):502-509. DOI: 10.13224/j.cnki.jasp.2020.03.006

      摘要 (1650) HTML (0) PDF 3.16 M (554) 评论 (0) 收藏

      摘要:为减少径向预旋系统的流动损失,运用数值模拟方法对不同盘腔进气位置的径向预旋系统进行分析,结果表明:随着盘腔进气径向位置的增加,预旋喷嘴出口气流旋流比随之逐渐减小,径向预旋系统的温降系数及总压损失系数均随之逐渐增大。当旋转雷诺数等于7.9×106,盘腔进气位置由低位向高位变化时温降系数最大可增加525%,同时总压损失系数增加3.93%。径向预旋系统内比熵增主要发生在预旋喷嘴和共转腔,约占系统总体比熵增的80%。随着盘腔进气径向位置的增加,径向预旋系统总体比熵增降低,预旋喷嘴比熵增占比逐渐增大,共转腔比熵增占比逐渐减小。

    • >叶轮机械
    • 某压气机静叶基元叶型优化设计及实验

      2020, 35(3):510-518. DOI: 10.13224/j.cnki.jasp.2020.03.007

      摘要 (542) HTML (0) PDF 3.45 M (749) 评论 (0) 收藏

      摘要:为了得到更加适合压气机静叶的叶型以降低气动损失,提取了静叶中径处的叶型,通过平面叶栅实验获得了原叶型的损失特性,发现原叶型气动损失较高,需要通过合理匹配设计参数来降低损失。为此,搭建数值优化平台在约束空间内搜寻气动损失更低的叶型,目标函数的构建综合考虑了多个冲角下的总压损失系数以提升叶片的变工况性能。优化结果显示:目标函数值降低了约9%,进一步实验研究发现,在实验涉及的整个马赫数和冲角范围内优化叶型比原叶型具有更低的总压损失系数,设计工况总压损失系数较原型叶型下降了31.3%,提升了叶型在正冲角边界附近的抗失速能力,设计进口马赫数正4°冲角下气流折转角增加1°。通过对实验结果的深入分析,解释了叶型性能提升的机理,对工作在相似环境的叶型设计及多目标优化方向给出了建议。

    • 涡轮叶片来流角扰动不确定性量化分析

      2020, 35(3):519-531. DOI: 10.13224/j.cnki.jasp.2020.03.008

      摘要 (428) HTML (0) PDF 5.24 M (713) 评论 (0) 收藏

      摘要:为了高效精确量化来流条件变化对涡轮叶片性能的影响,发展了一种基于自适应抽样技术的非嵌入式多项式混沌(NIPC)模型,并通过函数试验对自适应NIPC模型进行了验证。在不同尺度来流角变化下,用该模型量化涡轮叶栅能量损失系数的不确定性变化,同时与直接蒙特卡洛模拟和灵敏度分析进行对比。通过流场统计分析揭示来流角变化对气动不确定性的影响机理。结果表明:自适应NIPC在不同尺度来流角扰动下均能快速准确量化能量损失系数变化;能量损失系数不确定性变化对来流角扰动是非线性依赖关系;叶片通道内的激波对来流角变化最为敏感,是引起能量损失系数不确定性变化的主要因素。

    • 非均匀间隙下进气畸变对风扇性能的影响

      2020, 35(3):532-539. DOI: 10.13224/j.cnki.jasp.2020.03.009

      摘要 (560) HTML (0) PDF 5.48 M (643) 评论 (0) 收藏

      摘要:引入真实进气道的进气总压畸变,分析了偏心形式周向非均匀间隙下,畸变区相对位置对风扇流场的影响。结果表明,畸变流的影响总滞后于进口畸变位置,滞后角约为90°,在原有非均匀间隙流场中最恶劣的流动位置附加畸变流,对风扇性能的影响是最大的,且旋转失速与周向间隙收敛扇区的流动密切相关,推断一旦从大间隙到小间隙过程中无法将叶尖高熵溢流耗散,风扇转子将进入失速。

    • 特定涡旋流畸变对跨声速压气机性能的影响

      2020, 35(3):540-551. DOI: 10.13224/j.cnki.jasp.2020.03.010

      摘要 (459) HTML (0) PDF 7.47 M (663) 评论 (0) 收藏

      摘要:为研究S弯进气道出口和翼身融合体(HWB)飞机发动机进口两种特定涡旋流对压气机气动性能和稳定性的影响,将设计的新型叶片式旋流畸变发生器与跨声速压气机Stage 67进行联合数值仿真,获得进口旋流条件下的压气机特性线和流场分布,并与转子进口为均匀来流时进行对比分析。结果表明:在100%换算转速,S弯进气道产生的旋流使得压气机压比和效率在近峰值效率点分别下降0.06%、0.85%,稳定工作流量范围减小5.97%,压气机稳定裕度降低1.13%,压比和效率特性线均向左下方移动。HWB飞机产生的旋流使得压气机在近峰值效率点的压比增大6.61%,效率下降6.25%,稳定工作流量范围减小26.88%,压气机稳定裕度降低2.84%,压比特性线向右上方移动,效率特性线向右下方移动。

    • >结构、强度、振动
    • 考虑离心力与热应力非连续界面的微动磨损

      2020, 35(3):552-559. DOI: 10.13224/j.cnki.jasp.2020.03.011

      摘要 (507) HTML (0) PDF 2.68 M (675) 评论 (0) 收藏

      摘要:航空发动机、重型燃气轮机等旋转机械的组合转子轮盘间多采用非连续界面连接,通过建立两轮盘的环形平面接触有限元模型,并进行冷态预紧、升速工况下的结构分析及稳态温度场下的热-结构耦合分析。研究了环形接触界面的接触应力、切应力、磨损深度等微动参量的演变规律。结果表明:预紧力单独作用下接触应力从内径到外径逐渐减小,内径处接触应力和切应力最大,磨损最严重;随转速升高,离心力增大,接触界面的最大接触应力增大,最小接触应力减小,造成接触应力分布不均匀程度及磨损加剧;稳态温度场导致盘间接触应力均值大幅度上升,严重加剧轮盘间磨损;从盘间接触应力分布均匀性与界面微动磨损情况来说,热载荷是造成组合转子非连续界面微动磨损的主要因素。

    • 极低压下两相流挤压油膜阻尼器油膜参数特性

      2020, 35(3):560-568. DOI: 10.13224/j.cnki.jasp.2020.03.012

      摘要 (624) HTML (0) PDF 3.11 M (725) 评论 (0) 收藏

      摘要:针对极低供油压力工况开展实验研究,以考察气液两相流对挤压油膜阻尼器(SFD)油膜参数特性的影响。结果表明:当SFD入口气体体积分数小于0.9时,油膜阻尼随着入口气体体积分数的增加而减小,直至气体体积分数增大到0.9时,油膜阻尼是纯油状态时阻尼的60%;当气体体积分数大于0.9时,油膜阻尼大幅减小至几乎可以忽略不计。现有理论模型并不适用于极低供油压力工况。基于实验结果,找到了最符合SFD两相流动的等效黏度模型,其理论预测的油膜阻尼与不同供油压力下的实验数据吻合较好,为SFD两相流研究的模型选择提供了依据。

    • 基于非接触式测量的旋转叶片动应变重构方法

      2020, 35(3):569-580. DOI: 10.13224/j.cnki.jasp.2020.03.013

      摘要 (517) HTML (0) PDF 3.40 M (789) 评论 (0) 收藏

      摘要:基于叶端定时非接触式测量和振动响应传递比的概念,开展高速旋转叶片动应变重构方法的研究。在频域内推导了叶片任意测点位移与任意测点动应变的传递比,给出了单模态共振下响应传递比关于位移和应变模态振型的解析表达式;建立旋转叶片的三维(3D)有限元模型,开展考虑旋转预应力效应的叶片模态分析,提取位移和应变模态振型,获得任意转速下叶端位移与叶根关键点动应变的传递比。开展高速旋转叶片叶端定时非接触式测量实验,采用周向傅里叶算法对叶端定时信号进行处理,获得叶片在不同转速单模态共振下的叶端位移,结合响应传递比,重构5个旋转叶片的关键点动应变。结果表明:旋转叶片在9000r/min和13000r/min转速下发生1阶共振时,与应变片实测结果相比,叶根处应力最大点、次大点和边缘点3个关键点的动应变平均重构误差均小于15%,验证了旋转叶片动应变重构方法的有效性。

    • >火箭发动机
    • Experiment of gas plasma plume deflection by magnetic control

      2020, 35(3):581-588. DOI: 10.13224/j.cnki.jasp.2020.03.014

      摘要 (370) HTML (0) PDF 2.90 M (502) 评论 (0) 收藏

      摘要:To investigate the deflection effect of gas plasma plume controlled by magnetic field, a novel experimental scheme was presented. The Cs2CO3 catalytic ionization seeds were injected into the combustion chamber to obtain gas plasma on a high temperature magneto hydrodynamic (MHD) experiment rig. The plasma jet was deflected under the action of an external magnetic field, resulting in a thrust-vector effect. Particle image velocimetry (PIV) collected two-dimensional images of jet flow field. Through image processing and velocity vector analysis, the jet deflection angle can be obtained quantitatively. At 1800-2500K, the jet deflection was verified experimentally under the condition of 0.45T magnetic field strength. The results indicate that the jet deflection angle increases gradually with the increase of gas temperature, and above 2200K, the jet deflection angle increase obviously. In the process of gas plasma jet, it is feasible to realize the jet deflection controlled by MHD by adding an external magnetic field.

    • 导流槽构型对火箭起飞时流场和声场的影响

      2020, 35(3):589-596. DOI: 10.13224/j.cnki.jasp.2020.03.015

      摘要 (458) HTML (0) PDF 2.92 M (664) 评论 (0) 收藏

      摘要:为了研究四喷管运载火箭起飞时火箭周围噪声环境问题,建立燃气/空气双组分的可压缩流动模型,采用2阶Roe格式、SAS(scale-adaptive simulation)湍流模型和声学类比积分法Ffowcs-Williams Hawkings(FW-H)求解三维Navier-Stokes方程。以单机火箭的噪声问题为对象开展数值模拟,并将噪声数值计算结果与试验数据对比,误差在3dB以内(相对误差小于1.6%),验证预测噪声方法的有效性,进而研究四喷管运载火箭起飞阶段采用不同导流槽构型对箭体舱段区域噪声环境的影响,结果表明:在相同噪声接收点处,单侧导流槽对应的总声压级比双侧导流槽大,两者的总声压级(OASPL)之差最大为10.7dB。另外,当采用单侧导流槽时,沿周向接收点的噪声总声压相对单侧导流槽中心截面呈对称分布,而且沿导流出口方向逐渐增大。所建立的噪声数值方法为大推力捆绑运载火箭舱段的噪声环境预测及其控制提供一定参考。

    • 面向一体化计算的整体式液体冲压发动机建模及性能评估

      2020, 35(3):597-610. DOI: 10.13224/j.cnki.jasp.2020.03.016

      摘要 (428) HTML (0) PDF 5.90 M (646) 评论 (0) 收藏

      摘要:为发展冲压发动机性能工程预估方法,建立了面向一体化计算的整体式液体冲压发动机性能计算模型,并使其兼容基团贡献算法,提高了拓展性。利用模型分析了冲压发动机在攻角0°~6°、高度0~18km、马赫数2.0~3.5、余气系数1.0~2.9范围工况多维度连续变化下,比冲、推力系数和燃油质量流量的速度-高度特性、高度-节流特性、节流-速度特性以及攻角特性。研究结果表明:性能计算模型可在一体化计算条件下、基团贡献算法允许范围内,不依赖试验数据对冲压发动机性能预估,计算结果与技术参考值相比相对误差均小于14%。推力系数和比冲具有基本一致的速度-高度特性和相似的高度-节流特性,受燃气组分影响,推力系数和比冲的节流-速度特性差异明显。燃油质量流量的变化规律不同于比冲和推力系数,在进入平流层后呈现折缓趋势,而呈现连续性,速度越大、高度越低、余气系数越小,燃油质量流量越高,反之则越小。引入攻角以后,攻角越大,比冲和推力系数越低,进气道起动马赫数越高;攻角-起动马赫数曲线小范围内近似线性,攻角超过5.6°非线性加剧。

    • >气动热力学与总体设计
    • 螺旋桨安装效应对无人机气动特性影响

      2020, 35(3):611-619. DOI: 10.13224/j.cnki.jasp.2020.03.017

      摘要 (483) HTML (0) PDF 4.06 M (687) 评论 (0) 收藏

      摘要:为解决某型飞翼布局无人机(UAV)带动力构型风洞试验最大升阻比相对无动力状态大幅下降的问题,采用计算流体动力学(CFD)方法对无人机无动力与带动力构型进行了数值模拟,数值模拟结果分别与无动力以及带动力风洞试验数据吻合良好,在此基础上深入研究了螺旋桨安装效应对无人机气动特性的影响。结果表明:推力螺旋桨与机身之间气动干扰产生的低压区致使阻力增加,从而导致飞机最大升阻比相比无动力状态下降了30.7%。针对无人机在推力螺旋桨影响下出现的最大升阻比下降问题,采用增大螺旋桨与机身之间距离的方法可以有效地消除机身后部出现的低压区,减小了阻力,提升了无人机最大升阻比。桨毂拉长方案在8°和9°迎角下最大升阻比分别提升了17.3%和15.4%。

    • Influence of blade negative twist on aerodynamic performance of cycloidal propeller

      2020, 35(3):620-632. DOI: 10.13224/j.cnki.jasp.2020.03.018

      摘要 (448) HTML (0) PDF 7.65 M (592) 评论 (0) 收藏

      摘要:In order to improve the analysis accuracy of the influence of shape parameters on the aerodynamic performance of cycloidal propeller, a set of CFD methods suitable for the simulation of the cycloidal propeller flow field based on the dynamic overset mesh and Unsteady Reynolds time average method equation were established. The characteristics of the spanwise lift distribution and surface pressure distribution in case of different negative twist were provided. The influence of the negative twist on the aerodynamic performance of cycloidal propeller was emphasized. The results showed that, the change of negative twist caused the change of azimuth for cycloidal propeller lift. The larger the absolute value of negative twist, the lower the figure of merit of the cycloidal propeller; the smaller the negative twist, the more gentle the instantaneous aerodynamic fluctuation; the influence mechanism of negative twist on the aerodynamic performance of cycloidal propeller was generated from changing the non-uniform distribution of the cycloidal propeller's effective angle of

    • 基于实战使用的涡轴发动机空中起动飞行试验

      2020, 35(3):633-640. DOI: 10.13224/j.cnki.jasp.2020.03.019

      摘要 (441) HTML (0) PDF 3.60 M (727) 评论 (0) 收藏

      摘要:以涡轴发动机实战使用中空中紧急起动相关技术指标和特情处置措施验证为目的,设计了贴近实战使用特征的发动机空中起动试飞方案,并以某涡轴发动机设计定型试飞为依托,进行了不同飞行工况下的飞行试验验证和数据分析研究。结果表明:该空中起动试飞方案合理可行,能够最大程度的贴近部队实战使用的技术特征,满足对发动机起动高度包线、起动时间等研制总要求规定的技术指标验证需求;试飞结果表明该型涡轴发动机4000m以上高度停车至再起动期间直升机平均下滑高度约为670m,平均停车时间约为185s。该试飞方法和数据结果可为部队实战使用过程中空中遭遇停车的紧急处置提供支持和参考。

    • >自动控制
    • 某型涡扇发动机测量数据融合模型

      2020, 35(3):641-650. DOI: 10.13224/j.cnki.jasp.2020.03.020

      摘要 (420) HTML (0) PDF 2.13 M (530) 评论 (0) 收藏

      摘要:为提高传感器测量数据的有效性、可信性,基于共同工作方程,采用原始对偶内点法对某型涡扇发动机测量数据进行了融合分析。通过案例研究了当可测量参数存在误差时融合模型对流量、温度及压力等参数的预测效果;使用两种初值更新法,测试了不同工况下的物理运算时间;研究了约束违反程度对数据融合效果的影响;测试了关键传感器失效情况下,动态数据融合模型对缺失参数的预测效果,提出了进一步加快计算速度的方法。结果表明:通过数据融合,测量数据和未测量参数的不确定度下降到1%以内;通过算法优化,运算时间减小到5s,为发动机状态监控、传感器维护和传感器失效情况下发动机控制策略的制定提供了支持。

    • 航空发动机涡轮特性转换方法

      2020, 35(3):651-657. DOI: 10.13224/j.cnki.jasp.2020.03.021

      摘要 (591) HTML (0) PDF 2.51 M (856) 评论 (0) 收藏

      摘要:针对以落压比为横坐标,流量与效率为纵坐标,并按不同等转速线区分所表示的涡轮特性在涡轮处于临界状态时换算流量几乎保持不变,使得在航空发动机数学模型中应用该形式涡轮特性通过插值求解共同工作点时存在计算效率降低的不足,根据相似理论及等熵条件下涡轮膨胀功与落压比对应关系,推导了涡轮落压比、涡轮效率与涡轮当量功的关系,并结合抛物线插值方法给出了向以转速为横坐标,涡轮当量功与效率为纵坐标,并按不同等流量线区分所表示的涡轮特性转换方法。实例转换计算表明:采用所提出的转换计算方法,可有效解决以往涡轮特性插值流量基本不变的局限,并且两种格式的涡轮特性转换相对误差小于0.65%,满足工程要求。

    • >动力传输
    • 基于整机的中介轴承外圈剥落故障振动分析

      2020, 35(3):658-672. DOI: 10.13224/j.cnki.jasp.2020.03.022

      摘要 (484) HTML (0) PDF 4.95 M (589) 评论 (0) 收藏

      摘要:针对某典型双转子航空发动机中介轴承外圈疲劳剥落故障诊断问题,基于整机振动耦合动力学模型,导入中介轴承外圈早期剥落故障模型,通过数值积分方法获取故障激励下的整机振动响应。提取并分析了中介轴承外圈剥落故障特征,从信号分析中发现:①出现了4倍外圈故障特征频率及其两侧以外圈旋转频率为间隔的调制边频带;②随着不平衡量增加,特征频率分量基本不变,边频带变得更加突出;③随着轴承游隙的增加,外圈故障引起的冲击更为强烈,特征频率及其调制边频均显著增加。通过比较从轴承座到机匣各测点信号的故障特征,结果表明中介轴承早期疲劳剥落产生的冲击通过轴承座传递到机匣将产生很大程度的衰减,在机匣测点信号中轴承故障特征基本消失。

当期目录


文章目录

过刊浏览

年份

刊期

联系方式
  • 《航空动力学报》
  • 1986年创刊
  • 主办单位:
    中国航空学会  北京航空航天大学
  •  
  • 联系方式:
    北京市昌平区沙河高教园南三街9号 《航空动力学报》编辑部
  • 邮编:102206
  • 电话:010-61716749
  • 电子邮箱:JAP@buaa.edu.cn
  • 网址:
    http://www.jasp.com.cn
  •  
  • 刊号:
    1000-8055
  • CN11-2297/V
    公开发行
  • 国内定价:60元/期
您是第位访问者
中文版 ® 2024 版权所有
技术支持:北京勤云科技发展有限公司