Micro/nano-thin films are widely used in the fields of micro/nano-electromechanical system(MEMS/NEMS)and flexible electronics,and their mechanical properties have an important impact on the stability and reliability o...Micro/nano-thin films are widely used in the fields of micro/nano-electromechanical system(MEMS/NEMS)and flexible electronics,and their mechanical properties have an important impact on the stability and reliability of components.However,accurate characterization of the mechanical properties of thin films still faces challenges due to the complexity of film-substrate structure,and the characterization efficiency of traditional techniques is insufficient.In this paper,a high-throughput determination method of the elastic modulus of thin films is proposed based on the strain variance method,the feasibility of which is analyzed by the finite element method(FEM),and the specific tensile configuration with array-distributed thin films is designed and optimized.Based on the strain difference between the film-substrate region and the uncoated region,the elastic modulus of multiple films is obtained simultaneously,and the influences of film width,spacing,thickness,and distribution on the measurement of elastic modulus are elucidated.The results show that the change in film width has a more obvious effect on the elastic modulus determination than film spacing and thickness,i.e.,the larger the film width is,the closer the calculation results are to the theoretical value,and the change in calculation results tends to be stabilized when the film width increases to a certain length.Specifically,the simultaneous measurement of the elastic modulus of eight metal films on a polyimide(PI)substrate with a length of 110 mm and a width of 30 mm can be realized,and the testing throughput can be further increased with the extension of the substrate length.This study provides an efficient and low-cost method for measuring the elastic modulus of thin films,which is expected to accelerate the development of new thin film materials.展开更多
Taking a three-cable flexible photovoltaic(PV)support structure as the research subject,a finite element model was established.Utilizing a full-order flutter analysis method,the flutter critical wind speed and flutter...Taking a three-cable flexible photovoltaic(PV)support structure as the research subject,a finite element model was established.Utilizing a full-order flutter analysis method,the flutter critical wind speed and flutter frequency of the flexible PV support structure at a tilt angle of 0°were calculated.The results showed good agreement with wind tunnel test data.Further analysis examined the pretension effects in the load-bearing and stabilizing cables on the natural frequency and flutter critical wind speed of the flexible PV support structure.The research findings indicate increasing the pretension in the load-bearing cables significantly raises the natural frequencies of the first four modes.Specifically,as the pretension in the load-bearing cables increases from 22 to 102 kN,the flutter critical wind speed rises from 17.1 to 21.6 m/s.By contrast,the pretension in the stabilizing cable has a smaller effect on the natural frequency and flutter critical wind speed of the flexible PV support structure.When the pretension in the stabilizing cable increased from 22 to 102 kN,the flutter critical wind speed increased from 17.1 to 17.7 m/s.For wind-resistant design of flexible PV support structures,it is recommended to prioritize increasing the pretension in the load-bearing cables to enhance the structural flutter performance.展开更多
In this paper,we investigate the gravitational quasi-normal modes(QNMs)of Horndeski hairy black holes.By using the ordinary finite element method,the relations between the QNMs and the parameters of black holes are mi...In this paper,we investigate the gravitational quasi-normal modes(QNMs)of Horndeski hairy black holes.By using the ordinary finite element method,the relations between the QNMs and the parameters of black holes are minutely shown.Both the oscillating frequency and the decay rate increase with the increase of the scalar charge Q.While the oscillating frequency and the decay rate decrease with the black hole mass M increasing.And the quasi-period of the oscillation decreases when the angular momentum l of the perturbation field increases.However,the decay rate has barely changed for different values of 1.Thus,the numerical results indicate that Horndeski black holes are stable.展开更多
An innovative approach was introduced for the development of a AA6063 recrystallization model.This method incorporated a regression-based technique for the determination of material constants and introduced novel equa...An innovative approach was introduced for the development of a AA6063 recrystallization model.This method incorporated a regression-based technique for the determination of material constants and introduced novel equations for assessing the grain size evolution.Calibration and validation of this methodology involved a combination of experimentally acquired microstructural data from the extrusion of three different AA6063 profiles and results from the simulation using the Qform Extrusion UK finite element code.The outcomes proved the agreement between experimental findings and numerical prediction of the microstructural evolution.The trend of the grain size variation based on different process parameters was accurately simulated,both after dynamic and static recrystallization,with an error of less than 25% in almost the whole sampling computations.展开更多
A new three-dimensional laterolog array sonde（3D-LS） is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferentia...A new three-dimensional laterolog array sonde（3D-LS） is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferential detection abilities. Six investigation modes are designed using the 3D finite-element method and different investigation depths are simulated based on the pseudo-geometrical factor of the six modes. The invasion profile is described using multi-array radial logs. From the analysis of the pseudo-geometrical factor, the investigation depth of the 3D-LS is about 1.5 m for conductive invasion, which is close to that of the dual laterolog tool but greater than that of the highly integrated azimuthal laterolog sonde. The vertical and azimuthal resolution is also analyzed with the same method. The 3DLS can detect low-resistivity anomalies of 0.5 m thickness and 15？ around the borehole for infinitely thick formations. This study lays the foundation for more work on 3D laterolog array sonde for evaluating low-resistivity anomalies.展开更多
To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direc...To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.展开更多
In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is es...In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.展开更多
A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite e...A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite element method （FEM） and the bulk deformation of flexible-die was analyzed with element-free Galerkin method （EFGM）. The frictional contact between sheet metal and flexible-die was treated by the penalty function method. The sheet elastic flexible-die bulging process was analyzed with the FEM-EFGM program for coupled deformation between sheet metal and bulk flexible-die, called CDSB-FEM-EFGM for short. Compared with finite element code DEFORM-2D and experiment results, the CDSB-FEM-EFGM program is feasible. This method provides a suitable numerical method to analyze sheet flexible-die forming.展开更多
The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, differen...A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.展开更多
In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement wa...In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.展开更多
In order to study the dynamic response of high-voltage transmission lines under mechanical failure, a finite element model of a domestic 500-kV high-voltage transmission line system is established. The initial equilib...In order to study the dynamic response of high-voltage transmission lines under mechanical failure, a finite element model of a domestic 500-kV high-voltage transmission line system is established. The initial equilibrium condition of the coupling system model is verified by nonlinear static analysis. The transient dynamic analysis method is proposed to analyze the variation law of dynamic response under cable or insulator rupture, and the dynamic response of structural elements next to the broken span is calculated. The results show that upper crossarm cable rupture has no effect on cable tension at adjacent suspension points, but it has a significant influence on tension in the insulator and the tower component of the upper crossarm next to the broken span. The peak tension in the conductor of the upper crossarm at the suspension point exceeds the design value under insulator rupture. Insulator rupture has no effect on the tower component of the upper crossarm, but it has a significant influence on insulator tension of the upper crossarm. Insulator rupture should be taken into account in the design of overhead transmission lines. The research results can provide a theoretical basis for the design of transmission lines.展开更多
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton...Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.展开更多
In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was estab...In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.展开更多
In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the sof...In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.展开更多
Identifying suitable processing window is necessary but difficult for achieving favorable microstructure and performance in extrusion of large thick-walled pipe with difficult-to-deform Inconel 625 alloy. In this work...Identifying suitable processing window is necessary but difficult for achieving favorable microstructure and performance in extrusion of large thick-walled pipe with difficult-to-deform Inconel 625 alloy. In this work, a method was established for identifying the extrusion process window considering temperature control using response surface methodology. Firstly, the response surface models, which correlate temperature rise and peak temperature to key extrusion parameters, have been developed by orthogonal regression based on finite element calculated data. Secondly, the coupled effects of the key extrusion parameters on the temperature rise and peak temperature have been disclosed based on the regression models. Lastly, suitable extrusion processing windows, which are described by contour map of peak temperature in the space of extrusion speed and initial billet temperature, have been established for different extrusion ratios. Using the identified process window, a suitable combination of the key extrusion parameters can be determined conveniently and quickly.展开更多
To get the quantitive value of abnormal biological tissues, an inverse algorithm about the Young＇s modulus based on the boundary extraction and the image registration technologies is proposed. With the known displace...To get the quantitive value of abnormal biological tissues, an inverse algorithm about the Young＇s modulus based on the boundary extraction and the image registration technologies is proposed. With the known displacements of boundary tissues and the force distribution, the Young＇s modulus is calculated by constructing the unit system and the inverse finite element method （IFEM）. Then a tough range of the modulus for the whole tissue is estimated referring the value obtained before. The improved particle swarm optimizer （PSO） method is adopted to calculate the whole Yong＇s modulus distribution. The presented algorithm overcomes some limitations in other Young＇s modulus reconstruction methods and relaxes the displacements and force boundary condition requirements. The repetitious numerical simulation shows that errors in boundary displacement are not very sensitive to the estimation of next process; a final feasible solution is obtained by the improved PSO method which is close to the theoretical values obtained during searching in an extensive range.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
文摘Micro/nano-thin films are widely used in the fields of micro/nano-electromechanical system(MEMS/NEMS)and flexible electronics,and their mechanical properties have an important impact on the stability and reliability of components.However,accurate characterization of the mechanical properties of thin films still faces challenges due to the complexity of film-substrate structure,and the characterization efficiency of traditional techniques is insufficient.In this paper,a high-throughput determination method of the elastic modulus of thin films is proposed based on the strain variance method,the feasibility of which is analyzed by the finite element method(FEM),and the specific tensile configuration with array-distributed thin films is designed and optimized.Based on the strain difference between the film-substrate region and the uncoated region,the elastic modulus of multiple films is obtained simultaneously,and the influences of film width,spacing,thickness,and distribution on the measurement of elastic modulus are elucidated.The results show that the change in film width has a more obvious effect on the elastic modulus determination than film spacing and thickness,i.e.,the larger the film width is,the closer the calculation results are to the theoretical value,and the change in calculation results tends to be stabilized when the film width increases to a certain length.Specifically,the simultaneous measurement of the elastic modulus of eight metal films on a polyimide(PI)substrate with a length of 110 mm and a width of 30 mm can be realized,and the testing throughput can be further increased with the extension of the substrate length.This study provides an efficient and low-cost method for measuring the elastic modulus of thin films,which is expected to accelerate the development of new thin film materials.
基金The National Natural Science Foundation of China(No.52338011,52208481),China Postdoctoral Science Foundation(No.2023M730581).
文摘Taking a three-cable flexible photovoltaic(PV)support structure as the research subject,a finite element model was established.Utilizing a full-order flutter analysis method,the flutter critical wind speed and flutter frequency of the flexible PV support structure at a tilt angle of 0°were calculated.The results showed good agreement with wind tunnel test data.Further analysis examined the pretension effects in the load-bearing and stabilizing cables on the natural frequency and flutter critical wind speed of the flexible PV support structure.The research findings indicate increasing the pretension in the load-bearing cables significantly raises the natural frequencies of the first four modes.Specifically,as the pretension in the load-bearing cables increases from 22 to 102 kN,the flutter critical wind speed rises from 17.1 to 21.6 m/s.By contrast,the pretension in the stabilizing cable has a smaller effect on the natural frequency and flutter critical wind speed of the flexible PV support structure.When the pretension in the stabilizing cable increased from 22 to 102 kN,the flutter critical wind speed increased from 17.1 to 17.7 m/s.For wind-resistant design of flexible PV support structures,it is recommended to prioritize increasing the pretension in the load-bearing cables to enhance the structural flutter performance.
文摘In this paper,we investigate the gravitational quasi-normal modes(QNMs)of Horndeski hairy black holes.By using the ordinary finite element method,the relations between the QNMs and the parameters of black holes are minutely shown.Both the oscillating frequency and the decay rate increase with the increase of the scalar charge Q.While the oscillating frequency and the decay rate decrease with the black hole mass M increasing.And the quasi-period of the oscillation decreases when the angular momentum l of the perturbation field increases.However,the decay rate has barely changed for different values of 1.Thus,the numerical results indicate that Horndeski black holes are stable.
文摘An innovative approach was introduced for the development of a AA6063 recrystallization model.This method incorporated a regression-based technique for the determination of material constants and introduced novel equations for assessing the grain size evolution.Calibration and validation of this methodology involved a combination of experimentally acquired microstructural data from the extrusion of three different AA6063 profiles and results from the simulation using the Qform Extrusion UK finite element code.The outcomes proved the agreement between experimental findings and numerical prediction of the microstructural evolution.The trend of the grain size variation based on different process parameters was accurately simulated,both after dynamic and static recrystallization,with an error of less than 25% in almost the whole sampling computations.
基金sponsored by the National Oil and Gas Major Projects(No.2011ZX05020-009)
文摘A new three-dimensional laterolog array sonde（3D-LS） is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferential detection abilities. Six investigation modes are designed using the 3D finite-element method and different investigation depths are simulated based on the pseudo-geometrical factor of the six modes. The invasion profile is described using multi-array radial logs. From the analysis of the pseudo-geometrical factor, the investigation depth of the 3D-LS is about 1.5 m for conductive invasion, which is close to that of the dual laterolog tool but greater than that of the highly integrated azimuthal laterolog sonde. The vertical and azimuthal resolution is also analyzed with the same method. The 3DLS can detect low-resistivity anomalies of 0.5 m thickness and 15？ around the borehole for infinitely thick formations. This study lays the foundation for more work on 3D laterolog array sonde for evaluating low-resistivity anomalies.
基金financially supported by the National Natural Science Foundation of China(No.41574127 and 41174104)the National Key Technology R&D Program for the 13th five-year plan(No.2016ZX05018006-006)
文摘To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.
基金The Open Research Fund of Key Laboratory of Highway Engineering of Sichuan Province of Southw est Jiaotong University (No.LHTE002201102)
文摘In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.
基金Project(51275130)supported by the National Natural Science Foundation of China
文摘A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite element method （FEM） and the bulk deformation of flexible-die was analyzed with element-free Galerkin method （EFGM）. The frictional contact between sheet metal and flexible-die was treated by the penalty function method. The sheet elastic flexible-die bulging process was analyzed with the FEM-EFGM program for coupled deformation between sheet metal and bulk flexible-die, called CDSB-FEM-EFGM for short. Compared with finite element code DEFORM-2D and experiment results, the CDSB-FEM-EFGM program is feasible. This method provides a suitable numerical method to analyze sheet flexible-die forming.
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
文摘A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.
基金The National Natural Science Foundation of China(No.51378121)the Fok Ying Tung Education Foundation(No.141076)the Scientific Innovation Research of College Graduates in Jiangsu Province(No.KYLX_0164)
文摘In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.
基金The National Natural Science Foundation of China (No.50578038)the Science and Technology Project of the State Grid Corporation of China(No.SGKJ[2007]116)
文摘In order to study the dynamic response of high-voltage transmission lines under mechanical failure, a finite element model of a domestic 500-kV high-voltage transmission line system is established. The initial equilibrium condition of the coupling system model is verified by nonlinear static analysis. The transient dynamic analysis method is proposed to analyze the variation law of dynamic response under cable or insulator rupture, and the dynamic response of structural elements next to the broken span is calculated. The results show that upper crossarm cable rupture has no effect on cable tension at adjacent suspension points, but it has a significant influence on tension in the insulator and the tower component of the upper crossarm next to the broken span. The peak tension in the conductor of the upper crossarm at the suspension point exceeds the design value under insulator rupture. Insulator rupture has no effect on the tower component of the upper crossarm, but it has a significant influence on insulator tension of the upper crossarm. Insulator rupture should be taken into account in the design of overhead transmission lines. The research results can provide a theoretical basis for the design of transmission lines.
基金supported by the 973 Program (Grant No.2007CB209505)the National Natural Science Foundation of China (Grant No.40674061,40704019)
文摘Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.
基金Projects (11202125, 61175038) supported by the National Natural Science Foundation of China
文摘In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.
文摘In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.
基金Project(2009ZX04005-031-11)supported by the Major National Science and Technology Special Project of ChinaProject(318968)supported by the Marie Curie International Research Staff Exchange Scheme(IRSES,Mat Pro Future)within the 7th EC Framework Program(FP7)Project(B08040)supported by the 111 Plan,China
文摘Identifying suitable processing window is necessary but difficult for achieving favorable microstructure and performance in extrusion of large thick-walled pipe with difficult-to-deform Inconel 625 alloy. In this work, a method was established for identifying the extrusion process window considering temperature control using response surface methodology. Firstly, the response surface models, which correlate temperature rise and peak temperature to key extrusion parameters, have been developed by orthogonal regression based on finite element calculated data. Secondly, the coupled effects of the key extrusion parameters on the temperature rise and peak temperature have been disclosed based on the regression models. Lastly, suitable extrusion processing windows, which are described by contour map of peak temperature in the space of extrusion speed and initial billet temperature, have been established for different extrusion ratios. Using the identified process window, a suitable combination of the key extrusion parameters can be determined conveniently and quickly.
文摘To get the quantitive value of abnormal biological tissues, an inverse algorithm about the Young＇s modulus based on the boundary extraction and the image registration technologies is proposed. With the known displacements of boundary tissues and the force distribution, the Young＇s modulus is calculated by constructing the unit system and the inverse finite element method （IFEM）. Then a tough range of the modulus for the whole tissue is estimated referring the value obtained before. The improved particle swarm optimizer （PSO） method is adopted to calculate the whole Yong＇s modulus distribution. The presented algorithm overcomes some limitations in other Young＇s modulus reconstruction methods and relaxes the displacements and force boundary condition requirements. The repetitious numerical simulation shows that errors in boundary displacement are not very sensitive to the estimation of next process; a final feasible solution is obtained by the improved PSO method which is close to the theoretical values obtained during searching in an extensive range.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
