At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The who...At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.展开更多
The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at dif...The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing...Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing the number of design variables. In this study, to achieve lighter weight and higher stiffness, the simplified model of BIP was developed and combined with an optimization procedure;consequently, optimal designs of automotive body B-pillar were produced. B-pillar was divided into four quarters and each quarter was modelled by one simplified beam. In the optimization procedure, depth, width, and thickness of the simplified beams were considered as the design variables.Weight, bending and torsional stiffness were also considered as objective functions. The optimization procedure is composed of six stages: designing the experiments, calculating grey relational grade, calculating signal-to noise ratio,finding an optimum design using Taguchi grey relational analysis, performing sensitivity analysis using analysis of variance(ANOVA) and performing non-dominated sorting and multi-criteria decision making. The results show that the width of lower B-pillar has the highest effect(about 55%) and the obtained optimum design point could reduce the weight of B-pillar by about 40% without reducing the BIP stiffness by more than 1.47%.展开更多
Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction a...Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction and wear during exercise.The vital mechanical function展开更多
Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the sp...Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the spinal cord,which could be result in life threatening展开更多
On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary resul...On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary results show that many strong earthquakes occurred in the area of increased stress in the model. Though the results are preliminary, the quasi 3D finite element model is meaningful for strong earthquake prediction.展开更多
Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each oth...Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each other. The billet suffers plastic deformation (essentially, localized compressio n) during its rotation between the rotating dies. Compared to other numerical si mulation methods, the finite element method (FEM) has advantages in solving gene ral problems with complex shapes of the formed parts. In cross-wedge rolling, t here are four stages in the workpiece deformation process, namely knifing, guidi ng, stretching and sizing stage. It is time-consuming and expensive to design t he CWR process by trial and error method. The application of numerical simul ation for the CWR process will help engineers to efficiently improve the process development. Tselikov, Hayama, Jain and Kobayashi, and Higashimo applied the sl ip-line theory in study of CWR process analysis. Zb.pater studied CWR process i ncluding upsetting by upper-bound method. The above numerical simulation were b ased on the two-dimensional plain-strain assumption ignored the metal flow in workpiece axial direction. Therefore, the complex three-dimensional stress and deformation involved in CWR processes were not presented. Compared to other nume rical simulation methods, the finite element method (FEM) has advantages in solv ing general problems with complex shapes of the formed parts. As yet, a few 3-D finite element simulation studies on CWR process have been reported in literatu res. In this paper, the process of cross wedge rolling (CWR) has been simulated and analyzed by 3D rigid-plastic finite element method. Considering the charact eristic of CWR, the static implicit FEM program is selected. The models proposed in this study uses the commercial code DEFORM 3D to simulate the CWR process. T his is an implicit Lagrangian finite element code, which includes many new enhan cements functions. A new method of utilizing multiple processors using the MPI s tandard has been implemented. Automatic switching between the two different defo rmation solvers (Sparse Solver and Conjugate Gradient Solver) has also been impl emented in order to increase the speed of simulations. In this paper, all stages in CWR process are simulated to be able to closely understand and analyze the a ctual CWR process. For simulating all forming stages in CWR process, the dynam ic adaptive remeshing technology for tetrahedral solid elements was applied. T he stress distributions in cross section of forming workpiece are analyzed to in terpret fracture or rarefaction in the center of workpiece. Authors also analyze d the time-torque curve and the laws of load changing.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
In the process of thin-wall parts assembly for an antenna,the parts assembly deformation deviation is occurring due to the riveting assembly.In view of the riveting assembly deformation problems,it can be analyzed thr...In the process of thin-wall parts assembly for an antenna,the parts assembly deformation deviation is occurring due to the riveting assembly.In view of the riveting assembly deformation problems,it can be analyzed through transient and static simulation.In this work,the theoretical deformation model for riveting assembly is established with round head rivet.The simulation analysis for riveting deformation is carried out with the riveting assembly piece including four rivets,which comparing with the measuring points experiment results of riveting test piece through dealing with the experimental data using the point coordinate transform method and the space line fitting method.Simultaneously,the deformation deviation of the overall thin-wall parts assembly structure is analyzed through finite element simulation;and its results are verified by the measuring experiment for riveting assembly with the deformation deviation of some key points on the thin-wall parts.Through the comparison analysis,it is shown that the simulation results agree well with the experimental results,which proves the correctness and effectiveness of the theoretical analysis,simulation results and the given experiment data processing method.Through the study on the riveting assembly for thin-wall parts,it will provide a theoretical foundation for improving thin-wall parts assembly quality of large antenna in future.展开更多
The aim of this work is to analyze the stress distributions on a crown-luting cement-substrate system with a finite-element method in order to predict the likelihood of interfacial micro cracks, radial or circumferent...The aim of this work is to analyze the stress distributions on a crown-luting cement-substrate system with a finite-element method in order to predict the likelihood of interfacial micro cracks, radial or circumferential cracks, delamination, fracture and delamination with torsion. The contact and layer interface stresses in elastic layered half-space indented by an elastic sphere were examined using finite element method. The model consists of crown, luting cement and substrate. The solutions were carried out for three different elastic moduli of luting cement. It was placed between the cement and the substrate as a middle layer and its elastic module was chosen lower than the elastic module of crown and higher than the elastic module of dentin. An axisymmetric finite element mesh was set up for the stress analysis. Stress distributions on the contact surface and the interfaces of crown-luting cement and luting cement-dentin have been investigated for three different values of luting cement by using ANSYS. The effects of the luting cement which has three different elastic moduli on the pressure distribution and the location of interfacial stresses of the multi-layer model have been examined. The mechanism of crack initiation in the interfaces and interracial delamination was also studied quantitatively. For each luting cement, the pressure distribution is similar at the contact zone. Stress discontinuities occur at the perfect bonding interfaces of the crown-luting cement and the substrate-luting cement. The maximum stress jumps are obtained for the highest and the lowest elastic module of the luting cement. In the crown-luting cement-substrate system, failures may initiate at crown-luting cement region for luting cement with the lowest elastic module value. In addition, failures at luting cement-substrate region may occur for luting cement with the highest elastic module. In the luting cement, the medium elastic module value is more suitable for stress distribution in crown-luting cement-substrate interfaces.展开更多
高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框...高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。展开更多
Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studi...Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.展开更多
剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部...剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部宏-微观损伤模型的剪切型裂纹动态开裂模拟方法,定义了基于偏应变概念的物质点对的正伸长量,可作为预测剪切型裂纹扩展行为的动态开裂准则,一点的损伤定义为该点影响域范围内连接的物质键损伤的加权平均值,而物质键的损伤则与基于偏应变概念的物质点对的正伸长量相关联,并引入能量退化函数建立结构域几何拓扑损伤与能量损失之间的关系,将拓扑损伤与应力应变联系起来,通过能量退化函数修正了SBFEM的刚度系数矩阵,得到了子域在损伤状态下的刚度矩阵,推导了考虑结构损伤的SBFEM动力控制方程,采用Newmark隐式算法对控制方程进行时间离散.最后,通过3个典型算例验证了建议的模型可较好地模拟剪切型断裂问题,能够很好地捕捉剪切型裂纹的扩展路径,并得到较为准确的载荷-位移曲线.展开更多
文摘At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.
文摘The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.
文摘Optimization of an automotive body structure faces the difficulty of having too many design variables and a too large design search space. A simplified model of body-in-prime(BIP) can solve this difficulty by reducing the number of design variables. In this study, to achieve lighter weight and higher stiffness, the simplified model of BIP was developed and combined with an optimization procedure;consequently, optimal designs of automotive body B-pillar were produced. B-pillar was divided into four quarters and each quarter was modelled by one simplified beam. In the optimization procedure, depth, width, and thickness of the simplified beams were considered as the design variables.Weight, bending and torsional stiffness were also considered as objective functions. The optimization procedure is composed of six stages: designing the experiments, calculating grey relational grade, calculating signal-to noise ratio,finding an optimum design using Taguchi grey relational analysis, performing sensitivity analysis using analysis of variance(ANOVA) and performing non-dominated sorting and multi-criteria decision making. The results show that the width of lower B-pillar has the highest effect(about 55%) and the obtained optimum design point could reduce the weight of B-pillar by about 40% without reducing the BIP stiffness by more than 1.47%.
基金National Natural Science Foundation of China,10872147Natural Science Foundation of Tianjin,09JCYBJC1400
文摘Articular cartilage is a layer of low-friction,load-bearing soft hydrated tissue covering bone-ends in diarthrosis,which plays an important role in spreading the load,reducing the joint contact stress,joint friction and wear during exercise.The vital mechanical function
文摘Introduction-The cervical spine is subjected to injury frequently,especially among pilots who are usually on the condition of high acceleration.Injuries of the cervical spine will be potential risk of damage to the spinal cord,which could be result in life threatening
文摘On ths basis of interaction between faults, a finite element model for Southwest China is constructed,and the stress adjustment due to the strong earthquake occurrence in this region was studied. The preliminary results show that many strong earthquakes occurred in the area of increased stress in the model. Though the results are preliminary, the quasi 3D finite element model is meaningful for strong earthquake prediction.
文摘Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each other. The billet suffers plastic deformation (essentially, localized compressio n) during its rotation between the rotating dies. Compared to other numerical si mulation methods, the finite element method (FEM) has advantages in solving gene ral problems with complex shapes of the formed parts. In cross-wedge rolling, t here are four stages in the workpiece deformation process, namely knifing, guidi ng, stretching and sizing stage. It is time-consuming and expensive to design t he CWR process by trial and error method. The application of numerical simul ation for the CWR process will help engineers to efficiently improve the process development. Tselikov, Hayama, Jain and Kobayashi, and Higashimo applied the sl ip-line theory in study of CWR process analysis. Zb.pater studied CWR process i ncluding upsetting by upper-bound method. The above numerical simulation were b ased on the two-dimensional plain-strain assumption ignored the metal flow in workpiece axial direction. Therefore, the complex three-dimensional stress and deformation involved in CWR processes were not presented. Compared to other nume rical simulation methods, the finite element method (FEM) has advantages in solv ing general problems with complex shapes of the formed parts. As yet, a few 3-D finite element simulation studies on CWR process have been reported in literatu res. In this paper, the process of cross wedge rolling (CWR) has been simulated and analyzed by 3D rigid-plastic finite element method. Considering the charact eristic of CWR, the static implicit FEM program is selected. The models proposed in this study uses the commercial code DEFORM 3D to simulate the CWR process. T his is an implicit Lagrangian finite element code, which includes many new enhan cements functions. A new method of utilizing multiple processors using the MPI s tandard has been implemented. Automatic switching between the two different defo rmation solvers (Sparse Solver and Conjugate Gradient Solver) has also been impl emented in order to increase the speed of simulations. In this paper, all stages in CWR process are simulated to be able to closely understand and analyze the a ctual CWR process. For simulating all forming stages in CWR process, the dynam ic adaptive remeshing technology for tetrahedral solid elements was applied. T he stress distributions in cross section of forming workpiece are analyzed to in terpret fracture or rarefaction in the center of workpiece. Authors also analyze d the time-torque curve and the laws of load changing.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
基金Project(51675100)supported by the National Natural Science Foundation of ChinaProject(2016ZX04004008)supported by the National Numerical Control Equipment Major Project of ChinaProject(6902002116)supported by the Foundation of Certain Ministry of China
文摘In the process of thin-wall parts assembly for an antenna,the parts assembly deformation deviation is occurring due to the riveting assembly.In view of the riveting assembly deformation problems,it can be analyzed through transient and static simulation.In this work,the theoretical deformation model for riveting assembly is established with round head rivet.The simulation analysis for riveting deformation is carried out with the riveting assembly piece including four rivets,which comparing with the measuring points experiment results of riveting test piece through dealing with the experimental data using the point coordinate transform method and the space line fitting method.Simultaneously,the deformation deviation of the overall thin-wall parts assembly structure is analyzed through finite element simulation;and its results are verified by the measuring experiment for riveting assembly with the deformation deviation of some key points on the thin-wall parts.Through the comparison analysis,it is shown that the simulation results agree well with the experimental results,which proves the correctness and effectiveness of the theoretical analysis,simulation results and the given experiment data processing method.Through the study on the riveting assembly for thin-wall parts,it will provide a theoretical foundation for improving thin-wall parts assembly quality of large antenna in future.
文摘The aim of this work is to analyze the stress distributions on a crown-luting cement-substrate system with a finite-element method in order to predict the likelihood of interfacial micro cracks, radial or circumferential cracks, delamination, fracture and delamination with torsion. The contact and layer interface stresses in elastic layered half-space indented by an elastic sphere were examined using finite element method. The model consists of crown, luting cement and substrate. The solutions were carried out for three different elastic moduli of luting cement. It was placed between the cement and the substrate as a middle layer and its elastic module was chosen lower than the elastic module of crown and higher than the elastic module of dentin. An axisymmetric finite element mesh was set up for the stress analysis. Stress distributions on the contact surface and the interfaces of crown-luting cement and luting cement-dentin have been investigated for three different values of luting cement by using ANSYS. The effects of the luting cement which has three different elastic moduli on the pressure distribution and the location of interfacial stresses of the multi-layer model have been examined. The mechanism of crack initiation in the interfaces and interracial delamination was also studied quantitatively. For each luting cement, the pressure distribution is similar at the contact zone. Stress discontinuities occur at the perfect bonding interfaces of the crown-luting cement and the substrate-luting cement. The maximum stress jumps are obtained for the highest and the lowest elastic module of the luting cement. In the crown-luting cement-substrate system, failures may initiate at crown-luting cement region for luting cement with the lowest elastic module value. In addition, failures at luting cement-substrate region may occur for luting cement with the highest elastic module. In the luting cement, the medium elastic module value is more suitable for stress distribution in crown-luting cement-substrate interfaces.
文摘高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。
基金Otokar Otomotiv ve Savunma Sanayi A.S. for the financial support
文摘Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.
文摘剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部宏-微观损伤模型的剪切型裂纹动态开裂模拟方法,定义了基于偏应变概念的物质点对的正伸长量,可作为预测剪切型裂纹扩展行为的动态开裂准则,一点的损伤定义为该点影响域范围内连接的物质键损伤的加权平均值,而物质键的损伤则与基于偏应变概念的物质点对的正伸长量相关联,并引入能量退化函数建立结构域几何拓扑损伤与能量损失之间的关系,将拓扑损伤与应力应变联系起来,通过能量退化函数修正了SBFEM的刚度系数矩阵,得到了子域在损伤状态下的刚度矩阵,推导了考虑结构损伤的SBFEM动力控制方程,采用Newmark隐式算法对控制方程进行时间离散.最后,通过3个典型算例验证了建议的模型可较好地模拟剪切型断裂问题,能够很好地捕捉剪切型裂纹的扩展路径,并得到较为准确的载荷-位移曲线.
