This paper is devoted to the application of branch mode method in the critical speed ana-lysis of compound rotating systems, in which the distributed inertia including gyroscopic effectand distributed elastic support ...This paper is devoted to the application of branch mode method in the critical speed ana-lysis of compound rotating systems, in which the distributed inertia including gyroscopic effectand distributed elastic support are taken into account. Finally, the method introduced in this paper is used to calculate the critical speeds of anew-type spindle on the spinning frame. The first three critical Speeds are calculated and com-pared with the values obtained from the experimental approach and other theoretical methods.The results show that they are in good agreement with each other.展开更多
An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the rea...An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the real model superposition approximation method to the seismic response calculation of industrial buildings.The analysis using lumped mass-and-shear spring models indicates that for the equipment-to-structure frequency ratiosγf>1.1 orγf<0.9,the non-classical damping effect is limited,and the real mode superposition approximation method provides accurate estimates.For 0.9<γf<1.1,the system may have a pair of closely spaced frequency modes,and the non-zero off-diagonal damping terms have a non-negligible effect on the damping ratios and mode shape vectors of these modes.For 0.9<γf<1.1 and the equipment-to-structure mass ratiosγm<0.07,the real mode superposition approximation method results in large errors,while the approximation method can provide an accurate estimation for 0.9<γf<1.1 andγm>0.07.Furthermore,extensive parametric analyses are conducted,where both steel structures and reinforced concrete structures with equipment with various damping ratios are considered.Finally,the finite element analysis of a five-story industrial building is adopted to validate the proposed range of applicability.展开更多
In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thi...In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thing is to decompose a dynamic system into two sub-systems,a small-scale low-frequency system and a high-frequency system.The former can be analytically and efficiently solved with the mode superposition method,and the latter is dealt with through a time integration method such as the Newmark method.The summation of the responses of these two sub-systems is the responses of the original dynamic system.It is concluded that,with little sacrifice of efficiency,the combination method based on the strategy is more accurate than the combined time integration method,but it has the same accuracy order as that of the combined method.Numerical experiments validate the effectiveness of the proposed strategy.展开更多
The neutral grounding mode of medium-voltage distribution network decides the reliability, overvoltage, relay protection and electrical safety. Therefore, a comprehensive consideration of the reliability, safety and e...The neutral grounding mode of medium-voltage distribution network decides the reliability, overvoltage, relay protection and electrical safety. Therefore, a comprehensive consideration of the reliability, safety and economy is particularly important for the decision of neutral grounding mode. This paper proposes a new decision method of neutral point grounding mode for mediumvoltage distribution network. The objective function is constructed for the decision according the life cycle cost. The reliability of the neutral point grounding mode is taken into account through treating the outage cost as an operating cost. The safety condition of the neutral point grounding mode is preserved as the constraint condition of decision models, so the decision method can generate the most economical and reliable scheme of neutral point grounding mode within a safe limit. The example is used to verify the feasibility and effectiveness of the decision method.展开更多
In this letter, a new moment method using helical segments is presented to model Normal Mode Helical Antenna (NMHA). Using this method, the NMHA can be modeled by a few segments. The current distributions and radiatio...In this letter, a new moment method using helical segments is presented to model Normal Mode Helical Antenna (NMHA). Using this method, the NMHA can be modeled by a few segments. The current distributions and radiation patterns of some NMHAs are calculated.A comparison is made between results obtained using this helical segment algorithm and a linear segment algorithm, and the results of the two algorithms agree fairly well. When calculating the impedance matrix [Z], all the elements of the matrix can be obtained by only calculating a few elements with the application of the symmetric and periodic characteristics of the NMHA.Therefore, the CPU time and the memory storage are significantly reduced, with the accuracy and speed enhanced.展开更多
We compare three different methods to extract coherent modes from Doppler backscattering （DBS）, which are center of gravity （COG） of the complex amplitude spectrum, spectrum of DBS phase derivative （phase derivat...We compare three different methods to extract coherent modes from Doppler backscattering （DBS）, which are center of gravity （COG） of the complex amplitude spectrum, spectrum of DBS phase derivative （phase derivative method）, and phase spectrum, respectively. These three methods are all feasible to extract coherent modes, for example, geodesic acoustic mode oscillation. However, there are stilI differences between dealing with high frequency modes （several hundred kHz） and low frequency modes, （several kHz） hiding in DBS signal. There is a significant amount of power at low frequencies in the phase spectrum, which can be removed by using the phase derivative method and COG. High frequency modes are clearer by using the COO and the phase derivative method than the phase spectrum. The spectrum of DBS amplitude does not show the coherent modes detected by using COG, phase derivative method and phase spectrum. When two Doppler shifted peaks exist, coherent modes and their harmonics appear in the spectrum of DBS amplitude, which are introduced by the DBS phase.展开更多
A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are ...A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.展开更多
This paper presents a hybrid Trefftz （HT） boundary element method （BEM） by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a...This paper presents a hybrid Trefftz （HT） boundary element method （BEM） by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor （SIF） is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by usingthe Galerkin and the collocation techniques is effective for mode III fracture problems.展开更多
The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. Thi...The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.展开更多
<div style="text-align:justify;"> Currently, coupled mode theory (CMT) is widely used for calculating the coupling coefficient of twin-core fibers (TCFs) that are used in a broad range of important app...<div style="text-align:justify;"> Currently, coupled mode theory (CMT) is widely used for calculating the coupling coefficient of twin-core fibers (TCFs) that are used in a broad range of important applications. This approach is highly accurate for scenarios with weak coupling between the cores but shows significant errors in the strong coupling scenarios, necessitating the use of a more accurate method for coupling coefficient calculations. Therefore, in this work, we calculate the coupling coefficients of TCFs using the supermode theory with finite element method (FEM) that has higher accuracy than CMT, particularly for the strong coupling TCF. To investigate the origin of the differences between the results obtained by these two methods, the modal field distributions of the supermodes of TCF are simulated and analyzed in detail. </div>展开更多
Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements ...Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements appear since the presence of shear locking and spurious zero energy modes for thin plate problems. To overcome these shortcomings, we employ the large increment method(LIM) for the analyses of the plate bending problems, and propose a force-based 8-node quadrilateral plate(8NQP) element which is based on MindlinReissner plate theory and has no extra spurious zero energy mode. Several benchmark plate bending problems are presented to illustrate the accuracy and convergence of the plate element by comparing with the analytical solutions and displacement-based plate elements. The results show that the 8-node plate element produces fast convergence and accurate stress distributions in both the moderately thick and thin plate bending problems. The plate element is insensitive to mesh distortion and it can avoid the shear locking for thin plate analysis.展开更多
There are the application scope limits for single differential-mode current injection test method, so in order to carry out injection susceptibility test for two-pieces equipment interconnected with both ends of a cab...There are the application scope limits for single differential-mode current injection test method, so in order to carry out injection susceptibility test for two-pieces equipment interconnected with both ends of a cable simultaneously, a double differential-mode current in- jection test method (DDMCI) is proposed. The method adopted the equivalence source wave theorem and Baum-Liu-Tesche(BLT) equation as its theory foundation. The equivalent corresponding relation between injection voltage and radiation electric field intensity is derived, and the phase relation between the two injection voltage sources is confirmed. The results indicate that the amplitude and phase of the equivalent injection voltage source is closely related to the S parameter of directional coupling device, the transmission line length, and the source vector in BLT equation, but has nothing to do with the reflection coefficient between the two equipment pieces. Therefore, by choosing the right amplitude and phase of the double injection voltage sources, the DDMCI test is equivalent to the radiation test for two interconnected equipment of a system.展开更多
Combining the advantages of the stratified sampling and the importance sampling, a stratified importance sampling method （SISM） is presented to analyze the reliability sensitivity for structure with multiple failure...Combining the advantages of the stratified sampling and the importance sampling, a stratified importance sampling method （SISM） is presented to analyze the reliability sensitivity for structure with multiple failure modes. In the presented method, the variable space is divided into several disjoint subspace by n-dimensional coordinate planes at the mean point of the random vec- tor, and the importance sampling functions in the subspaces are constructed by keeping the sampling center at the mean point and augmenting the standard deviation by a factor of 2. The sample size generated from the importance sampling function in each subspace is determined by the contribution of the subspace to the reliability sensitivity, which can be estimated by iterative simulation in the sampling process. The formulae of the reliability sensitivity estimation, the variance and the coefficient of variation are derived for the presented SISM. Comparing with the Monte Carlo method, the stratified sampling method and the importance sampling method, the presented SISM has wider applicability and higher calculation efficiency, which is demonstrated by numerical examples. Finally, the reliability sensitivity analysis of flap structure is illustrated that the SISM can be applied to engineering structure.展开更多
We propose an improved finite temperature Lanczos method using the stochastic state selection method. In the finite temperature Lanczos method, we generate Lanczos states and calculate the eigenvalues. In addition we ...We propose an improved finite temperature Lanczos method using the stochastic state selection method. In the finite temperature Lanczos method, we generate Lanczos states and calculate the eigenvalues. In addition we have to calculate matrix elements that are the values of an operator between two Lanczos states. In the calculations of the matrix elements we have to keep the set of Lanczos states on the computer memory. Therefore the memory limits the system size in the calculations. Here we propose an application of the stochastic state selection method in order to weaken this limitation. This method is to select some parts of basis states stochastically and to abandon other basis state. Only by the selected basis states we calculate the inner product. After making the statistical average, we can obtain the correct value of the inner product. By the stochastic state selection method we can reduce the number of the basis states for calculations. As a result we can relax the limitation on the computer memory. In order to study the Higgs mode at finite temperature, we calculate the dynamical correlations of the two spin operators in the spin-1/2 Heisenberg antiferromagnet on the square lattice using the improved finite temperature Lanczos method. Our results on the lattices of up to 32 sites show that the Higgs mode exists at low temperature and it disappears gradually when the temperature becomes large. At high temperature we do not find this mode in the dynamical correlations.展开更多
为解决实振型分解反应谱法采用强迫解耦处理调谐质量阻尼器(tuned mass damper,TMD)减震结构具有非比例的阻尼矩阵及计算精度较差的问题,引入复模态理论,利用复振型在状态空间中的正交性和叠加性,并基于复完全平方组合(complex complete...为解决实振型分解反应谱法采用强迫解耦处理调谐质量阻尼器(tuned mass damper,TMD)减震结构具有非比例的阻尼矩阵及计算精度较差的问题,引入复模态理论,利用复振型在状态空间中的正交性和叠加性,并基于复完全平方组合(complex complete quadratic combination,CCQC)法和反应谱,提出一套适用于剪切型TMD减震结构的复振型分解反应谱方法。为匹配抗震规范,基于复振型的复共轭特性以及参数矩阵与状态方程系数矩阵间的关系,推导出实数形式的标准地震作用表达。通过数值分析比较实振型和复振型2种反应谱法的误差发现,实振型误差大于10%,复振型误差在5%以内。研究表明:复振型方法考虑了非比例阻尼和振型速度的影响,比实振型方法具有更高的计算精度和稳定性。匹配规范形式的标准地震作用表达式,可使用设计反应谱,形式简单易用。展开更多
文摘This paper is devoted to the application of branch mode method in the critical speed ana-lysis of compound rotating systems, in which the distributed inertia including gyroscopic effectand distributed elastic support are taken into account. Finally, the method introduced in this paper is used to calculate the critical speeds of anew-type spindle on the spinning frame. The first three critical Speeds are calculated and com-pared with the values obtained from the experimental approach and other theoretical methods.The results show that they are in good agreement with each other.
基金Fund of China National Industrial Building Diagnosis and Reconstruction Engineering Technology Research Center under Grant No.YZA2017Ky03the Beijing Natural Science Foundation under Grant No.JQ18029the National Natural Science Foundation of China under Grant No.52078277。
文摘An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the real model superposition approximation method to the seismic response calculation of industrial buildings.The analysis using lumped mass-and-shear spring models indicates that for the equipment-to-structure frequency ratiosγf>1.1 orγf<0.9,the non-classical damping effect is limited,and the real mode superposition approximation method provides accurate estimates.For 0.9<γf<1.1,the system may have a pair of closely spaced frequency modes,and the non-zero off-diagonal damping terms have a non-negligible effect on the damping ratios and mode shape vectors of these modes.For 0.9<γf<1.1 and the equipment-to-structure mass ratiosγm<0.07,the real mode superposition approximation method results in large errors,while the approximation method can provide an accurate estimation for 0.9<γf<1.1 andγm>0.07.Furthermore,extensive parametric analyses are conducted,where both steel structures and reinforced concrete structures with equipment with various damping ratios are considered.Finally,the finite element analysis of a five-story industrial building is adopted to validate the proposed range of applicability.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11872090 and 12172023).
文摘In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thing is to decompose a dynamic system into two sub-systems,a small-scale low-frequency system and a high-frequency system.The former can be analytically and efficiently solved with the mode superposition method,and the latter is dealt with through a time integration method such as the Newmark method.The summation of the responses of these two sub-systems is the responses of the original dynamic system.It is concluded that,with little sacrifice of efficiency,the combination method based on the strategy is more accurate than the combined time integration method,but it has the same accuracy order as that of the combined method.Numerical experiments validate the effectiveness of the proposed strategy.
文摘The neutral grounding mode of medium-voltage distribution network decides the reliability, overvoltage, relay protection and electrical safety. Therefore, a comprehensive consideration of the reliability, safety and economy is particularly important for the decision of neutral grounding mode. This paper proposes a new decision method of neutral point grounding mode for mediumvoltage distribution network. The objective function is constructed for the decision according the life cycle cost. The reliability of the neutral point grounding mode is taken into account through treating the outage cost as an operating cost. The safety condition of the neutral point grounding mode is preserved as the constraint condition of decision models, so the decision method can generate the most economical and reliable scheme of neutral point grounding mode within a safe limit. The example is used to verify the feasibility and effectiveness of the decision method.
文摘In this letter, a new moment method using helical segments is presented to model Normal Mode Helical Antenna (NMHA). Using this method, the NMHA can be modeled by a few segments. The current distributions and radiation patterns of some NMHAs are calculated.A comparison is made between results obtained using this helical segment algorithm and a linear segment algorithm, and the results of the two algorithms agree fairly well. When calculating the impedance matrix [Z], all the elements of the matrix can be obtained by only calculating a few elements with the application of the symmetric and periodic characteristics of the NMHA.Therefore, the CPU time and the memory storage are significantly reduced, with the accuracy and speed enhanced.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10990210 and 11475173the National Magnetic Confinement Fusion Energy Program of China under Grant Nos 2013GB106002 and 2014GB109002
文摘We compare three different methods to extract coherent modes from Doppler backscattering （DBS）, which are center of gravity （COG） of the complex amplitude spectrum, spectrum of DBS phase derivative （phase derivative method）, and phase spectrum, respectively. These three methods are all feasible to extract coherent modes, for example, geodesic acoustic mode oscillation. However, there are stilI differences between dealing with high frequency modes （several hundred kHz） and low frequency modes, （several kHz） hiding in DBS signal. There is a significant amount of power at low frequencies in the phase spectrum, which can be removed by using the phase derivative method and COG. High frequency modes are clearer by using the COO and the phase derivative method than the phase spectrum. The spectrum of DBS amplitude does not show the coherent modes detected by using COG, phase derivative method and phase spectrum. When two Doppler shifted peaks exist, coherent modes and their harmonics appear in the spectrum of DBS amplitude, which are introduced by the DBS phase.
基金Project supported by the National Natural Science Foundation of China(Nos.59825117 and 50175060).
文摘A meshless method integrated with linear elastic fracture mechanics(LEFM)is presented for 2D mixed-mode crack propagation analysis.The domain is divided automatically into sub-domains based on Voronoi cells,which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path.For each increment of the crack extension,the meshless method is applied to carry out a stress analysis of the cracked structure.The J-integral,which can be decomposed into mode Ⅰ and mode Ⅱ for mixed-mode crack,is used for the evaluation of the stress intensity factors(SIFs).The crack-propagation direction,predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method.The meshless results agree well with the experimental ones,which validates the accuracy and efficiency of the method.
基金the National Natural Science Foundation of China(10472082).
文摘This paper presents a hybrid Trefftz （HT） boundary element method （BEM） by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor （SIF） is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by usingthe Galerkin and the collocation techniques is effective for mode III fracture problems.
基金supported by the National Natural Science Foundation of China(No.11002054)
文摘The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.
文摘<div style="text-align:justify;"> Currently, coupled mode theory (CMT) is widely used for calculating the coupling coefficient of twin-core fibers (TCFs) that are used in a broad range of important applications. This approach is highly accurate for scenarios with weak coupling between the cores but shows significant errors in the strong coupling scenarios, necessitating the use of a more accurate method for coupling coefficient calculations. Therefore, in this work, we calculate the coupling coefficients of TCFs using the supermode theory with finite element method (FEM) that has higher accuracy than CMT, particularly for the strong coupling TCF. To investigate the origin of the differences between the results obtained by these two methods, the modal field distributions of the supermodes of TCF are simulated and analyzed in detail. </div>
基金the National Natural Science Foundation of China(No.10872128)
文摘Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements appear since the presence of shear locking and spurious zero energy modes for thin plate problems. To overcome these shortcomings, we employ the large increment method(LIM) for the analyses of the plate bending problems, and propose a force-based 8-node quadrilateral plate(8NQP) element which is based on MindlinReissner plate theory and has no extra spurious zero energy mode. Several benchmark plate bending problems are presented to illustrate the accuracy and convergence of the plate element by comparing with the analytical solutions and displacement-based plate elements. The results show that the 8-node plate element produces fast convergence and accurate stress distributions in both the moderately thick and thin plate bending problems. The plate element is insensitive to mesh distortion and it can avoid the shear locking for thin plate analysis.
基金Project supported by Arm Pre-research Program （51333040101）, National Defense 973 Program （6131380301 ）, National Natural Science Foundation of China （61040003）.
文摘There are the application scope limits for single differential-mode current injection test method, so in order to carry out injection susceptibility test for two-pieces equipment interconnected with both ends of a cable simultaneously, a double differential-mode current in- jection test method (DDMCI) is proposed. The method adopted the equivalence source wave theorem and Baum-Liu-Tesche(BLT) equation as its theory foundation. The equivalent corresponding relation between injection voltage and radiation electric field intensity is derived, and the phase relation between the two injection voltage sources is confirmed. The results indicate that the amplitude and phase of the equivalent injection voltage source is closely related to the S parameter of directional coupling device, the transmission line length, and the source vector in BLT equation, but has nothing to do with the reflection coefficient between the two equipment pieces. Therefore, by choosing the right amplitude and phase of the double injection voltage sources, the DDMCI test is equivalent to the radiation test for two interconnected equipment of a system.
基金National Natural Science Foundation of China (10572117,10802063,50875213)Aeronautical Science Foundation of China (2007ZA53012)+1 种基金New Century Program For Excellent Talents of Ministry of Education of China (NCET-05-0868)National High-tech Research and Development Program (2007AA04Z401)
文摘Combining the advantages of the stratified sampling and the importance sampling, a stratified importance sampling method （SISM） is presented to analyze the reliability sensitivity for structure with multiple failure modes. In the presented method, the variable space is divided into several disjoint subspace by n-dimensional coordinate planes at the mean point of the random vec- tor, and the importance sampling functions in the subspaces are constructed by keeping the sampling center at the mean point and augmenting the standard deviation by a factor of 2. The sample size generated from the importance sampling function in each subspace is determined by the contribution of the subspace to the reliability sensitivity, which can be estimated by iterative simulation in the sampling process. The formulae of the reliability sensitivity estimation, the variance and the coefficient of variation are derived for the presented SISM. Comparing with the Monte Carlo method, the stratified sampling method and the importance sampling method, the presented SISM has wider applicability and higher calculation efficiency, which is demonstrated by numerical examples. Finally, the reliability sensitivity analysis of flap structure is illustrated that the SISM can be applied to engineering structure.
文摘We propose an improved finite temperature Lanczos method using the stochastic state selection method. In the finite temperature Lanczos method, we generate Lanczos states and calculate the eigenvalues. In addition we have to calculate matrix elements that are the values of an operator between two Lanczos states. In the calculations of the matrix elements we have to keep the set of Lanczos states on the computer memory. Therefore the memory limits the system size in the calculations. Here we propose an application of the stochastic state selection method in order to weaken this limitation. This method is to select some parts of basis states stochastically and to abandon other basis state. Only by the selected basis states we calculate the inner product. After making the statistical average, we can obtain the correct value of the inner product. By the stochastic state selection method we can reduce the number of the basis states for calculations. As a result we can relax the limitation on the computer memory. In order to study the Higgs mode at finite temperature, we calculate the dynamical correlations of the two spin operators in the spin-1/2 Heisenberg antiferromagnet on the square lattice using the improved finite temperature Lanczos method. Our results on the lattices of up to 32 sites show that the Higgs mode exists at low temperature and it disappears gradually when the temperature becomes large. At high temperature we do not find this mode in the dynamical correlations.
