The initial stresses widely exist in elastic materials.While achieving a continuum stress-free configuration through compatible unloading is desirable,mechanical unloading alone frequently proves insufficient,posing c...The initial stresses widely exist in elastic materials.While achieving a continuum stress-free configuration through compatible unloading is desirable,mechanical unloading alone frequently proves insufficient,posing challenges in avoiding virtual stress-free configurations.In this paper,we introduce a novel concept of equivalent temperature variation to counteract the incompatible initial strain.Our focus is on initially stressed cylindrical and spherical elastomers,where we first derive the Saint-Venant,Beltrami-Michell,and Volterra integral conditions in orthogonal curvilinear coordinates using the exterior differential form theory.It is shown that for any given axially or spherically distributed initial stress,an equivalent temperature variation always exists.Furthermore,we propose two innovative initial stress forms based on the steady-state heat conduction.By introducing an equivalent temperature variation,the initial stress can be released through a compatible thermo-mechanical unloading process,offering valuable insights into the constitutive theory of initially stressed elastic materials.展开更多
Aim To study the Lie symmetries and the consered quantities of the holonomic systems with remainder coordinates. Methods Using the invariance of the ordinary differential equations under the infinitesimal transformati...Aim To study the Lie symmetries and the consered quantities of the holonomic systems with remainder coordinates. Methods Using the invariance of the ordinary differential equations under the infinitesimal transformations to establish the determining equations and the restriction equations of the Lie symmetries of the systems. Results and Conclusion the structure equation and the form of conserved quantities were obtained. An example was given to illustrate the application of the result.展开更多
The concise and informative representation of hyperspectral imagery is achieved via the introduced diffusion geometric coordinates derived from nonlinear dimension reduction maps - diffusion maps. The huge-volume high...The concise and informative representation of hyperspectral imagery is achieved via the introduced diffusion geometric coordinates derived from nonlinear dimension reduction maps - diffusion maps. The huge-volume high- dimensional spectral measurements are organized by the affinity graph where each node in this graph only connects to its local neighbors and each edge in this graph represents local similarity information. By normalizing the affinity graph appropriately, the diffusion operator of the underlying hyperspectral imagery is well-defined, which means that the Markov random walk can be simulated on the hyperspectral imagery. Therefore, the diffusion geometric coordinates, derived from the eigenfunctions and the associated eigenvalues of the diffusion operator, can capture the intrinsic geometric information of the hyperspectral imagery well, which gives more enhanced representation results than traditional linear methods, such as principal component analysis based methods. For large-scale full scene hyperspectral imagery, by exploiting the backbone approach, the computation complexity and the memory requirements are acceptable. Experiments also show that selecting suitable symmetrization normalization techniques while forming the diffusion operator is important to hyperspectral imagery representation.展开更多
The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccur...The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.展开更多
A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of th...A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of the model is based on an energy-conservation scheme of finite difference. The most important new feature of the model is that it is a truly compressible ocean model and it is free of the Boussinesq approximations. Thus, the new model is quite different from many existing models in the following ways: 1) the exact form of mass conservation, 2) the in-situ instantaneous pressure and the UNESCO equation of state to calculate density, 3) the in-situ density in the momentum. equations, 4) finite difference schemes that conserve the total energy. Initial tests showed that the model code runs smoothly, and it is quite stable. The quasi-steady circulation patterns generated by the new model compare well with existing models, but the time evolution of the new model seems different from some existing models. Thus, the non-Boussinesq models may provide more accurate information for climate study and satellite observations.展开更多
The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody sys...The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody systems. However, the formulation has so many principles in choosing the generalized coordinates that it hinders the implementation of modeling automation, A first order direct sensitivity analysis approach to multibody systems formulated with novel natural coordinates is presented. Firstly, a new selection method for natural coordinate is developed. The method introduces 12 coordinates to describe the position and orientation of a spatial object. On the basis of the proposed natural coordinates, rigid constraint conditions, the basic constraint elements as well as the initial conditions for the governing equations are derived. Considering the characteristics of the governing equations, the newly proposed generalized-ct integration method is used and the corresponding algorithm flowchart is discussed. The objective function, the detailed analysis process of first order direct sensitivity analysis and related solving strategy are provided based on the previous modeling system Finally, in order to verify the validity and accuracy of the method presented, the sensitivity analysis of a planar spinner-slider mechanism and a spatial crank-slider mechanism are conducted. The test results agree well with that of the finite difference method, and the maximum absolute deviation of the results is less than 3%. The proposed approach is not only convenient for automatic modeling, but also helpful for the reduction of the complexity of sensitivity analysis, which provides a practical and effective way to obtain sensitivity for the optimization problems of multibody systems.展开更多
Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with...Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.展开更多
A reconstructed edge-based smoothed triangular element, which is incorporated with the discrete shear gap (DSG) method, is formulated based on the global coordinate for analysis of Reissner-Mindlin plates. A symbolic ...A reconstructed edge-based smoothed triangular element, which is incorporated with the discrete shear gap (DSG) method, is formulated based on the global coordinate for analysis of Reissner-Mindlin plates. A symbolic integration combined with the smoothing technique is implemented to calculate the smoothed finite element matrices, which is integrated along the boundaries of each smoothing cell. Numerical results show that the proposed element is free from shear locking, and its results are in good agreement with the exact solutions, even for very thin plates with extremely distorted elements. The proposed element gives more accurate results than the original DSG element without smoothing, and it can be taken as an alternative element for analysis of Reissner-Mindlin plates. The prominent feature of the present element is that the integration scheme is unified in the smoothed form for all of the finite element matrices.展开更多
The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the...The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the feasibility of reproducing these velocity fields with computer models. If computer modeling were to prove successful in this context, it could become a useful prediction technique and research tool to enhance our understanding of natural river dynamics. A 3-D k-E turbulence hydrodynamic model in curvilinear coordinates is established to simulate the overhank flow. The bodyfitted coordinate is adopted in the horizontal plane, the part grid is adopted in the vertical direction, and the wall-function method is employed to simulate the bed resistance. The model is applied to the simulation of the meandering channel with straight flood plain banks, and the main velocities and secondary velocities for both the longitudinal and cross sections are presented. Comparison and analysis show that the results of simulation are fit to reflect the results of experiment. These results show the application value of the model to 3D overhank flow.展开更多
A<110>/2 screw dislocation is commonly dissociated into two <112>/6 Shockley partial dislocations on{111} planes in face-centered cubic metals.As the two partials are not purely screw,different mechanisms ...A<110>/2 screw dislocation is commonly dissociated into two <112>/6 Shockley partial dislocations on{111} planes in face-centered cubic metals.As the two partials are not purely screw,different mechanisms of cross-slip could take place,depending on the stacking fault energy,applied stress and tempe rature.It is crucial to classify the mechanisms of cross-slip because each mechanism possesses its own reaction path with a special activation process.In this work,molecular dynamics simulations have been performed systematically to explore the cross-slip mechanism under different stresses and temperatures in three different metals Ag,Cu and Ni that have different stacking fault energies of 17.8,44.4 and 126.8 mJ/m^2,re spectively.In Ag and Cu with low stacking fault energy,it is observed that the cross-slip mechanism of screw dislocations changes from the Fleischer obtuse angle(FLOA),to the Friedel-Escaig(FE),and then to the FL acute angle(FLAA) at low temperatures,with increasing the applied stress.However,when the temperature increases,the FE mechanism gradually becomes dominant,while the FLAA only occurs at the high stress region.In particular,the FLOA has not been observed in Ni because of its high stacking fault energy.展开更多
This investigation is intended to develop a computer procedure for the integration of NURBS geometry and the rational absolute nodal coordinate formulation （RANCF） finite element analysis. A linear transformation is...This investigation is intended to develop a computer procedure for the integration of NURBS geometry and the rational absolute nodal coordinate formulation （RANCF） finite element analysis. A linear transformation is given that can be used to convert the NURBS curve to RANCF cable element mesh retaining the same geometry and the same degree of continuity, including the discussion of continuity control and mesh refinement. The green strain tensor is used to establish the nonlinear dynamic equations with numerical examples to demonstrate the use of the procedure in the dynamic analysis of flexible bodies.展开更多
A crystal is a highly organized arrangement of atoms in a solid, wherein a unit cell is periodically repeated to form the crystal pattern. A unit cell is composed of atoms that are connected to some of their first nei...A crystal is a highly organized arrangement of atoms in a solid, wherein a unit cell is periodically repeated to form the crystal pattern. A unit cell is composed of atoms that are connected to some of their first neighbors by chemical bonds. A recent rule, entitled the even-odd rule, introduced a new way to calculate the number of covalent bonds around an atom. It states that around an uncharged atom, the number of bonds and the number of electrons have the same parity. In the case of a charged atom on the contrary, both numbers have different parity. The aim of the present paper is to challenge the even-odd rule on chemical bonds in well-known crystal structures. According to the rule, atoms are supposed to be bonded exclusively through single-covalent bonds. A distinctive criterion, only applicable to crystals, states that atoms cannot build more than 8 chemical bonds, as opposed to the classical model, where each atom in a crystal is connected to every first neighbor without limitation. Electrical charges can be assigned to specific atoms in order to compensate for extra or missing bonds. More specifically the article considers di-atomic body-centered-cubic, tetra-atomic and dodeca-atomic single-face-centered-cubic crystals. In body-centered crystals, atoms are interconnected by 8 covalent bonds. In face-centered crystal, the unit cell contains 4 or 12 atoms. For di-element crystals, the total number of bonds for both elements is found to be identical. The neutrality of the unit cell is obtained with an opposite charge on the nearest or second-nearest neighbor. To conclude, the even-odd rule is applicable to a wide number of compounds in known cubic structures and the number of chemical bonds per atom is not related to the valence of the elements in the periodic table.展开更多
A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
In this paper, unit moving trihedron is first constructed for a point on the surface of a revolution ellipsoid. Via translation, the origin of the trihedron coincides with that of Cartesian coordinates established at ...In this paper, unit moving trihedron is first constructed for a point on the surface of a revolution ellipsoid. Via translation, the origin of the trihedron coincides with that of Cartesian coordinates established at the center of the ellipsoid, and then through two coordinate rotations, the trihedron completely coincides with the Cartesian coordinates. Transformation formulae between the moving trihedron and unit Cartesian coordinate frameworks as well as transformation of point displacement between two unit coordinate frameworks are presented. Based on the above transformation formulae between two different coordinate frameworks, due to the fact that the displacement and moving trihedron of the point are both functions of the geodetic coordinates, components in the corresponding axis for differential of displacement vector and geodetic curves arc differential at the point in geodetic system can be obtained through complicated derivation. Displacement gradient matrix at the point in geodetic system is also given. Finally, expressions of strain and rotation tensor in geodetic coordinates are presented. Geometric meanings of the rotation tensor are explained in detail. The intrinsic relationship between strain tensors of sphere and ellipsoid are also discussed.展开更多
This paper introduces the canonical coordinates method to obtain the first integral of a single-degree freedom constraint mechanical system that contains conserva-tive and non-conservative constraint homonomic systems...This paper introduces the canonical coordinates method to obtain the first integral of a single-degree freedom constraint mechanical system that contains conserva-tive and non-conservative constraint homonomic systems. The definition and properties of canonical coordinates are introduced. The relation between Lie point symmetries and the canonical coordinates of the constraint mechanical system are expressed. By this re-lation, the canonical coordinates can be obtained. Properties of the canonical coordinates and the Lie symmetry theory are used to seek the first integrals of constraint mechanical system. Three examples are used to show applications of the results.展开更多
For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries, a numerical model is developed in curvilinear coordinates. In the model, the Boussinesq-type equations including the dissi...For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries, a numerical model is developed in curvilinear coordinates. In the model, the Boussinesq-type equations including the dissipation terms are em- ployed as the governing equations. In the present model, the dependent variables of the transformed equations are the free surface elevation and the utility velocity variables, instead of the usual primitive velocity variables. The introduction of utility velocity variables which are the products of the contravariant components of the velocity vector and the Jacobi ma- trix can make the transformed equations relatively concise, the treatment of lateral boundary conditions easier and the de- velopment of the program simpler. The predictor-corrector method and five-point finite-difference scheme are employed to discretize the time derivatives and the spatial ones, respectively. The numerical model is tested for three cases. It is found that the numerical results are in good agreement with the analytical results and experimental data.展开更多
An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-...An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.展开更多
Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of p...Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of perfectly matched layer (PML), MWI PML absorbing boundary condition (ABC) algorithm was deduced in 2D cylindrical coordinates. Numerical experiments were done to investigate the validity of MWI and its application in cylindrical coordinates FDTD algorithm. The results showed that MWI in cylindrical coordinates can be used to accurately calculate the numerical reflection error caused by different mesh increments in non uniform FDTD. MWI can also provide theoretical criterion to define the permitted variable range of mesh dimension. MWI PML ABC is easy to be applied and reduces low numerical reflection, which only causes a little higher reflection error compared with Teixeira's PML.展开更多
A scheme of space time conservation (STC) based on the method of space time conservation element and solution element (CE/SE) is represented in a nonorthogonal curvilinear coordinate system. The corresponding initia...A scheme of space time conservation (STC) based on the method of space time conservation element and solution element (CE/SE) is represented in a nonorthogonal curvilinear coordinate system. The corresponding initial and boundary conditions are discussed. It is seen that in the nonorthogonal coordinates the scheme maintains the advantages of the STC method, and is noted for its simple structure, clear physical meaning, rapid calculation and high accuracy. It is easy to extend to the multidimensional flow. The numerical results for a 2D Euler equation show good agreement with those from other computational methods and the experiment.展开更多
This article shows that in spherical polar coordinates, some noncentral separable potentials have super-symmetry and shape invariance in the r and θ dimensions, we choose Hartmann potential and ring-shaped oscillator...This article shows that in spherical polar coordinates, some noncentral separable potentials have super-symmetry and shape invariance in the r and θ dimensions, we choose Hartmann potential and ring-shaped oscillator astwo important examples, thus in principle the energy eigenvalues and energy eigenfunctions of such the potentials in ther and θ dimensions can be obtained by the method of supersymmetric quantum mechanics. Here we use an alternativemethod to get the required results.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.12241205 and 12032019)the National Key Research and Development Program of China(No.2022YFA1203200)the Strategic Priority Research Program of Chinese Academy of Sciences(Nos.XDB0620101 and XDB0620103)。
文摘The initial stresses widely exist in elastic materials.While achieving a continuum stress-free configuration through compatible unloading is desirable,mechanical unloading alone frequently proves insufficient,posing challenges in avoiding virtual stress-free configurations.In this paper,we introduce a novel concept of equivalent temperature variation to counteract the incompatible initial strain.Our focus is on initially stressed cylindrical and spherical elastomers,where we first derive the Saint-Venant,Beltrami-Michell,and Volterra integral conditions in orthogonal curvilinear coordinates using the exterior differential form theory.It is shown that for any given axially or spherically distributed initial stress,an equivalent temperature variation always exists.Furthermore,we propose two innovative initial stress forms based on the steady-state heat conduction.By introducing an equivalent temperature variation,the initial stress can be released through a compatible thermo-mechanical unloading process,offering valuable insights into the constitutive theory of initially stressed elastic materials.
文摘Aim To study the Lie symmetries and the consered quantities of the holonomic systems with remainder coordinates. Methods Using the invariance of the ordinary differential equations under the infinitesimal transformations to establish the determining equations and the restriction equations of the Lie symmetries of the systems. Results and Conclusion the structure equation and the form of conserved quantities were obtained. An example was given to illustrate the application of the result.
基金The National Key Technologies R & D Program during the 11th Five-Year Plan Period (No.2006BAB15B01)
文摘The concise and informative representation of hyperspectral imagery is achieved via the introduced diffusion geometric coordinates derived from nonlinear dimension reduction maps - diffusion maps. The huge-volume high- dimensional spectral measurements are organized by the affinity graph where each node in this graph only connects to its local neighbors and each edge in this graph represents local similarity information. By normalizing the affinity graph appropriately, the diffusion operator of the underlying hyperspectral imagery is well-defined, which means that the Markov random walk can be simulated on the hyperspectral imagery. Therefore, the diffusion geometric coordinates, derived from the eigenfunctions and the associated eigenvalues of the diffusion operator, can capture the intrinsic geometric information of the hyperspectral imagery well, which gives more enhanced representation results than traditional linear methods, such as principal component analysis based methods. For large-scale full scene hyperspectral imagery, by exploiting the backbone approach, the computation complexity and the memory requirements are acceptable. Experiments also show that selecting suitable symmetrization normalization techniques while forming the diffusion operator is important to hyperspectral imagery representation.
基金supported by the National Basic Research Program of China ( Grant No.2006CB403302)the National Natural Science Foundation of China (Grant Nos .50839001 and 50709004)the Scientific Research Foundation of the Higher Education Institutions of Liaoning Province (Grant No.2006T018)
文摘The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.
文摘A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of the model is based on an energy-conservation scheme of finite difference. The most important new feature of the model is that it is a truly compressible ocean model and it is free of the Boussinesq approximations. Thus, the new model is quite different from many existing models in the following ways: 1) the exact form of mass conservation, 2) the in-situ instantaneous pressure and the UNESCO equation of state to calculate density, 3) the in-situ density in the momentum. equations, 4) finite difference schemes that conserve the total energy. Initial tests showed that the model code runs smoothly, and it is quite stable. The quasi-steady circulation patterns generated by the new model compare well with existing models, but the time evolution of the new model seems different from some existing models. Thus, the non-Boussinesq models may provide more accurate information for climate study and satellite observations.
基金supported by National Defense Pre-research Foundation of China during the 12th Five-Year Plan Period(Grant No.51036050107)
文摘The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody systems. However, the formulation has so many principles in choosing the generalized coordinates that it hinders the implementation of modeling automation, A first order direct sensitivity analysis approach to multibody systems formulated with novel natural coordinates is presented. Firstly, a new selection method for natural coordinate is developed. The method introduces 12 coordinates to describe the position and orientation of a spatial object. On the basis of the proposed natural coordinates, rigid constraint conditions, the basic constraint elements as well as the initial conditions for the governing equations are derived. Considering the characteristics of the governing equations, the newly proposed generalized-ct integration method is used and the corresponding algorithm flowchart is discussed. The objective function, the detailed analysis process of first order direct sensitivity analysis and related solving strategy are provided based on the previous modeling system Finally, in order to verify the validity and accuracy of the method presented, the sensitivity analysis of a planar spinner-slider mechanism and a spatial crank-slider mechanism are conducted. The test results agree well with that of the finite difference method, and the maximum absolute deviation of the results is less than 3%. The proposed approach is not only convenient for automatic modeling, but also helpful for the reduction of the complexity of sensitivity analysis, which provides a practical and effective way to obtain sensitivity for the optimization problems of multibody systems.
基金supported by the National Natural Science Foundation of China (Grant Nos. 50839001 and 50979036)
文摘Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.
基金supported by the National Natural Science Foundation of China (Grants 11272118, 11372106)Fundamental Research Fund of the Central Universities (Grant 227201401203)
文摘A reconstructed edge-based smoothed triangular element, which is incorporated with the discrete shear gap (DSG) method, is formulated based on the global coordinate for analysis of Reissner-Mindlin plates. A symbolic integration combined with the smoothing technique is implemented to calculate the smoothed finite element matrices, which is integrated along the boundaries of each smoothing cell. Numerical results show that the proposed element is free from shear locking, and its results are in good agreement with the exact solutions, even for very thin plates with extremely distorted elements. The proposed element gives more accurate results than the original DSG element without smoothing, and it can be taken as an alternative element for analysis of Reissner-Mindlin plates. The prominent feature of the present element is that the integration scheme is unified in the smoothed form for all of the finite element matrices.
文摘The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the feasibility of reproducing these velocity fields with computer models. If computer modeling were to prove successful in this context, it could become a useful prediction technique and research tool to enhance our understanding of natural river dynamics. A 3-D k-E turbulence hydrodynamic model in curvilinear coordinates is established to simulate the overhank flow. The bodyfitted coordinate is adopted in the horizontal plane, the part grid is adopted in the vertical direction, and the wall-function method is employed to simulate the bed resistance. The model is applied to the simulation of the meandering channel with straight flood plain banks, and the main velocities and secondary velocities for both the longitudinal and cross sections are presented. Comparison and analysis show that the results of simulation are fit to reflect the results of experiment. These results show the application value of the model to 3D overhank flow.
基金financially supported by the Program of “One Hundred Talented People” of the Chinese Academy of Sciences (JBY)the National Natural Science Foundation of China (Nos. 51871223, 51771206, and 51790482)。
文摘A<110>/2 screw dislocation is commonly dissociated into two <112>/6 Shockley partial dislocations on{111} planes in face-centered cubic metals.As the two partials are not purely screw,different mechanisms of cross-slip could take place,depending on the stacking fault energy,applied stress and tempe rature.It is crucial to classify the mechanisms of cross-slip because each mechanism possesses its own reaction path with a special activation process.In this work,molecular dynamics simulations have been performed systematically to explore the cross-slip mechanism under different stresses and temperatures in three different metals Ag,Cu and Ni that have different stacking fault energies of 17.8,44.4 and 126.8 mJ/m^2,re spectively.In Ag and Cu with low stacking fault energy,it is observed that the cross-slip mechanism of screw dislocations changes from the Fleischer obtuse angle(FLOA),to the Friedel-Escaig(FE),and then to the FL acute angle(FLAA) at low temperatures,with increasing the applied stress.However,when the temperature increases,the FE mechanism gradually becomes dominant,while the FLAA only occurs at the high stress region.In particular,the FLOA has not been observed in Ni because of its high stacking fault energy.
基金supported by the National Natural Science Foundation of China(No.11172076)the Science and Technology Innovation Talent Foundation of Harbin(No.2012RFLXG020)
文摘This investigation is intended to develop a computer procedure for the integration of NURBS geometry and the rational absolute nodal coordinate formulation （RANCF） finite element analysis. A linear transformation is given that can be used to convert the NURBS curve to RANCF cable element mesh retaining the same geometry and the same degree of continuity, including the discussion of continuity control and mesh refinement. The green strain tensor is used to establish the nonlinear dynamic equations with numerical examples to demonstrate the use of the procedure in the dynamic analysis of flexible bodies.
文摘A crystal is a highly organized arrangement of atoms in a solid, wherein a unit cell is periodically repeated to form the crystal pattern. A unit cell is composed of atoms that are connected to some of their first neighbors by chemical bonds. A recent rule, entitled the even-odd rule, introduced a new way to calculate the number of covalent bonds around an atom. It states that around an uncharged atom, the number of bonds and the number of electrons have the same parity. In the case of a charged atom on the contrary, both numbers have different parity. The aim of the present paper is to challenge the even-odd rule on chemical bonds in well-known crystal structures. According to the rule, atoms are supposed to be bonded exclusively through single-covalent bonds. A distinctive criterion, only applicable to crystals, states that atoms cannot build more than 8 chemical bonds, as opposed to the classical model, where each atom in a crystal is connected to every first neighbor without limitation. Electrical charges can be assigned to specific atoms in order to compensate for extra or missing bonds. More specifically the article considers di-atomic body-centered-cubic, tetra-atomic and dodeca-atomic single-face-centered-cubic crystals. In body-centered crystals, atoms are interconnected by 8 covalent bonds. In face-centered crystal, the unit cell contains 4 or 12 atoms. For di-element crystals, the total number of bonds for both elements is found to be identical. The neutrality of the unit cell is obtained with an opposite charge on the nearest or second-nearest neighbor. To conclude, the even-odd rule is applicable to a wide number of compounds in known cubic structures and the number of chemical bonds per atom is not related to the valence of the elements in the periodic table.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
文摘In this paper, unit moving trihedron is first constructed for a point on the surface of a revolution ellipsoid. Via translation, the origin of the trihedron coincides with that of Cartesian coordinates established at the center of the ellipsoid, and then through two coordinate rotations, the trihedron completely coincides with the Cartesian coordinates. Transformation formulae between the moving trihedron and unit Cartesian coordinate frameworks as well as transformation of point displacement between two unit coordinate frameworks are presented. Based on the above transformation formulae between two different coordinate frameworks, due to the fact that the displacement and moving trihedron of the point are both functions of the geodetic coordinates, components in the corresponding axis for differential of displacement vector and geodetic curves arc differential at the point in geodetic system can be obtained through complicated derivation. Displacement gradient matrix at the point in geodetic system is also given. Finally, expressions of strain and rotation tensor in geodetic coordinates are presented. Geometric meanings of the rotation tensor are explained in detail. The intrinsic relationship between strain tensors of sphere and ellipsoid are also discussed.
基金Project supported by the National Natural Science Foundation of China(Nos.11072218 and 11272287)the Program for Changjiang Scholars and Innovative Research Team in University(PCSIRT)(No.IRT13097)
文摘This paper introduces the canonical coordinates method to obtain the first integral of a single-degree freedom constraint mechanical system that contains conserva-tive and non-conservative constraint homonomic systems. The definition and properties of canonical coordinates are introduced. The relation between Lie point symmetries and the canonical coordinates of the constraint mechanical system are expressed. By this re-lation, the canonical coordinates can be obtained. Properties of the canonical coordinates and the Lie symmetry theory are used to seek the first integrals of constraint mechanical system. Three examples are used to show applications of the results.
基金supported by the National Natural Science Foundation of China (Grant Nos .51079082 and 40676053)State Key Laboratory of Ocean Engineering ( Grant Nos . GKZD010012, GP010818 and GKZD010024)
文摘For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries, a numerical model is developed in curvilinear coordinates. In the model, the Boussinesq-type equations including the dissipation terms are em- ployed as the governing equations. In the present model, the dependent variables of the transformed equations are the free surface elevation and the utility velocity variables, instead of the usual primitive velocity variables. The introduction of utility velocity variables which are the products of the contravariant components of the velocity vector and the Jacobi ma- trix can make the transformed equations relatively concise, the treatment of lateral boundary conditions easier and the de- velopment of the program simpler. The predictor-corrector method and five-point finite-difference scheme are employed to discretize the time derivatives and the spatial ones, respectively. The numerical model is tested for three cases. It is found that the numerical results are in good agreement with the analytical results and experimental data.
文摘An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.
文摘Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of perfectly matched layer (PML), MWI PML absorbing boundary condition (ABC) algorithm was deduced in 2D cylindrical coordinates. Numerical experiments were done to investigate the validity of MWI and its application in cylindrical coordinates FDTD algorithm. The results showed that MWI in cylindrical coordinates can be used to accurately calculate the numerical reflection error caused by different mesh increments in non uniform FDTD. MWI can also provide theoretical criterion to define the permitted variable range of mesh dimension. MWI PML ABC is easy to be applied and reduces low numerical reflection, which only causes a little higher reflection error compared with Teixeira's PML.
文摘A scheme of space time conservation (STC) based on the method of space time conservation element and solution element (CE/SE) is represented in a nonorthogonal curvilinear coordinate system. The corresponding initial and boundary conditions are discussed. It is seen that in the nonorthogonal coordinates the scheme maintains the advantages of the STC method, and is noted for its simple structure, clear physical meaning, rapid calculation and high accuracy. It is easy to extend to the multidimensional flow. The numerical results for a 2D Euler equation show good agreement with those from other computational methods and the experiment.
文摘This article shows that in spherical polar coordinates, some noncentral separable potentials have super-symmetry and shape invariance in the r and θ dimensions, we choose Hartmann potential and ring-shaped oscillator astwo important examples, thus in principle the energy eigenvalues and energy eigenfunctions of such the potentials in ther and θ dimensions can be obtained by the method of supersymmetric quantum mechanics. Here we use an alternativemethod to get the required results.
