CN111488684A - 一种载荷平衡计算方法 - Google Patents

一种载荷平衡计算方法 Download PDF

Info

Publication number
CN111488684A
CN111488684A CN202010282730.0A CN202010282730A CN111488684A CN 111488684 A CN111488684 A CN 111488684A CN 202010282730 A CN202010282730 A CN 202010282730A CN 111488684 A CN111488684 A CN 111488684A
Authority
CN
China
Prior art keywords
load
equation
calculation method
balance calculation
creating
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.)
Pending
Application number
CN202010282730.0A
Other languages
English (en)
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.)
AVIC Aircraft Strength Research Institute
Original Assignee
AVIC Aircraft Strength Research Institute
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 AVIC Aircraft Strength Research Institute filed Critical AVIC Aircraft Strength Research Institute
Priority to CN202010282730.0A priority Critical patent/CN111488684A/zh
Publication of CN111488684A publication Critical patent/CN111488684A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Operations Research (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本申请属于飞机结构强度试验中载荷平衡领域,特别涉及一种载荷平衡计算方法。包括:步骤一:创建等效载荷方程;步骤二:基于最小二乘数学模型创建所述等效载荷方程的约束方程;步骤三:根据拉格朗日乘子法求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型;步骤四:根据所述载荷处理问题模型计算载荷谱。本申请的载荷平衡计算方法,将任务载荷根据力系等效原则分配到各个相应的载荷施加设备上,便于工程实施。本申请在计算载荷谱的过程中,简化了程序设计难度,提高了计算效率,实现自动化;提高加载点设计效率,降低设计员的劳动强度,优化试验实施布局。

Description

一种载荷平衡计算方法
技术领域
本申请属于飞机结构强度试验中载荷平衡领域,特别涉及一种载荷平衡计算方法。
背景技术
飞机结构强度试验是通过对飞机结构体施加模拟载荷,测量应力、应变及位移等力学参数,从而对结构体的承载能力和结构寿命作出正确的评价和估计,并为验证和优化结构设计提供可靠的依据。在静强度验证试验中,载荷处理是试验设计过程中的必要步骤。载荷处理的过程即是通过对已知的任务载荷的等效计算,求解出能够满足工程应用的可施加的载荷。目前现有技术中的计算方法还存在试验设计的效率低和工程试验实施的难度较高等缺点。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种载荷平衡计算方法,以解决现有级数存在的至少一个问题。
本申请的技术方案是:
一种载荷平衡计算方法,包括:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建所述等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型;
步骤四:根据所述载荷处理问题模型计算载荷谱。
可选地,步骤一中,所述创建等效载荷方程包括:
Figure BDA0002447331900000021
其中,Fx、Fy、Fz分别为任务载荷合力在x、y、z三个方向的力分量,Mx、My、Mz分别为任务载荷合力矩在x、y、z三个方向的力矩分量,fi为所求施加载荷
Figure BDA0002447331900000022
的大小,(ai,bi,ci)为
Figure BDA0002447331900000023
的方向向量,(xi,yi,zi)为
Figure BDA0002447331900000024
的作用点坐标。
可选地,步骤二中,所述创建所述等效载荷方程的约束方程包括:
构建等式约束的最小二乘数学模型如下:
Figure BDA0002447331900000025
其中,K为载荷等效约束矩阵,k1i=ai,k2i=bi,k3i=ci,k4i=(yici-zibi),k5i=(ziai-xici),k6i=(xibi-yiai),因此矩阵K为6×n的矩阵;G=[Fx Fy Fz Mx My Mz]T表示所有任务载荷在x、y、z三个方向的合力及合力矩。
可选地,步骤三中,所述求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型包括:
利用拉格朗日乘子法求解的条件极值问题,得到:
X=Xo-KT((KKT)-1(KXo-G))
定义载荷限制运算“[[]]”,设X=[x1,…,xi,…,xk]T,Y=[y1,…,yi,…,yk]T,则Y=[[X]],存在关系:
Figure BDA0002447331900000026
根据所述条件极值构建载荷处理问题模型,如下:
Figure BDA0002447331900000027
可选地,步骤三中,还包括获取所述载荷处理问题模型的迭代式:
对所述载荷处理问题模型进行展开可得到载荷处理问题的迭代式为:
Figure BDA0002447331900000031
发明至少存在以下有益技术效果:
本申请的载荷平衡计算方法,提高加载点设计效率,降低设计员的劳动强度,优化试验实施布局。
附图说明
图1是本申请一个实施方式的载荷平衡计算方法流程图;
图2是本申请一个实施方式的载荷平衡计算方法对载荷处理问题模型进行迭代计算的流程图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1至图2对本申请做进一步详细说明。
本申请提供了一种载荷平衡计算方法,包括以下步骤:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解约束方程的条件极值,并根据条件极值构建载荷处理问题模型;
步骤四:根据载荷处理问题模型计算载荷谱。
本申请的载荷平衡计算方法,载荷处理要求满足力系等效原则,并且必须在载荷限定范围之内,过大载荷会导致设备无法施加或超过试验件局部结构的承受能力。
在本申请的一个实施方式中,根据上述原则,首先创建等效载荷方程包括:
Figure BDA0002447331900000041
其中,Fx、Fy、Fz分别为任务载荷合力在x、y、z三个方向的力分量,Mx、My、Mz分别为任务载荷合力矩在x、y、z三个方向的力矩分量,fi为所求施加载荷
Figure BDA0002447331900000042
的大小,(ai,bi,ci)为
Figure BDA0002447331900000043
的方向向量,(xi,yi,zi)为
Figure BDA0002447331900000044
的作用点坐标。
其次,在创建等效载荷方程之后,创建等效载荷方程的约束方程包括:
构建等式约束的最小二乘数学模型如下:
Figure BDA0002447331900000045
其中,K为载荷等效约束矩阵,k1i=ai,k2i=bi,k3i=ci,k4i=(yici-zibi),k5i=(ziai-xici),k6i=(xibi-yiai),因此矩阵K为6×n的矩阵;G=[Fx Fy Fz Mx My Mz]T表示所有任务载荷在x、y、z三个方向的合力及合力矩。
本实施例中,步骤三中,求解约束方程的条件极值,并根据条件极值构建载荷处理问题模型包括:
利用拉格朗日乘子法求解的条件极值问题,得到:
X=Xo-KT((KKT)-1(KXo-G))
定义载荷限制运算“[[]]”,设X=[x1,…,xi,…,xk]T,Y=[y1,…,yi,…,yk]T,则Y=[[X]],存在关系:
Figure BDA0002447331900000051
根据条件极值构建载荷处理问题模型,如下:
Figure BDA0002447331900000052
在步骤三中,构建完成载荷处理问题模型后,还包括获取载荷处理问题模型的迭代式,便于通过软件实现载荷谱的计算:
对载荷处理问题模型进行展开可得到载荷处理问题的迭代式为:
Figure BDA0002447331900000053
在本申请的一个实施方式中,步骤四中,根据上述的载荷处理问题的迭代式计算载荷谱,程序设计流程如下:
(a)首先创建初始变量,一般设置为0;
(b)根据公式,建立迭代循环函数,并设置循环次数最大限制;
(c)检查最终变量值是否满足范围要求,如果满足则输出结果,不满足则再次进入迭代循环函数。
本申请的载荷平衡计算方法,将任务载荷根据力系等效原则分配到各个相应的载荷施加设备上,便于工程实施。本申请在计算载荷谱的过程中,简化了程序设计难度,提高了计算效率,实现自动化;提高加载点设计效率,降低设计员的劳动强度,优化试验实施布局。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。

Claims (5)

1.一种载荷平衡计算方法,其特征在于,包括:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建所述等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型;
步骤四:根据所述载荷处理问题模型计算载荷谱。
2.根据权利要求1所述的载荷平衡计算方法,其特征在于,步骤一中,所述创建等效载荷方程包括:
Figure FDA0002447331890000011
其中,Fx、Fy、Fz分别为任务载荷合力在x、y、z三个方向的力分量,Mx、My、Mz分别为任务载荷合力矩在x、y、z三个方向的力矩分量,fi为所求施加载荷
Figure FDA0002447331890000012
的大小,(ai,bi,ci)为
Figure FDA0002447331890000013
的方向向量,(xi,yi,zi)为
Figure FDA0002447331890000014
的作用点坐标。
3.根据权利要求2所述的载荷平衡计算方法,其特征在于,步骤二中,所述创建所述等效载荷方程的约束方程包括:
构建等式约束的最小二乘数学模型如下:
Figure FDA0002447331890000015
其中,K为载荷等效约束矩阵,k1i=ai,k2i=bi,k3i=ci,k4i=(yici-zibi),k5i=(ziai-xici),k6i=(xibi-yiai),因此矩阵K为6×n的矩阵;G=[Fx Fy Fz Mx My Mz]T表示所有任务载荷在x、y、z三个方向的合力及合力矩。
4.根据权利要求3所述的载荷平衡计算方法,其特征在于,步骤三中,所述求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型包括:
利用拉格朗日乘子法求解的条件极值问题,得到:
X=Xo-KT((KKT)-1(KXo-G))
定义载荷限制运算“[[]]”,设X=[x1,…,xi,…,xk]T,Y=[y1,…,yi,…,yk]T,则Y=[[X]],存在关系:
Figure FDA0002447331890000021
根据所述条件极值构建载荷处理问题模型,如下:
Figure FDA0002447331890000022
5.根据权利要求4所述的载荷平衡计算方法,其特征在于,步骤三中,还包括获取所述载荷处理问题模型的迭代式:
对所述载荷处理问题模型进行展开可得到载荷处理问题的迭代式为:
Figure FDA0002447331890000023
CN202010282730.0A 2020-04-12 2020-04-12 一种载荷平衡计算方法 Pending CN111488684A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010282730.0A CN111488684A (zh) 2020-04-12 2020-04-12 一种载荷平衡计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010282730.0A CN111488684A (zh) 2020-04-12 2020-04-12 一种载荷平衡计算方法

Publications (1)

Publication Number Publication Date
CN111488684A true CN111488684A (zh) 2020-08-04

Family

ID=71811792

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010282730.0A Pending CN111488684A (zh) 2020-04-12 2020-04-12 一种载荷平衡计算方法

Country Status (1)

Country Link
CN (1) CN111488684A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113536469A (zh) * 2021-08-03 2021-10-22 中国航空工业集团公司沈阳飞机设计研究所 一种驱动装置的输出力矩计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995918A (zh) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 一种机翼变形和振动对飞机传递对准影响的分析方法
CN104217072A (zh) * 2014-08-26 2014-12-17 中国直升机设计研究所 一种起落架机轮非线性动力学建模方法
CA2895895A1 (en) * 2014-07-15 2016-01-15 Dassault Systemes Simulia Corp. System and method of recovering lagrange multipliers in modal dynamic analysis
CN109583044A (zh) * 2018-11-09 2019-04-05 中国直升机设计研究所 基于rbf神经网络的直升机旋翼飞行载荷预测方法
CN109733638A (zh) * 2018-09-07 2019-05-10 中国飞机强度研究所 一种长支柱起落架大变形情况下载荷施加方法
CN109977556A (zh) * 2019-03-28 2019-07-05 中国飞机强度研究所 基于最小二乘法的载荷优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995918A (zh) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 一种机翼变形和振动对飞机传递对准影响的分析方法
CA2895895A1 (en) * 2014-07-15 2016-01-15 Dassault Systemes Simulia Corp. System and method of recovering lagrange multipliers in modal dynamic analysis
CN104217072A (zh) * 2014-08-26 2014-12-17 中国直升机设计研究所 一种起落架机轮非线性动力学建模方法
CN109733638A (zh) * 2018-09-07 2019-05-10 中国飞机强度研究所 一种长支柱起落架大变形情况下载荷施加方法
CN109583044A (zh) * 2018-11-09 2019-04-05 中国直升机设计研究所 基于rbf神经网络的直升机旋翼飞行载荷预测方法
CN109977556A (zh) * 2019-03-28 2019-07-05 中国飞机强度研究所 基于最小二乘法的载荷优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孙文瑜等: "《最优化方法》", 北京:高等教育出版社, pages: 148 - 151 *
高建等: "《基于遗传算法的飞机静强度试验载荷处理算法研究与应用》", 《工程与试验》, pages 20 - 21 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113536469A (zh) * 2021-08-03 2021-10-22 中国航空工业集团公司沈阳飞机设计研究所 一种驱动装置的输出力矩计算方法
CN113536469B (zh) * 2021-08-03 2023-10-10 中国航空工业集团公司沈阳飞机设计研究所 一种驱动装置的输出力矩计算方法

Similar Documents

Publication Publication Date Title
CN104915483B (zh) 一种地震作用下建筑物基础稳定性验算方法
CN108121865B (zh) 一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法
CN111159943B (zh) 一种动翼面封严结构的屈曲处理方法
CN114707390B (zh) 基于分层壳材料本构参数更新混合试验方法及分析方法
CN108763810A (zh) 一种桥梁静载试验荷载加载布置与调整方法
CN111488684A (zh) 一种载荷平衡计算方法
US20150073730A1 (en) Mechanical strain gauge simulation
CN107655650B (zh) 用于结构试验的杠杆加载系统碰撞检测方法
CN112818462A (zh) 车轮参数模型的生成方法、装置、存储介质和计算机设备
CN104834795B (zh) 包带连接结构接触摩擦非线性特性模拟方法及系统
CN107609234B (zh) 基于不确定性的稳健性分析方法及系统
CN113656891A (zh) 液体火箭动特性建模分析方法及终端设备
CN113158510A (zh) 一种索网结构的优化找形方法
CN115374667B (zh) 基于钢结构施工变形对结构设计性能影响的逆分析方法
CN109703781B (zh) 一种杠杆系统的计算设计方法
CN117094071A (zh) 一种卫星sar天线动力学快速建模方法及系统
Sayed et al. Inertia properties of indefinite quadratic forms
CN112528404A (zh) 适用于船舶舾装平台的格栅设计方法、装置、介质及终端
CN112818502A (zh) 光学镜面面形计算方法
KR20160077902A (ko) 시뮬레이션 인터페이스 생성방법 및 실시간 시뮬레이션 장치
CN106777650B (zh) 一种集中传载金属接头有限元网格尺寸确定方法
CN112528528A (zh) 一种混合结构的抗震计算模态叠加方法
CN115546210B (zh) 一种数字图像相关中的局部应变加权平均方法
CN110108544A (zh) 放射状杠杆加载系统空间初始位置确定方法
Ren et al. Evaluation of APMP. MM-K1 inter-comparisons by Monte Carlol method

Legal Events

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