274 results about "Stiffness matrix" patented technology

The invention discloses a jacket platform real-time correction and state prediction method based on a reduced-order model, and relates to the technical field of ocean engineering structure health monitoring and safety early warning. Comprising the steps of establishing a jacket platform finite element model; obtaining a sample matrix; constructing a reduced order model; reducing the order of the reduced-order model until the calculation precision of the reduced-order model meets a preset precision requirement, and obtaining a global reduced-order basis matrix and the reduced-order model meeting the preset precision requirement; reconstructing a stiffness matrix and a mass matrix by using the global reduced-order basis matrix; obtaining a corrected finite element model; inverting the marine environment load based on the on-site platform structure strain response data, and deducing the marine environment load; inputting the marine environment load into a pre-trained LSTM model, and outputting time series data including future load through the LSTM model; and predicting the future operation state trend of the jacket platform. The operation safety and reliability of the jacket platform are improved.

The invention discloses a rapid calculation method based on a non-coordination and reanalysis theory and a computer device. The invention relates to the technical field of structural design optimization, and solves the problem that it is difficult to realize efficient structural response analysis while ensuring the calculation precision in the prior art. The method comprises the following steps: selecting a structural member of a to-be-calculated structural response and discretizing the structural member; preprocessing the data, and determining a coarse and fine grid index set corresponding to interface fine grid nodes; solving a unit shape function, and embedding the shape function into a body stiffness matrix according to a degree-of-freedom position to serve as an interface constraint condition; calculating a variable element stiffness matrix, and completing the construction of a modified model stiffness matrix; carrying out iterative calculation to obtain a basis vector; and calculating and solving the structural response of the modified model. According to the method, the efficient response analysis of the non-coordinated variable structure is realized on the premise of ensuring the calculation precision.

The invention discloses a nonlinear beam unit analysis method considering prestressed tendon bonding slippage, which comprises the following steps: establishing an overall structure analysis model according to the geometric dimension of a prestressed concrete beam, dividing calculation units, and assembling unit stiffness matrixes into a structural overall stiffness matrix; applying load and boundary conditions to the overall structure analysis model, and solving the overall analysis model by adopting a Newton iteration method; and extracting calculation results of structural deformation, stress, section internal force, prestressed tendon relative slippage and the like, drawing a load-deformation curve, and calculating structural rigidity and bearing capacity according to the slope and peak value of the curve. The calculation unit provided by the invention is an 8-degree-of-freedom beam unit, not only is the stress characteristic of the beam structure considered, but also the bond slip characteristic of the prestressed tendon is considered, the problem that a traditional beam unit analysis model cannot consider the bond slip effect of the prestressed tendon is solved, and compared with an entity unit analysis model, the convergence is better, and the calculation efficiency is higher.

The invention relates to the technical field of force sensors, and discloses a nonlinear decoupling method and system for a three-axis force sensor, and the method comprises the steps: collecting multi-dimensional original data, and carrying out the adaptive noise reduction, dynamic selection of a wavelet basis, and construction of a filtering matrix, and obtaining a noise reduction signal through calculation; constructing a dynamic PINN proxy model, fusing material parameters and environment parameters, outputting a stiffness matrix, establishing an online parameter updating mechanism, calculating the obtained optimization parameters, and constructing an enhanced dynamic compensation matrix; the output stiffness matrix and the enhanced dynamic compensation matrix are fused, and a main singular value is reserved; reconstructing to obtain a reconstructed strain field, establishing a multi-physics field coupling model to compensate quantum-thermal coupling errors, and compensating wire inductive coupling under high frequency; constructing a multi-scale graph structure, strengthening a loss function through physical constraint, performing crystal plasticity hysteresis compensation, fusing multi-scale features, and calculating and outputting a decoupling force value; and an intelligent calibration decision is established, the drift trend is predicted, and the calibration period is dynamically adjusted.

The invention provides a crustal stress prediction method and device, electronic equipment and a storage medium, and relates to the technical field of seismic survey. The method comprises the following steps: obtaining observation seismic data of seismic wavelets in a strong VTI medium, and constructing a to-be-inverted parameter matrix based on a PP wave reflection coefficient corresponding to the strong VTI medium; constructing a posterior probability function obeyed by an inversion parameter matrix corresponding to the to-be-inverted parameter matrix based on a Bayesian inversion theory and observation seismic data, and determining a target functional based on a prior probability function and a likelihood function corresponding to the posterior probability function; determining medium density and each stiffness matrix coefficient based on an inversion parameter matrix solving result of the target functional, and determining a flexibility matrix of the strong VTI medium based on each stiffness matrix coefficient; and predicting the ground stress distribution of the target profile in the strong VTI medium based on the medium density and the positive strain matrix and the flexibility matrix corresponding to the strong VTI. Therefore, the crustal stress prediction accuracy under the strong VTI medium is improved.

The invention discloses a method and a system for acquiring deformation stress increment of extremely-thin lamellar weak surrounding rock. Relates to the technical field of tunnel engineering. Obtaining a strain increment matrix of the extremely-thin laminar weak surrounding rock under the overall coordinate system, calculating a stress tensor matrix of deformation of the extremely-thin laminar weak surrounding rock under the overall coordinate system by combining the rotation matrix and the element stiffness matrix under an elastic assumption theory, and calculating a deformation stress increment according to the stress tensor matrix after characteristic correction; according to the scheme, the method is improved on the basis of a traditional extremely-thin lamellar surrounding rock mechanical property simulation technology, and a stress increment matrix of deformation of the extremely-thin lamellar weak surrounding rock under an overall coordinate system is calculated on the basis of a strain increment matrix, a rotation matrix and an element stiffness matrix under an elastic assumption theory; after the characteristics of the stress increment matrix are corrected, the deformation stress increment is calculated, the stress deformation characteristics of the extremely-thin lamellar weak surrounding rock are accurately and reasonably described, and the engineering analysis efficiency of the extremely-thin lamellar weak surrounding rock is improved.

The invention discloses an oral cavity finite element simulation calculation method based on collision detection. The method comprises the following steps: 1) acquiring oral scanning data of teeth, generating surface grid data of the invisible appliance according to the oral scanning data, and then generating a tetrahedral grid from the surface grid of the invisible appliance for subsequent simulation calculation; 2) on the basis of collision detection, calculating a pushing vector of teeth to the invisible appliance, and taking the pushing vector as a displacement boundary condition of subsequent calculation; 3) constructing a three-dimensional finite element model taking voxels as units, calculating an element stiffness matrix, assembling an overall stiffness matrix, loading displacement boundary conditions, forming a finite element solution equation set, and obtaining a model node displacement result; and 4) calculating stress and strain of the finite element model according to a displacement result. Aiming at the orthodontic mechanical analysis of the invisible orthodontic appliance, the invention provides a modeling method for the interaction between the teeth and the invisible orthodontic appliance based on a collision detection technology, completes the subsequent simulation solution of a finite element model, and gives consideration to the solution precision and speed.

The invention discloses a bridge anti-fatigue performance analysis and optimization method and system based on finite element analysis. The method comprises the following steps: constructing a bridge finite element basic model and a random fleet time-varying load spectrum, and executing damage-rigidity dynamic coupling iterative calculation; a fatigue sensitive area is automatically identified through the strain energy density gradient, a high-precision sub-model is generated, a basic model stiffness matrix is updated in real time through a sub-model calculation result, and a fatigue evolution cloud picture containing a stiffness degradation track is generated; and finally, identifying a damage hotspot based on the cloud picture, and generating an optimization strategy containing a traffic current limiting threshold and a structure reinforcement opportunity by using a genetic algorithm. According to the method, the problem of prediction deviation of traditional static analysis is solved by simulating the reverse influence of damage on rigidity and internal force redistribution, and accurate evaluation and scientific decision-making of the whole life cycle of the bridge are achieved.

The invention discloses a vibration solving method for an elliptical perforated plate in any shape based on an elliptical energy element, and belongs to the technical field of thin-wall structures, and the method comprises the steps: carrying out the overall mapping of the elliptical perforated plate in any shape, building an elliptical equation in a dimensionless coordinate system, and carrying out the Boolean operation modeling of geometry; constructing a boundary condition function simulation constraint; simulating board deformation according to a first-order shear deformation theory, and constructing a global trial function in combination with a Legendre polynomial and a boundary condition function; establishing a total energy functional based on a minimum potential energy principle and a Ritz method, solving a total energy functional stationary value, and obtaining element calculation formulas in the stiffness matrix and the mass matrix; and constructing an elliptical energy element based on local mapping, and performing numerical integration on each element of the matrix based on the elliptical energy element and completing vibration problem solving. According to the method, an energy functional and a numerical integration formula are both constructed on a square solution domain in an overall mapping stage, so that the calculation solution process of the elliptical perforated plate in any shape is completely standard.

The invention relates to a BIM (Building Information Modeling)-based building structure modeling method, device, equipment and medium, and the method comprises the steps: aiming at a curvature-mechanics semantic fault problem existing in a process of converting a special-shaped curved surface structure from a BIM geometric model to finite element analysis, generating a curvature tensor field with parameter coordinates by analyzing an NURBS curved surface; integrating along the maximum curvature change direction to construct an internal force flow function, and forming a guide field for controlling the grid density; shell / beam units are divided based on the guide field gradient, space coordinates are marked, and a curvature-stiffness matrix and an internal force flow driving variable cross-section are combined for calculation to generate a structural analysis model; the gradient distribution of the guide field is dynamically corrected through displacement error feedback, and lossless transmission of geometric topology and mechanical properties is achieved. According to the method, the distortion risk of special-shaped structure analysis is remarkably reduced, artificial experience correction is avoided, and high-precision automatic modeling support is provided for a complex building structure.

The invention relates to a train-rail coupling vibration analysis method and device, computer equipment, a storage medium and a computer program product. The method comprises the following steps: dividing a plurality of track units, constructing a motion differential equation of a steel rail in the track units in a frequency domain, and solving to obtain a relation between longitudinal displacement and frequency on nodes at two ends of the steel rail and a relation between stress and frequency on the nodes at the two ends of the steel rail; constructing a steel rail dynamic stiffness matrix for the steel rails in the rail units, and constructing a rail dynamic stiffness matrix for the rail units according to the steel rail dynamic stiffness matrix; assembling to obtain an integral dynamic stiffness matrix of the train-rail system, and constructing a transfer matrix of a rail; and determining a track irregularity excitation power spectrum, and solving to obtain a dynamic response index of the train-track system on each degree of freedom under the track irregularity excitation power spectrum according to the transfer matrix and the overall dynamic stiffness matrix. By adopting the method, the calculation efficiency of the dynamic response index can be improved.

The present invention relates to the technical field of the intersection of computational mechanics and deep learning, and specifically to a deep learning method that combines a physical information neural network with a finite element method; specifically, partial differential equations are embedded into a neural network, comprising the following steps: S1, meshing the solution domain, and calculating the stiffness matrix and the node load array; S2, collecting observation data and boundary conditions; S3, using the energy represented by the stiffness matrix and the fitting error of the observation data as a loss function to train the neural network; S4, executing an optimization algorithm to optimize the neural network; S5, and repeating the training and optimization steps until convergence; the deep finite element method proposed in the present invention combines a physical information neural network with a finite element method, and exhibits significant advantages and beneficial effects in terms of the fusion of physical knowledge and observation data, modeling process and loss function design, model training stability and simplicity, computational efficiency and accuracy, pre- and post-processing technology and transfer learning capability, and applicability.

The invention belongs to the technical field of robot cooperative control, and particularly relates to a distributed self-adaptive control method and system for cooperation of multiple mechanical arms, and the method comprises the steps: decoupling a pure internal force vector from the data of a six-dimensional force sensor at the tail end through pre-identified kinetic parameters; calculating a cooperative confrontation index according to the pure internal force vector and the relative position vector, calculating a self-adaptive internal force penalty weight through nonlinear mapping based on the cooperative confrontation index, establishing a state-space equation of the system based on a double-integral model, introducing an equivalent environment stiffness matrix, mapping a prediction error vector in a prediction time domain into a prediction internal force vector, and obtaining a prediction internal force vector; and a cost function containing the product of the self-adaptive internal force penalty weight and the predicted internal force vector is constructed, an optimal control increment sequence is obtained by minimizing the cost function, and a speed instruction is updated to drive the mechanical arm to move. According to the method, the robustness and the execution stability of the system are effectively improved, and the cooperation safety is guaranteed.

The invention discloses a method for predicting fatigue fracture of a high-temperature discontinuous structure based on a phase field method under a variable-temperature working condition, and belongs to the technical field of fracture behavior prediction. Acquiring a stress-strain field of a coupling temperature field of the viscous-plastic material under any cyclic load; based on a stress strain field of a coupling temperature field, a phase field residual error and a phase field stiffness matrix are obtained by obtaining a fracture phase field of a viscous-plastic material under any cyclic load, a phase field variable is obtained through iterative solution, and a corresponding UEL subprogram is compiled; and correcting and updating a stress-strain field, a temperature field and a phase field by adopting a data interaction process of a UMAT subprogram, a UMATHT subprogram and a UEL subprogram, and predicting the fracture behavior of the viscous-plastic material under the temperature gradient by identifying a phase field variable and an evolution path thereof under a cyclic load.

The invention provides a volume coefficient matrix rapid solving method for residual stress field inference and application, and relates to the technical field of material stress measurement and characterization. The method comprises the following steps: S1, calculating an inverse matrix of a non-constraint node sub-stiffness matrix of an initial state of a part; s2, calculating a non-constraint node sub-stiffness matrix variable quantity caused by part material removal in each process step; s3, calculating an inverse matrix of a non-constraint node sub-stiffness matrix of the part in each process step; s4, calculating a volume coefficient matrix of the part in each process step; and S5, assembling the volume coefficient matrixes of the parts in each process step into a volume coefficient matrix. The method solves residual stress in residual stress field inference by utilizing a volume coefficient matrix. According to the method, the calculation efficiency of the volume coefficient matrix is remarkably improved, and the problem of low calculation efficiency of the volume coefficient matrix in the prior art is effectively solved.

The invention discloses a gas turbine pull rod rotor wheel disc vibration analysis method, system and device, and relates to the technical field of gas turbines.The method comprises the steps that a three-dimensional finite element kinetic equation of a rotating part is constructed; the contact stiffness of adjacent contact units is deduced by considering the rough contact characteristic of the contact interface and combining a microscopic contact mechanical model of the contact interface and a rough consolidation interface contact mechanical model based on the statistical theory; a stress strengthening effect is considered, and a stress stiffness matrix and a contact stiffness matrix are iteratively calculated; a central pull rod end face tooth structure considering local contact nonlinearity is constructed, and an end face tooth contact model considering collision separation is constructed; carrying out vibration response solution on the three-dimensional finite element kinetic equation by adopting a degree-of-freedom reduction method; according to the method, the vibration analysis model is established, the pull rod rotor does not need to be simplified, and geometric features, contact surface features and load effects of an original model can be reserved to a great extent.

The invention discloses a foundation pit supporting structure deformation and internal force calculation method. In the calculation method, an out-pit constraint stiffness matrix is introduced as a calculation factor, and a stiffness coefficient is iteratively calculated; the calculation method comprises the following steps: S1, constructing a foundation pit supporting structure calculation model; s2, a displacement-load total balance equation of the supporting structure under the undetermined supporting working condition of the foundation pit is determined; s3, a servo support is installed, prestress is loaded, a constraint stiffness matrix is introduced, a stiffness coefficient is calculated in an iterative mode, a displacement-load total balance equation of the supporting structure is established step by step, and the deformation condition of the supporting structure is analyzed and calculated; s4-5, the steps S2 to S3 are repeated until installation of all servo supports and prestress loading are completed, and displacement and internal force of the supporting structure are calculated; and S6, the foundation pit continues to be excavated to the pit bottom, deformation and internal force under all working conditions are calculated, and regulation and control calculation can be conducted in the process. By means of the calculation method, accurate calculation of deformation and internal force of the foundation pit supporting structure under the servo supporting effect can be achieved.

The invention is suitable for the field of portal frames, and discloses a portal frame dynamics analysis method which comprises the following steps: establishing a full-order finite element model of a portal frame, and obtaining a mass matrix and a stiffness matrix of the full-order finite element model; dividing the full-order finite element model into a plurality of substructures based on the modular structure characteristics of the portal frame, and decomposing the degree of freedom of each substructure into an internal degree of freedom and a boundary degree of freedom; constructing a fixed boundary main modal matrix and a constraint modal matrix of each substructure, and constructing a reduced-order basis matrix; performing modal polycondensation on a mass matrix and a stiffness matrix of the full-order finite element model by using the reduced-order basis matrix to generate a reduced-order mass matrix and a reduced-order stiffness matrix, and constructing a reduced-order kinetic equation; and solving the reduced-order kinetic equation to obtain reduced-order modal coordinates, carrying out inversion reduction in combination with the reduced-order basis matrix and the reduced-order modal coordinates to obtain a displacement response vector of the full-order finite element model, and retaining key kinetic information of the portal frame to the greatest extent through reduced-order efficiency improvement.

The invention belongs to the technical field of aircraft thin-wall structure design, and particularly relates to a method for correcting a quadrilateral shell unit in aircraft thin-wall structure finite element analysis, which comprises the following steps of: 1, establishing a unit coordinate system for the quadrilateral shell unit; step 2, under the unit coordinate system, establishing an equilibrium equation of a plane stress unit related to in-plane deformation, and calculating a rigidity matrix of the plane stress unit; step 3, under the unit coordinate system, establishing a balance equation of a bending plate unit related to out-of-plane deformation, and calculating a rigidity matrix of the bending plate unit; 4, applying an in-plane rotation penalty term on the rotational stiffness of the quadrilateral shell unit, and calculating the influence of the in-plane rotation penalty term on the strain energy function of the quadrilateral shell unit under the unit coordinate system to obtain an in-plane rotation stiffness matrix; 5, calculating a stiffness matrix of the quadrilateral shell unit under the unit coordinate system; and step 6, coordinate transformation is carried out to obtain a stiffness matrix of the quadrilateral shell unit under the overall coordinate system.

The invention discloses a robot rigidity modeling method and system based on joint module nonlinear rigidity parameter fitting, and the method comprises the steps: modeling the torsional rigidity of a joint module, and expressing the relation between the torsional angle of the joint module and the joint torque through a cubic polynomial function; modeling a robot rigidity model, and expressing a relation between a joint space rigidity matrix and a robot tail end rigidity matrix through robot kinematics; identifying a torsion angle and a torque of a joint module, and acquiring the torque of the joint module through a torsion sensor or based on joint current identification; a static flexibility experiment is conducted, and the corresponding relation between the tail end force vector and the tail end elastic deformation vector of the robot is obtained; non-linear stiffness parameter fitting identification is carried out, and non-linear stiffness parameter identification is realized through a least square problem model; according to the method, non-linear rigidity parameter identification of the joint module is achieved, whole machine rigidity modeling is completed on the basis, and a basis is provided for rigidity compensation of a robot and improvement of positioning precision.

The invention provides a robot variable impedance control method and system based on model predictive control and control equipment, and the method comprises the steps: constructing an iteration vector based on a current state vector and a control vector, and executing the iteration of the iteration vector under a plurality of set constraint conditions through employing an optimization model based on a non-linear objective function and the iteration vector, the plurality of constraint conditions comprise impedance parameters, upper and lower bound constraints of rigidity parameters, change rate constraints and passivity constraints. And extracting a control vector in the iteration vector, and obtaining an impedance matrix and a stiffness matrix based on the extracted control vector. And inputting the impedance matrix and the stiffness matrix into an impedance controller, and predicting to obtain a torque for executing robot control. According to the scheme, follow-up behaviors of the system can be predicted in a rolling mode in each control period, control input is optimized under the constraint condition, a systematic optimization strategy is provided for parameter adjustment, multi-target performance indexes and hardware constraints can be considered, and control robustness is improved.

The invention discloses an earthquake risk evaluation method and system, and relates to the technical field of earthquake numerical simulation, and the method comprises the steps: receiving a three-dimensional geometric model of a grid fault, and setting an initial condition; constructing a simulation operation equation to solve a shear stress loading rate, a normal stress loading rate and a state variable, and constructing an ordinary differential equation system; based on the three-dimensional geometric model and the initial conditions, accelerated simulation operation of the shear stress loading rate, the normal stress loading rate and the state variables is achieved through the ordinary differential equation system, and the accelerated simulation operation method comprises the steps that accelerated simulation operation is achieved through the method of combining an H-matrix and MPI, comprising the steps that a layered matrix is obtained based on a Green function under the half-space condition in combination with an adaptive cross-approximation algorithm, the layered matrix comprises a normal stiffness H-matrix and a shear stiffness H-matrix, and MPI parallel operation matrix-vector multiplication is adopted. According to the invention, the hierarchical matrix is introduced and the MPI parallel computing framework is fused, so that the simulation period is greatly shortened.

The invention relates to a method for predicting stress deformation of a confined space foundation pit support structure, which comprises the following steps of: S1, inputting soil parameters and foundation pit design parameters, and generating a stiffness matrix K; s2, determining an initial soil pressure parameter P, and calculating lateral deformation u of the retaining structure; s3: according to u, considering a foundation soil limited space, and updating P; calculating u again based on the updated P; s4, judging whether the error delta u meets the precision requirement or not, if so, turning to S5, otherwise, repeating S3, and iteratively calculating u until the error delta u meets the precision requirement; s5, after u is determined, soil shear stress is calculated, and then shear strain and soil shear rigidity are calculated; s6, calculating the elastic modulus k of the soil spring, and updating K; step S7, recording uT; then returning to S1, and entering the next excavation stage for calculation; and S8, after calculation is conducted to the designed excavation depth, P, u and the bending moment M are output, and the stress deformation condition of the foundation pit support structure is evaluated according to the P, u and the bending moment M. The method is beneficial for improving the precision of stress deformation prediction of the foundation pit support structure in the confined space.

The invention provides a calculation method and system for a herringbone micro anti-slide pile structure, the herringbone micro anti-slide pile structure is dispersed, and the interaction between a micro pile and a stratum is simplified into transverse and axial non-linear springs acting on nodes; determining nonlinear spring parameters; calculating the soil displacement and direction of each micro pile structure according to the landslide mass motion characteristics; establishing a structural overall stiffness matrix; establishing an overall increment balance equation set including the spring stiffness; the landslide mass displacement is divided into a plurality of displacement increments to be solved respectively, the spring stiffness and the pile bending stiffness are updated after each increment step, and the internal force and deformation of the herringbone micro pile structure under the corresponding landslide mass displacement effect are obtained through accumulation. According to the calculation method and system for the herringbone miniature anti-slide pile structure, material and geometric nonlinearity of the pile and interaction nonlinearity of the pile and the stratum can be considered at the same time, and calculation precision and efficiency are high.

The invention discloses a linear variation stiffness matrix method for suspension bridge space cable forming form finding, relates to the technical field of bridge numerical simulation, and solves the defects of a stiffness matrix in the process of solving the suspension bridge forming form finding through a segmented catenary method. The method comprises the steps that a suspension bridge is simplified into a space cable system composed of a main span, a side span and an anchor span, and the main span, the side span and the anchor span are sequentially subjected to a bridge forming and shape finding process; the linear change stiffness matrix solved in the bridge forming and shape finding process of the main span is used for representing the linear relation between the displacement change of the control point and the initial force change of the main cable; according to the method, an analytical expression of a linear change stiffness matrix is deduced, and an efficient iterative algorithm is established. According to the algorithm, the initial force correction is directly calculated according to the displacement error of the control point in each iteration.

The invention discloses a complex soft foundation waterscape miniature steel pipe pile reinforcing method which comprises the following steps: accurately predicting a deformation value of a potential settlement area by fusing thickness distribution parameters, soil compression modulus and numerical simulation, and comparing the deformation value with an engineering safety threshold. When deformation exceeds the standard, an optimization algorithm is adopted to adjust pile foundation arrangement, optimized pile foundation position parameters are generated, then connection node coordinates are extracted from the bottom plate structure model, and stress transmission efficiency and load distribution uniformity are calculated. And if the uniformity is lower than a threshold value, a reinforcement configuration scheme is finally generated by iteratively updating a connection node stiffness matrix, optimizing connection mechanism parameters, extracting mechanical transmission path data, simulating long-term stability and evaluating the durability improvement degree. According to the method, multi-level parameter optimization and numerical simulation are combined, so that the foundation stability and the load distribution uniformity are ensured, and the engineering durability and safety are remarkably improved.

The invention discloses a structural topology optimization design method for additive manufacturing and a related device, which can be used in the field of structural design, in the method, based on the shape of a structural design domain and a parameterized density field function including the structural design domain, a structural density field is determined, and a suspension angle constraint function at a structural boundary is constructed; constructing a boundary gradient constraint function based on a normal vector pointing to the interior of the structural design domain on the boundary of the structural design domain; discretizing the structural design domain into a plurality of finite element units, and constructing a material stiffness matrix to calculate an overall stiffness matrix; and constructing a topological optimization mathematical model according to the finite element balance constraint, the volume constraint, the suspension angle constraint function and the boundary gradient constraint function based on the overall stiffness matrix, and solving to obtain a structural topological optimization result. Therefore, through multi-constraint collaborative optimization, the conflict problem of suspended constraint and material anisotropy on the forming direction is solved, and a structural topological optimization result considering manufacturability and structural performance is obtained.

The invention relates to the technical field of computer-aided analysis, in particular to a floating foundation structure dynamic response real-time analysis system and method.The system comprises a residual vector calculation module for generating a residual force vector based on a full-order external excitation force vector, an increment base generation module for calculating a modal amplitude and constructing an increment base vector, and an increment base generation module for generating an increment base vector; the reduced-order matrix updating module constructs an enhanced reduced-order model by using the updated projection matrix, and the transient response solving module solves and obtains a corrected generalized coordinate value based on the enhanced reduced-order model and maps the corrected generalized coordinate value into a dynamic displacement field. According to the method, a residual vector is generated by calculating a difference value between full-order external excitation and internal force, residual force is mapped into static deformation distribution by using a stiffness matrix inverse operator, an incremental base is constructed through orthogonalization processing, adaptive updating of a projection space of a reduced-order model is realized, the system dimension is reduced, and the precision is guaranteed. And the calculation efficiency and the real-time performance of dynamic response analysis of the floating structure are improved.

The invention discloses a method for calculating broadband consistent damping force in time-domain dynamic analysis, which comprises the following steps of: constructing a finite element model of a hydroelectric plant, endowing the model with material attributes, and then obtaining an initial stiffness matrix of the hydroelectric plant by adopting a finite element method; idealizing the finite element model as a main spring, and connecting a plurality of combined elements consisting of sticky kettles and springs in parallel to form an idealized model; according to the initial stiffness matrix, the dynamic stiffness model of the combined element and a damping ratio in time domain dynamic analysis, calculating a stiffness coefficient of each combined element by adopting a least square method; when the dominant frequency of the combined element is uniformly distributed in the broadband range, calculating the viscosity coefficient of the combined element according to the dominant frequency of the combined element and the target damping ratio; and according to the rigidity coefficient and the viscosity coefficient, a time domain damping force model is adopted to calculate the consistent damping force of the hydroelectric plant in the broadband range.

The invention relates to a dam safety state intelligent prediction method and system based on regression modeling, and the method comprises the steps: carrying out the self-adaptive wavelet packet decomposition of structural data collected by a sensor, so as to extract a time sequence feature, and carrying out the differential coding of an environment variable, so as to generate a covariable feature; the method comprises the following steps: mapping a text event into vector features through a knowledge graph, inputting the features into a mixed primary function library formed by a physical primary function and a data primary function, and constructing a sparse regression model through a dynamic gating mechanism; a physical constraint term based on a stiffness matrix and an external load matching error is introduced into the model, a causal graph is dynamically updated through an incremental PC algorithm, and a causal regularization term is added into a loss function to improve the adaptability of the model. The method is suitable for dam structure health monitoring, intelligent early warning, long-term trend analysis and other application scenes, and has good real-time performance, interpretability and small sample generalization ability.