WO2014201984A1 - Super-unit construction method based structure-function analysis method and device therefor - Google Patents

Super-unit construction method based structure-function analysis method and device therefor Download PDF

Info

Publication number
WO2014201984A1
WO2014201984A1 PCT/CN2014/079933 CN2014079933W WO2014201984A1 WO 2014201984 A1 WO2014201984 A1 WO 2014201984A1 CN 2014079933 W CN2014079933 W CN 2014079933W WO 2014201984 A1 WO2014201984 A1 WO 2014201984A1
Authority
WO
WIPO (PCT)
Prior art keywords
super
displacement
matrix
structural
deformation
Prior art date
Application number
PCT/CN2014/079933
Other languages
French (fr)
Chinese (zh)
Inventor
丁桦
Original Assignee
深圳市网蓝实业有限公司
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 深圳市网蓝实业有限公司 filed Critical 深圳市网蓝实业有限公司
Publication of WO2014201984A1 publication Critical patent/WO2014201984A1/en

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
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Definitions

  • the present invention relates to the field of structural dynamics, and more particularly to a structural function analysis method based on a super-cell construction method and an apparatus therefor.
  • the structural function parameters are generally determined by a certain assumption of the deformation of the structure, or by a certain approximation, establishing the relationship between the equivalent generalized force and the generalized deformation.
  • the substructure method can be used to define the structural function parameters: establish the relationship between the external nodes of the substructure part and the external node forces.
  • the corresponding structural function parameter is the stiffness matrix of the substructure.
  • the traditional substructure method is difficult to give the relationship between the external nodes and the nodal forces of the substructure, so it is difficult to establish a structure similar to the structural static problem. Function definition.
  • model simplification is also a very critical technique in structural dynamic analysis.
  • the fundamental purpose of the structural dynamics model simplification is to obtain a low-order, efficient computational model that meets the engineering accuracy requirements, so that the simplified model can be used for performance analysis and simulation of the original complex model.
  • the existing structural dynamic model simplification methods can be mainly divided into the following three categories.
  • Dynamic substructure methods which are methods that directly obtain low-order models. Firstly, the low-order dynamic characteristics of each substructure are obtained, and then the comprehensive vibration equation of the whole structure expressed by low-order modal coordinates is obtained by the double coordination condition of displacement and force between substructures.
  • the new dynamic model simplification technology based on the deformation-corrected power reduction method, can better solve the various limitations and inconveniences in the traditional model simplification technology while maintaining high efficiency. Under certain assumptions, this method can be used to define and simplify the structural dynamics, stiffness and quality of substructures.
  • the invention name of the application No. 200810102136.8 is "power reduction algorithm based on deformation correction and super unit construction method”.
  • Chinese patent application discloses a power reduction algorithm based on deformation correction and a super unit construction method, which is hereby incorporated by reference in its entirety. Although this method can simplify the model of complex structure well, the computational complexity is greatly reduced. But right In terms of structural parts, it contains relatively many parameters, so that it is rather cumbersome to characterize the structure.
  • an object of the present invention is to provide a structural function analysis method based on a super-cell construction method, which simplifies parameters.
  • Step 2 According to the displacement characteristics of the boundary of the structural member, the displacement of the boundary node is further simplified, namely:
  • g is a super-unit generalized displacement vector
  • 7j and ⁇ are respectively a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member; using ⁇ to construct a structural member stiffness similar to the super-cell method, Mass and damping matrices, which are used as mechanical functions or characteristic parameters of structural members.
  • S is a deformation feature matrix that can directly select the component; when further precision is required, it can be obtained by optimization similar to the method of constructing ⁇ .
  • the structure is divided into several parts, each part is regarded as a super unit, and then divided into several motion synchronization areas according to the motion characteristics of each part, the part boundary nodes are reserved and the feature points of each part of the synchronization area are Add a super node.
  • the super unit construction method is then used to generate a simplified super unit model for the structural unit.
  • the present invention further simplifies the system based on the super-element construction method using the deformation feature matrix.
  • the calculation processing of the structural unit by this method can greatly simplify the parameters of the structural features and better maintain the accuracy of the calculation.
  • the simplest case is to directly use the aforementioned dynamic substructure method, the mass matrix of the substructure, the damping matrix of the substructure ( ', the stiffness matrix A ' of the substructure, directly as the functional parameters of the substructure;
  • the substructure In the case that the substructure is more complicated, the substructure can be divided into several smaller substructures.
  • the structural function parameters of each substructure are obtained by the above method, and the overall matrix integration is performed by splicing, that is, similar to finite element. Method, obtaining functional parameters of the substructure;
  • the present invention also provides an apparatus using the above method, the structural function analysis apparatus based on the super unit construction method comprising the following structural unit: a super unit model forming unit, which simplifies the structural unit generation according to the super unit construction method Super element model, an optimized approximation of the displacement approximating the original finite element model
  • is the generalized displacement vector of the superelement
  • is the optimal displacement mode matrix
  • the corresponding stiffness, mass and damping matrix can be obtained by transformation
  • the node displacement simplification unit the displacement of the boundary node according to the displacement feature of the structural member boundary
  • is a super-unit generalized displacement vector
  • is a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member
  • the structural member can be constructed by using a super-cell-like method Stiffness, mass and damping matrix, which are used as mechanical functions or characteristic parameters of structural members.
  • s is a deformation feature matrix of the component directly selected; when it is required to further improve the precision, the structure and the tau can be similarly constructed.
  • the method is obtained by optimization.
  • the present invention further simplifies the system based on the superelement construction method using the deformation feature matrix.
  • the calculation processing of the structural unit by this method can greatly simplify the parameters of the structural features and better maintain the accuracy of the calculation.
  • the simplest case is to directly use the aforementioned dynamic substructure method, the mass matrix of the substructure, the damping matrix of the substructure ( ', the stiffness matrix A ' of the substructure, directly as the functional parameters of the substructure;
  • the substructure can be divided into several smaller substructures, and the structural function parameters of each substructure are obtained by the foregoing method, and the overall matrix integration is performed by splicing, that is, similar finite element. Method, obtaining functional parameters of the substructure;
  • the invention will utilize a new dynamic model simplification technique, based on the dynamic reduction method of deformation correction, and give a dynamic substructure method, thereby realizing the characterization of structural function parameters.
  • FIG. 1 is a substructure mass condensation model obtained by a super unit analysis method
  • FIG. 2 is a frequency error curve of the embodiment of FIG. 1;
  • Figure 3 is a super unit obtained by a super unit construction method
  • Figure 4 is an overall structure composed of super units
  • Figure 5 is a structural unit obtained by simplification;
  • Figure 6 is an overall structure composed of structural units.
  • the more compact lung type of the present invention is a structural function analysis method of the present invention, specifically:
  • the simplified super-cell model is generated by the structural unit, and an optimized approximation of the displacement of the original finite element model is obtained, where g is the generalized displacement vector of the super element, and T is the optimized displacement mode matrix.
  • the corresponding stiffness, mass and damping matrix can be obtained by transformation;
  • q Bq , where g is a super-unit generalized displacement vector, 7 ⁇ , ⁇ are respectively a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member; using ⁇ to construct a structural member stiffness similar to the super-cell method, Mass and damping matrices, which are used as mechanical functions or characteristic parameters of structural members.
  • the deformed feature matrix of the component can be directly selected; when the precision needs to be further improved, the method can be similarly optimized by the method of constructing ⁇ .
  • the structural function analysis device based on the super unit construction method includes the following structural unit:
  • the super-cell model forming unit generates a simplified super-cell model from the structural unit according to the super-cell construction method, and obtains an optimized approximation ⁇ which approximates the displacement of the original finite element model, where is the generalized displacement vector of the super-unit, and ⁇ is optimized Displacement mode matrix, by which the corresponding stiffness, mass and damping matrix can be obtained;
  • the node displacement simplification unit further simplifies the boundary node displacement according to the displacement feature of the structural member boundary, that is, C, wherein the super-unit generalized displacement vector J and s are respectively the generalized displacement vector and the component deformation determined according to the deformation characteristics of the member.
  • Characteristic moment Array the structural stiffness, mass and damping matrix of structural members can be constructed by using S-like super-cell method, and they are used as mechanical functions or characteristic parameters of structural members.
  • is the deformation feature matrix of the component directly selected; when the precision needs to be further improved, it can be obtained by optimization similar to the method of constructing T.
  • the quality of each dynamic substructure is concentrated in the internal super node.
  • the displacement of both sides of the substructure shown in Fig. 2 is the translational displacement and rotation angle of the two midpoints (u, w, ⁇ ) (al )
  • the present invention further simplifies the system based on the super-element construction method using the deformation feature matrix. It can better solve various limitations and inconveniences in the traditional model simplification technology. Under certain assumptions, the structural dynamics function, stiffness and quality of the substructure can be performed by calculating the structural unit by this method. It is defined and greatly simplified, and it is better able to maintain the accuracy of the calculation while maintaining high efficiency.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The present invention relates to a structure-function analysis method and device therefor, based on a super-unit construction method. The present method and device, based on a super-unit construction method, use a deformation characteristic matrix to further simplify a system. Computational processing of structural units greatly simplifies the parameters of the structural characteristics, and better maintains computational accuracy.

Description

一种基于超单元构造方法的结构功能分析方法及其装置  Structural function analysis method based on super unit construction method and device thereof
技术领域 本发明涉及结构动力学领域,尤其是一种基于超单元构造方法的 结构功能分析方法及其装置。 TECHNICAL FIELD The present invention relates to the field of structural dynamics, and more particularly to a structural function analysis method based on a super-cell construction method and an apparatus therefor.
背景技术 结构力学领域所涉及的任何一个结构都具有它特定的力学功能。 长期以来人们对简单结构通常都已其特有的力学功能来描述它们,如: 弹簧,在使用方向上用其整体的弹簧系数来描述;杆、 梁均有相应的 力学特征参数来描述等。这些描述使问题变得简单、 明了 ,将结构的 功能用较直接的参数表述出来,为应用带来极大的便利和高效。但对 于较复杂的结构要对其功能进行简洁、有效的描述将是一件非常困难 的事。其中的关键,是在保证一定精度的前提下如何用尽可能少的参 数来描述结构的功能。 BACKGROUND OF THE INVENTION Any structure involved in the field of structural mechanics has its specific mechanical function. For a long time, simple structures have been described by their unique mechanical functions, such as: springs, which are described by their overall spring coefficient in the direction of use; rods and beams have corresponding mechanical characteristic parameters to describe. These descriptions make the problem simple and straightforward, and express the function of the structure with relatively straightforward parameters, which brings great convenience and efficiency to the application. But a simple and effective description of its function for more complex structures would be very difficult. The key to this is how to describe the function of the structure with as few parameters as possible while maintaining a certain degree of precision.
对于较简单的结构,在特定载荷作用下,其结构功能参数一般是 通过对结构的变形进行一定的假设限制,或者说进行一定的近似,建 立等效广义力和广义变形之间的关系来得到的 ,如,梁、 板、 壳等。 由于假设的存在,这些结构在使用中要受到诸多的限制。从另一个角 度看,当结构和载荷情况与假设相差较大时,结果的精度就满足不了 要求了。  For a simpler structure, under the action of a specific load, the structural function parameters are generally determined by a certain assumption of the deformation of the structure, or by a certain approximation, establishing the relationship between the equivalent generalized force and the generalized deformation. Such as, beams, plates, shells, etc. Due to the existence of assumptions, these structures are subject to many limitations in their use. From another perspective, when the structure and load conditions differ greatly from the assumptions, the accuracy of the results does not meet the requirements.
对于结构静力学问题,在一定的条件下,可以利用子结构方法来 进行结构功能参数的定义:建立子结构部分外部节点和外部节点力之 间的关系。 对应的结构功能参数就是子结构的刚度矩阵。  For the structural statics problem, under certain conditions, the substructure method can be used to define the structural function parameters: establish the relationship between the external nodes of the substructure part and the external node forces. The corresponding structural function parameter is the stiffness matrix of the substructure.
对于结构动力学问题,传统的子结构方法很难给出子结构部分外 部节点和节点力之间的关系,故很难建立类似结构静力学问题的结构 功能定义。 For the structural dynamics problem, the traditional substructure method is difficult to give the relationship between the external nodes and the nodal forces of the substructure, so it is difficult to establish a structure similar to the structural static problem. Function definition.
动态系统模型简化方法的研究和应用 ,一直是大型动态系统理论 研究和结构设计中的重要课题。 同样,在结构动力学领域,模型简化 也是结构动力分析中非常关键的技术。结构动力学模型简化的根本目 的是为了获得一个满足工程精度要求的低阶、有效的计算模型,从而 可以用简化模型对原始复杂模型进行性能分析以及模拟仿真。现有的 结构动力模型简化方法主要可分成以下三类。  The research and application of the simplified method of dynamic system model has always been an important topic in the theoretical research and structural design of large dynamic systems. Similarly, in the field of structural dynamics, model simplification is also a very critical technique in structural dynamic analysis. The fundamental purpose of the structural dynamics model simplification is to obtain a low-order, efficient computational model that meets the engineering accuracy requirements, so that the simplified model can be used for performance analysis and simulation of the original complex model. The existing structural dynamic model simplification methods can be mainly divided into the following three categories.
( I ) 自由度减缩方法,基本思路是从结构一般运动方程或特征 方程出发,用保留自由度表示縮聚掉的自由度,从而实现对模型的简 化。典型的此类方法有 Guyan-Irons法、 Kuhar法、 IRS法和模态縮聚 法等。  (I) The method of reducing the degree of freedom, the basic idea is to start from the general motion equation or characteristic equation of the structure, and use the degree of freedom to express the degree of freedom of the polycondensation, so as to simplify the model. Typical such methods are the Guyan-Irons method, the Kuhar method, the IRS method, and the mode polycondensation method.
( 2 )动态子结构方法,这类方法是直接得到低阶模型的方法。 首先得到各个子结构的低阶动力特性,然后通过子结构间位移和力的 双协调条件得到整体结构以低阶模态坐标表示的综合振动方程。  (2) Dynamic substructure methods, which are methods that directly obtain low-order models. Firstly, the low-order dynamic characteristics of each substructure are obtained, and then the comprehensive vibration equation of the whole structure expressed by low-order modal coordinates is obtained by the double coordination condition of displacement and force between substructures.
( 3 )结构等效方法,此类方法从结构力学分析出发 ,针对某种 特定的结构,用简单的结构等效复杂结构的主要特征,得到简化的力 学模型。  (3) Structural equivalence method. From the structural mechanics analysis, such a method obtains a simplified force model for a specific structure and uses the main features of a simple structural equivalent complex structure.
新的动态模型简化技术,基于变形修正的动力减縮方法,能够较 好地解决传统模型简化技术中的各种限制和不便,同时又能够保持较 高的效率。在一定的假设下,利用这种方法就可以对子结构的结构动 力学功能,刚度和质量,进行定义和简化。  The new dynamic model simplification technology, based on the deformation-corrected power reduction method, can better solve the various limitations and inconveniences in the traditional model simplification technology while maintaining high efficiency. Under certain assumptions, this method can be used to define and simplify the structural dynamics, stiffness and quality of substructures.
一旦完成结构功能的定义,就可以利用它来进行结构设计、分析 和制定验证实验方案。在开发复杂结构和进行替代结构设计是是非常 有用的。  Once the definition of the structural function is completed, it can be used to perform structural design, analysis, and validation of the experimental protocol. It is very useful to develop complex structures and make alternative structural designs.
申请号 200810102136.8的发明名称为 "基于变形修正的动力减縮 算法及超单元构造方法" 中国专利申请,公开了一种基于变形修正的 动力减縮算法及超单元构造方法,在此全文引用。虽然该方法能够很 好地对复杂结构的模型进行简化,从而使得运算量大大降低。但是对 于结构件来讲,它所包含的参数还是相对较多 ,这样的话 ,对于刻画 结构功能还是比较繁琐。 The invention name of the application No. 200810102136.8 is "power reduction algorithm based on deformation correction and super unit construction method". Chinese patent application discloses a power reduction algorithm based on deformation correction and a super unit construction method, which is hereby incorporated by reference in its entirety. Although this method can simplify the model of complex structure well, the computational complexity is greatly reduced. But right In terms of structural parts, it contains relatively many parameters, so that it is rather cumbersome to characterize the structure.
基于上述缺陷,本发明提出一种基于超单元构造方法的结构功能 分析方法及其采用该结构功能分析方法的装置。 发明内容 针对上述现有技术中的问题,本发明的目的在于基于超单元构造 方法,提供一种结构功能分析方法 ,简化参数。  Based on the above drawbacks, the present invention proposes a structural function analysis method based on a super-cell construction method and an apparatus using the same. SUMMARY OF THE INVENTION In view of the above problems in the prior art, an object of the present invention is to provide a structural function analysis method based on a super-cell construction method, which simplifies parameters.
为了实现上述目的 ,本发明提供了以下技术方案: 一种基于超单元构造方法的结构功能分析方法 ,具体为 : 步骤一 根据超单元构造方法将结构单元生成简化的超单元模型 , 得到逼近原始有限元模型的位移的一个优化的逼近 = ,其中 g是 超单元的广义位移向量, T是优化的位移模式矩阵,通过变换可以得 到对应于 g的刚度、 质量和阻尼矩阵;  In order to achieve the above object, the present invention provides the following technical solutions: A structural function analysis method based on a super-cell construction method, specifically: Step 1: According to a super-cell construction method, a simplified super-cell model is generated by a structural unit, and an approximating original limited is obtained. An optimized approximation of the displacement of the metamodel = , where g is the generalized displacement vector of the superelement, T is the optimized displacement mode matrix, and the stiffness, mass and damping matrix corresponding to g can be obtained by transformation;
步骤二 根据结构件边界的位移特征对边界节点位移做进一步简 化,即:  Step 2 According to the displacement characteristics of the boundary of the structural member, the displacement of the boundary node is further simplified, namely:
q = Bq ,其中 ,所述 g为超单元广义位移向量, 7j、 β分别为根据 构件变形特征确定的广义位移向量和构件变形特征矩阵;利用 β就能 够类似超单元的方法构造结构件刚度、质量和阻尼矩阵,将它们作为 结构件的力学功能或特征参数。  q = Bq , where g is a super-unit generalized displacement vector, 7j and β are respectively a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member; using β to construct a structural member stiffness similar to the super-cell method, Mass and damping matrices, which are used as mechanical functions or characteristic parameters of structural members.
优选地,当变形矩阵能够充分逼近真实位移时, S为可直接选用 构件的变形特征矩阵;当需要进一步提高精度时,可以类似与构造 Τ 的方法通过优化来得到 。  Preferably, when the deformation matrix can sufficiently approximate the true displacement, S is a deformation feature matrix that can directly select the component; when further precision is required, it can be obtained by optimization similar to the method of constructing Τ.
根据原模型的几何形状,把结构划分为若干部件,每一部件作为 一个超单元,再根据各部件的运动特征划分为若干运动同步性区域 , 保留部件边界节点并且在部件各同步性区域特征点增设超节点。然后 采用超单元构造方法将结构单元生成简化的超单元模型。 本发明基于超单元构造方法利用变形特征矩阵对系统进行进一 步的简化。通过这种方法对结构单元的计算处理,能够使得结构特征 的参数得到大大简化,并且能够更好地保持计算的精度。 最简单的情况是直接利用前述动态子结构方法,将子结构的质量 矩阵 、 子结构的阻尼矩阵 ('、 子结构的刚度矩阵 A '、直接作为子结 构的功能参数; According to the geometry of the original model, the structure is divided into several parts, each part is regarded as a super unit, and then divided into several motion synchronization areas according to the motion characteristics of each part, the part boundary nodes are reserved and the feature points of each part of the synchronization area are Add a super node. The super unit construction method is then used to generate a simplified super unit model for the structural unit. The present invention further simplifies the system based on the super-element construction method using the deformation feature matrix. The calculation processing of the structural unit by this method can greatly simplify the parameters of the structural features and better maintain the accuracy of the calculation. The simplest case is to directly use the aforementioned dynamic substructure method, the mass matrix of the substructure, the damping matrix of the substructure ( ', the stiffness matrix A ' of the substructure, directly as the functional parameters of the substructure;
( I )在子结构较复杂的情况下 ,可以将子结构分为若干个更小 的子结构 ,利用前述方法得到各个小子结构的结构功能参数,通过拼 接 ,即类似有限元进行总体矩阵集成的方法,得到该子结构的功能参 数;  (I) In the case that the substructure is more complicated, the substructure can be divided into several smaller substructures. The structural function parameters of each substructure are obtained by the above method, and the overall matrix integration is performed by splicing, that is, similar to finite element. Method, obtaining functional parameters of the substructure;
( 2 )对几何和载荷特征比较突出的子结构 ,可以对外部载荷和 节点位移进行进一步模型化,减少自由度。即对外部边界位移进行进 一步的同步限制 ,用边界位移同步性模型进行边界位移模型化,得到 广义位移和对应的广义力 ,如梁端部的中性面位移、端面转角和与之 对应的力与力矩。 基于上述方法,本发明还提供了一种采用上述方法的装置,所述 基于超单元构造方法的结构功能分析装置包括以下结构单元: 超单元模型形成单元,根据超单元构造方法将结构单元生成简化 的超单元模型,得到逼近原始有限元模型的位移的一个优化的逼近 (2) For substructures with prominent geometric and load characteristics, external loads and nodal displacements can be further modeled to reduce degrees of freedom. That is, the external boundary displacement is further limited, and the boundary displacement synchronization model is used to model the boundary displacement, and the generalized displacement and the corresponding generalized force are obtained, such as the neutral surface displacement, the end surface rotation angle and the corresponding force of the beam end. With torque. Based on the above method, the present invention also provides an apparatus using the above method, the structural function analysis apparatus based on the super unit construction method comprising the following structural unit: a super unit model forming unit, which simplifies the structural unit generation according to the super unit construction method Super element model, an optimized approximation of the displacement approximating the original finite element model
Ή, ,其中 是超单元的广义位移向量, τ是优化的位移模式矩阵, 通过变换可以得到对应于 的刚度、 质量和阻尼矩阵; 节点位移简化单元,根据结构件边界的位移特征对边界节点位移 做进一步简化,即 '. q ,其中 ,所述 ^为超单元广义位移向量 、 ^分别为根据构件变形特征确定的广义位移向量和构件变形特征矩 阵;利用 就能够类似超单元的方法构造结构件刚度、 质量和阻尼矩 阵 ,将它们作为结构件的力学功能或特征参数。 Ή , , where is the generalized displacement vector of the superelement, τ is the optimal displacement mode matrix, and the corresponding stiffness, mass and damping matrix can be obtained by transformation; the node displacement simplification unit, the displacement of the boundary node according to the displacement feature of the structural member boundary Further simplification, that is, '.q, where ^ is a super-unit generalized displacement vector, ^ is a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member; and the structural member can be constructed by using a super-cell-like method Stiffness, mass and damping matrix, which are used as mechanical functions or characteristic parameters of structural members.
优选的,当变形矩阵能够充分逼近真实位移时, s为直接选用构 件的变形特征矩阵; 当需要进一步提高精度时,能够类似与构造 τ 的方法通过优化来得到 。 Preferably, when the deformation matrix can sufficiently approximate the true displacement, s is a deformation feature matrix of the component directly selected; when it is required to further improve the precision, the structure and the tau can be similarly constructed. The method is obtained by optimization.
本发明基于超单元构造方法利用变形特征矩阵对系统进行进一 步的简化。通过这种方法对结构单元的计算处理,能够使得结构特征 的参数得到大大简化 ,并且能够更好地保持计算的精度。  The present invention further simplifies the system based on the superelement construction method using the deformation feature matrix. The calculation processing of the structural unit by this method can greatly simplify the parameters of the structural features and better maintain the accuracy of the calculation.
最简单的情况是直接利用前述动态子结构方法,将子结构的质量 矩阵 、 子结构的阻尼矩阵 ('、 子结构的刚度矩阵 A '、直接作为子结 构的功能参数; The simplest case is to directly use the aforementioned dynamic substructure method, the mass matrix of the substructure, the damping matrix of the substructure ( ', the stiffness matrix A ' of the substructure, directly as the functional parameters of the substructure;
( 1 )在子结构较复杂的情况下 ,可以将子结构分为若干个更小 的子结构,利用前述方法得到各个小子结构的结构功能参数,通过拼 接 ,即类似有限元进行总体矩阵集成的方法,得到该子结构的功能参 数;  (1) In the case that the substructure is more complicated, the substructure can be divided into several smaller substructures, and the structural function parameters of each substructure are obtained by the foregoing method, and the overall matrix integration is performed by splicing, that is, similar finite element. Method, obtaining functional parameters of the substructure;
( 2 )对几何和载荷特征比较突出的子结构 ,可以对外部载荷和 节点位移进行进一步模型化,减少自由度。即对外部边界位移进行进 一步的同步限制,用边界位移同步性模型进行边界位移模型化,得到 广义位移和对应的广义力 ,如梁端部的中性面位移、端面转角和与之 对应的力与力矩。  (2) For substructures with prominent geometric and load characteristics, external loads and nodal displacements can be further modeled to reduce degrees of freedom. That is, the external boundary displacement is further limited, and the boundary displacement synchronization model is used to model the boundary displacement, and the generalized displacement and the corresponding generalized force are obtained, such as the neutral surface displacement, the end surface rotation angle and the corresponding force of the beam end. With torque.
有益效果  Beneficial effect
本发明将利用一种新的动态模型简化技术,基于变形修正的动力 减縮方法,给出一种动态子结构方法,进而实现结构功能参数的刻画。  The invention will utilize a new dynamic model simplification technique, based on the dynamic reduction method of deformation correction, and give a dynamic substructure method, thereby realizing the characterization of structural function parameters.
附图说明 图 1为经过超单元分析方法得到的子结构质量凝聚模型; 图 2为图 1的实施例的频率误差曲线; BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a substructure mass condensation model obtained by a super unit analysis method; FIG. 2 is a frequency error curve of the embodiment of FIG. 1;
图 3为通过超单元构造方法得到的超单元;  Figure 3 is a super unit obtained by a super unit construction method;
图 4为由超单元组成的整体结构;  Figure 4 is an overall structure composed of super units;
图 5为通过简化得到的结构单元; 图 6为由结构单元组成的整体结构。 Figure 5 is a structural unit obtained by simplification; Figure 6 is an overall structure composed of structural units.
本发明的较隹实肺式 本发明的一种结构功能分析方法 ,具体为 : The more compact lung type of the present invention is a structural function analysis method of the present invention, specifically:
1 )根据超单元构造方法将结构单元生成简化的超单元模型 ,得 到逼近原始有限元模型的位移的一个优化的逼近 = ,其中 g是超 单元的广义位移向量, T是优化的位移模式矩阵,通过变换可以得到 对应于 的刚度、 质量和阻尼矩阵;  1) According to the super-cell construction method, the simplified super-cell model is generated by the structural unit, and an optimized approximation of the displacement of the original finite element model is obtained, where g is the generalized displacement vector of the super element, and T is the optimized displacement mode matrix. The corresponding stiffness, mass and damping matrix can be obtained by transformation;
2 )根据结构件边界的位移特征对边界节点位移做进一步简化 , 即:  2) The displacement of the boundary node is further simplified according to the displacement characteristics of the boundary of the structural member, namely:
q = Bq ,其中 ,所述 g为超单元广义位移向量, 7Ϊ、 β分别为根据 构件变形特征确定的广义位移向量和构件变形特征矩阵;利用 β就可 以类似超单元的方法构造结构件刚度、质量和阻尼矩阵,将它们作为 结构件的力学功能或特征参数。  q = Bq , where g is a super-unit generalized displacement vector, 7 Ϊ, β are respectively a generalized displacement vector and a component deformation characteristic matrix determined according to the deformation characteristics of the member; using β to construct a structural member stiffness similar to the super-cell method, Mass and damping matrices, which are used as mechanical functions or characteristic parameters of structural members.
当变形矩阵能够充分逼近真实位移时, 为可直接选用构件的变 形特征矩阵;当需要进一步提高精度时,可以类似与构造 Τ的方法通 过优化来得到 Β。  When the deformation matrix can fully approximate the real displacement, the deformed feature matrix of the component can be directly selected; when the precision needs to be further improved, the method can be similarly optimized by the method of constructing Τ.
所述基于超单元构造方法的结构功能分析装置包括以下结构单 元:  The structural function analysis device based on the super unit construction method includes the following structural unit:
超单元模型形成单元,根据超单元构造方法将结构单元生成简化 的超单元模型,得到逼近原始有限元模型的位移的一个优化的逼近 Έ, ,其中 是超单元的广义位移向量, Τ是优化的位移模式矩阵, 通过变换可以得到对应于 的刚度、 质量和阻尼矩阵; The super-cell model forming unit generates a simplified super-cell model from the structural unit according to the super-cell construction method, and obtains an optimized approximation which approximates the displacement of the original finite element model, where is the generalized displacement vector of the super-unit, and Τ is optimized Displacement mode matrix, by which the corresponding stiffness, mass and damping matrix can be obtained;
节点位移简化单元,根据结构件边界的位移特征对边界节点位移 做进一步简化,即 C ,其中 ,所述 为超单元广义位移向量 J、 s分别为根据构件变形特征确定的广义位移向量和构件变形特征矩 阵;利用 S就能够类似超单元的方法构造结构件刚度、 质量和阻尼矩 阵 ,将它们作为结构件的力学功能或特征参数。 The node displacement simplification unit further simplifies the boundary node displacement according to the displacement feature of the structural member boundary, that is, C, wherein the super-unit generalized displacement vector J and s are respectively the generalized displacement vector and the component deformation determined according to the deformation characteristics of the member. Characteristic moment Array; the structural stiffness, mass and damping matrix of structural members can be constructed by using S-like super-cell method, and they are used as mechanical functions or characteristic parameters of structural members.
优选的,当变形矩阵能够充分逼近真实位移时, ^为直接选用构 件的变形特征矩阵; 当需要进一步提高精度时,能够类似与构造 T 的方法通过优化来得到 。  Preferably, when the deformation matrix can sufficiently approximate the true displacement, ^ is the deformation feature matrix of the component directly selected; when the precision needs to be further improved, it can be obtained by optimization similar to the method of constructing T.
下面通过具体实例来进一步说明。  The following is further illustrated by specific examples.
平面悬臂梁  Plane cantilever beam
各个动态子结构的质量都集中在内部超节点处,模型示意图如图 The quality of each dynamic substructure is concentrated in the internal super node.
1所示。 频率求解结果如表 I所示 ,频率误差如图 2所示。 在前面模型的基础上,类似前面的动态子结构实例,将单个子结 构上下边的节点也作为内部节点处理,就可以得到如图 3所示的子结 构 ,拼接后就得到如图 2所示是整体结构。 1 is shown. The frequency solution results are shown in Table I, and the frequency error is shown in Figure 2. On the basis of the previous model, similar to the previous dynamic substructure example, the nodes above and below the single substructure are also treated as internal nodes, and the substructure shown in Fig. 3 can be obtained. After splicing, it is obtained as shown in Fig. 2. It is the overall structure.
进一步,假设子结构两端符合(梁的截面 )平面假设(直法线 ) , 将图 2 所示子结构两边的位移用两边中点的平动位移和转角 (u, w, Θ) ( a l )
Figure imgf000008_0001
Further, assuming that the two ends of the substructure conform to the (beam cross section) plane hypothesis (straight normal line), the displacement of both sides of the substructure shown in Fig. 2 is the translational displacement and rotation angle of the two midpoints (u, w, Θ) (al )
Figure imgf000008_0001
(每一边都用上式近似) 其中 y,是节点 i现对于边界中点的平行于 边界方向的坐标。 利用 ( a2 )就可以构造图 1所示子结构两边边界节 点位移向量和两边边界中 和转角的关系 (Each side is approximated by the above formula) where y is the coordinate of the node i now parallel to the boundary direction of the midpoint of the boundary. Using ( a2 ), we can construct the relationship between the displacement vector of the boundary points on both sides of the substructure shown in Fig. 1 and the sum of the corners of the two sides.
Figure imgf000008_0002
采用专利 2008 1 0102 1 36.8类似公式 ( 23 )至 ij ( 26 ) ,就可以给 出简化后的子结构系统 (图 5简化子结构,图 6为对应整体结构 )。 表 I为和 ANSYS结果比较,说明基本一致。 表 1 图 6实施例的频率结果与 ANSYS求解的结果对比
Figure imgf000008_0002
Using the patents 2008 1 0102 1 36.8 similar to equations (23) to ij (26), a simplified substructure system can be given (Fig. 5 simplified substructure, Fig. 6 corresponds to the overall structure). Table I compares with ANSYS results, and the explanation is basically the same. Table 1 Comparison of the frequency results of the embodiment of Figure 6 with the results of ANSYS solution
 Order
ANSYS求解结果 (Hz) 平面假设简化频率结果 相 对 误 差 次 (Hz) (%)  ANSYS Solution Results (Hz) Plane Simplified Frequency Results Relative Errors Times (Hz) (%)
1 4.7858 4.7992 0.28%1 4.7858 4.7992 0.28%
2 28.792 28.802 0.03%2 28.792 28.802 0.03%
3 74.608 74.667 0.08%3 74.608 74.667 0.08%
4 76.135 75.826 -0.41%4 76.135 75.826 -0.41%
5 139 1 37.9 -0.79%5 139 1 37.9 -0.79%
6 213.02 2 1 1 .04 -0.93%6 213.02 2 1 1 .04 -0.93%
7 223.79 224.58 0.35%7 223.79 224.58 0.35%
8 294.7 1 291.55 -1.07%8 294.7 1 291.55 -1.07%
9 372.85 374.93 0.56%9 372.85 374.93 0.56%
10 381.78 377.86 -1.03%10 381.78 377.86 -1.03%
1 1 472.68 468.27 -0.93%1 1 472.68 468.27 -0.93%
12 521.61 532.03 2.00%12 521.61 532.03 2.00%
13 566.33 554.32 -2.12%13 566.33 554.32 -2.12%
14 661.9 647.85 -2.12%14 661.9 647.85 -2.12%
15 669.75 690.05 3.03% 15 669.75 690.05 3.03%
最后应说明的是:显然,上述实施例仅仅是为清楚地说明本申请 所作的举例,而并非对实施方式的限定。对于所厲领域的普通技术人 员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。 这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而 易见的变化或变动仍处于本申请型的保护范围之中。 工业实用性 本发明基于超单元构造方法利用变形特征矩阵对系统进行进一 步的简化。 能够较好地解决传统模型简化技术中的各种限制和不便, 在一定的假设下 ,通过这种方法对结构单元的计算处理,就可以对子 结构的结构动力学功能,刚度和质量,进行定义和大大简化,并且能 够更好地保持计算的精度,同时又能够保持较高的效率。 It should be noted that the above-described embodiments are merely illustrative of the present invention and are not intended to limit the embodiments. Other variations or modifications of the various forms may be made by those of ordinary skill in the art. There is no need and no way to exhaust all of the implementations. Obvious changes or variations resulting therefrom are still within the scope of this application. Industrial Applicability The present invention further simplifies the system based on the super-element construction method using the deformation feature matrix. It can better solve various limitations and inconveniences in the traditional model simplification technology. Under certain assumptions, the structural dynamics function, stiffness and quality of the substructure can be performed by calculating the structural unit by this method. It is defined and greatly simplified, and it is better able to maintain the accuracy of the calculation while maintaining high efficiency.

Claims

权 利 要 求 书 Claim
1、 一种基于超单元构造方法的结构功能分析方法,其特征在 于:步骤如下 1. A structural function analysis method based on a super-cell construction method, characterized in that: the steps are as follows
步骤一 根据超单元构造方法将结构单元生成简化的超单元模型 , 得到逼近原始有限元模型的位移的一个优化的逼近 = ,其中 g是 超单元的广义位移向量, T是优化的位移模式矩阵,通过变换可以得 到对应于 g的刚度、 质量和阻尼矩阵;  Step 1 generates a simplified super element model based on the super element construction method to obtain a simplified approximation of the displacement of the original finite element model, where g is the generalized displacement vector of the super element and T is the optimized displacement mode matrix. The stiffness, mass and damping matrix corresponding to g can be obtained by transformation;
步骤二 根据结构件边界的位移特征对边界节点位移做进一步简 化,即: q - B f ,其中 ,所述 g为超单元广义位移向量 , Ίί、 β分别为 根据构件变形特征确定的广义位移向量和构件变形特征矩阵;利用 Β 就能够类似超单元的方法构造结构件刚度、质量和阻尼矩阵,将它们 作为结构件的力学功能或特征参数。  Step 2 further simplifies the displacement of the boundary node according to the displacement characteristics of the boundary of the structural member, namely: q - B f , where g is a super-unit generalized displacement vector, and Ίί, β are respectively generalized displacement vectors determined according to the deformation characteristics of the member And the deformation characteristic matrix of the member; the stiffness, mass and damping matrix of the structural member can be constructed by using the method of super-element, which is used as the mechanical function or characteristic parameter of the structural member.
2、 根据权利要求 1 所述的方法 ,其特征在于:当变形矩阵能 够充分逼近真实位移时, 为直接选用构件的变形特征矩阵;当需要 进一步提高精度时,能够类似与构造 Τ的方法通过优化来得到 S。 2. The method according to claim 1, wherein when the deformation matrix can sufficiently approximate the true displacement, the deformation characteristic matrix of the component is directly selected; when the precision needs to be further improved, the method capable of similarly constructing the flaw can be optimized. Come to get S.
3、 基于超单元构造方法的结构功能分析装置,其特征在于, 包括: 超单元模型形成单元 ,根据超单元构造方法将结构单元生成简 化的超单元模型,得到逼近原始有限元模型的位移的一个优化的逼近 u = Tq ,其中 g是超单元的广义位移向量, T是优化的位移模式矩阵, 通过变换可以得到对应于 g的刚度、 质量和阻尼矩阵; 节点位移简化单元 ,根据结构件边界的位移特征对边界节点位 移做进一步简化,即 -. q - Bci ,其中, 超单元广义位移向量 , q β分别为根据构件变形特征确定的广义位移向量和构件变形特征 矩阵;利用 β就能够类似超单元的方法构造结构件刚度、 质量和阻尼 矩阵 ,将它们作为结构件的力学功能或特征参数。 3. A structural function analysis device based on a super-cell construction method, comprising: a super-cell model forming unit, which generates a simplified super-cell model according to a super-cell construction method, and obtains a displacement approximating the displacement of the original finite element model The approximate approximation u = Tq , where g is the generalized displacement vector of the superelement, T is the optimal displacement mode matrix, and the stiffness, mass and damping matrix corresponding to g can be obtained by transformation; the node displacement simplification unit, according to the boundary of the structural member The displacement feature further simplifies the displacement of the boundary node, namely -. q - Bci , where the super-unit generalized displacement vector, q β is the generalized displacement vector and the component deformation characteristic matrix determined according to the deformation characteristics of the member respectively; Unit method to construct structural members for stiffness, mass and damping Matrices, which are used as mechanical functions or characteristic parameters of structural members.
4、 根据权利要求 3所述的装置 ,其特征在于:当变形矩阵能 够充分逼近真实位移时, 为直接选用构件的变形特征矩阵;当需要 进一步提高精度时,能够类似与构造 Τ的方法通过优化来得到 S。 4. The device according to claim 3, wherein when the deformation matrix can sufficiently approximate the true displacement, the deformation characteristic matrix of the component is directly selected; when the precision needs to be further improved, the method can be optimized similarly to the construction flaw. Come to get S.
PCT/CN2014/079933 2013-06-18 2014-06-16 Super-unit construction method based structure-function analysis method and device therefor WO2014201984A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201310242125.0 2013-06-18
CN201310242125.0A CN104239588B (en) 2013-06-18 2013-06-18 A kind of structural-functional analysis method based on ultra-unit construction method

Publications (1)

Publication Number Publication Date
WO2014201984A1 true WO2014201984A1 (en) 2014-12-24

Family

ID=52103954

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2014/079933 WO2014201984A1 (en) 2013-06-18 2014-06-16 Super-unit construction method based structure-function analysis method and device therefor

Country Status (2)

Country Link
CN (1) CN104239588B (en)
WO (1) WO2014201984A1 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105512361A (en) * 2015-11-22 2016-04-20 沈阳黎明航空发动机(集团)有限责任公司 Substructure analysis method based on ANSYS software
CN106855896A (en) * 2016-11-04 2017-06-16 南方科技大学 Method for designing based on hyperelement and spherical geometry structure
CN109918712B (en) * 2019-01-23 2022-07-05 昆明理工大学 Calculation method for solving composite stress intensity factor based on p-type finite element method and contour integral method
CN112580239B (en) * 2020-12-15 2022-07-05 中南大学 Substructure response reconstruction method and system, and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006171867A (en) * 2004-12-13 2006-06-29 Canon Inc Numerical value calculation method
CN101334803A (en) * 2008-03-18 2008-12-31 中国科学院力学研究所 Dynamic cut down algorithm based on distortion correction and ultra-unit construction method
CN102855363A (en) * 2012-09-26 2013-01-02 中国人民解放军装甲兵工程学院 Method for simulating and calculating dynamic stress of a vehicle body structure of high-speed tracked vehicle
US20130151217A1 (en) * 2009-03-24 2013-06-13 Landmark Graphics Corporation Systems and Methods for Modeling Drillstring Trajectories

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5594651A (en) * 1995-02-14 1997-01-14 St. Ville; James A. Method and apparatus for manufacturing objects having optimized response characteristics
CN100587696C (en) * 2008-03-19 2010-02-03 中国科学院力学研究所 Parameterization static state ultra-unit structuring method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006171867A (en) * 2004-12-13 2006-06-29 Canon Inc Numerical value calculation method
CN101334803A (en) * 2008-03-18 2008-12-31 中国科学院力学研究所 Dynamic cut down algorithm based on distortion correction and ultra-unit construction method
US20130151217A1 (en) * 2009-03-24 2013-06-13 Landmark Graphics Corporation Systems and Methods for Modeling Drillstring Trajectories
CN102855363A (en) * 2012-09-26 2013-01-02 中国人民解放军装甲兵工程学院 Method for simulating and calculating dynamic stress of a vehicle body structure of high-speed tracked vehicle

Also Published As

Publication number Publication date
CN104239588A (en) 2014-12-24
CN104239588B (en) 2018-02-16

Similar Documents

Publication Publication Date Title
Le et al. A gradient-based, parameter-free approach to shape optimization
Huang et al. Interactive shape interpolation through controllable dynamic deformation
Du et al. A convergent adaptive finite element algorithm for nonlocal diffusion and peridynamic models
Yamasaki et al. A structural optimization method based on the level set method using a new geometry‐based re‐initialization scheme
da Silva et al. Reliability-based topology optimization of continuum structures subject to local stress constraints
WO2014201984A1 (en) Super-unit construction method based structure-function analysis method and device therefor
Pietroni et al. Statics aware grid shells
Ghorbel et al. Dynamic stiffness formulation for free orthotropic plates
Sen et al. A fast and robust hybrid method for block‐structured mesh deformation with emphasis on FSI‐LES applications
Fan et al. A novel numerical manifold method with derivative degrees of freedom and without linear dependence
Narvydas et al. Applications of sub-modeling in structural mechanics
López et al. Simultaneous untangling and smoothing of moving grids
Vaxman Modeling polyhedral meshes with affine maps
de Frias et al. A multiscale mass scaling approach for explicit time integration using proper orthogonal decomposition
Mundis et al. An efficient flexible GMRES solver for the fully-coupled time-spectral aeroelastic system
Lu et al. Optimal placement of FBG sensors for reconstruction of flexible plate structures using modal approach
Amoozandeh Nobaveh et al. Symmetric kinetostatic behavior from asymmetric spatially curved beams
Shimoda Free-form optimization method for designing automotive shell structures
JP2007328596A (en) Method for vibration analysis of construct model having symmetry
CN105701287A (en) Method for designing three-way equal-rigidity rubber shock absorber of platform type inertial navigation system
Wang et al. Quasi-conforming triangular Reissner-Mindlin shell elements by using Timoshenko's Beam Function
WO2013078912A1 (en) Numerical method for simulating wake field of ship propeller
Ilanko et al. Bipenalty method from a frequency domain perspective
Shimoda et al. Multi-objective free-form optimization of shell structures
Benner et al. Model order reduction of large-scale dynamical systems with Jacobi-Davidson style eigensolvers

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: 14813372

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: 14813372

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 31.05.2016)

122 Ep: pct application non-entry in european phase

Ref document number: 14813372

Country of ref document: EP

Kind code of ref document: A1