CN106709174A - 一种基于有限元模型的飞机活动面偏转方法 - Google Patents

一种基于有限元模型的飞机活动面偏转方法 Download PDF

Info

Publication number
CN106709174A
CN106709174A CN201611161877.4A CN201611161877A CN106709174A CN 106709174 A CN106709174 A CN 106709174A CN 201611161877 A CN201611161877 A CN 201611161877A CN 106709174 A CN106709174 A CN 106709174A
Authority
CN
China
Prior art keywords
deflection
coordinate system
active face
aircraft
axis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201611161877.4A
Other languages
English (en)
Other versions
CN106709174B (zh
Inventor
李苗
陈军
李军贵
郑洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201611161877.4A priority Critical patent/CN106709174B/zh
Publication of CN106709174A publication Critical patent/CN106709174A/zh
Application granted granted Critical
Publication of CN106709174B publication Critical patent/CN106709174B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明属于飞机强度计算领域,涉及一种基于有限元模型的飞机活动面偏转方法。所述的方法包括如下步骤:步骤一,在基本坐标系下建立中立位置的活动面有限元模型;步骤二,在中立位置的活动面有限元模型中建立基准偏转坐标系,坐标系原点为某一悬挂支座的耳片孔心,X轴垂直于活动面前梁平面,前缘方向为正,Z轴为活动面转轴,翼尖方向为正,Y轴由右手法则确定;步骤三,将中立位置活动面有限元模型所有节点的坐标定义在基准偏转坐标系下;步骤四,将基准偏转坐标系绕其Z轴偏转要求的角度生成偏转坐标系,实现了飞机活动面的偏转。提供了一种基于有限元模型的飞机活动面偏转方法。

Description

一种基于有限元模型的飞机活动面偏转方法
技术领域
本发明属于飞机强度计算领域,涉及一种基于有限元模型的飞机活动面偏转方法。
背景技术
飞机强度分析中,针对活动面具有多偏度受载的问题,传统的解决办法是只建立活动面的中立位置有限元模型,对于带偏角的载荷情况采用作动器相对活动面偏转的方法来模拟因活动面偏转带来的操纵力臂的变化。对于舵面的单独求解该方法可行,但是需要将活动面与盒段连接求解来模拟更真实的支持刚度时,该方法无法实现多载荷工况下活动面不同偏转角度的受载状态。
发明内容
本发明的目的:提供一种基于有限元模型的飞机活动面偏转方法。
本发明的技术方案:一种基于有限元模型的飞机活动面偏转方法,其特征在于所述的方法包括如下步骤:
步骤一,在基本坐标系下建立中立位置的活动面有限元模型;
步骤二,在中立位置的活动面有限元模型中建立基准偏转坐标系,坐标系原点为某一悬挂支座的耳片孔心,X轴垂直于活动面前梁平面,前缘方向为正,Z轴为活动面转轴,翼尖方向为正,Y轴由右手法则确定;
步骤三,将中立位置活动面有限元模型所有节点的坐标定义在基准偏转坐标系下;
步骤四,将基准偏转坐标系绕其Z轴偏转要求的角度生成偏转坐标系,实现了飞机活动面的偏转。
优选地,基准偏转坐标系绕其Z轴偏转的角度对应于飞机的载荷情况。
优选地,采用PATRAN软件的PCL语言实现飞机活动面的偏转。
优选地,在基于PATRAN软件实现飞机活动面的偏转系统内预置有每种飞机载荷情况下活动面对应的偏转角度。
优选地,对于多种载荷情况,各载荷情况偏转坐标系输出到不同的偏转坐标系文件,文件名称采用载荷情况号。
本发明的有益效果:针对活动面具有多偏度受载的问题,传统的解决办法是只建立活动面的中立位置有限元模型,对于带偏角的载荷情况采用作动器相对活动面偏转的方法来模拟因活动面偏转带来的操纵力臂的变化。对于舵面的单独求解该方法可行,但是需要将活动面与盒段连接求解来模拟更真实的支持刚度时,该方法无法实现多载荷工况下活动面不同偏转角度的受载状态。本专利建立了以活动面转轴为坐标轴的基准偏转坐标系,活动面的坐标都是在基准偏转坐标系下给出,活动面需要偏转角度时将基准偏转坐标系绕转轴转动相应角度来得到偏转角度下的坐标系,便实现了整个活动面的偏转,符合实际结构受力状态,使得活动面的有限元分析结果更准确。同时采用PATRAN软件的PCL语言,编写了活动面结构参数化建模程序,实现了活动面结构建模的参数化和自动化以及模型偏转的批处理,实现了活动面与盒段连接求解时所有偏角工况下的模型快速偏转。
附图说明
图1升降舵有限元模型示意图;
图2升降舵基准偏转坐标系定义示意图;
图3升降舵参数化建模程序流程图;
图4基准偏转坐标系下定义的节点坐标格式示意图;
图5升降舵基准偏转坐标系偏转前后对比示意图;
图6升降舵有限元模型中立位置时与平尾相对位置示意图;
图7升降舵有限元模型偏转后与平尾相对位置示意图;
图8升降舵各载荷情况对应偏转角度文件格式示意图;
图9导入程序文件示意图;
图10程序化生成偏转坐标系文件格式示意图。
具体实施方式
以某型飞机升降舵为例来介绍一下飞机活动面偏转的具体操作流程:
步骤1、在基本坐标系即飞机总体坐标系下,建立升降舵中立位置的有限元模型,如图1所示基本坐标系1,左升降舵中立位置有限元模型2、右升降舵中立位置有限元模型3,升降舵中立位置即升降舵偏转角度为0度时的位置。该步骤中将常规有限元建模方法进行程序化,避免了手动生成大量节点和单元的繁琐工作。如图3所示为升降参数化建模程序流程图,其中的升降舵偏转模块为接下来重点介绍的部分。
步骤2、建立升降舵基准偏转坐标系,坐标系编号定义为73000,如图2所示。选取升降舵靠近翼根的1#悬挂支座6耳孔中心为坐标系原点,X轴垂直于升降舵前梁平面7,指向前缘方向5为正,Z轴为升降舵转轴8,指向翼尖方向4为正。
步骤3、将升降舵有限元模型的节点坐标定义在基准偏转坐标系下,即所有节点坐标值均是相对于基准偏转坐标系原点的位置来给出的。如图4所示为升降舵有限元模型节点在基准偏转坐标系下定义的输出文件格式示意图。
步骤4、将基准偏转坐标系绕其Z轴(即升降舵转轴)偏转一定的角度后生成偏转坐标系,因升降舵有限元模型节点坐标相对于基准偏转坐标系位置固定不变,因此坐标系偏转后,升降舵模型也跟着偏转了相同的角度。如图5所示为升降舵基准偏转坐标系偏转前后对比示意图,如图6、图7给出了偏转前后升降舵有限元模型相对于平尾盒段的位置示意图。升降舵在各载荷情况下偏转角度不同,在有限元分析计算时需建立各偏转角度下升降舵的有限元模型。而且结构设计往往进行多轮次载荷求解,且各轮载荷情况多达上百种、偏转角度也不同,因此对应每轮次载荷的求解要重新建立偏转模型,工作量大。因此本步骤利用PCL语言进行批处理偏转坐标系文件,即快速输出多个偏转坐标系文件,每个偏转坐标系文件与模型文件组合起来便为一个偏角下的模型。当载荷情况变化时,仅改变偏转坐标系文件即可,模型文件不变,大大减小了建模的工作量。
PCL命令语言是集成于MSC.PATRAN中的一个模块化高级编程语言和用户自定义工具,可供用户对PATRAN进行二次开发或编制专用程序,它的主要功能有:生成能直接从PATRAN调用的函数;生成窗体和控键;调用PATRAN的内部函数和操作其数据库;调用PATRAN外的可执行文件。这些特点使用户能灵活地调配PATRAN现有的功能和添加新的函数进行参数化建模。PCL语言类似于C语言,具有标准C语言的几乎全部功能,简单易用,适合于工程技术人员进行程序设计。
PCL程序可在文本文件中编写,文件后缀名为.PCL,其输入文件为升降舵各载荷情况对应偏转角度文件,如图8所示输入文件格式。其运行方法是在PATRAN软件的命令窗口中导入程序文件,如图9所示,导入命令为!!I NPUT。运行完成后即可生成不同载荷情况对应的偏转坐标系文件,每个载荷情况对应一个坐标系文件,文件名称以载荷情况号命名,输出偏转坐标系文件格式如图10所示。

Claims (5)

1.一种基于有限元模型的飞机活动面偏转方法,其特征在于所述的方法包括如下步骤:
步骤一,在基本坐标系下建立中立位置的活动面有限元模型;
步骤二,在中立位置的活动面有限元模型中建立基准偏转坐标系,坐标系原点为某一悬挂支座的耳片孔心,X轴垂直于活动面前梁平面,前缘方向为正,Z轴为活动面转轴,翼尖方向为正,Y轴由右手法则确定;
步骤三,将中立位置活动面有限元模型所有节点的坐标定义在基准偏转坐标系下;
步骤四,将基准偏转坐标系绕其Z轴偏转要求的角度生成偏转坐标系,实现了飞机活动面的偏转。
2.根据权利要求1所述的一种基于有限元模型的飞机活动面偏转方法,其特征为:基准偏转坐标系绕其Z轴偏转的角度对应于飞机的载荷情况。
3.根据权利要求1所述的一种基于有限元模型的飞机活动面偏转方法,其特征为:采用PATRAN软件的PCL语言实现飞机活动面的偏转。
4.根据权利要求3所述的一种基于有限元模型的飞机活动面偏转方法,其特征为:在基于PATRAN软件实现飞机活动面的偏转系统内预置有每种飞机载荷情况下活动面对应的偏转角度。
5.根据权利要求3所述的一种基于有限元模型的飞机活动面偏转方法,其特征为:对于多种载荷情况,各载荷情况偏转坐标系输出到不同的偏转坐标系文件,文件名称采用载荷情况号。
CN201611161877.4A 2016-12-15 2016-12-15 一种基于有限元模型的飞机活动面偏转方法 Active CN106709174B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611161877.4A CN106709174B (zh) 2016-12-15 2016-12-15 一种基于有限元模型的飞机活动面偏转方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611161877.4A CN106709174B (zh) 2016-12-15 2016-12-15 一种基于有限元模型的飞机活动面偏转方法

Publications (2)

Publication Number Publication Date
CN106709174A true CN106709174A (zh) 2017-05-24
CN106709174B CN106709174B (zh) 2019-12-24

Family

ID=58937870

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611161877.4A Active CN106709174B (zh) 2016-12-15 2016-12-15 一种基于有限元模型的飞机活动面偏转方法

Country Status (1)

Country Link
CN (1) CN106709174B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451337A (zh) * 2017-07-07 2017-12-08 中国航空工业集团公司西安飞机设计研究所 一种襟翼偏转坐标系建立方法
CN107526876A (zh) * 2017-08-01 2017-12-29 中国航空工业集团公司西安飞机设计研究所 一种三缝富勒襟翼多姿态快速建模方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915563A (zh) * 2010-07-20 2010-12-15 中国航空工业集团公司西安飞机设计研究所 一种飞机舵偏转角测量方法
CN101963499A (zh) * 2010-07-21 2011-02-02 中国航空工业集团公司西安飞机设计研究所 一种飞机舵面偏转角测量工具及测量方法
WO2016079987A1 (ja) * 2014-11-19 2016-05-26 パナソニックIpマネジメント株式会社 入出力操作装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915563A (zh) * 2010-07-20 2010-12-15 中国航空工业集团公司西安飞机设计研究所 一种飞机舵偏转角测量方法
CN101963499A (zh) * 2010-07-21 2011-02-02 中国航空工业集团公司西安飞机设计研究所 一种飞机舵面偏转角测量工具及测量方法
WO2016079987A1 (ja) * 2014-11-19 2016-05-26 パナソニックIpマネジメント株式会社 入出力操作装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451337A (zh) * 2017-07-07 2017-12-08 中国航空工业集团公司西安飞机设计研究所 一种襟翼偏转坐标系建立方法
CN107526876A (zh) * 2017-08-01 2017-12-29 中国航空工业集团公司西安飞机设计研究所 一种三缝富勒襟翼多姿态快速建模方法
CN107526876B (zh) * 2017-08-01 2020-08-11 中国航空工业集团公司西安飞机设计研究所 一种三缝富勒襟翼多姿态快速建模方法

Also Published As

Publication number Publication date
CN106709174B (zh) 2019-12-24

Similar Documents

Publication Publication Date Title
US9671777B1 (en) Training robots to execute actions in physics-based virtual environment
CN102073270B (zh) 单输入单输出时滞系统的分数阶pid控制方法
US10198422B2 (en) Information-processing equipment based on a spreadsheet
KR102335814B1 (ko) 스마트팩토리용 데이터 모델 생성 장치 및 그 방법
KR101862221B1 (ko) 항공기 비행제어법칙 시뮬레이션 방법 및 장치
CN106777689B (zh) 一种基于有限元模型的飞机双铰链舵面偏转方法
CN108595888A (zh) 面向工业机器人设计与验证的仿真平台与方法
CN107341847A (zh) 一种基于bim技术的钢结构拼装构件三维模型数据处理方法
CN104875204A (zh) 一种等离子空间切割机器人的离线编程模块及应用方法
WO2021205776A1 (ja) シミュレーションのためのプログラムおよび装置
CN106709174A (zh) 一种基于有限元模型的飞机活动面偏转方法
CN106914896A (zh) 一种机器人离线编程的构建方法
Du Val et al. FLIGHTLAB TM modeling for real-time simulation applications
CN103942380A (zh) 图形化控制系统设计与仿真工具
CN109986109B (zh) 钢构建材的钻孔加工设计方法
CN106227944B (zh) 环形零件的三维模型径向力的施加方法
Vosniakos et al. Structured design of flexibly automated manufacturing cells through semantic models and petri nets in a virtual reality environment
CN110334459B (zh) 一种输电塔塔线体系快速精细化建模系统及方法
Laguionie et al. A multi-process manufacturing approach based on STEP-NC data model
Stecken et al. Creating and using digital twins within simulation environments
CN112131707A (zh) 一种基于mockup的装配公差对臂架末端精度影响的仿真分析方法
CN207172078U (zh) 一种多异构工业机器人控制装置
CN105717871B (zh) 可快速配置通用多轴运动控制系统的方法
CN105739958B (zh) 一种整体框架有限元模型快速建立方法
Laguionie et al. Simulation and optimization in a multi-process environment using STEP-NC

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant