722 results about "Stress field" patented technology

The invention discloses a tunnel blasting simulation method and system for a karst area, and relates to the field of tunnel blasting simulation, and the method comprises the steps: building a three-dimensional geologic model of the karst area, and generating an initial ground stress distribution rule through numerical simulation; generating refined initial crustal stress field data; calculating stress wave intensity and time sequence data generated by explosion; acquiring stress wave propagation data; calculating the surrounding rock damage degree and the crack initial position and direction; simulating the surrounding rock crack propagation process, and calculating the penetration probability and the propagation speed of the crack near the karst cave in combination with the karst cave characteristic parameters; analyzing the surrounding rock stability and damage evolution trend of the karst cave area; and optimizing explosive source item parameters and stress wave propagation boundary conditions, and generating a karst area tunnel blasting simulation result. According to the method, the influence of the karst cave on stress wave propagation is accurately simulated, powerful support is provided for optimizing a blasting scheme and improving construction safety and efficiency, and the method has important guiding significance on tunnel engineering under complex geological conditions.

A method for quality inspection of an etching paste material based on a spectral response difference, including: identifying a position of density mutation boundary through spectral collection and density gradient analysis. Applying a specific frequency beam to obtain the resonance response parameter, generating stress field reconstruction data. Performing a reverse optical tracking to obtain a defect formation path, generating degradation prediction data. Establishing a self-organizing optical monitoring grid, to form a hierarchical spectral fingerprint library. Constructing a multi-point linked defect blocking network, to form an interconnected optical energy field. Monitoring a cooperative response to obtain network stability data, Identifying abnormal position coordinates based on a material quality grade and the network stability data, obtaining the abnormal position coordinates. Performing a phase adjustment to obtain a phase difference spectral set, conducting a phase comparison marking on the abnormal position coordinates, completing the quality inspection of the etching paste material.

The invention discloses a powder metallurgy gear hardness detection method and system, and belongs to the technical field of metallurgy gears, and the method comprises the steps: obtaining a to-be-detected gear three-dimensional model, and dividing a hardness evaluation region; a thermal response image is collected, a thermal response time gradient parameter TRTG is extracted, and a thermal diffusion characteristic matrix T is formed; establishing a hardness initial model H1 based on the TRTG and compactness; constructing a finite element stress field model, extracting a stress distribution modulation coefficient SDMC, and forming a stress characteristic matrix S; performing weighted correction on H1 based on T and S, and establishing a comprehensive hardness prediction model H2; outputting a hardness prediction result of the whole gear and the key area, and identifying a potential weak area; according to the method, thermal response and stress modulation parameters are fused, non-destructive and high-resolution hardness distribution prediction is achieved, and the method is suitable for quality evaluation and failure early warning of the powder metallurgy gear of a complex structure.

The invention relates to the technical field of material testing, and discloses a coating performance simulation method based on high and low temperature circulation and prestress loading. The method comprises the following steps: constructing a coating matrix composite geometric model containing a microporous structure by acquiring actually measured working condition data of a pipeline or a part, and quantifying surface roughness and interface curvature characteristics in combination with reverse engineering; a dynamic temperature-pressure coupling equation is generated based on the thermal expansion effect of a high-pressure medium, a micropore optical load is converted into equivalent heat flux density through a ray tracing algorithm, transient boundary conditions of a temperature field and a stress field are associated, a crack initiation path is simulated through an extended finite element method, and a life attenuation curve is predicted in combination with a Paris formula. And material parameters are dynamically corrected based on experimental data, and a multi-parameter collaborative optimization model is generated. According to the method, the problems of interface stress evaluation distortion and life prediction deviation under the combined load are solved, and multi-physics field co-simulation support is provided for coating compression resistance optimization.

The invention provides a multi-physics field cooperative regulation and control method in coal mine rock burst prevention and control, and relates to the technical field of coal mine rock burst prevention and control, and the method comprises the steps: synchronously collecting stress field, fracture field and vibration wave field data in real time through a distributed optical fiber sensor network and a micro-seismic monitoring system, and constructing a dynamically updated three-dimensional geomechanical model; based on an adaptive Kalman filtering algorithm, noise reduction multi-source data are fused, a dynamic stress concentration factor and an energy accumulation critical index are extracted, and a rock burst risk probability model is established; a multi-objective optimization algorithm is adopted to generate graded regulation and control instructions of grouting reinforcement, mining speed adjustment and pressure relief drilling, and the graded regulation and control instructions are executed in real time; and iteratively optimizing the model weight by using a transfer learning algorithm in combination with a historical case library to form a closed-loop feedback adaptive prevention and control system. According to the method, the analysis precision of rock mass fracture evolution under complex geological conditions is improved through multi-physics field collaborative perception and dynamic modeling, and risk quantification and real-time regulation and control are realized through multi-algorithm coupling analysis.

The invention discloses a casting detection method and system based on machine vision, and relates to the technical field of intersection of machine vision and material science, and the method comprises the steps: extracting the boundary coordinates of a region with a risk value exceeding a risk threshold according to a defect risk weight map, and forming a dynamic region of interest scanning frame; adjusting the incident angle of a polarized light source according to the track of the dynamic attention area scanning frame, collecting a local image, and obtaining a high-contrast image sequence for eliminating the interference of a surface oxide layer; extracting a defect edge feature map from the high-contrast image sequence through a deformable convolution kernel, and obtaining a defect candidate region mask; the defect candidate region mask is mapped to a von Mises stress field, and a crack germination probability is predicted in combination with a material fatigue limit parameter. According to the method, the accuracy and reliability of detection are improved, and the applicability under complex working conditions is remarkably improved, so that the casting quality control level is effectively improved, and the repair and quality control cost is reduced.

The invention discloses a coal and rock mass fracture seismic source and stress field joint inversion method and system, and the method comprises the following steps: arranging a multi-channel sonic sensor in a coal and rock mass, and collecting original acoustic emission / microseismic signal data; the head wave arrival time / amplitude of the filtered effective waveform data is automatically picked up through an improved AIC (Akaike Information Code), quantitative evaluation is carried out on the quality of a picked-up result, and high-precision positioning is carried out on waveform signals with a high signal-to-noise ratio and a low signal-to-noise ratio; constructing a tension fracture seismic source model based on a displacement discontinuous tensor, obtaining a moment tensor component of a fracture event, and calculating to obtain a spatial orientation, a tension-shear attribute, a deviation angle alpha, energy release and a moment magnitude of the fracture; establishing a stress field inversion model of the tension-shear fracture seismic source, and realizing the inversion of the stress direction and size of the tension-shear fracture seismic source; and a slip direction included angle minimization objective function is established for iteration until convergence conditions are met, and joint iteration inversion of the fracture and the stress field is realized.

The invention relates to a stamping die service life dynamic evaluation system and method based on multi-mode sensing. According to the system, a feature extraction module obtains stress distribution data and crack position information and inputs the data into a trained three-dimensional stress field model to obtain dynamic feature parameters of a stress concentration area under high-frequency cyclic loads; the crack prediction module judges the dynamic characteristic parameters based on a preset threshold value, and once the dynamic characteristic parameters exceed the preset threshold value, analysis is carried out according to the mapping relation between the stress level and the crack rate and the microstructure evolution rule at the same time, and a crack propagation prediction value is generated. And the service life evaluation module inputs the crack propagation prediction value into the trained health state evaluation model, and finally the service life evaluation value of the stamping die is obtained. The service life evaluation value of the stamping die under the existing working condition and the operation condition is accurately estimated, the production efficiency is effectively improved, the production cost is reduced, the potential crack risk of the die is perceived in advance, and accurate and dynamic die service life evaluation is achieved.

The invention discloses a 3D printing path planning method based on electric arc additive anisotropy and a stress field, and the method comprises the steps: constructing a CAD three-dimensional model of a part, and obtaining the stress field of the part through finite element simulation; determining a slice plane; mapping the stress field to a slice plane to form a force flow line; according to the obtained force flow line, rotation transformation regeneration is carried out according to the anisotropy of the used electric arc additive, and a reference trajectory considering the anisotropy of the material is obtained; under the principle of alternate arc starting and extinguishing, the corresponding printing sequence in the layers and between the layers is further planned, and the reference trajectory lines are connected end to end to be planned into a continuous path; and generating a code file, and printing according to the printing sequence. The comprehensive mechanical property of a printed piece is improved.

The present invention discloses a method and apparatus for simulating thermal stress of a casting mold during a service process, and a storage medium. The method comprises: importing a pre-built three-dimensional geometric model of the casting mold, processing the geometric model, and then performing grid division to obtain a casting simulation finite element physical model; assigning material parameters, interface parameters, process parameters, and boundary conditions to the finite element physical model to obtain a casting simulation finite element calculation model; performing casting process simulation calculations to obtain casting simulation results such as mold temperature field and stress field, and exporting result data; interpolating a temperature-stress correction factor; multiplying all correction factor data under the same node numbers with mold stress data to obtain final thermal stress distribution data of the mold during the service process.

The invention discloses a stamping die health state assessment method and system based on digital twinning, and relates to the technical field of die health state assessment, and the method comprises the following steps: a physical sensor network deployed on a stamping die collects die stamping process data in real time; constructing a finite element analysis simulation FEA model to simulate the working condition of the stamping die based on the geometric structure, the material attribute and the stamping process parameters of the die; according to the invention, the data of the die stamping process are collected in real time through the physical sensor network, and real-time calculation is carried out in combination with the finite element analysis simulation model, so that transient stress field, strain field and temperature field data of the die can be output in a short time, and real-time monitoring of the health state of the die is realized; by considering the degradation of the mold material performance along with the use time and the dynamic process of quantitative damage accumulation, the damage accumulation value is accurately calculated through the dynamic material performance database and the continuous damage mechanical model, and the accuracy of the evaluation result is improved.

The invention discloses an electric arc additive manufacturing temperature field and stress field real-time reconstruction method, which comprises the following steps of: 1, constructing a sequential coupling finite element model in which heat conduction analysis and thermoelastic-plastic mechanical analysis are combined on the basis of process parameters of an electric arc additive manufacturing process and three-dimensional geometry of parts, generating a multi-physical field coupling data set among heat flux density, temperature and stress; step 2, constructing and training a neural network model based on a multi-physics field coupling data set, and realizing mapping from heat flux density to temperature and stress; 3, dynamically optimizing heat source model parameters in combination with real-time monitoring data and process parameters in the electric arc additive manufacturing process; and 4, based on the optimized heat source model parameters and the trained neural network model, real-time reconstruction of a temperature field and a stress field in the electric arc additive manufacturing process is achieved. According to the invention, the requirements of rapid sensing and intelligent regulation and control of temperature and stress states on an additive manufacturing site are met.

The invention discloses a tunnel monitoring method and system based on a vibrating wire sensor, and relates to the technical field of tunnel engineering construction safety monitoring, and the method comprises the steps: collecting the vibration frequency data and spatial position coordinates of the vibrating wire sensor in a monitoring region; the method comprises the following steps: calculating a time difference and an energy attenuation ratio between adjacent sensors by extracting spectrum energy density distribution, and constructing a stress wave propagation field; dividing grids in the propagation field to calculate node stress values, generating a stress field distribution cloud picture and extracting a stress concentration area; calculating a gravity center position coordinate change rate and a stress time sequence curve of the stress concentration area, and determining fracture extension parameters; and calculating a surrounding rock bearing capacity attenuation curve based on the parameters, and outputting an early warning signal when the slope exceeds a threshold value. According to the invention, real-time monitoring and early warning of the tunnel surrounding rock instability risk are realized.

The invention discloses a simulation method and system based on multi-field coupling and dynamic grid processing, and belongs to the field of computational simulation, and the method comprises the steps: building and solving a coupling physical field model: carrying out finite element spatial discretization on a solving domain, building and solving a coupling partial differential equation set which at least comprises a temperature field, a displacement field and a stress field, and carrying out the finite element spatial discretization; obtaining numerical solutions of the physical fields in the current state; carrying out progressive unit removal based on a physical criterion; and finally, grid data structure updating and solver synchronization are executed. According to the method, the metal material can be gradually removed when the physical quantity reaches the critical condition, and the grid change is dynamically and robustly processed in the process.

The invention provides a method and a device for evaluating and measuring microscopic mechanical properties of TC4 alloy, relates to the technical field of microscopic mechanical properties, and realizes accurate evaluation of material properties through combination of multi-scale characterization and numerical simulation. Firstly, a standard sample is prepared through selective laser melting, and residual stress field distribution is obtained through X-ray diffraction; dividing a stress field into sub-regions, and establishing a quantitative relationship between microscopic structure parameters and residual stress by combining SEM and EBSD analysis; the method comprises the following steps: acquiring local mechanical properties through a nanoindentation test, and constructing a stress-performance correlation model; and finally, establishing a finite element model, transmitting stress and tissue parameters by adopting trilinear interpolation and an inverse distance weighting method, verifying data reliability through a stress difference criterion, and calculating mechanical properties of each unit. According to the method, multi-scale correlation from microstructure characteristics to macroscopic mechanical properties is innovatively realized, so that the mechanical property distribution of the material is accurately evaluated.

The invention discloses a PCB circuit microdefect real-time detection algorithm optimization method based on artificial intelligence, and relates to the field of printed circuit board monitoring, and the method comprises the steps: collecting multi-physical field data including a temperature field, a stress field, current density distribution and corrosion medium concentration in real time; generating time-space synchronous multi-modal fusion data; quantifying a thermal-mechanical-electric-chemical field interaction relationship through an inter-field coupling coefficient matrix; generating a four-dimensional synergistic effect matrix through tensor fusion; and based on the four-dimensional synergistic effect matrix, in combination with a spatio-temporal evolution diagram output by the carbonization path prediction model, establishing a three-dimensional mapping relationship among the crack size, the carbonization path and the material residual strength, and predicting a residual life prediction value. The method has the advantages that through real-time synchronous fusion of multi-physical field data and tensor coupling modeling, accurate tracking of a crack propagation path and dynamic evaluation of the residual life are achieved, the detection sensitivity is improved, and the life prediction error is reduced.

The invention discloses a landslide construction simulation method based on finite elements, particularly relates to the field of construction simulation, and realizes whole-process dynamic simulation through GTS-NX software. Firstly, a three-dimensional finite element model is established based on geological exploration data, a Mohr-Coulomb elastic-plastic constitutive model is adopted, and the bulk modulus and the shear modulus are calculated through the elastic modulus and the Poisson's ratio. The initial ground stress is simulated by adopting a self-weight stress field method, and the initial displacement influence is eliminated by adopting a displacement difference method. In construction process simulation, the unit passivation / activation technology is adopted to accurately reproduce the excavation and support time sequence, and cyclic calculation of excavation-stress release-support-balance is achieved. In nonlinear analysis, superposition of an elastic stress field and a plastic stress field is considered, and an unloading influence factor is introduced to correct anisotropy of a deformation modulus. A safety coefficient is calculated through a strength reduction method, a multi-working-condition evaluation system is established, a monitoring data fusion and parameter inversion mechanism is developed, and a three-level early warning system is formed.

The invention provides a deep tunnel excavation physical simulation test three-dimensional stress field real-time prediction method, and relates to the technical field of tunnel engineering. Comprising two stages of model construction and training and model use. In the model construction and training stage, a physical constraint reference model is constructed, and training optimization is completed; and in the model use stage, real-time prediction of a stress field of a deep engineering physical simulation test is realized through the trained and optimized model. The method has breakthrough advantages in the aspects of physical law constraint, real-time monitoring data fusion, training loss stability, boundary stress condition accuracy and the like, and is a high-reliability three-dimensional stress field prediction method oriented to a deep complex environment.

The invention relates to the technical field of hydraulic test of nickel alloy pipes, and discloses a hydraulic test system for a thin-wall nickel alloy pipe. According to the system, in the axial loading state of the thin-wall nickel alloy pipe, a closed water pressure cavity is established through an annular sealing clamp, and then strain distribution data on the surface of the pipe body is collected through a distributed optical fiber sensor array. And constructing a three-dimensional stress field reconstruction model based on the acquired strain distribution data, calculating an equivalent plastic deformation threshold value of the pipe body under different water pressure gradients, adjusting a water pressure loading curve slope parameter according to the equivalent plastic deformation threshold value, and generating a dynamic water pressure loading instruction. And the high-frequency pressure pulse generator executes the instruction, synchronously triggers the X-ray diffractometer to collect lattice distortion data of the tube body, inputs the data into the material microscopic damage evolution model, and outputs a microscopic defect spatial distribution map. And finally, material constitutive parameters of the three-dimensional stress field reconstruction model are corrected based on the atlas, closed-loop feedback control is formed, and multi-dimensional data acquisition and accurate management and control of the thin-wall nickel alloy pipe hydrostatic test are achieved.

The invention relates to a strong mining roadway thick-layer transboundary supporting method, and belongs to the technical field of roadway supporting, and the method comprises the steps that geometric parameters, rock mass mechanical parameters and stress field parameters of a strong mining roadway are obtained; based on the geometric parameters, the rock mass mechanical parameters and the stress field parameters, the deformation process and the damage process of the strong mining roadway surrounding rock are simulated, and the incidence relation between each deformation stage of the strong mining roadway surrounding rock and the roof partition damage characteristic is obtained; determining a key stage of irreversible deformation of the strong mining roadway surrounding rock based on the association relationship; determining the critical time of transboundary anchoring of the thick layer based on the key stage; determining the support parameters of the anchor rod and the critical anchoring length of the thick layer transboundary; and supporting the top plate of the strong mining roadway based on the critical time, the supporting parameters and the thick-layer transboundary anchoring length. The method has the effect of improving the stability of the strong mining roadway surrounding rock.

The invention provides a space measurement and analysis method and system for a neutron diffraction rocking curve of a single crystal material, and the method comprises the steps: S1, building a laboratory coordinate system and a sample coordinate system, and setting a test point grid; s2, collecting a two-dimensional diffraction image of each test point to obtain angle information of the test point; s3, integrating the angle information of the test points based on the test point grid, and preprocessing the two-dimensional diffraction image; s4, extracting key structure parameters of each test point based on the preprocessed two-dimensional diffraction image; s5, constructing a stress field according to the key structure parameters; and S6, constructing a distribution map of crystal internal lattice parameters in a three-dimensional space according to the stress field. The method provided by the invention effectively improves the capability of identifying the non-uniformity and local distortion of the crystal structure, can be widely applied to the fields of residual stress evaluation, defect evolution monitoring, material quality control and the like, and has a remarkable engineering application value.

The invention relates to a 3D point cloud-based concrete crack early recognition method. The method comprises the following steps of modeling a fine curvature change generated by a concrete surface at the initial stage of a crack into a high-order local nonlinear offset function; curvature derivative tensor fitting is carried out by using the third-order neighborhood of the point so as to distinguish a real crack sign and noise disturbance; a space stress projection field of the candidate area is formed in the crack germination stage, so that early deformation cracks are identified; based on a stress field guided point cloud reconstruction mechanism, enhancing the point density of the potential crack region identified by the perturbation model by adopting a local reversible interpolation algorithm; according to the interpolation, a crack space consistency constraint graph is constructed according to the continuity and curvature gradient of a surrounding normal, and the structural integrity is improved while the original crack form is reserved; and constructing a crack region into a vector flow field form, and fitting a growth path prediction trajectory of the crack region through the minimum principal curvature direction of crack end points in the point cloud.

The invention relates to the technical field of welding stress detection, and discloses a vehicle welding stress detection method and device, and the method comprises the steps: applying a preset load to a welding part, carrying out the nonlinear finite element analysis and data driving mapping processing of an initial deformation field matrix, carrying out the characteristic decomposition of a stress field tensor to form a characteristic tensor, and carrying out the detection of the welding stress. Carrying out threshold segmentation on stress concentration factor distribution, inputting an abnormal region coordinate set into a damage evaluation model to carry out life and damage quantification processing, carrying out data integration and visualization processing on damage indexes, and generating a welding stress detection report; by applying the preset load, collecting the three-dimensional strain data of the welding part, constructing the initial deformation field matrix and combining nonlinear finite element analysis and data driving mapping processing, the stress field tensor is directly obtained, the error caused by deducing stress distribution only depending on surface deformation is avoided, and the accuracy of stress calculation is improved.

The invention relates to the technical field of engineering mechanical structure optimization, in particular to an excavator structure topological optimization method and system based on real-time stress field reconstruction. The method comprises the following steps of: designing finite element analysis software, modeling and simulating in the finite element analysis software to obtain virtual data, designing a Kalman filter to carry out error correction on the virtual data and actual test data, training the Kalman filter by utilizing the corrected data to obtain prediction of real-time stress field distribution, and further, establishing a topological optimization model of an excavator structure, according to the method, material attributes are coded into a weight matrix in a neural network, and a stress constraint partial differential equation is used as a constraint condition, so that high-precision reconstruction of the stress field of the movable arm of the excavator under an actual working condition is realized, and an accurate boundary condition is provided for topological optimization.

The invention relates to the technical field of metal processing detection, in particular to a copper bar drawing residual stress detection method and system, and the method comprises the steps: collecting production line process state parameters in real time, and mapping the parameters to a predefined multi-dimensional process parameter space to generate a dynamic state vector; adjusting a key coefficient of a material constitutive equation based on the dynamic state vector, and generating a reference stress field; synchronously acquiring multi-physical field data on the surface of the copper bar; constructing a physically constrained machine learning model, and generating a cross-scale fusion stress field; an adjoint optimization equation is constructed based on dynamic process constraints, and model weight parameters are updated through real-time back propagation gradient; and outputting a new cross-scale fusion stress field according to the updated machine learning model, and obtaining related indexes of residual stress detection. According to the invention, the defects of static stiffness and insufficient data fusion of the model in the prior art are effectively overcome, and the real-time detection of the residual stress with adaptive technological parameters is realized.

The invention relates to the technical field of 3D printing, and discloses an environment-friendly packaging box rapid forming method and system based on biodegradable materials, and the method comprises the steps that a thermal-structure coupling finite element model is established; based on the thermal-structural coupling model, a multi-objective optimization model is constructed, and structural design parameters and initial printing path configuration are obtained; training a printing control strategy by utilizing a strategy optimization algorithm; according to the structural design parameters, the printing path configuration and the printing control strategy, a layer-by-layer printing path and a supporting structure configuration scheme matched with the printing object are generated; and the 3D printing process is executed, and closed-loop adjustment is conducted according to the control strategy. According to the method, the thermal-structural coupling finite element model based on the biodegradable material is constructed, and changes of a temperature field, a stress field and a thermal strain field in the printing process are accurately simulated, so that a printing path and supporting structure configuration are optimized, the printing precision is improved, structural instability caused by thermal deformation and interlayer stress is avoided, and the strength and stability of an object are ensured.

The invention discloses a landslide stability evaluation method and system combining a material point method and a transfer coefficient method, and the method comprises the steps: building a three-dimensional material point model for a landslide body based on the material point method, and obtaining the stress field distribution, deformation characteristics and fluid-solid coupling effect of the landslide body; dividing the landslide mass into a plurality of block areas, and determining the block area and the position of each material point; constructing a stress transfer matrix taking the positions and physical information of the material points as elements, calculating normal stress in the vertical direction between the block areas of the landslide mass and tangential stress parallel to the contact surface between the block areas based on the stress transfer matrix, and calculating a real-time stability coefficient; based on the displacement, the pore water pressure and stress and the real-time stability coefficient, the dynamic safety coefficient of the landslide mass is calculated, and stability evaluation is conducted on the landslide mass. Through combination of the material point method and the transfer coefficient method, the mechanical behavior and the material propagation process of the landslide mass are considered, and the evaluation precision of the landslide stability is improved.

The invention relates to a dynamic grid-based foundation pit temperature field-seepage-stress field coupling calculation method, which comprises the following steps of: establishing a nonlinear dynamic bidirectional coupling constitutive relationship of a temperature field, a seepage field and a stress field containing a two-dimensional foundation pit cofferdam, and obtaining a three-field full coupling equation set; in combination with a GPU acceleration solver, solving the three-field full-coupling equation set by adopting a sub-domain decoupling-asynchronous iteration strategy; and performing real-time early warning on the cofferdam construction period according to a solving result. The method has the advantages that compared with a traditional linear model, the precision of the established three-field full-coupling equation set is improved, and the interaction influence relation of three fields can be well simulated.

The invention relates to the technical field of data processing, and discloses a temperature effect modeling method and system for a concrete tower drum. The method comprises the steps of establishing a finite element model considering joint grading thermal resistance and prestress correction material parameters, adaptively adjusting a time step length according to a solar altitude change rate to calculate a radiation load, obtaining a convergence temperature field through bidirectional coupling iteration of a temperature field and a stress field, and extracting characteristic temperature to train a neural network to realize rapid prediction. According to the invention, the precision and efficiency of temperature effect analysis of the fabricated concrete tower drum are improved.

The invention belongs to the technical field of tunnel engineering, and particularly relates to a cold region tunnel water thermal coupling numerical modeling method and system. The method comprises the following steps: establishing an equal-proportion tunnel geometric model; defining equation variables, and setting various parameters, boundaries and initial conditions of materials; performing grid division to obtain a tunnel numerical calculation model; based on the tunnel numerical calculation model, solving a pre-constructed tunnel hydrothermal-thermal coupling mathematical model to obtain tunnel temperature field, moisture field and stress field distribution data; the tunnel water thermal-mechanical coupling mathematical model comprises a temperature field control equation, a moisture field control equation and a stress field control equation; and visualizing the distribution data of the temperature field, the moisture field and the stress field of the tunnel to obtain the change conditions of the temperature field, the moisture field and the stress field of the tunnel in winter. The method can effectively simulate and analyze the distribution conditions of moisture, temperature and stress of the tunnel in the cold region under the multi-field coupling condition, and provides scientific basis for design and maintenance of tunnel engineering in the cold region.