Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite struc...Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.展开更多
Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale pr...Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.展开更多
Dielectric elastomers(DEs)require balanced electric actuation performance and mechanical integrity under applied voltages.Incorporating high dielectric particles as fillers provides extensive design space to optimize ...Dielectric elastomers(DEs)require balanced electric actuation performance and mechanical integrity under applied voltages.Incorporating high dielectric particles as fillers provides extensive design space to optimize concentration,morphology,and distribution for improved actuation performance and material modulus.This study presents an integrated framework combining finite element modeling(FEM)and deep learning to optimize the microstructure of DE composites.FEM first calculates actuation performance and the effective modulus across varied filler combinations,with these data used to train a convolutional neural network(CNN).Integrating the CNN into a multi-objective genetic algorithm generates designs with enhanced actuation performance and material modulus compared to the conventional optimization approach based on FEM approach within the same time.This framework harnesses artificial intelligence to navigate vast design possibilities,enabling optimized microstructures for high-performance DE composites.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoreti...Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.展开更多
The precise microscopic feature of carbon-carbon(C/C)composites is essential for an accurate prediction of their mechanical behavior.After fabrication,actual microscopic feature differs from simple ideal spatial model...The precise microscopic feature of carbon-carbon(C/C)composites is essential for an accurate prediction of their mechanical behavior.After fabrication,actual microscopic feature differs from simple ideal spatial model.Micro-computed-tomography(CT)scan can well describe internal microstructures of composites.Therefore,a reconstructed model is developed based on mirco-CT,by a series of prodcedures including extracting components,generating new binary images and establishing a finite element(FE)model.Compared with the model designed by reconstructed commercial software MIMICS,the presented reconstructed FE model is superior in terms of high mesh quality and controllable mesh quantity.The precision of the model is verified by experiment.展开更多
The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-sc...The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-scale finite element methods for solutions to the electric potential and the displacement for composite material in periodic struc- ture under the coupled piezoelectricity are derived. The coupled two-scale relation of the electric potential and the displacement is set up, and some finite element approximate estimates and numerical examples which show the effectiveness of the method are presented.展开更多
This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforc...This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforced composite actuator under electromechanical loading. The four-variable refined plate theory is a simple and efficient higher-order shear deformation theory, which predicts parabolic variation of transverse shear stresses across the plate thickness and satisfies zero traction conditions on the plate free surfaces. The weak form of governing equations is derived using the principle of minimum potential energy, and a 4-node non-conforming rectangular plate element with 8 degrees of freedom per node is introduced for discretizing the domain. Several benchmark problems are solved by the developed MATLAB code and the obtained results are compared with those from exact and other numerical solutions, showing good agreement.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
A composite structure with frame and skin based on cabin structure in a large space telescope is studied in this paper.The frame is composed of longitudinal and transverse beams with hybrid bonded/bolted joints,and th...A composite structure with frame and skin based on cabin structure in a large space telescope is studied in this paper.The frame is composed of longitudinal and transverse beams with hybrid bonded/bolted joints,and the skin is connected to the frame by bolts.Tensile tests are conducted on the structure by a set of test stand.It is observed that residual deformation occurs in the first test of the structure,which is attributed to the relative sliding between the skin and frame because of bolt-hole clearances.The high tightening torque and the increased number of the skin-frame bolts contribute to the high stiffness of the structure.A finite element model(FEM)of this composite structure is established,and the simulation model is verified by the experimental results.The FEM is contrastively analyzed with different frame joints,and it is found that adhesive joining in the hybrid bonded/bolted joints enhances the stiffness of the structure significantly.Given that adhesive plays a leading role in the stiffness of the hybrid joints,Tie contact in FEM is proposed to simulate bonded or hybrid joints when studying the stiffness performance of undamaged structure.展开更多
In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. ...In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.展开更多
In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single clos...In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.展开更多
The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plastici...The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young's modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young's modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. ine Iormer results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.展开更多
To investigate the application of meta-model for finite element( FE) model updating of structures,the performance of two popular meta-model,i. e.,Kriging model and response surface model( RSM),were compared in detail....To investigate the application of meta-model for finite element( FE) model updating of structures,the performance of two popular meta-model,i. e.,Kriging model and response surface model( RSM),were compared in detail. Firstly,above two kinds of meta-model were introduced briefly. Secondly,some key issues of the application of meta-model to FE model updating of structures were proposed and discussed,and then some advices were presented in order to select a reasonable meta-model for the purpose of updating the FE model of structures. Finally,the procedure of FE model updating based on meta-model was implemented by updating the FE model of a truss bridge model with the measured modal parameters. The results showed that the Kriging model was more proper for FE model updating of complex structures.展开更多
The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems tha...The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems that are beyond the capabilities of classical finite element method(FEM).In PFEM,the computational domain is reconfigured for optimal solution by frequent remeshing and boundary updating.PFEM inherits many concepts,such as a Lagrangian description of continuum,from classic geomechanical FEM.This familiarity with more popular numerical methods facilitates learning and application.This work focuses on G-PFEM,a code specifically developed for the use of PFEM in geotechnical problems.The article has two purposes.The first is to give the reader an overview of the capabilities and main features of the current version of the G-PFEM and the second is to illustrate some of the newer developments of the code.G-PFEM can solve coupled hydro-mechanical static and dynamic problems involving the interaction of solid and/or deformable bodies.Realistic constitutive models for geomaterials are available,including features,such as structure and destructuration,which result in brittle response.The solutions are robust,solidly underpinned by numerical technology including mixedfield formulations,robust and mesh-independent integration of elastoplastic constitutive models and a rigorous and flexible treatment of contact interactions.The novel features presented in this work include the contact domain technique,a natural way to capture contact interactions and impose contact constraints between different continuum bodies,as well as a new simplified formulation for dynamic impact problems.The code performance is showcased by the simulation of several soil-structure interaction problems selected to highlight the novel code features:a rigid footing insertion in soft rock,pipeline insertion and subsequent lateral displacement on over-consolidated clay,screw-pile pull-out and the dynamic impact of a free-falling spherical penetrometer into clay.展开更多
Fiber reinforced polymer (FRP) composite materials having advantages such as higher strength to weight than conventional engineering materials, non-corrosiveness and modularization, which should help engineers to ob...Fiber reinforced polymer (FRP) composite materials having advantages such as higher strength to weight than conventional engineering materials, non-corrosiveness and modularization, which should help engineers to obtain more efficient and cost effective structural materials and systems. Currently, FRP composites are becoming more popular in civil engineering applications. The objectives of this research are to study performance and behavior of light weight multi-cellular FRP composite bridge decks (both module and system levels) under various loading conditions through finite element modeling, and to validate analytical response of FRP composite bridge decks with data from laboratory evaluations. The relative deflection, equivalent flexural rigidity, failure load (mode) and load distribution factors (LDF) based on FE results have been compared with experimental data and discussed in detail. The finite element results showing good correlations with experimental data are presented in this work.展开更多
The nonlinear finite element(FE) analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities(or gradients) to the model parameters are of significant i...The nonlinear finite element(FE) analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities(or gradients) to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method(DDM) has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction(SFSI) models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method(FFD).The relative importance(RI) of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM(e.g.,accurate,efficient,insensitive to numerical noise).Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.展开更多
Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of c...Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of complexity, and thus time, to the stress and deformation analyses of such components, even using numerical approaches such as finite elements. The analysis of composite laminate beams subjected to uniform extension, bending, and/or twisting loads was performed by a novel implementation of the usual finite element method. Due to the symmetric features of the deformations, only a thin slice of the beam to be analysed needs to be modelled. Conventional threedimensional solid finite elements were used for the structural discretization. The accurate deformation relationships were formulated and implemented through the coupling of nodal translational degrees of freedom in the numerical analysis. A sample solution for a rectangular composite laminate beam is presented to show the validity and accuracy of the proposed method.展开更多
A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiat...A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.展开更多
基金This work was supported by National Natural Science Foundation of China under Grant 11672266.
文摘Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.
基金Supported by Science Center for Gas Turbine Project of China (Grant No.P2022-B-IV-014-001)Frontier Leading Technology Basic Research Special Project of Jiangsu Province of China (Grant No.BK20212007)the BIT Research and Innovation Promoting Project of China (Grant No.2022YCXZ019)。
文摘Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFB3707803)the National Natural Science Foundation of China(Grant Nos.12072179 and 11672168)+1 种基金the Key Research Project of Zhejiang Lab(Grant No.2021PE0AC02)Shanghai Engineering Research Center for Inte-grated Circuits and Advanced Display Materials.
文摘Dielectric elastomers(DEs)require balanced electric actuation performance and mechanical integrity under applied voltages.Incorporating high dielectric particles as fillers provides extensive design space to optimize concentration,morphology,and distribution for improved actuation performance and material modulus.This study presents an integrated framework combining finite element modeling(FEM)and deep learning to optimize the microstructure of DE composites.FEM first calculates actuation performance and the effective modulus across varied filler combinations,with these data used to train a convolutional neural network(CNN).Integrating the CNN into a multi-objective genetic algorithm generates designs with enhanced actuation performance and material modulus compared to the conventional optimization approach based on FEM approach within the same time.This framework harnesses artificial intelligence to navigate vast design possibilities,enabling optimized microstructures for high-performance DE composites.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
基金This work is supported by the National Natural Science Foundation of China under the Grant 19772037 and 19902014
文摘Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.
基金supported by the National Natural Science Foundation of China (Nos.11272147,10772078)the Aviation Science Foundation (No.2013ZF52074)+1 种基金the State Key Laboratory of Mechanical Structural Mechanics and Control (No.0214G02)the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)
文摘The precise microscopic feature of carbon-carbon(C/C)composites is essential for an accurate prediction of their mechanical behavior.After fabrication,actual microscopic feature differs from simple ideal spatial model.Micro-computed-tomography(CT)scan can well describe internal microstructures of composites.Therefore,a reconstructed model is developed based on mirco-CT,by a series of prodcedures including extracting components,generating new binary images and establishing a finite element(FE)model.Compared with the model designed by reconstructed commercial software MIMICS,the presented reconstructed FE model is superior in terms of high mesh quality and controllable mesh quantity.The precision of the model is verified by experiment.
基金supported by the National Natural Science Foundation of China(Nos.10801042 and 11171257)the Specialized Research Fund for the Doctoral Program of Higher Education of China(No.20104410120001)
文摘The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-scale finite element methods for solutions to the electric potential and the displacement for composite material in periodic struc- ture under the coupled piezoelectricity are derived. The coupled two-scale relation of the electric potential and the displacement is set up, and some finite element approximate estimates and numerical examples which show the effectiveness of the method are presented.
文摘This research presents a finite element formulation based on four-variable refined plate theory for bending analysis of cross-ply and angle-ply laminated composite plates integrated with a piezoelectric fiber-reinforced composite actuator under electromechanical loading. The four-variable refined plate theory is a simple and efficient higher-order shear deformation theory, which predicts parabolic variation of transverse shear stresses across the plate thickness and satisfies zero traction conditions on the plate free surfaces. The weak form of governing equations is derived using the principle of minimum potential energy, and a 4-node non-conforming rectangular plate element with 8 degrees of freedom per node is introduced for discretizing the domain. Several benchmark problems are solved by the developed MATLAB code and the obtained results are compared with those from exact and other numerical solutions, showing good agreement.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
基金Supported by the National Natural Science Foundation of China(No.51805510)Science and Technology Development Plan Project of Jilin Province(No.20200201294JC)。
文摘A composite structure with frame and skin based on cabin structure in a large space telescope is studied in this paper.The frame is composed of longitudinal and transverse beams with hybrid bonded/bolted joints,and the skin is connected to the frame by bolts.Tensile tests are conducted on the structure by a set of test stand.It is observed that residual deformation occurs in the first test of the structure,which is attributed to the relative sliding between the skin and frame because of bolt-hole clearances.The high tightening torque and the increased number of the skin-frame bolts contribute to the high stiffness of the structure.A finite element model(FEM)of this composite structure is established,and the simulation model is verified by the experimental results.The FEM is contrastively analyzed with different frame joints,and it is found that adhesive joining in the hybrid bonded/bolted joints enhances the stiffness of the structure significantly.Given that adhesive plays a leading role in the stiffness of the hybrid joints,Tie contact in FEM is proposed to simulate bonded or hybrid joints when studying the stiffness performance of undamaged structure.
文摘In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.
基金The project supported by the National Natural Science Foundation of China (19932030)
文摘In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.
基金Project supported by the Research Committee of the Hong Kong Polytechnic University (No.G-U042).
文摘The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young's modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young's modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. ine Iormer results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.
基金Sponsored by the National Key Technology Research and Development Program of China(Grant No.2011BAK02B02)
文摘To investigate the application of meta-model for finite element( FE) model updating of structures,the performance of two popular meta-model,i. e.,Kriging model and response surface model( RSM),were compared in detail. Firstly,above two kinds of meta-model were introduced briefly. Secondly,some key issues of the application of meta-model to FE model updating of structures were proposed and discussed,and then some advices were presented in order to select a reasonable meta-model for the purpose of updating the FE model of structures. Finally,the procedure of FE model updating based on meta-model was implemented by updating the FE model of a truss bridge model with the measured modal parameters. The results showed that the Kriging model was more proper for FE model updating of complex structures.
基金financial support by Severo Ochoa Centre of Excellence (2019-2023) Grant No. CEX2018-000797-Sfunded by MCIN/AEI/10.13039/501100011033+1 种基金research projects BIA2017-84752-RPID2020-119598RB-I00
文摘The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems that are beyond the capabilities of classical finite element method(FEM).In PFEM,the computational domain is reconfigured for optimal solution by frequent remeshing and boundary updating.PFEM inherits many concepts,such as a Lagrangian description of continuum,from classic geomechanical FEM.This familiarity with more popular numerical methods facilitates learning and application.This work focuses on G-PFEM,a code specifically developed for the use of PFEM in geotechnical problems.The article has two purposes.The first is to give the reader an overview of the capabilities and main features of the current version of the G-PFEM and the second is to illustrate some of the newer developments of the code.G-PFEM can solve coupled hydro-mechanical static and dynamic problems involving the interaction of solid and/or deformable bodies.Realistic constitutive models for geomaterials are available,including features,such as structure and destructuration,which result in brittle response.The solutions are robust,solidly underpinned by numerical technology including mixedfield formulations,robust and mesh-independent integration of elastoplastic constitutive models and a rigorous and flexible treatment of contact interactions.The novel features presented in this work include the contact domain technique,a natural way to capture contact interactions and impose contact constraints between different continuum bodies,as well as a new simplified formulation for dynamic impact problems.The code performance is showcased by the simulation of several soil-structure interaction problems selected to highlight the novel code features:a rigid footing insertion in soft rock,pipeline insertion and subsequent lateral displacement on over-consolidated clay,screw-pile pull-out and the dynamic impact of a free-falling spherical penetrometer into clay.
基金Funded by Structural Engineering and Applied Mechanics (STREAM) Research Group(No.ENG-51-2-7-11-022-S),Faculty of Engineering,Prince of Songkla University,Hatyai Songkhla,Thailand
文摘Fiber reinforced polymer (FRP) composite materials having advantages such as higher strength to weight than conventional engineering materials, non-corrosiveness and modularization, which should help engineers to obtain more efficient and cost effective structural materials and systems. Currently, FRP composites are becoming more popular in civil engineering applications. The objectives of this research are to study performance and behavior of light weight multi-cellular FRP composite bridge decks (both module and system levels) under various loading conditions through finite element modeling, and to validate analytical response of FRP composite bridge decks with data from laboratory evaluations. The relative deflection, equivalent flexural rigidity, failure load (mode) and load distribution factors (LDF) based on FE results have been compared with experimental data and discussed in detail. The finite element results showing good correlations with experimental data are presented in this work.
基金National Key Research and Development Program of China under Grant No.2016YFC0701106Natural Sciences and Engineering Research Council of Canada via Discovery under Grant No.NSERC RGPIN-2017-05556 Li
文摘The nonlinear finite element(FE) analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities(or gradients) to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method(DDM) has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction(SFSI) models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method(FFD).The relative importance(RI) of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM(e.g.,accurate,efficient,insensitive to numerical noise).Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.
文摘Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of complexity, and thus time, to the stress and deformation analyses of such components, even using numerical approaches such as finite elements. The analysis of composite laminate beams subjected to uniform extension, bending, and/or twisting loads was performed by a novel implementation of the usual finite element method. Due to the symmetric features of the deformations, only a thin slice of the beam to be analysed needs to be modelled. Conventional threedimensional solid finite elements were used for the structural discretization. The accurate deformation relationships were formulated and implemented through the coupling of nodal translational degrees of freedom in the numerical analysis. A sample solution for a rectangular composite laminate beam is presented to show the validity and accuracy of the proposed method.
基金The project supported by the Research Grant Council of Hong Kong (HKUST 722196E, 6039197E)the National Natural Science Foundation of China(59809003)the Foundation of University Key Teacher by the Chinese Ministry of Education
文摘A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.