77 results about "Finite element program" patented technology

The invention relates to a probabilistic finite element method (PFEM)-based steel-bride fatigue reliability evaluation method, in particular to a fatigue reliability evaluation method by combining dynamic weighing data (WIM) of a car and the analysis of the PFEM, which aims at solving the problems that the bridge structure health monitoring cost is higher and the stress measuring point distribution is limited, and utilizes the real WIM data to statistically analyze car load and establishes a probabilistic distribution model of a car type, road distribution, axle weight and wheel-base. Finite element method analysis software is used for establishing a numeric model of a bridge. A probabilistic finite element program is compiled, stress amplitude and stress cyclic quantity of a key part under the random car load effect are extracted through the sampling, loading and finite element analysis to be statistically counted. The fatigue reliability calculation is performed according to the statistical result, and the reliability attenuation trend of the evaluated part in the subsequent application process is predicted. By adopting the method, a feasible analysis tool is provided for the fine fatigue evaluation of complicated steel-bridge details.

A method for determining the time-varying reliability of prestressed concrete box girder bridges based on the stochastic finite element method. For prestressed concrete box girders that are widely used in highway bridge systems and have common early-stage diseases, by combining degenerate shell elements, material constitutive relations and The crack model simulates the nonlinear characteristics of the box girder in the cracking stage, and the embedded concrete creep, shrinkage and steel strand stress relaxation model considers the deformation, prestress loss and internal force redistribution of the box girder during its entire life, through the reinforcement The uniform corrosion model and the pitting corrosion model calculate the section loss of steel bars (steel strands) in carbonization environment and chloride ion environment respectively. Finally, through the definition of random variables, sampling, calling finite element program calculation and extracting stress and strain results, the reliability calculation of the box girder bridge is realized and used to guide the subsequent maintenance and reinforcement work.

The invention discloses a testing method of I-type crack dynamic fracture toughness under explosive load. Detonators and a cement mortar sample with a crack are utilized to perform bursting testing research, the dynamic fracture toughness of the sample is determined by an experiment-numerical method, the sample bearing load and the crack initiation time are determined by a strain signal obtained through testing, the obtained time travel curve is input into a finite element program Ansys, the near field displacement of a crack tip is calculated by a 1 / 4 node unit, then the time travel curve of the I-type dynamic fracture stress intensity factor of the sample is obtained by a displacement extrapolation method, and the dynamic fracture toughness of the material is the stress intensity factor at the crack initiation moment, therefore, the new testing method of the I-type crack dynamic fracture toughness under explosive load is given. The method is simple, the operation is convenient, and the testing method has the practical application significance to the research of the dynamic fracture character of the rock under explosive load.

The invention relates to a finite-element analysis-based method for designing a profile of an autoclave molding fixture of a composite material member, belonging to the technical field of autoclave molding fixtures of composite material members. The finite-element analysis-based method comprises the steps of: firstly, using a member design profile as an initial profile of the fixture; developing a finite element program for realizing curing deformation analysis of the composite material member to obtain a profile of the composite material member after cured and deformed; under a global coordinate system, calculating a displacement L of each node, judging whether a size precision of the composite material member meets the requirement, if yes, completing the design of the profile of the fixture; if not, compensating a displacement of the node of the composite material member along the direction of a coordinate axis to the node of the profile of the fixture to obtain a new profile of the fixture; and repeating the processes of analyzing and judging a finite element until the profile of the fixture, meeting the precision requirement, is obtained. According to the finite-element analysis-based method, the problems that the profile of the fixture is directly extracted from the profile of the member and then the deformed member is repeatedly repaired are solved, the member forming precision is improved, and the development cost is reduced.

The invention provides a finite element analysis method for thermal-coupled vibration characteristics of an impeller structure. The technical scheme is that the finite element analysis method comprises the following steps: establishing a three-dimensional model of the impeller structure by utilizing a solidworks parameter modelling method; importing the three-dimensional model of the impeller structure into finite element analysis software ANSYS, setting the attributes, parameters, loads and constraints of materials according to the actual working conditions in an ANSYS software pre-processing module, and dividing finite element grids; calculating the first five orders of fixed frequencies of the impeller structure at different rotation speeds through an ANSYS software modal analysis module; calculating the first five orders of fixed frequencies of the impeller structure at different rotation speeds in high-temperature environment through considering the influences of a temperature field; comparing the fixed frequencies in above analysis to obtain the influences, caused by the temperature field, on the vibration characteristics of the impeller structure, so as to provide theoretical basis for the failure mechanism research of the impeller structure under the high-temperature environment.

The invention relates to a comprehensive analysis method for steam explosion of nuclear reactors. The method mainly comprises the following steps of: 1, obtaining initial parameters of the steam explosion; 2, calculating the change of melt in the steam explosion process; 3, calculating the energy exchange on a phase interface; 4, calculating the mass transfer on the phase interface; 5, calculatingthe momentum exchange on the phase interface; 6, semi-implicitly solving energy and momentum equations; 7, full-implicitly solving a mass conservation equation and outputting distribution of pressurewaves; 8, carrying out geometric modeling on main components by using a finite element program ABAQUS; 9, creating material physical properties and analyzing and inputting by taking the pressure wavedistribution in the step 7 as a load; 10, carrying out grid dividing on the main components established in the step 8; and 11, submitting analysis and outputting the stress change of the main components. The comprehensive analysis method disclosed by the invention has the advantages that the distribution situation of the pressure waves and the impact effect of the pressure waves on the main components under steam explosion accidents of the nuclear reactors can be comprehensively analyzed, so that bases are provided for more comprehensively and effectively evaluating the safety of the reactors.

The invention, a testing method for mode I crack dynamic fracture overall process parameters under blast loading, discloses the testing method for the dynamic crack initiation toughness, dynamic expansion toughness and dynamic crack arrest toughness of mode I crack under blast loading, and belongs to the technical field of geotechnical engineering. The configuration of an internal single crack disc specimen with a large size is an autonomous design; testing results cannot be influenced by reflected extension waves and transmission waves; and the method is suitable for researching the crack growth rate and growth behavior of the whole process. Blast loading adopts an 8# detonator; numeral calculation adopts a finite element difference software AUTODYN; and a result basically accords with atest. A load stress time history curve can be input into a finite element program ABAQUS; the time history curve of a dynamic stress intensity factor can be obtained; and the dynamic fracture toughness of a material can be determined through the fracture time of a crack. The specimen is simple in configuration and easy in processing; the test is concise and clear in process and convenient in operation; and researches on crack dynamic growth behavior and tests on the dynamic fracture toughness can be realized.

The invention belongs to the field of structural damage tolerance design of composite materials, and specifically relates to a method for analyzing residual compression strength of composite materials after impact damage; and the method is used for determining the residual compression strength of composite materials after impact damage. The method comprises three steps of 1, selecting a Hanshin failure criterion as a damage failure criterion for low-speed impact of a laminated plate according to failure characteristics of a composite material; 2, importing the Hanshin failure criterion by adopting a large-scale dynamic finite element program DYTRAN, so as to calculate a damage area under the low-speed impact of the laminated plate; and 3, carrying out rigidity attenuation on a damage area after the low-speed impact according to the damage area determined in the second step, and calculating the residual compression strength of the laminated plate after the low-speed impact by adopting a global-local model analysis method. The invention discloses a brand new analysis method, so that the residual strength of typical components of composite materials after impact damage can be effectively predicted, and basis is provided for the structural design, analysis and verification of the composite materials of airplanes.

The invention relates to an energy dissipation and shock absorption optimization design method which is a joint optimization design method developed based on an MATLAB (matrix laboratory) genetic algorithm toolbox and an SAP2000 finite element program platform. The method includes a preprocessing part, a cycle call optimization part and a post-processing part; the cycle call optimization part is genetic algorithm optimization and anti-seismic checking calculation, and the preprocessing part is used for determining related parameters and codes and establishing an original model and initialized populations for the cycle call optimization part; the cycle call optimization part is realized through the genetic algorithm toolbox GATBX and SAP2000API function compilation, optimum individual are selected, and the post-processing part is used for extracting and processing optimization results to generate a corresponding calculation sheet. The characteristics of the genetic algorithm used in multi-objective problem optimization and a finite element program used in structural energy dissipation and shock absorption reinforcement analysis are integrated, and the method has bigger advantages than a conventional energy dissipation and shock absorption design method.

The invention relates to a multi-block slide calculation method for a giant landslide, and belongs to the field of geological disaster engineering. The method is used for landslide disaster stability evaluation and engineering prevention and control. The method comprises the steps of: (1) finding out a landslide region range, deformation features and hydrogeology and engineering geology conditions; (2) based on a finite element program, establishing a finite element numerical model of a side slope; (3) inputting formation mechanical parameters of the side slope; (4) inputting reduction parameters and physical and mechanical characters of landslide massif for performing finite element calculation; (5) determining a side slope failure region; (6) continuously reducing failure strength parameters until slide surfaces are communicated; (7) determining a tension failure zone of each stage of slide surface; (8) setting the tension failure zone as a null unit in the model; and (9) re-calculating a new numerical model and searching out multiple stages of slide surfaces. According to the multi-block slide calculation method for the giant landslide, provided by the invention, the multiple stages of slide surfaces of the giant landslide can be accurately searched out and an important basis is provided for giant multi-stage landslide stability evaluation and engineering prevention and control.

The invention belongs to the technical field of building structure design and particularly relates to a checking calculation method for the stability of a building foundation under the earthquake action. According to the method, at first, a general finite element program ANSYS is selected for setting up a three-dimensional space finite element model of a building, loading, calculating and deriving are carried out according to the mechanics principle to obtain a concentrated mass bar unit model with the identical dynamic characteristics, the obtained concentrated mass bar model is guided into an ACS SASSI program, a foundation model is supplemented and set up in the ACS SASSI program according to the factory site foundation condition, the ground motion time history is input, and a calculation kernel of the ACS SASSI program is utilized for achieving calculation functions such as soil-structure interaction (SSI) and foundation nonlinearity. By the adoption of the method, an analysis model can conform to actuality better, the calculation result is more accurate and reliable, and the method is simple, applicable, easy to operate, and convenient to realize in actual engineering.

The invention discloses a welding nut finite element modeling method based on automobile fatigue simulation. The method includes the following steps that welding nut modeling is conducted, wherein a solid element network of a welding nut and a base metal is established through point adding and projection; modeling of a welding-nut and base-metal connection surface is conducted, wherein contact modeling is conducted on a bonding wire, the welding nut and the base metal; modeling of a bolt is conducted, wherein the bolt matched with the welding nut is subjected to modeling; modeling of assemblyof the welding nut and the bolt is conducted, wherein modeling of connection of the bolt and the welding nut is conducted, and pretightening force of the bolt is set; finite element calculation is conducted, the established model is moved into a fatigue finite element calculation model, and through the simulation of a finite element program, the force bearing mode of the finite element welding nutis obtained. The method not only considers the connection structure of a welding place, but also considers the contact effect of the nut to a plate, the deformation behavior of the welding nut in thefatigue simulation can be accurately simulated, and accurate plate stress distribution is obtained.

The invention discloses a testing method for the dynamic crack arrest toughness of an I-type crack under the impact load, and belongs to the technical field of geotechnical engineering. The method includes the steps of conducting impact test research through a large-diameter separated type Hopkinson pressing rod and a test configuration designed by an inventor, determining the dynamic crack toughness of a sample through an experiment-numerical value method, calculating the load borne by the sample and the moment of crack arrest through strain signals obtained through a test, inputting an obtained time travel curve to a finite element program Ansys, calculating the near-field displacement of a crack tip through a 1 / 4 node unit, obtaining the time travel curve of the I-type dynamic rupture stress strength factor of the sample through a displacement extrapolation method, obtaining the strain strength factor value corresponding to the crack arrest moment as the dynamic crack arrest toughness of the material, and then obtaining the pure I-type crack dynamic crack arrest toughness under the effect of the impact load. The method is simple and convenient to operate and should have more practical application meaning in research of the dynamic crack arrest characteristic of rock under the impact dynamic load.

The invention discloses an unconventional double-medium reservoir volume fracturing numerical simulation and parameter optimization method, which comprises the following steps of: developing a hydraulic fracture viscoelastic-plastic damage model based on explicit time integration by considering a friction effect between fractures on the basis of a double-medium pore elasticity theory; simulating random intersection and bifurcation of the hydraulic fractures meeting natural fractures by adopting a method of embedding zero-thickness fracture units into the inner boundaries of the computational model grids, and establishing a mathematical model of volume fracturing hydraulic fracture expansion of the unconventional double-medium reservoir; compiling a finite element program of unconventionalreservoir volume fracturing complex multi-crack fracturing and competitive expansion, and establishing a hydraulic fracturing finite element model for a casing-cement ring-cluster perforation-reservoir matrix with natural fractures; carrying out research on the crack initiation and expansion law of the complex crack in the near-wellbore zone in the cluster i, carrying out the research on the synchronous expansion law of the multi-cluster hydrodynamic main crack and optimizing the volume fracturing construction process. Simulation factors are comprehensive, and an oil and gas exploitation construction scheme is optimized.

The invention discloses a testing method for dynamic crack arrest toughness of an I-type crack under a middle-low speed impact load, and belongs to the technical field of rock and concrete engineering. A test piece is a V-shaped side open single crack half-bordered test piece designed by an inventor. The test piece has a V-shaped bottom boundary, and generates an oblique-upward compression wave under an impact load, and the horizontal component of the compression wave has a compression effect on an expansion crack, thereby further preventing cracking of the expansion crack. The testing methodadopts a middle-low speed drop hammer impact test device, finite element difference software AUTODYN is used for simulation calculation, and results are basically consistent with test results. An obtained load stress time history curve is inputted into finite element program ABAQUS, and a time-history curve of a dynamic stress intensity factor of the I-type crack is obtained. Finally, the dynamiccrack arrest toughness of a material is determined by the crack arrest time of the crack. The testing method for the dynamic crack arrest toughness of the pure I-type crack under the middle-low speedimpact load has the advantages of simple test piece configuration, easy processing, simple and clear test procedure and convenient operation, and is suitable for the research of crack dynamic expansion behaviors and dynamic fracture toughness tests.

The invention relates to a method for calculating and processing inrushing destruction of an excavation and depressurization coupling effect of a foundation pit with confined water. The method comprises the following steps of: (1) selecting a calculation range for the depressurization seepage field and the excavation stress field model of the foundation pit; (2) defining the attribute of a soil material around the pit within the calculation range; (3) creating a geometric solid model by using a finite element program, and carrying out element discretization; (4) setting boundary conditions for the depressurization seepage field and the excavation stress field model of the foundation pit; (5) carrying out finite element analog calculation on the excavation and depressurization coupling effect of a foundation pit model; and (6) obtaining the inrushing plastic failure distribution diagram for the soil around the foundation pit with confined water according to the calculated result, and judging whether the foundation pit can inrush or not. Compared with the prior art, according to the invention, the action of the coupling of a foundation pit excavation effect and a depressurization seepage effect is considered, and the shear strength, the heterogeneity and the elastoplasticity of the soil around the pit are considered, thus the defects of the traditional foundation pit inrushing judgment method are made up, and the accuracy and the reliability for the judgment of the inrushing destruction of the foundation pit are enhanced.

The invention provides a two-stage partitioned two-time polycondensation parallel computing system development method and a corresponding parallel computing system. The two-stage partitioned two-time polycondensation parallel computing system development method comprises the steps of 1 producing all program source codes required by finite element analysis parallel computing through a finite element program automatic generation system according to a mathematical equation expression for solving problems, 2 removing solution linear equation systems in the program source codes and conducting parallelized modification on remaining modules by means of a solution linear equation, 3 developing parallel solution linear equation systems based on two-stage partitioning and two-time polycondensation strategies and conducting seamless integration on the parallel solution linear equation systems and the modules subjected to parallelized modification, and 4 compiling and debugging the integrated system and generating a finite element analysis parallel computing executable program. By means of the two-stage partitioned two-time polycondensation parallel computing system development method and the corresponding parallel computing system, solutions can be provided for the finite element analysis problems in all disciplines and fields, and the quality and efficiency of finite element analysis parallel computing software development can be improved.

The invention discloses a finite element prediction method, under the action of force-electricity-heat multi-field coupling, of electrode-containing ferroelectric single crystal based on phase-field method analysis. The method comprises the following steps of establishing a new phase-field model with the multi-field coupling based on a Landau theory and an expansion interface model; obtaining the model establishment and the grid division of the electrode-containing ferroelectric single crystal by using ANSYS finite element software, and extracting node information and unit information; importing the obtained node information, the obtained unit information and considered different boundary condition information, and realizing a solution to the force-electricity-heat nonlinear coupling problem of the electrode-containing ferroelectric single crystal by using an FEAP fast solver; importing field variables which are obtained through FEAP finite element solution into a postprocessing program which is written by Fortran, and obtaining visual field variable results by using TECPLOT software. According to the method provided by the invention, not only is the process simple, but also the adaptability and the stability are good; a commercial finite element program can be developed to flexibly adapt to problem changes.

The invention relates to a method for constructing a slab-girder structure of a train body in the analysis of a train structure, which is used for solving the static force problem of a slab- girder combined structure in a train body steel structure. The method is characterized by comprising the following steps: constructing the slab-girder combined structure at a connection part of a girder unit and a slab unit, constructing an eccentric displacement relation of eccentric nodes when the slab unit is eccentrically connected with the girder unit; then formulating finite element programs constructing the slab-girder combined structure based on a penalty element method and a Lagrange multiplicator method respectively, wherein the finite element programs are respectively used for solving the static force problem of the slab- girder combined structure and the static force problem when the slab unit or the girder unit independently exists. The method can solve the contribution problem of theslab unit to the combined structure total rigidity matrix because of eccentric connection, has less total elements and nodes of the slab-girder combined structure, little occupation machine time of computation, low computation cost and results more accordant with practical situation, and is an effective method for constructing the slab-girder structure of the train.

The invention discloses an overground-underground bidirectional ultra-long concrete seamless design method and belongs to the field of constructional engineering. The method is characterized in that home and abroad large finite element programs, such as ANSYS, ABAQUS, SAP2000 and PMSAP, are used for calculation and analysis, and indoor and outdoor temperature data actually measured in practical projects is applied in the calculation so as to actually measure the regional average temperature; the maximum temperature difference of concrete shrinkage is calculated; and the division of a post-pouring strip or an expansion reinforcing strip and the restrained expansion rate of the used crack resistant waterproofing agent and expanding agent are provided in the design according to the conditions of individual buildings. By utilizing the method, the problem of difficulty in overground-underground bidirectional ultra-long concrete seamless design is solved successfully, requirements on entire use functions of the buildings are satisfied, the usable area is increased, and simultaneously, the construction speed is accelerated; problems in aspects of durability, fire resistance, water tightness, maintenance and the like and building defects caused by arrangement of a deformation joint are solved; and the arrangement of equipment pipelines can be facilitated.

The invention relates to a method for calculating a continuous wall with supporting legs under the action of horizontal loads. According to the specific constructing mode of the underground continuous wall with the supporting legs, firstly, a finite element program of a common elastic foundation beam method is modified and improved, a finite element program of a variable stiffness elastic foundation beam method is made, an engineering calculation method of internal force and horizontal displacement under the action of the horizontal loads is formed in combination of a depth judgment method on the interface of the supporting legs, and therefore when the continuous wall is used as a basement exterior wall and a foundation pit support structure, the lateral pressure of a soil body can be borne. The method can make engineer design personnel fast and effectively design and calculate the underground continuous wall with the supporting legs.

The invention discloses a cable force measuring method for a suspender cable of a suspension bridge. The cable force measuring method sequentially comprises the following steps of: step I, measuring the shape of a main cable of the suspension bridge; step II, according to physical parameters of the main cable, computing to obtain dead load intensity q0 of the main cable; step II, determining horizontal force H0 of the main cable; and step IV, constructing a finite element program and identifying node load of a suspension point, thereby obtaining cable force of the suspender cable. The cable force measuring method has relatively high practicability and precision, can be used for simultaneously measuring the shape of the main cable and the cable force of the suspender cable, and is higher in efficiency in comparison with a normal frequency method.

The invention relates to a parameter optimization method for a powder laser 3D (Three-Dimensional) printing forming technology. A computer finite element program is used for carrying out analysis andcalculation on temperature, stress and a strain field when metal powder is formed under different technical parameters, stress and strain values in a part printing process can be tracked, deformationand defects after parts are formed under different technical parameters are predicted, and then, an existing technical parameter is quickly optimized according to an analysis result to quickly obtainthe optimal technical parameter of powder forming. Compared with an existing method for testing through an experiment technology, the method has the advantages of short test period, low test cost andwide powder applicable range and is suitable for printing various types of metal powder.

The invention discloses a high-precision numerical value method for solving an interval heat convection diffusion problem. The high-precision numerical value method comprises the following steps: establishing a heat transfer control equation of the heat convection diffusion problem; introducing an interval variable to represent the uncertainty of an input parameter in a heat transfer model, and establishing an interval control equation of the heat convection diffusion problem; utilizing a Legendre polynomial expansion to carry out approximate representation on temperature response in the interval control equation; according to tensor product operation, determining a collocation point set in a high-dimension space; utilizing a finite element program to calculate the temperature response of all collocation points, establishing a linear equation set in regard to an expansion coefficient in a temperature response approximate expression, and adopting a least square method to solve; and on the basis of the smoothness of a polynomial function, determining an extreme point of the temperature response approximate expression, and obtaining the upper bound and the lower bound of interval temperature response. The high-precision numerical value method can systematically solve the heat convection diffusion problem with interval uncertainty parameters and effectively improves the calculation precision of the interval numerical value method.

The invention belongs to the technical field of nuclear building design, and relates to a floor response spectrum simplified calculation method considering an SSSI effect. The simplified calculation method is combined with two kinds of finite element programs of the general finite element program and the SSI analysis program; on one hand, part of the structure units in a building group are assembled into a unit by means of super-elements in the general finite element program; on the other hand, by means of format conversion of a dynamic characteristic matrix, the SSI analysis program can directly read unit matrix information of the general finite element program, and floor response spectrum simplified calculation considering the SSI effect is conducted by means of a computing kernel of theSSI analysis program. According to the floor response spectrum simplified calculation method considering the SSSI effect, the SSSI effect can be considered in a real, efficient and rapid mode, calculation results conform to reality on the basis of ensuring consistency and accuracy of a model, and the computing cost is low.

The invention provides a generator stator end part winding structure optimization method based on particle swarm algorithm and support vector machine, which comprises the following steps of: (1) carrying out finite element analysis, recording design variables and structural responses of each analysis as input variables and output variables respectively to obtain a sample point required for establishing an approximate model of a support vector machine; (2) using the obtained sample point, taking SVM parameters as design variables and carrying out initialization; (3) carrying out SVM parameter selection to establish a support vector machine approximate model with optimal parameters to predict structural response; (4) carrying out optimization combined with the approximate model of the support vector machine by using PSO as an optimization algorithm to obtain an optimal solution of the approximate model; (5) substituting the obtained optimal design variable values into the support vectormachine approximate model to call a finite element program for calculation, and comparing the calculation result with the finite element accurate calculation result, and judging the optimal solution.According to the invention, an accurate approximate model is established by using the SVM instead of the finite element analysis, and the optimization efficiency is improved.

A specific force capacity thrust bearing and a method for producing such a design for a particular application are disclosed. The magnetic flux density in the stator material is maximized by varying cross-sectional area normal to the flux path. After a set of initial parameters are chosen, the design can be improved upon by changing the design variables and then verifying the force capacity using a finite element program. By linking the finite element program to a model of the geometry and using some basic algorithms, it is possible to automatically iterate until an optimal design is reached. The resulting design has a much higher force capacity than designs typical of the prior art.

The invention discloses a structural parameter sensitivity analysis method based on a complex variable function method. The structural parameter sensitivity analysis method comprises the following steps that (1), a structural parameter pn to be calculated is selected for small perturbation; (2), a structural finite element calculation program is written, the complex variable function method is implanted into the finite element program, and after the virtual step length is added to the structural parameter, an objective function f(pn+ih) of a structure is solved in the complex space; (3), basedon the objective function f(pn+ih), a first order partial derivative which is the sensitivity of the objective function is solved; (4), by parity of reasoning, when n is smaller than m, n is equal tothe sum of n and 1, and a sensitivity value of the (n+1) structural parameter is obtained; when n is equal to m, the m structural parameter variable is set as a complex variable, and throughthe finite element analysis program, the objective function f(pn+ih) of the structure is solved in the complex space; finally, the sensitivity of the m parameter is solved by the first order partial derivative of the objective function in step (3). The method can utilize the complex variable function method to obtain a high-precision sensitivity matrix of the objective function to structural parameters.

The invention provides a suspension type waterproof curtain deep foundation pit dewatering scheme simulation optimization method, and the method comprises the following steps: S1, establishing a deep foundation pit dewatering hydrogeological conceptual model according to geological and hydrogeological characteristics of a research domain and hydraulic connection between aquifers; S2, according to the hydrogeological conceptual model, establishing a three-dimensional unstable seepage mathematical model and a land subsidence mathematical model of the water-bearing layer soil compression amount, coupling the three-dimensional unstable seepage mathematical model and the land subsidence mathematical model, and developing and solving through a finite element program; S3, carrying out parameter inversion through water pumping test data, and obtaining an optimized three-dimensional coupling numerical model of unstable seepage and land subsidence; S4, carrying out deep foundation pit dewatering and ground subsidence simulation calculation of different combination depths of the suspension type waterproof curtain and the dewatering well filter, and obtaining the optimal deep foundation pit efficient dewatering technical scheme. The coupling model constructed by the method meets certain precision and confidence, and the calculation speed and stability of the model are improved.

The invention provides a solid deformation interface calculation method considering liquid surface tension. A displacement-surface force decoupling numerical method is directly considered, discrete points of a fluid-solid boundary surface in a finite element model are fitted by utilizing a cubic spline curve, the curvature of each point is calculated, overall curvature data of a liquid-solid interface is obtained, and the fluid-solid coupling load problem with surface tension is solved. According to the method, boundary curvature calculation and finite element analysis are alternately carriedout, and the functions of preprocessing, postprocessing, nonlinear constitutive relation and the like of an existing finite element program can be fully utilized.