Method for implementing ultimate strength analysis of plate frame structure based on isogeometric analysis

A method for implementing ultimate strength analysis of plate frame structure based on isogeometric analysis, comprising:

• according to the plate frame structure, constructing an IGA model, wherein the IGA model comprises several plane shell elements;
• carrying out a coordinate transformation on the plane shell elements;
• carrying out a weak coupling on boundaries of transformed plane shell elements by a Nitsche method and carrying out a mesh refinement;
• setting material properties, loads and boundary conditions of refined plane shell elements, and constructing a nonlinear equation according to setting results;
• and
• solving the nonlinear equation by a Newton-Raphson iteration method, and carrying out a simulation analysis by an NX software to obtain an ultimate strength of the plate frame structure;
• wherein the coordinate transformation is carried out through four kinds of coordinate systems, comprising global Cartesian coordinate system, local Cartesian coordinate system from control points, parametric coordinate system and local coordinate system at Gaussian integration;
• a process of carrying out the simulation analysis by the NX software comprises:
• according to the plate frame structure, obtaining element information and design information of a geometric model, wherein the element information of the geometric model comprises an order of basic functions, node vectors, control point coordinates and weights;
• the design information comprises node vectors, element number information, topological relationship between elements and nodes and normal vectors of the control points;
• constructing the geometric model through the element information and the design information, and traversing the elements corresponding to all the node vectors of the geometric model circularly, and calculating an element stiffness matrix in a corresponding local Cartesian coordinate system;
• obtaining a local coordinate system of Gaussian integral points of the nodes, calculating a global elastic coefficient matrix, obtaining a Jacobi matrix corresponding to the nodes, calculating a strain displacement matrix according to the global elastic coefficient matrix and the Jacobi matrix, and calculating and obtaining an element stiffness matrix in the global Cartesian coordinate system according to the strain displacement matrix;
• connecting and assembling the element stiffness matrix in the global Cartesian coordinate system into an overall stiffness matrix;
• and
• according to a material and the nonlinear equation of the plate frame structure, constructing a nonlinear solution algorithm by using an incremental step length dynamic calculation mode, and obtaining the ultimate strength of the plate frame structure by solving the nonlinear equation through the nonlinear solution algorithm;
• the incremental step length dynamic calculation mode is:
• according to an incremental step length, carrying out an incremental step length loop, then initializing related variables and matrices in an operation calculation of the NX software, and entering an iterative loop, wherein the incremental step length is obtained by presetting or calculating;
• in the iterative loop, calculating the element stiffness matrix and an initial stress node force matrix in the global Cartesian coordinate system one by one in the geometric model, and connecting and assembling into the overall stiffness matrix and a global initial stress node force matrix respectively;
• wherein in a calculation process of the element stiffness matrix and the initial stress node force matrix, obtaining the element stiffness matrix and the initial stress node force matrix in the global Cartesian coordinate system by using reduced integral and based on the global elastic coefficient matrix and the strain displacement matrix obtained by the calculation;
• according to the overall stiffness matrix and the global initial stress node force matrix, calculating displacement increments;
• carrying out an isogeometric nonlinear iteration based on the nonlinear equation to judge whether convergence conditions are met, and if the convergence conditions are met, outputting all the displacement increments, namely the ultimate strength of the plate frame structure;
• wherein the convergence conditions are:
• the incremental step length is equal to a maximum load, an iterations is less than a maximum iterations, and an incremental step length split time is less than a maximum split time.