Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficul...Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.展开更多
Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the resul...Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the results should be given together when the measurement result is given. Nowadays most researches on straightness focus on error calculation and only several research projects evaluate the measurement uncertainty based on ＂The Guide to the Expression of Uncertainty in Measurement（GUM）＂. In order to compute spatial straightness error（SSE） accurately and rapidly and overcome the limitations of GUM, a quasi particle swarm optimization（QPSO） is proposed to solve the minimum zone SSE and Monte Carlo Method（MCM） is developed to estimate the measurement uncertainty. The mathematical model of minimum zone SSE is formulated. In QPSO quasi-random sequences are applied to the generation of the initial position and velocity of particles and their velocities are modified by the constriction factor approach. The flow of measurement uncertainty evaluation based on MCM is proposed, where the heart is repeatedly sampling from the probability density function（PDF） for every input quantity and evaluating the model in each case. The minimum zone SSE of a shaft measured on a Coordinate Measuring Machine（CMM） is calculated by QPSO and the measurement uncertainty is evaluated by MCM on the basis of analyzing the uncertainty contributors. The results show that the uncertainty directly influences the product judgment result. Therefore it is scientific and reasonable to consider the influence of the uncertainty in judging whether the parts are accepted or rejected, especially for those located in the uncertainty zone. The proposed method is especially suitable when the PDF of the measurand cannot adequately be approximated by a Gaussian distribution or a scaled and shifted t-distribution and the measurement model is non-linear.展开更多
A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structur...A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structural strength in engineering applications.First,a topology optimization model is established for a lightweight structure with element stress as constraints.Second,the stress globalization method is adopted to convert local stress constraints into strain energy constraints,which overcomes the difficulties caused by local stress constraints,such as model establishment,sensitivity analysis,and massive solution calculations.Third,the sensitivity of the objective function and constraint function is analyzed,and the method of moving asymptotes is employed to solve the optimization model.In addition,the additive hyperelasticity technique is utilized to solve the numerical instability induced by structures undergoing large deformation.Numerical examples are given to validate the feasibility of the proposed method.The method provides a significant reference for geometrically nonlinear optimization design.展开更多
This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically so...This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization （SIMP） scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes （MMA） belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.展开更多
A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of ...A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of revolutions. All of the interface integral items participating in forming the potential and complementary energy functionals of the composite materials are expressed in terms of intrinsic quantities of the ellipsoidal of revolutions. Based on this, the upper and the lower bound for the effective elastic moduli of the composite materials with inclusions described above have been derived. Under three limiting conditions of sphere, disk and needle shaped inclusions, the results of this paper will return to the bounds obtained by Hashin([6]) (1992).展开更多
Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterativ...Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.展开更多
Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume cont...Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.展开更多
The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading.Meanwhile,the application of critical plane method in the prediction of multiaxial fatigu...The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading.Meanwhile,the application of critical plane method in the prediction of multiaxial fatigue life has made certain progress.According to the law of thermodynamics,a new damage evolution equation is developed in the present study to predict the fatigue life of geometrically discontinuous structure under tension-torsion loading based on damage mechanics and the critical plane method.The essence of this approach is tha t the st rain parame ter of the uniaxial nonlinear fatigue damage model is replaced with the equivalent strain,which consists of the releva nt parame ters of the critical plane.However,it is difficult to calculate the stress-strain status and the critical plane position of geometrically dis?continuous structure by theoretical methods because of the existence of stress concentration and the multiaxial nonproportional characteristics.Therefore,a new numerical simulation method is proposed to determine the critical plane of geometrically discontinuous structure under multiaxial loading by means of the finite element method and MATLAB software.The fatigue life of notched specimens subjected to combined bending and torsion is predicted using the proposed met hod,and the result is compared with t hose from the experimen ts and the Manson-Cfiffin law.The comparisons show that the proposed method is superior to the Manson-Coffin law and is capable of reproducing the experimental results reasonably when the geometry of the structure is complex.It completely meets the needs of engineering practice.展开更多
The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was e...The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first_order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C_shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.展开更多
This paper describes geometrical essentials of some iteration methods （e.g. Newton iteration, secant line method, etc.） for solving nonlinear equations and advances some geometrical methods of iteration that are fle...This paper describes geometrical essentials of some iteration methods （e.g. Newton iteration, secant line method, etc.） for solving nonlinear equations and advances some geometrical methods of iteration that are flexible and efficient.展开更多
In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical featu...In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.展开更多
This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element couplin...This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.展开更多
Micro cracking in solid material is a phenomenon of solid plane bifurcation with incompatible micro-rotation. Therefore it is important to analysis the micro-rotation field in a strained material in fractute mechanics...Micro cracking in solid material is a phenomenon of solid plane bifurcation with incompatible micro-rotation. Therefore it is important to analysis the micro-rotation field in a strained material in fractute mechanics. On the basis of strain-rotation (S-R) decomposition theorem, the geometric criterion of cracking has been established. In order to verify its validity, using grating experimental method the thin compact tension specimens were investigated. The analysis results show a good agreement between the experimental and geometric criterion.展开更多
This paper introduced a bootstrap method called truncated geometric bootstrap method for time series stationary process. We estimate the parameters of a geometric distribution which has been truncated as a probability...This paper introduced a bootstrap method called truncated geometric bootstrap method for time series stationary process. We estimate the parameters of a geometric distribution which has been truncated as a probability model for the bootstrap algorithm. This probability model was used in resampling blocks of random length, where the length of each blocks has a truncated geometric distribution. The method was able to determine the block sizes b and probability p attached to its random selections. The mean and variance were estimated for the truncated geometric distribution and the bootstrap algorithm developed based on the proposed probability model.展开更多
In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending ...In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.展开更多
The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consis...The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.展开更多
We present in this paper a numerical study of the validity limit of the geometrical optics approximation compared with a differential method which is established according to rigorous formalisms based on the electroma...We present in this paper a numerical study of the validity limit of the geometrical optics approximation compared with a differential method which is established according to rigorous formalisms based on the electromagnetic theory. The precedent studies show that this method is adapted to the study of diffraction by periodic rough surfaces. We determine by two methods the emissivity of gold and tungsten for surfaces with a rectangular or sinusoidal profile, for a wavelength equal to 0.55 microns. The monochromatic directional emissivity of these surfaces clearly depends on the angle of incidence, the surface profile, height, period and the nature of the material. We perform our calculations by a method of coupled wave analysis (CWA) and a geometric optics method (GOA). The latter method is theoretically valid only when the dimensions of the cavities are very large compared to the wavelength, while the CWA is theoretically correct whatever these dimensions. The main purpose of this work is to investigate the validity limit of GOA compared with CWA. The obtained results for a fixed height of the grating, allowed us to delimit the validity domain of the optic geometrical approximation for the treated cases. Finally, the agreement between the emissivity calculated by the differential method and that given on the basis of the homogenization theory is satisfactory when the period is much smaller than the wavelength.展开更多
This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were di...This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were discretized by Finite Volume Method(FVM)with second-order approximation scheme and deferred correction.Moreover,the algebraic method and the differential method were used to generate the non-orthogonal structured curvilinear grids.Furthermore,the influence of some parameters of geometric multigrid method,as well as lexicographical Gauss–Seidel(Lex-GS),η-line Gauss–Seidel(η-line-GS),Modified Strongly Implicit(MSI)and modified incomplete LU decomposition(MILU)solvers on the Central Processing Unit(CPU)time was investigated.Therefore,several parameters of multigrid method and solvers were tested for the problem,with the use of nonorthogonal structured curvilinear grids and multigrid method,resulting in an algorithm with the combination that achieved the best results and CPU time.The geometric multigrid method with Full Approximation Scheme(FAS),V-cycle and standard coarsening ratio for this problem were utilized.This article shows how to calculate the coordinates transformation metrics in the coarser grids.Results show that the MSI and MILU solvers are the most efficient.Moreover,theMSI solver is faster thanMILU for both grids generators;and the solutions are more accurate for the Burgers problem with grids generated using elliptic equations.展开更多
The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quad...The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.展开更多
A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) n...A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) not only the stress fields of a circular disk can be displayed, but the u -and v -displacement fields can also be produced.展开更多
基金Project supported by the National 973 Program (No.2004CB719402), the National Natural Science Foundation of China (No. 10372030)the Open Research Projects supported by the Project Fund of the Hubei Province Key Lab of Mechanical Transmission & Manufacturing Engineering Wuhan University of Science & Technology (No.2003A16).
文摘Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.
基金supported by National Natural Science Foundation of China (Grant No. 51075198)Jiangsu Provincial Natural Science Foundation of China (Grant No. BK2010479)+2 种基金Innovation Research of Nanjing Institute of Technology, China (Grant No. CKJ20100008)Jiangsu Provincial Foundation of 333 Talents Engineering of ChinaJiangsu Provincial Foundation of Six Talented Peak of China
文摘Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the results should be given together when the measurement result is given. Nowadays most researches on straightness focus on error calculation and only several research projects evaluate the measurement uncertainty based on ＂The Guide to the Expression of Uncertainty in Measurement（GUM）＂. In order to compute spatial straightness error（SSE） accurately and rapidly and overcome the limitations of GUM, a quasi particle swarm optimization（QPSO） is proposed to solve the minimum zone SSE and Monte Carlo Method（MCM） is developed to estimate the measurement uncertainty. The mathematical model of minimum zone SSE is formulated. In QPSO quasi-random sequences are applied to the generation of the initial position and velocity of particles and their velocities are modified by the constriction factor approach. The flow of measurement uncertainty evaluation based on MCM is proposed, where the heart is repeatedly sampling from the probability density function（PDF） for every input quantity and evaluating the model in each case. The minimum zone SSE of a shaft measured on a Coordinate Measuring Machine（CMM） is calculated by QPSO and the measurement uncertainty is evaluated by MCM on the basis of analyzing the uncertainty contributors. The results show that the uncertainty directly influences the product judgment result. Therefore it is scientific and reasonable to consider the influence of the uncertainty in judging whether the parts are accepted or rejected, especially for those located in the uncertainty zone. The proposed method is especially suitable when the PDF of the measurand cannot adequately be approximated by a Gaussian distribution or a scaled and shifted t-distribution and the measurement model is non-linear.
基金This work was supported by the National Natural Science Foundation of China(11872080)Beijing Natural Science Foundation(3192005)。
文摘A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structural strength in engineering applications.First,a topology optimization model is established for a lightweight structure with element stress as constraints.Second,the stress globalization method is adopted to convert local stress constraints into strain energy constraints,which overcomes the difficulties caused by local stress constraints,such as model establishment,sensitivity analysis,and massive solution calculations.Third,the sensitivity of the objective function and constraint function is analyzed,and the method of moving asymptotes is employed to solve the optimization model.In addition,the additive hyperelasticity technique is utilized to solve the numerical instability induced by structures undergoing large deformation.Numerical examples are given to validate the feasibility of the proposed method.The method provides a significant reference for geometrically nonlinear optimization design.
基金the National '973' Key Fundamental Research Projects of China(No.2003CB716207)the National '863' High-Tech Development Projects of China(No.2006AA04Z162)also the Australian Research Council(No.ARC-DP0666683).
文摘This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization （SIMP） scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes （MMA） belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.
文摘A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of revolutions. All of the interface integral items participating in forming the potential and complementary energy functionals of the composite materials are expressed in terms of intrinsic quantities of the ellipsoidal of revolutions. Based on this, the upper and the lower bound for the effective elastic moduli of the composite materials with inclusions described above have been derived. Under three limiting conditions of sphere, disk and needle shaped inclusions, the results of this paper will return to the bounds obtained by Hashin([6]) (1992).
基金Supported by the National Natural Science Foundation of China(61272300)
文摘Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.
基金Projects(51621006,51779251)supported by the National Natural Science Foundation of China。
文摘Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.
基金the National Natural Science Foundation of China(Grant No.51605212)the Natural Science Foundation of Gansu Province(Grant No.17JR5RA122)the Project of Hongliu First-class Disciplines Development Program of Lanzhou University of Technology.
文摘The fatigue damage model based on theory of damage mechanics is capable of predicting the fatigue life under multiaxial loading.Meanwhile,the application of critical plane method in the prediction of multiaxial fatigue life has made certain progress.According to the law of thermodynamics,a new damage evolution equation is developed in the present study to predict the fatigue life of geometrically discontinuous structure under tension-torsion loading based on damage mechanics and the critical plane method.The essence of this approach is tha t the st rain parame ter of the uniaxial nonlinear fatigue damage model is replaced with the equivalent strain,which consists of the releva nt parame ters of the critical plane.However,it is difficult to calculate the stress-strain status and the critical plane position of geometrically dis?continuous structure by theoretical methods because of the existence of stress concentration and the multiaxial nonproportional characteristics.Therefore,a new numerical simulation method is proposed to determine the critical plane of geometrically discontinuous structure under multiaxial loading by means of the finite element method and MATLAB software.The fatigue life of notched specimens subjected to combined bending and torsion is predicted using the proposed met hod,and the result is compared with t hose from the experimen ts and the Manson-Cfiffin law.The comparisons show that the proposed method is superior to the Manson-Coffin law and is capable of reproducing the experimental results reasonably when the geometry of the structure is complex.It completely meets the needs of engineering practice.
文摘The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first_order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C_shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.
文摘This paper describes geometrical essentials of some iteration methods （e.g. Newton iteration, secant line method, etc.） for solving nonlinear equations and advances some geometrical methods of iteration that are flexible and efficient.
基金This project is supported by Provincial Project Foundation of Science and Technology of Guangdong, China(No.2002104040101).
文摘In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.
基金supported by grants from the National Natural Science Foundation of China (51478130)the Guangzhou Municipal Education Bureau’s Scientific Research Project, China (2024312217)+1 种基金the China Scholarship Council (201808440070)the 111 Project of China (D21021).
文摘This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.
文摘Micro cracking in solid material is a phenomenon of solid plane bifurcation with incompatible micro-rotation. Therefore it is important to analysis the micro-rotation field in a strained material in fractute mechanics. On the basis of strain-rotation (S-R) decomposition theorem, the geometric criterion of cracking has been established. In order to verify its validity, using grating experimental method the thin compact tension specimens were investigated. The analysis results show a good agreement between the experimental and geometric criterion.
文摘This paper introduced a bootstrap method called truncated geometric bootstrap method for time series stationary process. We estimate the parameters of a geometric distribution which has been truncated as a probability model for the bootstrap algorithm. This probability model was used in resampling blocks of random length, where the length of each blocks has a truncated geometric distribution. The method was able to determine the block sizes b and probability p attached to its random selections. The mean and variance were estimated for the truncated geometric distribution and the bootstrap algorithm developed based on the proposed probability model.
文摘In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.
基金supported by the Natural Science Foundation of Hubei Province(CN)(Grant No.2019CFB693)the Research Foundation of the Education Department of Hubei Province(CN)(Grant No.B2019003)the open Foundation of the Key Laboratory of Metallurgical Equipment and Control of Education Ministry(CN)(Grant No.2015B14).
文摘The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.
文摘We present in this paper a numerical study of the validity limit of the geometrical optics approximation compared with a differential method which is established according to rigorous formalisms based on the electromagnetic theory. The precedent studies show that this method is adapted to the study of diffraction by periodic rough surfaces. We determine by two methods the emissivity of gold and tungsten for surfaces with a rectangular or sinusoidal profile, for a wavelength equal to 0.55 microns. The monochromatic directional emissivity of these surfaces clearly depends on the angle of incidence, the surface profile, height, period and the nature of the material. We perform our calculations by a method of coupled wave analysis (CWA) and a geometric optics method (GOA). The latter method is theoretically valid only when the dimensions of the cavities are very large compared to the wavelength, while the CWA is theoretically correct whatever these dimensions. The main purpose of this work is to investigate the validity limit of GOA compared with CWA. The obtained results for a fixed height of the grating, allowed us to delimit the validity domain of the optic geometrical approximation for the treated cases. Finally, the agreement between the emissivity calculated by the differential method and that given on the basis of the homogenization theory is satisfactory when the period is much smaller than the wavelength.
文摘This paper seeks to develop an efficient multigrid algorithm for solving the Burgers problem with the use of non-orthogonal structured curvilinear grids in L-shaped geometry.For this,the differential equations were discretized by Finite Volume Method(FVM)with second-order approximation scheme and deferred correction.Moreover,the algebraic method and the differential method were used to generate the non-orthogonal structured curvilinear grids.Furthermore,the influence of some parameters of geometric multigrid method,as well as lexicographical Gauss–Seidel(Lex-GS),η-line Gauss–Seidel(η-line-GS),Modified Strongly Implicit(MSI)and modified incomplete LU decomposition(MILU)solvers on the Central Processing Unit(CPU)time was investigated.Therefore,several parameters of multigrid method and solvers were tested for the problem,with the use of nonorthogonal structured curvilinear grids and multigrid method,resulting in an algorithm with the combination that achieved the best results and CPU time.The geometric multigrid method with Full Approximation Scheme(FAS),V-cycle and standard coarsening ratio for this problem were utilized.This article shows how to calculate the coordinates transformation metrics in the coarser grids.Results show that the MSI and MILU solvers are the most efficient.Moreover,theMSI solver is faster thanMILU for both grids generators;and the solutions are more accurate for the Burgers problem with grids generated using elliptic equations.
文摘The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.
文摘A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) not only the stress fields of a circular disk can be displayed, but the u -and v -displacement fields can also be produced.
