Owing to their global search capabilities and gradient-free operation,metaheuristic algorithms are widely applied to a wide range of optimization problems.However,their computational demands become prohibitive when ta...Owing to their global search capabilities and gradient-free operation,metaheuristic algorithms are widely applied to a wide range of optimization problems.However,their computational demands become prohibitive when tackling high-dimensional optimization challenges.To effectively address these challenges,this study introduces cooperative metaheuristics integrating dynamic dimension reduction(DR).Building upon particle swarm optimization(PSO)and differential evolution(DE),the proposed cooperative methods C-PSO and C-DE are developed.In the proposed methods,the modified principal components analysis(PCA)is utilized to reduce the dimension of design variables,thereby decreasing computational costs.The dynamic DR strategy implements periodic execution of modified PCA after a fixed number of iterations,resulting in the important dimensions being dynamically identified.Compared with the static one,the dynamic DR strategy can achieve precise identification of important dimensions,thereby enabling accelerated convergence toward optimal solutions.Furthermore,the influence of cumulative contribution rate thresholds on optimization problems with different dimensions is investigated.Metaheuristic algorithms(PSO,DE)and cooperative metaheuristics(C-PSO,C-DE)are examined by 15 benchmark functions and two engineering design problems(speed reducer and composite pressure vessel).Comparative results demonstrate that the cooperative methods achieve significantly superior performance compared to standard methods in both solution accuracy and computational efficiency.Compared to standard metaheuristic algorithms,cooperative metaheuristics achieve a reduction in computational cost of at least 40%.The cooperative metaheuristics can be effectively used to tackle both high-dimensional unconstrained and constrained optimization problems.展开更多
This paper extends the original 2D discontinuous deformation analysis(DDA)method proposed by Shi to 3D cases,and presents the formulations of the 3D DDA.The formulations maintain the characteristics of the original 2D...This paper extends the original 2D discontinuous deformation analysis(DDA)method proposed by Shi to 3D cases,and presents the formulations of the 3D DDA.The formulations maintain the characteristics of the original 2D DDA approach.Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact,such as of vertex-to-vertex,vertex- to-edge and edge-to-edge types,can be handled easily based on the C-P method.The matrices of equilibrium equations have been given in detail for programming purposes.The C program codes for the 3D DDA are developed.The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples.The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes.Finally,implications and future extensions are discussed.展开更多
In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a ...In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.展开更多
Shanghai has experienced the greatest land subsidence in China in the past sixty years and produced undesirable environmental impact. However, horizontal ground deformation has not been understood yet. Therefore groun...Shanghai has experienced the greatest land subsidence in China in the past sixty years and produced undesirable environmental impact. However, horizontal ground deformation has not been understood yet. Therefore ground deformation monitoring together with the analysis of its driving forces are critical for geo-hazards early-warning, city planning and sustainable urbanization in Shanghai. In this paper, two-dimensional ground deformation monitoring was performed in Shanghai with SBAS and MSBAS InSAR methods. Twenty-nine Multi-Look Fine 6 （MF6） Radarsat-2 SLC data acquired during 2011-2013 were used to derive vertical ground deformation. Meanwhile, six descending Multi-Look Fine 6 （MF6） and four ascending Multi-Look Fine 2 （MF2） spanning April to August, 2008, were used to derive vertical and horizontal ground deformation during the observation period. The results indicate that vertical and horizontal deformations in 2008 were not homogeneously distributed in different districts ranging from 0-2 cm/year. Vertical deformation rate during 2011-2013 were decreased to less than 1 cm/year in most district of Shanghai area. Activities from groundwater exploitation and rapid urbanization are responsible for most of the ground deformation in Shanghai. Thus, future ground deformation in vertical and horizontal directions should be warranted.展开更多
Based on mesh regeneration and stress interpolation from an old mesh to a new one, a large deformation finite element model is developed for the study of the behaviour of circular plate anchors subjected to uplift loa...Based on mesh regeneration and stress interpolation from an old mesh to a new one, a large deformation finite element model is developed for the study of the behaviour of circular plate anchors subjected to uplift loading. For the deterruination of the distributions of stress components across a clay foundation, the Recovery by Equilibrium in Patches is extended to plastic analyses. ABAQUS, a commercial finite element package, is customized and linked into our program so as to keep automatic and efficient running of large deformation calculation. The quality of stress interpolation is testified by evaluations of Tresca stress and nodal reaction forces. The complete pulling-up processes of plate anchors buried in homogeneous clay arc simulated, and typical pulling force-displacement responses of a deep anchor and a shallow anchor are compared. Different from the results of previous studies, large deformation analysis is of the capability of estimating the breakaway between the anchor bottom and soils. For deep anchors, the variation of mobilized uplift resistance with anchor settlement is composed of three stages, and the initial buried depths of anchors affect the separation embedment slightly. The uplift bearing capacity of deep anchors is usually higher than that of shallow anchors.展开更多
A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream su...A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.展开更多
A deformation monitoring network that covers part of North China area and takes the Beijing region as the center was measured for two times with high precision GPS in 1995 and 1996 respectively. The results from reme...A deformation monitoring network that covers part of North China area and takes the Beijing region as the center was measured for two times with high precision GPS in 1995 and 1996 respectively. The results from remeasurement indicate that present horizontal movement in the monitored area is characterized by relative motion among several main tectonic blocks. Considering the spatial distribution features obtained from geological survey and results on seismic wave and activity in the area, and stratified features of crustal medium in depth, a three dimensional finite element medium model is designed. And under the conditions of taking and not taking the action manner of the background stress field in the studied area into account, the relative motion between tectonic blocks is calculated and modeled. Based on the results from the analysis and calculations the dynamic mechanism for the present horizontal deformation in the area is discussed.展开更多
Quasicrystals （QCs） are sensitive to the piezoelectric （PE） effect. This paper studies static deformation of a multilayered one-dimensional （1D） hexagonal QC plate with the PE effect. The exact closed-form solut...Quasicrystals （QCs） are sensitive to the piezoelectric （PE） effect. This paper studies static deformation of a multilayered one-dimensional （1D） hexagonal QC plate with the PE effect. The exact closed-form solutions of the extended displacement and traction for a homogeneous piezoelectric quasicrystal （PQC）plate are derived from an eigensystem. The general solutions for multilayered PQC plates are then obtained using the propagator matrix method when mechanical and electrical loads are applied on the top surface of the plate. Numerical examples for several sandwich plates made up of PQC, PE, and QC materials are provided to show the effect of stacking sequence on phonon, phason, and electric fields under mechanical and electrical loads, which is useful in designing new composites for engineering structures.展开更多
A three-dimensional （3-D） transient model has been developed to investigate plasma deformation driven by a magnetic field and its influence on arc stability in a circuit breaker. The 3-D distribution of electric cur...A three-dimensional （3-D） transient model has been developed to investigate plasma deformation driven by a magnetic field and its influence on arc stability in a circuit breaker. The 3-D distribution of electric current density is obtained from a current continuity equation along with the generalized Ohm＇s law; while the magnetic field induced by the current flowing through the arc column is calculated by the magnetic vector potential equation. When gas interacts with an arc column, fundamental factors, such as Ampere＇s law, Ohm＇s law, the turbulence model, transport equations of mass, momentum and energy of plasma flow, have to be coupled for aria- lyzing the phenomenon. The coupled interactions between arc and plasma flow are described in the fl＇amework of time-dependent magnetohydrodynamic （MHD） equations in conjunction with a K-~ turbulence model. Simulations have been focused on sausage and kink instabilities in plasma （these phenomena are related tO pinch effects and electromagnetic fields）. The 3-D sjm- ulation reveals the relation between plasma deformation and instability phenomena, which affect arc stability during circuit breaker operation. Plasma deformation is the consequence of coupled interactions between the electromagnetic force and plasma flow described in simulations.展开更多
Based on the finite element numerical algorithm, the coseismic displacements of the Wenchuan Ms8.0 earthquake are calculated with the rupture slip vectors derived by Ji and Hayes as well as Nishimura and Yaji. Except ...Based on the finite element numerical algorithm, the coseismic displacements of the Wenchuan Ms8.0 earthquake are calculated with the rupture slip vectors derived by Ji and Hayes as well as Nishimura and Yaji. Except in a narrow strip around the rupture zone, the coseismic displacements are consistent with those from GPS observation and InSAR interpretation. Numerical results show that rupture slip vectors and elastic properties have profound influences on the surface coseismic deformation, Results from models with different elastic parameters indicate that： （1） in homogeneous elastic medium, the surface displacements are weakly dependent on Poisson＇s ratio and independent of the elastic modulus; （2） in horizontally homogeneous medium with a weak zone at its middle, the thickness of the weak zone plays a significant role on calculating the surface displacements; （3） in horizontally and vertically heterogeneous medium, the surface displacements depend on both Poisson＇s ratio and elastic modulus. Calculations of eoseismic deformation should take account of the spatial variation of the elastic properties. The misfit of the numerical results with that from the GPS observations in the narrow strip around the rupture zone suggests that a much more complicated rupture model of the Wenchuan earthquake needs to be established in future study.展开更多
A new dimensionless number is proposed for dynamic plastic deformation analysis of clamped circular plates under underwater explosion loads by introducing dimensional analysis method to the basic dynamical governing e...A new dimensionless number is proposed for dynamic plastic deformation analysis of clamped circular plates under underwater explosion loads by introducing dimensional analysis method to the basic dynamical governing equations of circular plates.The relation between dimensionless final plastic deformation of circular plates and the new dimensionless number is established based on massive underwater explosion test data.Meanwhile,comparative analysis was discussed with two other published dimensionless parameters which indicated the new dimensionless number proposed in this paper is more effective and extensive to predict the dynamic plastic response of circular plates under underwater explosion condition.展开更多
We investigate the effect of disorder and mechanical deformation on a two-dimensional photonic crystal waveguide. The dispersion characteristics and transmittance of the waveguide are studied using the finite element ...We investigate the effect of disorder and mechanical deformation on a two-dimensional photonic crystal waveguide. The dispersion characteristics and transmittance of the waveguide are studied using the finite element method. Results show that the geometric change of the dielectric material perpendicular to the light propagation direction has a larger influence on the waveguide characteristics than that parallel to the light propagation direction. Mechanical deformation has an obvious influence on the performance of the waveguide. In particular, longitudinal deformed structure exhibits distinct optical characteristics from the ideal one. Studies on this work will provide useful guideline to the fabrication and practical applications based on photonic crystal waveguides.展开更多
Line heating process is a very complex phenomenon as a variety of factors affects the amount of residual deformations. Numerical thermal and mechanical analysis of line heating for prediction of residual deformation i...Line heating process is a very complex phenomenon as a variety of factors affects the amount of residual deformations. Numerical thermal and mechanical analysis of line heating for prediction of residual deformation is time consuming. In the present work dimensional analysis has been presented to obtain a new relationship between input parameters and resulting residual deformations during line heating process. The temperature distribution and residual deformations for 6 mm, 8 mm, 10 mm and 12 mm thick steel plates were numerically estimated and compared with experimental and published results. Extensive data generated through a validated FE model were used to find co-relationship between the input parameters and the resulting residual deformation by multiple regression analysis. The results obtained from the deformation equations developed in this work compared well with those of the FE analysis with a drop in the computation time in the order of 100 (computational time required for FE analysis is around 7 200 second to 9 000 seconds and where the time required for getting the residual deformation by developed equations is only 60 to 90 seconds).展开更多
Using self-researched gas drainage borehole stability dynamic monitoring device, three-dimensional deformation characteristics of borehole under steady vertical load were researched experimentally and systematically. ...Using self-researched gas drainage borehole stability dynamic monitoring device, three-dimensional deformation characteristics of borehole under steady vertical load were researched experimentally and systematically. This research indicated that under the action of steady loading, the mechanical deformation path of the simulated gas drainage borehole is gradually complicated, and the decay of the borehole circumferential strain is an important characterization of the prediction and early warning of borehole instability and collapse. The horizontal position of borehole occurs compressive strain, and the vertical of which occurs tensile strain under the action of vertical stress. At the initial stage of loading, the vertical strain is more sensitive than that in the horizontal direction. After a certain period of time, the horizontal strain is gradually higher than the vertical one, and the intersection of the borehole horizontal diameter and the hole wall is the stress concentration point. With the increase of the depth of hole, the strain shows a gradual decay trend as a whole, and the vertical strain decays more observably, but there is no absolute position correlation between the amount of strain decay and the increase in borehole depth,and the area within 1.5 times the orifice size is the borehole stress concentration zone.展开更多
Diffusion-induced deformation during electrochemical cycling plays an important role in determining structural durability of the electrodes in lithium-ion batteries. In this work, we investigate the coupling between d...Diffusion-induced deformation during electrochemical cycling plays an important role in determining structural durability of the electrodes in lithium-ion batteries. In this work, we investigate the coupling between diffusion and stress in the boundary conditions of a bilayer electrode, and analyze the evolution of the lithium concentration and stress. Numerical simulations are performed under four different combinations of the boundary conditions between diffusion and mechanical deformation. The stress distributes uniformly in the bilayer electrode for all four cases. The concentration of lithium at the interface is discontinuous for the cases with fixed boundary conditions and is continuous for the cases with a surface at stress-free state. For the bilayer electrode fixed at both surfaces, the magnitude of the stress in the bilayer electrode increases with the increase of the diffusion time. This study reveals the importance of incorporating the coupling between diffusion and stress in the boundary condition in the analysis of the structural durability of lithium-ion batteries and in the design of multilayered and/or gradient electrodes.展开更多
The Laji Shan—Jishi Shan tectonic belt(LJTB),located in the southern part of the northeastern Tibetan Plateau(NETP),is a tectonic window to reveal regional tectonic deformation in the NETP.However,its kinematics in t...The Laji Shan—Jishi Shan tectonic belt(LJTB),located in the southern part of the northeastern Tibetan Plateau(NETP),is a tectonic window to reveal regional tectonic deformation in the NETP.However,its kinematics in the Holocene remains controversial.We obtain the latest and dense horizontal velocity field based on data collected from our newly constructed and existing GNSS stations.Combined with fault kinematics from geologic observations,we analyze the crustal deformation characteristics along the LJTB.The results show that:(1)The Laji Shan fault(LJF)is inactive,and the northwest-oriented Jishi Shan fault(JSF)exhibits a significant dextral and thrust slip.(2)The transpression along the arc-shaped LJTB accommodates deformation transformation between the dextral Riyue Shan fault and the sinistral west Qinling fault.(3)With the continuous pushing of the Indian plate,internal strains in the Tibetan Plateau are continuously transferred in the northeast via the LJTB as they are gradually dissipated near the LJTB and translated into significant crustal uplift in these regions.展开更多
Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on ...Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.展开更多
The derivation of nonlinear integrable evolution partial differential equations in higher dimensions has always been the holy grail in the field of integrability.The well-known modified Kd V equation is a prototypical...The derivation of nonlinear integrable evolution partial differential equations in higher dimensions has always been the holy grail in the field of integrability.The well-known modified Kd V equation is a prototypical example of an integrable evolution equation in one spatial dimension.Do there exist integrable analogs of the modified Kd V equation in higher spatial dimensions?In what follows,we present a positive answer to this question.In particular,rewriting the(1+1)-dimensional integrable modified Kd V equation in conservation forms and adding deformation mappings during the process allows one to construct higher-dimensional integrable equations.Further,we illustrate this idea with examples from the modified Kd V hierarchy and also present the Lax pairs of these higher-dimensional integrable evolution equations.展开更多
he convergence between India and Eurasia is partly accommodated by motion of a few large blocks along strike\|slip faults. About 1800km long Altyn Tagh fault strikes in N60～80°E and allows the northeastward disp...he convergence between India and Eurasia is partly accommodated by motion of a few large blocks along strike\|slip faults. About 1800km long Altyn Tagh fault strikes in N60～80°E and allows the northeastward displacement of the Tibet plateau relative to the Tarim. The Altyn Tagh fault zone is a typical transpressional fault zone, characterized by blocks rotation and crustal shortening and vertical extrusion of blocks within the Altyn Tagh strike\|slip system. Differences of three\|dimensional deformation and configuration of the active structures are recognized at different segment of the Altyn Tagh fault zone.1\ Structural configuration of the Altyn Tagh fault zone\;In the Altyn Tagh strike\|slip fault zone, the assemblage pattern of the (active) faults is in the form of parallel plumes, especially in the eastern and the western segments of the Altyn Tagh fault zone. These plumes structures in the eastern segment are assembled by string\|like left lateral strike\|slip fault and broom\|like thrusting faults, and in the western segment by arc\|like left lateral strike\|slip faults along with thrusting faults and normal faults. In the middle segment of the Altyn Tagh fault zone, the structures are characterized by the string\|like left lateral strike\|slip faults in the center and reverse thrusting faults on the two sides.展开更多
Existing experimental results have shown that four types of physical mechanisms, namely, martensite transformation, martensite reorientation, magnetic domain wall motion and magnetization vector rotation, can be activ...Existing experimental results have shown that four types of physical mechanisms, namely, martensite transformation, martensite reorientation, magnetic domain wall motion and magnetization vector rotation, can be activated during the magneto-mechanical deformation of NiMnGa ferromagnetic shape memory alloy (FSMA) single crystals. In this work, based on irreversible thermodynamics, a three-dimensional (3D) single crystal constitutive model is constructed by considering the aforementioned four mechanisms simultaneously. Three types of internal variables, i.e., the volume fraction of each martensite variant, the volume fraction of magnetic domain in each variant and the deviation angle between the magnetization vector, and easy axis are introduced to characterize the magneto-mechanical state of the single crystals. The thermodynamic driving force of each mechanism and the thermodynamic constraints on the constitutive model are obtained from Clausius's dissipative inequality and constructed Gibbs free energy. Then, thermodynamically consistent kinetic equations for the four mechanisms are proposed, respectively. Finally, the ability of the proposed model to describe the magneto-mechanical deformation of NiMnGa FSMA single crystals is verified by comparing the predictions with corresponding experimental results. It is shown that the proposed model can quantitatively capture the main experimental phenomena. Further, the proposed model is used to predict the deformations of the single crystals under the non-proportional mechanical loading conditions.展开更多
基金funded by National Natural Science Foundation of China(Nos.12402142,11832013 and 11572134)Natural Science Foundation of Hubei Province(No.2024AFB235)+1 种基金Hubei Provincial Department of Education Science and Technology Research Project(No.Q20221714)the Opening Foundation of Hubei Key Laboratory of Digital Textile Equipment(Nos.DTL2023019 and DTL2022012).
文摘Owing to their global search capabilities and gradient-free operation,metaheuristic algorithms are widely applied to a wide range of optimization problems.However,their computational demands become prohibitive when tackling high-dimensional optimization challenges.To effectively address these challenges,this study introduces cooperative metaheuristics integrating dynamic dimension reduction(DR).Building upon particle swarm optimization(PSO)and differential evolution(DE),the proposed cooperative methods C-PSO and C-DE are developed.In the proposed methods,the modified principal components analysis(PCA)is utilized to reduce the dimension of design variables,thereby decreasing computational costs.The dynamic DR strategy implements periodic execution of modified PCA after a fixed number of iterations,resulting in the important dimensions being dynamically identified.Compared with the static one,the dynamic DR strategy can achieve precise identification of important dimensions,thereby enabling accelerated convergence toward optimal solutions.Furthermore,the influence of cumulative contribution rate thresholds on optimization problems with different dimensions is investigated.Metaheuristic algorithms(PSO,DE)and cooperative metaheuristics(C-PSO,C-DE)are examined by 15 benchmark functions and two engineering design problems(speed reducer and composite pressure vessel).Comparative results demonstrate that the cooperative methods achieve significantly superior performance compared to standard methods in both solution accuracy and computational efficiency.Compared to standard metaheuristic algorithms,cooperative metaheuristics achieve a reduction in computational cost of at least 40%.The cooperative metaheuristics can be effectively used to tackle both high-dimensional unconstrained and constrained optimization problems.
基金The project supported by the National Natural Science Foundation of China (50139010)
文摘This paper extends the original 2D discontinuous deformation analysis(DDA)method proposed by Shi to 3D cases,and presents the formulations of the 3D DDA.The formulations maintain the characteristics of the original 2D DDA approach.Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact,such as of vertex-to-vertex,vertex- to-edge and edge-to-edge types,can be handled easily based on the C-P method.The matrices of equilibrium equations have been given in detail for programming purposes.The C program codes for the 3D DDA are developed.The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples.The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes.Finally,implications and future extensions are discussed.
基金supported by the Key Project of Chinese National Programs for Fundamental Research and Development(2010CB731502)the National Natural Science Foundation of China(50978745)
文摘In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.
基金supported by the China Science National Foundation (No. 41372353)
文摘Shanghai has experienced the greatest land subsidence in China in the past sixty years and produced undesirable environmental impact. However, horizontal ground deformation has not been understood yet. Therefore ground deformation monitoring together with the analysis of its driving forces are critical for geo-hazards early-warning, city planning and sustainable urbanization in Shanghai. In this paper, two-dimensional ground deformation monitoring was performed in Shanghai with SBAS and MSBAS InSAR methods. Twenty-nine Multi-Look Fine 6 （MF6） Radarsat-2 SLC data acquired during 2011-2013 were used to derive vertical ground deformation. Meanwhile, six descending Multi-Look Fine 6 （MF6） and four ascending Multi-Look Fine 2 （MF2） spanning April to August, 2008, were used to derive vertical and horizontal ground deformation during the observation period. The results indicate that vertical and horizontal deformations in 2008 were not homogeneously distributed in different districts ranging from 0-2 cm/year. Vertical deformation rate during 2011-2013 were decreased to less than 1 cm/year in most district of Shanghai area. Activities from groundwater exploitation and rapid urbanization are responsible for most of the ground deformation in Shanghai. Thus, future ground deformation in vertical and horizontal directions should be warranted.
文摘Based on mesh regeneration and stress interpolation from an old mesh to a new one, a large deformation finite element model is developed for the study of the behaviour of circular plate anchors subjected to uplift loading. For the deterruination of the distributions of stress components across a clay foundation, the Recovery by Equilibrium in Patches is extended to plastic analyses. ABAQUS, a commercial finite element package, is customized and linked into our program so as to keep automatic and efficient running of large deformation calculation. The quality of stress interpolation is testified by evaluations of Tresca stress and nodal reaction forces. The complete pulling-up processes of plate anchors buried in homogeneous clay arc simulated, and typical pulling force-displacement responses of a deep anchor and a shallow anchor are compared. Different from the results of previous studies, large deformation analysis is of the capability of estimating the breakaway between the anchor bottom and soils. For deep anchors, the variation of mobilized uplift resistance with anchor settlement is composed of three stages, and the initial buried depths of anchors affect the separation embedment slightly. The uplift bearing capacity of deep anchors is usually higher than that of shallow anchors.
基金Sponsored by National Natural Science Foundation of China(50175095)Provincial Natural Science Foundation of Hebei of China(502173)
文摘A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.
文摘A deformation monitoring network that covers part of North China area and takes the Beijing region as the center was measured for two times with high precision GPS in 1995 and 1996 respectively. The results from remeasurement indicate that present horizontal movement in the monitored area is characterized by relative motion among several main tectonic blocks. Considering the spatial distribution features obtained from geological survey and results on seismic wave and activity in the area, and stratified features of crustal medium in depth, a three dimensional finite element medium model is designed. And under the conditions of taking and not taking the action manner of the background stress field in the studied area into account, the relative motion between tectonic blocks is calculated and modeled. Based on the results from the analysis and calculations the dynamic mechanism for the present horizontal deformation in the area is discussed.
基金Project supported by the National Natural Science Foundation of China(Nos.11502123 and11262012)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(No.2015JQ01)
文摘Quasicrystals （QCs） are sensitive to the piezoelectric （PE） effect. This paper studies static deformation of a multilayered one-dimensional （1D） hexagonal QC plate with the PE effect. The exact closed-form solutions of the extended displacement and traction for a homogeneous piezoelectric quasicrystal （PQC）plate are derived from an eigensystem. The general solutions for multilayered PQC plates are then obtained using the propagator matrix method when mechanical and electrical loads are applied on the top surface of the plate. Numerical examples for several sandwich plates made up of PQC, PE, and QC materials are provided to show the effect of stacking sequence on phonon, phason, and electric fields under mechanical and electrical loads, which is useful in designing new composites for engineering structures.
文摘A three-dimensional （3-D） transient model has been developed to investigate plasma deformation driven by a magnetic field and its influence on arc stability in a circuit breaker. The 3-D distribution of electric current density is obtained from a current continuity equation along with the generalized Ohm＇s law; while the magnetic field induced by the current flowing through the arc column is calculated by the magnetic vector potential equation. When gas interacts with an arc column, fundamental factors, such as Ampere＇s law, Ohm＇s law, the turbulence model, transport equations of mass, momentum and energy of plasma flow, have to be coupled for aria- lyzing the phenomenon. The coupled interactions between arc and plasma flow are described in the fl＇amework of time-dependent magnetohydrodynamic （MHD） equations in conjunction with a K-~ turbulence model. Simulations have been focused on sausage and kink instabilities in plasma （these phenomena are related tO pinch effects and electromagnetic fields）. The 3-D sjm- ulation reveals the relation between plasma deformation and instability phenomena, which affect arc stability during circuit breaker operation. Plasma deformation is the consequence of coupled interactions between the electromagnetic force and plasma flow described in simulations.
基金supported by the National Natural Science Foundation of China under grant No.40774045supported by the program from Chinese Academy of Sciences under grant No. KZCX2-YW-QN507
文摘Based on the finite element numerical algorithm, the coseismic displacements of the Wenchuan Ms8.0 earthquake are calculated with the rupture slip vectors derived by Ji and Hayes as well as Nishimura and Yaji. Except in a narrow strip around the rupture zone, the coseismic displacements are consistent with those from GPS observation and InSAR interpretation. Numerical results show that rupture slip vectors and elastic properties have profound influences on the surface coseismic deformation, Results from models with different elastic parameters indicate that： （1） in homogeneous elastic medium, the surface displacements are weakly dependent on Poisson＇s ratio and independent of the elastic modulus; （2） in horizontally homogeneous medium with a weak zone at its middle, the thickness of the weak zone plays a significant role on calculating the surface displacements; （3） in horizontally and vertically heterogeneous medium, the surface displacements depend on both Poisson＇s ratio and elastic modulus. Calculations of eoseismic deformation should take account of the spatial variation of the elastic properties. The misfit of the numerical results with that from the GPS observations in the narrow strip around the rupture zone suggests that a much more complicated rupture model of the Wenchuan earthquake needs to be established in future study.
基金supported by the National Natural Science Foundation of China(12402444)。
文摘A new dimensionless number is proposed for dynamic plastic deformation analysis of clamped circular plates under underwater explosion loads by introducing dimensional analysis method to the basic dynamical governing equations of circular plates.The relation between dimensionless final plastic deformation of circular plates and the new dimensionless number is established based on massive underwater explosion test data.Meanwhile,comparative analysis was discussed with two other published dimensionless parameters which indicated the new dimensionless number proposed in this paper is more effective and extensive to predict the dynamic plastic response of circular plates under underwater explosion condition.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 60908028, 60971068, 10979065, and 61275201)the Fundamental Research Funds for the Central Universities of China (Grant Nos. 2011RC0402 and 2012RC0402)the Program for New Century Excellent Talents in University of China (Grant No. NCET-10-0261)
文摘We investigate the effect of disorder and mechanical deformation on a two-dimensional photonic crystal waveguide. The dispersion characteristics and transmittance of the waveguide are studied using the finite element method. Results show that the geometric change of the dielectric material perpendicular to the light propagation direction has a larger influence on the waveguide characteristics than that parallel to the light propagation direction. Mechanical deformation has an obvious influence on the performance of the waveguide. In particular, longitudinal deformed structure exhibits distinct optical characteristics from the ideal one. Studies on this work will provide useful guideline to the fabrication and practical applications based on photonic crystal waveguides.
文摘Line heating process is a very complex phenomenon as a variety of factors affects the amount of residual deformations. Numerical thermal and mechanical analysis of line heating for prediction of residual deformation is time consuming. In the present work dimensional analysis has been presented to obtain a new relationship between input parameters and resulting residual deformations during line heating process. The temperature distribution and residual deformations for 6 mm, 8 mm, 10 mm and 12 mm thick steel plates were numerically estimated and compared with experimental and published results. Extensive data generated through a validated FE model were used to find co-relationship between the input parameters and the resulting residual deformation by multiple regression analysis. The results obtained from the deformation equations developed in this work compared well with those of the FE analysis with a drop in the computation time in the order of 100 (computational time required for FE analysis is around 7 200 second to 9 000 seconds and where the time required for getting the residual deformation by developed equations is only 60 to 90 seconds).
基金financial support of Distinguished scholars of yueqi (NO. 800015Z1179)National Science Fund subsidized project (51474220)Basic scientific research project of the CPC Central Committee (NO. 2009QZ03)
文摘Using self-researched gas drainage borehole stability dynamic monitoring device, three-dimensional deformation characteristics of borehole under steady vertical load were researched experimentally and systematically. This research indicated that under the action of steady loading, the mechanical deformation path of the simulated gas drainage borehole is gradually complicated, and the decay of the borehole circumferential strain is an important characterization of the prediction and early warning of borehole instability and collapse. The horizontal position of borehole occurs compressive strain, and the vertical of which occurs tensile strain under the action of vertical stress. At the initial stage of loading, the vertical strain is more sensitive than that in the horizontal direction. After a certain period of time, the horizontal strain is gradually higher than the vertical one, and the intersection of the borehole horizontal diameter and the hole wall is the stress concentration point. With the increase of the depth of hole, the strain shows a gradual decay trend as a whole, and the vertical strain decays more observably, but there is no absolute position correlation between the amount of strain decay and the increase in borehole depth,and the area within 1.5 times the orifice size is the borehole stress concentration zone.
基金the National Natural Science Foundation of China (Grant 11402054), the Natural Science Foundation of Fujian Provincial (Grant 2018J01663)2016 Open Projects of Key Laboratory for Strength and Vibration of Mechanical Structures (Grant SV2016-KF-18)the National Science Foundation (Grant CMMI-1634540, monitored by Dr. Khershed Cooper).
文摘Diffusion-induced deformation during electrochemical cycling plays an important role in determining structural durability of the electrodes in lithium-ion batteries. In this work, we investigate the coupling between diffusion and stress in the boundary conditions of a bilayer electrode, and analyze the evolution of the lithium concentration and stress. Numerical simulations are performed under four different combinations of the boundary conditions between diffusion and mechanical deformation. The stress distributes uniformly in the bilayer electrode for all four cases. The concentration of lithium at the interface is discontinuous for the cases with fixed boundary conditions and is continuous for the cases with a surface at stress-free state. For the bilayer electrode fixed at both surfaces, the magnitude of the stress in the bilayer electrode increases with the increase of the diffusion time. This study reveals the importance of incorporating the coupling between diffusion and stress in the boundary condition in the analysis of the structural durability of lithium-ion batteries and in the design of multilayered and/or gradient electrodes.
基金supported by the National Science Foundation of China(41874117)the Second Tibetan Plateau Scientific Expedition and Research Program(SETP)(2019QZKK0901)Natural Science Basic Research Program of Shaanxi(Program No.2023-JC-ON-0309)。
文摘The Laji Shan—Jishi Shan tectonic belt(LJTB),located in the southern part of the northeastern Tibetan Plateau(NETP),is a tectonic window to reveal regional tectonic deformation in the NETP.However,its kinematics in the Holocene remains controversial.We obtain the latest and dense horizontal velocity field based on data collected from our newly constructed and existing GNSS stations.Combined with fault kinematics from geologic observations,we analyze the crustal deformation characteristics along the LJTB.The results show that:(1)The Laji Shan fault(LJF)is inactive,and the northwest-oriented Jishi Shan fault(JSF)exhibits a significant dextral and thrust slip.(2)The transpression along the arc-shaped LJTB accommodates deformation transformation between the dextral Riyue Shan fault and the sinistral west Qinling fault.(3)With the continuous pushing of the Indian plate,internal strains in the Tibetan Plateau are continuously transferred in the northeast via the LJTB as they are gradually dissipated near the LJTB and translated into significant crustal uplift in these regions.
基金Project(51008188)supported by National Natural Science Foundation of ChinaProject(KLE-TJGE-B1302)supported by Key Laboratory Fund of Geotechnical and Underground Engineering of Ministry of Education,ChinaProject(SKLGDUEK1205)supported by Open Program of State Key Laboratory for Geomechanics and Deep Underground Engineering,China
文摘Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.
基金sponsored by the National Natural Science Foundations of China(Nos.12235007,11975131,11435005,12275144,11975204)KC Wong Magna Fund in Ningbo UniversityNatural Science Foundation of Zhejiang Province No.LQ20A010009。
文摘The derivation of nonlinear integrable evolution partial differential equations in higher dimensions has always been the holy grail in the field of integrability.The well-known modified Kd V equation is a prototypical example of an integrable evolution equation in one spatial dimension.Do there exist integrable analogs of the modified Kd V equation in higher spatial dimensions?In what follows,we present a positive answer to this question.In particular,rewriting the(1+1)-dimensional integrable modified Kd V equation in conservation forms and adding deformation mappings during the process allows one to construct higher-dimensional integrable equations.Further,we illustrate this idea with examples from the modified Kd V hierarchy and also present the Lax pairs of these higher-dimensional integrable evolution equations.
文摘he convergence between India and Eurasia is partly accommodated by motion of a few large blocks along strike\|slip faults. About 1800km long Altyn Tagh fault strikes in N60～80°E and allows the northeastward displacement of the Tibet plateau relative to the Tarim. The Altyn Tagh fault zone is a typical transpressional fault zone, characterized by blocks rotation and crustal shortening and vertical extrusion of blocks within the Altyn Tagh strike\|slip system. Differences of three\|dimensional deformation and configuration of the active structures are recognized at different segment of the Altyn Tagh fault zone.1\ Structural configuration of the Altyn Tagh fault zone\;In the Altyn Tagh strike\|slip fault zone, the assemblage pattern of the (active) faults is in the form of parallel plumes, especially in the eastern and the western segments of the Altyn Tagh fault zone. These plumes structures in the eastern segment are assembled by string\|like left lateral strike\|slip fault and broom\|like thrusting faults, and in the western segment by arc\|like left lateral strike\|slip faults along with thrusting faults and normal faults. In the middle segment of the Altyn Tagh fault zone, the structures are characterized by the string\|like left lateral strike\|slip faults in the center and reverse thrusting faults on the two sides.
基金the National Natural Science Foundation of China (Grant 11602203)Young Elite Scientist Sponsorship Program by the China Association for Science and Technology (Grant 2016QNRC001)Fundamental Research Funds for the Central Universities (Grant 2682018CX43).
文摘Existing experimental results have shown that four types of physical mechanisms, namely, martensite transformation, martensite reorientation, magnetic domain wall motion and magnetization vector rotation, can be activated during the magneto-mechanical deformation of NiMnGa ferromagnetic shape memory alloy (FSMA) single crystals. In this work, based on irreversible thermodynamics, a three-dimensional (3D) single crystal constitutive model is constructed by considering the aforementioned four mechanisms simultaneously. Three types of internal variables, i.e., the volume fraction of each martensite variant, the volume fraction of magnetic domain in each variant and the deviation angle between the magnetization vector, and easy axis are introduced to characterize the magneto-mechanical state of the single crystals. The thermodynamic driving force of each mechanism and the thermodynamic constraints on the constitutive model are obtained from Clausius's dissipative inequality and constructed Gibbs free energy. Then, thermodynamically consistent kinetic equations for the four mechanisms are proposed, respectively. Finally, the ability of the proposed model to describe the magneto-mechanical deformation of NiMnGa FSMA single crystals is verified by comparing the predictions with corresponding experimental results. It is shown that the proposed model can quantitatively capture the main experimental phenomena. Further, the proposed model is used to predict the deformations of the single crystals under the non-proportional mechanical loading conditions.
