Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Multi-scale Simulation Method with Coupled Finite/Discrete Element Model and Its Application

The existing research on continuous structure is usually analyzed with finite element method (FEM) and granular medium with discrete element method (DEM), but there are few researches on the coupling interaction between continuous structure and discrete medium. To the issue of this coupling interaction, a multi-scale simulation method with coupled finite/discrete element model is put forward, in their respective domains of discrete and finite elements, the nodes follow force law and motion law of their own method, and on the their interaction interface, the touch type between discrete and finite elements is distinguished as two types: full touch and partial touch, the interaction force between them is calculated with linear elastic model. For full touch, the contact force is proportional to the overlap distance between discrete element and finite element patch. For partial touch, first the finite element patch is extended on all sides indefinitely to be a complete plane, the full contact force can be obtained with the touch type between discrete element and plane being viewed as full touch, then the full overlap area between them and the actual overlap area between discrete element and finite element patch are computed, the actual contact force is obtained by scaling the full contact force with a factor which is determined by the ratio of the actual overlap area to the full overlap area. The contact force is equivalent to the finite element nodes and the force and displacement on the nodes can be computed, so the ideal simulation results can be got. This method has been used to simulate the cutter disk of the earth pressure balance shield machine (EPBSM) made in North Heavy Industry (NHI) with its excavation diameter of 6.28 m cutting and digging the sandy clay layer. The simulation results show that as the gradual increase of excavating depth of the cutter head, the maximum stress occurs at the roots of cutters on the cutter head, while for the soil, the largest stress is distributed at the region which directly contacted with the cutters. The proposed research can provide good solutions for correct design and installation of cutters, and it is necessary to design mounting bracket to fix cutters on cutter head.

Two-scale finite element method for piezoelectric problem in periodic structure

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.

High-efciency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.

Study on spline wavelet finite-element method in multi-scale analysis for foundation

A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.

Prediction of grain scale plasticity of NiTi shape memory alloy based on crystal plasticity finite element method

Grain scale plasticity of NiTi shape memory alloy(SMA)during uniaxial compression deformation at 400℃was investigated through two-dimensional crystal plasticity finite element simulation and corresponding analysis based on the obtained orientation data.Stress and strain distributions of the deformed NiTi SMA samples confirm that there exhibits a heterogeneous plastic deformation at grain scale.Statistically stored dislocation(SSD)density and geometrically necessary dislocation(GND)density were further used in order to illuminate the microstructure evolution during uniaxial compression.SSD is responsible for sustaining plastic deformation and it increases along with the increase of plastic strain.GND plays an important role in accommodating compatible deformation between individual grains and thus it is correlated with the misorientation between neighboring grains,namely,a high GND density corresponds to large misorientation between grains and a low GND density corresponds to small misorientation between grains.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.

Dynamic Crack Propagation Analysis Using Scaled Boundary Finite Element Method

The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.

Modeling Technology in Finite Element Analysis of Electrostatic Proximity Fuze Problem

In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.

Finite Element Simulation of the Polymeric Liquid Flow in Coat-hanger Die

A three - dimenslonal finite element model is developed to deal with the polymeric liquid flow in coat - hanger die. This model is used to predict the flow behavior of the 2％ CMC/watsr solution in the coat - hanger die with linearly taper manifolds and its validity is experimentally verified quantitatively and qualitatively by using Laser Doppler Velocimetry and Particle Image Velocimetry respectively.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

Parametric study on single shot peening by dimensional analysis method incorporated with finite element method

Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.

常泰长江大桥组合索塔锚固结构钢-混传剪构造足尺模型试验研究

常泰长江大桥索塔锚固结构采用钢箱-核芯混凝土组合结构,为研究该新型组合索塔锚固结构钢-混传剪构造的受力特性,进行钢-混传剪构造足尺模型试验研究。制作2个锚固结构足尺节段试验模型,通过压剪试验研究锚固结构的荷载~滑移曲线及应力、应变分布等受力特性,并通过有限元模型分析锚固结构的传力机理和各组件的内力分配比例,推导剪力钉剪力计算方法。结果表明:在2.14倍单索最大索力荷载作用下,锚固结构保持弹性状态,钢壁板未产生明显滑移,钢-混界面最大滑移不超过0.25 mm,该锚固结构中钢-混传剪构造至少具有2.14倍的安全系数;荷载作用下,剪力钉剪力从上至下逐渐增大,锚腹板附近底部3排剪力钉剪力较大,钢-混传剪构造至少存在剪力钉和界面摩擦力2种传剪机制,钢-混传剪构造的承载能力显著提高;钢-混传剪构造受力过程分为粘结力传力阶段和局部滑移阶段,剪力钉剪力分布不仅与沿剪切方向长度分布有关,也与荷载的大小线性相关。

平钢板-UHPC组合桥面板纵向抗负弯性能试验研究

某大桥次边跨及中跨主梁为钢-UHPC组合梁,钢-UHPC组合梁的桥面板首次采用一种通过PBL剪力键将8 mm平钢板与15 cm的UHPC层连接起来的新型组合桥面体系。为了探究该桥新型钢-UHPC组合桥面板抗负弯矩性能与安全性,完成了2块足尺模型的静力性能试验,进行了桥梁整体受力和桥面板局部受力计算,以及桥面板抗负弯矩极限承载力构成分析。结果表明:UHPC名义开裂强度达8.90 MPa以上,其值远大于实桥荷载作用下UHPC层上缘的纵桥向最大拉应力,组合桥面板负弯矩抗裂性能良好,能满足工程需求;组合桥面板抗负弯承载力的计算,UHPC受拉本构宜采用三折线模型,抗负弯承载力简化计算的UHPC受拉区等效均布应力折减系数可取0.76,抗负弯承载力状态的受压钢板受力仍处于线弹性阶段,抗负弯破坏模式与抗正弯明显不同,为受拉钢筋拉断;PBL剪力键能保证了钢板与UHPC板整体共同受力;组合桥面板在负弯矩作用下的延性良好,裂宽增长缓慢,UHPC的名义极限强度达名义开裂强度的4.1倍以上,到达极限荷载后,组合板承载力没有明显下降。

正交异性钢-钢纤维混凝土组合桥面板疲劳极限状态研究

为了解正交异性钢-钢纤维混凝土(Steel Fiber Reinforced Concrete,SFRC)组合桥面板的疲劳性能及失效机理,以川南城际铁路临港长江大桥为背景,设计、制作正交异性钢-SFRC组合桥面板足尺模型进行疲劳试验,采用ANSYS软件建立试件有限元模型,研究关键细节的疲劳应力和开裂等情况。基于有限元模型,分析不同设计参数(钢顶板厚度、栓钉布置、SFRC抗拉强度及层厚)下组合桥面板疲劳极限状态的失效模式,并提出了主要控制参数取值建议。结果表明:200万次疲劳加载后足尺模型的实测最大裂缝宽度为0.136 mm,出现在SFRC层上缘,钢结构未开裂,组合桥面板疲劳性能良好;组合桥面板疲劳极限状态主要由SFRC层开裂控制,钢顶板厚度、栓钉布置对组合桥面板疲劳性能的影响较小,SFRC抗拉强度和层厚对组合桥面板疲劳性能影响较大,为主要控制参数;在常规正交异性钢桥面板上铺设薄层(厚度不超过50 mm)SFRC时,SFRC抗拉强度不应小于5 MPa,在常规正交异性钢桥面板上铺设普通SFRC(钢纤维体积含量不高于1%,抗拉强度不高于3 MPa)时,SFRC层厚不宜低于100 mm。

Research on Analytical Method of Fatigue Characteristics of Soft Yoke Mooring System Based on Full-Scale Measurement

By focusing on the vulnerability of the structure of marine equipments, together with considering the randomness of meta-ocean load in statistics, a kind of analytical method of fatigue characteristics of marine structure based on full- scale and actual measurement of prototype is proposed. On the basis of short-term field measurement results of structural response, research is carried out on the fatigue analysis of hinge joints at the upper part of the Soft Yoke single point Mooring System （SYMS） by simultaneously monitoring the environmental load and considering the design criteria of offshore structure. Through analysis of finite element modeling, the time-histories of typical stress response are obtained, and then the assessment of fatigue damage at key hinge joints is conducted. The simulation results indicate that the proposed method can accurately analyze the fatigue damage of offshore engineering structure caused by the effect of wave load. The present analytical method of fatigue characteristics can be extended on other offshore engineering structures subjected to meta-ocean load.