PHC管桩静压沉桩有限元分析

预应力高强度混凝土管桩(PHC管桩)具有应用范围广、单桩承载力高、工程造价低等鲜明特点,所以在我国发展迅速。根据南昌地区特殊的土层特点,依托该地区具体工程项目,采用有ANSYS有限元软件分析管桩在沉桩过程中土中应力应变分布及隆起特点。

Numerical simulation of two-dimensional magnetotelluric forward modeling of Heifangtai loess landslide in Gansu,China

Loess soils are widely distributed worldwide and typical in northwest China,and excessive agricultural irrigation has caused landslides in the area,specifi cally in the Heifangtai loess region in Lanzhou,Gansu,China.Geophysical exploration is an essential method in landslide engineering geological surveys,and geological surveying,drilling,geophysical prospecting,monitoring,and other methods are used for performing engineering geological evaluation and obtaining comprehensive basic data for landslide protection design and construction.The theoretical feasibility of using geophysical methods in loess landslide detection is essential.On the basis of the shallow geological structure of the Heifangtai landslide region in Lanzhou,Gansu,China,a typical geoelectric model of the magnetotelluric method was established,and the loess landslide area was modeled through a two-dimensional fi nite element method,forward numerical simulation,and engineering geological analysis.The distribution characteristics of the magnetotelluric field were determined.This is a typical application of the geological process analysis method in geophysical exploration.This study provides the typical stratigraphic structure and electrical characteristics of different groundwater distributions in Heifangtai,Gansu,China,verifi es the accuracy of forward modeling and calculation results,and provides a detailed theoretical basis for landslide detection through magnetotelluric methods.Through the numerical simulation of the forward modeling of the Heifangtai landslide region in Lanzhou,Gansu,China,this study can provide a detailed geophysical basis for landslide investigation,corroborate results of geological investigation and landslide design,and facilitate the sustainable development of agriculture in Heifangtai.

沉降控制复合桩基数值模拟分析

沉降控制复合桩基在上海软土地区建筑工程中得到了广泛应用,根据沉降控制复合桩基自身的特点,提出了简化的三维有限元分析模型,模拟了桩土及其承台之间复杂的相互作用,通过计算分析,得到了若干可供工程设计参考的结论。

Finite element numerical simulation of 2.5D direct current method based on mesh refinement and recoarsement

To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.

Dynamic analysis of pavement on long span steel bridge decks

In order to analyze the dynamic response of pavement on long-span steel bridge decks under random dynamic loads, the irregularities of the pavement surface is simulated with the power spectrum density function, and the random load is calculated according to a vehicle vibration equation of vehicle model. The mechanical responses of three different cases are compared by using a transient dynamic analysis method, i. e., under random dynamic load, constant moving load and dead load respectively. The results indicate that the mid-span of two adjacent transversal diaphragms is the worst load position. The maximum vertical displacement and the maximum transversal tensile stress of the pavement are 1.33 times and 1.39 times as much as those when only considering the impact coefficients. This study not only provides a theoretical basis for the mixture design and structural design of pavement, but also puts forward higher demand on the construction and maintenance for steel deck pavement.

Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data

The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.

Two-dimensional frequency-domain acoustic fullwaveform inversion with rugged topography

We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to solve the cutoff botmdary problem as well as to consider the requirement of using the same subdivision grid in joint multifrequency inversion. The proposed method introduces the attenuation factor, and by adjusting it, acoustic waves are sufficiently attenuated in the attenuation layer to minimize the cutoff boundary effect. Based on the law of exponential attenuation, expressions for computing the attenuation factor and the thickness of attenuation layers are derived for different frequencies. In multifrequency-domain FWI, the conjugate gradient method is used to solve equations in the Gauss-Newton algorithm and thus minimize the computation cost in calculating the Hessian matrix. In addition, the effect of initial model selection and frequency combination on FWI is analyzed. Examples using numerical simulations and FWI calculations are used to verify the efficiency of the proposed method.

Nonlinear dynamical response of high-voltage transmission lines based on cable dropping

In order to study the dynamic response of high-voltage transmission lines under mechanical failure, a finite element model of a domestic 500-kV high-voltage transmission line system is established. The initial equilibrium condition of the coupling system model is verified by nonlinear static analysis. The transient dynamic analysis method is proposed to analyze the variation law of dynamic response under cable or insulator rupture, and the dynamic response of structural elements next to the broken span is calculated. The results show that upper crossarm cable rupture has no effect on cable tension at adjacent suspension points, but it has a significant influence on tension in the insulator and the tower component of the upper crossarm next to the broken span. The peak tension in the conductor of the upper crossarm at the suspension point exceeds the design value under insulator rupture. Insulator rupture has no effect on the tower component of the upper crossarm, but it has a significant influence on insulator tension of the upper crossarm. Insulator rupture should be taken into account in the design of overhead transmission lines. The research results can provide a theoretical basis for the design of transmission lines.

NEW LIMITED MEMORY SYMMETRIC RANK ONE ALGORITHM FOR LARGE-SCALE UNCONSTRAINED OPTIMIZATION

A new limited memory symmetric rank one algorithm is proposed. It combines a modified self-scaled symmetric rank one （SSR1） update with the limited memory and nonmonotone line search technique. In this algorithm, the descent search direction is generated by inverse limited memory SSR1 update, thus simplifying the computation. Numerical comparison of the algorithm and the famous limited memory BFGS algorithm is given. Comparison results indicate that the new algorithm can process a kind of large-scale unconstrained optimization problems.

Residual stress and welding distortion of Al/steel butt joint by arc-assisted laser welding-brazing

The thermo-elastic.plastic finite element method(FEM)is used to simulate the thermo-mechanical behavior of Al/steel tungsten inert gas(TIG)arc-assisted laser welding-brazing(A-LWB)butt joint.The influence of material nonlinearity,geometrical nonlinearity and work hardening on the welding process is studied,and the differences in the welding temperature field,residual stress and welding distortion by A-LWB and by single laser welding-brazing(SLWB)are analyzed.The results show that the thermal cycle,residual stress distribution and welding distortion by the numerical simulation are in good agreement with the measured data by experiments,which verifies the effectiveness of FEM.Compared with the SLWB,A-LWB can make the high-temperature distribution zone of weld in width direction wider,decrease the transverse tensile stress in the weld and reduce the distribution range of longitudinal tensile stress.And the welding deformation also decreases to some extent.

Numerical modeling of anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods

Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal axes symmetry,respectively.We then provide a simple derivation of the anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods.Analogous to the mean resistivity,we propose a formulation for deriving the mean polarizability.We also present a three-dimensional finite element algorithm for modeling the direct current resistivity and time-domain induced polarization using an unstructured tetrahedral grid.Finally,we provide the apparent resistivity and apparent chargeability curves of a tilted,transversely isotropic medium with diff erent angles,respectively.The subsequent results illustrate the anisotropy paradoxes of direct current resistivity and time-domain induced polarization.

Dynamic analysis of a flexible rotor supported by ball bearings with damping rings based on FEM and lumped mass theory

A dynamic model of a flexible rotor supported by ball bearings with rubber damping rings was proposed by combining the finite element and the mass-centralized method.In the proposed model,the rotor was built with the Timoshenko beam element,while the supports and bearing outer rings were modelled by the mass-centralized method.Meanwhile,the influences of the rotor’s gravity,unbalanced force and nonlinear bearing force were considered.The governing equations were solved by precise integration and the Runge-Kutta hybrid numerical algorithm.To verify the correctness of the modelling method,theoretical and experimental analysis is carried out by a rotor-bearing test platform,where the error rate between the theoretical and experimental studies is less than 10%.Besides that,the influence of the rubber damping ring on the dynamic properties of the rotor-bearing coupling system is also analyzed.The conclusions obtained are in agreement with the real-world deployment.On this basis,the bifurcation and chaos behaviors of the coupling system were carried out with rotational speed and rubber damping ring’s stiffness.The results reveal that as rotational speed increases,the system enters into chaos by routes of crisis,quasi-periodic and intermittent bifurcation.However,the paths of crisis,quasi-periodic bifurcation,and Hopf bifurcation to chaos were detected under the parameter of rubber damping ring’s stiffness.Additionally,the bearing gap affects the rotor system’s dynamic characteristics.Moreover,the excessive bearing gap will make the system’s periodic motion change into chaos,and the rubber damping ring’s stiffness has a substantial impact on the system motion.

Improved Particle Swarm Optimization for Solving Transient Nonlinear Inverse Heat Conduction Problem in Complex Structure

Accurately solving transient nonlinear inverse heat conduction problems in complex structures is of great importance to provide key parameters for modeling coupled heat transfer process and the structure’s optimization design.The finite element method in ABAQUS is employed to solve the direct transient nonlinear heat conduction problem.Improved particle swarm optimization(PSO)method is developed and used to solve the transient nonlinear inverse problem.To investigate the inverse performances,some numerical tests are provided.Boundary conditions at inaccessible surfaces of a scramjet combustor with the regenerative cooling system are inversely identified.The results show that the new methodology can accurately and efficiently determine the boundary conditions in the scramjet combustor with the regenerative cooling system.By solving the transient nonlinear inverse problem,the improved particle swarm optimization for solving the transient nonlinear inverse heat conduction problem in a complex structure is verified.

Analytic and Experimental Studies of the Errors in Numerical Methods for the Valuation of Options

The value of a European option satisfies the Black-Scholes equation with appropriately specified final and boundary conditions.We transform the problem to an initial boundary value problem in dimensionless form.There are two parameters in the coefficients of the resulting linear parabolic partial differential equation.For a range of values of these parameters,the solution of the problem has a boundary or an initial layer.The initial function has a discontinuity in the first-order derivative,which leads to the appearance of an interior layer.We construct analytically the asymptotic solution of the equation in a finite domain.Based on the asymptotic solution we can determine the size of the artificial boundary such that the required solution in a finite domain in x and at the final time is not affected by the boundary.Also,we study computationally the behaviour in the maximum norm of the errors in numerical solutions in cases such that one of the parameters varies from finite (or pretty large) to small values,while the other parameter is fixed and takes either finite (or pretty large) or small values. Crank-Nicolson explicit and implicit schemes using centered or upwind approximations to the derivative are studied.We present numerical computations,which determine experimentally the parameter-uniform rates of convergence.We note that this rate is rather weak,due probably to mixed sources of error such as initial and boundary layers and the discontinuity in the derivative of the solution.

Finite Element Analysis of Polygon Shaped Shell Roof

A new spherical triangular finite element based on shallow shell formulation is developed in this paper. The element has six degrees of freedom at each comer node, five of which are the essential external degrees of freedom and the additional sixth is associated with the in-plane shell rotation. The displacement fields of the element satisfy the exact requirement of rigid body modes of motion. The element is based on independent strain assumption insofar as it is allowed by the compatibility equations. The element developed herein is first validated by applying it to the analysis of a benchmark problem involving a standard spherical shell with simply supported edges. The results of the analysis showed that reasonably accurate results were obtained even when modeling the shells using fewer elements compared to other shell element types. The element is then used in a finite element model to analyze polygon shaped spherical roof structures. The distribution of the various components of deflection and stress is obtained. Furthermore, the effect of introducing circular arched beams as stiffeners spanning the two diagonally opposite end comers is investigated. It is found that the stiffeners reduced the deflections and the stresses in the roof structure by considerable value.

Comparison on construction of strut-and-tie models for reinforced concrete deep beams

With consideration of the differences between concrete and steel,three solutions using genetic evolutionary structural optimization algorithm were presented to automatically develop optimal strut-and-tie model for deep beams.In the finite element analysis of the first method,the concrete and steel rebar are modeled by a plane element and a bar element,respectively.In the second method,the concrete and steel are assigned to two different plane elements,whereas in the third method only one kind of plane element is used with no consideration of the differences of the two materials.A simply supported beam under two point loads was presented as an example to verify the validity of the three proposed methods.The results indicates that all the three methods can generate optimal strut-and-tie models and the third algorithm has powerful capability in searching more optimal results with less computational effort.The effectiveness of the proposed algorithm III has also been demonstrated by other two examples.

Rock Cutting Analysis Employing Finite and Discrete Element Methods

The petroleum industry has shown great interest in the study of drilling optimization on pre-salt formations given the low rates of penetration observed so far. Rate of penetration is the key to economically drill the pre-salt carbonate rock. This work presents the results of numerical modeling through finite element method and discrete element method for single cutter drilling in carbonate samples. The work is relevant to understand the mechanics of drill bit-rock interaction while drilling deep wells and the results were validated with experimental data raised under simulated downhole conditions. The numerical models were carried out under different geometrical configurations, varying the cutter chamfer size and back-rake angles. The forces generated on the cutter are translated into mechanical specific energy as this parameter is often used to measure drilling efficiency. Results indicate that the chamfer size does not change significantly the mechanical specific energy values, characteristic. Results also show there is a significant increase although the cutter aggressiveness is influenced by this geometrical in drilling resistance for larger values of back-rake angle.

Torsional inertia moment of beam element with complex section analysis based on FEM

Currently, for the analysis of complex bridge based on beam element, the calculation of cross-section torsional inertia moment is still an unresolved technical problem. Most current calculation of section torsional inertia moment is an approximate analytic method for two-dimensional cross-section, which is not fully consistent with the actual situation, and do not consider the effects of diaphragm in bridge. In order to analyze accurately cable-stayed bridge, suspension bridge and other complex bridge structures based on beam element, the calculation method of section torsional inertia moment based on finite element method (FEM) is invented in this paper. Firstly, setting up local cantilever fine model with solid element or shell element and applying torsion on the end of cantilever. Secondly, calculating the torsion angle of cantilever by FEM method and then the torsional moment through equivalent beam method. Finally, the examples of the section torsional moment calculation of concrete model with solid element with diaphragm and steel girder box model with shell element with diaphragm are used to verify the calculation method, which is applied to the suspension bridge design and construction control special software SBNA developed by Research Institute of Highway Ministry of Transport. Taizhou Bridge under construction is one of the examples.

Prediction Model for Cooling of an Electrical Unit with Time-Dependent Heat Generation

The satisfactory performance of electrical equipments depends on their operating temperature. In order to maintain these devices within the safe temperature limits, an effective cooling is needed. High heat transfer rate of compact in size and reliable operation are the challenges of a thermal design engineer of electronic equipment. Then, it has been simulated the transient a three-dimensional model to study the heating phenomenon with two assumption values of heat generation. To control for the working of this equipment, cooling process was modeled by choosing one from different cooling technique. Constant low speed fan at one direction of air flow was used for cooling to predict the reducing of heating temperature through working of this equipment. Numerical Solution of finite difference time domain method （FDTD） has been utilized to simulate the temporal and spatial temperature profiles through two processes, which would minimize the solution errors.