CN111488684A - 一种载荷平衡计算方法 - Google Patents
一种载荷平衡计算方法 Download PDFInfo
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012360 testing method Methods 0.000 abstract description 11
- 238000013461 design Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force 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
技术领域
本申请属于飞机结构强度试验中载荷平衡领域,特别涉及一种载荷平衡计算方法。
背景技术
飞机结构强度试验是通过对飞机结构体施加模拟载荷,测量应力、应变及位移等力学参数,从而对结构体的承载能力和结构寿命作出正确的评价和估计,并为验证和优化结构设计提供可靠的依据。在静强度验证试验中,载荷处理是试验设计过程中的必要步骤。载荷处理的过程即是通过对已知的任务载荷的等效计算,求解出能够满足工程应用的可施加的载荷。目前现有技术中的计算方法还存在试验设计的效率低和工程试验实施的难度较高等缺点。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种载荷平衡计算方法,以解决现有级数存在的至少一个问题。
本申请的技术方案是:
一种载荷平衡计算方法,包括:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建所述等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型;
步骤四:根据所述载荷处理问题模型计算载荷谱。
可选地,步骤一中,所述创建等效载荷方程包括:
其中,Fx、Fy、Fz分别为任务载荷合力在x、y、z三个方向的力分量,Mx、My、Mz分别为任务载荷合力矩在x、y、z三个方向的力矩分量,fi为所求施加载荷的大小,(ai,bi,ci)为的方向向量,(xi,yi,zi)为的作用点坐标。
可选地,步骤二中,所述创建所述等效载荷方程的约束方程包括:
构建等式约束的最小二乘数学模型如下:
其中,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]],存在关系:
根据所述条件极值构建载荷处理问题模型,如下:
可选地,步骤三中,还包括获取所述载荷处理问题模型的迭代式:
对所述载荷处理问题模型进行展开可得到载荷处理问题的迭代式为:
发明至少存在以下有益技术效果:
本申请的载荷平衡计算方法,提高加载点设计效率,降低设计员的劳动强度,优化试验实施布局。
附图说明
图1是本申请一个实施方式的载荷平衡计算方法流程图;
图2是本申请一个实施方式的载荷平衡计算方法对载荷处理问题模型进行迭代计算的流程图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1至图2对本申请做进一步详细说明。
本申请提供了一种载荷平衡计算方法,包括以下步骤:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解约束方程的条件极值,并根据条件极值构建载荷处理问题模型;
步骤四:根据载荷处理问题模型计算载荷谱。
本申请的载荷平衡计算方法,载荷处理要求满足力系等效原则,并且必须在载荷限定范围之内,过大载荷会导致设备无法施加或超过试验件局部结构的承受能力。
在本申请的一个实施方式中,根据上述原则,首先创建等效载荷方程包括:
其中,Fx、Fy、Fz分别为任务载荷合力在x、y、z三个方向的力分量,Mx、My、Mz分别为任务载荷合力矩在x、y、z三个方向的力矩分量,fi为所求施加载荷的大小,(ai,bi,ci)为的方向向量,(xi,yi,zi)为的作用点坐标。
其次,在创建等效载荷方程之后,创建等效载荷方程的约束方程包括:
构建等式约束的最小二乘数学模型如下:
其中,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]],存在关系:
根据条件极值构建载荷处理问题模型,如下:
在步骤三中,构建完成载荷处理问题模型后,还包括获取载荷处理问题模型的迭代式,便于通过软件实现载荷谱的计算:
对载荷处理问题模型进行展开可得到载荷处理问题的迭代式为:
在本申请的一个实施方式中,步骤四中,根据上述的载荷处理问题的迭代式计算载荷谱,程序设计流程如下:
(a)首先创建初始变量,一般设置为0;
(b)根据公式,建立迭代循环函数,并设置循环次数最大限制;
(c)检查最终变量值是否满足范围要求,如果满足则输出结果,不满足则再次进入迭代循环函数。
本申请的载荷平衡计算方法,将任务载荷根据力系等效原则分配到各个相应的载荷施加设备上,便于工程实施。本申请在计算载荷谱的过程中,简化了程序设计难度,提高了计算效率,实现自动化;提高加载点设计效率,降低设计员的劳动强度,优化试验实施布局。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (5)
1.一种载荷平衡计算方法,其特征在于,包括:
步骤一:创建等效载荷方程;
步骤二:基于最小二乘数学模型创建所述等效载荷方程的约束方程;
步骤三:根据拉格朗日乘子法求解所述约束方程的条件极值,并根据所述条件极值构建载荷处理问题模型;
步骤四:根据所述载荷处理问题模型计算载荷谱。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113536469A (zh) * | 2021-08-03 | 2021-10-22 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种驱动装置的输出力矩计算方法 |
Citations (6)
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 | 中国飞机强度研究所 | 基于最小二乘法的载荷优化方法 |
-
2020
- 2020-04-12 CN CN202010282730.0A patent/CN111488684A/zh active Pending
Patent Citations (6)
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)
Title |
---|
孙文瑜等: "《最优化方法》", 北京:高等教育出版社, pages: 148 - 151 * |
高建等: "《基于遗传算法的飞机静强度试验载荷处理算法研究与应用》", 《工程与试验》, pages 20 - 21 * |
Cited By (2)
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 |