Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explaine...Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explained by analyzing the interface potential, the tangential energy flux, and the normal energy flux. A rigorous criterion for tip-splitting instability is established analytically, i.e., the ratio of the cellular tip radius to the cellular width α 〉3/2/π≈ 0.3899, which is in good agreement with simulation results. This study also reveals that the cellular tip splitting instability is attributable to weak Gibbs–Thomson energy acting on the interface.展开更多
In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plas...In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plastic solutions of the crack lip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field, to calculate the plastic stress intensity factors, are also developed. Therefore, a complete analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.展开更多
在剖析了V.keilia-Borok研究小组近年来提出的强震预报(Time of increasedprobability)方法的基础上,本文简要地介绍了该方法的原理,采用1970-01-01至1991-12-31的全国地震目录,通过两种途径,用Tip方法对云南省和台湾省发生的7级以上(含...在剖析了V.keilia-Borok研究小组近年来提出的强震预报(Time of increasedprobability)方法的基础上,本文简要地介绍了该方法的原理,采用1970-01-01至1991-12-31的全国地震目录,通过两种途径,用Tip方法对云南省和台湾省发生的7级以上(含7级)强震做了检验地震目录,通过表明Tip方法对强震的预报成功率还是较高的.在此基础上,我们又试验将该方法推广应用于云南省发生的6.5级以上(含6.5级)地震,结果表明适当修改地震流函数后,也可以得到较满意的结果.展开更多
Because of the randomness of many impact factors influencing the dynamic assembly relationship of complex machinery, the reliability analysis of dynamic assembly relationship needs to be accomplished considering the r...Because of the randomness of many impact factors influencing the dynamic assembly relationship of complex machinery, the reliability analysis of dynamic assembly relationship needs to be accomplished considering the randomness from a probabilistic perspective. To improve the accuracy and efficiency of dynamic assembly relationship reliability analysis, the mechanical dynamic assembly reliability（MDAR） theory and a distributed collaborative response surface method（DCRSM） are proposed. The mathematic model of DCRSM is established based on the quadratic response surface function, and verified by the assembly relationship reliability analysis of aeroengine high pressure turbine（HPT） blade-tip radial running clearance（BTRRC）. Through the comparison of the DCRSM, traditional response surface method（RSM） and Monte Carlo Method（MCM）, the results show that the DCRSM is not able to accomplish the computational task which is impossible for the other methods when the number of simulation is more than 100 000 times, but also the computational precision for the DCRSM is basically consistent with the MCM and improved by 0.40-4.63% to the RSM, furthermore, the computational efficiency of DCRSM is up to about 188 times of the MCM and 55 times of the RSM under 10000 times simulations. The DCRSM is demonstrated to be a feasible and effective approach for markedly improving the computational efficiency and accuracy of MDAR analysis. Thus, the proposed research provides the promising theory and method for the MDAR design and optimization, and opens a novel research direction of probabilistic analysis for developing the high-performance and high-reliability of aeroengine.展开更多
Nanoscale sliding contacts of smooth surfaces or between a single asperity and a smooth surface have been widely investigated by molecular dynamics simulations, while there are few studies on the sliding contacts betw...Nanoscale sliding contacts of smooth surfaces or between a single asperity and a smooth surface have been widely investigated by molecular dynamics simulations, while there are few studies on the sliding contacts between two rough surfaces. Actually, the friction of two rough surfaces considering interactions between more asperities should be more realistic. By using multiscale method, friction characteristics of two dimensional nanoscale sliding contacts between rigid multi-asperity tips and elastic textured surfaces are investigated. Four nanoscale textured surfaces with different texture shapes are designed, and six multi-asperity tips composed of cylindrical asperities with different radii are used to slide on the textured surfaces. Friction forces are compared for different tips, and effects of the asperity radii on the friction characteristics are investigated. Average friction forces for all the cases are listed and compared, and effects of texture shapes of the textured surfaces are discussed. The results show that textured surface II has a better structure to reduce friction forces. The multi-asperity tips composed of asperities with R=20r0 （r0=0.227 7 nm） or R=30r0 get higher friction forces compared with other cases, and more atoms of the textured surfaces are taken away by these two tips, which are harmful to reduce friction or wear. For the case of R=10ro, friction forces are also high due to large contact areas, but the sliding processes are stable and few atoms are taken away by the tip. The proposed research considers interactions between more asperities to make the model approach to the real sliding contact problems. The results will help to vary or even control friction characteristics by textured surfaces, or provide references to the design of textured surfaces.展开更多
A compliant tower is modeled as a partially dry, partially tapered, damped Timoshenko beam with the superstructure modeled as an eccentric tip mass, and a non-classical damped boundary at the base. The foundation is m...A compliant tower is modeled as a partially dry, partially tapered, damped Timoshenko beam with the superstructure modeled as an eccentric tip mass, and a non-classical damped boundary at the base. The foundation is modeled as a combination of a linear spring and a torsional spring, along with parallel linear and torsional dampers(Kelvin-Voigt model). The superstructure adds to the kinetic energy of the system without affecting the potential energy, thereby reducing the natural frequencies. The weight of the superstructure acts as an axial compressive load on the beam, reducing its natural frequencies further. The empty space factor due to the truss-type structure of the tower is included. The effect of shear deformation and rotary inertia are included in the vibration analysis; with the non-uniform beam mode-shapes being a weighted sum of the uniform beam mode-shapes satisfying the end condition. The weights are evaluated by the Rayleigh-Ritz(RR) method, and verified using finite element method(FEM). The weight of the superstructure acts as an axial compressive load on the beam. Kelvin-Voigt model of structural damping is included.A part of the structure being underwater, the virtual added inertia is included to calculate the wet natural frequencies. A parametric study is done for various magnitudes of tip mass and various levels of submergence. The computational efficiency and accuracy of the Rayleigh-Ritz method, as compared to the FEA, has been demonstrated. The advantage of using closed-form trial functions is clearly seen in the efficacy of calculating the various energy components in the RR method.展开更多
The displacement fields u_x,u_y at growing crack tip of LYI2-M specimens with double edge cracks are measured using moire method.The experimental singularity fields are compared with GH theoretical field [12-14].The s...The displacement fields u_x,u_y at growing crack tip of LYI2-M specimens with double edge cracks are measured using moire method.The experimental singularity fields are compared with GH theoretical field [12-14].The size and shape of the experimental GH singularity fields are obtained.The error in both the experimental and theoretical evaluations is controlled within ±10%.The experiments show that there is (In A/r)^(+1)singularity dominant around a growing crack tip.The shape of this dominant region ranges from butterfly wing to oblate and circular.Inside GH-field,there is a 3-D deformed damage zone where no GH sin- gularity exists.展开更多
Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenien...Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenient for engineering analyses. A numerical example is presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.展开更多
A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discon...A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.展开更多
Reproducing kernel particle Method （RKPM） is a meshless technology which has proven very useful for solving problems of elastic plastic fracture mechanics. The mode I plastic zone shape at the crack-tip in a work-ha...Reproducing kernel particle Method （RKPM） is a meshless technology which has proven very useful for solving problems of elastic plastic fracture mechanics. The mode I plastic zone shape at the crack-tip in a work-hard ening material is obtained using RKPM. Ramberg-Osgood stress-strain relation is assumed and the crack-tip stress intensity factor （SIF） before and after formation of the plastic zone are examined. To impose the essential boundary conditions, penalty method is used. To construct the shape functions in the vicinity of the crack and crack-tip, both the diffraction and visibility criteria are employed. A comparison between two conventional treatments, visibility and diffraction, to crack discontinuity is conducted. The effects of different dilation parameters on SIF under plane-stress and plane-strain conditions are ~tudied. Results including plastic zone shape are compared with finite element method （FEM） to show the accuracy of RKPM. The main objective is to study the effects of different dilation parameters on SIF under plane stress and plane strain conditions and to obtain the mode I plastic zone shape at the crack-tip in a work hardening material using RKPM.展开更多
The distribution of hydrogen near the notch tip of the austenitic steel 21Cr9Ni9MnN under mode Ⅰ loading has been studid with ion microprobe mass analyzer.Two peaks of hydrogen accumulation,one at notch tip and other...The distribution of hydrogen near the notch tip of the austenitic steel 21Cr9Ni9MnN under mode Ⅰ loading has been studid with ion microprobe mass analyzer.Two peaks of hydrogen accumulation,one at notch tip and other far apart from it,were found.The distribution of hydrostatic stress and plastic strain in front of the notch tip under planar stress condition were calculated using ADINA non-linear finite element method.The calculated result in compari- son with data measured showed that the hydrogen concentration peak at the notch tip is caused by trapping effect of dislocation on dissolved hydrogen,while the other from the hydrostatic stress.The influence of stress intensity factor on the distribution of hydrogen con- centration and the process of hydrogen accumulation near notch tip was also discussed.展开更多
An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. I...An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.展开更多
The crack tip fields are investigated for a cracked functionally graded material （FGM） plate by Reissner＇s linear plate theory with the consideration of the transverse shear deformation generated by bending. The el...The crack tip fields are investigated for a cracked functionally graded material （FGM） plate by Reissner＇s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson＇s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner＇s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner＇s effect when the in-homogeneity parameter approaches zero.展开更多
基金Project supported by the National Basic Research Program of China(Grant No.2011CB610401)the National Natural Science Foundation of China(Grant No.51371151)the Free Research Fund of State Key Laboratory of Solidification Processing,China(Grant No.100-QP-2014)
文摘Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explained by analyzing the interface potential, the tangential energy flux, and the normal energy flux. A rigorous criterion for tip-splitting instability is established analytically, i.e., the ratio of the cellular tip radius to the cellular width α 〉3/2/π≈ 0.3899, which is in good agreement with simulation results. This study also reveals that the cellular tip splitting instability is attributable to weak Gibbs–Thomson energy acting on the interface.
文摘In this paper, a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample. The elastic-plastic solutions of the crack lip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field, to calculate the plastic stress intensity factors, are also developed. Therefore, a complete analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.
文摘在剖析了V.keilia-Borok研究小组近年来提出的强震预报(Time of increasedprobability)方法的基础上,本文简要地介绍了该方法的原理,采用1970-01-01至1991-12-31的全国地震目录,通过两种途径,用Tip方法对云南省和台湾省发生的7级以上(含7级)强震做了检验地震目录,通过表明Tip方法对强震的预报成功率还是较高的.在此基础上,我们又试验将该方法推广应用于云南省发生的6.5级以上(含6.5级)地震,结果表明适当修改地震流函数后,也可以得到较满意的结果.
基金supported by National Natural Science Foundation of China(Grant Nos.51175017,51245027)Innovation Foundation of Beihang University for PhD Graduates,China(Grant No.YWF-12-RBYJ008)Research Fund for the Doctoral Program of Higher Education of China(Grant No.20111102110011)
文摘Because of the randomness of many impact factors influencing the dynamic assembly relationship of complex machinery, the reliability analysis of dynamic assembly relationship needs to be accomplished considering the randomness from a probabilistic perspective. To improve the accuracy and efficiency of dynamic assembly relationship reliability analysis, the mechanical dynamic assembly reliability（MDAR） theory and a distributed collaborative response surface method（DCRSM） are proposed. The mathematic model of DCRSM is established based on the quadratic response surface function, and verified by the assembly relationship reliability analysis of aeroengine high pressure turbine（HPT） blade-tip radial running clearance（BTRRC）. Through the comparison of the DCRSM, traditional response surface method（RSM） and Monte Carlo Method（MCM）, the results show that the DCRSM is not able to accomplish the computational task which is impossible for the other methods when the number of simulation is more than 100 000 times, but also the computational precision for the DCRSM is basically consistent with the MCM and improved by 0.40-4.63% to the RSM, furthermore, the computational efficiency of DCRSM is up to about 188 times of the MCM and 55 times of the RSM under 10000 times simulations. The DCRSM is demonstrated to be a feasible and effective approach for markedly improving the computational efficiency and accuracy of MDAR analysis. Thus, the proposed research provides the promising theory and method for the MDAR design and optimization, and opens a novel research direction of probabilistic analysis for developing the high-performance and high-reliability of aeroengine.
基金supported by National Natural Science Foundation of China(Grant Nos.51205313,50975232)111 Project(Grant No.B13044)Northwestern Polytechnical University Foundation for Fundamental Research,China(Grant No.JC20110249)
文摘Nanoscale sliding contacts of smooth surfaces or between a single asperity and a smooth surface have been widely investigated by molecular dynamics simulations, while there are few studies on the sliding contacts between two rough surfaces. Actually, the friction of two rough surfaces considering interactions between more asperities should be more realistic. By using multiscale method, friction characteristics of two dimensional nanoscale sliding contacts between rigid multi-asperity tips and elastic textured surfaces are investigated. Four nanoscale textured surfaces with different texture shapes are designed, and six multi-asperity tips composed of cylindrical asperities with different radii are used to slide on the textured surfaces. Friction forces are compared for different tips, and effects of the asperity radii on the friction characteristics are investigated. Average friction forces for all the cases are listed and compared, and effects of texture shapes of the textured surfaces are discussed. The results show that textured surface II has a better structure to reduce friction forces. The multi-asperity tips composed of asperities with R=20r0 （r0=0.227 7 nm） or R=30r0 get higher friction forces compared with other cases, and more atoms of the textured surfaces are taken away by these two tips, which are harmful to reduce friction or wear. For the case of R=10ro, friction forces are also high due to large contact areas, but the sliding processes are stable and few atoms are taken away by the tip. The proposed research considers interactions between more asperities to make the model approach to the real sliding contact problems. The results will help to vary or even control friction characteristics by textured surfaces, or provide references to the design of textured surfaces.
文摘A compliant tower is modeled as a partially dry, partially tapered, damped Timoshenko beam with the superstructure modeled as an eccentric tip mass, and a non-classical damped boundary at the base. The foundation is modeled as a combination of a linear spring and a torsional spring, along with parallel linear and torsional dampers(Kelvin-Voigt model). The superstructure adds to the kinetic energy of the system without affecting the potential energy, thereby reducing the natural frequencies. The weight of the superstructure acts as an axial compressive load on the beam, reducing its natural frequencies further. The empty space factor due to the truss-type structure of the tower is included. The effect of shear deformation and rotary inertia are included in the vibration analysis; with the non-uniform beam mode-shapes being a weighted sum of the uniform beam mode-shapes satisfying the end condition. The weights are evaluated by the Rayleigh-Ritz(RR) method, and verified using finite element method(FEM). The weight of the superstructure acts as an axial compressive load on the beam. Kelvin-Voigt model of structural damping is included.A part of the structure being underwater, the virtual added inertia is included to calculate the wet natural frequencies. A parametric study is done for various magnitudes of tip mass and various levels of submergence. The computational efficiency and accuracy of the Rayleigh-Ritz method, as compared to the FEA, has been demonstrated. The advantage of using closed-form trial functions is clearly seen in the efficacy of calculating the various energy components in the RR method.
基金The project suppoted by National Natural Science Foundation of China
文摘The displacement fields u_x,u_y at growing crack tip of LYI2-M specimens with double edge cracks are measured using moire method.The experimental singularity fields are compared with GH theoretical field [12-14].The size and shape of the experimental GH singularity fields are obtained.The error in both the experimental and theoretical evaluations is controlled within ±10%.The experiments show that there is (In A/r)^(+1)singularity dominant around a growing crack tip.The shape of this dominant region ranges from butterfly wing to oblate and circular.Inside GH-field,there is a 3-D deformed damage zone where no GH sin- gularity exists.
文摘Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenient for engineering analyses. A numerical example is presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.
基金Project supported by the National Natural Science Foundation of China (No. 10272037)
文摘A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.
文摘Reproducing kernel particle Method （RKPM） is a meshless technology which has proven very useful for solving problems of elastic plastic fracture mechanics. The mode I plastic zone shape at the crack-tip in a work-hard ening material is obtained using RKPM. Ramberg-Osgood stress-strain relation is assumed and the crack-tip stress intensity factor （SIF） before and after formation of the plastic zone are examined. To impose the essential boundary conditions, penalty method is used. To construct the shape functions in the vicinity of the crack and crack-tip, both the diffraction and visibility criteria are employed. A comparison between two conventional treatments, visibility and diffraction, to crack discontinuity is conducted. The effects of different dilation parameters on SIF under plane-stress and plane-strain conditions are ~tudied. Results including plastic zone shape are compared with finite element method （FEM） to show the accuracy of RKPM. The main objective is to study the effects of different dilation parameters on SIF under plane stress and plane strain conditions and to obtain the mode I plastic zone shape at the crack-tip in a work hardening material using RKPM.
文摘The distribution of hydrogen near the notch tip of the austenitic steel 21Cr9Ni9MnN under mode Ⅰ loading has been studid with ion microprobe mass analyzer.Two peaks of hydrogen accumulation,one at notch tip and other far apart from it,were found.The distribution of hydrostatic stress and plastic strain in front of the notch tip under planar stress condition were calculated using ADINA non-linear finite element method.The calculated result in compari- son with data measured showed that the hydrogen concentration peak at the notch tip is caused by trapping effect of dislocation on dissolved hydrogen,while the other from the hydrostatic stress.The influence of stress intensity factor on the distribution of hydrogen con- centration and the process of hydrogen accumulation near notch tip was also discussed.
基金co-supported by the National High Technology Research and Development Program of China(No.2015AA015303)the National Natural Science Foundation of China(No.11272152)+1 种基金the Aeronautical Science Foundation of China(No.20152752033)the Open Project of Key Laboratory of Aerodynamic Noise Control
文摘An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.
基金supported by the National Natural Science Foundation of China(Nos.90305023 and 11172332)
文摘The crack tip fields are investigated for a cracked functionally graded material （FGM） plate by Reissner＇s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson＇s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner＇s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner＇s effect when the in-homogeneity parameter approaches zero.
