Fatigue crack growth is a critical phenomenon in engineering structures,accounting for a significant percentage of structural failures across various industries.Accurate prediction of crack initiation,propagation path...Fatigue crack growth is a critical phenomenon in engineering structures,accounting for a significant percentage of structural failures across various industries.Accurate prediction of crack initiation,propagation paths,and fatigue life is essential for ensuring structural integrity and optimizing maintenance schedules.This paper presents a comprehensive finite element approach for simulating two-dimensional fatigue crack growth under linear elastic conditionswith adaptivemesh generation.The source code for the programwas developed in Fortran 95 and compiled with Visual Fortran.To achieve high-fidelity simulations,the methodology integrates several key features:it employs an automatic,adaptive meshing technique that selectively refines the element density near the crack front and areas of significant stress concentration.Specialized singular elements are used at the crack tip to ensure precise stress field representation.The direction of crack advancement is predicted using the maximum tangential stress criterion,while stress intensity factors are determined through either the displacement extrapolation technique or the J-integral method.The simulation models crack growth as a series of linear increments,with solution stability maintained by a consistent transfer algorithm and a crack relaxation method.The framework’s effectiveness is demonstrated across various geometries and loading scenarios.Through rigorous validation against both experimental data and established numerical benchmarks,the approach is proven to accurately forecast crack trajectories and fatigue life.Furthermore,the detailed description of the program’s architecture offers a foundational blueprint,serving as a valuable guide for researchers aiming to develop their specialized software for fracture mechanics analysis.展开更多
A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, an...A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, and the transition probability matrix is directly calculated by the results of a discrete element method （DEM） simulation. The Markov property of the BFB is discussed by the comparison results calculated from both static and dynamic transition probability matrices. The static matrix is calculated based on the Markov chain while the dynamic matrix is calculated based on the memory property of the particle movement. Results show that the difference in the trends of particle movement between the static and dynamic matrix calculation is very small. Besides, the particle mixing curves of the MCM and DEM have the same trend and similar numerical values, and the details show the time averaged characteristic of the MCM and also expose its shortcoming in describing the instantaneous particle dynamics in the BFB.展开更多
Using the finite element code ANSYS/LS-DYNA, a dynamic finite element modelwith an elastic-linear-kinematic-hardening plastic material is established to analyzeelastic-plastic stresses in the railhead in the impact pr...Using the finite element code ANSYS/LS-DYNA, a dynamic finite element modelwith an elastic-linear-kinematic-hardening plastic material is established to analyzeelastic-plastic stresses in the railhead in the impact process of wheel and rail occurring at thegap of rail joint. The model is based on the discrete elastic support condition of the rails, whichis suitable for the actual situation of wheel/track rolling contact. In the analysis the influencesof axle load, yield stress and tangent modulus of rail material on the stresses and strains areinvestigated in detail. The distribution of stresses and strains in the jointed railhead are given.It is found that the axle load, yield stress and tangent modulus of rail material greatly affect thestresses and strains in the railhead during impacting. The study provides a reliable method anduseful datum for the further research on fatigue and wear of railhead and improving the rail jointmode.展开更多
A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation ...A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.展开更多
Numerical simulation for fluid flow over an attached rigid body with a deformable ring bubble is analyzed based on the velocity potential theory together with the boundary element method （BEM）. The analysis is focus...Numerical simulation for fluid flow over an attached rigid body with a deformable ring bubble is analyzed based on the velocity potential theory together with the boundary element method （BEM）. The analysis is focused on the axisymmetric case. The bubble surface is treated as a well defined air-liquid interface and is tracked by a mixed Eulerian-Lagrangian method. The points of intersection between the bubble and body are treated, specially in the numerical procedure. The auxiliary function method is adopted to calculate the pressure on the body surface and in the flow field. The convergence study is undertaken to assess the developed numerical method and the computation code. Some case studies are undertaken in which the interactions between the bubble/body and the incoming flow field are simulated. The effects of various physical parameters on the interactions are investigated.展开更多
This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforc...This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforced composite actuator under electromechanical loading. The four-variable refined plate theory is a simple and efficient higher-order shear deformation theory, which predicts parabolic variation of transverse shear stresses across the plate thickness and satisfies zero traction conditions on the plate free surfaces. The weak form of governing equations is derived using the principle of minimum potential energy, and a 4-node non-conforming rectangular plate element with 8 degrees of freedom per node is introduced for discretizing the domain. Several benchmark problems are solved by the developed MATLAB code and the obtained results are compared with those from exact and other numerical solutions, showing good agreement.展开更多
Crustal tectonic activities are essentially the consequences of the accumulation and release of in situ stress. Therefore, studying the stress state near active faults is important for understanding crustal dynamics a...Crustal tectonic activities are essentially the consequences of the accumulation and release of in situ stress. Therefore, studying the stress state near active faults is important for understanding crustal dynamics and earthquake occurrences. In this paper, using in situ stress measurement results obtained by hydraulic fracturing in the vicinity of the Longmenshan fault zone before and after the Wenchuan Ms 8.0 earthquake and finite element modeling, the variation of stress state before and after the Wenchuan M. 8.0 earthquake is investigated. The results show that the shear stress, which is proportional to the difference between principal stresses, increases with depth and distance from the active fault in the calm period or after the earthquakes, and tends to approach to the regional stress level outside the zone influenced by the fault. This distribution appears to gradually reverse with time and the change of fault properties such as frictional strength. With an increase in friction coefficient, low stress areas are reduced and areas with increased stress accumulation are more obvious near the fault. In sections of the fault with high frictional strengths, in situ stress clearly increases in the fault. Stress accumulates more rapidly in the fault zone relative to the surrounding areas, eventually leading to a stress field that peaks at the fault zone. Such a reversal in the stress field between the fault zone and surrounding areas in the magnitude of the stress field is a potential indicator for the occurrence of strong earthquakes.展开更多
The symplectic algorithm and the energy conservation algorithm are two important kinds of algorithms to solve Hamiltonian systems. The symplectic Runge- Kutta （RK） method is an important part of the former, and the ...The symplectic algorithm and the energy conservation algorithm are two important kinds of algorithms to solve Hamiltonian systems. The symplectic Runge- Kutta （RK） method is an important part of the former, and the continuous finite element method （CFEM） belongs to the later. We find and prove the equivalence of one kind of the implicit RK method and the CFEM, give the coefficient table of the CFEM to simplify its computation, propose a new standard to measure algorithms for Hamiltonian systems, and define another class of algorithms --the regular method. Finally, numerical experiments are given to verify the theoretical results.展开更多
The distribution of saturated soft clay is greatly wide in China. The current main measures adopted to deal with soft soil foundations may lead to environmental pollution, even some engineering accidents may happen on...The distribution of saturated soft clay is greatly wide in China. The current main measures adopted to deal with soft soil foundations may lead to environmental pollution, even some engineering accidents may happen on soft soil foundations. In order to solve engineering problems of saturated soft soil foundations well, researches of mechanical properties of them are necessary. One of the most important mechanical characteristics of saturated soft clay is its cyclic accumulative deformation under cyclic loadings. For saturated soft clay, the cyclic accumulative deformation is similar to the creep behavior under static loadings. Therefore, the cyclic accumulative deformation is equivalent to the creep, the number of loading cycles is seen as the time, and this study develops a practical method for predicting the cyclic accumulative deformation of saturated soft clay with the creep theory. The method is a pseudostatic elasto-plastic finite element method implemented by ABAQUS software. A fitted equation between cyclic accumulative strain and number of loading cycles and the empirical relationship of parameters of fitted equation were established with aseries of cyclic triaxial compression tests. Then with this empirical relationship of parameters, the method developed by this study was employed to predict the cyclic accumulative deformation under cyclic triaxial tension tests. Predicted results were in good agreement with test results, and the effectiveness of this method was thus validated for different stress states. The method was then applied in analyzing the cyclic accumulative deformation for soft soil foundation of a pile-supported wharf structure.展开更多
In this paper, the necessary theoretical analysis for the approximation boundary element method to solve dynamical response of viscoelastic thin plate presented in [1] is.discussed. The theorem of existence and uniq...In this paper, the necessary theoretical analysis for the approximation boundary element method to solve dynamical response of viscoelastic thin plate presented in [1] is.discussed. The theorem of existence and uniqueness of the approximate solution andthe error estimation are also obtained. Based on these conclusions , the principle forchoosing the mesh size and the number of truncated terms in the fundamental solution are given. It isshown that the theoretical ana analysis in this paper are consistent with thenumerical results in [1].展开更多
It has been well recognized that, due to anisotropic packing structure of granular material, the true stress in a specimen is different from the applied stress. However, very few research efforts have been focused on ...It has been well recognized that, due to anisotropic packing structure of granular material, the true stress in a specimen is different from the applied stress. However, very few research efforts have been focused on quantifying the relationship between the true stress and applied stress. In this paper, we derive an explicit relationship among applied stress tensor, material-fabric tensor, and force-fabric tensor; and we propose a relationship between the true stress tensor and the applied stress tensor. The validity of this derived relationship is examined by using the discrete element simulation results for granular material under biaxial and triaxial loading con- ditions.展开更多
A novel single-mode single-polarization （SMSP） photonic crystal fibre has been proposed and analysed based on the polarization-dependent coupling and absorption effect via a full-vector finite element method with pe...A novel single-mode single-polarization （SMSP） photonic crystal fibre has been proposed and analysed based on the polarization-dependent coupling and absorption effect via a full-vector finite element method with perfectly matched layers. The numerical results predict that very efficient SMSP operation can be achieved with both high bandwidth and high extinction ratio at low loss penalty. Effects of the fibre structural parameters on the SMSP bandwidth and extinction ratio have been explored, which will provide useful guide for the design and fabrication of the fibre. The results obtained will be instructive for the realization of new SMSP fibres with high performance.展开更多
In this paper, we introduce a z-axis quartz gyroscope using a double-H tuning fork, which has a high sensitivity. However, it also causes a large mechanical quadrature error. The laser trimming method is used to suppr...In this paper, we introduce a z-axis quartz gyroscope using a double-H tuning fork, which has a high sensitivity. However, it also causes a large mechanical quadrature error. The laser trimming method is used to suppress this error at quartz level. The trimming law is obtained through the finite element method （FEM）. A femtosecond laser processing system is used to trim the gold balancing masses on the beams, and experimental results are basically consistent with the simulated ones. The mechanical quadrature error is suppressed by 96%, from 26.3° s-1 to 1.1° s-1. Nonlinearity changes from 1.48% to 0.30%, angular random walk （ARW） is reduced from 2.19° h-1/2 to 1.42° h-1/2, and bias instability is improved by a factor of 7.7, from 197.6° h-1 to 25.4° h-1.展开更多
Objective:This study aims to clarify how the stimulation of acupuncture points is achieved by needles with different surface texture during acupuncture;it also seeks to lessen injury at the insertion site and increase...Objective:This study aims to clarify how the stimulation of acupuncture points is achieved by needles with different surface texture during acupuncture;it also seeks to lessen injury at the insertion site and increase the therapeutic efficacy of acupuncture,by simulating the mechanical effects of various needle surface patterns on Zusanli(ST36)without changing the radius of acupuncture needles.Methods:Five acupuncture needle models with different surface patterns,including the smooth needle,the lined needle,the ringed needle,the left-hand threaded needle and the right-hand threaded needle,and a layered model of the Zusanli acupoint were used to investigate how to reduce tissue damage and increase stimulation during acupuncture treatment,Puncturing of the skin as well as liftinginserting and twisting needle manipulations were simulated using these models,and the degree of damage and force of stimulation caused by the acupuncture needles with different surface patterns during acupuncturewerecompared.Results:The smooth needle and the lined needle caused the least tissue damage during insertion,while the left-hand threaded and the right-hand threaded needles caused the most damage.The ringed needle,the left-hand threaded needle and the right-hand threaded needle stimulated the acupoint tissue more during lifting-inserting manipulations,while the lined needle and the smooth needle produced less stimulation.The stimulation of the lined needle on the acupoint tissue was the largest during twisting manipulation,whereas the left-hand threaded needle and the right-hand threaded needle had smaller effects.In lifting-inserting and twisting manipulations,both the left-hand threaded needle and right-hand threaded needle provided more stimulation,but the torsion direction in which they produced better stimulation was theopposite.Conclusion:According to the simulation results,the ringed pattern enhances stimulation best in the lifting-inserting manipulation,whereas the lined pattern enhances stimulation best in the twisting manipulation.Both the right-hand and left-hand thread patterns have certain enhancing effects in these two operations.Taking the geometric properties of the pattern into account,the left-hand thread pattern and the right-hand thread pattern have the geometric characteristics of both the lined pattern and the ringed pattern.To conclude,a pattern perpendicular to the movement direction during the acupuncture manipulation creates more stimulation.These results have significance for future needle design.展开更多
In this paper, a novel photonic crystal fiber （PCF） with high birefringence and nonlinearity is designed. The charac- teristics of birefringence, dispersion and nonlinearity are studied by using the full-vector fini...In this paper, a novel photonic crystal fiber （PCF） with high birefringence and nonlinearity is designed. The charac- teristics of birefringence, dispersion and nonlinearity are studied by using the full-vector finite element method （FVFEM）. The numerical results show that the phase birefringence and nonlinear coefficient of PCF can be up to 4.51× 10-3 and 32.8972 w-l.km-1 at 1.55 μm, respectively. The proposed PCF could be found to have important applications in the polarization-dependent nonlinear optics such as the pulse compress and reshaping in the C waveband.展开更多
To improve the bending load-carrying capacity （ BLCC） of under-matched butt joint under four-point bending load in the elastic stage, the shape design of the reinforcement is studied based on the theoretics of mecha...To improve the bending load-carrying capacity （ BLCC） of under-matched butt joint under four-point bending load in the elastic stage, the shape design of the reinforcement is studied based on the theoretics of mechanics of materials. The concept, criterion, realization condition and design proposal of equal bending load-carrying capacity （EBLCC） are put forward. The theoretical analysis results have been verified by the finite element method. The simulation results are coincident basically with the ones of theoretical analysis. The research results show that the shape design of the reinforcement of EBLCC can improve BLCC of under-matched butt joint and the unilateral-side type reinforcement can replace double-side symmetry展开更多
Using premium casing connections instead of API ones is one of the mosteffective technique to prevent casing failure. The factors contribute to the strength of premiumcasing connections are studied with FEA and full-s...Using premium casing connections instead of API ones is one of the mosteffective technique to prevent casing failure. The factors contribute to the strength of premiumcasing connections are studied with FEA and full-scale test. The criterions are presented thatensure the connection's strength higher than the pipe. At the same time, the method is given todecrease the peak stress of the connection so as to improve its anticorruption property. At last,full-scale tests are done to test the strength of the connections designed with the methoddescribed, the results show that the connection's strength is higher than the pipe. This indicatedthat the method described is effective in designing premium casing connection.展开更多
The adjustment and the regeneration of the flow field unstructuredmesh as well as the fi- nitelement solutions for the Euler equationsin the transient process of the separation o external stores form anaircraft at a p...The adjustment and the regeneration of the flow field unstructuredmesh as well as the fi- nitelement solutions for the Euler equationsin the transient process of the separation o external stores form anaircraft at a predetermined downward translation and rotation areinvestigated. It is shown that the adjustment or regeneration ofgrids needed for moving store can be confined within a small regionaround the store and the finite element numerical solution method ofEuler equations for fixed grids can also be used for moving grids ifthe conservative fluxes F_k of fixed grids are replaced by E_k ofmoving grids. The relation-ships between F_k and E_k are alsoestablished. The numerical results For practical examples of aircraftwith external stores are shown.展开更多
This paper describes an investigation of a high-strength concrete frame reinforced with high-strength rebars that was tested in the structure engineering laboratory at Shenyang Jianzhu University. The frame specimen w...This paper describes an investigation of a high-strength concrete frame reinforced with high-strength rebars that was tested in the structure engineering laboratory at Shenyang Jianzhu University. The frame specimen was pseudo- dynamically loaded to indicate three earthquake ground motions of different hazard levels, after which the test specimen was subjected to a pseudo-static loading. This paper focuses on the design, construction and experiment of the test frame and validation of the simulation models. Research shows that a high-strength concrete frame reinforced with high-strength rebars is more efficient and economical than a traditional reinforced concrete frame structure. In addition to the economies achieved by effective use of materials, research shows that the frame can provide enough strength to exceed conventional reinforced concrete frames and provide acceptable ductility. The test study provides evidence to validate the performance of a high- strength concrete frame designed according to current seismic code provisions. Based on previous test research, a nonlinear FEM analysis is completcd by using OpenSees software, The dynamic responses of the frame structure are numerically analyzed, The results of the numerical simulation show that the model can calculate the seismic responses of the frame by OpenSees. At the same time, the test provides additional opportunities to validate the performance of the simulation models.展开更多
Based on the Hamilton principle and the moderate deflection beam theory, discretizing the helicopter blade into a number of beam elements with 15 degrees of freedora, and using a quasi-steady aero-model, a nonlinear c...Based on the Hamilton principle and the moderate deflection beam theory, discretizing the helicopter blade into a number of beam elements with 15 degrees of freedora, and using a quasi-steady aero-model, a nonlinear coupled rotor/fuselage equation is established. A periodic solution of blades and fuselage is obtained through aeroelastic coupled trim using the temporal finite element method （TEM）. The Peters dynamic inflow model is used for vehicle stability. A program for computation is developed, which produces the blade responses, hub loads, and rotor pitch controls. The correlation between the analytical results and related literature is good. The converged solution simultaneously satisfies the blade and the vehicle equilibrium equations.展开更多
基金funding of the Deanship of Graduate Studies and Scientific Research,Jazan University,Saudi Arabia,through Project number:JU-20250230-DGSSR-RP-2025.
文摘Fatigue crack growth is a critical phenomenon in engineering structures,accounting for a significant percentage of structural failures across various industries.Accurate prediction of crack initiation,propagation paths,and fatigue life is essential for ensuring structural integrity and optimizing maintenance schedules.This paper presents a comprehensive finite element approach for simulating two-dimensional fatigue crack growth under linear elastic conditionswith adaptivemesh generation.The source code for the programwas developed in Fortran 95 and compiled with Visual Fortran.To achieve high-fidelity simulations,the methodology integrates several key features:it employs an automatic,adaptive meshing technique that selectively refines the element density near the crack front and areas of significant stress concentration.Specialized singular elements are used at the crack tip to ensure precise stress field representation.The direction of crack advancement is predicted using the maximum tangential stress criterion,while stress intensity factors are determined through either the displacement extrapolation technique or the J-integral method.The simulation models crack growth as a series of linear increments,with solution stability maintained by a consistent transfer algorithm and a crack relaxation method.The framework’s effectiveness is demonstrated across various geometries and loading scenarios.Through rigorous validation against both experimental data and established numerical benchmarks,the approach is proven to accurately forecast crack trajectories and fatigue life.Furthermore,the detailed description of the program’s architecture offers a foundational blueprint,serving as a valuable guide for researchers aiming to develop their specialized software for fracture mechanics analysis.
基金The National Science Foundation of China(No.51276036,51306035)the Fundamental Research Funds for the Central Universities(No.KYLX_0114)
文摘A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, and the transition probability matrix is directly calculated by the results of a discrete element method （DEM） simulation. The Markov property of the BFB is discussed by the comparison results calculated from both static and dynamic transition probability matrices. The static matrix is calculated based on the Markov chain while the dynamic matrix is calculated based on the memory property of the particle movement. Results show that the difference in the trends of particle movement between the static and dynamic matrix calculation is very small. Besides, the particle mixing curves of the MCM and DEM have the same trend and similar numerical values, and the details show the time averaged characteristic of the MCM and also expose its shortcoming in describing the instantaneous particle dynamics in the BFB.
基金National Natural Science Foundation of China(No.599355100)Foundation for Excellent PhD Thesis of University of Ministry of Education of China (No.200048)
文摘Using the finite element code ANSYS/LS-DYNA, a dynamic finite element modelwith an elastic-linear-kinematic-hardening plastic material is established to analyzeelastic-plastic stresses in the railhead in the impact process of wheel and rail occurring at thegap of rail joint. The model is based on the discrete elastic support condition of the rails, whichis suitable for the actual situation of wheel/track rolling contact. In the analysis the influencesof axle load, yield stress and tangent modulus of rail material on the stresses and strains areinvestigated in detail. The distribution of stresses and strains in the jointed railhead are given.It is found that the axle load, yield stress and tangent modulus of rail material greatly affect thestresses and strains in the railhead during impacting. The study provides a reliable method anduseful datum for the further research on fatigue and wear of railhead and improving the rail jointmode.
基金Project supported by the National Natural Science Foundation of China(Nos.10971203,11271340,and 11101381)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.
基金supported by the Excellent Young Scientists Fund(No.51222904)the National NaturalScience Foundation of China(Nos.11202057 and 11302056)
文摘Numerical simulation for fluid flow over an attached rigid body with a deformable ring bubble is analyzed based on the velocity potential theory together with the boundary element method （BEM）. The analysis is focused on the axisymmetric case. The bubble surface is treated as a well defined air-liquid interface and is tracked by a mixed Eulerian-Lagrangian method. The points of intersection between the bubble and body are treated, specially in the numerical procedure. The auxiliary function method is adopted to calculate the pressure on the body surface and in the flow field. The convergence study is undertaken to assess the developed numerical method and the computation code. Some case studies are undertaken in which the interactions between the bubble/body and the incoming flow field are simulated. The effects of various physical parameters on the interactions are investigated.
文摘This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforced composite actuator under electromechanical loading. The four-variable refined plate theory is a simple and efficient higher-order shear deformation theory, which predicts parabolic variation of transverse shear stresses across the plate thickness and satisfies zero traction conditions on the plate free surfaces. The weak form of governing equations is derived using the principle of minimum potential energy, and a 4-node non-conforming rectangular plate element with 8 degrees of freedom per node is introduced for discretizing the domain. Several benchmark problems are solved by the developed MATLAB code and the obtained results are compared with those from exact and other numerical solutions, showing good agreement.
基金supported by the research funds of the Institute of Geomechanics,Chinese Academy of Geological Science(Grant No.DZLXJK201404)the Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period(Grant No.2012BAK19B03–3)
文摘Crustal tectonic activities are essentially the consequences of the accumulation and release of in situ stress. Therefore, studying the stress state near active faults is important for understanding crustal dynamics and earthquake occurrences. In this paper, using in situ stress measurement results obtained by hydraulic fracturing in the vicinity of the Longmenshan fault zone before and after the Wenchuan Ms 8.0 earthquake and finite element modeling, the variation of stress state before and after the Wenchuan M. 8.0 earthquake is investigated. The results show that the shear stress, which is proportional to the difference between principal stresses, increases with depth and distance from the active fault in the calm period or after the earthquakes, and tends to approach to the regional stress level outside the zone influenced by the fault. This distribution appears to gradually reverse with time and the change of fault properties such as frictional strength. With an increase in friction coefficient, low stress areas are reduced and areas with increased stress accumulation are more obvious near the fault. In sections of the fault with high frictional strengths, in situ stress clearly increases in the fault. Stress accumulates more rapidly in the fault zone relative to the surrounding areas, eventually leading to a stress field that peaks at the fault zone. Such a reversal in the stress field between the fault zone and surrounding areas in the magnitude of the stress field is a potential indicator for the occurrence of strong earthquakes.
基金Project supported by the National Natural Science Foundation of China (No. 11071067)the Hunan Graduate Student Science and Technology Innovation Project (No. CX2011B184)
文摘The symplectic algorithm and the energy conservation algorithm are two important kinds of algorithms to solve Hamiltonian systems. The symplectic Runge- Kutta （RK） method is an important part of the former, and the continuous finite element method （CFEM） belongs to the later. We find and prove the equivalence of one kind of the implicit RK method and the CFEM, give the coefficient table of the CFEM to simplify its computation, propose a new standard to measure algorithms for Hamiltonian systems, and define another class of algorithms --the regular method. Finally, numerical experiments are given to verify the theoretical results.
基金funded by Science Commission Fund of Chongqing(Grant No.cstc2016jcyj A0123)Open Fund of Key Laboratory of regulation technology for inlandwaterway in transportation industry of Chongqing Jiaotong University (Grant No. NHHD-201506)
文摘The distribution of saturated soft clay is greatly wide in China. The current main measures adopted to deal with soft soil foundations may lead to environmental pollution, even some engineering accidents may happen on soft soil foundations. In order to solve engineering problems of saturated soft soil foundations well, researches of mechanical properties of them are necessary. One of the most important mechanical characteristics of saturated soft clay is its cyclic accumulative deformation under cyclic loadings. For saturated soft clay, the cyclic accumulative deformation is similar to the creep behavior under static loadings. Therefore, the cyclic accumulative deformation is equivalent to the creep, the number of loading cycles is seen as the time, and this study develops a practical method for predicting the cyclic accumulative deformation of saturated soft clay with the creep theory. The method is a pseudostatic elasto-plastic finite element method implemented by ABAQUS software. A fitted equation between cyclic accumulative strain and number of loading cycles and the empirical relationship of parameters of fitted equation were established with aseries of cyclic triaxial compression tests. Then with this empirical relationship of parameters, the method developed by this study was employed to predict the cyclic accumulative deformation under cyclic triaxial tension tests. Predicted results were in good agreement with test results, and the effectiveness of this method was thus validated for different stress states. The method was then applied in analyzing the cyclic accumulative deformation for soft soil foundation of a pile-supported wharf structure.
文摘In this paper, the necessary theoretical analysis for the approximation boundary element method to solve dynamical response of viscoelastic thin plate presented in [1] is.discussed. The theorem of existence and uniqueness of the approximate solution andthe error estimation are also obtained. Based on these conclusions , the principle forchoosing the mesh size and the number of truncated terms in the fundamental solution are given. It isshown that the theoretical ana analysis in this paper are consistent with thenumerical results in [1].
基金the financial support of the National Natural Science Foundation of China(51178044)Program for New Century Excellent Talents in University (NCET-11-0579)
文摘It has been well recognized that, due to anisotropic packing structure of granular material, the true stress in a specimen is different from the applied stress. However, very few research efforts have been focused on quantifying the relationship between the true stress and applied stress. In this paper, we derive an explicit relationship among applied stress tensor, material-fabric tensor, and force-fabric tensor; and we propose a relationship between the true stress tensor and the applied stress tensor. The validity of this derived relationship is examined by using the discrete element simulation results for granular material under biaxial and triaxial loading con- ditions.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 60372100 and 60507007）, and the State Key Development Program for Basic Research of China （Grant No 2003CB314906）.Acknowledgment The first author Liu Xiao-Yi would like to thank professor Yu Chong-Xiu for her constant encouragement and kind support.
文摘A novel single-mode single-polarization （SMSP） photonic crystal fibre has been proposed and analysed based on the polarization-dependent coupling and absorption effect via a full-vector finite element method with perfectly matched layers. The numerical results predict that very efficient SMSP operation can be achieved with both high bandwidth and high extinction ratio at low loss penalty. Effects of the fibre structural parameters on the SMSP bandwidth and extinction ratio have been explored, which will provide useful guide for the design and fabrication of the fibre. The results obtained will be instructive for the realization of new SMSP fibres with high performance.
基金Project supported by the National Natural Science Foundation of China(Grant No.61027007)
文摘In this paper, we introduce a z-axis quartz gyroscope using a double-H tuning fork, which has a high sensitivity. However, it also causes a large mechanical quadrature error. The laser trimming method is used to suppress this error at quartz level. The trimming law is obtained through the finite element method （FEM）. A femtosecond laser processing system is used to trim the gold balancing masses on the beams, and experimental results are basically consistent with the simulated ones. The mechanical quadrature error is suppressed by 96%, from 26.3° s-1 to 1.1° s-1. Nonlinearity changes from 1.48% to 0.30%, angular random walk （ARW） is reduced from 2.19° h-1/2 to 1.42° h-1/2, and bias instability is improved by a factor of 7.7, from 197.6° h-1 to 25.4° h-1.
基金National Natural Science Foundation of China(grant No:12172092,82174488)Shanghai Key Laboratory of Acupuncture Mechanism and Acupoint Function(grant No:21DZ2271800)+1 种基金Military Medical Research Special Project of"Xinhuo Chuancheng"of Naval Medical UniversityMilitary Medical Talent Project of"San Hang"of Naval Medical University(project No:2019-QH-12).
文摘Objective:This study aims to clarify how the stimulation of acupuncture points is achieved by needles with different surface texture during acupuncture;it also seeks to lessen injury at the insertion site and increase the therapeutic efficacy of acupuncture,by simulating the mechanical effects of various needle surface patterns on Zusanli(ST36)without changing the radius of acupuncture needles.Methods:Five acupuncture needle models with different surface patterns,including the smooth needle,the lined needle,the ringed needle,the left-hand threaded needle and the right-hand threaded needle,and a layered model of the Zusanli acupoint were used to investigate how to reduce tissue damage and increase stimulation during acupuncture treatment,Puncturing of the skin as well as liftinginserting and twisting needle manipulations were simulated using these models,and the degree of damage and force of stimulation caused by the acupuncture needles with different surface patterns during acupuncturewerecompared.Results:The smooth needle and the lined needle caused the least tissue damage during insertion,while the left-hand threaded and the right-hand threaded needles caused the most damage.The ringed needle,the left-hand threaded needle and the right-hand threaded needle stimulated the acupoint tissue more during lifting-inserting manipulations,while the lined needle and the smooth needle produced less stimulation.The stimulation of the lined needle on the acupoint tissue was the largest during twisting manipulation,whereas the left-hand threaded needle and the right-hand threaded needle had smaller effects.In lifting-inserting and twisting manipulations,both the left-hand threaded needle and right-hand threaded needle provided more stimulation,but the torsion direction in which they produced better stimulation was theopposite.Conclusion:According to the simulation results,the ringed pattern enhances stimulation best in the lifting-inserting manipulation,whereas the lined pattern enhances stimulation best in the twisting manipulation.Both the right-hand and left-hand thread patterns have certain enhancing effects in these two operations.Taking the geometric properties of the pattern into account,the left-hand thread pattern and the right-hand thread pattern have the geometric characteristics of both the lined pattern and the ringed pattern.To conclude,a pattern perpendicular to the movement direction during the acupuncture manipulation creates more stimulation.These results have significance for future needle design.
基金Project partly supported by the State Major Basic Research Development Program of China(Grant No.2010CB327604)the National Natural Science Foundation of China(Grant No.61377100)
文摘In this paper, a novel photonic crystal fiber （PCF） with high birefringence and nonlinearity is designed. The charac- teristics of birefringence, dispersion and nonlinearity are studied by using the full-vector finite element method （FVFEM）. The numerical results show that the phase birefringence and nonlinear coefficient of PCF can be up to 4.51× 10-3 and 32.8972 w-l.km-1 at 1.55 μm, respectively. The proposed PCF could be found to have important applications in the polarization-dependent nonlinear optics such as the pulse compress and reshaping in the C waveband.
文摘To improve the bending load-carrying capacity （ BLCC） of under-matched butt joint under four-point bending load in the elastic stage, the shape design of the reinforcement is studied based on the theoretics of mechanics of materials. The concept, criterion, realization condition and design proposal of equal bending load-carrying capacity （EBLCC） are put forward. The theoretical analysis results have been verified by the finite element method. The simulation results are coincident basically with the ones of theoretical analysis. The research results show that the shape design of the reinforcement of EBLCC can improve BLCC of under-matched butt joint and the unilateral-side type reinforcement can replace double-side symmetry
文摘Using premium casing connections instead of API ones is one of the mosteffective technique to prevent casing failure. The factors contribute to the strength of premiumcasing connections are studied with FEA and full-scale test. The criterions are presented thatensure the connection's strength higher than the pipe. At the same time, the method is given todecrease the peak stress of the connection so as to improve its anticorruption property. At last,full-scale tests are done to test the strength of the connections designed with the methoddescribed, the results show that the connection's strength is higher than the pipe. This indicatedthat the method described is effective in designing premium casing connection.
文摘The adjustment and the regeneration of the flow field unstructuredmesh as well as the fi- nitelement solutions for the Euler equationsin the transient process of the separation o external stores form anaircraft at a predetermined downward translation and rotation areinvestigated. It is shown that the adjustment or regeneration ofgrids needed for moving store can be confined within a small regionaround the store and the finite element numerical solution method ofEuler equations for fixed grids can also be used for moving grids ifthe conservative fluxes F_k of fixed grids are replaced by E_k ofmoving grids. The relation-ships between F_k and E_k are alsoestablished. The numerical results For practical examples of aircraftwith external stores are shown.
基金Liaoning Provincial Natural Science Foundation of China Under Grant No. 20052005
文摘This paper describes an investigation of a high-strength concrete frame reinforced with high-strength rebars that was tested in the structure engineering laboratory at Shenyang Jianzhu University. The frame specimen was pseudo- dynamically loaded to indicate three earthquake ground motions of different hazard levels, after which the test specimen was subjected to a pseudo-static loading. This paper focuses on the design, construction and experiment of the test frame and validation of the simulation models. Research shows that a high-strength concrete frame reinforced with high-strength rebars is more efficient and economical than a traditional reinforced concrete frame structure. In addition to the economies achieved by effective use of materials, research shows that the frame can provide enough strength to exceed conventional reinforced concrete frames and provide acceptable ductility. The test study provides evidence to validate the performance of a high- strength concrete frame designed according to current seismic code provisions. Based on previous test research, a nonlinear FEM analysis is completcd by using OpenSees software, The dynamic responses of the frame structure are numerically analyzed, The results of the numerical simulation show that the model can calculate the seismic responses of the frame by OpenSees. At the same time, the test provides additional opportunities to validate the performance of the simulation models.
基金Project supported by the National Natural Science Foundation of China (No. 10872089)
文摘Based on the Hamilton principle and the moderate deflection beam theory, discretizing the helicopter blade into a number of beam elements with 15 degrees of freedora, and using a quasi-steady aero-model, a nonlinear coupled rotor/fuselage equation is established. A periodic solution of blades and fuselage is obtained through aeroelastic coupled trim using the temporal finite element method （TEM）. The Peters dynamic inflow model is used for vehicle stability. A program for computation is developed, which produces the blade responses, hub loads, and rotor pitch controls. The correlation between the analytical results and related literature is good. The converged solution simultaneously satisfies the blade and the vehicle equilibrium equations.
