A GEOMETRIC FEATURE FOR FINITE ELEMENT SCHEMES

In this note we characterize the geometric facture of a (μ, r, k) - FES. Namely, for a C~μ triangular in- terpolation schcme with C^r vertex data, any angle of the macrotriangle must be divided into at least (μ + 1)/ (r + 1 - μ) parts.

Relationships between femoral strength evaluated by nonlinear finite element analysis and BMD,material distribution and geometric morphology

Precisely quantifying the strength of the proximal femur and accurately assessing hip fracture risk would enable those at high risk to be identified so that preventive interventions could be taken.Development of better measures of femoral strength using the clinically

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.

Frame-invariance in finite element formulations of geometrically exact rods

This article is concerned with finite element implementations of the three- dimensional geometrically exact rod. The special attention is paid to identifying the con- dition that ensures the frame invariance of the resulting discrete approximations. From the perspective of symmetry, this requirement is equivalent to the commutativity of the employed interpolation operator I with the action of the special Euclidean group SE（3）, or I is SE（3）-equivariant. This geometric criterion helps to clarify several subtle issues about the interpolation of finite rotation. It leads us to reexamine the finite element for- mulation first proposed by Simo in his work on energy-momentum conserving algorithms. That formulation is often mistakenly regarded as non-objective. However, we show that the obtained approximation is invariant under the superposed rigid body motions, and as a corollary, the objectivity of the continuum model is preserved. The key of this proof comes from the observation that since the numerical quadrature is used to compute the integrals, by storing the rotation field and its derivative at the Gauss points, the equiv- ariant conditions can be relaxed only at these points. Several numerical examples are presented to confirm the theoretical results and demonstrate the performance of this al- gorithm.

Nonlinear anisotropic finite element analysis of liquefiable tunnel–sand–pile interaction under seismic excitation

Nonlinear time‐history analysis can be used to determine the liquefiable behaviors of the tunnel-sand-pile interaction(TSPI)model with the consideration of sand anisotropy.This study presents the nonlinear response of the TSPI model with the existence of liquefaction under seismic excitation.The analysis reveals that tunnel and pile behave as isotropic elements,while sand shows isotropic,orthotropic,and anisotropic characteristics.Three constitutive models including UBC3D‐PLM(two yield surfaces associated with the hardening rule),NGI‐ADP(yielding with associated plastic potential function),and a user‐specified constitutive model are adopted to evaluate the isotropic,orthotropic,and anisotropic behaviors of sand.On this basis,two finite element‐based codes(ETABS 18.1.1 and Plaxis 3D)are used to evaluate sand behaviors and responses.Responses of the tunnel,sand,pile,and excess pore pressure ratio are recorded in the interaction zone by varying the pile diameter,tunnel diameter,and tunnel-pile clearance.Compared with the orthotropic and isotropic conditions,lower variations of results are found in the anisotropic condition,except for the case of generation of excess pore pressure.In addition,the present reanalysis results are in agreement with previous analytical and case study results,which further indicates the effectiveness of the finite element‐based numerical codes.

Vector form Intrinsic Finite Element Method for the Two-Dimensional Analysis of Marine Risers with Large Deformations

Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static and dynamic analyses of marine risers. The proposed approach uses the vector form intrinsic finite element(VFIFE) method, which is based on vector mechanics theory and numerical calculation. In this method, the risers are described by a set of particles directly governed by Newton's second law and are connected by weightless elements that can only resist internal forces. The method does not require the integration of the stiffness matrix, nor does it need iterations to solve the governing equations. Due to these advantages, the method can easily increase or decrease the element and change the boundary conditions, thus representing an innovative concept of solving nonlinear behaviors, such as large deformation and large displacement. To prove the feasibility of the VFIFE method in the analysis of the risers, rigid and flexible risers belonging to two different categories of marine risers, which usually have differences in modeling and solving methods, are employed in the present study. In the analysis, the plane beam element is adopted in the simulation of interaction forces between the particles and the axial force, shear force, and bending moment are also considered. The results are compared with the conventional finite element method(FEM) and those reported in the related literature. The findings revealed that both the rigid and flexible risers could be modeled in a similar unified analysis model and that the VFIFE method is feasible for solving problems related to the complex behaviors of marine risers.

A lumped mass nonconforming finite element method for nonlinear parabolic integro-differential equations on anisotropic meshes

A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.

ALE FINITE ELEMENT ANALYSIS OF NONLINEAR SLOSHING PROBLEMS BASED ON F LUID VELOCITY POTENTIAL

Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.

Assessment of natural frequency of installed offshore wind turbines using nonlinear finite element model considering soil-monopile interaction

A nonlinear finite element model is developed to examine the lateral behaviors of monopiles, which support offshore wind turbines(OWTs) chosen from five different offshore wind farms in Europe. The simulation is using this model to accurately estimate the natural frequency of these slender structures, as a function of the interaction of the foundations with the subsoil. After a brief introduction to the wind power energy as a reliable alternative in comparison to fossil fuel, the paper focuses on concept of natural frequency as a primary indicator in designing the foundations of OWTs. Then the range of natural frequencies is provided for a safe design purpose. Next, an analytical expression of an OWT natural frequency is presented as a function of soil-monopile interaction through monopile head springs characterized by lateral stiffness KL, rotational stiffness KRand cross-coupling stiffness KLRof which the differences are discussed. The nonlinear pseudo three-dimensional finite element vertical slices model has been used to analyze the lateral behaviors of monopiles supporting the OWTs of different wind farm sites considered. Through the monopiles head movements(displacements and rotations), the values of KL, KRand KLRwere obtained and substituted in the analytical expression of natural frequency for comparison. The comparison results between computed and measured natural frequencies showed an excellent agreement for most cases. This confirms the convenience of the finite element model used for the accurate estimation of the monopile head stiffness.

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.

Analytical study of performance evaluation for seismic retrofitting of reinforced concrete building using 3D dynamic nonlinear finite element analysis

This paper presents three-dimensional finite element （FE） analyses of an all-frame model of a three-story reinforced concrete （RC） building damaged in the 1999 Taiwan Chi-Chi Earthquake. Non-structural brick walls of the building acted as a seismic resistant element although their contributions were neglected in the design. Hence, the entire structure of a typical frame was modeled and static and dynamic nonlinear analyses were conducted to evaluate the contributions of the brick walls. However, the results of the analyses were considerably overestimated due to coarse mesh discretizations, which were unavoidable due to limited computer resources. This study corrects the overestimations by modifying （1） the tensile strengths and （2） shear stiffness reduction factors of concrete and brick. The results indicate that brick walls improve frame strength although shear failures are caused in columns shortened by spandrel walls. Then, the effectiveness of three types of seismic retrofits is evaluated. The maximum drift of the first floor is reduced by 89.3%, 94.8%, and 27.5% by Steel-confined, FuI1-RC, and Full-brick models, respectively. Finally, feasibility analyses of models with soils were conducted. The analyses indicated that the soils elongate the natural period of building models although no significant differences were observed.

Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

A review of flexibility-based finite element method for beam-column elements

For material nonlinear problem,elements derived with the flexibility-based method are more accurate than classical elements derived with the stiffness-based method. A review of the current state of the art of the flexibility-based finite element method is provided to enhance the robustness of structure analysis. The research on beam-column elements is the mainstream in the research on flexibility-based finite element method at present. The original development of flexibility-based finite element method is reviewed, and the further development of this method is then presented in several specific aspects, such as geometrically nonlinear analysis and dynamic analysis. The further research needed to be carried out in the future is finally discussed.

PRINCIPAL COMPONENT DECOMPOSITION BASED FINITE ELEMENT MODEL UPDATING FOR STRAIN-RATE-DEPENDENCE NONLINEAR DYNAMIC PROBLEMS

Thin wail component is utilized to absorb impact energy of a structure. However, the dynamic behavior of such thin-walled structure is highly non-linear with material, geometry and boundary non-linearity. A model updating and validation procedure is proposed to build accurate finite element model of a frame structure with a non-linear thin-walled component for dynamic analysis. Design of experiments （DOE） and principal component decomposition （PCD） approach are applied to extract dynamic feature from nonlinear impact response for correlation of impact test result and FE model of the non-linear structure. A strain-rate-dependent non-linear model updating method is then developed to build accurate FE model of the structure. Computer simulation and a real frame structure with a highly non-linear thin-walled component are employed to demonstrate the feasibility and effectiveness of the proposed approach.

Seismic Deformation and Seismic Resistance Analysis of Shapai Roller Compacted Concrete Arch Dam Based on Field Monitoring and Dynamic Finite Element Method

Shapai Roller Compacted Concrete(RCC) Arch Dam is the highest RCC arch dam of the 20th century in the world with a maximum height of 132m,and it is the only concrete arch dam near the epicentre of Wenchuan earthquake on May 12th,2008.The seismic damage to the dam and the resistance of the dam has drawn great attention.This paper analyzed the response and resistance of the dam to the seismic wave using numerical simulations with comparison to the monitored data.The field investigation after the earthquake and analysis of insitu data record showed that there was only little variation in the opening size at the dam and foundation interface,transverse joints and inducing joints before and after the earthquake.The overall stability of the dam abutment resistance body was quite good except a little relaxation was observed.The results of the dynamic finite element method(FEM) showed that the sizes of the openings obtained from the numerical modeling are comparable with the monitored values,and the change of the opening size is in millimeter range.This study revealed that Shapai arch dam exhibited high seismic resistance and overload capacity in the Wenchuan earthquake event.The comparison of the monitored and simulated results showed that the numerical method applied in this paper well simulated the seismic response of the dam.The method could be useful in the future application on the safety evaluation of RCC dams.

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

RANDOM MICROSTRUCTURE FINITE ELEMENT METHOD FOR EFFECTIVE NONLINEAR PROPERTIES OF COMPOSITE MATERIALS

Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.

NONLINEAR QUASI-CONFORMING FINITE ELEMENT METHOD

The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given.

Nonlinear Finite Element Analysis of Steel Fiber Reinforced Concrete Deep Beams

By the nonlinear finite element analysis （FEA） method, the mechanical properties of the steel fiber reinforced concrete （SFRC） deep beams were discussed in terms of the crack load and ultimate bearing capacity. In the simulation process, the ANSYS parametric design language （APDL） was used to set up the finite element model; the model of bond stress-slip relationship between steel bar and concrete was established. The nonlinear FEA results and test results demonstrated that the steel fiber can not only significantly improve the cracking load and ultimate bearing capacity of the concrete but also repress the development of the cracks. Meanwhile, good agreement was found between the experimental data and FEA results, if the unit type, the parameter model and the failure criterion are selected reasonably.

Nonlinear Finite Element Analysis of the Structure of Door Seals

In order to evaluate the influence of the seal structure on door dosing force, nonlinear finite dement methed is introduced to analyze compression deformation of a door seal for SANTANA （name of the car made by Shanghai Volkswagen Co. Ltd）. MSC. Marc software is used to analyze the large deformation of the seal and the compression test is done to prove the computational results. The results show that the compression loads of the door seal are larger than the standard value of Shanghai Volkswagen Co. Ltd and the seal structure needs to be optimized. There are consistent relationships between calculating results and experimental results and the simulation method is effective.