This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationa...This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.展开更多
The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimiz...The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimizing their design.However,conventional numerical methods can face challenges with the non-linearities inherent in hyperelasticity and the complex spatial variations in FGMs.This paper presents a novel hybrid numerical approach combining Physics-Informed Neural Networks(PINNs)with Finite Element Method(FEM)derived data for the robust analysis of thick-walled,axisymmetric,heterogeneous,hyperelastic pressure vessels with elliptical geometries.A PINN framework incorporating neo-Hookean constitutive relations is developed in MATLAB.To enhance training efficiency and accuracy,the PINN’s loss function is augmented with displacement data obtained from high-fidelity FEM simulations performed in ANSYS.The methodology is rigorously validated by comparing PINN-predicted displacement and von Mises stress fields against ANSYS benchmarks for various scenarios of FGMconfigurations(with material properties varying according to a power law)subjected to internal and external pressurization.The results demonstrate excellent agreement between the proposed hybrid PINN-FEMapproach and conventional FEMsolutions across all test cases,accurately capturing complex deformation patterns and stress concentrations.This study highlights the potential of data-augmented PINNs as an effective and accurate computational tool for tackling complex solid mechanics problems involving non-linearmaterials and significant heterogeneity,offering a promising avenue for future research in engineering design and analysis.展开更多
Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a l...Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.展开更多
In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Conver...In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.展开更多
One of the practical approaches in identifying structures is the non-linear resonant decay method which identifies a non-linear dynamic system utilizing a model based on linear modal space containing the underlying li...One of the practical approaches in identifying structures is the non-linear resonant decay method which identifies a non-linear dynamic system utilizing a model based on linear modal space containing the underlying linear system and a small number of extra terms that exhibit the non-linear effects.In this paper,the method is illustrated in a simulated system and an experimental structure.The main objective of the non-linear resonant decay method is to identify the non-linear dynamic systems based on the use of a multi-shaker excitation using appropriated excitation which is obtained from the force appropriation approach.The experimental application of the method is indicated to provide suitable estimates of modal parameters for the identification of non-linear models of structures.展开更多
We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to ...We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).展开更多
In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Trans...In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.展开更多
This paper presents a microcomputer program system for the three-dimensional non-linear finite element method (FEM), which is developed for applications in mining engineering, such as excavation, filling, caving, damm...This paper presents a microcomputer program system for the three-dimensional non-linear finite element method (FEM), which is developed for applications in mining engineering, such as excavation, filling, caving, damming, etc. Our attention is focused on some of the key aspects of the programing techniques for an applicable microcomputer program, such as the structure of the program, data processing in the core, communication between the core and mass storage, input and output systems, etc. In the last part of this paper, two test examples are presented verifying the computational accuracy of this programing system and illustrating its applicability in practical engineering.展开更多
Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynam...Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.展开更多
高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(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%,动能快速下降与内能急剧增加是建筑破坏的关键特征。展开更多
In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length wa...In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length was set as the resistance in the limit state. The nonlinear FEM was used to obtain the crack length of the foundation surface of the gravity dam, and the linear response surface method based on the orthogonal test design method was used to calculate the reliability, providing a reasonable and simple method for calculating the reliability of the serviceability limit state. The Longtan RCC gravity dam was chosen as an example. An orthogonal test, including eleven factors and two levels, was conducted, and the tensile reliability was calculated. The analysis shows that this method is reasonable.展开更多
In Europe, computation of displacement demand for seismic assessment of existing buildings is essentially based on a simplified formulation of the N2 method as prescribed by Eurocode 8(EC8). However, a lack of accurac...In Europe, computation of displacement demand for seismic assessment of existing buildings is essentially based on a simplified formulation of the N2 method as prescribed by Eurocode 8(EC8). However, a lack of accuracy of the N2 method in certain conditions has been pointed out by several studies. This paper addresses the assessment of effectiveness of the N2 method in seismic displacement demand determination in non-linear domain. The objective of this work is to investigate the accuracy of the N2 method through comparison with displacement demands computed using non-linear timehistory analysis(NLTHA). Results show that the original N2 method may lead to overestimation or underestimation of displacement demand predictions. This may affect results of mechanical model-based assessment of seismic vulnerability at an urban scale. Hence, the second part of this paper addresses an improvement of the N2 method formula by empirical evaluation of NLTHA results based on EC8 ground-classes. This task is formulated as a mathematical programming problem in which coefficients are obtained by minimizing the overall discrepancy between NLTHA and modified formula results. Various settings of the mathematical programming problem have been solved using a global optimization metaheuristic. An extensive comparison between the original N2 method formulation and optimized formulae highlights benefits of the strategy.展开更多
This paper demonstrates the approximate analytical solution to a non-linear singular two-point boundary-value problem which describes oxygen diffusion in a planar cell. The model is based on diffusion equation contain...This paper demonstrates the approximate analytical solution to a non-linear singular two-point boundary-value problem which describes oxygen diffusion in a planar cell. The model is based on diffusion equation containing a non-linear term related to Michaelis-Menten kinetics of enzymatic reaction. Approximate analytical expression of concentration of oxygen is derived using new Homotopy perturbation method for various boundary conditions. The validity of the obtained solutions is verified by the numerical results.展开更多
For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describ...For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describes, an adaptation of Newton-Raphson method was used for solving the highly nonlinear system of equations describing the formation of equilibrium products in reacting of fuel-additive-air mixtures. This study also shows what possible of the results. In this paper, to be present the efficient numerical algorithms for. solving the combustion problem, to be used nonlinear equations based on the iteration method and high order of the Taylor series. The modified Adomian decomposition method was applied to construct the numerical algorithms. Some numerical illustrations are given to show the efficiency of algorithms. Comparisons of results by the new Matlab routines and previous routines, the result data indicate that the new Matlab routines are reliable, typical deviations from previous results are less than 0.05%.展开更多
An iterative algorithm for modeling of non-linear joint by the displacement discontinuity method (DDM) was described, and the effect of the non-linear joint on the in-situ stress field was investigated in this paper. ...An iterative algorithm for modeling of non-linear joint by the displacement discontinuity method (DDM) was described, and the effect of the non-linear joint on the in-situ stress field was investigated in this paper. The Barton-Bandis (BB) non-linear joint model and failure criterion were adopted in the new DDM program. Using this program, the stress field around the non-linear joint was obtained, the parameters analysis of the joint was carried out, and the deformation and stress distribution of the joint were studied. The simulation results show that: (1)the in-situ stress is significantly affected by the joint; (2)the increase of stiffness, friction angle, and thickness of the joint affect the stress concentration in different ways; (3)the influence distance of the joint changes with the angle of the joint; (4)the deformation and stress of the joint change with the point position.展开更多
文摘This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.
文摘The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimizing their design.However,conventional numerical methods can face challenges with the non-linearities inherent in hyperelasticity and the complex spatial variations in FGMs.This paper presents a novel hybrid numerical approach combining Physics-Informed Neural Networks(PINNs)with Finite Element Method(FEM)derived data for the robust analysis of thick-walled,axisymmetric,heterogeneous,hyperelastic pressure vessels with elliptical geometries.A PINN framework incorporating neo-Hookean constitutive relations is developed in MATLAB.To enhance training efficiency and accuracy,the PINN’s loss function is augmented with displacement data obtained from high-fidelity FEM simulations performed in ANSYS.The methodology is rigorously validated by comparing PINN-predicted displacement and von Mises stress fields against ANSYS benchmarks for various scenarios of FGMconfigurations(with material properties varying according to a power law)subjected to internal and external pressurization.The results demonstrate excellent agreement between the proposed hybrid PINN-FEMapproach and conventional FEMsolutions across all test cases,accurately capturing complex deformation patterns and stress concentrations.This study highlights the potential of data-augmented PINNs as an effective and accurate computational tool for tackling complex solid mechanics problems involving non-linearmaterials and significant heterogeneity,offering a promising avenue for future research in engineering design and analysis.
文摘泥石流作为高破坏性混合流体,其携带的块石对框架结构的冲击机制尚未得到充分研究。为揭示含块石泥石流冲击框架结构的动力响应与损伤机制,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合数值方法,构建流体-块石-结构的多耦合数值模型,模拟不同冲击速度和角度下框架结构的损伤过程,并结合两相溃坝试验来验证模型有效性。研究表明:泥石流冲击导致结构损伤经历“接触-扩散-反弹-堆积/冲击”四个阶段,当冲击速度超过6.00 m/s时,结构产生不可恢复损伤,且块石阻隔作用使流体上部冲击力大于下部,引发结构中部最先发生集中损伤;此外,冲击力峰值随速度和角度增大呈非线性增长,10.00 m/s与90°工况下框架柱底冲击力达497.17 kN,超过结构抗冲击承载力;最后,数值模拟与经验公式所得冲击力结果误差在13.95%~29.00%,数量级一致,验证了模型可靠性。研究结果可为泥石流高发区框架结构的抗冲击设计提供参考。
文摘Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.
基金the Natural Science Foundation of China(Grant Nos.61673169,11301127,11701176,11626101,and 11601485)The Natural Science Foundation of Huzhou City(Grant No.2018YZ07).
文摘In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.
文摘One of the practical approaches in identifying structures is the non-linear resonant decay method which identifies a non-linear dynamic system utilizing a model based on linear modal space containing the underlying linear system and a small number of extra terms that exhibit the non-linear effects.In this paper,the method is illustrated in a simulated system and an experimental structure.The main objective of the non-linear resonant decay method is to identify the non-linear dynamic systems based on the use of a multi-shaker excitation using appropriated excitation which is obtained from the force appropriation approach.The experimental application of the method is indicated to provide suitable estimates of modal parameters for the identification of non-linear models of structures.
基金supported by the National Natural Science Foundation of China (41721003, 41974022, 41774024, 41874001)Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, China(20-02-05)
文摘We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).
文摘In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.
文摘This paper presents a microcomputer program system for the three-dimensional non-linear finite element method (FEM), which is developed for applications in mining engineering, such as excavation, filling, caving, damming, etc. Our attention is focused on some of the key aspects of the programing techniques for an applicable microcomputer program, such as the structure of the program, data processing in the core, communication between the core and mass storage, input and output systems, etc. In the last part of this paper, two test examples are presented verifying the computational accuracy of this programing system and illustrating its applicability in practical engineering.
基金supported by the National Natural Science Foundation of China(51879159,51490675,11432009,and 51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for doctoral student
文摘Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.
文摘高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。
文摘In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length was set as the resistance in the limit state. The nonlinear FEM was used to obtain the crack length of the foundation surface of the gravity dam, and the linear response surface method based on the orthogonal test design method was used to calculate the reliability, providing a reasonable and simple method for calculating the reliability of the serviceability limit state. The Longtan RCC gravity dam was chosen as an example. An orthogonal test, including eleven factors and two levels, was conducted, and the tensile reliability was calculated. The analysis shows that this method is reasonable.
文摘In Europe, computation of displacement demand for seismic assessment of existing buildings is essentially based on a simplified formulation of the N2 method as prescribed by Eurocode 8(EC8). However, a lack of accuracy of the N2 method in certain conditions has been pointed out by several studies. This paper addresses the assessment of effectiveness of the N2 method in seismic displacement demand determination in non-linear domain. The objective of this work is to investigate the accuracy of the N2 method through comparison with displacement demands computed using non-linear timehistory analysis(NLTHA). Results show that the original N2 method may lead to overestimation or underestimation of displacement demand predictions. This may affect results of mechanical model-based assessment of seismic vulnerability at an urban scale. Hence, the second part of this paper addresses an improvement of the N2 method formula by empirical evaluation of NLTHA results based on EC8 ground-classes. This task is formulated as a mathematical programming problem in which coefficients are obtained by minimizing the overall discrepancy between NLTHA and modified formula results. Various settings of the mathematical programming problem have been solved using a global optimization metaheuristic. An extensive comparison between the original N2 method formulation and optimized formulae highlights benefits of the strategy.
文摘This paper demonstrates the approximate analytical solution to a non-linear singular two-point boundary-value problem which describes oxygen diffusion in a planar cell. The model is based on diffusion equation containing a non-linear term related to Michaelis-Menten kinetics of enzymatic reaction. Approximate analytical expression of concentration of oxygen is derived using new Homotopy perturbation method for various boundary conditions. The validity of the obtained solutions is verified by the numerical results.
文摘For a specific combustion problem involving calculations of several species at the equilibrium state, it is simpler to write a general computer program and calculate the combustion concentration. Original work describes, an adaptation of Newton-Raphson method was used for solving the highly nonlinear system of equations describing the formation of equilibrium products in reacting of fuel-additive-air mixtures. This study also shows what possible of the results. In this paper, to be present the efficient numerical algorithms for. solving the combustion problem, to be used nonlinear equations based on the iteration method and high order of the Taylor series. The modified Adomian decomposition method was applied to construct the numerical algorithms. Some numerical illustrations are given to show the efficiency of algorithms. Comparisons of results by the new Matlab routines and previous routines, the result data indicate that the new Matlab routines are reliable, typical deviations from previous results are less than 0.05%.
基金Western Transport Construction Science and Technology Project of the Ministry of Transport of China ( No. 2009318000046)
文摘An iterative algorithm for modeling of non-linear joint by the displacement discontinuity method (DDM) was described, and the effect of the non-linear joint on the in-situ stress field was investigated in this paper. The Barton-Bandis (BB) non-linear joint model and failure criterion were adopted in the new DDM program. Using this program, the stress field around the non-linear joint was obtained, the parameters analysis of the joint was carried out, and the deformation and stress distribution of the joint were studied. The simulation results show that: (1)the in-situ stress is significantly affected by the joint; (2)the increase of stiffness, friction angle, and thickness of the joint affect the stress concentration in different ways; (3)the influence distance of the joint changes with the angle of the joint; (4)the deformation and stress of the joint change with the point position.
