CN113076677A - Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic Hermite shape function - Google Patents
Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic Hermite shape function Download PDFInfo
- Publication number
- CN113076677A CN113076677A CN202110406021.3A CN202110406021A CN113076677A CN 113076677 A CN113076677 A CN 113076677A CN 202110406021 A CN202110406021 A CN 202110406021A CN 113076677 A CN113076677 A CN 113076677A
- Authority
- CN
- China
- Prior art keywords
- flexible body
- body structure
- structure unit
- equation
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Complex Calculations (AREA)
Abstract
The invention discloses a novel method for numerical simulation of a high-order nonlinear finite element of a flexible body structure based on a quintic Hermite shape function, which is used for deducing the quintic Hermite interpolation shape function and an attribute equation based on a quintic Hermite high-order interpolation theory; determining a flexible body structure unit and parameters, a physical model and simplifying assumptions which need to be solved by the flexible body structure unit; establishing a mathematical motion partial differential equation of the flexible body structure unit in space and time and a stretching limiting condition along the axial direction, namely solving a mathematical model in a domain; thereby establishing a system equation of the flexible body structure unit and adopting a Newton iteration method to carry out iteration solution on the parameters. Has the advantages that: the simulation precision of the space deformation, the corner and the curvature of the slender flexible body structure is improved; the continuity of the curvature at the finite element node is ensured; the fast convergence is realized, and particularly, the prediction precision of the high-order vibration frequency of the flexible body is higher; when the problems of transverse large gradient load and large bending of the flexible body are solved, higher precision guarantee is provided.
Description
Technical Field
The invention relates to the technical field of flexible body (slider structures) simulation, in particular to a high-order nonlinear finite element numerical simulation method of a flexible body structure based on a quintic Hermite shape function.
Background
Up to now, the most commonly used finite element method for flexible body structural elements (such as beam elements) in the world is still a finite element method based on cubic Hermite shape functions (cubic Hermite shape functions). The method can only ensure the continuity of the space configuration (spatial configuration) of the flexible body and the first derivative (namely the corner slope angle) of the configuration in the discrete nodes of the finite element.
Since the spatial form of the flexible body unit structure is simulated and segmented into a cubic polynomial function, the first derivative is a spatial quadratic function, namely, quadratic distribution, and the second derivative is a spatial linear function, namely, linear distribution. Although this method has achieved good numerical simulation accuracy, it has some inherent drawbacks: the simulated result of the method can only ensure the continuity of the first derivative function of the space form of the space flexible structure at the finite element node, but cannot ensure the continuity of the second derivative function; and the Dirichlet (Dirichlet) boundary condition cannot be applied for the flexible body cell boundary condition with zero curvature (curvature), resulting in deviation (deviation) of the physical model and the numerical model boundary condition; when the flexible body is subjected to a large gradient load (namely, when the distribution of bending moment and shearing force inside the structure body has large mutation), the algorithm is difficult to accurately capture the positions and the sizes of extreme values of the bending moment and the shearing force (position and mass of the critical bending and shear force), and a very fine finite element unit is required to capture the extreme values.
When the flexible body structure is designed, the natural vibration frequency of the flexible body often needs to be accurately known, so that the natural frequency range of main external excitation load possibly encountered by the structure in the use process is avoided, and the negative influence caused by resonance is prevented or reduced. By adopting a cubic Hermite finite element method, the natural frequency of high-order vibration of the flexible body structure is difficult to be measured accurately, a large error exists, and the larger the unit is, the larger the error is.
Disclosure of Invention
Aiming at the problems of the traditional method, the invention provides a high-order nonlinear finite element numerical simulation method of a flexible body structure based on a quintic Hermite shape function. The quintic Hermite special-shaped function finite element method is introduced into the computer finite element numerical simulation of the flexible body structure, and the computer numerical simulation precision of the space deformation, the corner and the curvature of the slender flexible body structure is improved.
In order to achieve the purpose, the invention adopts the following specific technical scheme:
a flexible body structure high-order nonlinear finite element numerical simulation method based on a quintic Hermite shape function is characterized by comprising the following specific steps:
the method comprises the following steps: based on a Quintic Hermite (Quintic Hermite) high-order interpolation theory, a Quintic Hermite interpolation shape function is derived under a one-dimensional coordinate system, and the Quintic Hermite interpolation shape function comprises 6 Quintic functions:
wherein the content of the first and second substances,xaand xbRespectively as the start point coordinate and the end point coordinate of the space finite element unit node, and the coordinate interval range under the one-dimensional overall coordinate system x is x ∈ [ ]a,xb](ii) a The coordinate interval range of xi under the local coordinate system is xi epsilon [0,1];
The attribute equation of the quintic Hermite interpolation shape function is as follows:
step two: determining a flexible body structure unit and parameters to be solved by the flexible body structure unit, determining a physical model of the flexible body structure unit according to size parameters of the flexible body structure unit, and establishing a mathematical motion partial differential equation of the flexible body structure unit in space and time on the basis of adopting a reasonable simplifying assumption; analyzing to obtain the deformation characteristic of the flexible body structure unit in the axial direction, and acquiring the stretching limiting condition of the flexible body structure unit along the axial direction;
the mathematical motion partial differential equation of the flexible body structure unit in space and time is concretely as follows:
wherein M ismIs a quality matrix;
r、q is a variable which changes with the arc length s and the time t under an arc coordinate system;
wherein r is a space instantaneous position vector of the flexible body structure unit after deformation;
b is the bending rigidity of the flexible body structure unit, B is EI, E is the Young modulus of the material of the flexible body structure unit, and I is the moment of inertia of the flexible body structure unit to the neutral axis;
q is the load distribution of the external transverse load along the length direction of the flexible body structure unit;
the black point at the top of the symbol represents the derivative with respect to time, and the two points represent the second derivative with respect to time; the "'" in the upper right hand corner of the symbol indicates the first derivative to spatial arc length; "" in the upper right hand corner of the symbol indicates the second derivative to spatial arc length; the "'" in the upper right hand corner of the symbol indicates the third derivative to the spatial arc length; the "" in the upper right hand corner of the symbol indicates the fourth derivative to spatial arc length; by default, bold black indicates that the physical quantity is a vector or tensor, and non-bold indicates a scalar quantity.
Step three: substituting the quintic Hermite interpolation shape function (1) obtained in the first step into a mathematical motion partial differential equation (3) of the flexible body structure unit in space and time, integrating the whole flexible body structure unit by adopting partial integration, and establishing a nonlinear system equation on the space unit of the flexible body structure unit by combining the quintic Hermite interpolation shape function attribute equation (2) in the first step:
wherein L represents the length of the flexible body structural unit;
ai(s) represents the expression of the quintic hermitian interpolation function in the overall arc coordinate system;
ai(s)=Ji·φi(xi) (i ═ 1 to 6); establishing a conversion relation of the overall coordinate system and the form function expressions under the unit coordinate system; j. the design is a squareiRepresenting Jacobian transformation coefficients under different coordinate systems;
the relationship between subscript i and length L of the flexible body structural unit is as follows:
step four: applying boundary conditions at two ends of the flexible body structure unit, and defining the equal sign right side of equation (5) as the generalized force F of the flexible body structure unitiGeneralized force F of the flexible body structural unitiThe expression is as follows:
step five: discretizing the parameters in the nonlinear system equation (5) in the third step by using a tensor algorithm:
obtaining the spatial instantaneous position distribution r (s, t) and the high-order derivative r of the flexible body structural unit by the quintic Hermite interpolation shape function approximate simulation(n)(s,t);
By a quadratic interpolation function pm(s) approximate simulation to obtain effective tension distribution inside the flexible body structural unitMass distribution M (s, t) inside the flexible body structure unit, and bending rigidity distribution B(s) inside the flexible body structure unit;
the specific simulation physical quantity equation is as follows:
step six: substituting the numerical simulation physical quantity equation (7) obtained in the step five into a nonlinear system equation (5) on a space unit of the flexible body structure unit, and combining the generalized force F of the flexible body structure unit in the step fiveiAnd (3) deriving a discrete equation in the structural unit of the flexible body by using an expression (6) and an attribute equation (2) of a quintic Hermite interpolation shape function:
αikman integral coefficient matrix related to a bending stiffness term generated by the Galerkin method;
βikman integral coefficient matrix related to a tension term generated by the Galerkin method;
γikman integral coefficient matrix related to an inertial force term generated by the Galerkin method;
qinfor the direction of external transversely-distributed loads acting on the flexible body unitAn amount;
Finis the generalized force vector at the unit node;
step seven: performing coordinate system transformation on discrete equations (8) within the flexible body structural units: the coordinate system of the quintestic Hermite interpolation shape function is converted from a space overall coordinate system through a Jacobian conversion coefficient JiAnd (3) converting the data into a unit local coordinate system to obtain a discrete equation in the flexible body structure unit under the local coordinate system, so that the rigidity coefficient matrix in the equation (8) can be conveniently solved.
Wherein, the discrete system equation (9) in the flexible body structure unit under the local coordinate system is:
all stiffness coefficient matrices in equation (8) solve the relationship:
step eight: establishing a stretching limiting condition of the flexible body structure unit along the axial direction in a simultaneous step II, repeating the steps III to V based on a similar discrete method, discretizing the stretching limiting condition of the flexible body structure unit along the axial direction to obtain a discrete equation (11) corresponding to the stretching limiting condition of the flexible body structure unit along the axial direction;
step nine: and solving the nonlinear coupling system equation on the space unit of the flexible body structure unit by a Newton iteration method according to the discrete system equation (9) in the flexible body structure unit under the local coordinate system and the discrete equation (11) corresponding to the stretching limit condition of the equation flexible body structure unit along the axial direction.
By the design, a finite element algorithm of a quintic Hermite special function is adopted,
the curvature is brought into generalized displacement to be directly solved, so that the continuity of the curvature at a finite element node can be guaranteed, and the continuity of the bending moment of the flexible body and the bending moment induced stress in a space solution domain is guaranteed; the algorithm can solve a smaller total stiffness matrix under the condition of dividing fewer units, and compared with the traditional finite element method, the method is higher in precision and low in calculation cost; under the condition that the same flexible body structure unit is divided into the same thickness, the precision of the new algorithm for solving the natural frequency of the flexible body is higher, particularly the precision of the high-order natural frequency is higher, so that the more accurate dynamic characteristic of the flexible body can be obtained, and the dynamic and fatigue prediction precision is improved. For solving the problem of the flexible body with large curvature, the precision of the new algorithm can be well guaranteed no matter the thickness of the divided units, and the traditional method needs extremely fine units to obtain higher precision, so that the numerical finite element of the flexible body is solved, the reliability of the numerical simulation result is better than that of the traditional finite element method no matter the thickness of the units of the calculation result of the new method, and the confidence of a decision maker on the calculation result and engineering design can be increased.
The flexible body structure unit has a high utilization ratio in engineering and shows a trend of a higher utilization ratio (because the engineering structure meets the strength condition, the use of construction materials can be effectively reduced by adopting a large number of flexible body slender structures, the structure is light and convenient to install, and the total construction cost can be greatly reduced), and the flexible body structure unit is particularly suitable for various ocean engineering projects. Such as various marine risers, pipelines, anchor chains, umbilical cables, pipe-in-pipe, multi-layer pipes, novel composite pipes, and the like; cables, high-voltage wires, optical cables, ropes of various materials; engineering beam structures of various shapes, such as bearing beams, upright beams, L-shaped/T-shaped/U-shaped/rectangular beams for buildings, long-span bridges, wings of airplanes, fan blades for wind power generation, vertical supporting uprights and the like, belong to flexible body structures.
Further, in the second step, the deformation characteristics of the flexible body structure unit in the axial direction at least include a case that the flexible body structure unit is not stretchable along the axial direction and a case that the flexible body structure unit is deformed and stretched along the axial direction;
for the case that the flexible body structure unit is not stretchable along the axial direction, the axial stretching limiting conditions are as follows:
r′·r′=1 (4-a)
for the flexible body structure unit along the axial deformation stretching condition, the axial stretching limiting condition is as follows:
r′·r′=(1+ε)2 (4-b)
in equations (4-a) and (4-b), the middle black dot represents the inner product or dot product of the vectors;
epsilon represents the axial strain of the flexible body structural unit,wherein T is the local axial tangential tension, A represents the cross-sectional area of the flexible body structure unit, and E is the Young modulus of the flexible body structure unit.
The limit of the stretch of the flexible body in the axial direction generally includes the different situations of non-stretch in the axial direction, small deformation stretch, large deformation stretch, etc. In the present application, a large deformation tensile flexible body is not an object of study, and examples thereof include a rope having a large axial strain.
In a further technical solution, in step eight, for the case that the flexible body structural unit is not stretchable along the axial direction, the discrete equation corresponding to the stretching limitation condition of the axial direction is:
for the case that the flexible body structure unit has deformation and stretching along the axial direction, the discrete equation corresponding to the stretching limiting condition of the axial direction is as follows:
for example, the flexible body structure is a marine riser, and the relational expression of the axial strain of the flexible body structure unit, the axial effective tension and the internal and external hydrostatic pressures of the riser is as follows:
Afthe cross section area is corresponding to the outer diameter of the flexible body structure unit;
Acthe cross section area is corresponding to the inner diameter of the flexible body structure unit;
ρfgyAf-ρcgyActhe pressure difference of the static water external pressure and the static water internal pressure at the position of a vertical coordinate y is corresponding to the external fluid and the internal fluid of the flexible body structure unit;
when the flexible body structure unit is not in the fluid, i.e. the environment has a small density such as air, ρ can be approximately consideredfgyAf-ρcgyAc=0;
The coefficient matrix in the equation can be solved by the following integral relation:
further, when the flexible body structure unit is in a static state, and a newton iteration method is used for solving the nonlinear system equation on the space unit of the flexible body structure unit, the first term, i.e., the unsteady term, of the equation (9) can be ignored, and the problem can be simplified to a greater extent. If the letter r is used to represent the number of iteration steps, the equations of the iterative algorithm for the r-th step and the r-1-th step are detailed in equations (14) and (15):
in a further technical scheme, when the flexible body structure unit is in a dynamic state, a Newton iteration method is adopted to solve a nonlinear system equation on a space unit of the flexible body structure unit, and a letter r is used for representing iteration steps. The non-constant term of equation (9) is not negligible, and the total time to be simulated in the time domain needs to be divided into a plurality of time steps s, so that the system updates the displacement, velocity and acceleration of the flexible body after balancing at each time step, and performs iteration at the next time step by using the updated information. Is a process that moves forward with time steps and iterates continuously. For the current time step, an iterative relationship algorithm can be established to solve the axial small-deformation flexible body motion equation and solve the power problem, and the size of the time step is generally determined according to the convergence condition. The iteration is carried out at each time step (indicated by letter s), and the iterative algorithm formula of the r step and the r-1 step is detailed in equations (16) and (17):
compared with the traditional cubic Hermite special function finite element method, the method has the beneficial effects that:
the curvature of the flexible body is brought into generalized displacement by adopting a new quintic hermite interpolation shape function finite element algorithm (1) to be directly solved, and the continuity of the curvature at a finite element node can be ensured, so that the continuity of the bending moment of the flexible body and the continuity of bending moment induced stress in a space solution domain are ensured
(2) The new algorithm can solve a smaller total stiffness matrix under the condition of dividing fewer units, and can obtain higher precision compared with the traditional cubic Hermite shape function finite element method, so that the computer computing cost can be saved.
(3) Under the condition of the same unit thickness, the precision of the new algorithm for solving the natural frequency of the flexible body, particularly the precision of the high-order natural frequency, can be higher, so that the more accurate dynamic characteristic of the flexible body can be obtained, and the dynamic strength calculation and fatigue life prediction precision of the flexible body are improved.
(4) For solving the problem of the flexible body with large curvature, the precision of the new algorithm can be well guaranteed no matter the thickness of the divided units, while the traditional method (such as a common cubic Hermite shape function finite element method) needs extremely fine units to obtain higher precision, so that the finite element numerical value of the flexible body is solved, the calculated result of the new method is better than that of the traditional finite element method no matter the thickness of the units, and the confidence of a decision maker on the structure calculated result and the engineering design can be increased.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is an image of a quintic Hermite shape function;
FIG. 3 is a flow chart of a flexible body structure element dynamics state solution for a quintic hermitian function finite element method;
FIG. 4 is a schematic view of a simple beam model according to example 1;
FIG. 5 is a schematic diagram of the spatial configuration of the riser of example 3 before deformation;
FIG. 6 is a schematic view of the deformed riser in the embodiment 3.
Detailed Description
The following provides a more detailed description of the embodiments and the operation of the present invention with reference to the accompanying drawings.
A method for simulating a high-order nonlinear finite element numerical value of a flexible body structure based on a quintic Hermite shape function can be seen by referring to FIG. 1, and comprises the following specific steps:
the method comprises the following steps: based on a Quintic Hermite (Quintic Hermite) high-order interpolation theory, a Quintic Hermite interpolation shape function is deduced under a one-dimensional coordinate system, and a shape function image is shown in detail in figure 2; the quintic hermite interpolation shape function comprises 6 quintic shape functions, specifically:
wherein the content of the first and second substances,xaand xbRespectively as the start point coordinate and the end point coordinate of the space finite element unit node, and the coordinate interval range under the one-dimensional overall coordinate system x is x ∈ [ ]a,xb](ii) a The coordinate interval range of xi under the local coordinate system is xi epsilon [0,1];
The attribute equation of the quintic Hermite interpolation shape function is as follows:
step two: determining a flexible body structure unit and parameters to be solved by the flexible body structure unit, determining a physical model of the flexible body structure unit according to size parameters of the flexible body structure unit, and establishing a mathematical motion partial differential equation of the flexible body structure unit in space and time on the basis of adopting a reasonable model hypothesis; analyzing to obtain the deformation characteristic of the flexible body structure unit in the axial direction, and acquiring the stretching limiting condition of the flexible body structure unit along the axial direction;
the mathematical motion partial differential equation of the flexible body structure unit in space and time is concretely as follows:
wherein M ismIs a quality matrix;
r、q is a variable which changes with the arc length s and the time t under an arc coordinate system;
wherein r is a space instantaneous position vector of the flexible body structure unit after deformation;
b is the bending rigidity of the flexible body structure unit, B is EI, E is the Young modulus of the material of the flexible body structure unit, and I is the moment of inertia of the flexible body structure unit to the neutral axis;
q is the load distribution of the external transverse load along the length direction of the flexible body structure unit;
the black point at the top of the symbol represents the derivative with respect to time, and the two points represent the second derivative with respect to time; the "'" in the upper right hand corner of the symbol indicates the first derivative to spatial arc length; "" in the upper right hand corner of the symbol indicates the second derivative to spatial arc length; the "'" in the upper right hand corner of the symbol indicates the third derivative to the spatial arc length; the "" in the upper right hand corner of the symbol indicates the fourth derivative to spatial arc length; by default, bold black indicates that the physical quantity is a vector or tensor, and non-bold indicates a scalar quantity.
In the second step, the deformation characteristics of the flexible body structure unit in the axial direction at least comprise the condition that the flexible body structure unit is not stretched along the axial direction and the condition that the flexible body structure unit is deformed and stretched along the axial direction;
for the case that the flexible body structure unit is not stretchable along the axial direction, the axial stretching limiting conditions are as follows:
r′·r′=1 (4-a)
for the flexible body structure unit along the axial deformation stretching condition, the axial stretching limiting condition is as follows:
r′·r′=(1+ε)2 (4-b)
in equations (4-a) and (4-b), the middle black dot represents the inner product or dot product of the vectors;
epsilon represents the axial strain of the flexible body structural unit,wherein T is the local axial tangential tension, A represents the cross-sectional area of the flexible body structure unit, and E is the Young modulus of the flexible body structure unit.
Step three: substituting the quintic Hermite interpolation shape function (1) obtained in the first step into a mathematical motion partial differential equation (3) of the flexible body structure unit in space and time, integrating the whole flexible body structure unit by adopting partial integration, and establishing a nonlinear system equation on the space unit of the flexible body structure unit by combining the quintic Hermite interpolation shape function attribute equation (2) in the first step:
wherein L represents the length of the flexible body structural unit;
ai(s) represents the expression of the quintic hermitian interpolation function in the overall arc coordinate system;
ai(s)=Ji·φi(xi) (i ═ 1 to 6); establishing a conversion relation of the overall coordinate system and the form function expressions under the unit coordinate system;
Jirepresenting Jacobian transformation coefficients under different coordinate systems;
the relationship between subscript i and length L of the flexible body structural unit is as follows:
step four: applying boundary conditions at two ends of the flexible body structure unit, and defining the equal sign right side of equation (5) as the generalized force F of the flexible body structure unitiGeneralized force F of the flexible body structural unitiThe expression is as follows:
step five: discretizing the parameters in the nonlinear system equation (5) in the third step by using a tensor algorithm:
obtaining the spatial instantaneous position distribution r (s, t) and the high-order derivative r of the flexible body structural unit by the quintic Hermite interpolation shape function approximate simulation(n)(s,t);
By a quadratic interpolation function pm(s) approximate simulation to obtain effective tension distribution inside the flexible body structural unitMass distribution M (s, t) inside the flexible body structure unit, and bending rigidity distribution B(s) inside the flexible body structure unit;
the specific simulation physical quantity equation is as follows:
step six: substituting the numerical simulation physical quantity equation (7) obtained in the step five into a nonlinear system equation (5) on a space unit of the flexible body structure unit, and combining the generalized force F of the flexible body structure unit in the step fiveiAnd (3) deriving a discrete equation in the structural unit of the flexible body by using an expression (6) and an attribute equation (2) of a quintic Hermite interpolation shape function:
αikman integral coefficient matrix related to a bending stiffness term generated by the Galerkin method;
βikman integral coefficient matrix related to a tension term generated by the Galerkin method;
γikman integral coefficient matrix related to an inertial force term generated by the Galerkin method;
qinvectors generated for external laterally distributed loads applied to the flexible body elements;
Finis the generalized force vector at the unit node;
step seven: performing coordinate system transformation on discrete equations (8) within the flexible body structural units: the coordinate system of the quintestic Hermite interpolation shape function is converted from a space overall coordinate system through a Jacobian conversion coefficient JiConverting the data into a unit local coordinate system to obtain a discrete equation in the flexible body structure unit under the local coordinate system, and conveniently solving a rigidity coefficient matrix in the equation (8);
wherein, the discrete system equation (9) in the flexible body structure unit under the local coordinate system is:
all stiffness coefficient matrices in equation (8) solve the relationship:
step eight: establishing a stretching limiting condition of the flexible body structure unit along the axial direction in a simultaneous step II, repeating the steps III to V based on a similar discrete method, discretizing the stretching limiting condition of the flexible body structure unit along the axial direction to obtain a discrete equation (11) corresponding to the stretching limiting condition of the flexible body structure unit along the axial direction;
in step eight, for the case that the flexible body structure unit is not stretchable along the axial direction, the discrete equation corresponding to the stretching limitation condition of the axial direction is as follows:
for the case of the flexible body structure unit deforming and stretching along the axial direction, the stretching limit condition of the axial direction corresponds to a discrete equation:
wherein, the relational expression of flexible body structure unit axial strain and effective pulling force of axial and riser inside and outside hydrostatic pressure does:
Afthe cross section area is corresponding to the outer diameter of the flexible body structure unit; for flexible body structural units as risers AfThe cross section area is a circular cross section area corresponding to the outer diameter of the vertical pipe structural unit;
Acthe cross section area is corresponding to the inner diameter of the flexible body structure unit; for flexible body structural units as risers AcIs a circular cross-sectional area corresponding to the inner diameter of the stand pipe.
ρfgyAf-ρcgyAcIs a corresponding flexible bodyThe pressure difference between the static water external pressure and the internal pressure of the external fluid and the internal fluid of the structural unit at the vertical coordinate y;
when the flexible body structure unit is not in the fluid, the correction term rhofgyAf-ρcgyAc=0;
The coefficient matrix in equation (11) can be solved by the following integral relation:
step nine: and solving the nonlinear coupling system equation on the space unit of the flexible body structure unit by a Newton iteration method according to the discrete system equation (9) in the flexible body structure unit under the local coordinate system and the discrete equation (11) corresponding to the stretching limit condition of the equation flexible body structure unit along the axial direction.
When the flexible body structure unit is in a static state, a Newton iteration method is adopted to solve a nonlinear system equation on a space unit of the flexible body structure unit, the number of iteration steps is represented by a letter r, and the iterative algorithm formula of the step r and the step r-1 is detailed in equations (14) and (15):
and when the flexible body structure unit is in a dynamic state, solving a nonlinear system equation on the space unit of the flexible body structure unit by adopting a Newton iteration method, and expressing iteration steps by using a letter r. The non-constant term of equation (9) is not negligible, and the total time to be simulated in the time domain needs to be divided into a plurality of time steps s, so that the system updates the displacement, velocity and acceleration of the flexible body after balancing at each time step, and performs iteration at the next time step by using the updated information. Is a process that moves forward with time steps and iterates continuously. For the current time step, an iterative relationship algorithm can be established to solve the axial small-deformation flexible body motion equation and solve the power problem, and the size of the time step is generally determined according to the convergence condition. The iteration is carried out at each time step (indicated by letter s), and the iterative algorithm formula of the r step and the r-1 step is detailed in equations (16) and (17):
specifically, referring to fig. 3, when the flexible body structure unit is in a dynamic state, the step of solving the nonlinear system equation on the space unit of the flexible body structure unit by using a newton iteration method includes:
s1: collecting all input data required by finite element numerical simulation calculation of the flexible body structure unit;
s2: initializing displacement, speed, acceleration, effective tension and all convergence control parameters of the flexible body structure unit, determining time step length and total time step, and setting the initialization time step s to be 0
S3: the letter r is used to represent the number of iteration steps, the initialization iteration step r is 0, and s is s +1
S4: r is r + 1; establishing a finite element unit stiffness array of the flexible body structure unit;
s5: assembling all the unit stiffness arrays into a total stiffness array;
s6: applying boundary conditions and loads and numerically solving a total stiffness array;
s7: updating the displacement, speed and acceleration of the current iteration step;
s8: judging whether the displacement increment is smaller than the preset convergence precision or not; if yes, converging, and outputting a convergence result of the step s; if not, go to step S9;
s9: judging whether the current r is smaller than the maximum allowable iteration number; if so, the process returns to step S4, otherwise, the process diverges and returns to step S2.
The algorithm provided by the invention has a wide application range, wherein the flexible structure has a high utilization ratio in engineering and shows a trend of higher utilization ratio (because the engineering structure meets the strength condition, the use of construction materials can be effectively reduced by adopting a large number of flexible body slender structures, and the structure is light and convenient to install, so that the total construction cost can be greatly reduced), and particularly in ocean engineering projects. Such as various marine risers, pipelines, anchor chains, umbilical cables, pipe-in-pipe, multi-layer pipes, novel composite pipes, and the like; cables, high-voltage wires, optical cables, signal transmission cables, ropes and the like made of various materials; various engineering beam structures, such as bearing beams, upright beams, L-shaped/T-shaped/U-shaped/rectangular beams for buildings, long-span bridges, wings of airplanes, fan blades and uprights for wind power generation and the like, all belong to slender flexible structures. By applying the novel high-order interpolation algorithm of the patent, the computer numerical simulation precision of various elongated flexible bodies can be improved, and the continuity of the curvature and bending induced stress distribution in the structure in the problem solving domain is guaranteed. Meanwhile, due to the rapid convergence of the new algorithm, the calculation result is not sensitive to the division thickness of the finite element unit, so that a thicker grid can be divided under the condition of ensuring the convergence of the algorithm, and the result and the precision of the result meeting the engineering requirement can be obtained under the condition of reducing the total freedom degree of the result solution.
The invention is described below by way of three examples:
the first embodiment is as follows: as can be seen from fig. 4, the flexible body structure unit is: a simply supported beam which is elastically supported; and solving the position and the value of the maximum bending moment of the simply supported beam with the elastic support under the action of seven distributed loads.
The parameters of the simply supported beam of the elastic support are as follows: l8 m, EI 6.9e3Nm2Q0 is 50N/m, and the rigidity of the uniformly distributed elastic support is Ks 1.0e 4N/m/m;
based on the above steps, the results of calculating the bending moment by the quintic hermite finite element new method (QH) are shown in table 1, and the results of calculating the bending moment by the conventional cubic hermite finite element method (CH) are also shown in table 1:
table 1: comparison of results of calculation of bending moment in working example 1
In the table, the Mesh Size table divides the cell thickness, and the cell length is 1m as shown in 1; MBM represents the value of the maximum bending moment; MBM-best represents the most accurate maximum bending moment value which can be obtained; pos represents the location where the maximum bending moment occurs on the beam; rdif represents the relative error;
as can be seen from Table 1, the calculation result of the traditional cubic Hermite finite element method (CH) is sensitive to the unit division thickness, and the position of the maximum bending moment is not very accurate; the maximum bending moment value obtained by the quintic Hermite finite element new method (QH) is insensitive to the thickness of the unit division, and the position of the maximum bending moment is very accurate under the condition of the coarsest unit division. Similar conclusions can be applied to the calculation of the maximum shear force in table 2 below.
The results of calculating the maximum shear using the conventional cubic hermite finite element method (CH) and the quintic hermite finite element new method (QH) proposed by the present invention are shown in table 2:
table 2: comparison of maximum shear calculation results of example 1
In the table: the Mesh Size table divides the thickness of the unit, and the unit length is 1m as 1; MSF denotes the value of maximum shear; MSF-best represents the most accurate maximum shear value that can be obtained; rdif represents the relative error; gdops represents the total number of degrees of freedom required for solution.
Example two: a journal paper was published in 1976 by Dareing and Huang for natural frequency calculations for a drilling riser. The specific literature names are: dareing, d.w., Huang, t.,1976.Natural frequency of marine driling rise.j.pet.tech 28, 813-818.
The flexible body structure unit is: a drilling riser; the flexible body structural unit size parameters are detailed in table 3;
table 3: natural frequency calculation input parameters for drilling risers
The results of the calculations to obtain the natural frequency of the drilling riser according to the above calculation steps are detailed in table 4:
table 4: comparison of results of natural frequency calculation in example 2
Considering that the natural frequency of the drilling RISER has no mathematical analytic solution, if the result of the quintic Hermite finite element new method (RISER3D _ QH) proposed by the present invention in the case of 14 element divisions is taken as the reference standard result, the results of the old finite element method (cubic Hermite FEM) in the case of 7 and 14 element divisions and the new algorithm (quintic Hermite FEM) in the case of 7 element divisions can be compared with the reference result, and the calculated relative error is shown in table 4. As can be seen from Table 4, after the new method (RISER3D _ QH) refines the unit, the natural frequency of the drilling RISER has little variation with respect to the error, and the percentage variation of the natural frequency of the highest order 12 is not more than 1%. For the traditional finite element algorithm (CABLE3D _ CH), when the unit is changed from thick to thin, the change of the natural frequency of 12 th order is as high as about 30 percent, which shows that the result of the traditional finite element algorithm is very sensitive to the thickness division of the unit and the convergence speed is not good; the quintic Hermite finite element new method (QH) provided by the invention has considerable precision no matter the thickness of the unit, and is insensitive to the thickness division of the unit.
Example three: predicting the bending moment of a shallow water Jumper after a large transverse translation is carried out on a top platform;
this arithmetic example refers to the international journal paper published in Connaire, et al (2015), which is collectively referred to as: advances in the subsea roller analysis using quasi-ratios and the Newton-Raphson method; from the documents Advances in the subsea roller analysis using ration-ratios and the Newton-Raphson method, it can be seen that the main input parameters are shown in tables 5 and 6:
TABLE 5 input parameters for simulating deformation of flexible pipe when large transverse translation of top platform occurs
TABLE 6 hypothetical parameters for simulating deformation of a flexible pipe during large lateral translations of the top end platform
TABLE 7 comparison of bending moment calculations
Referring to fig. 5, the present example mainly simulates the large bending deformation of a flexible riser after the top end of the flexible riser (i.e., the origin of coordinates of fig. 5) moves horizontally to the right 15 meters along the x direction, and compares how the maximum bending moment inside the riser changes when different unit thicknesses are divided by using the conventional finite element algorithm (CH) and the quintessential finite element new method (QH) proposed by the present invention.
Table 7 shows the maximum bending moment comparison of the two algorithms, the conventional cubic hermite finite element method (CH) and the quinhermite finite element method (QH), for different element thicknesses when the riser is moved to the position of FIG. 6.
As can be seen from the results shown in table 7, if the calculation result of the quintic hermitian finite element new method (QH) in the case of the 0.125-meter cell division is taken as the benchmarking result, the result accuracy in the case of the other cell division can be analyzed. Wherein BMQH-Dif represents the quintic Hermite finite element New method (QH) relative error; BMCH-Dif represents the relative error of the conventional cubic Hermite finite element method (CH); it can be seen that the result of the quintic hermite finite element new method (QH) does not change much regardless of the unit thickness, whereas the conventional cubic hermite finite element method (CH) necessarily requires extremely fine meshing to obtain a more accurate result, and the coarser the mesh, the larger the error. BM _ FDM in table 7 indicates that the calculation is analyzed by using the conventional finite element method (CH), but the post-processing uses a 5-point difference method to obtain information of the second derivative, i.e. curvature, by using the obtained information difference, and then calculates the bending moment value. With the above information, the Quintic Hermite shape function (finite element) finite element method adopted by the patent is mainly applicable to the following working conditions or examples:
and (3) accurately predicting the numerical finite element of the maximum bending moment/shearing force of various flexible bodies under the action of any load.
The accurate numerical finite element accurate prediction of the high-order natural vibration frequency of various flexible bodies (including flexible bodies which cannot resist bending, such as anchor chains and ropes with extremely small bending rigidity) is carried out.
When the change gradient of the transverse load born by the flexible body is large, the prediction of the internal force of the flexible body can be obviously improved.
The internal force prediction with higher precision is carried out on the condition that a flexible body with rigidity bears larger bending (the problem of large curvature);
based on the algorithm, the internal force calculated by the flexible body static force and dynamic force problem is more accurate, and the method can be concluded that the method is not only suitable for the flexible body strength problem, but also has remarkable help for the prediction accuracy of the flexible body fatigue problem.
It should be noted that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make variations, modifications, additions or substitutions within the spirit and scope of the present invention.
Claims (5)
1. A flexible body structure high-order nonlinear finite element numerical simulation method based on a quintic Hermite shape function is characterized by comprising the following specific steps:
the method comprises the following steps: based on a Quintic Hermite (Quintic Hermite) high-order interpolation theory, a Quintic Hermite interpolation shape function is deduced, and the Quintic Hermite interpolation shape function comprises 6 Quintic functions:
wherein the content of the first and second substances,xaand xbRespectively as the start point coordinate and the end point coordinate of the space finite element node, and the coordinate interval range under the overall coordinate system x is x ∈ [ x [ ]a,xb](ii) a The coordinate interval range of xi under the local coordinate system is xi epsilon [0,1];
The attribute equation of the quintic Hermite interpolation shape function is as follows:
step two: determining a flexible body structure unit and parameters to be solved by the flexible body structure unit, determining a physical model of the flexible body structure unit according to size parameters of the flexible body structure unit, and establishing a mathematical motion partial differential equation of the flexible body structure unit in space and time; analyzing to obtain the deformation characteristic of the flexible body structure unit in the axial direction, and acquiring the stretching limiting condition of the flexible body structure unit along the axial direction;
the mathematical motion partial differential equation of the flexible body structure unit in space and time is concretely as follows:
wherein M ismIs a quality matrix;
r、q is a variable which changes with the arc length s and the time t under an arc coordinate system;
wherein r is a space instantaneous position vector of the flexible body structure unit after deformation;
b is the bending rigidity of the flexible body structure unit, B is EI, E is the Young modulus of the material of the flexible body structure unit, and I is the moment of inertia of the flexible body structure unit to the neutral axis;
q is the load distribution of the external transverse load along the length direction of the flexible body structure unit;
the black point at the top of the symbol represents the derivative with respect to time, and the two points represent the second derivative with respect to time; the "'" in the upper right hand corner of the symbol indicates the first derivative to spatial arc length; the symbol "", in the upper right hand corner, represents the second derivative to spatial arc length; the ""' in the upper right hand corner of the symbol indicates the third derivative to the spatial arc length; the ""' in the upper right hand corner of the symbol indicates the fourth derivative to the spatial arc length;
step three: substituting the quintic Hermite interpolation shape function (1) obtained in the first step into a mathematical motion partial differential equation (3) of the flexible body structure unit in space and time, integrating the whole flexible body structure unit by adopting partial integration, and establishing a nonlinear system equation on the space unit of the flexible body structure unit by combining the quintic Hermite interpolation shape function attribute equation (2) in the first step:
wherein L represents the length of the flexible body structural unit;
ai(s) represents the expression of the quintic hermitian interpolation function in the overall arc coordinate system;
ai(s)=Ji·φi(xi) (i is 1-6), and establishing a conversion relation of the overall coordinate system and the unit coordinate system lower form function expressions;
jirepresenting Jacobian transformation coefficients under different coordinate systems;
the relationship between subscript i and length L of the flexible body structural unit is as follows:
step four: applying boundary conditions at two ends of the flexible body structure unit, and defining the equal sign right side of equation (5) as the generalized force F of the flexible body structure unitiGeneralized force F of the flexible body structural unitiThe expression is as follows:
step five: discretizing the parameters in the nonlinear system equation (5) in the third step by using a tensor algorithm:
obtaining the spatial instantaneous position distribution r (s, t) and the high-order derivative r of the flexible body structural unit through quintic Hermite interpolation shape function simulation(n)(s,t);
By a quadratic interpolation function pm(s) simulation to obtain effective tension distribution inside the structure unit of the flexible bodyMass distribution M (s, t) inside the flexible body structure unit, and bending rigidity distribution B(s) inside the flexible body structure unit;
the specific simulation physical quantity equation is as follows:
step six: substituting the numerical simulation physical quantity equation (7) obtained in the step five into a flexible body knotCombining the nonlinear system equation (5) on the space unit of the structural unit with the generalized force F of the flexible body structural unit in the step fiveiAnd (3) deriving a discrete equation in the structural unit of the flexible body by using an expression (6) and an attribute equation (2) of a quintic Hermite interpolation shape function:
αikman integral coefficient matrix related to a bending stiffness term generated by the Galerkin method;
βikman integral coefficient matrix related to a tension term generated by the Galerkin method;
γikman integral coefficient matrix related to an inertial force term generated by the Galerkin method;
qinvectors generated for external laterally distributed loads applied to the flexible body elements;
Finis the generalized force vector at the unit node;
step seven: performing coordinate system transformation on discrete equations (8) within the flexible body structural units: the coordinate system of the quintestic Hermite interpolation shape function is converted from a space overall coordinate system through a Jacobian conversion coefficient JiConverting the data into a unit local coordinate system to obtain a discrete equation in the flexible body structure unit under the local coordinate system, and solving a rigidity coefficient matrix in an equation (8);
wherein, the discrete system equation (9) in the flexible body structure unit under the local coordinate system is:
all stiffness coefficient matrices in equation (8) solve the relationship:
step eight: establishing a stretching limiting condition of the flexible body structure unit along the axial direction in a simultaneous step II, repeating the step III to the step V, discretizing the stretching limiting condition of the flexible body structure unit along the axial direction to obtain a discrete equation (11) corresponding to the stretching limiting condition of the flexible body structure unit along the axial direction;
step nine: and solving a nonlinear system coupling equation on the space unit of the flexible body structure unit by a Newton iteration method according to a discrete system equation (9) in the flexible body structure unit under the local coordinate system and a discrete equation (11) corresponding to the stretching limit condition of the equation flexible body structure unit along the axial direction.
2. The fifth order hermitian function based high order nonlinear finite element numerical simulation method of a flexible body structure according to claim 1, wherein in the second step, the deformation characteristics of the flexible body structure unit in the axial direction at least include a condition that the flexible body structure unit is not stretchable in the axial direction and a condition that the flexible body structure unit is stretched in the axial direction;
for the case that the flexible body structure unit is not stretchable along the axial direction, the axial stretching limiting conditions are as follows:
r′·r′=1 (4-a)
for the flexible body structure unit along the axial deformation stretching condition, the axial stretching limiting condition is as follows:
r′·r′=(1+ε)2 (4-b)
in equations (4-a) and (4-b), the middle black dot represents the inner product or dot product of the vectors;
3. The fifth order hermitian function-based numerical simulation of a flexible body structure high order nonlinear finite element according to claim 2, wherein in the eighth step, for the case that the flexible body structure unit is not stretchable along the axial direction, the discrete equation corresponding to the stretching constraint condition along the axial direction is:
for the case that the flexible body structure unit has deformation and stretching along the axial direction, the discrete equation corresponding to the stretching limiting condition of the axial direction is as follows:
wherein, the relational expression of flexible body structure unit axial strain and effective pulling force of axial and riser inside and outside hydrostatic pressure does:
Afthe cross section area is corresponding to the outer diameter of the flexible body structure unit;
Acthe cross section area is corresponding to the inner diameter of the flexible body structure unit;
ρfgyAf-ρcgyActhe pressure difference of the static water external pressure and the static water internal pressure at the position of a vertical coordinate y is corresponding to the external fluid and the internal fluid of the flexible body structure unit;
when the flexible body structure unit is not in the fluid, then rhofgyAf-ρcgyAc=0;
The coefficient matrix in the equation can be solved by the following integral relation:
4. the fifth order hermitian function based flexible body structure high order nonlinear finite element numerical simulation method of claim 3, wherein:
when the flexible body structure unit is in a static state, a Newton iteration method is adopted to solve a nonlinear system equation on a space unit of the flexible body structure unit, the number of iteration steps is represented by a letter r, and the iterative algorithm formula of the step r and the step r-1 is detailed in formulas (14) and (15):
5. the fifth order hermitian function based flexible body structure high order nonlinear finite element numerical simulation method of claim 3, wherein:
when the flexible body structure unit is in a dynamic state, performing numerical integration on time by adopting a Newmark-beta numerical calculation method, performing time step division on finite element simulation time, and expressing the time step by using s; for the current s time step, when a Newton iteration method is adopted to solve the nonlinear system equation on the space unit of the flexible body structure unit, the letter r represents the number of iteration steps, and the iterative algorithm formula of the r step and the r-1 step is detailed in formulas (16) and (17):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406021.3A CN113076677B (en) | 2021-04-15 | 2021-04-15 | Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic hermite function |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406021.3A CN113076677B (en) | 2021-04-15 | 2021-04-15 | Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic hermite function |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113076677A true CN113076677A (en) | 2021-07-06 |
CN113076677B CN113076677B (en) | 2023-03-17 |
Family
ID=76617606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110406021.3A Active CN113076677B (en) | 2021-04-15 | 2021-04-15 | Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic hermite function |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113076677B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361183A (en) * | 2021-07-15 | 2021-09-07 | 中国第一汽车股份有限公司 | Method for forecasting ultimate bearing capacity of automobile half shaft |
CN115146377A (en) * | 2022-06-23 | 2022-10-04 | 西北工业大学 | Coordinated design method for structural rigidity of aircraft with layout of connecting wings |
CN117216911A (en) * | 2023-11-07 | 2023-12-12 | 天津大学 | Single-column offshore wind turbine structure response calculation method based on inertia release theory |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008234039A (en) * | 2007-03-16 | 2008-10-02 | Furukawa Electric Co Ltd:The | Simulation system and simulation method for long flexible structure, program and information recording medium |
US20100169060A1 (en) * | 2008-12-31 | 2010-07-01 | Zhenhai Zhu | Method, system, and computer program prodcut for lithography simulation in electronic design automation |
US20190392099A1 (en) * | 2018-06-25 | 2019-12-26 | Worcester Polytechnic Institute | Method and apparatus for modeling electromagnetic fields using hermite finite elements |
CN111159636A (en) * | 2019-12-04 | 2020-05-15 | 大连理工大学 | Flexible multi-body system dynamics semi-analytic sensitivity analysis method based on absolute node coordinate description |
-
2021
- 2021-04-15 CN CN202110406021.3A patent/CN113076677B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008234039A (en) * | 2007-03-16 | 2008-10-02 | Furukawa Electric Co Ltd:The | Simulation system and simulation method for long flexible structure, program and information recording medium |
US20100169060A1 (en) * | 2008-12-31 | 2010-07-01 | Zhenhai Zhu | Method, system, and computer program prodcut for lithography simulation in electronic design automation |
US20190392099A1 (en) * | 2018-06-25 | 2019-12-26 | Worcester Polytechnic Institute | Method and apparatus for modeling electromagnetic fields using hermite finite elements |
CN111159636A (en) * | 2019-12-04 | 2020-05-15 | 大连理工大学 | Flexible multi-body system dynamics semi-analytic sensitivity analysis method based on absolute node coordinate description |
Non-Patent Citations (6)
Title |
---|
F. A. RADU等: "Hermite finite elements for convection-diffusion equations", 《ARXIV》 * |
余志刚等: "基于高阶有限元方法的裂纹斜梁振动特性分析", 《振动与冲击》 * |
刘承范等: "用配置法结合适应埃尔米特元族求解运移方程", 《世界地质》 * |
夏逸鸣等: "区间五次Hermite样条多小波Euler-Bernoulli梁单元", 《固体力学学报》 * |
朱礼云: "深水海底管道S型铺管法安装分析", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》 * |
靳红玲等: "空间柔性梁的动态可靠性分析", 《西安电子科技大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361183A (en) * | 2021-07-15 | 2021-09-07 | 中国第一汽车股份有限公司 | Method for forecasting ultimate bearing capacity of automobile half shaft |
CN113361183B (en) * | 2021-07-15 | 2023-02-21 | 中国第一汽车股份有限公司 | Method for forecasting ultimate bearing capacity of automobile half shaft |
CN115146377A (en) * | 2022-06-23 | 2022-10-04 | 西北工业大学 | Coordinated design method for structural rigidity of aircraft with layout of connecting wings |
CN115146377B (en) * | 2022-06-23 | 2023-02-24 | 西北工业大学 | Coordinated design method for structural rigidity of aircraft with layout of connecting wings |
CN117216911A (en) * | 2023-11-07 | 2023-12-12 | 天津大学 | Single-column offshore wind turbine structure response calculation method based on inertia release theory |
CN117216911B (en) * | 2023-11-07 | 2024-02-02 | 天津大学 | Single-column offshore wind turbine structure response calculation method based on inertia release theory |
Also Published As
Publication number | Publication date |
---|---|
CN113076677B (en) | 2023-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113076677B (en) | Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic hermite function | |
Kan et al. | Nonlinear dynamic and deployment analysis of clustered tensegrity structures using a positional formulation FEM | |
Li et al. | Geometrically nonlinear analysis of thin-shell structures based on an isogeometric-meshfree coupling approach | |
Tessler et al. | Inverse FEM for full-field reconstruction of elastic deformations in shear deformable plates and shells | |
Langley | Application of the dynamic stiffness method to the free and forced vibrations of aircraft panels | |
Singh et al. | Analysis of cracked plate using higher-order shear deformation theory: Asymptotic crack-tip fields and XIGA implementation | |
Baiz et al. | Linear buckling analysis of cracked plates by SFEM and XFEM | |
Heltai et al. | A natural framework for isogeometric fluid–structure interaction based on BEM–shell coupling | |
Craugh et al. | Coupled finite element and cellular automata methods for analysis of composite structures with fluid–structure interaction | |
Fernandes et al. | Multi-scale modelling for bending analysis of heterogeneous plates by coupling BEM and FEM | |
CN106777691A (en) | For the rubber o-ring finite element modeling method of Structural Dynamics emulation | |
Duarte et al. | Mesh‐independent p‐orthotropic enrichment using the generalized finite element method | |
Bel Hadj Ali et al. | On static analysis of tensile structures with sliding cables: the frictional sliding case | |
CN111651907A (en) | Modeling analysis method of complex cable network system | |
CN108491591A (en) | A kind of hot environment lower curve stiffened panel finite element method | |
Stanford et al. | Formulation of analytical design derivatives for nonlinear unsteady aeroelasticity | |
Al-Dojayli et al. | Accurate modeling of contact using cubic splines | |
Thari et al. | Adaptive reduced-order modeling for non-linear fluid–structure interaction | |
Jonsson et al. | Development of flutter constraints for high-fidelity aerostructural optimization | |
Bitaraf et al. | Large deflection analysis of flexible plates by the meshless finite point method | |
CN110147571B (en) | Topology optimization method and device of component structure | |
Binesh et al. | Elasto-plastic analysis of reinforced soils using mesh-free method | |
CN115292953A (en) | Mechanical simulation analysis method for analyzing two-dimensional periodic heterogeneous structure | |
Anderson et al. | Adjoint-Based High-Fidelity Aeroelastic Optimization of Wind Turbine Blade for Load Stress Minimization | |
Asok | Interface adapted LBB-stable finite elements on fluid structure interaction problems in fully Eulerian framework |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |