Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction a...Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction and wear during exercise.The vital mechanical function展开更多
In this work,the effect of ultrasonic vibration modes on the mechanical properties and relaxation of residual stress in 6061-T6 aluminum alloy was studied.A new ultrasonic vibration Johnson-Cook model was proposed,and...In this work,the effect of ultrasonic vibration modes on the mechanical properties and relaxation of residual stress in 6061-T6 aluminum alloy was studied.A new ultrasonic vibration Johnson-Cook model was proposed,and the relaxation and distribution of residual stress under ultrasonic vibration were predicted and analyzed using the finite element method(FEM).The mechanical properties of 6061-T6 aluminum alloy under different ultrasonic vibration modes were analyzed through experiments involving notched specimen tensile testing and scanning electron microscopy(SEM)analysis.The findings indicate that ultrasonic vibration treatment during deformation,unloading,and load-holding,as well as treatment with its natural ultrasonic frequency,can effectively release residual stress;however,treatment with its natural frequency has the highest rate of release up to 65.4%.Ultrasonic vibration treatment during deformation better inhibits fracture under the same conditions.The FEM results are in good agreement with the experimental results,and it can be used as a valid tool for predicting residual stress release under ultrasonic vibration.展开更多
At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The who...At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing...Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing the number of design variables. In this study, to achieve lighter weight and higher stiffness, the simplified model of BIP was developed and combined with an optimization procedure;consequently, optimal designs of automotive body B-pillar were produced. B-pillar was divided into four quarters and each quarter was modelled by one simplified beam. In the optimization procedure, depth, width, and thickness of the simplified beams were considered as the design variables.Weight, bending and torsional stiffness were also considered as objective functions. The optimization procedure is composed of six stages: designing the experiments, calculating grey relational grade, calculating signal-to noise ratio,finding an optimum design using Taguchi grey relational analysis, performing sensitivity analysis using analysis of variance(ANOVA) and performing non-dominated sorting and multi-criteria decision making. The results show that the width of lower B-pillar has the highest effect(about 55%) and the obtained optimum design point could reduce the weight of B-pillar by about 40% without reducing the BIP stiffness by more than 1.47%.展开更多
Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the sp...Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the spinal cord,which could be result in life threatening展开更多
Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by ...Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
The Cr-plated coating inside a gun barrel can effectively improve the barrel’s erosion resistance and thus increase the service life.However,due to the cyclic thermal load caused by high-temperature gunpowder,micro-e...The Cr-plated coating inside a gun barrel can effectively improve the barrel’s erosion resistance and thus increase the service life.However,due to the cyclic thermal load caused by high-temperature gunpowder,micro-element damage tends to occur within the Cr coating/steel substrate interface,leading to a gradual deterioration in macro-mechanical properties for the material in the related region.In order to mimic this cyclic thermal load and,thereby,study the thermal erosion behavior of the Cr coating on the barrel’s inner wall,a laser emitter is utilized in the current study.With the help of in-situ tensile test and finite element simulation results,a shear stress distribution law of the Cr coating/steel substrate and a change law of the interface ultimate shear strength are identified.Studies have shown that the Cr coating/steel substrate interface’s ultimate shear strength has a significant weakening effect due to increasing temperature.In this study,the interfacial ultimate shear strength decreases from 2.57 GPa(no erosion)to 1.02 GPa(laser power is 160 W).The data from this experiment is employed to establish a Cr coating/steel substrate interface shear damage model.And this model is used to predict the flaking process of Cr coating by finite element method.The simulation results show that the increase of coating crack spacing and coating thickness will increase the service life of gun barrel.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite ...In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.展开更多
We present an inverse methodology for deriving viscoplasticity constitutive model parameters for use in explicit finite element simulations of dynamic processes using functional experiments, i.e., those which provide ...We present an inverse methodology for deriving viscoplasticity constitutive model parameters for use in explicit finite element simulations of dynamic processes using functional experiments, i.e., those which provide value beyond that of constitutive model development. The developed methodology utilises Bayesian optimisation to minimise the error between experimental measurements and numerical simulations performed in LS-DYNA. We demonstrate the optimisation methodology using high hardness armour steels across three types of experiments that induce a wide range of loading conditions: ballistic penetration, rod-on-anvil, and near-field blast deformation. By utilising such a broad range of conditions for the optimisation, the resulting constitutive model parameters are generalised, i.e., applicable across the range of loading conditions encompassed the by those experiments(e.g., stress states, plastic strain magnitudes, strain rates, etc.). Model constants identified using this methodology are demonstrated to provide a generalisable model with superior predictive accuracy than those derived from conventional mechanical characterisation experiments or optimised from a single experimental condition.展开更多
This article established groundwater flows differential equation mathematical model of San iiang Plain on the hydrology theory foundation, and used the analysis finite element method to liner change the differential e...This article established groundwater flows differential equation mathematical model of San iiang Plain on the hydrology theory foundation, and used the analysis finite element method to liner change the differential equation into the large-scale system of linear equations, it took linear equations as a part of constraint conditions of the optimized model, carried on the groundwater flow status equation and the optimized model the coupling, and carries on the solution with the Lingo software. The results indicated that this local shallow layer groundwater resources were rich and have the big development potential. But recent years water resources disposition was unreasonable and ground water mining quantity was oversized, these caused the region water flux to assume the drop tendency.展开更多
On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary resul...On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary results show that many strong earthquakes occurred in the area of increased stress in the model. Though the results are preliminary, the quasi 3D finite element model is meaningful for strong earthquake prediction.展开更多
To obtain the optimal process parameters of stamping forming, finite element analysis and optimization technique were integrated via transforming multi-objective issue into a single-objective issue. A Pareto-based gen...To obtain the optimal process parameters of stamping forming, finite element analysis and optimization technique were integrated via transforming multi-objective issue into a single-objective issue. A Pareto-based genetic algorithm was applied to optimizing the head stamping forming process. In the proposed optimal model, fracture, wrinkle and thickness varying are a function of several factors, such as fillet radius, draw-bead position, blank size and blank-holding force. Hence, it is necessary to investigate the relationship between the objective functions and the variables in order to make objective functions varying minimized simultaneously. Firstly, the central composite experimental(CCD) with four factors and five levels was applied, and the experimental data based on the central composite experimental were acquired. Then, the response surface model(RSM) was set up and the results of the analysis of variance(ANOVA) show that it is reliable to predict the fracture, wrinkle and thickness varying functions by the response surface model. Finally, a Pareto-based genetic algorithm was used to find out a set of Pareto front, which makes fracture, wrinkle and thickness varying minimized integrally. A head stamping case indicates that the present method has higher precision and practicability compared with the "trial and error" procedure.展开更多
A three-dimensional finite element model was established for the milling of thin-walled parts. The physical model of the milling of the part was established using the AdvantEdge FEM software as the platform. The alumi...A three-dimensional finite element model was established for the milling of thin-walled parts. The physical model of the milling of the part was established using the AdvantEdge FEM software as the platform. The aluminum alloy impeller was designated as the object to be processed and the boundary conditions which met the actual machining were set. Through the solution, the physical quantities such as the three-way cutting force, the tool temperature, and the tool stress were obtained, and the calculation of the elastic deformation of the thin-walled blade of the free-form surface at the contact points between the tool and the workpiece was realized. The elastic deformation law of the thin-walled blade was then predicted. The results show that the maximum deviation between the predicted value and the actual measured machining value of the elastic deformation was 26.055 μm; the minimum deviation was 2.011 μm, with the average deviation being 10.154 μm. This shows that the prediction is in close agreement with the actual result.展开更多
基金National Natural Science Foundation of China,10872147Natural Science Foundation of Tianjin,09JCYBJC1400
文摘Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction and wear during exercise.The vital mechanical function
基金Project(51775480)supported by the National Natural Science Foundation of ChinaProjects(E2018203143,E2022203050)supported by the Natural Science Foundation of Hebei Province,China。
文摘In this work,the effect of ultrasonic vibration modes on the mechanical properties and relaxation of residual stress in 6061-T6 aluminum alloy was studied.A new ultrasonic vibration Johnson-Cook model was proposed,and the relaxation and distribution of residual stress under ultrasonic vibration were predicted and analyzed using the finite element method(FEM).The mechanical properties of 6061-T6 aluminum alloy under different ultrasonic vibration modes were analyzed through experiments involving notched specimen tensile testing and scanning electron microscopy(SEM)analysis.The findings indicate that ultrasonic vibration treatment during deformation,unloading,and load-holding,as well as treatment with its natural ultrasonic frequency,can effectively release residual stress;however,treatment with its natural frequency has the highest rate of release up to 65.4%.Ultrasonic vibration treatment during deformation better inhibits fracture under the same conditions.The FEM results are in good agreement with the experimental results,and it can be used as a valid tool for predicting residual stress release under ultrasonic vibration.
文摘At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.
文摘Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing the number of design variables. In this study, to achieve lighter weight and higher stiffness, the simplified model of BIP was developed and combined with an optimization procedure;consequently, optimal designs of automotive body B-pillar were produced. B-pillar was divided into four quarters and each quarter was modelled by one simplified beam. In the optimization procedure, depth, width, and thickness of the simplified beams were considered as the design variables.Weight, bending and torsional stiffness were also considered as objective functions. The optimization procedure is composed of six stages: designing the experiments, calculating grey relational grade, calculating signal-to noise ratio,finding an optimum design using Taguchi grey relational analysis, performing sensitivity analysis using analysis of variance(ANOVA) and performing non-dominated sorting and multi-criteria decision making. The results show that the width of lower B-pillar has the highest effect(about 55%) and the obtained optimum design point could reduce the weight of B-pillar by about 40% without reducing the BIP stiffness by more than 1.47%.
文摘Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the spinal cord,which could be result in life threatening
基金National Natural Science Foundation of China (40521002 and 40474013).
文摘Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
文摘The Cr-plated coating inside a gun barrel can effectively improve the barrel’s erosion resistance and thus increase the service life.However,due to the cyclic thermal load caused by high-temperature gunpowder,micro-element damage tends to occur within the Cr coating/steel substrate interface,leading to a gradual deterioration in macro-mechanical properties for the material in the related region.In order to mimic this cyclic thermal load and,thereby,study the thermal erosion behavior of the Cr coating on the barrel’s inner wall,a laser emitter is utilized in the current study.With the help of in-situ tensile test and finite element simulation results,a shear stress distribution law of the Cr coating/steel substrate and a change law of the interface ultimate shear strength are identified.Studies have shown that the Cr coating/steel substrate interface’s ultimate shear strength has a significant weakening effect due to increasing temperature.In this study,the interfacial ultimate shear strength decreases from 2.57 GPa(no erosion)to 1.02 GPa(laser power is 160 W).The data from this experiment is employed to establish a Cr coating/steel substrate interface shear damage model.And this model is used to predict the flaking process of Cr coating by finite element method.The simulation results show that the increase of coating crack spacing and coating thickness will increase the service life of gun barrel.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金the National Natural Science Foundation of China (5963140) Doctor Point Fund of National Education Committee Parent Company Fund of Aviation Industry
文摘In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.
文摘We present an inverse methodology for deriving viscoplasticity constitutive model parameters for use in explicit finite element simulations of dynamic processes using functional experiments, i.e., those which provide value beyond that of constitutive model development. The developed methodology utilises Bayesian optimisation to minimise the error between experimental measurements and numerical simulations performed in LS-DYNA. We demonstrate the optimisation methodology using high hardness armour steels across three types of experiments that induce a wide range of loading conditions: ballistic penetration, rod-on-anvil, and near-field blast deformation. By utilising such a broad range of conditions for the optimisation, the resulting constitutive model parameters are generalised, i.e., applicable across the range of loading conditions encompassed the by those experiments(e.g., stress states, plastic strain magnitudes, strain rates, etc.). Model constants identified using this methodology are demonstrated to provide a generalisable model with superior predictive accuracy than those derived from conventional mechanical characterisation experiments or optimised from a single experimental condition.
文摘This article established groundwater flows differential equation mathematical model of San iiang Plain on the hydrology theory foundation, and used the analysis finite element method to liner change the differential equation into the large-scale system of linear equations, it took linear equations as a part of constraint conditions of the optimized model, carried on the groundwater flow status equation and the optimized model the coupling, and carries on the solution with the Lingo software. The results indicated that this local shallow layer groundwater resources were rich and have the big development potential. But recent years water resources disposition was unreasonable and ground water mining quantity was oversized, these caused the region water flux to assume the drop tendency.
文摘On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary results show that many strong earthquakes occurred in the area of increased stress in the model. Though the results are preliminary, the quasi 3D finite element model is meaningful for strong earthquake prediction.
基金Project(2012ZX04010-081) supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China
文摘To obtain the optimal process parameters of stamping forming, finite element analysis and optimization technique were integrated via transforming multi-objective issue into a single-objective issue. A Pareto-based genetic algorithm was applied to optimizing the head stamping forming process. In the proposed optimal model, fracture, wrinkle and thickness varying are a function of several factors, such as fillet radius, draw-bead position, blank size and blank-holding force. Hence, it is necessary to investigate the relationship between the objective functions and the variables in order to make objective functions varying minimized simultaneously. Firstly, the central composite experimental(CCD) with four factors and five levels was applied, and the experimental data based on the central composite experimental were acquired. Then, the response surface model(RSM) was set up and the results of the analysis of variance(ANOVA) show that it is reliable to predict the fracture, wrinkle and thickness varying functions by the response surface model. Finally, a Pareto-based genetic algorithm was used to find out a set of Pareto front, which makes fracture, wrinkle and thickness varying minimized integrally. A head stamping case indicates that the present method has higher precision and practicability compared with the "trial and error" procedure.
基金Project(U1530138)supported by the National Natural Science Foundation of ChinaProject(A1-8903-17-0103)supported by the Natural Science Foundation of Shanghai Municipal Education Commission,China
文摘A three-dimensional finite element model was established for the milling of thin-walled parts. The physical model of the milling of the part was established using the AdvantEdge FEM software as the platform. The aluminum alloy impeller was designated as the object to be processed and the boundary conditions which met the actual machining were set. Through the solution, the physical quantities such as the three-way cutting force, the tool temperature, and the tool stress were obtained, and the calculation of the elastic deformation of the thin-walled blade of the free-form surface at the contact points between the tool and the workpiece was realized. The elastic deformation law of the thin-walled blade was then predicted. The results show that the maximum deviation between the predicted value and the actual measured machining value of the elastic deformation was 26.055 μm; the minimum deviation was 2.011 μm, with the average deviation being 10.154 μm. This shows that the prediction is in close agreement with the actual result.
