The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear li...The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.展开更多
This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffract...This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.展开更多
This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the d...This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green's fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are deriv...This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are derived based on Biot's theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by: (1) checking the satisfaction extent of the boundary conditions, (2) comparing the degenerated solutions of a single-phased case with well- known solutions, and (3) examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional a...This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.展开更多
Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant s...Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant sources, including force dipole, electric dipole, moment, force dilatation and rotation. Such force and charge sources may model defects like vacancies, foreign particles and dislocations. The locations and orientations of the stress and charge sources with respect to the crack are arbitrary.展开更多
The dynamic soil-tunnel interaction is studied by indirect boundary element method (IBEM), using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subj...The dynamic soil-tunnel interaction is studied by indirect boundary element method (IBEM), using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.展开更多
An exact and complete solution of the problem of a half-planecrack in an infinite transversely isotropic piezoelectric body ispresented. The upper and lower crack faces are assumed to be loadedantisym- metrically by a...An exact and complete solution of the problem of a half-planecrack in an infinite transversely isotropic piezoelectric body ispresented. The upper and lower crack faces are assumed to be loadedantisym- metrically by a couple of tangential point forces inopposite directions. The solution is derived through a lim- itingprocedure from that of a penny-shaped crack. The expressions for theelectroelastic field are given in terms of elementary functions.Finally, the numerical results of the second and third mode stressintensity factors k_2 and k_3 of piezoelectric materials and elasticmaterials are compared in figures.展开更多
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line sourc...Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.展开更多
The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform...The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley's axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.展开更多
Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its face...Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results.展开更多
In this paper,we connsider the plane crack problem of the compound orthotropic material for loads applied symmetrically with respect to the crack plane by means of method of complex variable.
The dynamic stress intensity factors in a half-plane weakened by several finite moving cracks are investigated by employing the Fourier complex transformation. Stress analysis is performed in a half-plane containing a...The dynamic stress intensity factors in a half-plane weakened by several finite moving cracks are investigated by employing the Fourier complex transformation. Stress analysis is performed in a half-plane containing a single dislocation and without dislocation. An exact solution in a closed form to the stress fields and displacement is ob- tained. The Galilean transformation is used to transform between coordinates connected to the cracks. The stress components are of the Cauchy singular kind at the location of dislocation and the point of application of the the influence of crack length and crack running force. Numerical examples demonstrate velocity on the stress intensity factor.展开更多
This paper derives explicit expressions for the propagation of Gaussian beams carrying two vortices of equal charges m = ±1diffracted at a half-plane screen, which enables the study of the dynamic evolution of vo...This paper derives explicit expressions for the propagation of Gaussian beams carrying two vortices of equal charges m = ±1diffracted at a half-plane screen, which enables the study of the dynamic evolution of vortices in the diffraction field. It shows that there may be no vortices, a pair or several pairs of vortices of opposite charges m -=±, -1 in the diffraction field. Pair creation, annihilation and motion of vortices may appear upon propagation. The off-axis distance additionally affects the evolutionary behaviour. In the process the total topological charge is equal to zero, which is unequal to that of the vortex beam at the source plane. A comparison with the free-space propagation of two vortices of equal charges and a further extension are made.展开更多
We employ fundamental equations of non-homogeneous elasticity and Fourierintegral transformations to obtain the general solutions of the stress function.On thebasis of these points of view and when the forces on the b...We employ fundamental equations of non-homogeneous elasticity and Fourierintegral transformations to obtain the general solutions of the stress function.On thebasis of these points of view and when the forces on the boundary are arbityary for nonhomogeneous half-plane problems with the Young’s modulus E(x)-E_0θxp[βx].accurate solutions are obtained At last with the degeneracy it is obtained that thefamous Boussnesq solution and this method is successful.展开更多
Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based o...Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green's functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green's function can be divided into two parts: primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.展开更多
Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "impro...Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "improved cosine half- range expansion" algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.展开更多
Two kinds of symplectic schemes are proposed for the equations of vortex system in half plane.One is the Reflecting-method which is based on the symplectic schemes for the system in the whole plane.The other is constr...Two kinds of symplectic schemes are proposed for the equations of vortex system in half plane.One is the Reflecting-method which is based on the symplectic schemes for the system in the whole plane.The other is constructed directly according to the feature of the system.Aseries of numerical results are presented to show the effectiveness of our schemes.展开更多
This paper investigates the effect of atomic disorder on the electronic structure, magnetism, and half-metallicity of full-Heusler Co2FeSi alloy by using the full-potential linearized augmented plane wave method withi...This paper investigates the effect of atomic disorder on the electronic structure, magnetism, and half-metallicity of full-Heusler Co2FeSi alloy by using the full-potential linearized augmented plane wave method within the generalized gradient approximation (GGA) and GGA-kU schemes. It considers three types of atomic disorders in Co2FeSi alloy: the Co-Fe, Co-Si, and Fe-Si disorders. Total energy calculations show that of the three types of disorders, the Fe-Si disorder is more likely to occur. It finds that for the Co Si disorder, additional states appear in the minority band-gap at the EF and the half-metallcity is substantially destroyed, regardless of the disorder level. On the other hand, the Co-Fe and Fe-Si disorders have little effect on the half-metallicity at a low disorder level. When increasing the disorder levels, the half-metallcity is destroyed at about 9 % of the Co-Fe disorder level, while that stays at 25 % of the Fe-Si disorder level.展开更多
基金The project supported by the National Natural Science Foundation of China
文摘The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.
基金National Natural Science Foundation of China Under Grant No.50908156 and 50978183Tianjin Natural Science Foundation Under Grant No. 07JCZDJC10100
文摘This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green's fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.
基金Program for New Century Excellent Talents in University Under Grant No. NCET-05-0248the Key Program for Applied Basic Research of Tianjin Municipality Under Grant No. 07JCZDJC10100
文摘This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green's functions of compressive and shear wave sources are derived based on Biot's theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by: (1) checking the satisfaction extent of the boundary conditions, (2) comparing the degenerated solutions of a single-phased case with well- known solutions, and (3) examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.
基金Project supported by the National Natural Science Foundation of China (No. 10172075)the Yu-Ying Foundation of Hunan University.
文摘Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant sources, including force dipole, electric dipole, moment, force dilatation and rotation. Such force and charge sources may model defects like vacancies, foreign particles and dislocations. The locations and orientations of the stress and charge sources with respect to the crack are arbitrary.
基金National Natural Science Foundation of China under Grant 51378384Key Project of Natural Science Foundation of Tianjin Municipality under Grant 12JCZDJC29000
文摘The dynamic soil-tunnel interaction is studied by indirect boundary element method (IBEM), using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.
基金the National Natural Science Foundation of China(No.19872060 and 69982009)the Postdoctoral Foundation of China
文摘An exact and complete solution of the problem of a half-planecrack in an infinite transversely isotropic piezoelectric body ispresented. The upper and lower crack faces are assumed to be loadedantisym- metrically by a couple of tangential point forces inopposite directions. The solution is derived through a lim- itingprocedure from that of a penny-shaped crack. The expressions for theelectroelastic field are given in terms of elementary functions.Finally, the numerical results of the second and third mode stressintensity factors k_2 and k_3 of piezoelectric materials and elasticmaterials are compared in figures.
基金supported by National Natural Science Foundation of China (50978183)
文摘Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.
基金supported by the National Natural Science Foundation of China (50908156,50978183)Tianjin Research Program of Application Foundation and Advanced Technology (12JCQNJC04700)
文摘The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley's axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.
基金The project supported by the Guangdong Provincial Natural Science Foundationthe Science Foundation of Shantou University
文摘Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results.
文摘In this paper,we connsider the plane crack problem of the compound orthotropic material for loads applied symmetrically with respect to the crack plane by means of method of complex variable.
文摘The dynamic stress intensity factors in a half-plane weakened by several finite moving cracks are investigated by employing the Fourier complex transformation. Stress analysis is performed in a half-plane containing a single dislocation and without dislocation. An exact solution in a closed form to the stress fields and displacement is ob- tained. The Galilean transformation is used to transform between coordinates connected to the cracks. The stress components are of the Cauchy singular kind at the location of dislocation and the point of application of the the influence of crack length and crack running force. Numerical examples demonstrate velocity on the stress intensity factor.
基金supported by the National Natural Science Foundation of China (Grant No. 10874125)the Foundation of Education Department of Sichuan Province of China (Grant No. 10ZA063)
文摘This paper derives explicit expressions for the propagation of Gaussian beams carrying two vortices of equal charges m = ±1diffracted at a half-plane screen, which enables the study of the dynamic evolution of vortices in the diffraction field. It shows that there may be no vortices, a pair or several pairs of vortices of opposite charges m -=±, -1 in the diffraction field. Pair creation, annihilation and motion of vortices may appear upon propagation. The off-axis distance additionally affects the evolutionary behaviour. In the process the total topological charge is equal to zero, which is unequal to that of the vortex beam at the source plane. A comparison with the free-space propagation of two vortices of equal charges and a further extension are made.
文摘We employ fundamental equations of non-homogeneous elasticity and Fourierintegral transformations to obtain the general solutions of the stress function.On thebasis of these points of view and when the forces on the boundary are arbityary for nonhomogeneous half-plane problems with the Young’s modulus E(x)-E_0θxp[βx].accurate solutions are obtained At last with the degeneracy it is obtained that thefamous Boussnesq solution and this method is successful.
文摘Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green's functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green's function can be divided into two parts: primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.
文摘Diffraction of a two-dimensional (2D) semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called "improved cosine half- range expansion" algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.
文摘Two kinds of symplectic schemes are proposed for the equations of vortex system in half plane.One is the Reflecting-method which is based on the symplectic schemes for the system in the whole plane.The other is constructed directly according to the feature of the system.Aseries of numerical results are presented to show the effectiveness of our schemes.
基金supported by the National Natural Science Foundation of China (Grant No 10664005)
文摘This paper investigates the effect of atomic disorder on the electronic structure, magnetism, and half-metallicity of full-Heusler Co2FeSi alloy by using the full-potential linearized augmented plane wave method within the generalized gradient approximation (GGA) and GGA-kU schemes. It considers three types of atomic disorders in Co2FeSi alloy: the Co-Fe, Co-Si, and Fe-Si disorders. Total energy calculations show that of the three types of disorders, the Fe-Si disorder is more likely to occur. It finds that for the Co Si disorder, additional states appear in the minority band-gap at the EF and the half-metallcity is substantially destroyed, regardless of the disorder level. On the other hand, the Co-Fe and Fe-Si disorders have little effect on the half-metallicity at a low disorder level. When increasing the disorder levels, the half-metallcity is destroyed at about 9 % of the Co-Fe disorder level, while that stays at 25 % of the Fe-Si disorder level.