The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes libra...The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes library of finite element analysis program and Windows-type Graphical User Interfaces by VC + + and its MFC are developed. The reliability, reusability and extensibility of program are enhanced. It is a reference to develop the large-scale, versatile and powerful systems of object-oriented finite element software.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess v...The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.展开更多
The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric va...The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.展开更多
We present the formulation and applications of JAX-CPFEM,an open-source,GPU-accelerated,and differentiable 3-D crystal plasticity finite element method(CPFEM)software package.Leveraging the modern computing architectu...We present the formulation and applications of JAX-CPFEM,an open-source,GPU-accelerated,and differentiable 3-D crystal plasticity finite element method(CPFEM)software package.Leveraging the modern computing architecture JAX,JAX-CPFEM features high performance through array programming and GPU acceleration,achieving a 39×speedup in a polycrystal case with~52,000 degrees of freedom compared to MOOSE with MPI(8 cores).Furthermore,JAX-CPFEM utilizes the automatic differentiation technique,enabling users to handle complex,non-linear constitutive materials laws without manually deriving the case-specific Jacobian matrix.Beyond solving forward problems,JAX-CPFEM demonstrates its potential in an inverse design pipeline,where initial crystallographic orientations of polycrystal copper are optimized to achieve targeted mechanical properties under deformations.The end-to-end differentiability of JAX-CPFEM allows automatic sensitivity calculations and high-dimensional inverse design using gradient-based optimization.The concept of differentiable JAX-CPFEM provides an affordable,flexible,and multi-purpose tool,advancing efficient and accessible computational tools for inverse design in smart manufacturing.展开更多
The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iter...The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iteration algorithm is used to solve this formulation, which distinguishes the integration points of the rigid zones and the plastic zones and solves a series of the quadratic programming to overcome the difficulties caused by the nonsmoothness and the nonlinearity of the objective function. This method has been used to carry out the rigid-plastic FEM analysis. An example is given to demonstrate the effectiveness of this method.展开更多
This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound ...This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.展开更多
This paper shows that the alternating direction method can be used to solve the structured inverse quadratic eigenvalue problem with symmetry, positive semi-definiteness and sparsity requirements. The results of numer...This paper shows that the alternating direction method can be used to solve the structured inverse quadratic eigenvalue problem with symmetry, positive semi-definiteness and sparsity requirements. The results of numerical examples show that the proposed method works well.展开更多
Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software desig...Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software designed using an object-oriented approach can be significantly more robust than traditional codes. This paper describes a special kind of implementation of object-oriented programming which is rather hybrid in nature, in the development of a finite element code for engineering analysis of metal working problems using C++, and discusses the advantages of this approach.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class ...This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class hierarchy structure of object-oriented finite element structural analysis, and gives a part C+ + code description.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for th...The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.展开更多
A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The ...A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.展开更多
In the present study,the performance of reinforced concrete tunnel(RCT)under internal water pressure is evaluated by using nonlinear finite element analysis and surrogate models.Several parameters,including the compre...In the present study,the performance of reinforced concrete tunnel(RCT)under internal water pressure is evaluated by using nonlinear finite element analysis and surrogate models.Several parameters,including the compressive and tensile strength of concrete,the size of the longitudinal reinforcement bar,the transverse bar diameter,and the internal water pre ssure,are considered as the input variables.Based on the levels of variables,36 mix designs are selected by the Taguchi method,and 12 mix designs are proposed in this study.Carbon fiber reinforced concrete(CFRC)or glass fiber reinforced concrete(GFRC)is considered for simulating these 12 samples.Principal component regression(PCR),Multi Ln equation regression(MLnER),and gene expression programming(GEP)are employed for predicting the percentage of damaged surfaces(PDS)of the RCT,the effective tensile plastic strain(ETPS),the maximum deflection of the RCT,and the deflection of crown of RCT.The error terms and statistical parameters,including the maximum positive and negative errors,mean absolute percentage error(MAPE),root mean square error(RMSE),coefficient of determination,and normalized square error(NMSE),are utili zed to evaluate the accuracy of the models.Based on the results,GEP performs better than other models in predicting the outputs.The results sh ow that the internal water pressure and the mechanical properties of concrete have the most effect on the damag e and deflection of the RCT.展开更多
A mathematical model of optimization for cargo hold structure of small and medium-sized duplex stainless steel (DSS) chemical tankers based on the rule is described.The weight of DSS is defined as the objective functi...A mathematical model of optimization for cargo hold structure of small and medium-sized duplex stainless steel (DSS) chemical tankers based on the rule is described.The weight of DSS is defined as the objective function.The explicit formulas of the relevant rules,the explicit geometric relations formed by the structural arrangement and the implicit constraints such as layout,lines,expert experience,workmanship and so on are taken as the constraints.The design variables and the algorithm are determined by analyzing the logical relations among the factors.Besides,the efficiency of the optimum algorithm is improved by utilizing the distribution law of the optimal corrugations.Finally,the comprehensive optimization of the cargo hold DSS is accomplished by C++ programming.The weight of the optimized scheme of the example is 4.93% lower than that of the original scheme.The finite element analysis results show that the optimized scheme is satisfied with the requirements of the rules and has a good structure performance.The optimization method based on the rules and accomplished by the programming is practical for structural preliminary design.展开更多
文摘The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes library of finite element analysis program and Windows-type Graphical User Interfaces by VC + + and its MFC are developed. The reliability, reusability and extensibility of program are enhanced. It is a reference to develop the large-scale, versatile and powerful systems of object-oriented finite element software.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金part of the TPS projecta Vied-Newton PhD scholarship+1 种基金a Dixon scholarship from Imperial College London,UKthe Dean’s Fund from Imperial College London for financial support(2017-2020)。
文摘The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.
基金Project supported by the National Natural Science Foundation of China(Nos.10225212, 10421002 and 10332010)the NCET Program provided by the Ministry of Education and the National Key Basic Research Special Foundation of China (No.2005CB321704)
文摘The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.
基金support from the Department of Defense Vannevar Bush Faculty Fellowship,USA N00014-19-1-2642from the NSF Engineering Research Center for Hybrid Autonomous Manufacturing Moving from Evolution to Revolution(ERC‐HAMMER)under Award Number EEC-2133630.
文摘We present the formulation and applications of JAX-CPFEM,an open-source,GPU-accelerated,and differentiable 3-D crystal plasticity finite element method(CPFEM)software package.Leveraging the modern computing architecture JAX,JAX-CPFEM features high performance through array programming and GPU acceleration,achieving a 39×speedup in a polycrystal case with~52,000 degrees of freedom compared to MOOSE with MPI(8 cores).Furthermore,JAX-CPFEM utilizes the automatic differentiation technique,enabling users to handle complex,non-linear constitutive materials laws without manually deriving the case-specific Jacobian matrix.Beyond solving forward problems,JAX-CPFEM demonstrates its potential in an inverse design pipeline,where initial crystallographic orientations of polycrystal copper are optimized to achieve targeted mechanical properties under deformations.The end-to-end differentiability of JAX-CPFEM allows automatic sensitivity calculations and high-dimensional inverse design using gradient-based optimization.The concept of differentiable JAX-CPFEM provides an affordable,flexible,and multi-purpose tool,advancing efficient and accessible computational tools for inverse design in smart manufacturing.
文摘The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iteration algorithm is used to solve this formulation, which distinguishes the integration points of the rigid zones and the plastic zones and solves a series of the quadratic programming to overcome the difficulties caused by the nonsmoothness and the nonlinearity of the objective function. This method has been used to carry out the rigid-plastic FEM analysis. An example is given to demonstrate the effectiveness of this method.
基金The project supported by National Natural Science Foundation of China.
文摘This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.
基金Supported by Youth Teacher Education and Research Funds of Fujian(Grant No.JAT170911).
文摘This paper shows that the alternating direction method can be used to solve the structured inverse quadratic eigenvalue problem with symmetry, positive semi-definiteness and sparsity requirements. The results of numerical examples show that the proposed method works well.
文摘Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software designed using an object-oriented approach can be significantly more robust than traditional codes. This paper describes a special kind of implementation of object-oriented programming which is rather hybrid in nature, in the development of a finite element code for engineering analysis of metal working problems using C++, and discusses the advantages of this approach.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
文摘This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class hierarchy structure of object-oriented finite element structural analysis, and gives a part C+ + code description.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
基金theNationalKeyBasicResearchSpecialFoundation (G1 9990 3 2 80 5 ) the FoundationforUniversityKeyTeacherbytheMinistryofEducationo
文摘The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.
基金This is part of the TPS projecta Vied-Newton PhD scholarship and a Dixon scholarship from Imperial College London, UK, for supporting his studies at Imperial College Londonthe Dean’s Fund from Imperial College London for financial support (2017-2020).
文摘A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.
文摘In the present study,the performance of reinforced concrete tunnel(RCT)under internal water pressure is evaluated by using nonlinear finite element analysis and surrogate models.Several parameters,including the compressive and tensile strength of concrete,the size of the longitudinal reinforcement bar,the transverse bar diameter,and the internal water pre ssure,are considered as the input variables.Based on the levels of variables,36 mix designs are selected by the Taguchi method,and 12 mix designs are proposed in this study.Carbon fiber reinforced concrete(CFRC)or glass fiber reinforced concrete(GFRC)is considered for simulating these 12 samples.Principal component regression(PCR),Multi Ln equation regression(MLnER),and gene expression programming(GEP)are employed for predicting the percentage of damaged surfaces(PDS)of the RCT,the effective tensile plastic strain(ETPS),the maximum deflection of the RCT,and the deflection of crown of RCT.The error terms and statistical parameters,including the maximum positive and negative errors,mean absolute percentage error(MAPE),root mean square error(RMSE),coefficient of determination,and normalized square error(NMSE),are utili zed to evaluate the accuracy of the models.Based on the results,GEP performs better than other models in predicting the outputs.The results sh ow that the internal water pressure and the mechanical properties of concrete have the most effect on the damag e and deflection of the RCT.
文摘A mathematical model of optimization for cargo hold structure of small and medium-sized duplex stainless steel (DSS) chemical tankers based on the rule is described.The weight of DSS is defined as the objective function.The explicit formulas of the relevant rules,the explicit geometric relations formed by the structural arrangement and the implicit constraints such as layout,lines,expert experience,workmanship and so on are taken as the constraints.The design variables and the algorithm are determined by analyzing the logical relations among the factors.Besides,the efficiency of the optimum algorithm is improved by utilizing the distribution law of the optimal corrugations.Finally,the comprehensive optimization of the cargo hold DSS is accomplished by C++ programming.The weight of the optimized scheme of the example is 4.93% lower than that of the original scheme.The finite element analysis results show that the optimized scheme is satisfied with the requirements of the rules and has a good structure performance.The optimization method based on the rules and accomplished by the programming is practical for structural preliminary design.
