划船器项目可行性研究报告征地用(专业版)

划船器项目

可行性研究报告

（征地用/专业版）

普慧投资研究中心

划船器项目

可行性研究报告

（征地用/专业版）

项目负责人：齐宪臣注册咨询工程师

参加人员：郑西芳注册咨询工程师

胡冰月注册咨询工程师

王子奇高级经济师

杜翔宇高级工程师

项目审核人：张子宏注册咨询工程师

普慧投资研究中心

目录

划船器项目可行性研究报告常见问题解答 ............ 错误!未定义书签。

1、划船器项目应该在经信委还是发改委立项？ (1)

2、编制划船器项目可行性研究报告企业需提供的资料清单 (1)

一、总论 (2)

（一）项目背景 (2)

1、项目名称 (2)

2、建设单位概况 (2)

3、可行性研究报告编制依据 (2)

4、项目提出的理由与过程 (3)

（二）项目概况 (3)

1、拟建项目 (3)

2、建设规模与目标 (3)

3、主要建设条件 (3)

4、项目投入总资金及效益情况 (4)

5、主要技术经济指标 (4)

（三）主要问题说明 (6)

1、项目资金来源问题 (6)

2、项目技术设备问题 (6)

3、项目供电供水保障问题 (6)

二、市场预测 (7)

（一）划船器市场分析 (7)

1、国际市场 (7)

2、国内市场 (7)

（二）主要竞争企业分析（略） (8)

（三）目标市场分析 (9)

1、目标市场调查 (9)

2、价格现状与预测 (10)

（四）营销策略 (10)

1、销售队伍建设 (10)

2、销售网络建设 (10)

3、销售策略 (10)

三、建设规模与产品方案 (12)

（一）建设规模 (12)

（二）产品方案 (12)

四、场址选择 (13)

（一）场址所在位置现状 (13)

1、地点与地理位置 (13)

2、场址土地权属类别及占地面积 (13)

3、土地利用现状 (14)

（二）场址建设条件 (14)

1、地理环境位置 (14)

2、地形、地貌 (14)

3、气候、水文 (14)

4、交通运输条件 (14)

5、公用设施社会依托条件 (14)

6、环境保护条件 (15)

7、法律支持条件 (15)

8、征地、拆迁、移民安置条件 (15)

9、施工条件 (15)

五、技术方案、设备方案和工程方案 (16)

（一）技术方案 (16)

1、生产方法 (16)

2、工艺流程 (17)

（二）主要设备方案 (18)

1、设备选配原则 (18)

2、设备选型表 (19)

（三）工程方案 (20)

1、土建工程设计方案 (20)

2、主要建、构筑物的建筑特征、结构及面积方案 (21)

3、建筑及安装工程量及造价 (22)

六、主要原材料、燃料供应 (23)

（一）主要原料材料供应 (23)

（二）燃料及动力供应 (23)

（三）主要原材料、燃料及动力价格 (23)

（四）主要原材料、燃料年需要量表 (24)

七、总图运输与公用辅助工程 (25)

（一）总图布置 (25)

1、平面布置 (25)

2、竖向布置及道路 (25)

3、总平面图 (25)

4、总平面布置主要指标表 (28)

（二）场内外运输 (28)

1、场外运输量及运输方式 (28)

2、场内运输量及运输方式 (28)

3、场外运输设施及设备 (29)

（三）公共辅助工程 (29)

1、供水工程 (29)

2、供电工程 (30)

3、通信系统设计方案 (35)

4、通风采暖工程 (36)

5、防雷设计 (37)

6、防尘设计 (37)

7、维修及仓储设施 (38)

八、节能措施 (39)

（一）节能措施 (39)

1、节能规范 (39)

2、设计原则 (39)

3、节能方案 (39)

（二）能耗指标分析 (42)

1、用能标准与能耗计算方法 (42)

2、能耗状况和能耗指标分析 (43)

九、节水措施 (44)

（一）节水措施 (44)

（二）水耗指标分析 (44)

十、环境影响评价 (45)

（一）场址环境条件 (45)

（二）项目建设和生产对环境的影响 (45)

1、项目建设对环境的影响 (45)

2、项目生产对环境的影响 (46)

（三）环境保护措施方案 (47)

1、设计依据 (47)

2、环保措施 (47)

（四）环境保护投资 (49)

（五）环境影响评价 (49)

十一、劳动安全卫生与消防 (50)

（一）劳动安全与职业卫生 (50)

1、设计依据 (50)

2、设计执行的主要标准 (50)

3、设计内容及原则 (50)

4、职业安全 (50)

5、职业卫生 (51)

6、辅助卫生用室 (51)

7、职业安全卫生机构 (51)

（二）消防 (51)

1、设计依据 (51)

2、总平面布置 (52)

3、建筑部分 (52)

4、电气部分 (52)

5、给排水部分 (52)

十二、组织机构与人力资源配置 (53)

（一）组织机构 (53)

1、项目法人组建方案 (53)

2、管理机构组织方案 (53)

（二）人力资源配置 (53)

1、生产作业班次 (53)

2、项目劳动定员 (53)

3、职工工资福利 (53)

4、员工来源及招聘方案 (54)

5、员工培训 (54)

十三、项目实施进度 (55)

（一）建设工期 (55)

（二）项目实施进度安排 (55)

（三）项目实施进度表 (55)

十四、招标方案 (56)

（一）编制招标计划的依据 (56)

（二）招标内容 (56)

十五、投资估算 (58)

（一）投资估算依据 (58)

（二）建设投资估算 (58)

1、建筑工程费 (58)

2、设备及工器具购置费 (58)

3、安装及装修工程费 (58)

4、土地购置及整理费 (59)

5、工程建设其他费用 (59)

6、基本预备费 (59)

7、涨价预备费 (59)

8、建设期利息 (59)

（三）流动资金估算 (59)

（四）项目投入总资金 (59)

（六）投资使用计划 (59)

十六、融资方案 (60)

（一）资本金筹措 (60)

（二）债务资金筹措 (60)

（三）融资方案分析 (60)

十七、财务评价 (61)

（一）计算依据及相关说明 (61)

1、项目测算参考依据 (61)

2、项目测算基本设定 (61)

（二）销售收入、销售税金及附加和增值税估算 (62)

1、销售收入 (62)

2、销售税金及附加费用 (62)

（三）总成本费用估算 (62)

1、直接成本 (62)

2、工资及福利费用 (62)

3、折旧及摊销 (62)

4、修理费 (62)

5、财务费用 (63)

6、其它费用 (63)

7、总成本费用 (63)

（四）财务评价报表 (63)

1、项目损益及利润分配表 (63)

2、项目财务现金流量表 (63)

（五）财务评价指标 (63)

1、投资利润率，投资利税率 (63)

2、财务内部收益率、财务净现值、投资回收期 (64)

（七）不确定性分析 (64)

1、敏感性分析 (64)

2、盈亏平衡分析 (64)

（八）财务评价结论 (65)

十八、项目经济效益与社会效益 (66)

（一）经济效益 (66)

（二）社会效益 (66)

十九、风险分析 (67)

（一）项目风险因素识别 (67)

1、法律及政策风险 (67)

2、市场风险 (67)

3、建设风险 (67)

4、环保风险 (67)

（二）项目风险防控措施 (67)

1、法律及政策风险防控措施 (67)

2、市场风险防控措施 (67)

3、建设风险防控措施 (68)

4、环保风险防控措施 (68)

二十、结论与建议 (69)

（一）结论 (69)

（二）建议 (69)

二十一、附件 (70)

（一）附表 (70)

（二）附图 (78)

附表：

1、附表1 项目建筑工程费估算表

2、附表2 项目设备及工器具购置费估算表

3、附表3 工程建设其他费用估算表

4、附表4 流动资金估算表（万元）

5、附表5 项目投入总资金估算表（万元）

6、附表6 项目投入总资金使用计划表（万元）

7、附表7 项目销售税金及附加费用（万元）

8、附表8 项目直接成本表（万元）

9、附表9 项目摊销估算表（万元）

10、附表10 项目折旧估算表（万元）

11、附表11 项目总成本费用估算表（万元）

12、附表12 项目损益及利润分配表（万元）

13、附表13 项目财务现金流量表（万元）

附图：

1、建设项目地理位置图

2、项目厂区平面布置图

附件：

1、企业法人营业执照

2、项目备案请示

划船器项目可行性研究报告常见问题解答

1

2

3

4

5

6

7

8

9

飞行器系统仿真

《飞行器系统仿真与CAD 》学习报告 第一部分仿真(40) 题目1：给定导弹相对于目标的运动学方程组为 q k q V q V q r q V q V r m m ,sin )sin(),cos(cos r(0) = 5km, q(0) = 60deg, (0) = 30deg,V = , V m = , 1Ma = 340m/s, k = 2 （1） 建立系统的方框图模型； （2） 用MATLAB 语言编写S —函数 （3） 用窗口菜单对(1), (2)进行仿真，动态显示结果; （4）用命令行对(1), (2)进行仿真，以图形显示结果 答： （1） （2）用MATLAB 语言编写S 函数 function [sys,x0,str,ts]=CAD1_sfun(t,x,u,flag) switch flag case 0 [sys,x0,str,ts]=mdlInitializeSizes; case 1 sys = mdlDerivatives(t,x,u); case 3 sys = mdlOutputs(t,x,u); case {2,4,9} sys = []; otherwise

error('unhandled flag=',num2str(flag)) end function [sys,x0,str,ts]=mdlInitializeSizes sizes=simsizes; =3; =0; =3; =0; =1; =1; sys=simsizes(sizes); str=[]; x0=[5000,pi/3,pi/6]; ts=[0 0]; function sys=mdlDerivatives(t,x,u) vm=*340; v=*340; k=2; dx(1)=vm*cos(x(2))-v*cos(x(2)-x(3)); dx(2)=(v*sin(x(2)-x(3))-vm*sin(x(2)))/x(1); dx(3)=k*dx(2); sys=dx; function sys=mdlOutputs(t,x,u) sys=x; 调用S函数的模型框图 (3)框图仿真结果：

飞行器仿真原理

在无风、无侧滑的情况下，飞行器航迹坐标系下的运动学方程[2]为： cos()sin sin()cos cos cos cos sin()sin sin s v p t p s s t p s s t d m p Q mg d d m p Y mg d d mv p Y d θψαθαγγθθαγγ?=+Φ--????=+Φ+-????-=+Φ+?? （1） 其中m 为飞行器质量；v 为速度；p 为发动机动力；α为迎角；p Φ为发动机安装角；Q 为空间阻力；θ为俯仰角；s γ为滚转角；Y 为升力；s ψ为偏航角。 图1显示了机体坐标系下的飞行器受力情况；图2表示了地面坐标系和航迹坐标系的关系；其中，d d d Ox y z 表示地面坐标系，地面坐标系固定于地面，原点选在 地面的某一点，d y 铅直向上，d x 和d z 在水平面内。h h h Ox y z 表示航迹坐标系，航 迹坐标系原点在飞行器质心，h x 沿飞行器速度向量v ，即飞行器飞行方向，h y 在包含v 的铅垂直平面内，h z 垂直于铅垂平面。 图1 飞行器受力分析 图2 地面坐标系和轨迹坐标系 为了更清晰、简练地描述这些运动学的量，我们令 cos()/sin()cos cos /sin()cos sin /x p y p s s z p s s n p Q mg n p Y mg n p Y mg ααγγαγγ???=+Φ-??????=+Φ+??????=+Φ+???? （2） 称x n 、y n 、z n 为过载，把（2）式代入（1）式得到 []()()()()sin ()()()cos ()()cos ()()s v t z t t y t t x t d n t t g d d v t n t t g d d v t t n t g d θψθθθ?=-??????=-??????-=?? （3） 从式（3）可以明确看出：x n 、y n 、z n 反映了飞行器因主动运动而产生的加速度，而sin θ和cos θ则是由于飞行器的重力产生的加速度。

四旋翼飞行器建模与仿真Matlab

四轴飞行器的建模与仿真 摘要 四旋翼飞行器是一种能够垂直起降的多旋翼飞行器,它非常适合近地侦察、监视的任务,具有广泛的军事和民事应用前景。本文根据对四旋翼飞行器的机架结构和动力学特性做详尽的分析和研究,在此基础上建立四旋翼飞行器的动力学模型。四旋翼飞行器有各种的运行状态,比如:爬升、下降、悬停、滚转运动、俯仰运动、偏航运动等。本文采用动力学模型来描述四旋翼飞行器的飞行姿态。在上述研究和分析的基础上,进行飞行器的建模。动力学建模是通过对飞行器的飞行原理和各种运动状态下的受力关系以及参考牛顿-欧拉模型建立的仿真模型，模型建立后在Matlab/simulink软件中进行仿真。 关键字:四旋翼飞行器,动力学模型，Matlab/simulink Modeling and Simulating for a quad-rotor aircraft ABSTRACT The quad-rotor is a VTOL multi-rotor aircraft. It is very fit for the kind of reconnaissance mission and monitoring task of near-Earth, so it can be used in a wide range of military and civilian applications. In the dissertation, the detailed analysis and research on the rack structure and dynamic characteristics of the laboratory four-rotor aircraft is showed in the dissertation. The dynamic model of the four-rotor aircraft areestablished. It also studies on the force in the four-rotor aircraft flight principles and course of the campaign to make the research and analysis. The four-rotor aircraft has many operating status, such as climbing, downing, hovering and rolling movement, pitching movement and yawing movement. The dynamic model is used to describe the four-rotor aircraft in flight in the dissertation. On the basis of the above analysis, modeling of the aircraft can be made. Dynamics modeling is to build models under the principles of flight of the aircraft and a variety of state of motion, and Newton - Euler model with reference

飞行器六自由度仿真

1引言 现在的战争已不是过去大刀长矛的时代，他早已成为国家综合实力的体现，这很大程度取决于军事高科技。这其中导弹作为精确打击的利器关乎国家的战略安全。而研究其包括导弹在内的飞行器精确制导与控制便显得十分的重要。 飞行器最优制导规律研究是进行武器系统总体方案论证和提高制导性能及精度的关键技术之一。而要进行制导规律最优性研究一方面需要研究合适的制导规律,另一方面需要进行接近实际情况的全面的大量的仿真研究。 仿真验证包括建模与仿真两个方面。在大型工程的方案论证阶段甚至包括实际研制的各个阶段，都要进行仿真检验以论证可行性、合理性和最优性。仿真技术在工程应用特别在高端武器系统总体设计和方案论证中具有极为重要的作用。对制导问题的研究在国外倍受重视。在公开发表的文献中,专门讨论制导规律方面的研究论文很多，可见制导规律的研究是非常重要的。但是仅有理论研究是远远不够的，因为设计的所谓最优制导规律大都是把实际系统进行了大量简化情况下推导出来的，因而与实际情况差别较大。也就是说理论上是最优的制导规律或参数在实际系统中不一定是最优的。因此,必须建立接近实际状态的数学模型和仿真软件。通过仿真计算确定出制导系统的最优参数及制导规律的控制效果，才能最终确定制导规律的最优性。 目前国内外这类问题研究主要存在下列三个问题:其一是模型被简化。从众多公开发表的文献资料看,大都是把控制系统简化为一阶、二阶、或三阶等根模型来推导制导规律,并据此模型进行仿真计算。其二是把飞行器的六自由度空间运动状态简化为平面运动状态进行仿真研究，以这种把飞行器运动限制在平面范围内进行仿真计算是有局限性的。其三是在全弹道数字仿真中仅选取几个特征点参数来代表全弹道的气动力参数，这些参数要表征全弹道动态过程是比较片面的，因而仿真结论的可信度是比较差的。 若把飞行器看成一个刚体，则它在空间的运动，可以看做是质心的移动和绕质心的转动的合成运动。质心的移动取决于作用在飞行器上的力，绕质心的转动则取决于作用在飞行器上相对于质心的力矩。在飞行中，作用在飞行器上的力主要有：总空气动力、发动机的推力和重力等。作用在飞行器上的力矩有：空气动力引起的空气动力矩，由发动机推力(若推力作用线不通过飞行器质心时)引起的推力矩等。作用在导弹上的空气动

四旋翼飞行器仿真-实验报告

动态系统建模仿真实验报告（2） 四旋翼飞行器仿真 姓名： 学号： 指导教师： 院系： 2014.12.28

1实验容 基于Simulink建立四旋翼飞行器的悬停控制回路，实现飞行器的悬停控制； 建立GUI界面，能够输入参数并绘制运动轨迹； 基于VR Toolbox建立3D动画场景，能够模拟飞行器的运动轨迹。 2实验目的 通过在 Matlab 环境中对四旋翼飞行器进行系统建模，使掌握以下容： 四旋翼飞行器的建模和控制方法 在Matlab下快速建立虚拟可视化环境的方法。 3实验器材 硬件：PC机。 工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。 4实验原理 4.1四旋翼飞行器 四旋翼飞行器通过四个螺旋桨产生的升力实现飞行，原理与直升机类似。四个旋翼位于一个几何对称的十字支架前，后，左，右四端，如图 1 所示。旋翼由电机控制；整个飞行器依靠改变每个电机的转速来实现飞行姿态控制。 图1四旋翼飞行器旋转方向示意图

在图 1 中， 前端旋翼 1 和后端旋翼 3 逆时针旋转， 而左端旋翼 2 和右端的旋翼 4 顺时针旋转， 以平衡旋翼旋转所产生的反扭转矩。 由此可知， 悬停时， 四只旋翼的转速应该相等，以相互抵消反扭力矩；同时等量地增大或减小四只旋翼的转速，会引起上升或下降运动；增大某一只旋翼的转速，同时等量地减小同组另一只旋翼的转速，则产生俯仰、横滚运动；增大某一组旋翼的转速，同时等量减小另一组旋翼的转速，将产生偏航运动。 4.2建模分析 四旋翼飞行器受力分析,如图 2 所示 图2四旋翼飞行器受力分析示意图 旋翼机体所受外力和力矩为： 重力mg , 机体受到重力沿w z -方向； 四个旋翼旋转所产生的升力i F (i= 1 , 2 , 3 , 4)，旋翼升力沿b z 方向； 旋翼旋转会产生扭转力矩i M (i= 1 , 2 , 3 , 4)。i M 垂直于叶片的旋翼平面，与旋转矢量相反。 力模型为：2i F i F k ω= ，旋翼通过螺旋桨产生升力。F k 是电机转动力系数， 可取826.1110/N rpm -?，i ω为电机转速。旋翼旋转产生旋转力矩Mi(i=1,2,3,4)，

四旋翼飞行器仿真 实验报告

动态系统建模仿真 实验报告（2）四旋翼飞行器仿真 姓名： 学号： 指导教师： 院系： 2014.12.28

1实验内容 基于Simulink建立四旋翼飞行器的悬停控制回路，实现飞行器的悬停控制； 建立GUI界面，能够输入参数并绘制运动轨迹； 基于VR Toolbox建立3D动画场景，能够模拟飞行器的运动轨迹。 2实验目的 通过在Matlab 环境中对四旋翼飞行器进行系统建模，使掌握以下内容：四旋翼飞行器的建模和控制方法 在Matlab下快速建立虚拟可视化环境的方法。 3实验器材 硬件：PC机。 工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。 4实验原理 4.1四旋翼飞行器 四旋翼飞行器通过四个螺旋桨产生的升力实现飞行，原理与直升机类似。四个旋翼位于一个几何对称的十字支架前，后，左，右四端，如图 1 所示。旋翼由电机控制；整个飞行器依靠改变每个电机的转速来实现飞行姿态控制。 图1四旋翼飞行器旋转方向示意图

在图 1 中， 前端旋翼 1 和后端旋翼 3 逆时针旋转， 而左端旋翼 2 和右端的旋翼 4 顺时针旋转， 以平衡旋翼旋转所产生的反扭转矩。 由此可知， 悬停时， 四只旋翼的转速应该相等，以相互抵消反扭力矩；同时等量地增大或减小四只旋翼的转速，会引起上升或下降运动；增大某一只旋翼的转速，同时等量地减小同组另一只旋翼的转速，则产生俯仰、横滚运动；增大某一组旋翼的转速，同时等量减小另一组旋翼的转速，将产生偏航运动。 4.2建模分析 四旋翼飞行器受力分析,如图 2 所示 图2四旋翼飞行器受力分析示意图 旋翼机体所受外力和力矩为： 重力mg , 机体受到重力沿w z -方向； 四个旋翼旋转所产生的升力i F (i= 1 , 2 , 3 , 4)，旋翼升力沿b z 方向； 旋翼旋转会产生扭转力矩i M (i= 1 , 2 , 3 , 4)。i M 垂直于叶片的旋翼平面，与旋转矢量相反。 力模型为：2i F i F k ω= ，旋翼通过螺旋桨产生升力。F k 是电机转动力系数， 可取826.1110/N rpm -?，i ω为电机转速。旋翼旋转产生旋转力矩Mi(i=1,2,3,4)，

飞行器系统仿真

《飞行器系统仿真与C A D》学习报告 第一部分仿真(40) 题目1：给定导弹相对于目标的运动学方程组为 r(0) = 5km, q(0) = 60deg, (0) = 30deg,V = , V m= , 1Ma = 340m/s, k = 2 （1）建立系统的方框图模型； （2）用MATLAB语言编写S—函数 （3）用窗口菜单对(1), (2)进行仿真，动态显示结果; （4）用命令行对(1), (2)进行仿真，以图形显示结果 答： （1） （2）用MATLAB语言编写S函数 function [sys,x0,str,ts]=CAD1_sfun(t,x,u,flag) switch flag case 0 [sys,x0,str,ts]=mdlInitializeSizes; case 1 sys = mdlDerivatives(t,x,u); case 3 sys = mdlOutputs(t,x,u); case {2,4,9} sys = []; otherwise error('unhandled flag=',num2str(flag)) end function [sys,x0,str,ts]=mdlInitializeSizes sizes=simsizes; =3; =0; =3;

=0; =1; =1; sys=simsizes(sizes); str=[]; x0=[5000,pi/3,pi/6]; ts=[0 0]; function sys=mdlDerivatives(t,x,u) vm=*340; v=*340; k=2; dx(1)=vm*cos(x(2))-v*cos(x(2)-x(3)); dx(2)=(v*sin(x(2)-x(3))-vm*sin(x(2)))/x(1); dx(3)=k*dx(2); sys=dx; function sys=mdlOutputs(t,x,u) sys=x; 调用S 函数的模型框图 (3)框图仿真结果： S 函数仿真结果： （4）命令输入 clear;clc [t x ] = sim('CAD1'); hSimulink = figure(); subplot(3, 1, 1); plot(t,x(:,1)); grid; ylabel('r'); subplot(3, 1, 2); plot(t,x(:,2)); grid; ylabel('q'); subplot(3, 1, 3); plot(t,x(:,3)); grid; ylabel('sigma'); [t x ] = sim('CAD1_S'); hSFun = figure(); subplot(3, 1, 1); plot(t,x(:,1)); grid; ylabel('r'); subplot(3, 1, 2); plot(t,x(:,2)); grid; ylabel('q'); subplot(3, 1, 3); plot(t,x(:,3)); grid; ylabel('sigma'); 模型仿真结果： S 函数仿真结果： 题目2：给出动态方程 0)0(,1)0(;1)1(2===+-+? ???x x tx x x x ; (1) 用MATLAB 语言编写S —函数; (2) 用命令行gear/adams 法对(1)进行仿真，显示曲线x(t=0：100); (3) 建立方框图,用RK45仿真50秒,显示曲线

飞行器建模与仿真

飞行器建模与仿真作业 班级：飞设12 姓名：潘周周 学号：2110702015

作业一 一、题目 已知微分方程组 x x 2 1=? x x 32=? x x 4 3=? x x 5 4=? x x 6 5 =? ??? ?????------=? w x a w x a w x a w x a w x a x x a w x 5 6 54054304320320 211660 6 用四阶龙哥库塔法求解该方程组。 二、程序 1.编写rk4f.m 函数 function r=rk4f(x1,x2,x3,x4,x5,x6,xs,ts,tf,N) %四阶龙格库塔法计算微分方程组程序 %x1,x2,x3,x4,x5,x6分别是初值 %xs 是方程组中的x %ts ，tf 代表代表左端点和右端点 %N 是迭代步数 w0=50; a1=3.86; a2=7.46; a3=9.13; a4=7.46; a5=3.86; a6=1; aa=zeros(1,N+1); ab=zeros(1,N+1); ac=zeros(1,N+1); ad=zeros(1,N+1); ae=zeros(1,N+1); af=zeros(1,N+1);

aa(1)=x1; ab(1)=x2; ac(1)=x3; xd(1)=x4; xe(1)=x5; af(1)=x6; h=(tf-ts)/N; t=ts:h:tf; for i=1:1:N; %ki1(i=1,2,...6) k11=x2; k21=x3; k31=x4; k41=x5; k51=x6; k61=w0^6/a6*(xs-x1-a1*x2/w0-a2*x3/w0^2-a3*x4/w0^3-a4*x5/w0^4-a5*x6/w0 ^5); %ki2(i=1,2,...6) k12=x2+h/2*k21; k22=x3+h/2*k31; k32=x4+h/2*k41; k42=x5+h/2*k51; k52=x6+h/2*k61; k62=w0^6/a6*(xs-(x1+h/2*k11)-a1*k12/w0-a2*k22/w0^2-a3*k32/w0^3-a4*k42 /w0^4-a5*k52/w0^5); %%ki3(i=1,2,...6) k13=x2+h/2*k22; k23=x3+h/2*k32; k33=x4+h/2*k42; k43=x5+h/2*k52; k53=x6+h/2*k62; k63=w0^6/a6*(xs-(x1+h/2*k12)-a1*k13/w0-a2*k23/w0^2-a3*k33/w0^3-a4*k43 /w0^4-a5*k53/w0^5); %ki4(i=1,2,...6) k14=x2+h*k23; k24=x3+h*k33; k34=x4+h*k43;

飞行器轨迹仿真大作业

飞行器轨迹仿真大作业 [文档副标题]

1. 积分方法 1） 欧拉积分 1(,)n n n n x x h f x t +=+? 2） 二阶龙格库塔 112121() 2 (,) (,)n n n n n n h x x k k k f t x k f t h x hk +=++==++ 3） 四阶龙格库塔 1123461212232224 3(22)(,) (,) (,) (,) h i i i i h h i i h h i i i i x x K K K K K f t x K f t x K K f t x K K f t h x hK +=++++==++=++= ++ 2. 计算例子 实现自由落体运动的轨迹求解。 (,)h v D v h v g m ?=??=--?? 00100/,10000,5v m s h m m kg ===，阻力系数0.02d C =，小球半径0.1r m =, 209.81/g m s =，重力加速度取常数，大气密度采用插值计算。求落地速度和落地时间，观察速度的变化规律。 大气密度变化规律： [高度/M] 0.0 1000 3000 5000 7000 9000 11000 15000 20000 [密度/kg/m^3] 1.225 1.112 0.909 0.736 0.589 0.466 0.364 0.194 0.088

3. 图形 3.1 用Euler 法画出的图 图1 速度随时间的变化图 图2 高度随时间的变化 3.2 用Runge-Kutta 法画的图形 时间t 速度v 01020 30405060 时间t 高度h

飞行器CFD仿真分析专题培训

飞行器CFD仿真分析专题培训 应群里兄弟对飞行器的热爱本部组织举行飞行器空气动力学流体CFD专题培训，这是飞机设计中密不可少，鉴于目前关于此培训市场上没有，而专题的针对型培训价格昂贵，所以组织此培训来满足大家的需求。 培训老师：张永立，原ANSYS中国西部大区流体技术经理，10年以上CFD 工程仿真经验。 课程分三大节，主要内容是飞行器CFD仿真分析，包括模型处理、网格划分、计算设置与求解、后处理等内容。课程时间为8月11日到8月15日，晚上8点开始。 通过这次培训可以让初学者基本掌握Fluent飞行器外流场气动分析的基本知识，为后续高级课程打好基础。 为了确保培训效果，提供培训讲解视频。 第一大节课： 1.CFD理论基础简介 2.CFD仿真分析流程及要点说明 3.演示案例1：FLUENT翼型外流场绕流流体动力学分析（2D） a)DesignModeler几何建模 b)Meshing网格划分 c)FLUENT物理前处理设置 d)FLUENT求解及监控 e)CFD-Post计算结果后处理 第二大节课： 1.FLUENT物理前处理设置技术详解 a)材料定义 b)计算区域/计算模型定义 c)边界条件定义 2.FLUENT求解技术讲解 a)求解器类型 b)初始条件 c)数据监控/收敛准则 3.CFD-POST结果后处理技术详解 4.演示案例2：FLUENT翼身结构外流场气动分析 a)导入已有网格 b)定义计算模型和边界条件

c)求解监控和收敛准则 d)初始化与求解 e)后处理（气动载荷分布、载荷积分、力矩计算、翼型升阻力特性处理等） 第三大节课： 1.提高飞行器气动分析精度的主要途径 2.ICEM CFD结构化网格讲解 d)ICEM CFD基本介绍 e)六面体网格block的操作理念 f)六面体网格划分操作讲解 3.案例演示3：弹体外流场气动分析 d)ICEM CFD结构化网格划分和导出 e)FLUENT物理前处理设置 f)FLUENT求解及监控 g)CFD-Post计算结果后处理