• 2007年第22卷第8期文章目次
    全 选
    显示方式: |
    • >首届试验仿真与数据挖掘测试与分析会议(624)
    • 磁控等离子体对尾喷管壁传热特性的影响

      2007, 22(8):1209-1215.

      摘要 (1473) HTML (0) PDF 1002.89 K (400) 评论 (0) 收藏

      摘要:根据磁场作用下等离子体的湍流和传热能力将受到抑制的现象, 提出利用磁场控制低温等离子体隔离高温燃气与喷管壁的方法, 以减少高温燃气对壁面的传热, 从而达到降低壁面温度的目的.分别建立诱导磁场方程求解洛伦兹力和磁场作用下的k-ε湍流模型求解湍流粘度, 数值模拟了不同强度磁场作用下的磁控等离子体流动和传热特性.结果表明, 磁场能够有效地抑制湍流强度, 降低传热能力, 从而有效地降低壁面温度;并且磁场越强, 效果越明显.

    • 高位预旋进气转静盘腔换热实验

      2007, 22(8):1216-1221.

      摘要 (1550) HTML (0) PDF 685.62 K (460) 评论 (0) 收藏

      摘要:将某型发动机的涡轮盘腔结构简化成高位预旋进气的转静盘腔模型, 结合热色液晶先进测温显示技术, 通过实验的方法研究了转静盘间距、进气流量、旋转雷诺数以及旋流比等参数的变化对冷气降温效果和转盘换热效果的影响.实验结果表明:高位预旋进气对转盘外缘的冷却效果较好, 转盘上的温度分布呈现为同心圆形状;随着旋转雷诺数和转静盘间距的增大或者进气流量的减小, 冷气的总温降减小, 气流对转盘的换热效果变差.

    • 一种新颖的燃烧室出口温度场调试方法

      2007, 22(8):1222-1226.

      摘要 (1533) HTML (0) PDF 1.09 M (422) 评论 (0) 收藏

      摘要:以某小型回流燃烧室为基础验证平台, 试验研究改变燃油喷嘴插入火焰筒涡流器组件深度对燃烧室出口温度场的影响, 并对试验结果进行了比较和分析.结果表明:对同一燃烧室, 燃油喷嘴与火焰筒涡流器组件之间有一个最佳的轴向位置, 此时的燃烧室的出口温度分布系数最小;在相同的燃烧室工作状态, 采用不同喷雾锥角的燃油喷嘴, 这个最佳位置不同.该研究为燃烧室出口温度场的调试及燃烧室的研制提供了一条有效的技术途径.

    • 预旋进气位置对转静盘腔换热影响的数值研究

      2007, 22(8):1227-1232.

      摘要 (1593) HTML (0) PDF 939.21 K (489) 评论 (0) 收藏

      摘要:对具有不同预旋进气位置的转静盘腔内的换热进行了数值模拟.在进气流量恒定的条件下, 改变预旋孔在静盘上的径向位置, 研究了旋转雷诺数、预旋角、转静盘间距等参数对进出口气流的无量纲总温降和转盘换热效果的影响.研究表明:旋转雷诺数的增大和转静盘间距的减小都能够增强气流对转盘的换热效果;当预旋孔的径向位置增大时, 无量纲总温降和转盘表面的平均努塞尔数都是先增大后减小, 预旋孔存在一个最优的径向位置.

    • 数值分析二级涡流器环形燃烧室的燃烧性能

      2007, 22(8):1233-1240.

      摘要 (1503) HTML (0) PDF 2.03 M (419) 评论 (0) 收藏

      摘要:采用多维经验分析法在任意曲线坐标系下对包括二级突扩扩压器、双级轴向涡流器及火焰筒在内的单头部环形燃烧室的燃烧性能进行计算.并对RNGk-ε紊流模型加以改进, 使其更适用于数值模拟非交错贴体坐标系下的三维两相化学反应流.采用EBU-Arrhenius紊流燃烧模型、六通量热辐射模型以及颗粒群轨道模型对油雾燃烧进行模拟.同时, 数值分析了不同涡流器几何尺寸对燃烧室整体流场和燃烧性能的影响, 数值计算结果与实验数据相当一致, 表明计算方法合理, 计算程序可靠, 可用来估算环形燃烧室的燃烧性能.

    • 带隔热屏加力燃烧室热态流场计算

      2007, 22(8):1241-1246.

      摘要 (1496) HTML (0) PDF 903.78 K (466) 评论 (0) 收藏

      摘要:在任意曲线坐标系下对带有横向波纹隔热屏、外冷却通道和尾喷口的涡喷加力燃烧室的三维热态流场进行数值模拟.利用椭圆型微分方程和区域法生成三维贴体网格, k-ε模型预估紊流特性, EBU-Arre-henius紊流燃烧模型估算化学反应速率, 六通量模型预估辐射通量以考虑其对壁面和气流温度的影响.在非交错网格系统下应用混合差分格式离散控制方程, SIMPLE算法求解离散方程, 数值分析两种进口气流参数分布对加力室热态流场和壁温分布的影响.所得计算值与试验数据基本相符, 表明计算方法合理.

    • 周期电加热控制律对除冰表面温度的影响

      2007, 22(8):1247-1251.

      摘要 (1474) HTML (0) PDF 1.36 M (430) 评论 (0) 收藏

      摘要:研究周期性电热除冰系统加热控制律对表面温度的影响是进行系统设计和优化的前提.在采用焓法模型对二维电热除冰系统矩形单元的瞬态传热模型进行建模的基础上, 采用控制容积法全隐格式求解控制微分方程, 对周期性电热除冰表面的温度进行了计算.分析了影响表面温度的主要因素, 得到表面温度与加热控制律的关系, 结果表明:供电热流密度的大小和供电周期对系统除冰效果有重要影响, 需要合理匹配.

    • 竹节形扰流元对流动与换热特性影响的数值研究

      2007, 22(8):1252-1257.

      摘要 (1510) HTML (0) PDF 719.98 K (422) 评论 (0) 收藏

      摘要:采用数值计算方法, 研究了扰流柱形状对通道中流体流动和换热的影响, 模拟了带竹节形扰流柱通道的流场.结果表明:相对于圆柱形扰流柱, 装有圆柱竹节形和椭圆竹节形扰流柱阵列矩形通道的压力损失系数分别为62%~77%和25%~27%, 而恒热流壁面的平均奴塞尔数分别降低了9%22%和22%24%.随着扰流柱排列间距的减小, 压力损失增加, 而改变扰流柱排列横向间距引起的压力损失的变化要比改变流向间距显著.

    • Study on the breakup lengths of free round liquid jets

      2007, 22(8):1258-1263.

      摘要 (1295) HTML (0) PDF 636.13 K (445) 评论 (0) 收藏

      摘要:An experi ment was conducted to measure the breakup lengths of water jets with a high-speed video camera for Weber numbers from0 to about 1.1×103.The initial jet diameters are changed from0.3 mmto 1. 0 mm.The results indicate that at lowjet velocity the breakup lengths of the jets are increased linearly from 0 to a maxi mumvalue.At the Weber number about 20 the breakup length of the jet reaches its maxi mum value for different initial jet diameter.Acomputation based onthe dispersion equationis conducted to study the relationship between the growth rate of the jet surface wave and the maxi mumbreakup length.The computations showthat the maxi mumgrowth rate for the axisymmetric surface wave has a turning point at Weber number about 20, and that agrees well with the experi ments.

    • 减阻表面换热特性的实验

      2007, 22(8):1264-1267.

      摘要 (1415) HTML (0) PDF 536.83 K (459) 评论 (0) 收藏

      摘要:采用实验方法对带"V"型沟槽表面通道在不同雷诺数下的流动与换热特性进行了研究.研究结果表明:宽与深均为0.5 mm的沟槽板通道在实验范围内具有明显的减阻效果, 减阻效果随着雷诺数的增大而减小;相比平板通道流动, 低雷诺数范围下(10000<Re<30000)沟槽板的传热性能会减弱;高雷诺数范围下(30000<Re<80000)沟槽板的传热性能有所增强.且强化换热比在实验工况范围内随着雷诺数的增加而增加.

    • 气冷涡轮叶片三维换热问题计算

      2007, 22(8):1268-1272.

      摘要 (1499) HTML (0) PDF 618.20 K (456) 评论 (0) 收藏

      摘要:建立了三维气冷涡轮叶片换热计算的热-流耦合计算模型, 将辐射热量作为源项加入到方程中, 采用热-流耦合方法完成了某气冷涡轮叶片带辐射的三维换热计算.辐射热量的计算采用了P-1辐射模型, 比较了计算辐射和不计辐射时叶片表面的温度值及温度分布的差异.计算表明对于1600 K这样的高温燃气, 辐射热量不可忽略, 考虑辐射热量后, 叶片温度明显升高, 温度分布也产生了显著差异.

    • 微小“T”型蒸发管蒸发率实验

      2007, 22(8):1273-1278.

      摘要 (1445) HTML (0) PDF 714.30 K (439) 评论 (0) 收藏

      摘要:研究了微小尺寸"T"型蒸发管燃油蒸发率在不同工况下的变化规律.实验中蒸发管入口空气常压.实验研究了气油比、蒸发管入口空气温度、蒸发管壁温度和蒸发管入口空气流速对蒸发率的影响.试验结果表明, 气油比、蒸发管入口空气温度和蒸发管壁温是影响蒸发率的主要因素.蒸发管入口空气流速对蒸发率的影响较小.该蒸发管在管壁和进口空气同时加温, 并且气油比大于4.5的条件下, 蒸发率接近100%, 蒸发效果好.

    • >气动热力学与总体设计
    • 定几何二元倒置“X”型混压式超声速进气道数值仿真与实验验证

      2007, 22(8):1279-1284.

      摘要 (1444) HTML (0) PDF 666.64 K (424) 评论 (0) 收藏

      摘要:针对一种定几何二元倒置"X"型混压式超声速进气道开展了数值仿真研究, 结果表明:随着来流马赫数的增加, 进气道总压恢复系数下降, 而流量系数却先上升, 在设计点达到最大值后下降.当攻角变化时, 小攻角α<6°时迎背风两侧进气道总压恢复系数虽有下降但变化幅度不大, 对于流量系数, 在小攻角α<6°下背风侧进气道高于迎风侧进气道, 但两侧总的流量随攻角变化不大;在大攻角状态下(α=6°~9°), 背风侧进气道总压恢复系数和流量系数均下降剧烈, 而迎风侧进气道总压恢复系数虽然下降但流量系数却有上升.同时, 与实验结果对比表明, 两者规律趋势一致.

    • 扑翼微型飞行器飞行姿态模型研究

      2007, 22(8):1285-1289.

      摘要 (1604) HTML (0) PDF 1.68 M (483) 评论 (0) 收藏

      摘要:扑翼微型飞行器飞行质量主要取决于能否对其飞行姿态进行有效控制, 而建立准确的飞行姿态模型尤为重要.通过对鸟和昆虫的飞行机理尤其是其飞行过程中翅膀的运动规律进行研究, 并考虑机械设计方面的因素, 对扑翼微型飞行器的飞行姿态建立了较为完整的动力学模型和数学模型.由分析可知机身所受外力为空气动力、重力和机翼作用于机身的驱动力, 而采用扑动与扭转两个自由度飞行的机翼所产生的驱动力是由瞬时平移力和扭转循环力合成的瞬时空气动力.数值仿真证实了动力学模型的正确性.

    • 典型二元高超声速进气道设计方法研究

      2007, 22(8):1290-1296.

      摘要 (1612) HTML (0) PDF 614.28 K (431) 评论 (0) 收藏

      摘要:综合了一系列典型二元高超声速进气道的设计和性能估算方法, 给出了可行的设计原则.在满足流量、增压以及工作范围(起动性能和反压承受能力)的条件下, 给出了进气道进口、外压波系、内压缩通道、唇罩及隔离段的设计方法.采用此方法, 以H=22800 m、Ma0=6为设计点, 完成了一高超声速进气道的初步设计, 并估算得到了进气道性能参数、进气道的起动马赫数和反压承受能力, 对比CFD计算结果, 误差不大.通过该方法得到的进气道具有结构简单、流量系数大、压缩损失小的特点, 不通过优化即可得到性能较为良好的模型.

    • 高超声速飞行器-进气道一体化热流数值计算

      2007, 22(8):1297-1302.

      摘要 (1473) HTML (0) PDF 593.09 K (484) 评论 (0) 收藏

      摘要:采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.

    • 一种可用于声衬优化设计的管道声传播有限元模型

      2007, 22(8):1303-1308.

      摘要 (1520) HTML (0) PDF 473.97 K (523) 评论 (0) 收藏

      摘要:以航空发动机消声短舱声学设计为目标, 发展了一种管道声传播有限元模型, 它包含了壁面声阻抗、管道截面积和形状变化、以及管内非均匀流动等典型影响因素.应用该模型针对一种真实的高速转子声源进行了管道消声的数值模拟, 不仅得到了细致的管内声场分布情况, 更利用其快速的特点完成了声衬声阻抗的优化计算.结果表明该模型具有工程应用价值.

    • 基于CAIV的航空发动机性能与费用的综合权衡模型框架

      2007, 22(8):1309-1314.

      摘要 (1510) HTML (0) PDF 2.17 M (389) 评论 (0) 收藏

      摘要:运用"费用作为独立变量"方法, 对航空发动机立项论证中的性能与费用的综合权衡与决策方法进行了研究.首先建立了航空发动机总拥有费用的参数估算模型, 研究了发动机主要性能参数与其主要设计参数之间的关系并建立了性能模型.然后研究了如何确定设计参数的约束条件即"权衡空间".最后将以上研究成果进行了综合, 建立了发动机性能与总拥有费用的综合权衡模型框架, 可用于发动机的设计参数与费用的决策.

    • 微型扑翼飞行器的升力风洞试验

      2007, 22(8):1315-1319.

      摘要 (1450) HTML (0) PDF 931.58 K (477) 评论 (0) 收藏

      摘要:针对微型扑翼飞行器在专门的低雷诺数试验风洞中进行了扑翼的升力风洞试验, 测试出了微型扑翼飞行器在正常飞行以及向上爬升状态时的升力特性曲线, 并且探讨了扑翼频率、飞行速度以及扑翼的迎角对扑翼飞行升力的影响, 在此基础上, 还分析了扑翼升力的组成部分以及相互关系.实验结果为微型扑翼飞行器总体、气动的进一步设计提供了依据.

    • 计算涡扇发动机风车起动特性的辨识模型

      2007, 22(8):1320-1324.

      摘要 (1665) HTML (0) PDF 1.51 M (529) 评论 (0) 收藏

      摘要:针对高空风车条件下涡扇发动机起动性能估算问题, 提出了运用约简遗传规划算法进行模型结构自动选择的思路, 对基于约简遗传规划的风车起动模型进行了研究.以某型涡扇发动机为例, 运用该算法对给定条件下该型发动机的风车起动过程进行了建模仿真, 得到了起动过程的估算结果.将估算结果与试验数据对比, 验证了方法的正确性、可行性.

    • 复杂形式球型收敛调节片喷管内流场计算及分析

      2007, 22(8):1325-1329.

      摘要 (1596) HTML (0) PDF 669.40 K (460) 评论 (0) 收藏

      摘要:为了分析复杂形式球型收敛调节片喷管(SCFN)三维内流场特性, 对六角形和领结形的SCFN进行几何建模, 并利用N-S方程和S-A模型, 对两种喷管在不同落压比(NPR)下进行了内流场数值模拟.研究表明, 计算结果与实验结果吻合较好;在此基础上, 对两种喷管的流场特性及性能进行了比较分析.

    • 平板上有叶片绕流的源项模型

      2007, 22(8):1330-1334.

      摘要 (1511) HTML (0) PDF 724.71 K (421) 评论 (0) 收藏

      摘要:将源项模型应用于叶片涡流发生器的数值模拟中, 以叶片在当地流场中产生的升力作为源项一个分量原型, 并引入了法向力作为源项另一个分量原型.该源项由叶片的形状、大小、位置以及来流条件等参数决定.在数值模拟中, 叶片位于平板上, 其高度约为当地边界层厚度的1/5.在相同计算条件下, 对"升力-法向力"源项模型计算结果和全网格计算结果及实验数据进行比较.源项模型可以得到与全网格数值模拟基本一致的诱导速度场和涡量场, 而网格量减少了65%.

    • >大飞机动力
    • 成败型与指数型单元混合系统的可靠性熵法评定

      2007, 22(8):1335-1339.

      摘要 (1347) HTML (0) PDF 288.42 K (444) 评论 (0) 收藏

      摘要:利用信息量等效的基本原则, 研究了成败型单元与指数型单元混合组成系统的可靠性评定问题.给出了成败型单元与指数型单元混合组成成败型系统或指数型系统时, 单元试验信息熵法折合的基本模型, 导出了成败型信息与指数型信息相互折合的基本公式, 并通过算例及其评定结果的分析比较说明了方法的具体应用及其有效性.

    • 基于旋转坐标系转轴振动信号的滚动轴承故障诊断方法

      2007, 22(8):1340-1345.

      摘要 (1485) HTML (0) PDF 1.41 M (484) 评论 (0) 收藏

      摘要:实验验证了基于旋转坐标系转轴振动信号的滚动轴承局部故障诊断方法.建立起滚动轴承故障实验台, 采用压电晶体加速度传感器测取转轴的振动.通过滚动轴承内环、外环、滚动体和保持架四种典型故障的实验分析结果, 证明了该方法的可行性, 并发现转轴振动信号的传递路径相对简单, 对轴承内环、滚动体和保持架早期微小故障的发现与诊断可能是有利的.

    • 涡轮冷却叶片参数化建模及多学科设计优化

      2007, 22(8):1346-1351.

      摘要 (1533) HTML (0) PDF 925.17 K (520) 评论 (0) 收藏

      摘要:建立了一个涉及结构、气动、传热、振动、强度和寿命等学科的涡轮冷却叶片多学科设计优化系统, 进行了单孔薄壁冷却叶片的多学科设计优化.提出了单孔薄壁冷却叶片的参数化造型方法, 叶片叶型采用5次多项式构造, 气动与传热为三维耦合分析;叶片体积平均温度与最高温度为优化目标, 强度、振动和寿命等学科相关参数为约束, 模拟退火与序列二次规划组合算法进行叶片参数空间寻优, 在保持冷却气体流量不变的条件下, 优化提高了冷却效果, 降低了叶片材料的性能要求.

    • 发动机零部件稳态应力计算的状态函数法

      2007, 22(8):1352-1355.

      摘要 (1350) HTML (0) PDF 297.69 K (434) 评论 (0) 收藏

      摘要:借鉴发动机气动参数相似准则的基本思想, 提出了根据典型工作点的应力结果直接推算出其它状态点的状态函数法.利用该方法可以快速计算发动机任意状态点的稳态应力、位移、温度, 可以便捷地为构件寿命预估提供应力谱, 使发动机零部件寿命计算过程大大简化.该方法在发动机零部件寿命预估和管理等方面有较高的实用价值.

    • >叶轮机械
    • 发动机零部件稳态应力计算的状态函数法

      2007, 22(8):1356-1364.

      摘要 (1261) HTML (0) PDF 1.30 M (423) 评论 (0) 收藏

      摘要:An opti mization process is used to redesign blades of a high-pressure compressor.An artificial neural network (ANN) method is coupled to Navier-Stokes solvers and is applied to three different redesigns.A newrotor blade of a transonic compressor is designed by modifying thick, stacking line andinlet angle using a 3Dapproach, with a significant efficiencyi mprovement at the design point.The off-design behavior of this new compressor is also checked afterwards, which shows that the whole performance of the inlet stage is improved over a wide range of mass flow.The losses are reduced, proving the good performance of the opti mum.The whole results indicate that the opti mization method can find i mproved design and can be integrated in a design procedure.

    • 吸附式跨声速压气机叶栅流场数值模拟

      2007, 22(8):1365-1370.

      摘要 (1517) HTML (0) PDF 1.84 M (474) 评论 (0) 收藏

      摘要:使用MISES程序数值模拟了跨声速吸附式压气机叶栅流场, 重点研究了吸气量和吸气位置对跨声速压气机叶栅气动性能的影响.结果表明, 叶栅来流马赫数和方向一定时, 吸气位置和吸气量是相互关联的关键参数, 不同的吸气位置对应着不同的最佳吸气量, 且随着吸气位置向后缘远离激波, 最佳吸气量呈逐渐增大之势.从吸气对叶片吸力面边界层的影响效果分析, 理想的吸气位置应该是在靠近激波后附面层发展到一局部极大值即将进入过渡段的位置附近.

    • NASA67级非定常流场的频域分析

      2007, 22(8):1371-1377.

      摘要 (1527) HTML (0) PDF 1.25 M (406) 评论 (0) 收藏

      摘要:通过对NASA67级典型截面非定常流场的快速傅立叶分解, 在频域内对转静干扰效应、交界面典型脉动量的周向波形以及周期边条的相位延迟角进行了分析.同时可以看到一些在时域中很难直观看到的现象.最后通过对比采用前5阶谐波合成出来的确定应力和非定常计算合成出来的确定应力, 分析了非线性谐波法在确定应力建模方面的可行性及优势.

    • >自动控制
    • 考虑切换特性的航空发动机多回路控制系统设计

      2007, 22(8):1378-1383.

      摘要 (1821) HTML (0) PDF 1.61 M (454) 评论 (0) 收藏

      摘要:针对航空发动机多回路控制系统切换, 分析了多回路控制系统切换的性能指标要求, 基于极点配置与H鲁棒控制综合的方法, 构造了线性矩阵不等式(LMIs), 以此为充分条件给出了多回路切换时的H∞控制器设计方法, 并与单回路系统设计切换进行比较仿真, 进一步表明了航空发动机切换控制的必要性与有效性.

    • 涡扇发动机高空起动模型研究

      2007, 22(8):1384-1390.

      摘要 (1655) HTML (0) PDF 1.93 M (439) 评论 (0) 收藏

      摘要:基于慢车以上的部件级模型, 通过外推部件特性, 建立了双转子涡扇发动机部件级高空起动模型, 可执行高空风车起动或起动机起动仿真.模型考虑了起动过程中油气比对燃烧效率的影响, 燃烧室、高低压涡轮的热惯性及起动过程中总压恢复系数的变化, 采用容积动力学方法避免了迭代运算, 仿真结果表明该模型具有实时性, 仿真转速误差小于5%, 能够满足发动机高空起动性能仿真的需求.

    • 基于RBF网络辨识的涡扇发动机双变量神经网络PID解耦控制

      2007, 22(8):1391-1395.

      摘要 (1496) HTML (0) PDF 602.41 K (395) 评论 (0) 收藏

      摘要:根据神经网络与PID算法相结合的思想, 针对涡扇发动机双变量控制中变量之间的耦合问题, 提出基于径向基函数神经网络(RBF)辨识的发动机双变量神经网络PID解耦控制, 并给出控制系统的控制结构及原理.仿真结果表明, 该方法控制精度高、跟踪性能强、鲁棒性良好, 能够有效地减小各回路之间的耦合影响, 并保证控制系统具有良好的稳态和动态性能, 适合航空发动机控制.

    • 航空发动机数控系统中LVDT传感器信号处理及在线故障检测

      2007, 22(8):1396-1400.

      摘要 (1400) HTML (0) PDF 489.91 K (429) 评论 (0) 收藏

      摘要:通过对线位移差动变压器(LVDT)信号处理电路的深刻研究, 给出了芯片发热、温漂、信号精度等问题的一般解决方法, 提出了一种通过软、硬件结合, 在不增加外围器件的情况下对LVDT传感器进行在线故障检测的定值判断法.发动机台架试验结果表明, 所提出的方法是行之有效的.

当期目录


文章目录

过刊浏览

年份

刊期

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