Coastal port pile foundation structure load calculation method and system based on coupling mode

Abstract

The invention discloses a coastal port pile foundation structure load calculation method based on a coupling mode, and the method comprises the steps: a hydrodynamic module which evaluates the load and mechanical properties of a port pile foundation structure under different environment conditions; the water power module at least comprises an atmospheric mode module, a sea wave mode module, an ocean mode module and a coupler module; wherein the coupler module exchanges and synchronizes data among the atmosphere mode module, the sea wave mode module and the ocean mode module; the data conversion module is used for converting the wind field data, the wave field data and the flow field data obtained by the hydrodynamic force module into aerodynamic force and hydrodynamic force input files of the pile foundation structure; the finite element analysis module is used for receiving the aerodynamic force and hydrodynamic force input files, establishing a port pile foundation structure model and calculating stress distribution and displacement response of a pile foundation structure, and simulation results output by the finite element analysis module include but not limited to a stress nephogram, displacement field distribution and a buckling phenomenon of a pile foundation.

Description

Technical Field

[0001] The invention relates to the technical field of port pile foundation load model analysis, and in particular to a method and system for calculating the load of a coastal port pile foundation structure based on a coupling mode. Background Art

[0002] With the deepening of global marine development activities, the design and construction of ports and docks are gradually developing towards deep water and large-scale development, and the environmental loads and safety risks faced by engineering structures have also increased significantly. Traditional dock design methods often rely on empirical formulas and simplified models. These methods have many shortcomings when dealing with load calculations in complex marine environments, such as ignoring the coupling effects of wind, waves, and currents, and being unable to accurately predict the structural stress conditions under extreme weather conditions. Therefore, there is an urgent need for a new calculation method that can better combine coupled numerical models and pile group structural mechanics models to accurately simulate loads in complex marine environments and provide a scientific basis for the design, construction, and operation of port terminals.

[0003] Therefore, a method and system for calculating the load of coastal port pile foundation structures based on a coupling mode is provided, which can simulate and calculate the load of port pile foundations in a complex marine environment. Summary of the invention

[0004] Therefore, in order to solve the above defects in the prior art, the present invention provides a method for calculating the load of coastal port pile foundation structure based on a coupling mode, so as to realize the simulation calculation of the port pile foundation load in a complex marine environment.

[0005] The present invention discloses a method for calculating the load of a coastal port pile foundation structure based on a coupling mode, comprising: A hydrodynamic module for evaluating the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module at least includes an atmospheric model module for simulating wind field data of the offshore dock area, a wave model module for calculating wave field data based on wind field data, an ocean model module for calculating flow field data based on the wind field data and wave field data, and a coupler module; wherein the coupler module exchanges and synchronizes data among the atmospheric model module, the wave model module, and the ocean model module by integrating multiple data exchange protocols and synchronization algorithms; A data conversion module converts the wind field data, wave field data and flow field data obtained by the hydrodynamic module into aerodynamic and hydrodynamic input files for the pile foundation structure; and The finite element analysis module receives the aerodynamic and hydrodynamic input files, calculates the stress distribution and displacement response of the pile foundation structure, and outputs simulation results including but not limited to the stress cloud diagram, displacement field distribution and buckling phenomenon of the pile foundation.

[0006] In the hydrodynamics module, the method for configuring model parameters and boundary conditions is: S21. Prepare atmospheric initial conditions and boundary conditions for the atmospheric model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the WRF model, and define key parameters including the size of the simulation domain, time step, physical process options, and output frequency; after configuration, the WRF model simulates atmospheric conditions and generates high-precision simulation results of the regional wind field. The meteorological simulation data includes wind speed, wind direction, temperature, and humidity; S22. Prepare ocean initial conditions and boundary conditions for the ocean model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the ROMS model, and set physical parameters including simulation domain, time step, vertical mixing and bottom drag number; the ROMS model receives wind field data from the WRF model, simulates ocean dynamic processes, and generates simulation data including ocean current velocity, temperature and salinity; S23. According to the requirement for the accuracy of predicting the ocean dynamic conditions in the target sea area, prepare the initial wave conditions and boundary conditions for the wave model module; configure the SWAN model to define the wave domain settings, boundary conditions and wind field input; the SWAN model receives the wind field data of the WRF model and the ocean flow field data of the ROMS model, and simulates the wave characteristics of the nearshore area through the input data. The generated simulation data includes wave height, period and propagation direction.

[0007] In the hydrodynamic module, the MCT coupler in the coupler module manages the data exchange between the WRF model, the ROMS model and the SWAN model. The MCT coupler realizes the bidirectional transmission of physical quantities between the modes by defining the grid mapping and communication paths between the modes, including: S31. The two-way coupling between the atmospheric model module and the ocean model module is as follows: the atmospheric boundary conditions provided by the WRF model, including but not limited to wind speed, air pressure, air temperature, humidity and precipitation data, are transmitted to the ROMS model for calculating the sea surface heat flux, momentum flux and fresh water flux, which affect the dynamic and thermal processes of the ocean model module; the ROMS model feeds back the calculated data including sea surface temperature and sea surface current velocity to the WRF model, which affects the thermal processes and wind field distribution in the atmospheric model module; S32. The two-way coupling between the ocean model module and the wave model module is as follows: the wave parameters calculated by the SWAN model are transmitted to the ROMS model for correcting the sea surface stress and vertical mixing coefficient, which affects the mixing layer depth and surface flow field of the ocean model module by changing the dynamic process of the sea surface; the ROMS model feeds back the calculated ocean current data to the SWAN model for correcting the wave propagation path and energy distribution, so that the wave propagation direction and energy concentration area in the wave model module change; S33. The two-way coupling between the atmospheric model module and the ocean wave model module is as follows: the wind field data provided by the WRF model drives the wave generation and propagation process in the SWAN model, and determines the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module can affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure and momentum exchange of the atmospheric model module by changing the sea surface roughness.

[0008] Finite element analysis modules include: S51. Importing the marine environment data processed by the data conversion module into the Fluent model, including but not limited to wind speed, wave height and ocean current speed; S52. Use the boundary condition setting function of the Fluent model to convert wind speed and wave height data into dynamic loads and apply them to the corresponding parts of the wharf structure; the ocean current velocity data is used as the fluid dynamic load input; S53. Accurately map the ocean model grid data to the nodes of the Fluent model to ensure that the spatial distribution of the data is consistent with reality; S54: After completing the boundary conditions and load input, the simulation is started. The Fluent model calculates the response of the port pile foundation according to the above environmental loads, including load, settlement and deformation.

[0009] The finite element analysis module uses the Morison equation to calculate the hydrodynamic loads:

[0010] Among them, the inertia force is:

[0011] The resistance is:

[0012] in, F(t) is the total hydrodynamic load; ρ is the water density; D is the pile diameter; U(t) is the water flow velocity; is the time derivative of water velocity; C m is the additional mass coefficient, C d is the resistance coefficient.

[0013] The dynamic response equation of the Fluent model under the simultaneous action of aerodynamic loads and hydrodynamic loads is: , Where M is the mass matrix, representing the inertia of the structure; is the acceleration vector; C is the damping matrix; is the velocity vector; K is the stiffness matrix, which represents the elastic restoring force of the structure; u(t) is the displacement vector; is the aerodynamic load; is the hydrodynamic load.

[0014] The calculation formula for pile foundation settlement is: , in, S is the total settlement of the pile foundation, P is the pile top load, A is the cross-sectional area of ​​the pile, E is the elastic modulus of the pile material, L is the length of the pile, λ It is a parameter related to the stiffness of the soil around the pile.

[0015] The present invention also provides a coastal port pile foundation structure load calculation system based on a coupling mode, and the system adopts the coastal port pile foundation structure load calculation method based on a coupling mode as described above for calculation.

[0016] The technical solution of the present invention has the following advantages: In the present invention, a high-precision multi-physical field coupling model is adopted. By fully considering the nonlinear interactions between the physical fields, the simulation accuracy of the wind, wave and flow field coupling effects is significantly improved, and the error accumulation and systematic deviation of the traditional single model method are greatly reduced; at the same time, a modular integrated design architecture is adopted, which has extremely high scalability and flexibility, and improves the system's computing efficiency and resource utilization; the present invention not only improves the simulation precision and spatiotemporal resolution, but also greatly shortens the calculation time of large-scale complex ocean dynamic processes, meets the forecast requirements with high real-time performance and strict precision requirements, and realizes high-precision simulation calculation of port pile foundation loads in complex marine environments. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0018] Figure 1 It is a schematic diagram of the load calculation method of the present invention; Figure 2 This is a schematic diagram of the hydrodynamic module described in Example 3; Figure 3 This is a schematic diagram of the typhoon wind field distribution cloud map near the pile foundation described in Example 3; Figure 4 This is a schematic diagram of the cloud map of the ocean current field distribution near the pile foundation described in Example 3; Figure 5 This is a schematic diagram of the cloud map of the ocean wave field distribution near the pile foundation described in Example 3; Figure 6 Schematic diagram of the void pressure calculation results of pile foundation points A and B in Example 3. DETAILED DESCRIPTION

[0019] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0020] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0021] Example 1: Figure 1 As shown, this embodiment provides a method for calculating the load of a coastal port pile foundation structure based on a coupling mode, including: A hydrodynamic module for evaluating the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module at least includes an atmospheric model module for simulating wind field data of the offshore dock area, a wave model module for calculating wave field data based on wind field data, an ocean model module for calculating flow field data based on the wind field data and wave field data, and a coupler module; wherein the coupler module exchanges and synchronizes data among the atmospheric model module, the wave model module, and the ocean model module by integrating multiple data exchange protocols and synchronization algorithms; The data conversion module converts the wind field data, wave field data and flow field data obtained by the hydrodynamic module into aerodynamic and hydrodynamic input files of the pile foundation structure. In this embodiment, the data conversion module can output data files that conform to the corresponding formats according to the requirements of different finite element analysis software. The data conversion module uses an algorithm based on fluid dynamics theory to convert wind field, wave field and flow field data into aerodynamic and hydrodynamic loads. It comprehensively considers the geometric characteristics, material properties and environmental conditions of the pile foundation structure to ensure that the generated aerodynamic and hydrodynamic input files can truly reflect the stress conditions of the pile foundation in a complex marine environment. and a finite element analysis module, which receives the aerodynamic and hydrodynamic input files, calculates the stress distribution and displacement response of the pile foundation structure, and the output simulation results include but are not limited to the stress cloud map, displacement field distribution and buckling phenomenon of the pile foundation, i.e., load, settlement and deformation; in this embodiment, the finite element analysis module can establish a refined model of the port pile foundation structure based on the three-dimensional finite element modeling technology; it takes into account the interaction force between the pile foundation and the soil, the nonlinear mechanical properties of the pile foundation, and the deformation and stress distribution under various load conditions; at the same time, the finite element analysis module can also output simulation results including the stress cloud map, displacement field distribution and possible buckling phenomenon of the pile foundation, providing a scientific basis for the design and safety assessment of the port pile foundation structure; wherein, the generated mechanical response analysis report includes the spatial distribution map of the shear force, bending moment, axial force and torque of the pile foundation, vibration spectrum analysis, and overall stability assessment of the structure. At the same time, the three-dimensional model of the port pile foundation structure can dynamically adjust the arrangement of the pile foundation, including but not limited to the diameter, length, inclination, spacing and inclination direction of the pile foundation, to adapt to different types of dock design requirements and construction conditions.

[0022] Specifically, Figure 2 As shown in the figure, the atmospheric model module generates high-precision simulation results of the regional wind field by receiving and processing real-time meteorological data. The meteorological data includes but is not limited to wind speed, wind direction, temperature, air pressure and humidity. The atmospheric model module uses the WRF (Weather Research and Forecasting) mode for calculation. The WRF mode realizes the simulation of atmospheric dynamic processes on a regional scale through high-resolution grid division. It can capture the changes in wind fields at all levels in the atmosphere. Especially in the case of extreme weather events such as typhoons, the atmospheric model module can provide accurate wind field data to provide basic data support for subsequent wave model and ocean model calculations. The ocean model module uses the ROMS (Regional Ocean Modeling System) model for calculations. The ROMS model accurately simulates the flow field changes at different depths in the sea area by considering various ocean dynamic processes such as tidal currents, eddies, and ocean currents. The ocean model module can also provide detailed three-dimensional flow field distribution in combination with the initial conditions and boundary conditions of the actual sea area, especially simulating the dynamic changes of ocean currents on multi-level vertical sections, providing detailed fluid mechanics data support for the hydrodynamic analysis of pile foundations. The wave model module uses the SWAN (Simulating Waves Nearshore) model to simulate the wave propagation process. The SWAN model is based on wind field data and can calculate the distribution of wave energy, wave height, wavelength, wave direction, and nonlinear interaction of waves. The wave model module is particularly suitable for simulating wave characteristics under complex sea conditions, including shallow water areas, areas with significant seabed topography, and areas with strong currents. The wave field data output by the SWAN model will serve as the key input data for subsequent hydrodynamic calculations. The MCT (Model Coupling Toolkit) coupler of the coupler module integrates multiple data exchange protocols and synchronization algorithms to achieve seamless connection between modes, ensuring that the calculation results of the atmospheric model, wave model and ocean model can remain consistent in time and space; the coupler also has an error correction function, which can correct data deviations caused by model differences or calculation errors in real time during data transmission, thereby improving the accuracy of the overall forecast model.

[0023] Furthermore, in the hydrodynamic module, the method for configuring model parameters and boundary conditions is as follows: S21. Prepare atmospheric initial conditions and boundary conditions for the atmospheric model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the WRF model, and define key parameters including the size of the simulation domain, time step, physical process options, and output frequency; after configuration, the WRF model simulates atmospheric conditions and generates high-precision simulation results of the regional wind field. The meteorological simulation data includes wind speed, wind direction, temperature, and humidity; S22. Prepare ocean initial conditions and boundary conditions for the ocean model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the ROMS model, and set physical parameters including simulation domain, time step, vertical mixing and bottom drag number; the ROMS model receives wind field data from the WRF model, simulates ocean dynamic processes, and generates simulation data including ocean current velocity, temperature and salinity; S23. According to the requirement for the accuracy of predicting the ocean dynamic conditions in the target sea area, prepare the initial wave conditions and boundary conditions for the wave model module; configure the SWAN model to define the wave domain settings, boundary conditions and wind field input; the SWAN model receives the wind field data of the WRF model and the ocean flow field data of the ROMS model, and simulates the wave characteristics of the nearshore area through the input data. The generated simulation data includes wave height, period and propagation direction.

[0024] Specifically, in the hydrodynamic module, the MCT coupler in the coupler module manages the data exchange between the WRF model, the ROMS model and the SWAN model. The MCT coupler realizes the bidirectional transmission of physical quantities between the modes by defining the grid mapping and communication paths between the modes, specifically including: S31. The two-way coupling between the atmospheric model module and the ocean model module is as follows: the atmospheric boundary conditions provided by the WRF model, including but not limited to wind speed, air pressure, air temperature, humidity and precipitation data, are transmitted to the ROMS model for calculating the sea surface heat flux, momentum flux and fresh water flux, which affect the dynamic and thermal processes of the ocean model module; the ROMS model feeds back the calculated data including sea surface temperature and sea surface current velocity to the WRF model, which affects the thermal processes and wind field distribution in the atmospheric model module; S32. The two-way coupling between the ocean model module and the wave model module is as follows: the wave parameters calculated by the SWAN model are transmitted to the ROMS model for correcting the sea surface stress and vertical mixing coefficient, which affects the mixing layer depth and surface flow field of the ocean model module by changing the dynamic process of the sea surface; the ROMS model feeds back the calculated ocean current data to the SWAN model for correcting the wave propagation path and energy distribution, so that the wave propagation direction and energy concentration area in the wave model module change; S33. The two-way coupling between the atmospheric model module and the ocean wave model module is as follows: the wind field data provided by the WRF model drives the wave generation and propagation process in the SWAN model, and determines the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module can affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure and momentum exchange of the atmospheric model module by changing the sea surface roughness.

[0025] Furthermore, the data conversion module includes: S41. The mesoscale wind field data file output by the mesoscale atmosphere-wave-ocean coupled forecast model is used to regenerate the wind field of the area around the dock using the full-field turbulent wind simulation software TurbSim to create wind speed data of a two-dimensional rectangular grid; the wind speed data includes the u component, v component and w component of the wind speed, and is saved in an input format suitable for the Fluent model.

[0026] S42. Extract the time series data of wave height, velocity components (u, v, w), acceleration and dynamic pressure at specific points from the wave field and flow field data of the mesoscale model; use these data to generate the small-scale wave field data file required by Fluent, and configure the corresponding hydrodynamic effects in the Fluent model; S43. Extract flow velocity and flow direction data from the mesoscale flow field data. These data will be used to set fluid dynamic loads in Fluent. Flow velocity and flow direction data help simulate the stress conditions of port pile foundations in actual marine environments. S44. Input the processed data into the Fluent model for numerical simulation to calculate the load, settlement and deformation response of the port pile foundation under complex marine conditions.

[0027] Specifically, the finite element analysis module includes: S51. Importing the marine environment data processed by the data conversion module into the Fluent model, including but not limited to wind speed, wave height and ocean current speed; S52. Use the boundary condition setting function of the Fluent model to convert wind speed and wave height data into dynamic loads and apply them to the corresponding parts of the wharf structure; the ocean current velocity data is used as the fluid dynamic load input; S53. Accurately map the ocean model grid data to the nodes of the Fluent model to ensure that the spatial distribution of the data is consistent with reality; S54: After completing the boundary conditions and load input, the simulation is started. The Fluent model calculates the response of the port pile foundation according to the above environmental loads, including load, settlement and deformation.

[0028] The finite element analysis module uses the Morison equation to calculate the hydrodynamic loads: , Among them, the inertia force is: , The resistance is: , in, F(t) is the total hydrodynamic load; ρ is the water density; D is the pile diameter; U(t) is the water flow velocity; is the time derivative of water velocity; C m is the additional mass coefficient, C d is the resistance coefficient.

[0029] The dynamic response equation of the Fluent model under the simultaneous action of aerodynamic loads and hydrodynamic loads is: , Where M is the mass matrix, representing the inertia of the structure; is the acceleration vector; C is the damping matrix; is the velocity vector; K is the stiffness matrix, which represents the elastic restoring force of the structure; u(t) is the displacement vector; is the aerodynamic load; is the hydrodynamic load.

[0030] The calculation formula for pile foundation settlement is: , in, S is the total settlement of the pile foundation, P is the pile top load, A is the cross-sectional area of ​​the pile, E is the elastic modulus of the pile material, L is the length of the pile, λ It is a parameter related to the stiffness of the soil around the pile.

[0031] The calculation of pile foundation deformation is as follows: Under the action of horizontal load H, the horizontal displacement of the pile foundation is δ Calculated by Winkler elastic foundation beam theory. According to the elastic beam theory, the horizontal displacement of the pile under the action of the foundation reaction force is δ(x) Distance from pile top x The expression at is: , in, EI is the bending stiffness of the pile foundation, E is the elastic modulus, I is the moment of inertia of the cross section, k s is the horizontal resistance coefficient of soil to pile foundation, δ(x) For the distance x Horizontal displacement of pile foundation at .

[0032] Embodiment 2: Based on the above embodiment, this embodiment provides a coastal port pile foundation structure load calculation system based on coupling mode, which adopts the coastal port pile foundation structure load calculation method based on coupling mode as described above for calculation.

[0033] Example 3: This example takes a port pile foundation as an example to further explain the present application in detail.

[0034] S1: Use satellite remote sensing or data from nearby meteorological observation stations to obtain environmental information about the sea area where the offshore wind farm is located. These data can provide key data such as wind fields, wave fields, flow fields, sea surface temperature, salinity, etc. that have wide coverage and temporal continuity. After long-term observation accumulation, they have high temporal and spatial resolution and consistency. First, select appropriate satellite remote sensing products, such as scatterometer wind field data, wave height data, sea surface temperature data, etc., to ensure that the environmental data obtained has sufficient accuracy and availability. Next, perform necessary processing on these data, including noise removal, error correction, interpolation processing, and spatial matching, to ensure the accuracy and consistency of the data. Finally, select representative data that can accurately reflect the characteristics of the target sea area from the processed data as input for the subsequent hydrodynamic model.

[0035] S2: According to actual needs, select an appropriate target sea area, and then use the meteorological data with a spatial resolution of 0.25°×0.25° provided by the European Center for Medium-Range Weather Forecasts (ECMWF), the global high-resolution coastline data provided by the British Ocean Data Center (BODC), and the global seafloor topography data provided by GEBCO to extract the water depth information of the target sea area, as well as the harmonic constants of the main tidal components obtained from the Global Tide and Surge Model (GTSM). By integrating these data resources, complete the settings of the mesoscale atmospheric model WRF, the shallow water wave numerical model SWAN, and the regional ocean model ROMS to ensure the sensitivity and accuracy of the model and meet the simulation requirements of the target sea area.

[0036] like Figure 2 As shown, the hydrodynamic module of this embodiment includes an atmospheric model module, an ocean model module, a wave model module and a coupler module. That is, a coupled numerical model system is constructed, which integrates the WRF, ROMS and SWAN modes, and performs parallel calculations through a high-performance computing platform to simulate the interaction between wind, waves and currents in a complex ocean environment; among them, the WRF mode is used to simulate the atmospheric dynamic process and obtain high-resolution wind field, temperature field and humidity field data; the ROMS mode is used to simulate the ocean dynamic process and obtain the ocean current speed, temperature and salinity distribution; the SWAN mode is used to simulate the nearshore wave dynamic process and calculate the wave height, period and propagation direction. The three modes interact and exchange information through the data coupling interface (i.e., the running-in period module), thereby realizing a high-precision simulation of the coupling effect of wind, waves and currents in the ocean environment.

[0037] Atmospheric model module: The WRF (Weather Research and Forecasting) model is used for calculations; in the initialization phase, it is necessary to prepare the initial atmospheric conditions and boundary conditions; these conditions are usually provided by global or regional meteorological reanalysis data, such as the ERA5 data of the European Centre for Medium-Term Weather Forecasts (ECMWF), and the data must be converted to the NetCDF format that can be read by the WRF model through a format conversion tool. Among them, when configuring the WRF model, it is necessary to define key parameters such as the size of the simulation domain, time step, physical process options (such as cloud microphysical processes, radiation transfer, boundary layer parameterization, etc.), and output frequency; the size of the simulation domain should be determined according to the scope of the target sea area to ensure coverage of the dock area and its surrounding environment; the choice of time step needs to balance the calculation accuracy and computing resource consumption, and is usually set to a few minutes to a few hours. After the configuration is completed, the WRF model begins to simulate atmospheric conditions and generate meteorological data such as wind speed, wind direction, temperature, and humidity; this data will provide key inputs for subsequent ocean models and wave models. The results of this embodiment are shown in the figure below. Figure 3 shown.

[0038] Ocean model module: ROMS (Regional Ocean Modeling System) mode is used for calculation; in the initialization stage, it is necessary to prepare the initial and boundary conditions of the ocean, including data such as temperature, salinity and current velocity; these data are usually obtained from observational data or ocean reanalysis products, such as HYCOM (Hybrid Coordinate Ocean Model) global ocean analysis reanalysis data. Among them, when configuring the ROMS model, it is necessary to set physical parameters such as simulation domain, time step, vertical mixing and bottom drag; the ROMS model can take into account various ocean dynamic processes such as tides, eddies, and ocean currents, and accurately simulate the flow field changes at different depths in the sea area. When the ROMS model is running, it will receive wind field data from the WRF model, simulate ocean dynamic processes, and generate data such as ocean current velocity, temperature and salinity. The results of this embodiment are as follows. Figure 4 shown.

[0039] Wave model module: Use the SWAN (Simulating Waves Nearshore) model to simulate the wave propagation process; in the initialization stage, it is necessary to prepare the initial wave conditions and boundary conditions, including wave height, period, propagation direction and other data; these data are obtained from field observations, satellite remote sensing or historical wave data. Among them, when configuring the SWAN model, it is necessary to define the wave domain settings, boundary conditions and wind field inputs; the SWAN model will receive the wind field data of the WRF model and the ocean flow field data of the ROMS model, and use these inputs to simulate the wave characteristics of the nearshore area, including key data such as wave height, period and direction. The results of this embodiment are as follows: Figure 5 shown.

[0040] Coupler module: The MCT (Model Coupling Toolkit) coupler is used to manage data exchange and synchronization between the three modes of WRF, ROMS and SWAN. The MCT coupler integrates multiple data exchange protocols and synchronization algorithms to achieve seamless connection between the modes. In the configuration phase, the grid mapping and communication paths between the modes need to be defined to ensure that data can be transmitted accurately and efficiently. During the model operation, the MCT coupler is responsible for managing data exchange and synchronization between WRF, ROMS and SWAN; specifically, it includes: S31. The atmospheric boundary conditions (such as wind speed, air pressure, temperature, humidity, etc.) provided by the WRF model are transmitted to the ROMS model to calculate the sea surface heat flux, momentum flux and fresh water flux; these physical quantities directly affect the dynamic and thermal processes in the ocean model; the sea surface temperature (SST), sea surface current velocity and other data calculated by the ROMS model are fed back to the WRF model, affecting the thermal processes and wind field distribution in the atmospheric model; this two-way coupling ensures the mutual influence and consistency of physical processes between the atmosphere and the ocean; S32. The wave parameters (such as wave height, wavelength, wave energy, etc.) calculated by the SWAN model are transmitted to the ROMS model to correct the sea surface stress and vertical mixing coefficient; these wave parameters affect the mixed layer depth and surface flow field of the ocean model by changing the dynamic process of the sea surface; the ocean current data calculated by the ROMS model are fed back to the SWAN model to correct the wave propagation path and energy distribution; the influence of ocean current on wave propagation characteristics is particularly significant, especially in strong current areas such as tidal currents or typhoon paths; S33. The wind field data provided by the WRF model (such as wind speed and direction at a height of 10 meters) directly drives the wave generation and propagation process in the SWAN model, and determines the intensity and direction of the waves. At the same time, the wave parameters in the SWAN model can affect the offshore wind field of the WRF model. Waves affect wind speed, turbulent structure and momentum exchange by changing the roughness of the sea surface. During severe storms, high waves will increase the roughness of the sea surface, thereby weakening the wind speed and having a significant impact on the wind field distribution.

[0041] Furthermore, the data conversion module is responsible for converting the wind field, wave field and flow field data output by the atmospheric model, ocean model and wave model into the aerodynamic and hydrodynamic input files of the pile foundation structure; this step is the key to connecting the physical model with the engineering structure analysis; specifically: S41. Using the mesoscale wind field data file output by the mesoscale atmosphere-wave-ocean coupled forecast model, the wind field of the area around the dock is regenerated using the full-field turbulent wind simulation software TurbSim. Create wind speed data of a two-dimensional rectangular grid, including the u component, v component, and w component of the wind speed, and save it in a format suitable for Fluent input.

[0042] S42. Extract the time series data of wave elevation, velocity components (u, v, w), acceleration and dynamic pressure at specific points from the wave field and flow field data of the mesoscale model. These data are used to generate the small-scale wave field data file required by Fluent, and configure the corresponding hydrodynamic effects in Fluent. At the same time, extract the flow velocity and flow direction data from the mesoscale flow field data. These data are input into Fluent as fluid dynamic loads to help simulate the stress conditions of port pile foundations in actual marine environments.

[0043] Specifically, the finite element analysis module can accurately simulate and analyze the mechanical response of the port pile foundation structure in a complex marine environment; by receiving the aerodynamic and hydrodynamic input files, the module uses three-dimensional finite element modeling technology to establish a highly refined port pile foundation structure model. Among them, the Morison equation used in the finite element analysis module to calculate the hydrodynamic load is: , Among them, the inertia force is: , The resistance is: , In the formula, F(t) is the total hydrodynamic load; ρ is the water density; D is the pile diameter; U(t) is the water flow velocity; is the time derivative of water velocity; C m is the additional mass coefficient, ranging from 1.5 to 2.0; C d is the resistance coefficient, which depends on the roughness of the pile and the Reynolds number and ranges from 0.7 to 1.2.

[0044] Furthermore, Fluent implements the dynamic response equation of the wind turbine pile foundation structure under the simultaneous action of the above-mentioned aerodynamic loads and hydrodynamic loads as follows: , Where M is the mass matrix, representing the inertia of the structure; is the acceleration vector; C is the damping matrix; is the velocity vector; K is the stiffness matrix, which represents the elastic restoring force of the structure; u(t) is the displacement vector; is the aerodynamic load; is the hydrodynamic load.

[0045] At the same time, the load calculation model also takes into account the nonlinear interaction between the pile foundation and the surrounding soil. Especially under complex geological conditions, it can accurately simulate the nonlinear deformation and failure mechanism of the pile foundation. The calculation formula for pile foundation settlement is: , in, S is the total settlement of the pile foundation, P is the pile top load, A is the cross-sectional area of ​​the pile, E is the elastic modulus of the pile material, L is the length of the pile, λ It is a parameter related to the stiffness of the soil around the pile.

[0046] Here, it is necessary to add that: The interaction between the port pile foundation structure and the surrounding soil is a key factor affecting the stability and durability of the structure. To accurately describe this interaction, the contact mechanics theory is introduced into the finite element model, and parameters such as friction coefficient, contact stiffness and lateral earth pressure are combined. Specifically, the friction force ( τ ) is calculated by the following formula: , in, is the normal stress; is the friction coefficient.

[0047] At the same time, due to the nonlinear characteristics of the material, the pile foundation structure will show nonlinear deformation and stress distribution when subjected to complex loads. For this reason, nonlinear constitutive relations such as the Mohr-Coulomb criterion are used in the finite element model to describe the yield behavior of the pile foundation material. Specifically, the Mohr-Coulomb yield criterion is expressed as: , in, σ 1 and σ 3 are the maximum and minimum principal stresses, σ c is the cohesion of the material, ϕ is the internal friction angle.

[0048] In the finite element analysis, the influence of various load conditions on the pile foundation structure is considered, such as wind load, wave load and dynamic load caused by flow field; the superposition effect of these loads may lead to complex stress distribution and deformation mode of the pile foundation structure. To this end, the stress-strain relationship equation is used to calculate the stress state and displacement response of the structure under different working conditions; specifically, the stress-strain relationship is expressed by the following linear elastic equation: , in, σ is stress, E is the elastic modulus, ϵ For strain.

[0049] Furthermore, in the Fluent environment, establishing a detailed three-dimensional finite element model of the port pile foundation structure is the basis of the entire analysis. During the modeling process, it is necessary to fully consider the geometric dimensions, material properties, and contact relationship between the pile foundation and the surrounding soil. The details are as follows: First, the geometric dimensions of the pile foundation need to be accurately input based on actual engineering data, including pile diameter, pile length, and burial depth. These parameters directly affect the accuracy of the model and the reliability of the simulation results. Secondly, the definition of material properties includes the elastic modulus, Poisson's ratio, density, etc. of the pile foundation material. These parameters will affect the deformation and stress response of the structure under different loads. In order to improve the accuracy of the model, it is necessary to describe the mechanical properties of different materials in a refined manner, such as using a piecewise linear model to characterize the elastic-plastic behavior of the material. In addition, the accuracy of mesh division is also an important factor affecting the calculation results. The density, shape and distribution of the mesh must match the complexity of the structure to ensure the accuracy and convergence of the calculation results. Finally, the model boundary conditions include the contact relationship between the pile foundation and the soil, the fixed end conditions and the free end conditions. The setting of the boundary conditions not only determines the stress distribution of the model, but also affects the convergence of the calculation process. In actual operation, by adjusting the mesh division, boundary conditions and material parameters many times, a three-dimensional finite element model that meets the accuracy requirements is finally obtained.

[0050] The processed wind field, wave field and flow field data are imported into Fluent, and these data are converted into dynamic loads applied to the port pile foundation structure through the boundary condition setting function of the software; in this process, the environmental load must be accurately applied to the corresponding parts of the model to ensure that the direction, size and point of action of the load are consistent with the actual engineering conditions.

[0051] In order to realize the complex interaction between pile foundation and soil, Fluent software provides a variety of contact algorithms, such as Coulomb friction model, stick-slip contact model, etc. These models can accurately simulate the effects of friction, lateral soil pressure and end bearing force, especially considering the displacement of pile foundation in vertical and horizontal directions and the contact stress changes caused by it. In order to improve the simulation accuracy, the contact stiffness parameter is introduced in the contact algorithm. This parameter is closely related to the properties of the contact surface between pile foundation and soil, and can effectively reflect the contact mechanical behavior under different soil conditions.

[0052] After completing the model construction and load setting, start the Fluent simulation calculation process. During the simulation process, the software will calculate the mechanical response of the port pile foundation structure according to the input environmental load conditions, including the stress distribution, displacement response, settlement, and buckling phenomenon of the pile foundation.

[0053] In addition, to ensure the accuracy and convergence of the simulation results, the time step and calculation parameters need to be set reasonably during the calculation. A time step that is too large will cause distortion of the calculation results, while a time step that is too small may cause a significant increase in the calculation time, or even cause the calculation to not converge. To this end, by gradually adjusting the time step and parameter settings, a balance point can be found so that the simulation results are both accurate and engineering practical.

[0054] After the simulation calculation is completed, a detailed mechanical response analysis report is generated. Figure 6 As shown, it is a schematic diagram of the calculation results of the void pressure at two points A and B of the pile foundation calculated by this embodiment. Specifically, by viewing the stress cloud map and displacement field distribution map of the pile foundation, the stress condition and deformation characteristics of the pile foundation under different load conditions can be intuitively understood. In particular, the stress concentration area and displacement mutation point under extreme weather conditions are paid attention to. Through vibration spectrum analysis, the vibration characteristics and stability of the pile foundation under dynamic load can be evaluated, focusing on parameters such as resonance frequency and damping ratio to ensure the long-term stability and durability of the structure in the marine environment; according to the simulation results and mechanical response analysis report, targeted structural optimization suggestions are proposed. For example, the diameter, length, inclination and other parameters of the pile foundation are adjusted to improve the stress performance and stability of the structure; the contact relationship between the pile foundation and the soil is strengthened to increase friction and lateral earth pressure, etc. At the same time, in order to ensure the accuracy and reliability of the calculation results, the load calculation model needs to be verified and optimized. The model output is verified by comparing field observation data, historical statistical data and other independent verification sources, and the model parameters and structural settings are adjusted according to the verification results to improve the prediction accuracy and scope of application of the model.

[0055] In summary, the present application achieves at least the following effects compared to the prior art: (1) High-precision multi-physics coupling simulation. A high-precision multi-physics coupling model is used, integrating the WRF, ROMS and SWAN models; this integrated coupling strategy can capture the multi-scale complex dynamic processes under extreme weather such as typhoons, including the evolution of wind fields in the atmospheric boundary layer, the dynamic response of the ocean interior, and the generation and propagation of offshore waves. By fully considering the nonlinear interactions between the various physical fields, the simulation accuracy of the coupling effects of wind, waves and flow fields is significantly improved. Especially under strong typhoon conditions, the system can accurately reproduce the dynamic evolution process of the complex ocean environment, greatly reducing the error accumulation and systematic deviation of the traditional single model method.

[0056] (2) Advanced modular integrated design. The modular integrated design architecture enables independent development and flexible combination of various functional modules. For example, the AeroDyn module is responsible for aerodynamic calculations, the HydroDyn module handles hydrodynamics, and the SubDyn module realizes the coordinated coupling of multiple physical fields. Through this modular design, the overall architecture of the system has extremely high scalability and flexibility. Users can adjust module configurations or introduce new physical processes according to specific application requirements, ensuring the wide applicability and foresight of the model, while improving the system's computing efficiency and resource utilization.

[0057] (3) Parallel computing optimization. In terms of computing efficiency, an advanced parallel computing framework is used, which enables the system to adjust the grid resolution in real time according to the gradient changes of the physical field and optimize the allocation of computing resources. This not only improves the precision and spatiotemporal resolution of the simulation, but also greatly shortens the calculation time of large-scale complex ocean dynamic processes. This strategy of combining parallel processing with adaptive grids significantly improves the computing performance of the model and meets the forecasting requirements with high real-time performance and strict accuracy requirements.

[0058] Obviously, the above embodiments are merely examples for the purpose of clear explanation, and are not intended to limit the implementation methods. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the implementation methods here. The obvious changes or modifications derived therefrom are still within the scope of protection of the invention.

Claims

1. A method for calculating the load of pile foundation structures of coastal ports based on a coupling model, characterized in that: include: A hydrodynamic module for evaluating the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module at least includes an atmospheric model module for simulating wind field data of the offshore dock area, a wave model module for calculating wave field data based on wind field data, an ocean model module for calculating flow field data based on the wind field data and wave field data, and a coupler module; wherein the coupler module exchanges and synchronizes data among the atmospheric model module, the wave model module, and the ocean model module by integrating multiple data exchange protocols and synchronization algorithms; A data conversion module converts the wind field data, wave field data and flow field data obtained by the hydrodynamic module into aerodynamic and hydrodynamic input files for the pile foundation structure; and The finite element analysis module receives the aerodynamic and hydrodynamic input files, establishes a port pile foundation structure model, and calculates the stress distribution and displacement response of the pile foundation structure. The output simulation results include but are not limited to the stress cloud diagram, displacement field distribution and buckling phenomenon of the pile foundation.

2. The method for calculating the load of coastal port pile foundation structure based on coupling mode according to claim 1 is characterized in that: In the hydrodynamics module, the method for configuring model parameters and boundary conditions is: S21. Prepare atmospheric initial conditions and boundary conditions for the atmospheric model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the WRF model, and define key parameters including the size of the simulation domain, time step, physical process options, and output frequency; after configuration, the WRF model simulates atmospheric conditions and generates high-precision simulation results of the regional wind field. The meteorological simulation data includes wind speed, wind direction, temperature, and humidity; S22. Prepare ocean initial conditions and boundary conditions for the ocean model module according to the accuracy requirements for predicting ocean dynamic conditions in the target sea area; configure the ROMS model, and set physical parameters including simulation domain, time step, vertical mixing and bottom drag number; the ROMS model receives wind field data from the WRF model, simulates ocean dynamic processes, and generates simulation data including ocean current velocity, temperature and salinity; S23. According to the requirement for the accuracy of predicting the ocean dynamic conditions in the target sea area, prepare the initial wave conditions and boundary conditions for the wave model module; configure the SWAN model to define the wave domain settings, boundary conditions and wind field input; the SWAN model receives the wind field data of the WRF model and the ocean flow field data of the ROMS model, and simulates the wave characteristics of the nearshore area through the input data. The generated simulation data includes wave height, period and propagation direction.

3. The method for calculating the load of coastal port pile foundation structure based on coupling mode according to claim 2 is characterized in that: In the hydrodynamic module, the MCT coupler in the coupler module manages the data exchange between the WRF model, the ROMS model and the SWAN model. The MCT coupler realizes the bidirectional transmission of physical quantities between the modes by defining the grid mapping and communication paths between the modes, including: S31. The two-way coupling between the atmospheric model module and the ocean model module is as follows: the atmospheric boundary conditions provided by the WRF model, including but not limited to wind speed, air pressure, air temperature, humidity and precipitation data, are transmitted to the ROMS model for calculating the sea surface heat flux, momentum flux and fresh water flux, which affect the dynamic and thermal processes of the ocean model module; the ROMS model feeds back the calculated data including sea surface temperature and sea surface current velocity to the WRF model, which affects the thermal processes and wind field distribution in the atmospheric model module; S32. The two-way coupling between the ocean model module and the wave model module is as follows: the wave parameters calculated by the SWAN model are transmitted to the ROMS model for correcting the sea surface stress and vertical mixing coefficient, which affects the mixing layer depth and surface flow field of the ocean model module by changing the dynamic process of the sea surface; the ROMS model feeds back the calculated ocean current data to the SWAN model for correcting the wave propagation path and energy distribution, so that the wave propagation direction and energy concentration area in the wave model module change; S33. The two-way coupling between the atmospheric model module and the ocean wave model module is as follows: the wind field data provided by the WRF model drives the wave generation and propagation process in the SWAN model, and determines the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module can affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure and momentum exchange of the atmospheric model module by changing the sea surface roughness.

4. The method for calculating the load of coastal port pile foundation structure based on coupling mode according to claim 1 is characterized in that: Finite element analysis modules include: S51. Importing the marine environment data processed by the data conversion module into the Fluent model, including but not limited to wind speed, wave height and ocean current speed; S52. Use the boundary condition setting function of the Fluent model to convert wind speed and wave height data into dynamic loads and apply them to the corresponding parts of the wharf structure; the ocean current velocity data is used as the fluid dynamic load input; S53. Accurately map the ocean model grid data to the nodes of the Fluent model to ensure that the spatial distribution of the data is consistent with reality; S54: After completing the boundary conditions and load input, the simulation is started. The Fluent model calculates the response of the port pile foundation according to the above environmental loads, including load, settlement and deformation.

5. The method for calculating the load of pile foundation structure of coastal ports based on coupling mode according to claim 4 is characterized in that: The finite element analysis module uses the Morison equation to calculate the hydrodynamic loads: ， Among them, the inertia force is: , The resistance is: , in, F(t) is the total hydrodynamic load; ρ is the water density; D is the pile diameter; U(t) is the water flow velocity; is the time derivative of water velocity; C m is the additional mass coefficient, C d is the resistance coefficient.

6. The method for calculating the load of pile foundation structure of coastal ports based on coupling mode according to claim 5 is characterized in that: The dynamic response equation of the Fluent model under the simultaneous action of aerodynamic loads and hydrodynamic loads is: ， Where M is the mass matrix, representing the inertia of the structure; is the acceleration vector; C is the damping matrix; is the velocity vector; K is the stiffness matrix, which represents the elastic restoring force of the structure; u(t) is the displacement vector; is the aerodynamic load; is the hydrodynamic load.

7. The method for calculating the load of pile foundation structure of coastal ports based on coupling mode according to claim 6 is characterized in that: The calculation formula for pile foundation settlement is: ， in, S is the total settlement of the pile foundation, P is the pile top load, A is the cross-sectional area of ​​the pile, E is the elastic modulus of the pile material, L is the length of the pile, λ It is a parameter related to the stiffness of the soil around the pile.

8. The load calculation system for coastal port pile foundation structure based on coupling mode is characterized by: The system uses the coastal port pile foundation structure load calculation method based on the coupling mode described in any one of claims 1-7 for calculation.

Images