北航大型通用软件
- 格式:doc
- 大小:1002.50 KB
- 文档页数:16
一、国外软件1.Unigra phics(UG) UG是Uni graphi csSol utio n s公司的拳头产品。
该公司首次突破传统CA D/CAM模式,为用户提供一个全面的产品建模系统。
在U G中,优越的参数化和变量化技术与传统的实体、线框和表面功能结合在一起,这一结合被实践证明是强有力的,并被大多数C A D/CAM软件厂商所采用。
U G最早应用于美国麦道飞机公司。
它是从二维绘图、数控加工编程、曲面造型等功能发展起来的软件。
90年代初,美国通用汽车公司选中U G 作为全公司的CAD/CAE/CAM/CIM主导系统,这进一步推动了U G的发展。
1997年10月Uni graphi csSol utio n s公司与Intergraph公司签约,合并了后者的机械CAD产品,将微机版的S O LIDEDGE软件统一到Pa rasol i d平台上。
由此形成了一个从低端到高端,兼有Unix工作站版和W indo wsNT微机版的较完善的企业级CAD/CAE/CAM/PDM集成系统。
2.SOLIDE DGE SOLIDE DGE是真正Wind ows软件。
它不是将工作站软件生硬地搬到Wi ndows平台上,而是充分利用W i ndo ws基于组件对象模型(COM)的先进技术重写代码。
SOLIDE DGE与M i c rosoftOffi ce兼容,与Wind o ws的OL E技术兼容,这使得设计师们在使用C A D系统时,能够进行Wi ndows下字处理、电子报表、数据库操作等。
SOLIDE DGE具有友好的用户界面,它采用一种称为SmartRibb on的界面技术,用户只要按下一个命令按钮,既可以在Sm artRi bbon上看到该命令的具体的内容和详细的步骤,同时在状态条上提示用户下一步该做什么。
航空科学与工程学院《航空工程大型通用软件应用》大作业机翼结构设计与分析组号第3组小组成员11051090 赵雅甜11051093 廉佳11051100 王守财11051108 刘哲11051135 张雄健11051136 姜南6月目录一 CATIA部分....................................... 错误!未定义书签。
( 一) 作业要求..................................... 错误!未定义书签。
( 二) 作业报告..................................... 错误!未定义书签。
1、三维模型图................................... 错误!未定义书签。
2、工程图....................................... 错误!未定义书签。
二 FLUENT部分...................................... 错误!未定义书签。
( 一) 作业要求..................................... 错误!未定义书签。
( 二) 作业报告..................................... 错误!未定义书签。
1、计算方法和流程............................... 错误!未定义书签。
2、网格分布图................................... 错误!未定义书签。
3、气动力系数................................... 错误!未定义书签。
4、翼型表面压力曲线............................. 错误!未定义书签。
5、翼型周围压力云图............................. 错误!未定义书签。
北航MOOC系统Android客户端NABC北航MOOC⼿机客户端NABC分析1) N (Need 需求) MOOC是Massive Open Online Course的缩写,通常被译为⼤型开放式⽹络课程,它最早在08年的时候由⼀位加拿⼤的学者提出,经过了⼏年的沉寂,终于在2012年随着⼏⼤课程平台:Coursera、edX、Udacity的崛起⽽迅速成为全球最为热门的教育话题之⼀。
MOOC这是⼀种在线教育形式,任何⼈都能免费注册使⽤,它有与线下课程类似的作业评估体系和考核⽅式。
按时完成作业和考试的学习者还可能收获导师签署的课程证书。
可以说,MOOC技术是在⾼等教育资源分布不均的情况下,由计算机技术的成熟⽽推动的教育领域的⼀项重⼤变⾰。
这项技术的⽬的在于使更多的在现实⽣活中没有办法接受⾼等教育的⼈享受到和顶级⾼校学⽣⼀样的教育资源,同时也能够让⾼校的教育资源得到更充分的利⽤。
随着技术的发展和推⼴,这场由在线学习掀起的⾰命,可能将会重塑⾼等教育的版图。
⽽在中国,这个⾼等教育资源分配不均的问题尤为严重,接受到中国顶尖⾼等教育的⼈只是很少的⼀部分。
据统计,课程平台Coursera 的注册学⽣其中约41000⼈来⾃中国,占总⼈数的40%左右。
这些接受在线教育的⼈群中,有些⽤户是学⽣,希望利⽤外部资源补充他们的课堂作业。
有些⽤户是上班族,希望助⼒他们的职业发展。
有些⽤户是找⼯作的⼈,希望获得“成就报告”或“验证证书”以改善他们的简历。
还有些⽤户是终⾝学习者,只是喜欢增长⾃⼰的知识。
由此可见,国内对在线教育的需求还是相当⼤的,所以中国⼈发展⾃⼰的MOOC课程⽹络已经是⼤势所趋。
北京航空航天⼤学作为中国顶尖⾼校之⼀,⾃然有必要发展⾃⼰学校的在线教育,将北航校内的精品课程推⼴开来,既可以提升北航的学校知名度,⼜能使本校的教育资源被更多的⼈获得,何乐⽽不为呢?本次的MOOC⼿机客户端是基于学长在IOS平台上的开发基础,我们将程序移植到Android平台,同时对功能进⾏更进⼀步的完善,使更多的⽤户接触到这个客户端并从中获益。
北京航空航天大学科技成果——机器人离线编程系统
(Robot CAX)
成果简介Robot CAX利用计算机图形学,在计算机上建立机器人及其工作环境的模型,开发规划算法,通过对模型的控制和操作,对机器人进行轨迹规划,生成机器人控制程序。
Robot CAX使用通用CAD/CAM系统(如CATIA、Master CAM等)生成的G代码或APT代码作为加工轨迹进行输入。
获取加工轨迹之后,Robot CAX进行运动仿真、碰撞检查、代码优化等操作,以校验出机器人加工的正确性和可达性。
在这个工作程序中,省略了示教过程,机器人无需停工,因此可以节省大量的加工时间。
同时,该系统还可以自由定义末端执行器、工装、导轨、旋转台等其他外围设备。
仿真优化完成后,Robot CAX直接输出优化后的机器人控制代码,进而控制机器人进行实际加工。
可实现机器人及外围环境的三维可视化仿真,三维环境的运动路径仿真,机器人的运动学和动力学仿真,运动路径的友好人机交互,机器人关节轴超限检查,利用机器人冗余自由度进行运动优化,奇异点自动回避,基于B样条的运动路径插值技术,G代码、APT与机器人语言之间的转换,兼容多种工业机器人,兼容多种CAD/CAM系统。
Robot CAX软件集成了完全自主知识产权的三维内核,是一个独立于特定CAD/CAM系统的软件程序。
因此,用户并不受制于具体CAD/CAM系统或其版本,可以根据具体情况,配套使用任意的CAD/CAM系统。
CAM软件CAM (computer Aided Manufacturing,计算机辅助制造)的核心是计算机数值控制(简称数控),是将计算机应用于制造生产过程的过程或系统。
CAM软件是具有CAM功能的软件的统称,常见的有CAXA,UG,Hypermill等目录等多道工序,这些都是通过程序指令控制运作的,只要改变程序指令就可改变加工过程,数控的这种加工灵活性称之为“柔性”。
CAM(computer Aided Manufacturing,计算机辅助制造):利用计算机来进行生产设备管理控制和操作的过程。
它输入信息是零件的工艺路线和工序内容,输出信息是刀具加工时的运动轨迹(刀位文件)和数控程序。
CAM (computer Aided Manufacturing,计算机辅助制造)的核心是计算机数值控制(简称数控),是将计算机应用于制造生产过程的过程或系统。
1952年美国麻省理工学院首先研制成数控铣床。
数控的特征是由编码在穿孔纸带上的程序指令来控制机床。
此后发展了一系列的数控机床,包括称为“加工中心”的多功能机床,能从刀库中自动换刀和自动转换工作位置,能连续完成锐、钻、饺、攻丝等多道工序,这些都是通过程序指令控制运作的,只要改变程序指令就可改变加工过程,数控的这种加工灵活性称之为“柔性”。
加工程序的编制不但需要相当多的人工,而且容易出错,最早的CAM便是计算机辅助加工零件编程工作。
麻省理工学院于1950年研究开发数控机床的加工零件编程语言APT,它是类似FORTRAN的高级语言。
增强了几何定义、刀具运动等语句,应用APT使编写程序变得简单。
这种计算机辅助编程是批处理的。
CAM系统一般具有数据转换和过程自动化两方面的功能。
CAM所涉及的范围,包括计算机数控,计算机辅助过程设计。
市面上的CAM软件有:UG NX、Pro/NC、CATIA、MasterCAM、SurfCAM、SPACE-E、CAMWORKS、WorkNC、TEBIS、HyperMILL、Powermill、Gibbs CAM、FEATURECAM、topsolid、solidcam、cimtron、vx、esprit、gibbscam、Edgecam......等等数控除了在机床应用以外,还广泛地用于其它各种设备的控制,如冲压机、火焰或等离子弧切割、激光束加工、自动绘图仪、焊接机、装配机、检查机、自动编织机、电脑绣花和服装裁剪等,成为各个相应行业CAM的基础。
CAM软件分类及功能概述国产CAM软件CAXA制造工程师目前,国产CAM软件只有一款,就是CAXA。
是北京数码大方科技有限公司(CAXA),北航海尔自主研发的CAD/CAM 软件,主要以教育行业为主,并且是数控大赛指定参赛软件。
Computer Aided X Alliance -Always a stepAhead ”“X:technology,product,solution and service …”联盟合作的领先一步的计算机辅助技术与服务。
拥有完全自主知识产权的系列化,CAD、CAPP、CAM、DNC、EDM、PDM、MES、MPM等PLM软件产品和解决方案。
覆盖了制造业信息化设计、工艺、制造和管理四大领域。
产品广泛应用于装备制造、电子电器、汽车、国防军工、航空航天、工程建设、教育等各个行业。
国外CAM软件Ugnx(UG);Mastercam;Solidcam;Topsolid;Cimtron;Procam(ptc);Powermill(Delcam);Hypermill(openmind);Worknc;Esprit;Gibbscam;surfcam;Edgecam...............西门子UGS:UG NX最新版本:NX7强大的造型能力和数控编程能力,几乎所有的飞机发动机和汽车发动机都是UG设计,军工领域等高端工程领域运用广泛,高端领域与catia并驾齐驱;Procam多轴车削,冲加工,等离子/激光切割,及电火花加工等。
英国达尔康Delcam:Powermill,Featurecam,Partmaker,ArtcamPowermill 2-5轴(最新版本10版本);FeatureCAM - 产品加工、车铣复合、线切割;PartMaker - 瑞士型纵切机床、车削中心编程;ArtCAM - 立体艺术浮雕CAD/CAM系统;DentMILL - 牙科专业CAM加工系统。
德国Openmind:Hypermill致力于叶轮叶片的加工模块美国CNC:Mastercam实现功能:铣切,线切割,车削,实体,浮雕,木雕Mastercam模块:design设计,lathe车,mill铣,wire 线切割:模块构成:Mastercam 铣削;Mastercam 线切割;Mastercam 车削;Mastercam 木工铣削;Mastercam in Solidworks;Mastercam 三维雕刻;Mastercam实体设计;最新版本:mastercam X4 MU1(X就是10的意思)以色列:Cimtron,GibbsCAMcimtron模具用户运用广泛。
/m-selig/ads/coord_database.html (超级牛的网站)/m-selig/ads.html以下介绍航天相关的软件,主要是飞机设计领域的相关软件a. 飞机外型设计1. Raymer Simplified Aircraft Design Spreadsheet (Excel Shareware)网址:/ac-size.html2. Design Calc.xlshttp://www.hoppenbrouwer-home.nl/ikarus/software/design.htm3. Design Analysis Program网址:/articles/design/helmutlelkedesignanalysis.htmb. 气动力分析1. TornadoTornado是利用Vortex lattice method来模拟飞机外型的各种气动力参数。
输入飞机每一段机翼的翼形、与原点的相对位置等相关参数后,就可以计算出飞机所受到的力与力矩,以及相关的气动力参数。
不过有个蛮大的缺点是Tornado只能使用NACA系列的翼剖面,NACA以外的就无法计算。
另外,如果机身的截面不是NACA翼形的话,机身对气动力参数的影响也无法列入考虑,因此无法准确估计飞机的横向气动力参数首页:http://www.ave.kth.se/divisions/aero/software/tornado/2. AVLAerodynnamic analysis, trim calculation, dynamic stability analysis,aircraft configuraiotn analysis/drela/Public/web/avl/3. Software for Aerodynamics and Aircraft Design/~mason/Mason_f/MRsoft.html4. CFD codes list网址:/cfd/CFD_codes_p.htmlc. 翼形(airfoil)分析1. UIUC airfoil viewerUIUC是很有名的大站,提供各种翼形点的数据viewer 下载:/airfoil/index.html2. Xfoil有名的翼形分析软件,可以输入翼形的点资料以及飞行的雷诺数,算出升力系数等相关参数 DOS界面,比较不漂亮/drela/Public/web/xfoil/3. Javafoil功能同上,有完善的图形使用者接口,很容易使用首页:http://www.mh-aerotools.de/airfoils/javafoil.htm下载:http://www.mh-aerotools.de/airfoils/java/javafoilinstaller.msi说明书:http://www.mh-aerotools.de/airfoils/jf_users_manual.htm有的计算机必须安装Java虚拟机器/zh_TW/download/index.jsp安装完后,只要执行开始菜单内的Javafoil-VBScript就可以使用!4. Profili同样是翼形分析软件,不过有输出不同弦长翼形的功能,对制作飞机的人来说很有用注册版有比较多功能,不过免费版一般说来还很够用首页:/eng/default.htm下载:/download/profiliv2.exe其它相关的翼形分析软件可以到以下的网站找到/m-selig/ads.html翼形数据库1. UIUC/m-selig/ads/coord_database.html2. NASG这网站的参与人员参与了十年的鸟人大赛,致力于发展人力飞机,所以他们也有一系列的翼形数据库/afdb/list-airfoil-e.phtml3. Winfoil/airfoildb/airfoil.aspd. 推进系统1. Javaprop螺旋桨设计软件,可以指定翼形等相关参数,再计算出螺旋桨的性能首页:http://www.mh-aerotools.de/airfoils/javaprop.htm下载:http://www.mh-aerotools.de/airfoils/java/javapropinstaller.msi在线使用版:http://www.mh-aerotools.de/airfoils/javaprop.htm同样的可能要安装Java虚拟机器/安装完后,只要执行开始菜单内的Javaprop Script-VBScript就可以使用!写Javafoil和Javaprop的都是同一个人-Martin Hepperle,在他的网站http://www.mh-aerotools.de/airfoils/index.htm可以发现很多和模型飞机空气动力学有关的知识2. Propeller Performance Calculator简单好用的软件,但算出来的结果仅供参考/_frames/f_propcalc.shtml下载:/_frames/f_downloads.shtml3. Electric Motor Performance Calculator/_frames/f_emotor.shtml下载:/_frames/f_downloads.shtml4. Qprop分析螺旋桨和电动马达性能的程序,不容易上手/drela/Public/web/qprop/5. ThrustHP螺旋桨飞机推力分析软件,分析结果看看就好http://www.hoppenbrouwer-home.nl/ikarus/software/thrusthpv20d.htm6. PropellerScanner只要用扫描仪把螺旋桨的上视图及侧视图扫成图档,就可以利用Javaprop分析螺旋桨的外型及twist资料http://www.mh-aerotools.de/airfoils/index.htme. 飞机动态仿真软件1. AerosimAerosim是一个MATLAB / Simulink 的toolbox,可以对飞机的六自由度非线性模型进行模拟,还有很多控制器设计的相关模块。
航空科学与工程学院《航空工程大型通用软件应用》大作业组号47小组成员11051165 王海涛11051196 张彧11051210 杜酉昌11051217 郭皓11051218 张启冬2014年6月小组成员分工表目录1.Catia部分 (3)1.1.作业要求: (3)1.2.作业报告: .............................................. 错误!未定义书签。
2.Fluent部分 (3)2.1.作业要求: (4)2.2.作业报告: .............................................. 错误!未定义书签。
3.Ansys部分......................................................... 错误!未定义书签。
3.1.作业要求: (7)3.2.作业报告: .............................................. 错误!未定义书签。
4.Matlab部分 ....................................................... 错误!未定义书签。
4.1.作业要求: .............................................. 错误!未定义书签。
4.2.作业报告: .............................................. 错误!未定义书签。
1.Catia部分:1.1作业要求:1、用提供的翼型数据,生成翼型曲线;2、参考所给的机翼图例,使用参数化设计,参数为:翼根弦长=495mm、翼中弦长=270mm、翼尖弦长=150mm、内翼展长=525mm、外翼展长=800mm、内翼后掠角=10度、外翼后掠角=10度翼尖不要求倒圆,最后在Part设计环境中生成实心体机翼(不是曲面);3、在工程图环境中生成三视图,并标注尺寸;4、上述三维模型和工程图分别截图插入word文档的Catia部分;最后提交作业时,part文件和word报告一起提交。
航空科学与工程学院《航空工程大型通用软件应用》大作业组号47小组成员11051165 王海涛11051196 张彧11051210 杜酉昌11051217 郭皓11051218 张启冬2014年6月小组成员分工表目录1.Catia部分 (3)1.1.作业要求: (3)1.2.作业报告: .............................................. 错误!未定义书签。
2.Fluent部分 (3)2.1.作业要求: (4)2.2.作业报告: .............................................. 错误!未定义书签。
3.Ansys部分......................................................... 错误!未定义书签。
3.1.作业要求: (7)3.2.作业报告: .............................................. 错误!未定义书签。
4.Matlab部分 ....................................................... 错误!未定义书签。
4.1.作业要求: .............................................. 错误!未定义书签。
4.2.作业报告: .............................................. 错误!未定义书签。
1.Catia部分:1.1作业要求:1、用提供的翼型数据,生成翼型曲线;2、参考所给的机翼图例,使用参数化设计,参数为:翼根弦长=495mm、翼中弦长=270mm、翼尖弦长=150mm、内翼展长=525mm、外翼展长=800mm、内翼后掠角=10度、外翼后掠角=10度翼尖不要求倒圆,最后在Part设计环境中生成实心体机翼(不是曲面);3、在工程图环境中生成三视图,并标注尺寸;4、上述三维模型和工程图分别截图插入word文档的Catia部分;最后提交作业时,part文件和word报告一起提交。
5、生成的Catia实体机翼模型会在后续作业中使用。
1.2作业报告:1、作业截图如下:图错误!文档中没有指定样式的文字。
-1 三维模型截图图错误!文档中没有指定样式的文字。
-2 工程图截图2、part文件见电子版附件。
2.Fluent部分2.1作业要求:1、从Catia软件中输出igs文件,导入到Gambit软件中;2、用与机翼对称面相距Z的平面撕裂机翼表面生成计算所需的翼型曲线;3、采用分区的四边形Map网格或者应用尺寸函数的非结构网格生成计算网格;4、计算条件:压力远场边界条件、SA湍流模型;5、计算结果提取内容:网格分布图、气动力系数、翼型表面压力曲线、翼型周围的压力云图和速度云图、翼型周围的速度矢量图和流线图。
6、word报告:计算方法和流程、计算结果提取内容。
2.2计算流程1.导入GAMBIT所画的default_id7076网格数据;2.对读入的网格进行检查;3.选择求解器为压力基;4.选择物理模型,激活能量方程,并选择SA湍流模型;5.定义流体属性,将空气密度项改成由理想气体状态方程确定;6.定义操作条件;7.定义边界条件,定义各部分对应的边界条件,其中压力远场边界条件中马赫数为0.3,X 和Y方向流速分量分别为cos2°(0.999391)和sin2°(0.0348995);8.设置参考值,将其中的面积和长度项都改为1,计算起始位置选择边界条件为压力远场处;9.求解方法选择,库朗数及松弛因子均选用程序默认值;10.设置收敛监视,将残差各项收敛条件值均改为0.0001;11.定义初始条件,选择计算开始处为压力远场对应边界,然后点击初始化按钮;12.设置自动存储,设置每200布进行自动存储13.迭代计算;14.计算300步后收敛,保存结果。
2.3网格分布图图错误!文档中没有指定样式的文字。
-3 网格分布图2.4空气动力系数cl=0.22854028 cd=0.00581640192.5翼型表面压力曲线2.6 翼型周围的压力云图和速度云图2.7 翼型周围的速度矢量图2.8 迭代过程(相关文件见电子版附件)3.Ansys部分3.1.作业要求:1、将Catia模型导入Ansys生成实体模型;2、采用扫略网格划分计算模型;3、输入参数:机翼为各向同性线弹性材料,弹性模量90GPa,泊松比0.2,材料的拉伸强度为120MPa,压缩强度为80MPa;4、计算条件:翼根为固支边界条件;机翼上、下表面分布承受0.1MPa、0.5MPa均布气动压力,翼尖中点上承受700N集中力作用,方向与上表面压力方向相同;5、计算要求:计算机翼的全场应力和变形,并对机翼的强度进行校核,给出由翼根到翼尖任意路径上的应力和挠度分布规律;计算机翼的固有振动频率6、提交结果:计算报告和命令流。
3.2作业报告:3.2.2 第二强度理论应力云图(S2)3.2.4 应变分布云图3.2.5 机翼变形示意图3.2.6 应力沿某一路径分布示意图3.2.7 挠度曲线3.2.8 振动频率3.2.9命令流见附录4.Matlab部分:4.1作业要求:从Fluent计算结果中导出数据:1.在三维空间绘制翼型外形压力系数分布图;2.计算整个翼型沿X方向的阻力;3.计算整个翼型沿Y方向的升力;4.计算整个翼型相对于1/4弦线长位置的俯仰力矩5.提交计算报告:命令语句及计算结果4.2作业报告图错误!文档中没有指定样式的文字。
-4 翼型外形压力系数分布图命令语句function finalread;%读取数据graph;%在三维空间绘制翼型外形压力系数分布图returnfunction readglobal pre_up pre_down upperfoil lowerfoilpre_up=importdata('file_up.txt');pre_down=importdata('file_down.txt');upperfoil=importdata('upfoil.txt');lowerfoil=importdata('downfoil.txt');returnfunction graph%在k中分别存储翼型的x,y,压力值global k pre_up pre_down upperfoil lowerfoil[num_node,dummy] = size(pre_up) ;i=num_node;ii=i;[num_node,dummy] = size(pre_down) ;i=i+num_node;K=sparse(i,3) ;for ie=1:1:iif(ie<=ii)k(ie,1) = pre_up(ie,1);j=(k(ie,1)+0.242711)/0.41*100;[num_node,dummy] = size(upperfoil) ;for ig=1:1:num_nodeif(j<upperfoil(ig,1))k(ie,2) = (upperfoil(ig-1,2)+(j-upperfoil(ig-1,1))/(upperfoil(ig,1)-upperfoil(ig-1,1))*(upperfoil(ig,2)-upperfo il(ig-1,2)))/100*0.41;break;endendk(ie,3) = pre_up(ie,2)/2628.1;elsek(ie,1) = pre_down((ie-ii),1);j=(k(ie,1)+0.242711)/0.41*100;[num_node,dummy] = size(lowerfoil) ;for ig=1:1:num_nodeif(j<lowerfoil(ig,1))k(ie,2) = (lowerfoil(ig-1,2)+(j-lowerfoil(ig-1,1))/(lowerfoil(ig,1)-lowerfoil(ig-1,1))*(lowerfoil(ig,2)-lowerfo il(ig-1,2)))/100*0.41;break;endendk(ie,3) = pre_down((ie-ii),2)/2628.1;endend%在三维空间绘制翼型外形压力系数分布图for ie=1:1:iplot3(k(ie,1),k(ie,2),k(ie,3),'r-*')hold onplot3(k(ie,1),k(ie,2),0,'k.')hold online([k(ie,1) k(ie,1)],[k(ie,2) k(ie,2)],[k(ie,3) 0],'LineStyle','-')hold onendhold off%计算整个翼型沿X方向的阻力、升力及整个翼型相对于1/4弦线长位置的俯仰力矩dra=0lift=0;mom=0;for ie=1:1:(ii-1)if(k(ie,1)==k(ie+1,1)) continue;endp1 = [k(ie,1) k(ie,2)];p2 = [k(ie+1,1) k(ie+1,2)];force=(k(ie,3)+k(ie+1,3))/2*2628.1*norm(p1-p2);dra=dra+force*(k(ie,2)-k(ie+1,2))/norm(p1-p2);lift=lift+force*(k(ie,1)-k(ie+1,1))/norm(p1-p2);mom=mom+force*(k(ie,2)-k(ie+1,2))/norm(p1-p2)*(k(ie,2)+k(ie+1,2))/2+force*(k(ie,1)-k(ie+1,1) )/norm(p1-p2)*(-0.1302-(k(ie,1)+k(ie+1,1))/2);endfor ie=(i-ii):1:(i-1)if(k(ie,1)==k(ie+1,1)) continue;endp1 = [k(ie,1) k(ie,2)];p2 = [k(ie+1,1) k(ie+1,2)];force=(k(ie,3)+k(ie+1,3))/2*2628.1*norm(p1-p2);dra=dra+force*(k(ie,2)-k(ie+1,2))/norm(p1-p2);lift=lift+force*(k(ie,1)-k(ie+1,1))/norm(p1-p2);mom=mom+force*(k(ie,2)-k(ie+1,2))/norm(p1-p2)*(k(ie,2)+k(ie+1,2))/2+force*(k(ie,1)-k(ie+1,1) )/norm(p1-p2)*(-0.1302-(k(ie,1)+k(ie+1,1))/2);end%显示结果fprintf( ' 阻力升力力矩\n' ) ;fprintf( '%16.8e N %16.8e N %16.8e N/m\n',dra,lift,mom) ;return最终结果:曲线图见上截图X方向上的阻力:1.51073775e+02 NY方向上的升力:1.30671697e+03 N1/4弦线长处的俯仰力矩:2.5315138e+01 N/m。