Load calculation system for coastal port pile foundation structures based on coupling model

Abstract

The present invention discloses a load calculation system for coastal port pile foundation structures based on a coupling mode, comprising: a hydrodynamic module, which evaluates the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module at least includes an atmospheric mode module, a wave mode module, an ocean mode module and a coupler module; wherein the coupler module exchanges and synchronizes data between the atmospheric mode module, the wave mode module and the ocean mode module; a data conversion module, which converts 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; and a finite element analysis module, which 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, wherein the simulation results output include but are not limited to the stress cloud map, displacement field distribution and buckling phenomenon of the pile foundation.

Description

Technical Field

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

[0002] With the deepening of global ocean development activities, the design and construction of port terminals are gradually developing towards deeper water and larger scale, and the environmental loads and safety risks faced by engineering structures have also increased significantly. Traditional terminal design methods often rely on empirical formulas and simplified models. These methods have many shortcomings when dealing with load calculations in complex marine environments. For example, they ignore the coupling effects of wind, waves, and currents and cannot 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 coastal port pile foundation structure load calculation system based on a coupling mode is provided, which can simulate and calculate the port pile foundation load 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 coastal port pile foundation structure load calculation system based on a coupling mode to realize the simulation calculation of the port pile foundation load in a complex marine environment.

[0005] The present invention discloses a coastal port pile foundation structure load calculation system based on a coupling mode, comprising:

[0006] A hydrodynamic module, which evaluates the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module includes at least an atmospheric model module for simulating wind field data in the offshore terminal area, a wave model module for calculating wave field data based on the 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;

[0007] 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

[0008] The finite element analysis module receives the aerodynamic and hydrodynamic input files and calculates the stress distribution and displacement response of the pile foundation structure. The simulation results output include but are not limited to the stress cloud diagram, displacement field distribution and buckling phenomenon of the pile foundation.

[0009] In the hydrodynamics module, the method for configuring model parameters and boundary conditions is:

[0010] S21. Prepare atmospheric initial and boundary conditions for the atmospheric model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the WRF model, defining key parameters including the size of the simulation domain, time step, physical process options, and output frequency. Once configured, the WRF model simulates atmospheric conditions, including wind speed, direction, temperature, and humidity, to generate high-precision simulation results for the regional wind field.

[0011] S22. Prepare ocean initial and boundary conditions for the ocean model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the ROMS model, setting physical parameters including the simulation domain, time step, vertical mixing, and bottom drag number. The ROMS model receives wind field data from the WRF model to simulate ocean dynamic processes. The generated simulation data includes ocean current velocity, temperature, and salinity.

[0012] S23. Prepare initial wave conditions and boundary conditions for the wave model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the SWAN model to define the wave domain settings, boundary conditions, and wind field inputs. The SWAN model receives wind field data from the WRF model and ocean current data from the ROMS model, and uses this input data to simulate wave characteristics in the nearshore area. The generated simulation data includes wave height, period, and propagation direction.

[0013] In the hydrodynamic module, the MCT coupler in the coupler module manages data exchange between the WRF model, the ROMS model, and the SWAN model. The MCT coupler implements bidirectional transmission of physical quantities between the modes by defining grid mapping and communication paths between the modes. Specifically, the following are performed:

[0014] S31. The bidirectional coupling between the atmospheric model module and the ocean model module is as follows: 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 calculation of sea surface heat flux, momentum flux, and freshwater flux, which influence the dynamic and thermodynamic processes in the ocean model module. The ROMS model feeds back calculated data, including sea surface temperature and sea surface current velocity, to the WRF model, which influences the thermodynamic processes and wind field distribution in the atmospheric model module.

[0015] S32. The bidirectional 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 to modify the sea surface stress and vertical mixing coefficient. This affects the mixed layer depth and surface flow field of the ocean model module by changing the dynamic processes of the sea surface. The ROMS model feeds back the calculated ocean current data to the SWAN model to modify the wave propagation path and energy distribution, resulting in changes in the wave propagation direction and energy concentration area in the wave model module.

[0016] S33. The bidirectional 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, determining the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure, and momentum exchange in the atmospheric model module by changing the sea surface roughness.

[0017] Finite element analysis modules include:

[0018] 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;

[0019] 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; ocean current velocity data is used as the fluid dynamic load input;

[0020] S53. Accurately map ocean model grid data to nodes in the Fluent model to ensure that the spatial distribution of the data is consistent with reality;

[0021] S54: After completing the boundary conditions and load input, the simulation is started. The Fluent model calculates the response of the port pile foundation based on the above environmental loads, including load, settlement, and deformation.

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

[0023] ,

[0024] Among them, the inertial force is: ,

[0025] The resistance is: ,

[0026] 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 drag coefficient.

[0027] The dynamic response equation of the Fluent model under the simultaneous action of aerodynamic loads and hydrodynamic loads is:

[0028] ,

[0029] 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:

[0031] ,

[0032] 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.

[0033] The technical solution of the present invention has the following advantages:

[0034] In the present invention, a high-precision multi-physics field coupling model is adopted. By fully considering the nonlinear interactions between various 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 fineness and spatiotemporal resolution, but also greatly shortens the calculation time of large-scale complex ocean dynamic processes, meets the forecast requirements of high real-time and strict precision requirements, and realizes high-precision simulation calculation of port pile foundation loads in complex marine environments. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0036] Figure 1 Schematic diagram of the load calculation method of the present invention;

[0037] Figure 2 This is a schematic diagram of the hydrodynamic module described in Example 2;

[0038] Figure 3 Schematic diagram of the void pressure calculation results of pile foundation points A and B in Example 2. DETAILED DESCRIPTION

[0039] The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described are only some embodiments of the present invention, not all embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

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

[0041] Example 1: Figure 1 As shown, this embodiment provides a coastal port pile foundation structure load calculation system based on a coupling mode, including:

[0042] A hydrodynamic module, which evaluates the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module includes at least an atmospheric model module for simulating wind field data in the offshore terminal area, a wave model module for calculating wave field data based on the 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;

[0043] 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. In this embodiment, the data conversion module can output data files that conform to 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.

[0044] A finite element analysis module receives the aerodynamic and hydrodynamic input files and calculates the stress distribution and displacement response of the pile foundation structure. The output simulation results include but are not limited to the pile foundation stress cloud diagram, displacement field distribution, and buckling phenomenon, namely 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 three-dimensional finite element modeling technology. It takes into account the interaction 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 pile foundation stress cloud diagram, displacement field distribution, and possible buckling phenomenon, providing a scientific basis for the design and safety assessment of the port pile foundation structure. The generated mechanical response analysis report includes the spatial distribution diagram 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 pile foundation layout, including but not limited to the pile foundation diameter, length, inclination angle, spacing, and tilt direction, to adapt to different types of terminal design requirements and construction conditions.

[0045] Specifically, if 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) model for calculation. The WRF model simulates atmospheric dynamic processes on a regional scale through high-resolution grid division. It can capture wind field changes 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, providing basic data support for subsequent wave model and ocean model calculations.

[0046] The ocean model module uses the Regional Ocean Modeling System (ROMS) model for calculations. The ROMS model accurately simulates the flow field changes at different depths within 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 based on the initial and boundary conditions of the actual sea area. In particular, it can simulate the dynamic changes of ocean currents on multiple vertical sections, providing detailed fluid mechanics data support for the hydrodynamic analysis of pile foundations.

[0047] The ocean wave model module uses the SWAN (Simulating Waves Nearshore) model to simulate wave propagation. Based on wind field data, the SWAN model can calculate wave energy distribution, wave height, wavelength, wave direction, and nonlinear wave interactions. The ocean wave model module is particularly suitable for simulating wave characteristics in 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 key input data for subsequent hydrodynamic calculations.

[0048] The MCT (Model Coupling Toolkit) coupler in the coupler module integrates multiple data exchange protocols and synchronization algorithms to achieve seamless connection between models, ensuring that the calculation results of the atmospheric model, wave model, and ocean model are 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.

[0049] Furthermore, in the hydrodynamics module, the method for configuring model parameters and boundary conditions is as follows:

[0050] S21. Prepare atmospheric initial and boundary conditions for the atmospheric model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the WRF model, defining key parameters including the size of the simulation domain, time step, physical process options, and output frequency. Once configured, the WRF model simulates atmospheric conditions, including wind speed, direction, temperature, and humidity, to generate high-precision simulation results for the regional wind field.

[0051] S22. Prepare ocean initial and boundary conditions for the ocean model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the ROMS model, setting physical parameters including the simulation domain, time step, vertical mixing, and bottom drag number. The ROMS model receives wind field data from the WRF model to simulate ocean dynamic processes. The generated simulation data includes ocean current velocity, temperature, and salinity.

[0052] S23. Prepare initial wave conditions and boundary conditions for the wave model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the SWAN model to define the wave domain settings, boundary conditions, and wind field inputs. The SWAN model receives wind field data from the WRF model and ocean current data from the ROMS model, and uses this input data to simulate wave characteristics in the nearshore area. The generated simulation data includes wave height, period, and propagation direction.

[0053] Specifically, in the hydrodynamic module, the MCT coupler in the coupler module manages data exchange between the WRF model, the ROMS model, and the SWAN model. The MCT coupler implements bidirectional transmission of physical quantities between the modes by defining grid mapping and communication paths between the modes. Specifically, the following steps are performed:

[0054] S31. The bidirectional coupling between the atmospheric model module and the ocean model module is as follows: 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 calculation of sea surface heat flux, momentum flux, and freshwater flux, which influence the dynamic and thermodynamic processes in the ocean model module. The ROMS model feeds back calculated data, including sea surface temperature and sea surface current velocity, to the WRF model, which influences the thermodynamic processes and wind field distribution in the atmospheric model module.

[0055] S32. The bidirectional 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 to modify the sea surface stress and vertical mixing coefficient. This affects the mixed layer depth and surface flow field of the ocean model module by changing the dynamic processes of the sea surface. The ROMS model feeds back the calculated ocean current data to the SWAN model to modify the wave propagation path and energy distribution, resulting in changes in the wave propagation direction and energy concentration area in the wave model module.

[0056] S33. The bidirectional 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, determining the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure, and momentum exchange in the atmospheric model module by changing the sea surface roughness.

[0057] Furthermore, the data conversion module includes:

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

[0059] S42. Extract time series data of wave height, velocity components (u, v, w), acceleration, and dynamic pressure at specific locations from the wave and flow field data of the mesoscale model; use this data to generate the small-scale wave field data file required by Fluent, and configure the corresponding hydrodynamic effects in the Fluent model;

[0060] S43. Extract flow velocity and direction data from the mesoscale flow field data. This data will be used to set fluid dynamic loads in Fluent. This flow velocity and direction data helps simulate the stress conditions of port pile foundations in actual marine environments.

[0061] 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.

[0062] Specifically, the finite element analysis module includes:

[0063] 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;

[0064] 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; ocean current velocity data is used as the fluid dynamic load input;

[0065] S53. Accurately map ocean model grid data to nodes in the Fluent model to ensure that the spatial distribution of the data is consistent with reality;

[0066] S54: After completing the boundary conditions and load input, the simulation is started. The Fluent model calculates the response of the port pile foundation based on the above environmental loads, including load, settlement, and deformation.

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

[0068] ,

[0069] Among them, the inertial force is: ,

[0070] The resistance is: ,

[0071] 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 drag coefficient.

[0072] The dynamic response equation of the Fluent model under the simultaneous action of aerodynamic loads and hydrodynamic loads is:

[0073] ,

[0074] 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.

[0075] The calculation formula for pile foundation settlement is: ,

[0076] 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.

[0077] The calculation of pile foundation deformation is as follows:

[0078] 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:

[0079] ,

[0080] in, IT is the flexural 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 the pile foundation.

[0081] Example 2: This example takes a port pile foundation as an example to further explain this application in detail.

[0082] 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. with 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 have 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.

[0083] S2: Based on actual needs, an appropriate target sea area is selected. Then, using meteorological data with a spatial resolution of 0.25° × 0.25° provided by the European Centre for Medium-Range Weather Forecasts (ECMWF), high-resolution global coastline data from the UK Ocean Data Centre (BODC), and global seafloor topography data from GEBCO, the target sea area's water depth information is extracted. The harmonic constants of the main tidal components are obtained from the Global Tide and Surge Model (GTSM). By integrating these data resources, the settings of the WRF mesoscale atmospheric model, the SWAN shallow water wave numerical model, and the ROMS regional ocean model are completed to ensure model sensitivity and accuracy to meet the simulation requirements of the target sea area.

[0084] like Figure 2As shown, the hydrodynamic module of this embodiment includes an atmospheric model module, an ocean model module, a wave model module, and a coupler module. This constructs a coupled numerical model system integrating the WRF, ROMS, and SWAN models. Parallel computation is performed on a high-performance computing platform to simulate the interactions between wind, waves, and currents in complex ocean environments. The WRF model simulates atmospheric dynamic processes, acquiring high-resolution wind, temperature, and humidity field data; the ROMS model simulates ocean dynamic processes, acquiring current velocity, temperature, and salinity distribution; and the SWAN model simulates nearshore wave dynamic processes, calculating wave height, period, and propagation direction. These three models interact and exchange information through a data coupling interface (i.e., the run-in module), enabling high-precision simulation of the coupled wind, wave, and current interactions in the ocean environment.

[0085] The atmospheric model module uses the WRF (Weather Research and Forecasting) model for calculations. During the initialization phase, atmospheric initial and boundary conditions must be prepared. These conditions are typically provided by global or regional meteorological reanalysis data, such as the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 data. This data must be converted to the NetCDF format readable by the WRF model using a format conversion tool. When configuring the WRF model, key parameters must be defined, including the size of the simulation domain, time step, physical process options (such as cloud microphysics, radiative transfer, boundary layer parameterization), and output frequency. The size of the simulation domain should be determined based on the target sea area to ensure coverage of the port area and its surroundings. The time step should balance computational accuracy and resource consumption, typically ranging from several minutes to several hours. Once configured, the WRF model begins simulating atmospheric conditions, generating meteorological data such as wind speed, direction, temperature, and humidity. This data provides key input for subsequent ocean and wave models.

[0086] Ocean Model Module: Calculations are performed using the Regional Ocean Modeling System (ROMS). During the initialization phase, initial ocean and boundary conditions, including temperature, salinity, and current velocity, are required. These data are typically obtained from observational data or ocean reanalysis products, such as the HYCOM (Hybrid Coordinate Ocean Model) global ocean analysis reanalysis data. When configuring the ROMS model, physical parameters such as the simulation domain, time step, vertical mixing, and bottom drag must be set. The ROMS model can account for various ocean dynamic processes, such as tidal currents, eddies, and currents, accurately simulating flow field variations at different depths within the ocean. During operation, the ROMS model receives wind field data from the WRF model to simulate ocean dynamic processes and generate data such as ocean current velocity, temperature, and salinity.

[0087] The Ocean Wave Model Module uses the SWAN (Simulating Waves Nearshore) model to simulate wave propagation. During the initialization phase, initial wave and boundary conditions must be prepared, including data such as wave height, period, and propagation direction. This data is obtained from field observations, satellite remote sensing, or historical wave data. When configuring the SWAN model, the wave domain settings, boundary conditions, and wind field inputs must be defined. The SWAN model receives wind field data from the WRF model and ocean current data from the ROMS model, using these inputs to simulate nearshore wave characteristics, including key data such as wave height, period, and direction.

[0088] Coupler Module: The MCT (Model Coupling Toolkit) coupler is used to manage data exchange and synchronization between the WRF, ROMS, and SWAN modes. The MCT coupler integrates multiple data exchange protocols and synchronization algorithms to achieve seamless integration between the modes. During the configuration phase, grid mapping and communication paths between the modes must be defined to ensure accurate and efficient data transfer. During model execution, the MCT coupler manages data exchange and synchronization between WRF, ROMS, and SWAN. Specifically, it includes:

[0089] S31. The atmospheric boundary conditions (such as wind speed, pressure, temperature, and humidity) provided by the WRF model are transmitted to the ROMS model for calculation of sea surface heat flux, momentum flux, and freshwater flux. These physical quantities directly influence the dynamic and thermodynamic processes in the ocean model. Data such as sea surface temperature (SST) and sea surface current velocity calculated by the ROMS model are fed back to the WRF model, influencing the thermodynamic 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.

[0090] S32. Wave parameters calculated by the SWAN model (such as wave height, wavelength, and wave energy) are transmitted to the ROMS model to correct sea surface stress and vertical mixing coefficients. These wave parameters affect the mixed layer depth and surface flow field of the ocean model by changing the dynamic processes at the sea surface. Ocean current data calculated by the ROMS model are fed back to the SWAN model to correct wave propagation paths and energy distribution. Ocean currents have a particularly significant impact on wave propagation characteristics, especially in areas with strong currents such as tidal currents or typhoon paths.

[0091] S33. The wind field data provided by the WRF model (such as wind speed and direction at a height of 10 meters) directly drive the wave generation and propagation process in the SWAN model, determining the wave intensity and direction. At the same time, the wave parameters in the SWAN model can affect the offshore wind field in the WRF model. Waves affect wind speed, turbulence structure, and momentum exchange by changing the sea surface roughness. During severe storms, high waves increase sea surface roughness, thereby reducing wind speed and significantly affecting wind field distribution.

[0092] 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:

[0093] S41. Using the mesoscale wind field data file output by the mesoscale atmosphere-wave-ocean coupled forecast model, regenerate the wind field for the area around the pier using TurbSim, a full-field turbulent wind simulation software. Create wind speed data for a two-dimensional rectangular grid, including the u, v, and w components of the wind speed, and save it in a format suitable for Fluent input.

[0094] S42. Extract time series data of wave elevation, velocity components (u, v, w), acceleration, and dynamic pressure at specific points from the wave and flow field data of the mesoscale model. These data are used to generate the small-scale wave field data files required by Fluent and configure the corresponding hydrodynamic effects in Fluent. At the same time, extract flow velocity and 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.

[0095] Specifically, the finite element analysis module can accurately simulate and analyze the mechanical response of port pile foundation structures in complex marine environments. By receiving aerodynamic and hydrodynamic input files, the module uses three-dimensional finite element modeling technology to build a highly refined port pile foundation structure model. The Morison equation used in the finite element analysis module to calculate hydrodynamic loads is:

[0096] ,

[0097] Among them, the inertial force is: ,

[0098] The resistance is: ,

[0099] Where, 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.

[0100] 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:

[0101] ,

[0102] 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.

[0103] 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:

[0104] ,

[0105] 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.

[0106] Here, we need to add that:

[0107] 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:

[0108] ,

[0109] in, is the normal stress; is the friction coefficient.

[0110] At the same time, due to the nonlinear characteristics of the material, the pile foundation structure will exhibit nonlinear deformation and stress distribution when subjected to complex loads. To this end, 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:

[0111] ,

[0112] in, σ 1 and σ 3 are the maximum and minimum principal stresses, σ c is the cohesion of the material, ϕ is the internal friction angle.

[0113] In finite element analysis, the effects of various load conditions on pile foundation structures are considered, such as wind loads, wave loads, and dynamic loads caused by flow fields. The superposition of these loads may lead to complex stress distribution and deformation patterns in the pile foundation structure. To this end, the finite element model calculates the stress state and displacement response of the structure under different working conditions using the stress-strain relationship equation. Specifically, the stress-strain relationship is expressed using the following linear elastic equation:

[0114] ,

[0115] in, σ is stress, E is the elastic modulus, ϵ For strain.

[0116] Furthermore, within the Fluent environment, establishing a detailed 3D finite element model of the port pile foundation structure is the foundation of the entire analysis. During the modeling process, the pile foundation's geometric dimensions, material properties, and contact relationship with the surrounding soil must be fully considered. The details are as follows:

[0117] First, the geometric dimensions of the pile foundation need to be accurately input based on actual engineering data, including parameters such as pile diameter, pile length, and burial depth; these parameters directly affect the accuracy of the model and the reliability of the simulation results;

[0118] 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 detail, such as using a piecewise linear model to characterize the elastic-plastic behavior of the material. In addition, the accuracy of the 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.

[0119] 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, through multiple adjustments to the mesh division, boundary conditions, and material parameters, a three-dimensional finite element model that meets the accuracy requirements is finally obtained.

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

[0121] To account for the complex interactions between piles and soil, Fluent software provides a variety of contact algorithms, such as the Coulomb friction model and the stick-slip contact model. These models accurately simulate the effects of friction, lateral earth pressure, and end bearing forces, particularly accounting for vertical and horizontal displacements of the pile and the resulting changes in contact stress. To enhance simulation accuracy, a contact stiffness parameter is incorporated into the contact algorithm. This parameter is closely related to the properties of the pile-soil interface and effectively reflects the contact mechanical behavior under varying soil conditions.

[0122] After completing the model construction and load setting, the Fluent simulation process was initiated. During the simulation, the software calculated the mechanical response of the port pile foundation structure based on the input environmental load conditions, including the stress distribution, displacement response, settlement, and buckling of the pile foundation.

[0123] Furthermore, to ensure the accuracy and convergence of simulation results, the time step and computational parameters must be appropriately set during calculations. Excessively large time steps can distort the results, while too small can significantly increase the computation time or even lead to non-convergence. Therefore, by gradually adjusting the time step and parameter settings, a balance can be found that ensures both accurate and practical simulation results.

[0124] After the simulation calculation is completed, a detailed mechanical response analysis report is generated. Figure 3Figure 2 shows a schematic diagram of the void pressure calculation results at points A and B in the pile foundation, calculated using this embodiment. Specifically, by viewing the pile foundation's stress contours and displacement field distribution diagrams, one can intuitively understand the stress and deformation characteristics of the pile foundation under different load conditions. In particular, one can focus on stress concentration areas and displacement mutation points under extreme weather conditions. Vibration spectrum analysis can be used to evaluate the vibration characteristics and stability of the pile foundation under dynamic loads, focusing on parameters such as resonant frequency and damping ratio to ensure the long-term stability and durability of the structure in marine environments. Based on the simulation results and mechanical response analysis report, targeted structural optimization recommendations can be made. For example, parameters such as the pile foundation diameter, length, and inclination can be adjusted to improve the structure's stress resistance and stability; the contact relationship between the pile foundation and the soil can be strengthened to increase friction and lateral earth pressure. Furthermore, 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 it with field observation data, historical statistical data, and other independent verification sources. Based on the verification results, the model parameters and structural settings are adjusted to improve the model's prediction accuracy and applicability.

[0125] In summary, this application achieves at least the following effects compared to the prior art:

[0126] (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 conditions such as typhoons, including the wind field evolution in the atmospheric boundary layer, the internal dynamic response of the ocean, and the generation and propagation of offshore waves. By fully considering the nonlinear interactions between the various physical fields, the simulation accuracy of the coupled 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 traditional single model methods.

[0127] (2) Advanced modular integrated design. A modular integrated design architecture is adopted, which 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 collaborative 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.

[0128] (3) Parallel computing optimization. In terms of computational 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 simulation precision and spatiotemporal resolution, but also significantly shortens the calculation time of large-scale and complex ocean dynamic processes. This strategy of combining parallel processing with adaptive grids significantly improves the computational performance of the model and meets the forecasting requirements of high real-time and strict accuracy.

[0129] Obviously, the above embodiments are merely examples for clarity of explanation and are not intended to limit the implementation methods. Those skilled in the art will readily appreciate that other variations or modifications based on the above descriptions are possible. It is not necessary and impossible to enumerate all implementation methods here. Obvious variations or modifications arising therefrom remain within the scope of protection of the present invention.

Claims

1. The coastal port pile foundation structure load calculation system based on the coupling model is characterized by: include: A hydrodynamic module, which evaluates the load and mechanical properties of the port pile foundation structure under different environmental conditions; the hydrodynamic module includes at least an atmospheric model module for simulating wind field data in the offshore terminal area, a wave model module for calculating wave field data based on the 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 a finite element analysis module that 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; Among them, 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; ocean current velocity data is used as the fluid dynamic load input; S53. Accurately map ocean model grid data to nodes in 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, start the simulation. The Fluent model calculates the response of the port pile foundation based on the load, including load, settlement, and deformation.

2. The coupled mode-based coastal port pile foundation structure load calculation system 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 and boundary conditions for the atmospheric model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the WRF model, defining key parameters including the size of the simulation domain, time step, physical process options, and output frequency. Once configured, the WRF model simulates atmospheric conditions, including wind speed, direction, temperature, and humidity, to generate high-precision simulation results for the regional wind field. S22. Prepare ocean initial and boundary conditions for the ocean model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the ROMS model, setting physical parameters including the simulation domain, time step, vertical mixing, and bottom drag number. The ROMS model receives wind field data from the WRF model to simulate ocean dynamic processes. The generated simulation data includes ocean current velocity, temperature, and salinity. S23. Prepare initial wave conditions and boundary conditions for the wave model module based on the accuracy requirements for predicting ocean dynamic conditions in the target sea area. Configure the SWAN model to define the wave domain settings, boundary conditions, and wind field inputs. The SWAN model receives wind field data from the WRF model and ocean current data from the ROMS model, and uses the received data as input to simulate wave characteristics in the nearshore area. The generated simulation data includes wave height, period, and propagation direction.

3. The coupled mode-based coastal port pile foundation structure load calculation system according to claim 2 is characterized in that: In the hydrodynamic module, the MCT coupler in the coupler module manages data exchange between the WRF model, the ROMS model, and the SWAN model. The MCT coupler implements bidirectional transmission of physical quantities between the modes by defining grid mapping and communication paths between the modes. Specifically, the following are performed: S31. The bidirectional coupling between the atmospheric model module and the ocean model module is as follows: 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 calculation of sea surface heat flux, momentum flux, and freshwater flux, which influence the dynamic and thermodynamic processes in the ocean model module. The ROMS model feeds back calculated data, including sea surface temperature and sea surface current velocity, to the WRF model, which influences the thermodynamic processes and wind field distribution in the atmospheric model module. S32. The bidirectional 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 to modify the sea surface stress and vertical mixing coefficient. This affects the mixed layer depth and surface flow field of the ocean model module by changing the dynamic processes of the sea surface. The ROMS model feeds back the calculated ocean current data to the SWAN model to modify the wave propagation path and energy distribution, resulting in changes in the wave propagation direction and energy concentration area in the wave model module. S33. The bidirectional 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, determining the intensity and direction of the waves in the ocean wave model module; the wave parameters in the SWAN module affect the offshore wind field of the WRF model, and the waves affect the wind speed, turbulence structure, and momentum exchange in the atmospheric model module by changing the sea surface roughness.

4. The coupled mode-based coastal port pile foundation structure load calculation system according to claim 3 is characterized in that: The finite element analysis module uses the Morison equation to calculate hydrodynamic loads: ， Among them, the inertial 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 drag coefficient.

5. The coupled mode-based coastal port pile foundation structure load calculation system according to claim 4 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.

6. The coupled mode-based coastal port pile foundation structure load calculation system according to claim 5 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.

Images