33 results about "Finite element solution" patented technology

The invention relates to the technical field of electromagnetic scattering, in particular to a frequency domain response calculation method and device of an electromagnetic field, equipment and a storage medium. Performing recursive transformation on the time domain basic system equation according to the discrete time domain set and a backward Euler method to obtain a time domain transformation equation and an algebraic matrix equation; carrying out back substitution solving operation on the algebraic matrix equation according to a time domain stepping method and a time domain transformation equation; performing full-time-domain response calculation on the time-domain finite element solution based on a time-domain finite element method to obtain a time-domain base solution; calculating the time domain base solution according to a preset linear superposition principle to obtain a frequency domain response; a continuous problem is dispersed through a backward Euler method, rapid back substitution solution is carried out after factorization for one time, and redundant calculation is reduced; based on a linear superposition principle, phase weighting is carried out to synthesize a time domain base solution, frequency domain response of any frequency is rapidly obtained, frequency-by-frequency repeated calculation is avoided, efficiency is improved, and the method is suitable for broadband high-precision electromagnetic field analysis and has remarkable application value.

The invention belongs to the technical field of pressure vessel structure analysis, and particularly relates to an intelligent reduced basis finite element method for rapidly analyzing the deformation behavior of a pressure vessel. According to the method, a finite element method is improved by adopting a reduced basis method, an online and offline combined solution mode is developed, after design parameters such as a structure size, a load type and size, material attributes and the like are preset, a reduced basis space is obtained by adopting an offline solution mode, and then a mapping relation between a reduced basis solution and a finite element solution is established; according to the method, order reduction of a finite element model is achieved, the degree of freedom of a problem is reduced, the calculation amount needed by solving a linear algebraic equation set through a traditional finite element method is greatly reduced while the precision is guaranteed, a Gaussian process regression method is further introduced to describe the elastic-plastic deformation behavior of the pressure vessel, and therefore the deformation behavior of the pressure vessel is rapidly solved.

This invention belongs to the field of numerical electromagnetic field solution, specifically a mesh adaptive method based on reconstructive posterior error estimation. The invention first performs finite element modeling of the target microwave component, introducing boundary conditions and excitations to establish a corresponding electromagnetic simulation model. Then, a tetrahedral mesh is used to divide the solution domain, and finite element electromagnetic simulation analysis is performed on the model to calculate the electric field at the centroid of the four faces of the tetrahedral mesh. The least squares principle is then used to refit the electric field on the tetrahedral mesh to reconstruct a first-order linear fitting expression. Next, the reconstructed solution at the nodes is calculated using the reconstructed fitting expression, and posterior error estimation is performed based on the electric field at the nodes. Finally, the mesh is refined according to the error on the tetrahedral mesh, and finite element electromagnetic simulation analysis is performed again until the finite element solution meets the accuracy requirements. This invention has the advantages of simple engineering implementation, low computational cost, and high universality.

The application discloses a kind of T-shaped sheet piece angle joint welding's welding temperature field and deformation prediction method, it is related to welding temperature field and deformation simulation prediction technical field, method steps include obtaining T-shaped sheet piece angle joint welding's temperature field theoretical solution;T-shaped sheet piece grid division and heat source setting are carried out;Solving the temperature field of T-shaped sheet piece angle joint welding, and comparing and analyzing iterative finite element analysis setting with the temperature field theoretical solution obtained;Obtain the welding deformation and stress field of T-shaped sheet piece.The application solves the temperature field of T-shaped sheet piece angle joint welding by using finite element analysis method, and compares and analyzes iterative finite element analysis setting with the temperature field theoretical solution, improve the precision of finite element solution.

The invention relates to the technical field of building curtain wall design, in particular to a parameterized finite element simulation calculation method of a point-supported glass curtain wall for improving design efficiency and a storage medium, and the method comprises the following steps: defining an input parameter set, and setting the input parameter set as a global variable; obtaining a segmentation geometric model; unit and material characteristics are given to the segmentation geometric model; performing grid division on the segmented geometric model to generate a finite element grid, so that corresponding finite element nodes exist at the positions of the geometric key points; finite element nodes are recognized, displacement constraint is applied, and meanwhile evenly-distributed pressure loads are applied; mechanical response data are obtained, and when any variable is changed, modeling calculation is repeated based on the changed parameters; and extracting stress distribution, deformation and support counterforce data from a finite element solution calculation result. According to the invention, parameterized automatic simulation is realized, and stress concentration is accurately captured; the utilization of computing resources is optimized while the large deformation behavior of the brittle glass material is accurately simulated, and the design efficiency and the structural safety are improved.

The invention discloses an efficient simulation method for collision between long ellipsoid particles and a wall surface based on discrete elements, and belongs to the technical field of gas-solid two-phase flow numerical simulation. The method comprises the steps that a universal long ellipsoid particle-wall surface collision elastic force model is established, and the model represents a universal function relation between elastic force and particle geometry, material attributes and kinematics state parameters; based on the elastic force model and the momentum theorem, constructing a particle-wall surface collision dynamic model; based on the kinetic equation, a discrete time method and a numerical integration algorithm are adopted, and force, speed and energy parameters, changing along with time, of the long ellipsoid particles and the wall face in the whole collision period are solved. According to the method, the analytic model is embedded into a discrete element simulation framework to replace calculation-intensive finite element solution, so that second-level efficient solution is realized while the calculation precision is ensured, and the calculation efficiency is remarkably improved.

The invention belongs to the technical field of fuel tank simulation analysis, and discloses a fuel tank lightning environment direct effect analysis method and system, and the method comprises the following steps: building a fuel tank simulation model, and carrying out the material attribute setting of the fuel tank simulation model; finite element solution analysis is conducted on the fuel tank simulation model with the set material attributes, and a simulation analysis result is generated; and performing comparison verification based on the simulation analysis result and the actual lightning stroke test data of the fuel tank. The method solves the problem that the aircraft fuel tank lacks effective simulation analysis means under the requirements of lightning environment airworthiness certification terms and laws and regulations, different lightning stroke scenes can be quickly covered, and the analysis period is greatly shortened.

The present invention provides a method and device for optimizing the head shape and water entry strategy of a water-air cross-media robot, including determining a design domain; performing Latin hypercube sampling within the design domain to obtain a design parameter vector; drawing a NURBS curve to form a water-air cross-media robot model; resetting parameters; generating a simulation domain and meshing it; performing a finite element solution on the current mesh to extract solution data; constructing low-quality and high-quality data sets; obtaining the optimal hyperparameters of a proxy model through Bayesian optimization; training the proxy model based on the optimal hyperparameters; constructing an optimization target and employing a multi-objective optimization algorithm to obtain a solution that satisfies the biomimetic morphological and component constraints under the optimization target; determining whether the proxy model has converged and outputting the optimal design parameters. By constructing a fully automated design paradigm for parametric modeling, numerical simulation, and parameter optimization, the present invention avoids manual real-time intervention in the intermediate process, significantly improving process efficiency.

The invention discloses an intelligent bearing fault monitoring method and device for wind power operation and maintenance, and belongs to the technical field of wind power fault monitoring. A joint simulation analysis mode based on HyperMesh and ANSYS is introduced, and refined modeling and grid division are performed on a fan bearing structure, so that the bearing fault monitoring efficiency is improved; stress and vibration response characteristics of the bearing under actual working conditions are analyzed by combining finite element solution, bearing vibration energy distribution and key response areas are accurately reflected, and effective coverage of axial, radial and transverse vibration information is realized; according to the method, the arranged sensors can more effectively acquire vibration information closely related to the running state of the bearing, and compared with a unified or empirical measuring point arrangement mode in the prior art, the signal-to-noise ratio and identifiability of bearing fault characteristic signals can be improved, and the probability of occurrence of missing detection and misjudgment can be reduced.

The application belongs to the technical field of pressure vessel structure analysis, and specifically relates to an intelligent reduced basis finite element method for fast analysis of pressure vessel deformation behavior. The application improves the finite element method by using the reduced basis method, and develops an online and offline combined solving mode. After design parameters such as structure size, load type and size, material properties and the like are preset, the reduced basis space is obtained by using the offline solving mode, and then a mapping relationship between the reduced basis solution and the finite element solution is established, so that the finite element model is reduced, the degree of freedom of the problem is reduced, the calculation amount required for solving the linear algebraic equation set by the traditional finite element method is greatly reduced while the accuracy is ensured, and further, the Gaussian process regression method is introduced to describe the elastic-plastic deformation behavior of the pressure vessel, so that the fast solution of the pressure vessel deformation behavior is realized.

The invention discloses a pungent-keeping finite element solving method for dynamic response of a structure under a Hamiltonian system. According to the method, a finite element method is applied to a Hamiltonian equation, a Pocking finite element method is deduced, and it is proved that the Pocking finite element method is of a strict Pocking structure. And on the other hand, programmed solution of a high-order iteration format is realized by utilizing mathematical calculation software. In addition, the flexibility of the pungency-preserving finite element method to the linear Hamiltonian system is also discussed. A plurality of numerical examples verify the excellent performance of the Pocking finite element method in the structural dynamic response problem, and prove the capability of keeping energy conservation and excellent stability of the Pocking finite element method in the long-term dynamic simulation process, and compared with a traditional numerical method, the Pocking finite element method has remarkable advantages.

The embodiment of the invention provides a simulation method and device for craniocerebral growth and a storage medium, and belongs to the technical field of images. The method comprises the following steps: constructing three-dimensional brain models in one-to-one correspondence with standard brain CT sample images corresponding to different growth time periods; taking the three-dimensional brain model corresponding to an initial time period in different growth time periods as an initial simulation model for segmentation to obtain a plurality of brain model components; performing finite element solution of growth parameters on a gridding three-dimensional craniocerebral model obtained by gridding each craniocerebral model component to obtain a predicted craniocerebral model of the next growth time period; obtaining a first target growth parameter corresponding to the current growth time period according to the predicted craniocerebral model and the three-dimensional craniocerebral model in the same growth time period; and taking the predicted craniocerebral model as an initial simulation model to obtain a first target growth parameter of the next growth time period. According to the embodiment of the invention, the simulation efficiency and simulation precision of skull simulation can be considered.

The present application belongs to the technical field of welded structure reliability evaluation, and particularly relates to a calculation method of stress intensity factor of gap defect weld under tensile shear load, comprising: selecting an automobile sheet part as a research object, and constructing a strip two-dimensional model of a gap defect weld containing an initial gap; based on the linear elastic superposition principle, the load borne by the gap defect weld is decomposed into multiple load types; based on the strip two-dimensional model of the gap defect weld containing the initial gap, the analytical solution of the stress intensity factor under the load generation mode corresponding to different load types is calculated respectively, and the analytical solution of the stress intensity factor under different load generation modes is calculated; a finite element model of the gap defect weld is constructed and finite element analysis is carried out, and the finite element solution of the stress intensity factor under different load generation modes is calculated; the analytical solution and the finite element solution of the stress intensity factor under different load generation modes are compared, and the accuracy of the analytical solution of the stress intensity factor is verified.

The present invention relates to the fields of image processing and biomechanics, and specifically is a wrist joint mechanical detection method based on finite element analysis, including: acquiring medical imaging data, image preprocessing, assigning material properties, setting mechanical loads and applying boundary conditions, finite element solution and result analysis. Through finite element analysis, the present invention can accurately simulate the mechanical behavior of the wrist joint under different loads, and comprehensively analyze its stress and strain distribution. The present invention can not only detect the stress conditions of the bones, but also conduct a detailed analysis of the mechanical response of the cartilage tissue, providing richer and more accurate information for clinical diagnosis, and can also perform separate segmentation and assignment processing for the lesions. The present invention fully considers the impact of the lesion tissue on the mechanical properties of the wrist joint, conducts personalized analysis for patients, and can formulate a more accurate treatment plan for each patient to improve the treatment effect.

The application belongs to the technical field of three-dimensional electromagnetic radiation numerical solution, and relates to an electromagnetic radiation simulation analysis method based on local grid self-adaptive division. Firstly, the application carries out finite element modeling on a target antenna, introduces PML boundary and wave port excitation to establish a corresponding finite element calculation domain model, and adopts tetrahedral grid division to classify the calculation domain model; then finite element electromagnetic radiation simulation analysis is carried out to obtain finite element solution, and error estimation and encryption strategy are used to give the encryption size of the grid nodes in the non-PML region; finally, the grid is encrypted, finite element electromagnetic radiation simulation analysis is carried out again, and the process is repeated until the overall solution error is less than the set error limit, the final calculation grid is obtained, and finite element electromagnetic radiation simulation analysis is carried out. The application only encrypts the tetrahedral grid in the non-PML calculation domain, ensures the original characteristics of the tetrahedral grid in the PML calculation domain, and realizes the rapid simulation design of the antenna radiation characteristics while ensuring the calculation accuracy.

The present application relates to a kind of compressor blade grid parameterization joint call analysis method.The present application is with MATLAB software as total control hub, by the flow control and data interaction ability of MATLAB, the grid parameterization deformation technology of HyperMesh software and the parametric finite element solution technology of ANSYS software are organically integrated together, construct a full-automatic closed loop analysis method and system from variable input to result output.The method realizes the quick response and batch calculation of several kinds of uncertain variables such as geometry, material, load, and is especially suitable for the fatigue life reliability analysis of compressor blade.

The present application relates to the technical field of building curtain wall design, and more particularly to a parameterized finite element simulation calculation method of point support glass curtain wall for improving design efficiency and a storage medium, comprising: defining an input parameter set and setting it as a global variable; obtaining a segmented geometric model; giving the segmented geometric model unit and material properties; performing mesh division on the segmented geometric model to generate finite element meshes, so that corresponding finite element nodes exist at geometric key point positions; identifying the finite element nodes and applying displacement constraints, while applying uniform pressure loads; obtaining mechanical response data, when any variable is changed, repeating modeling calculation based on the changed parameters; extracting stress distribution, deformation and support reaction data from the finite element solution calculation results. The present application realizes parameterized automatic simulation, accurately captures stress concentration; optimizes the use of computing resources while accurately simulating the large deformation behavior of glass brittle materials, and improves design efficiency and structural safety.

This application provides a nonlinear finite element accelerated solution system and method. The system includes: a physical information neural network model construction module configured to: acquire configuration points in the solution domain of the nonlinear partial differential equation to be solved; input the configuration points into a fully connected neural network to obtain a predicted solution; and obtain a physical information neural network model based on the predicted solution; an input module configured to: construct a finite element model of the nonlinear partial differential equation to be solved; transform the nonlinear partial differential equation to be solved into a system of nonlinear algebraic equations; and generate a predicted value vector based on the finite element model; and a solution module configured to: solve the system of nonlinear algebraic equations using a nonlinear solver to obtain the target solution. This addresses the problem that in current finite element methods for parameterized research scenarios, each time a new parameter value is input, a complete and computationally expensive nonlinear finite element solution process needs to be started independently, leading to repetitive calculations and low efficiency.

The invention discloses an in-situ homogenization method and system for an electronic packaging structure and a medium, and belongs to the technical field of microelectronic integration. The method comprises the following steps: carrying out layering, pixelation and unitization preprocessing on a packaging structure, and dividing macroscopic units; an in-situ structure unit is constructed through in-layer boundary expansion and out-layer boundary expansion, and a boundary overrun problem is processed by combining a low-modulus false expansion region; applying multiple groups of displacement / thermal field boundary conditions, and extracting an equivalent stiffness matrix and an anisotropic thermal expansion coefficient of the macroscopic unit through finite element solution; and finally, replacing the macroscopic unit with an equivalent homogeneous material, and constructing a packaging macroscopic model to support high-efficiency-thermal simulation. According to the method, the cross-layer and cross-scale material interaction effect is truly simulated through in-situ boundary expansion, the limitation of periodic boundary condition assumption in traditional homogenization is broken through, the equivalent parameter calculation precision is remarkably improved, and an efficient solution is provided for reliability evaluation of complex electronic packaging.

The invention discloses an optimal design method and system for an offshore wind power floating foundation structure, and belongs to the field of ocean engineering structure design, and the method comprises the following steps: S1, building a three-dimensional finite element model of the floating foundation structure, and defining a design domain, a non-design domain, material attributes, volume ratio constraints and limit load working conditions for the three-dimensional finite element model; and S2, independently performing finite element solution for each limit load working condition, and calculating the flexibility sensitivity of each unit in the design domain, so that the overall performance and robustness of the final topological configuration are remarkably improved, the sensitivity of a high-conflict region is actively inhibited, and the reliability of the topological configuration is improved. The optimization process can avoid invalid swing between mutually contradictory rigidity requirements, so that a real global optimal solution can be explored and converged, unprecedented design insight and controllability are provided for design engineers, and convergence stagnation or divergence caused by target conflicts is effectively avoided.

The application discloses a design method of a heterogeneous superstructure material with zero Poisson's ratio and high bearing capacity. In the linear elasticity theory framework, first, two kinds of cell configurations with theoretically positive and negative Poisson's ratio values are selected based on a structural Poisson's ratio calculation formula of a classical beam theory; then, finite element solutions of the Poisson's ratio and the Young's modulus of the selected cells are calculated; and finally, the two kinds of cells with positive and negative Poisson's ratios are assembled in a column and row arrangement mode according to a given layout rule to obtain a heterogeneous lattice superstructure material with zero Poisson's ratio and high bearing capacity. The application can obtain a heterogeneous superstructure material with zero Poisson's ratio and high bearing capacity through cell assembly design. The heterogeneous lattice superstructure material formed by assembly has stronger deformation resistance and one-way deformation capacity of the zero Poisson's ratio structure, and can provide scheme support for actual engineering design.

The present application relates to the technical field of electromagnetic scattering, and particularly relates to a method and device for calculating frequency domain response of electromagnetic field, equipment and storage medium; a time domain basic system equation is recursively transformed according to a discrete time domain set and a backward Euler method to obtain a time domain transformation equation and an algebraic matrix equation; the algebraic matrix equation is back-substituted and solved according to a time domain stepping method and the time domain transformation equation; a time domain finite element solution is calculated based on a time domain finite element method to obtain a time domain base solution; the time domain base solution is calculated according to a preset linear superposition principle to obtain a frequency domain response; a continuous problem is discretized by the backward Euler method, and a fast back-substitution solution is obtained after one factor decomposition to reduce redundant calculation; the time domain base solution is phase-weighted and synthesized based on the linear superposition principle to quickly obtain an arbitrary frequency domain response, avoid repeated calculation for each frequency, improve efficiency, be suitable for wideband high-precision electromagnetic field analysis, and have significant application value.

A wallboard assembly performance evaluation method based on actual measurement data driving is characterized by comprising the steps that actual measurement data of key loading points of an assembly site are obtained firstly, and a real wallboard finite element model is generated on the premise that the volume mesh topology is kept unchanged based on actual measurement deviation; then carrying out assembly condition analysis and extracting performance indexes, and further constructing a training sample to realize online performance evaluation based on AI; and finally, returning a finite element analysis result to the twinborn visualization platform for state mapping and display. According to the method, a closed loop of rapid generation of a real wallboard, stable finite element solution, automatic construction of a training sample and twinborn visual display is realized.

The invention belongs to the technical field of reliability evaluation of welding structures, and particularly relates to a calculation method of a stress intensity factor of a gap defect welding spot under tension-shear load, which comprises the following steps: selecting an automobile thin plate as a research object, and constructing a strip-shaped two-dimensional model of the gap defect welding spot containing an initial gap; based on the linear elastic superposition principle, the load borne by the gap defect welding spot is decomposed into various load types; based on the strip-shaped two-dimensional model of the gap defect welding spot containing the initial gap, respectively calculating analytic solutions of stress intensity factors in load generation modes corresponding to different load types, and calculating analytic solutions of stress intensity factors in different load generation modes; constructing a finite element model of the gap defect welding spots, performing finite element analysis, and calculating finite element solutions of stress intensity factors in different load generation modes; and comparing the analytic solutions and the finite element solutions of the stress intensity factors in different load generation modes, and verifying the accuracy of the analytic solutions of the stress intensity factors.

The application provides a porous medium fluid-solid coupling iterative acceleration calculation method based on a Fourier neural operator, and relates to the technical field of porous medium fluid-solid coupling simulation. The method comprises the following steps: constructing a sample data set containing heterogeneous parameters and boundary conditions, and obtaining a pore pressure field and a volume strain field true value by using a traditional numerical simulation; training a Fourier neural operator mechanical proxy model based on the data set, wherein the Fourier neural operator mechanical proxy model takes a pressure field, a Young's modulus field, a Poisson's ratio and a boundary load as input, and takes a volume strain field as output. Subsequently, a mixed iterative calculation framework based on a constant stress splitting is constructed: a fluid module adopts a finite volume method to solve a pressure field, a solid module calls the trained proxy model to quickly infer a strain field, and a two-way coupling iteration is realized through a porosity evolution model until convergence. The method replaces the traditional finite element solution by using the proxy model, and significantly reduces the solution cost of a large linear system under the premise of maintaining physical consistency and calculation accuracy.

The invention relates to the technical field of data processing, in particular to a novel arbitrary polygon unit calculation method for structure simulation analysis, which comprises the following steps of: acquiring a geometric model, extracting node coordinate data, calculating coordinates of a central point and an edge midpoint of a polygon unit to form an extended coordinate data set; a quadrilateral subunit topological structure is generated through geometric connection operation, a strain displacement relation is established, equivalent strain matrix data is calculated through constraint equation fusion, numerical integration and stability correction are carried out by adopting a single-point integration algorithm and an hourglass control algorithm, and rigidity matrix data controlled by the hourglass are output. Integrating to a finite element solution system for numerical solution, and calculating stress-strain distribution based on node displacement data. According to the method, the problems of low efficiency and poor precision caused by dependence on subunit decomposition in the prior art are solved, and the calculation stability and the engineering applicability are improved.

This invention belongs to the field of electromagnetic simulation calculation and provides an electrically large-scale electromagnetic simulation method based on a hybrid algorithm for partitioned solution. First, the large-scale electromagnetic target is divided into solution regions according to component structure. For different algorithm solution domains, a domain decomposition method is used to decompose the hybrid algorithm solution matrix into subdomains, obtaining the hybrid algorithm solution matrix equation after domain decomposition. Then, for the hybrid algorithm solution matrix after domain decomposition, an outermost preprocessor is constructed using a block Jacobi preprocessing method, and an inner preprocessor is constructed using a multi-wavefront block incomplete decomposition method. Finally, multi-threaded parallel technology is used to iteratively solve the hybrid algorithm solution matrix equation, obtaining solutions for both the finite element solution domain and the integral equation solution domain. Based on this, this invention has advantages such as significant multi-algorithm fusion, high preprocessing efficiency, small memory footprint, and good parallel scalability, enabling efficient electromagnetic simulation analysis of electrically large-scale targets.

The invention discloses an ecological bank protection gabion finite element analysis and structure improvement method, and relates to the field of water conservancy and geotechnical engineering.The ecological bank protection gabion finite element analysis and structure improvement method comprises the steps that a gabion structure is divided into a flexible net cage layer, a contact coupling layer and a stone filling layer, the stone filling layer adopts equivalent particle domain unit modeling, and the contact coupling layer is used for establishing a mechanical coupling relationship between the two layers; normal contact rigidity, tangential rigidity and friction coefficient parameters are set in the contact coupling layer so as to describe the top pressure, the friction reaction force and the slippage behavior of the stone on the mesh cage; applying a dynamic water impact load with a time pulse characteristic to the outer surface of the flexible mesh cage; and supporting constraint is applied to the bottom of the model, transient finite element solution is carried out, and dynamic stress distribution of each layer of structure is obtained. The problem that discontinuous stress and pulse impact response of the gabion cannot be accurately simulated through traditional finite element analysis is solved.

The present application relates to the technical field of building engineering construction, in particular to a kind of assembly type component hoisting process simulation method and system based on BIM. By obtaining component model and reinforcement geometric coordinates containing virtual hook node, embedded lifting ring node and rigging;According to the position of hook and lifting ring node to determine the geometric length of rigging, nonlinear finite element solution is carried out combined with rigging length increment sequence under different working conditions, and sleeve hole deformation coordinates under each working condition are obtained;According to the position difference between deformation coordinates and corresponding reinforcement coordinates, determine the hole group deviation distribution set, extract the deviation distribution radius and distribution center offset vector;Determine the optimal working condition based on the numerical difference of radius, and generate control instructions according to the offset vector of optimal working condition and rigging length increment sequence.The present application realizes the quantitative decoupling of shape error and position error, improves the on-site alignment adjustment precision through step-by-step control of shape correction first and then alignment.

The invention provides an ABAQUS-based energy pile thermosetting coupling numerical simulation method, and relates to the technical field of energy geotechnical engineering.The method comprises the steps that a three-dimensional geometric model of a soil body, a pile body and a heat exchange pipe is constructed; the soil body is partitioned, and the pile body and the heat exchange tube are subjected to Boolean operation to generate an integrated pile body component with the common joint of an interface; assembling to form a three-dimensional solid model and defining material parameters; performing grid discretization by adopting a temperature-displacement coupling unit; establishing a multi-stage analysis step sequence, wherein the multi-stage analysis step sequence comprises two crustal stress balance and temperature-displacement coupling analysis steps; applying load and boundary conditions; and finite element solution is carried out, and temperature field, thermal stress and pile-soil interface mechanical response data are extracted. According to the method, the whole energy pile construction and service process can be truly reflected, accurate coupling calculation of the temperature field and the stress field is achieved, and a reliable basis is provided for energy pile design optimization.