[go: up one dir, main page]

WO2023240860A1 - 一种自动铺丝工装的设计方法、装置、设备及存储介质 - Google Patents

一种自动铺丝工装的设计方法、装置、设备及存储介质 Download PDF

Info

Publication number
WO2023240860A1
WO2023240860A1 PCT/CN2022/125592 CN2022125592W WO2023240860A1 WO 2023240860 A1 WO2023240860 A1 WO 2023240860A1 CN 2022125592 W CN2022125592 W CN 2022125592W WO 2023240860 A1 WO2023240860 A1 WO 2023240860A1
Authority
WO
WIPO (PCT)
Prior art keywords
section
tooling
beam structure
deformation
cross
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.)
Ceased
Application number
PCT/CN2022/125592
Other languages
English (en)
French (fr)
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.)
Chengdu Aircraft Industrial Group Co Ltd
Original Assignee
Chengdu Aircraft Industrial Group Co Ltd
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 Chengdu Aircraft Industrial Group Co Ltd filed Critical Chengdu Aircraft Industrial Group Co Ltd
Publication of WO2023240860A1 publication Critical patent/WO2023240860A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C70/00—Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts
    • B29C70/04—Shaping composites, i.e. plastics material comprising reinforcements, fillers or preformed parts, e.g. inserts comprising reinforcements only, e.g. self-reinforcing plastics
    • B29C70/28—Shaping operations therefor
    • B29C70/30—Shaping by lay-up, i.e. applying fibres, tape or broadsheet on a mould, former or core; Shaping by spray-up, i.e. spraying of fibres on a mould, former or core
    • B29C70/38—Automated lay-up, e.g. using robots, laying filaments according to predetermined patterns
    • B29C70/382—Automated fiber placement [AFP]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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

Definitions

  • This application relates to the technical field of aviation equipment, and specifically to a design method, device, equipment and storage medium for automatic wire laying equipment.
  • the cylindrical parts in some aircraft structures can be manufactured by automatic filament laying.
  • the tooling must rotate at a constant speed or at an accelerated speed and cannot block the filament laying head. Therefore, the main features of the automatic filament laying tool for cylindrical parts are
  • the structure of the tooling is cylindrical. Both ends of the tooling are connected to the rotary drive joints of the wire laying machine. The middle part of the tooling is suspended in the air without any support.
  • the automatic laying equipment for cylindrical parts will bear four types of loads: gravity, the pressure of the head of the laying machine, the support force at both ends, and the centrifugal force generated by the accelerated rotation of the laying machine. Therefore, when designing automatic wire laying tooling for cylindrical parts, it is necessary to accurately calculate the maximum deformation of the tooling during the automatic wire laying process, as a basis for measuring whether the tooling meets the technical requirements and optimizing the tooling structure.
  • the main purpose of this application is to provide a design method, device, equipment and storage medium, device, equipment and storage medium for automatic wire laying tooling, aiming to solve the shortcoming of the existing technology that the design cycle of the silk laying tooling is too long.
  • this application provides a design method for automatic wire laying equipment, including the following steps:
  • the deformation curve function of the equivalent beam structure is obtained according to the load function, and the deformation amount of the equivalent beam structure is calculated according to the deformation curve function;
  • the wire laying tooling is qualified, otherwise the load distribution of the equivalent beam structure is determined repeatedly. and cross-section size parameters until the deformation amount is qualified.
  • the wire laying tooling includes a main shaft and a number of auxiliary structures installed on the main shaft, then the wire laying tooling is determined to be a main shaft and its auxiliary structure type tooling, and the main shaft is converted into an equivalent beam structure, and each of the auxiliary structures is The structure is applied as an external load on the equivalent beam structure;
  • the wire-laying tooling is an integrated structure, it is determined that the wire-laying tooling is an integral structural tooling, and the wire-laying tooling is converted into an equivalent beam structure.
  • determining the load distribution and cross-sectional size parameters of the equivalent beam structure includes the following steps:
  • the structural form is a main axis and its auxiliary structural tooling
  • the main axis is a beam of equal cross-section
  • the mass of each of the auxiliary structures is evenly arranged along the axial direction of the wire laying tooling
  • the equivalent beam structure is For a cross-section beam, the cross-sectional size parameters of the equal-section beam are determined as the cross-sectional size parameters of the main axis;
  • the equivalent beam structure is For equal-section beams that bear non-uniformly distributed loads, the equivalent beam structure is divided into several segmented equal-section beams, and the cross-sectional size of the main axis in the section of each segmented equal-section beam is determined as each of the segmented equal-section beams.
  • the cross-sectional size parameters, the accessory mechanisms on each segmented cross-section beam are applied as external loads on the equivalent beam structure;
  • the equivalent beam structure is a variable-section beam, and the equivalent beam structure is divided into Several segmented equal-section beams, and the cross-sectional size of the main axis at the midpoint section of each segmented equal-section beam is determined as the cross-sectional size parameter of each of the segmented equal-section beams, and the accessory mechanism on each of the segmented equal-section beams is as External loads are applied to the equivalent beam structure;
  • the structural form is an overall structure type tooling, and the overall structure of the wire laying tooling is a beam of constant cross-section, then its equivalent beam structure is a beam of constant cross-section, and the cross-sectional size of the equivalent beam structure is determined as the cross-sectional size of the wire-laying tooling; Otherwise, if it is a variable cross-section beam, the equal cross-section is divided into several segmented equal-section beams, and the cross-sectional size corresponding to the midpoint section of each segmented cross-section beam is used as the cross-sectional size parameter of each segmented equal-section beam.
  • the length u of the segmented equal-section beam and the total length U of the main axis satisfy the following relationship: u ⁇ U(1-A)/10, where A represents the calculation accuracy, and its value range is 0-1.
  • calculating the load function of the equivalent beam according to the load distribution and size parameters includes the following steps:
  • the load function is the following expression:
  • Z 0 0
  • the P i is the concentrated force
  • Q j is the inertial force
  • Q k is the weight of the accessory structure
  • Z i is the coordinate of the concentrated force
  • Z j-1 and Z j are respectively the segmented beams of equal cross-section corresponding to the inertial force
  • Z k-1 and Z k are respectively the starting and ending coordinates of each segment of the segmented equal-section beam corresponding to the gravity Q k of the accessory structure
  • i, j, k, r, m, n are all Positive integers, all less than or equal to the number of segmented equal-section beams, i, j, k represent the number of each segment, r, m, n represent the number of segmented segments;
  • a represents the coordinates of the segmented plane;
  • z represents the point to be solved coordinate.
  • the load function is the following expression:
  • Z 0 0
  • P i is the concentrated force
  • Q j is the inertial force
  • Z i is the coordinate of the concentrated force
  • Z j-1 and Z j are respectively the starting and ending coordinates of each section of the segmented equal-section beam corresponding to the inertial force Q j
  • i, j, r, and m are all positive integers, and are less than or equal to the number of segmented equal-section beams.
  • i and j represent the numbers of each segment
  • r and m represent the number of segmented segments
  • a represents the coordinates of the segmented plane
  • z represents the coordinates of the segment to be Find the coordinates of the solution point.
  • calculating the maximum deformation of the equivalent beam structure according to the deformation curve function includes the following steps:
  • the hinge boundary conditions are brought into the deformation curve function to calculate the hinge support deformation, and the fixed support boundary conditions are brought into the deformation curve function to calculate the fixed support deformation;
  • hinge deformation and the fixed support deformation Compare the hinge deformation and the fixed support deformation, and use the larger of the hinge deformation and the fixed support deformation as the upper limit of the maximum deformation.
  • the hinge deformation and the fixed support deformation are The smaller one is used as the lower limit of the maximum deformation.
  • this application also provides a design device for automatic wire laying equipment, including:
  • the equivalent conversion module is used to determine the structural form of the wire-laying tooling, and convert the wire-laying tooling into an equivalent beam structure according to the structural form;
  • a data generation module used to determine the load distribution and cross-sectional size parameters of the equivalent beam structure
  • a first calculation module configured to calculate the load function of the equivalent beam structure according to the load distribution and size parameters
  • the second calculation module is used to derive the deformation curve function of the equivalent beam structure based on the load function, and calculate the deformation amount of the equivalent beam structure based on the deformation curve function;
  • the judgment module combines the deformation amount and whether it meets the deformation amount judgment standard of the wire laying tooling. If the deformation meets the deformation amount judgment standard, the wire laying tooling is qualified. Otherwise, the equivalent beam structure is determined repeatedly. load distribution and cross-sectional size parameters until the deformation amount is qualified.
  • the present application also provides an electronic device.
  • the electronic device includes a memory and a processor.
  • a computer program is stored in the memory, and the processor executes the computer program to implement the above method.
  • the present application also provides a computer-readable storage medium, a computer program is stored on the computer-readable storage medium, and the processor executes the computer program to implement the above method.
  • This application first determines the structural form of the wire laying tool, then converts the wire laying tool into an equivalent beam structure based on the structural form, and determines the load distribution and cross-sectional size of the equivalent beam structure based on the size of the wire laying tool. Finally, Calculate the deformation range of the wire-laying tooling by combining the load distribution of the equivalent beam structure and the deformation curve function; at the same time, determine whether the calculated deformation amount of the wire-laying tooling is qualified based on the deformation criterion;
  • this application converts the complex wire-laying tooling into an equivalent beam structure with a relatively simple structure, and then reflects the deformation of the wire-laying tooling through stress analysis and deformation analysis of the equivalent beam structure; It realizes the simplification of the calculation process of the deformation of the wire laying tooling while ensuring the calculation accuracy. Therefore, in the technical solution described in this application, only the key parameters such as the structural form, load distribution and cross-sectional size of the wire laying tooling need to be determined.
  • this application realizes the calculation of the deformation of the wire laying tooling before the detailed design of the tooling after the technical solution and key dimensions of the wirelaying tooling are determined, which not only reduces the time required for the wirelaying tooling
  • the detailed design cycle of the tooling and at the same time, it can quickly determine whether the structural form and key dimensions of the wire-laying tooling are reasonable during the detailed design process, and quickly determine whether the maximum deformation of the wire-laying tooling meets the technical requirements, and provide a basis for the subsequent detailed design process of the wire-laying tooling.
  • the iterative optimization of various parameters provides data support and basis, thereby further shortening the design cycle of the entire wire laying tooling.
  • the equivalent beam structure modeling and load calculation method used in the design method of automatic wire laying tooling described in this application includes the manufacturing and assembly errors of the tooling, and provides the corresponding tooling deformation when the error is minimum and maximum. value, which can effectively guide the actual engineering and accurately evaluate the extreme conditions of the tooling.
  • Figure 1 is a schematic structural diagram of an electronic device of the hardware operating environment involved in the present invention
  • Figure 2 is a flow chart of a design method for automatic wire laying equipment provided by an embodiment of the present application
  • Figure 3 is a schematic structural diagram of the silk laying tooling described in Embodiment 1 of the present application.
  • Figure 4 is a cross-sectional view of the intermediate mandrel in Embodiment 1 of the present application.
  • Figure 5 is a cross-sectional view of the connection joint described in Embodiment 1 of the present application.
  • Figure 6 is a stress analysis diagram of the state of the hinges at both ends of the wire laying tool described in Embodiment 1 of the present application;
  • Figure 7 is a stress analysis diagram of the fixed support state at both ends of the wire laying tooling described in Embodiment 1 of the present application;
  • Figure 8 is a functional module schematic diagram of the design device of the automatic wire laying tooling according to the present invention.
  • Figure 9 is a schematic structural diagram of the wire laying tooling described in Embodiment 2 of the present application.
  • Figure 10 is a stress analysis diagram of the hinge state at both ends of the wire laying tool described in Embodiment 2 of the present application;
  • Figure 11 is a stress analysis diagram of the fixed support state at both ends of the wire laying tooling described in Embodiment 2 of the present application;
  • connection can be a fixed connection, a detachable connection, or an integral body; It can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; it can be an internal connection between two elements or an interactive relationship between two elements, unless otherwise clearly limited.
  • fixing can be a fixed connection, a detachable connection, or an integral body; It can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; it can be an internal connection between two elements or an interactive relationship between two elements, unless otherwise clearly limited.
  • Figure 1 is a schematic diagram of the electronic device structure of the hardware operating environment involved in the embodiment of the present application.
  • the electronic device may include: a processor 1001, such as a central processing unit (Central Processing Unit, CPU), a communication bus 1002, a user interface 1003, a network interface 1004, and a memory 1005.
  • the communication bus 1002 is used to realize connection communication between these components.
  • the user interface 1003 may include a display screen (Display) and an input unit such as a keyboard (Keyboard).
  • the optional user interface 1003 may also include a standard wired interface and a wireless interface.
  • the network interface 1004 may optionally include a standard wired interface or a wireless interface (such as a wireless fidelity (WIreless-FIdelity, WI-FI) interface).
  • the memory 1005 can be a high-speed random access memory (Random Access Memory, RAM) memory or a stable non-volatile memory (Non-Volatile Memory, NVM), such as a disk memory.
  • RAM Random Access Memory
  • NVM Non-Volatile Memory
  • the memory 1005 may optionally be a storage device independent of the aforementioned processor 1001.
  • Figure 1 does not constitute a limitation of the electronic device, and may include more or fewer components than shown, or combine certain components, or arrange different components.
  • a memory as a storage medium may include an operating system, a data storage module, a network communication module, a user interface module and an electronic program.
  • the network interface 1004 is mainly used for data communication with the network server;
  • the user interface 1003 is mainly used for data interaction with the user;
  • the processor 1001 and the memory 1005 in the electronic device of the present invention can be provided in
  • the electronic device calls the automatic wire laying tool design device stored in the memory 1005 through the processor 1001, and executes the automatic wire laying tool design method provided by the embodiment of the present application.
  • this application discloses a design method for automatic wire laying equipment, which includes the following steps:
  • the wire laying tool includes a main shaft located in the center.
  • auxiliary structures are provided around the main shaft, and each auxiliary structure is connected to the main shaft; when in use, the wire laying tool
  • the spindle is connected to external power equipment; if it meets the above characteristics, it is determined that the wire laying tooling is the main spindle and its ancillary structure tooling;
  • the fiber laying tool described in this embodiment is a main spindle and ancillary structural tooling.
  • the two ends of the mandrel in the middle of the fiber laying tool are connected to the rotary drive joint of the fiber laying machine, and the other attachments are The structure is not connected to the rotary drive joint of the filament laying machine, so the component composed of the middle mandrel and the rotary drive joint of the filament laying machine is converted into an equivalent beam structure, while the accessory structure is equivalent to an external load;
  • the structural form is a main shaft and its auxiliary structure tooling
  • the main shaft is a beam of equal cross-section
  • the mass of each auxiliary structure is evenly arranged along the axial direction of the wire laying tooling
  • the equivalent beam structure It is a beam of constant cross-section
  • the cross-sectional size parameters of the constant-section beam are the cross-sectional size parameters of the main axis
  • the main axis located in the middle is a beam of equal cross-section.
  • the mass of each accessory structure is evenly distributed along the axial direction of the main axis. Therefore, this embodiment performs calculations in the form of a beam of equal cross-section. ;
  • the material of the automatic wire laying tooling for cylindrical parts is INVAR steel, and the material parameters are shown in Table 1.
  • the length of the central mandrel of the fiber laying tool is 3900mm, and the length of the rotary drive joint of the fiber laying machine is 500mm;
  • the equivalent beam structure of this method includes the self-weight of the central core shaft, the weight of each accessory core membrane, and the weight of the central core membrane.
  • the mass of the upper core mold, the right core mold, the lower core mold and the left core mold are all 1085.28kg.
  • the distance between the cross-section centroids of the upper core mold, the right core mold, the lower core mold and the left core mold is the rotation axis.
  • step S44 Establish the deformation determination standard of the silk laying tooling based on the design requirements of the silk laying tooling. If the deformation in step S44 all falls within the range of the deformation determination standard, then the silk laying tooling design is qualified, otherwise step S21 needs to be adjusted. load distribution and cross-sectional size parameters in the model, and repeat the above calculation steps until the deformation amount is qualified.
  • an aircraft component frame beam gap design device including:
  • the equivalent conversion module is used to determine the structural form of the wire-laying tooling, and convert the wire-laying tooling into an equivalent beam structure according to the structural form;
  • a data generation module used to determine the load distribution and cross-sectional size parameters of the equivalent beam structure
  • a first calculation module configured to calculate the load function of the equivalent beam structure according to the load distribution and size parameters
  • the second calculation module is used to derive the deformation curve function of the equivalent beam structure based on the load function, and calculate the deformation amount of the equivalent beam structure based on the deformation curve function;
  • the judgment module combines the deformation amount and whether it meets the deformation amount judgment standard of the wire laying tooling. If the deformation meets the deformation amount judgment standard, the wire laying tooling is qualified. Otherwise, the equivalent beam structure is determined repeatedly. load distribution and cross-sectional size parameters until the deformation amount is qualified.
  • the fiber laying tool does not have a main shaft and ancillary structures installed on the main shaft.
  • both ends of the fiber laying tool are connected to the rotary drive joints of the fiber laying machine; therefore, the fiber laying tool is determined to be an integral structure. work clothes;
  • the wire-laying tooling described in this embodiment is an integral structural tooling, which does not have ancillary structures. Therefore, the wire-laying tooling is converted into an equivalent beam structure as a whole;
  • the structural form is an overall structure type tooling, and the overall structure of the wire laying tooling is a beam of constant cross-section, then its equivalent beam structure is a beam of constant cross-section, and the cross-sectional size of the equivalent beam structure is the cross-sectional size of the wire-laying tooling; Otherwise, it is a variable cross-section beam, and the equal cross-section is divided into a number of mutually spliced segmented equal-section beams, and the cross-sectional size parameter of each segmented equal-section beam is the cross-sectional size corresponding to the midpoint section of each segmented cross-section beam;
  • each cross section on the main body The areas are not equal, so the equivalent beam structure is determined to be a variable cross-section beam;
  • the maximum pressure of the fiber laying machine indenter is 100N, and the position of the fiber laying machine indenter is located at the midpoint of the tooling; the required calculation accuracy is 60%.
  • the material of the automatic wire laying tooling for cylindrical parts is INVAR steel, and the material parameters are shown in Table 5.
  • the load borne by the converted equivalent beam structure includes the tool's own weight, the support force at both ends, and the concentrated force of the wire laying machine indenter;
  • M 1 ⁇ M 27 are the masses of each section of the 27 sections of equal section beams
  • the supporting force of the laying machine on the rotating drive joint of the laying machine is P 1 .
  • the supporting force of the laying machine on the rotating driving joint of the laying machine is P 2 ;
  • M 1 ⁇ M 27 are the weights of each section of beams of equal section, where a represents the coordinates of the dividing plane;
  • F(z) is the equivalent beam structure load function after expansion in step S3
  • I(z) is the cross-section about the y-axis Moment of inertia.
  • the cross-sectional moment of inertia of each segmented equal-section beam in this embodiment is shown in Table 8.
  • C and D are both integral constants.
  • the maximum deformation value ⁇ max of the equivalent beam structure for the main force transmission of the automatic wire laying tool for cylindrical parts is between the maximum deformation values of the equivalent beam structure with hinged supports at both ends and fixed supports at both ends.
  • step S43 Establish the deformation determination standard of the silk laying tooling based on the design requirements of the silk laying tooling. If the deformation in step S43 all falls within the range of the deformation determination standard, then the silk laying tooling design is qualified, otherwise step S21 needs to be adjusted. load distribution and cross-sectional size parameters in the model, and repeat the above calculation steps until the deformation amount is qualified.
  • the wire-laying tooling in this embodiment is an overall structural tooling and a variable-section beam. If the overall structural tooling is a constant-section beam, there is no need to divide the entire equivalent beam structure. It only needs to be divided into The final result can be obtained by calculating the load function and deformation deflection curve function of the equivalent beam structure as a whole, that is, the calculation is performed by removing the relevant expressions of the accessory structure based on the calculation method described in Embodiment 1.
  • the wire laying tooling is the main spindle and its ancillary structure tooling, and it is a beam of equal cross-section. If any structure in the structure is a variable-section structure, then the main shaft and its auxiliary structural tooling are variable-section beams.
  • the main shaft needs to be divided according to the method adopted in step S22 of Embodiment 2, and then the load function is calculated according to the relevant methods recorded in Embodiment 2. and the deformation deflection curve function to get the final result;
  • this application converts the complex wire-laying tooling into an equivalent beam structure with a relatively simple structure, and then reflects the deformation of the wire-laying tooling through stress analysis and deformation analysis of the equivalent beam structure; It realizes the simplification of the calculation process of the deformation of the wire laying tooling while ensuring the calculation accuracy. Therefore, in the technical solution described in this application, only the key parameters such as the structural form, load distribution and cross-sectional size of the wire laying tooling need to be determined.
  • this application realizes the calculation of the deformation of the wire laying tooling before the detailed design of the tooling after the technical solution and key dimensions of the wirelaying tooling are determined, which not only reduces the time required for the wirelaying tooling
  • the detailed design cycle of the tooling and at the same time, it can quickly determine whether the structural form and key dimensions of the wire-laying tooling are reasonable during the detailed design process, and quickly determine whether the maximum deformation of the wire-laying tooling meets the technical requirements, and provide a basis for the subsequent detailed design process of the wire-laying tooling.
  • the iterative optimization of various parameters provides data support and basis, thereby further shortening the design cycle of the entire wire laying tooling.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Robotics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Composite Materials (AREA)
  • Mechanical Engineering (AREA)
  • Rod-Shaped Construction Members (AREA)
  • Working Measures On Existing Buildindgs (AREA)
  • Bridges Or Land Bridges (AREA)

Abstract

本申请公开了一种自动铺丝工装的设计方法、装置、设备及存储介质，先确定铺丝工装的结构形式，根据结构形式将铺丝工装转换为等效梁结构，结合铺丝工装的尺寸确定等效梁结构的载荷分布情况和截面尺寸，最后结合等效梁结构的载荷分布和变形曲线函数计算铺丝工装的变形范围，结合判定基准判定铺丝工装的变形量是否合格；本申请将铺丝工装转换为结构相对简单的等效梁结构，进而通过对等效梁结构的受力分析和变形分析反应铺丝工装的变形；因此只需确定铺丝工装的结构形式、载荷分布和截面尺寸等参数即可进行计算，简化变形计算需要的参数；从而在铺丝工装技术方案和关键尺寸确定之后，详细设计之前实现变形量的计算，缩减铺丝工装的详细设计周期。

Description

一种自动铺丝工装的设计方法、装置、设备及存储介质 技术领域
本申请涉及航空设备技术领域，具体涉及一种自动铺丝工装的设计方法、装置、设备及存储介质。
背景技术
部分飞机结构中的筒形件可采用自动铺丝的方式进行制造，自动铺丝时工装要进行匀速或加速旋转运动，且不能对铺丝头造成遮挡，因此筒形件自动铺丝工装主要特征为：工装结构形式为筒形，工装两端与铺丝机旋转驱动接头连接，工装中间部位悬空无任何支撑。筒形件自动铺丝工装在自动铺丝工艺过程中会承受重力、铺丝机压头压力、两端支撑力、铺丝机加速旋转产生的离心力4类载荷。因此在进行筒形件自动铺丝工装设计时，需对工装在自动铺丝工艺过程中的最大变形量进行准确计算，作为衡量工装是否满足技术要求以及工装结构优化的依据。
目前，在获取飞机筒形件自动铺丝工装最大变形量时，可采用实物测量、有限元仿真计算等方法；但实物测量时必须将工装制造出来后才可进行，而有限元仿真计算必须在完成工装结构图纸或数模详细设计后才可进行；由此导致在铺丝工装设计前期缺乏有效的数据支撑，导致设计周期过长。
发明内容
本申请的主要目的在于提供一种自动铺丝工装的设计方法、装置、设备及存储介质、装置、设备及存储介质，旨在解决现有技术中铺丝工装设计周期过长的缺陷。
为实现上述目的，本申请提供一种自动铺丝工装的设计方法，包 括以下步骤：
确定铺丝工装的结构形式，并根据所述结构形式将所述铺丝工装转换为等效梁结构；
确定所述等效梁结构的载荷分布及截面尺寸参数；
根据所述载荷分布和尺寸参数计算所述等效梁结构的载荷函数；
根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
结合所述变形量和是否满足铺丝工装的变形量判定基准，若所述变形量满足所述变形量判定基准，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤，直至所述变形量合格。
可选的，确定铺丝工装的结构形式，并根据所述结构形式将所述铺丝工装转换为等效梁结构，包括以下步骤：
若铺丝工装包括主轴及若干安装于所述主轴上的附属结构，则判定所述铺丝工装为主轴及其附属结构类工装，并将所述主轴转换为等效梁结构，各所述附属结构作为外载荷施加于所述等效梁结构上；
若铺丝工装为一体式结构，则判定所述铺丝工装为整体结构类工装，并将所述铺丝工装转换为等效梁结构。
可选的，确定所述等效梁结构的载荷分布及截面尺寸参数，包括以下步骤：
若结构形式为主轴及其附属结构类工装，所述主轴为等截面梁，且各所述附属结构的质量沿所述铺丝工装的轴向方向均匀布置，则所述等效梁结构为等截面梁，将所述等截面梁的截面尺寸参数确定为所述主轴的截面尺寸参数；
若结构形式为主轴及其附属结构类工装，所述主轴为等截面梁， 且各所述附属结构的质量沿所述铺丝工装的轴向方向呈非均匀布置，所述等效梁结构为承受非均布载荷的等截面梁，将所述等效梁结构分割为若干分段等截面梁，并将各分段等截面梁截面中主轴的截面尺寸确定为各所述分段等截面梁的截面尺寸参数，各所述分段截面梁上的附属机构作为外载荷施加于所述等效梁结构上；
若结构形式为主轴及其附属结构类工装，所述主轴为变截面梁，则无论附属结构的质量如何分布，所述等效梁结构均为变截面梁，将所述等效梁结构分割为若干分段等截面梁，并将各分段等截面梁中点截面处主轴的截面尺寸确定为各所述分段等截面梁的截面尺寸参数，各所述分段截面梁上的附属机构作为外载荷施加于所述等效梁结构上；
若结构形式为整体结构类工装，且铺丝工装整体结构为等截面梁，则其等效梁结构为等截面梁，将所述等效梁结构的截面尺寸确定为铺丝工装的截面尺寸；否则为变截面梁，将所述等截面分割为若干分段等截面梁，将各分段截面梁中点截面对应的截面尺寸作为各所述分段等截面梁的截面尺寸参数。
可选的，分段等截面梁的长度u与主轴总长度U之间满足以下关系：u≤U(1-A)/10，其中A表示计算精度，其取值范围为0-1。
可选的，根据所述载荷分布和尺寸参数计算所述等效梁的载荷函数，包括以下步骤：
若结构形式为主轴及其附属结构类工装，则载荷函数为以下表达式：
Figure PCTCN2022125592-appb-000001
其中Z 0＝0，
Figure PCTCN2022125592-appb-000002
所述P i为集中力，Q j为惯性力，Q k为附属结构重量；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j对应的分段等截面梁各段起始和结束坐标；Z k-1、Z k分别为附属结构重力Q k对应的分段等截面梁各段的起始和结束坐标；i、j、k、r、m、n均为正整数，且均小于等于分段等截面梁的个数，i、j、k表示各分割段的编号，r、m、n表示分割段数量；a表示分割面坐标；z表示待求解点的坐标。
若等效梁为整体结构类工装，则载荷函数为以下表达式：
Figure PCTCN2022125592-appb-000003
其中，Z 0＝0，
Figure PCTCN2022125592-appb-000004
P i为集中力，Q j为惯性力；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j对应的分段等截面梁各段起始和结束坐标；i、j、r、m均为正整数，且均小于等于分段等截面梁的个数，i、j表示各分割段的编号，r、m表示分割段数量；a表示分割面坐标；z表示待求解点的坐标。
可选的，根据所述变形曲线函数计算等效梁结构的最大变形量，包括以下步骤：
确定等效梁结构两端处于铰支状态时的铰支边界条件；
确定等效梁结构两端处于固支状态是固支边界条件；
将所述铰支边界条件带入到变形曲线函数内计算铰支变形量，将固支边界条件带入到变形曲线函数内计算固支变形量；
比较铰支变形量和固支变形量，将所述铰支变形量和固支变形量中的较大者作为最大变形量的上限值，所述铰支变形量和固支变形量中的较小者作为最大变形量的下限值。
可选的，铰支边界条件为ω(0)＝ω(L)＝0；所述固支边界条件为ω(0)＝ω(L)＝0，d[ω(0)]/dz＝d[ω(L)]/dz＝0。
此外，为实现上述目的，本申请还提供一种自动铺丝工装的设计装置，包括：
等效转换模块，用于确定铺丝工装的结构形式，并根据所述结构形式将所述铺丝工装转换为等效梁结构；
数据生成模块，用于确定所述等效梁结构的载荷分布及截面尺寸参数；
第一计算模块，用于根据所述载荷分布和尺寸参数计算所述等效梁结构的载荷函数；
第二计算模块，用于根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
判断模块，结合所述变形量和是否满足铺丝工装的变形量判定基准，若所述变形量满足所述变形量判定基准，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤，直至所述变形量合格。
此外，为实现上述目的，本申请还提供一种电子设备，该电子设备包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器执行所述计算机程序，实现上述方法。
此外，为实现上述目的，本申请还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，处理器执行所述计算机程序，实现上述方法
与现有技术相比，本申请具有以下有益效果：
本申请首先确定铺丝工装的结构形式，再将根据所述结构形式将铺丝工装转换为等效梁结构，并结合铺丝工装的尺寸确定等效梁结构的载荷分布情况和截面尺寸，最后结合等效梁结构的载荷分布和变形曲线函数计算所述铺丝工装的变形范围；同时结合变形量的判定基准判定计算得到的铺丝工装的变形量是否合格；
与现有技术相比，本申请通过将复杂的铺丝工装转换为结构相对简单的等效梁结构，进而通过对等效梁结构的受力分析和变形分析反应所述铺丝工装的变形；其在保证计算精度的条件下，实现了对铺丝工装变形量计算过程的简化，因此在本申请所述技术方案中，只需要确定铺丝工装的结构形式、载荷分布和截面尺寸等关键参数即可进行计算，大大简化了进行变形计算需要的参数；因此，与现有技术中的实物测量法相比，本申请不需要实际制作铺丝工装的实物模型，其不但经济成本高，同时缩减了实物制作的时间，进一步缩短了铺丝工装的设计周期；
而与现有技术中的有限元仿真计算相比，本申请在铺丝工装技术方案和关键尺寸确定之后，在工装详细设计之前实现了对铺丝工装变形量的计算，其不但缩减了铺丝工装的详细设计周期，同时为详细设计过程中快速判断铺丝工装的结构形式和关键尺寸是否合理，快速判定铺丝工装的最大变形量等是否满足技术要求，为后续铺丝工装详细设计过程中各项参数的迭代优化提供了数据支撑和依据，从而进一步缩短了整个铺丝工装的设计周期。
本申请所述的自动铺丝工装的设计方法中所运用的等效梁结构建模及载荷计算方法，包含了工装的制造、装配误差，给出了误差最小和误差最大时各自对应的工装变形值，能够有效指导工程实际，对工装的极限情况进行准确评估。
附图说明
图1为本发明涉及的硬件运行环境的电子设备结构示意图；
图2为本申请实施例提供的一种自动铺丝工装的设计方法的流程图；
图3为本申请实施方式1所述铺丝工装结构示意图；
图4为本申请实施方式1中间芯轴剖视图；
图5为本申请实施方式1所述连接接头截面图；
图6为本申请实施方式1所述铺丝工装两端铰支状态受力分析图；
图7为本申请实施方式1所述铺丝工装两端固支状态受力分析图；
图8为本发明涉及的自动铺丝工装的设计装置的功能模块示意图；
图9为本申请实施方式2所述铺丝工装结构示意图；
图10为本申请实施方式2所述铺丝工装两端铰支状态受力分析图；
图11为本申请实施方式2所述铺丝工装两端固支状态受力分析图；
附图标记：1-接头，2-中间芯轴，3-上侧芯模，4-右侧芯模，5-下侧芯模，6-左侧芯模，7-工装，1001-处理器，1002-通信总线，1003-用户接口，1004-网络接口，1005-存储器。
本申请目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一 部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明，本发明实施例中所有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
在本发明中，除非另有明确的规定和限定，术语“连接”、“固定”等应做广义理解，例如，“固定”可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
另外，若本发明实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，全文中出现的“和/或”的含义，包括三个并列的方案，以“A和/或B”为例，包括A方案或B方案或A和B同时满足的方案。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
实施方式1
参照图1，图1为本申请实施例方案涉及的硬件运行环境的电子 设备结构示意图。
如图1所示，该电子设备可以包括：处理器1001，例如中央处理器(Central Processing Unit，CPU)，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(WIreless-FIdelity，WI-FI)接口)。存储器1005可以是高速的随机存取存储器(Random Access Memory，RAM)存储器，也可以是稳定的非易失性存储器(Non-Volatile Memory，NVM)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的结构并不构成对电子设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
参照图1，作为一种存储介质的存储器中可以包括操作系统、数据存储模块、网络通信模块、用户接口模块以及电子程序。
在图1所示的电子设备中，网络接口1004主要用于与网络服务器进行数据通信；用户接口1003主要用于与用户进行数据交互；本发明电子设备中的处理器1001、存储器1005可以设置在电子设备中，所述电子设备通过处理器1001调用存储器1005中存储的自动铺丝工装的设计装置，并执行本申请实施方式提供的自动铺丝工装设计方法。
如图2所示，本申请公开了一种自动铺丝工装的设计方法，包括以下步骤：
S1、确定铺丝工装的结构形式，并根据所述结构形式将所述铺丝 工装转换为等效梁结构；
如图3所示，铺丝工装包括位于中心位置的主轴，绕所述主轴的轴线，在所述主轴的四周均设置有附属结构，各附属结构与所述主轴相连；使用时，铺丝工装的主轴与外部动力设备相连；符合上述特征，则判定该铺丝工装为主轴及其附属结构类工装；
如上所述，本实施方式所述的铺丝工装为主轴及其附属结构类工装，从图3中可以看出，铺丝工装中部的芯轴两端与铺丝机旋转驱动接头相连，其余附属结构则不与铺丝机旋转驱动接头相连，因此将中部芯轴及铺丝机旋转驱动接头组成的构件转换为等效梁结构，而附属结构则等效为外载荷；
S2、确定所述等效梁结构的载荷分布及截面尺寸参数；
S21、若结构形式为主轴及其附属结构类工装，所述主轴为等截面梁，且各所述附属结构的质量沿所述铺丝工装的轴向方向均匀布置，则所述等效梁结构为等截面梁，等截面梁的截面尺寸参数为主轴的截面尺寸参数；
如图3所示，在本实施方式中，位于中部的主轴为等截面梁，同时沿所述主轴的轴线方向，各个附属结构的质量均匀分布，因此本实施方式以等截面梁的方式进行计算；
同时结合铺丝工装的工艺参数可知，该型铺丝工装加速旋转的最大角速度为A max＝π/6(rad/s)；铺丝机压头最大压力为100N，铺丝机压头位置位于工装中点处；要求计算精度为90％。
筒形件自动铺丝工装材料为INVAR钢，材料参数如表1所示。
表1 INVAR钢材料参数
Figure PCTCN2022125592-appb-000005
筒形件自动铺丝工装各部分的重量、重心如表2所示。
表2 筒形件自动铺丝工装各部分重量、重心
  重量(kg/个) 重心(mm)
接头 M1＝69.80 (0，0，250)(0，0，4650)
中间芯轴 M2＝619.04 (0，0，1950)
上侧芯模 M3＝1085.28 (200.33，-46.03，1950)
右侧芯模 M4＝1085.28 (46.03，200.33，1950)
下侧芯模 M5＝1085.28 (-200.33，46.03，1950)
左侧芯模 M6＝1085.28 (-46.03，-200.33，1950)
如图4和图5所示，其标注了筒形件自动铺丝工装主传力等效梁结构截面尺寸，同时上述尺寸的具体参数如表3所示。
表3 筒形件自动铺丝工装主传力等效梁结构截面尺寸
Figure PCTCN2022125592-appb-000006
同时所述铺丝工装的中部芯轴长为3900mm，铺丝机旋转驱动接头长度为500mm；
需要指出的是，如图3所示，以步骤S1中所述等效梁结构的一端为原点，等效梁结构的轴线方向为Z轴建立坐标系，其中所述X轴和Y轴均在所述等效梁结构的径向方向上；结合所述铺丝工装的中部芯轴长为3900mm，铺丝机旋转驱动接头长度为500mm等参数，可以得出表2中各坐标参数；
S3、根据所述载荷分布和尺寸参数计算所述等效梁结构的载荷函数；
由于本实施方式的铺丝工装类型为主轴及其附属结构类工装；同时结合受力分析可知，本事是方式的等效梁结构祝好包括中部芯轴的自重，各个附属芯膜的重量，中间芯轴两端受到的支撑力，铺丝机压头的集中力以及加速旋转的离心力；
因此所述载荷函数的表达式为：
Figure PCTCN2022125592-appb-000007
其中；Z 0＝0，
Figure PCTCN2022125592-appb-000008
同时结合步骤S21中的坐标系可知，在Z 1＝0mm，该处铺丝机对铺丝机旋转驱动接头的支撑力为P 1；等效梁结构另一端为Z 2＝4900mm，该处铺丝机对铺丝机旋转驱动接头的支撑力为P 2；铺丝机压头位于工装中点Z 3＝2450mm，压头压力为P 3＝100N；
两端铺丝机旋转驱动接头质量为M 1＝69.80kg，铺丝机旋转驱动接头截面形心距离旋转轴中心距离r 1＝0mm，其中一个连接接头的站 位坐标为Z 1＝0～Z 4＝500mm，另一个连接接头的站位坐标为Z 5＝4400mm～Z 2＝4900mm；
中间芯轴质量为M 2＝619.04kg，中间芯轴截面形心距离旋转轴中心距离r 2＝0mm；
上侧芯模、右侧芯模、下侧芯模、左侧芯模质量均为1085.28kg，上侧芯模、右侧芯模、下侧芯模、左侧芯模截面形心距离旋转轴中心距离均为r 3＝r 4＝r 5＝r 6＝205.55mm；中间芯轴及各外侧芯模的站位坐标均为Z 4＝500mm～Z 5＝4400mm；a表示分割面坐标；g为重力加速度，g＝9800mm/s 2。
结合上述参数可将上述展开为：
F(z)＝P 1<z-Z 1> -1+P 2<z-Z 2> -1-P 3<z-Z 3> -1
-M 1[g+r 1*A 2][<z-Z 1> 0-<z-Z 4> 0]/(Z 4-Z 1)
-M 1[g+r 1*A 2][<z-Z 5> 0-<z-Z 2> 0]/(Z 2-Z 5)
-M 2[g+r 2*A 2][<z-Z 4> 0-<z-Z 5> 0]/(Z 5-Z 4)
-M 3[g+r 3*A 2][<z-Z 4> 0-<z-Z 5> 0]/(Z 5-Z 4)
-M 4[g+r 4*A 2][<z-Z 4> 0-<z-Z 5> 0]/(Z 5-Z 4)
-M 5[g+r 5*A 2][<z-Z 4> 0-<z-Z 5> 0]/(Z 5-Z 4)
-M 6[g+6*A 2][<z-Z 4> 0-<z-Z 5> 0]/(Z 5-Z 4)
其中，Z 0＝0，
Figure PCTCN2022125592-appb-000009
S4、根据所述载荷函数得出等效梁结构的变形曲线函数，并根据 所述变形曲线函数计算等效梁结构的变形量；
S41、确定等效梁结构两端处于铰支状态时的铰支边界条件；
根据受力分析可得，当等效梁结构的两端处于铰支状态时，其边界条件为ω(0)＝ω(L)＝0；
S42、确定等效梁结构两端处于固支状态是固支边界条件；
根据受力分析可得，当等效梁结构的两端处于固支状态时，其边界条件ω(0)＝ω(L)＝0，d[ω(0)]/dz＝d[ω(L)]/dz＝0；
S43、将所述铰支边界条件带入到变形曲线函数内计算铰支变形量，将固支边界条件带入到变形曲线函数内计算固支变形量；
如图6所示，其标识了等效梁结构两端处于铰支状态时的受力状态，根据等效梁结构的变形挠度曲线函数：
Figure PCTCN2022125592-appb-000010
对步骤S3中的载荷函数进行计算可得：
Figure PCTCN2022125592-appb-000011
其中q 1＝M 2(g+r 2*A 2)/H+(M 3+M 4+M 5+M 6)*(g+r 3*A 2)/H＝12.4640N/mm，L 1＝3900mm；q 2＝M 1*(g+r 1*A 2)/h＝1.3681N/mm，L 2＝500mm；E 1＝E 2＝206000Mpa；I 1＝I(Z 4～Z 5)＝130307500mm 4，I 2＝I(Z 1～Z 4)＝I(Z 5～Z 2)＝53689000mm 4；
同时结合铰支边界条件ω(0)＝ω(L)＝0可以得出ω max＝3.2656mm；
如图7所示，其标识了等效梁结构两端处于固支状态时的受力状态，根据等效梁结构的变形挠度曲线函数：
Figure PCTCN2022125592-appb-000012
对步骤S3中的载荷函数进行计算可得，
Figure PCTCN2022125592-appb-000013
其中q 1＝M 2(g+r 2*A 2)/H+(M 3+M 4+M 5+M 6)*(g+r 3*A 2)/H＝12.4640N/mm，L 1＝3900mm；q 2＝M 1*(g+r 1*A 2)/h＝1.3681N/mm，L 2＝500mm；E 1＝E 2＝206000Mpa；I 1＝I(Z 4～Z 5)＝130307500mm 4，I 2＝I(Z 1～Z 4)＝I(Z 5～Z 2)＝53689000mm 4；
同时结合固支边界条件ω(0)＝ω(L)＝0，d[ω(0)]/dz＝d[ω(L)]/dz＝0，可以得出ω max＝1.1053mm。
S44、比较铰支变形量和固支变形量，其较大者作为最大变形量的上限值，较小这作为最大变形量的下限值。
根据上述计算，等效梁结构最大变形量的上限值为ω max＝3.2656mm，其下限位值为ω max＝1.1053mm，即1.1053mm<ω max<3.2656mm；
S5、结合所述变形量和铺丝工装的变形量判定基准，若所述变形量在判定基准允许范围内，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤直至所述变形量合格。
结合铺丝工装的设计要求制定所述铺丝工装的变形量判定基准， 取步骤S44中的变形量全部落入到变形量判定基准的范围内，则铺丝工装设计合格，否则需调整步骤S21中的载荷分布及截面尺寸参数，并重复上述计算步骤直至所述变形量合格。
参照图8，基于相同的发明思路，本申请的实施方式还提供一种飞机部件框梁间隙设计装置，包括：
等效转换模块，用于确定铺丝工装的结构形式，并根据所述结构形式将所述铺丝工装转换为等效梁结构；
数据生成模块，用于确定所述等效梁结构的载荷分布及截面尺寸参数；
第一计算模块，用于根据所述载荷分布和尺寸参数计算所述等效梁结构的载荷函数；
第二计算模块，用于根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
判断模块，结合所述变形量和是否满足铺丝工装的变形量判定基准，若所述变形量满足所述变形量判定基准，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤，直至所述变形量合格。
实施方式2
S1、确定铺丝工装的结构形式，根据所述结构形式将所述铺丝工装转换为等效梁结构；
如图9所示，铺丝工装没有主轴和安装于主轴上的附属结构，使用时，铺丝工装的两端分别与铺丝机旋转驱动接头相连；因此判定所述铺丝工装为整体结构类工装；
如上所述，本实施方式所述的铺丝工装为整体结构类工装，其不存在附属结构，因此将所述铺丝工装作为一个整体转换为等效梁结构；
S2、确定所述等效梁结构的载荷分布及截面尺寸参数；
S21、若结构形式为整体结构类工装，且铺丝工装整体结构为等截面梁，则其等效梁结构为等截面梁，所述等效梁结构的截面尺寸为铺丝工装的截面尺寸；否则为变截面梁，将所述等截面分割为若干相互拼接的分段等截面梁，各所述分段等截面梁的截面尺寸参数为各分段截面梁中点截面对应的截面尺寸；
如图9所示，在本实施方式中，沿所述铺丝工装的轴线方向，除所述铺丝工装的两端用于连接旋转驱动接头的部分为等截面结构外，主体上各个横截面面积不等，因此判定所述等效梁结构为变截面梁；
同时结合铺丝工装的工艺参数可知；铺丝机压头最大压力为100N，铺丝机压头位置位于工装中点处；要求计算精度为60％。
筒形件自动铺丝工装材料为INVAR钢，材料参数如表5所示。
表5 INVAR钢材料参数
Figure PCTCN2022125592-appb-000014
筒形件自动铺丝工装各部分的重量、重心如表6所示。
表6 筒形件自动铺丝工装各部分重量、重心
  重量(kg/个) 重心(mm)
接头 m 7＝54.263kg (0，0，250)
接头 m 8＝54.263kg (0，0，4650)
工装 m 9＝1507.168kg (0，0，2300)
需要指出的是，如图9所示，以所述等效梁结构的一端为原点，等效梁结构的轴线方向为Z轴建立坐标系，其中所述X轴和Y轴均在所述等效梁结构的径向方向上；结合所述铺丝工装的中部芯轴长为 3900mm，铺丝机旋转驱动接头长度为500mm等参数，可以得出表6中各坐标参数；
S22、将所述变截面等效梁结构分割为若干段等截面梁；
由于计算精度A＝60％，除两端的铺丝机旋转驱动接头外，工装主体部分的长度为3900mm，因此工装主体部分的分段等截面梁中各段的长度u与总长度U＝3900mm之间的关系满足：u≤U(1-A)/10＝156mm。取u＝156mm，将主传力等效梁结构分割为25段等截面梁，各段等截面梁的截面尺寸为各段中点截面所对应的截面尺寸，各段质量直接量取获得。
加上两端的铺丝机旋转驱动接头等截面梁，共为27段等截面梁，分段等截面梁各段截面尺寸及各段重量如表7所示。
表7 分段等截面梁各段截面
Figure PCTCN2022125592-appb-000015
S3、根据所述载荷分布和尺寸参数计算所述等效梁结构的载荷函数；
由于铺丝工装为整体结构类工装，因此其转换后的等效梁结构承受的载荷包含工装自重、两端支撑力、铺丝机压头集中力；
则载荷函数为以下表达式：
Figure PCTCN2022125592-appb-000016
其中Z 0＝0，
Figure PCTCN2022125592-appb-000017
由于等效梁结构被划分为27段，因此上述展开后的表达式为：
F(z)＝P 1<z-Z 0> -1+P 2<z-Z 27> -1-P 3<z-Z 14> -1-M 1g[<z-Z 0> 0-<z-Z 1> 0]/(Z 1-Z 0)
-M 2g[<z-Z 1> 0-<z-Z 2> 0]/(Z 2-Z 1)-M 3g[<z-Z 2> 0-<z-Z 3> 0]/(Z 3-Z 2)
-M 4g[<z-Z 3> 0-<z-Z 4> 0]/(Z 4-Z 3)-…-M 27g[<z-Z 26> 0-<z-Z 27> 0]/(Z 27-Z 26)
其中Z 0＝0，
Figure PCTCN2022125592-appb-000018
其中M 1～M 27为27段等截面梁中各段质量，Z 0～Z 27为各段等截面梁的起止点坐标；以主传力等效梁结构一端为原点，Z 0＝0mm，该处铺丝机对铺丝机旋转驱动接头的支撑力为P 1，主传力等效梁结构另一端为Z 27＝4900mm，该处铺丝机对铺丝机旋转驱动接头的支撑力为P 2；铺丝机压头位于工装中点Z 14＝2528mm处，压头压力为P 3＝100N；M 1～M 27为各段等截面梁的重量，其中a表示分割面坐标；g为重力加速度，g＝9800mm/s 2。
S4、根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
S41、确定等效梁结构两端处于铰支状态时的铰支边界条件；
根据受力分析可得，当等效梁结构的两端处于铰支状态时，其边界条件为ω(0)＝ω(Z 27)＝0，其中Z 27＝4900mm；
S42、确定等效梁结构两端处于固支状态是固支边界条件；
根据受力分析可得，当等效梁结构的两端处于固支状态时，其边 界条件ω(0)＝ω(Z 27)＝0，d[ω(0)]dz＝d[ω(Z 27)]dz＝0，其中Z 27＝4900mm；
S43、将所述铰支边界条件带入到变形曲线函数内计算铰支变形量，将固支边界条件带入到变形曲线函数内计算固支变形量；
根据材料力学的弯曲梁理论可知，等效梁结构的变形挠度曲线函数：
Figure PCTCN2022125592-appb-000019
其中F(z)为步骤S3中展开后的等效梁结构载荷函数，E(z)为材料弹性模量，本实施例中E(z)＝206000Mpa；I(z)为关于y轴的截面惯性矩，本实施例中各分段等截面梁的截面惯性矩如表8所示，C、D均为积分常数。
表8 各分段等截面梁的截面惯性矩
Figure PCTCN2022125592-appb-000020
Figure PCTCN2022125592-appb-000021
筒形件自动铺丝工装主传力等效梁结构的最大变形值ω max介于两端铰支和两端固支的等效梁结构的最大变形值之间，结合等效梁结构两端铰支的边界条件：ω(0)＝ω(Z 27)＝0，其中Z 27＝4900mm和等效梁结构两端固支的边界条件：ω(0)＝ω(Z 27)＝0，d[ω(0)]dz＝d[ω(Z 27)]dz＝0，其中Z 27＝4900mm；
同时由于本实施方式将等效梁结构划分对多段，因此为提高计算效率， 采用编程算法可得筒形件自动铺丝工装主传力等效梁结构的变形挠曲线最大值，如图10和图11所示，计算得到等效梁结构最大变形量的上限值为ω max＝0.645mm，其下限位值为ω max＝0.080mm，即0.080mm<ω max<0.645mm；
S5、结合所述变形量和铺丝工装的变形量判定基准，若所述变形量在判定基准允许范围内，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤直至所述变形量合格。
结合铺丝工装的设计要求制定所述铺丝工装的变形量判定基准，取步骤S43中的变形量全部落入到变形量判定基准的范围内，则铺丝工装设计合格，否则需调整步骤S21中的载荷分布及截面尺寸参数，并重复上述计算步骤直至所述变形量合格。
需要指出的是，本实施方式中的铺丝工装为整体结构类工装，且为变截面梁，如果整体结构类工装为等截面梁，则不需要将整个等效梁结构进行划分，只需要将等效梁结构作为一个整体进行载荷函数和变形挠度曲线函数的计算即可得出最终结果，即在实施方式1所述计算方法的基础上去掉附属结构的相关表达式进行计算。
同时在实施方式1中的铺丝工装为主轴及其附属结构类工装，且其为等截面梁，如果主轴及其附属结构类工装中，沿铺丝工装的轴向方向，中部芯轴或辅助结构中任意结构为变截面结构，则主轴及其附属结构类工装为变截面梁，需要按照实施方式2步骤S22中采用的方法对主轴进行划分，然后按照实施方式2记载的相关方法计算载荷函数和变形挠度曲线函数即可得到最终结果；
针对上述两种情况的具体计算过程，本说明书就不再进行详细说明；
与现有技术相比，本申请通过将复杂的铺丝工装转换为结构相对 简单的等效梁结构，进而通过对等效梁结构的受力分析和变形分析反应所述铺丝工装的变形；其在保证计算精度的条件下，实现了对铺丝工装变形量计算过程的简化，因此在本申请所述技术方案中，只需要确定铺丝工装的结构形式、载荷分布和截面尺寸等关键参数即可进行计算，大大简化了进行变形计算需要的参数；因此，与现有技术中的实物测量法相比，本申请不需要实际制作铺丝工装的实物模型，其不但经济成本高，同时缩减了实物制作的时间，进一步缩短了铺丝工装的设计周期；
而与现有技术中的有限元仿真计算相比，本申请在铺丝工装技术方案和关键尺寸确定之后，在工装详细设计之前实现了对铺丝工装变形量的计算，其不但缩减了铺丝工装的详细设计周期，同时为详细设计过程中快速判断铺丝工装的结构形式和关键尺寸是否合理，快速判定铺丝工装的最大变形量等是否满足技术要求，为后续铺丝工装详细设计过程中各项参数的迭代优化提供了数据支撑和依据，从而进一步缩短了整个铺丝工装的设计周期。
以上仅为本申请的优选实施例，并非因此限制本申请的专利范围，凡是利用本申请说明书及附图内容所做的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本申请的专利保护范围内。

Claims (8)

  1. 一种自动铺丝工装的设计方法，其特征在于，包括以下步骤：
    若铺丝工装包括主轴及若干安装于所述主轴上的附属结构，则判定所述铺丝工装为主轴及其附属结构类工装，并将所述主轴转换为等效梁结构，各所述附属结构作为外载荷施加于所述等效梁结构上；
    若铺丝工装为一体式结构，则判定所述铺丝工装为整体结构类工装，并将所述铺丝工装转换为等效梁结构；
    确定所述等效梁结构的载荷分布及截面尺寸参数；
    若结构形式为主轴及其附属结构类工装，则载荷函数为以下表达式：
    Figure PCTCN2022125592-appb-100001
    ，其中
    Figure PCTCN2022125592-appb-100002
    P i为集中力，Q j为惯性力，Q k为附属结构重量；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j对应的分段等截面梁各段起始和结束坐标；Z k-1、Z k分别为附属结构重力Q k对应的分段等截面梁各段的起始和结束坐标；i、j、k、r、m、n均为正整数，且均小于等于分段等截面梁的个数，i、j、k表示各分割段的编号，r、m、n表示分割段数量；a表示分割面坐标；
    若等效梁为整体结构类工装，则载荷函数为以下表达式：
    Figure PCTCN2022125592-appb-100003
    其中，
    Figure PCTCN2022125592-appb-100004
    P i为集中力，Q j为惯性力；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j 对应的分段等截面梁各段起始和结束坐标；i、j、r、m均为正整数，且均小于等于分段等截面梁的个数，i、j表示各分割段的编号，r、m表示分割段数量；a表示分割面坐标；
    根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
    结合所述变形量和是否满足铺丝工装的变形量判定基准，若所述变形量满足所述变形量判定基准，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤，直至所述变形量合格。
  2. 根据权利要求1所述的一种自动铺丝工装的设计方法，其特征在于，所述确定所述等效梁结构的载荷分布及截面尺寸参数，包括以下步骤：
    若结构形式为主轴及其附属结构类工装，所述主轴为等截面梁，且各所述附属结构的质量沿所述铺丝工装的轴向方向均匀布置，则所述等效梁结构为等截面梁，将所述等截面梁的截面尺寸参数确定为所述主轴的截面尺寸参数；
    若结构形式为主轴及其附属结构类工装，所述主轴为等截面梁，且各所述附属结构的质量沿所述铺丝工装的轴向方向呈非均匀布置，所述等效梁结构为承受非均布载荷的等截面梁，将所述等效梁结构分割为若干分段等截面梁，并将各分段等截面梁截面中主轴的截面尺寸确定为各所述分段等截面梁的截面尺寸参数，各所述分段截面梁上的附属机构作为外载荷施加于所述等效梁结构上；
    若结构形式为主轴及其附属结构类工装，所述主轴为变截面梁，则无论附属结构的质量如何分布，所述等效梁结构均为变截面梁，将所述等效梁结构分割为若干分段等截面梁，并将各分段等截面梁中点 截面处主轴的截面尺寸确定为各所述分段等截面梁的截面尺寸参数，各所述分段截面梁上的附属机构作为外载荷施加于所述等效梁结构上；
    若结构形式为整体结构类工装，且铺丝工装整体结构为等截面梁，则其等效梁结构为等截面梁，将所述等效梁结构的截面尺寸确定为铺丝工装的截面尺寸；否则为变截面梁，将所述等截面分割为若干分段等截面梁，将各分段截面梁中点截面对应的截面尺寸作为各所述分段等截面梁的截面尺寸参数。
  3. 根据权利要求2所述的一种自动铺丝工装的设计方法，其特征在于，所述分段等截面梁的长度与主轴总长度之间满足以下关系：u≤U(1-A)/10，其中A表示计算精度，其取值范围为0-1。
  4. 根据权利要求1所述的一种自动铺丝工装的设计方法，其特征在于，所述根据所述变形曲线函数计算等效梁结构的最大变形量，包括以下步骤：
    确定等效梁结构两端处于铰支状态时的铰支边界条件；
    确定等效梁结构两端处于固支状态是固支边界条件；
    将所述铰支边界条件带入到变形曲线函数内计算铰支变形量，将固支边界条件带入到变形曲线函数内计算固支变形量；
    比较铰支变形量和固支变形量，其较大者作为最大变形量的上限值，较小这作为最大变形量的下限值。
  5. 根据权利要求4所述的一种自动铺丝工装的设计方法，其特征在于，所述铰支边界条件为ω(0)＝ω(L)＝0；所述固支边界条件为ω(0)＝ω(L)＝0，d[ω(0)]/dz＝d[ω(L)]/dz＝0。
  6. 一种自动铺丝工装的设计装置，其特征在于，包括：
    等效转换模块，用于在铺丝工装包括主轴及若干安装于所述主轴 上的附属结构时，判定所述铺丝工装为主轴及其附属结构类工装，并将所述主轴转换为等效梁结构，各所述附属结构作为外载荷施加于所述等效梁结构上；
    或在铺丝工装为一体式结构时，判定所述铺丝工装为整体结构类工装，并将所述铺丝工装转换为等效梁结构；
    数据生成模块，用于确定所述等效梁结构的载荷分布及截面尺寸参数；
    第一计算模块，用于在结构形式为主轴及其附属结构类工装时，则确定载荷函数为以下表达式：
    Figure PCTCN2022125592-appb-100005
    ，其中
    Figure PCTCN2022125592-appb-100006
    P i为集中力，Q j为惯性力，Q k为附属结构重量；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j对应的分段等截面梁各段起始和结束坐标；Z k-1、Z k分别为附属结构重力Q k对应的分段等截面梁各段的起始和结束坐标；i、j、k、r、m、n均为正整数，且均小于等于分段等截面梁的个数，i、j、k表示各分割段的编号，r、m、n表示分割段数量；a表示分割面坐标；
    或在等效梁为整体结构类工装，则确定载荷函数为以下表达式：
    Figure PCTCN2022125592-appb-100007
    其中，
    Figure PCTCN2022125592-appb-100008
    P i为集中力，Q j为惯性力；Z i为集中力所在坐标；Z j-1、Z j分别为惯性力Q j对应的分段等截面梁各段起始和结束坐标；i、j、r、m均为正整数， 且均小于等于分段等截面梁的个数，i、j表示各分割段的编号，r、m表示分割段数量；a表示分割面坐标；
    第二计算模块，用于根据所述载荷函数得出等效梁结构的变形曲线函数，并根据所述变形曲线函数计算等效梁结构的变形量；
    判断模块，结合所述变形量和是否满足铺丝工装的变形量判定基准，若所述变形量满足所述变形量判定基准，则所述铺丝工装合格，否则重复确定所述等效梁结构的载荷分布及截面尺寸参数的步骤，直至所述变形量合格。
  7. 一种电子设备，其特征在于，该电子设备包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器执行所述计算机程序，实现如权利要求1-5中任一项所述的方法。
  8. 一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，处理器执行所述计算机程序，实现如权利要求1-5中任一项所述的方法。
PCT/CN2022/125592 2022-06-17 2022-10-17 一种自动铺丝工装的设计方法、装置、设备及存储介质 Ceased WO2023240860A1 (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210684191.2 2022-06-17
CN202210684191.2A CN114770977B (zh) 2022-06-17 2022-06-17 一种自动铺丝工装的设计方法、装置、设备及存储介质

Publications (1)

Publication Number Publication Date
WO2023240860A1 true WO2023240860A1 (zh) 2023-12-21

Family

ID=82421006

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/125592 Ceased WO2023240860A1 (zh) 2022-06-17 2022-10-17 一种自动铺丝工装的设计方法、装置、设备及存储介质

Country Status (2)

Country Link
CN (1) CN114770977B (zh)
WO (1) WO2023240860A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118024623A (zh) * 2024-04-15 2024-05-14 杭州艾美依航空制造装备有限公司 多功能自动铺丝系统
CN119952989A (zh) * 2025-04-01 2025-05-09 中建材(上海)航空技术有限公司 一种复合材料曲面结构的二维料片展开及成型方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114770977B (zh) * 2022-06-17 2022-10-25 成都飞机工业（集团）有限责任公司 一种自动铺丝工装的设计方法、装置、设备及存储介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000328645A (ja) * 1999-05-24 2000-11-28 Nippon Kaiser Kk 格子梁解析システム
CN104636543A (zh) * 2014-12-31 2015-05-20 哈尔滨工业大学 一种基于有限差分法的重型龙门铣床横梁重力变形预测方法
CN104714478A (zh) * 2014-12-31 2015-06-17 哈尔滨工业大学 基于有限差分法的重型双柱立车横梁重力变形预测方法
CN106273545A (zh) * 2016-09-29 2017-01-04 西安交通大学 一种实现预浸带输送和铺放速度一致的重送方法及装置
CN106339523A (zh) * 2016-06-24 2017-01-18 北京玻钢院复合材料有限公司 一种大扭矩复合材料传动轴结构的设计方法
WO2018151298A1 (ja) * 2017-02-17 2018-08-23 新日鐵住金株式会社 分析方法、設計方法、製造方法、及びプログラム
CN109484624A (zh) * 2018-11-29 2019-03-19 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种飞机复合材料机翼翼梁及翼根区连接结构
KR20190064992A (ko) * 2017-12-01 2019-06-11 김승원 합성보의 구조 설계 수행 방법 및 컴퓨터로 판독 가능한 저장매체
CN111460710A (zh) * 2020-03-30 2020-07-28 成都飞机工业（集团）有限责任公司 一种基于铺丝轨迹的复材固化变形仿真建模方法
CN114488814A (zh) * 2022-01-26 2022-05-13 中国航空制造技术研究院 一种消除重力扰动高精度铺放张力控制方法
CN114770977A (zh) * 2022-06-17 2022-07-22 成都飞机工业（集团）有限责任公司 一种自动铺丝工装的设计方法、装置、设备及存储介质

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE159476T1 (de) * 1990-08-09 1997-11-15 First Green Park Pty Ltd Verpackungssystem
US5904972A (en) * 1995-06-07 1999-05-18 Tpi Technology Inc. Large composite core structures formed by vacuum assisted resin transfer molding
JP4186334B2 (ja) * 1999-09-13 2008-11-26 日新イオン機器株式会社 イオン注入装置
US6655633B1 (en) * 2000-01-21 2003-12-02 W. Cullen Chapman, Jr. Tubular members integrated to form a structure
US20140372082A1 (en) * 2000-12-01 2014-12-18 Aleksandr I. KAMENOMOSTSKIY Tool for optimized thin wall profile member (tpm) and tpm-panel design and selection
US9278484B2 (en) * 2008-04-17 2016-03-08 The Boeing Company Method and apparatus for producing contoured composite structures and structures produced thereby
JP5368416B2 (ja) * 2010-11-18 2013-12-18 日本車輌製造株式会社 分岐装置
CN102479260A (zh) * 2010-11-23 2012-05-30 大连兆阳软件科技有限公司 一种对梁进行精确分析的智能软件
DE102010063165A1 (de) * 2010-12-15 2012-06-21 Hochschule Merseburg (Fh) Druckbasierte Flächenausgleichsmodule zur Aufnahme von auf Baukonstruktionen wirkenden Kräften mittels sogenannter Fluidkolben
CN102108946B (zh) * 2011-01-17 2013-01-09 南京航空航天大学 复合铺层式风力机叶片及其制造方法
US10169492B2 (en) * 2011-06-20 2019-01-01 The Boeing Company Fiber placement optimization for steered-fiber plies
CN105631073B (zh) * 2014-10-31 2019-05-21 中国航空工业集团公司西安飞机设计研究所 一种优化制造复合材料波纹梁的方法
US10195817B2 (en) * 2015-01-02 2019-02-05 The Boeing Company Skin-stringer design for composite wings
US10434726B1 (en) * 2015-07-13 2019-10-08 The Boeing Company Forming thermoplastic composite parts having steered fiber orientations
CN105760629A (zh) * 2016-03-25 2016-07-13 兰州理工大学 风力机叶片主梁铺层优化设计方法
CN106626434B (zh) * 2016-11-29 2019-03-08 中国航空工业集团公司沈阳飞机设计研究所 一种飞机复合材料隔框结构制造及设计方法
CN106599485B (zh) * 2016-11-29 2020-10-09 中国航空工业集团公司沈阳飞机设计研究所 一种框梁整体件优化设计方法
CN106515045B (zh) * 2016-12-18 2019-02-15 湖北三江航天江北机械工程有限公司 复合材料连接裙自动铺放成型方法
CN108268692A (zh) * 2017-11-21 2018-07-10 北京航空航天大学 自动铺丝复合材料固化变形预报方法
CN108647386A (zh) * 2018-04-04 2018-10-12 江苏理工学院 一种应用于frp防撞梁结构/材料的协同优化设计方法
US10974465B2 (en) * 2018-06-20 2021-04-13 The Boeing Company Method and system for generating a layup plan for forming a composite laminate
ES2857911T3 (es) * 2018-07-19 2021-09-29 Airbus Operations Slu Pieza de composite curva y su método de fabricación
US20200023569A1 (en) * 2018-07-20 2020-01-23 Valorbec Societe En Commandite Method and system for 4d printing of composites
US11084223B2 (en) * 2018-08-20 2021-08-10 Lawrence Livermore National Security, Llc Optimal toolpath generation system and method for additively manufactured composite materials
US20200173881A1 (en) * 2018-12-03 2020-06-04 The Boeing Company Method for Determining Noodle and Interface Residual Stresses and Properties Under Hygrothermal-Mechanical Loadings
CN109573476B (zh) * 2019-01-25 2024-01-05 深圳市泰克光电科技有限公司 一种灯丝检测分光机
US11807355B2 (en) * 2019-06-28 2023-11-07 The Boeing Company Landing gear system with composite flex beam
CN111639419B (zh) * 2020-05-09 2022-06-21 哈尔滨工程大学 一种船舱内爆小变形塑性毁伤模式问题的判断方法
CN112123811B (zh) * 2020-07-24 2022-08-12 北京卫星制造厂有限公司 一种分体工装复材支架的成型方法
CN111859690B (zh) * 2020-07-27 2022-10-18 吉林重通成飞新材料股份公司 一种梁单元模拟的风电叶片静力测试方法、装置及设备
CN112487676B (zh) * 2020-11-13 2022-09-06 福州大学 一种深海管道s型铺设动力模拟方法
CN113408105B (zh) * 2021-05-24 2022-10-18 沃飞长空科技(成都)有限公司 曲面结构上材料应力状态的确定方法及装置
CN113486417B (zh) * 2021-05-29 2024-01-19 中铁建工集团有限公司 一种控制主塔和钢梁施工仿真系统

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000328645A (ja) * 1999-05-24 2000-11-28 Nippon Kaiser Kk 格子梁解析システム
CN104636543A (zh) * 2014-12-31 2015-05-20 哈尔滨工业大学 一种基于有限差分法的重型龙门铣床横梁重力变形预测方法
CN104714478A (zh) * 2014-12-31 2015-06-17 哈尔滨工业大学 基于有限差分法的重型双柱立车横梁重力变形预测方法
CN106339523A (zh) * 2016-06-24 2017-01-18 北京玻钢院复合材料有限公司 一种大扭矩复合材料传动轴结构的设计方法
CN106273545A (zh) * 2016-09-29 2017-01-04 西安交通大学 一种实现预浸带输送和铺放速度一致的重送方法及装置
WO2018151298A1 (ja) * 2017-02-17 2018-08-23 新日鐵住金株式会社 分析方法、設計方法、製造方法、及びプログラム
KR20190064992A (ko) * 2017-12-01 2019-06-11 김승원 합성보의 구조 설계 수행 방법 및 컴퓨터로 판독 가능한 저장매체
CN109484624A (zh) * 2018-11-29 2019-03-19 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种飞机复合材料机翼翼梁及翼根区连接结构
CN111460710A (zh) * 2020-03-30 2020-07-28 成都飞机工业（集团）有限责任公司 一种基于铺丝轨迹的复材固化变形仿真建模方法
CN114488814A (zh) * 2022-01-26 2022-05-13 中国航空制造技术研究院 一种消除重力扰动高精度铺放张力控制方法
CN114770977A (zh) * 2022-06-17 2022-07-22 成都飞机工业（集团）有限责任公司 一种自动铺丝工装的设计方法、装置、设备及存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YANG HUIFANG, LI ZHONGKAI, JIN HUI: "Analysis of Key Component of Large Composites Fiber Placement Machine", AERONAUTICAL MANUFACTURING TECHNOLOGY, no. 22, 15 November 2010 (2010-11-15), pages 72 - 75, XP093117616, ISSN: 1671-0833, DOI: 10.16080/j.issn1671-833x.2010.22.015 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118024623A (zh) * 2024-04-15 2024-05-14 杭州艾美依航空制造装备有限公司 多功能自动铺丝系统
CN119952989A (zh) * 2025-04-01 2025-05-09 中建材(上海)航空技术有限公司 一种复合材料曲面结构的二维料片展开及成型方法

Also Published As

Publication number Publication date
CN114770977B (zh) 2022-10-25
CN114770977A (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
WO2023240860A1 (zh) 一种自动铺丝工装的设计方法、装置、设备及存储介质
CN111141518B (zh) 一种基于模型的非对称转子轴承系统不平衡量识别方法
CN114858407A (zh) 一种静气动弹性风洞试验机翼模型及复合材料蒙皮
CN110069822A (zh) 一种叶片动应变测量的传感器布置方法
CN114722686B (zh) 一种基于有限元分析的大型设备吊耳设计及优化方法
CN116754134A (zh) 基于试验与仿真数据融合的转子不平衡状态识别方法
CN117272639A (zh) 一种用于桥梁顶推施工受力分析方法
CN102063539B (zh) 一种基于有限元的惯性平台残余应力释放仿真方法
CN110807276A (zh) 导入残余应力分析铸件结构强度的方法及系统
CN113946994A (zh) 一种基于数字孪生体的光滑有限元数值计算方法
CN110160740A (zh) 一种一体化冲击气动力测量系统
CN117932719A (zh) 一种求解悬索结构的非线性分析方法
CN117828833A (zh) 一种cmc结构热梯度环境下的力学性能预测方法
CN103745038A (zh) 摩天轮结构有限元模型的建模方法
CN111881606B (zh) 一种考虑截面变形的伸缩臂臂架动力学建模方法
CN111444644A (zh) 一种基于有限元技术的工业机器人整机刚度优化方法
CN111639386A (zh) 一种航空发动机用无键相虚拟动平衡方法
CN114936483B (zh) 大型轴承静力学仿真分析的有限元等效建模方法
CN115618163B (zh) 一种加筋结构体固有频率计算方法
CN116305646A (zh) 一种基于柔性基础传递函数的整机转子动力学瞬态计算方法及装置
CN116150835B (zh) 直角坐标系下抛物线拱面内自由振动频率近似解析方法
CN115130299A (zh) 一种单支点柔壁喷管型面曲线设计方法
CN114218823A (zh) 基于三维模型的运载火箭疲劳载荷赋值系统及方法
CN109989890B (zh) 风力发电机组钢制筒形塔架设计方法、装置、设备及介质
CN119334788B (zh) 一种汽轮机转子蠕变弯曲预测方法及蠕变弯曲预测装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22946532

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 22946532

Country of ref document: EP

Kind code of ref document: A1