A finite element optimization technique to determine critical imperfections of shell structures

In this paper, the numerical model developed in the previous paper is used to study the seismic performance of elevated liquid-filled steel conical tanks. A number of conical tanks which are classified as tall or broad tanks according to the ratio of the tank radius to its height are considered. The consistent shell element is used to model the tank surfaces, while the coupled boundary-shell element formulation is employed to obtain the fluid added-mass which simulates the dynamic pressure resulting from a seismic motion. Linear springs are used to model the supporting towers. The natural frequencies of the liquid-filled tanks due to both horizontal and vertical excitations are evaluated. This is followed by a non-linear dynamic analysis, using an appropriately scaled real input ground motion, and which includes the effect of both geometric and material non-linearities. Thin-walled structures of this kind may exhibit inelastic behaviour and a tendency to develop localized buckles, thus diminishing stiffness. The consequence could lead to overall instability of the structure. In general, time-history analyses indicate that liquid-filled conical tanks, often possessing apparently adequate safety factors under hydrostatic loading, are shown to be very sensitive to seismic loading when ground motion frequencies contain those of the fundamental frequencies of the vessels themselves. © 1997 John Wiley & Sons, Ltd.

The cable tower is the bearing component of long-span suspension bridges, and its structure is very high and bear large force, which determines the stability and is the key of safety control. As for the height of the main tower of a long-span suspension bridge up to 195.3 m, the finite element software ANSYS is used to establish a three-dimensional finite element model (FEM), and the effects of geometric nonlinearity and material nonlinearity on the stability of the main tower are analyzed. The calculation results show that geometrical nonlinearity and material defects have significant influence on the main tower stability, and the nonlinear stability should be considered under wind load in the design calculation.

In this paper, the mathematical model of a jet pipe electrohydraulic flow control servo valve with a built-in mechanical feedback is developed. The actuator dynamics is also included. The solid and finite element (FE) models of the critical and precision components, feedback spring assembly and flexure tube of jet pipe assembly, were carried out first and simulated with appropriate boundary and loading conditions in a commercially available FE code, IDEAS, version 7.0, to ascertain the stiffnesses. To validate the obtained finite element results, experiments were carried out with a specially designed and fabricated test set-up. Moderately good agreement was observed between finite element and experimental values. The required simulation parameters for a jet pipe servo valve-actuator system, such as the moment of inertia, mass of the spool, mass of the piston, stiffness of the critical and precision components, and lever arm length, were obtained from the solid and finite element models. The developed model was simulated in commercially available software, MATLAB 6.5, Release 13. The step response of the valve was analysed. The simulation clearly showed the dynamics of the jet pipe, spool, and actuator displacement in achieving the equilibrium position when the applied torque and restoring torques balance, known as the steady state.

In the framework of finite elements (FEs) applications, this paper proposes the use of the node-dependent kinematics (NDK) concept to the large deflection and post-buckling analysis of thin-walled metallic one-dimensional (1D) structures. Thin-walled structures could easily exhibit local phenomena which would require refinement of the kinematics in parts of them. This fact is particularly true whenever these thin structures undergo large deflection and post-buckling. FEs with kinematics uniform in each node could prove inappropriate or computationally expensive to solve these locally dependent deformations. The concept of NDK allows kinematics to be independent in each element node; therefore, the theory of structures changes continuously over the structural domain. NDK has been successfully applied to solve linear problems by the authors in previous works. It is herein extended to analyze in a computationally efficient manner nonlinear problems of beam-like structures. The unified 1D FE model in the framework of the Carrera Unified Formulation (CUF) is referred to. CUF allows introducing, at the node level, any theory/kinematics for the evaluation of the cross-sectional deformations of the thin-walled beam. A total Lagrangian formulation along with full Green–Lagrange strains and 2nd Piola Kirchhoff stresses are used. The resulting geometrical nonlinear equations are solved with the Newton–Raphson linearization and the arc-length type constraint. Thin-walled metallic structures are analyzed, with symmetric and asymmetric C-sections, subjected to transverse and compression loadings. Results show how FE models with NDK behave as well as their convenience with respect to the classical FE analysis with the same kinematics for the whole nodes. In particular, zones which undergo remarkable deformations demand high-order theories of structures, whereas a lower-order theory can be employed if no local phenomena occur: this is easily accomplished by NDK analysis. Remarkable advantages are shown in the analysis of thin-walled structures with transverse stiffeners.

A composite finite element model for cold bend simulation of energy pipelines is presented in this paper. Four-node shell elements with material and geometric nonlinearity are used to model a pipe in straight condition. An elastic pipe, having the same nodal coordinates as the main pipe along with elastic radial links are used as a tool to prevent local buckling and ovalization of the main pipe during the cold bend process. By dividing the elastic pipe into series of rings along the axis of the pipe and by conducting a four-step procedure, residual curvature is developed in a specific segment of a pipe. Based on the proposed concept, different methods of cold bending are discussed and the results are presented. University of Alberta cold bend trials were used to validate the proposed finite element model. The moment-curvature response, pattern of imperfections, and distribution of maximum residual strains are obtained by the finite element model and compared with the test results.

To determine nonlinear material properties of passive, diastolic myocardium using magnetic resonance imaging (MRI) tissue-tagging, finite element analysis (FEA) and nonlinear optimization. Alterations in the diastolic material properties of myocardium may pre-date the onset of or exist exclusive of systolic ventricular dysfunction in disease states such as hypertrophy and heart failure. Accordingly, significant effort has been expended recently to characterize the material properties of myocardium in diastole. The present study defines a new technique for determining material properties of passive myocardium using finite element (FE) models of the heart, MRI tissue-tagging and nonlinear optimization. This material parameter estimation algorithm is employed to estimate nonlinear material parameter sin the in vivo canine heart and provides the necessary framework to study the full complexities of myocardial material behavior in health and disease. Material parameters for a proposed exponential strain energy function were determined by minimizing the least squares difference between FE model-predicted and MRI-measured diastolic strains. Six mongrel dogs underwent MRI imaging with radiofrequency (RF) tissue-tagging. Two-dimensional diastolic strains were measured from the deformations of the MRI tag lines. Finite element models were constructed from early diastolic images and were loaded with the mean early to late left ventricular and right ventricular diastolic change in pressure measured at the time of imaging. A nonlinear optimization algorithm was employed to solve the least squares objective function for hte material parameters. Average material parameters for the six dogs were E = 28,722 +/- 15984 dynes/cm2 and c = 0.00182 +/- 0.00232 cm2/dyne. This parameter estimation algorithm provides the necessary framework for estimating the nonlinear, anisotropic and non-homogeneous material properties of passive myocardium in health and disease in the in vivo beating heart.

In this paper, a finite element model updating (FEMU) method is proposed based on the response surface model (RSM) and genetic algorithm (GA) to establish a high-precision finite element (FE) model of space station scientific experiment racks. First, the fine solid and mixed FE models are established, respectively, and a comparison of the modal test results is conducted. Then, an orthogonal experimental design is used to analyze the significance of the parameters, and the variables to be modified are determined. The design parameters are sampled via the Latin hyperbolic method and are substituted into the FE model to obtain the modal parameters of the scientific experiment rack. The mapping relationship between the design and modal parameters is fitted by constructing the Kriging function, and the RSM is established. The design parameters of the scientific experiment rack are optimized via GA, and the initial FE model is updated, which has the advantage of improving the computing efficiency. Finally, the updated FE model of the experiment rack is verified by frequency sweep and random vibration tests. The experimental results show that the proposed approach has high precision and computing efficiency, and compared with the test results, the modal frequency errors of the updated model are within 5%, and the vibration response errors under random excitation of the updated model are within 7%.

Skull stiffness of children is of great significance to protect brain tissues during traffic accidents. In this paper, the effects of skull stiffness on traumatic brain injury (TBI) are investigated by using a validated 6-year-old child head finite element (FE) model according to Chinese standard GB/T 24550-2009 Results showed that the HIC value, intracranial pressure, Von Mises stress and shear strain of brain tissue increased with the decrease of skull stiffness during the impact of FE model with engine hood. Therefore, the injury risk of the child head decreases with the increase of skull stiffness. And the skull with higher stiffness can provide a better protection for brain tissues during traffic accidents.

There has been a growing demand for nonlinear model updating procedures in the structural dynamics community as advanced vehicles have started to incorporate nonlinearity into their designs. Finite element (FE) model updating is difficult for a nonlinear structure for several reasons: there may be many unknown parameters in the FE model so that multiple solutions may exist, and it is very expensive to compute the structure’s nonlinear normal modes from a full FE model that are an excellent metric to use for nonlinear updating. A recent work showed that some of these challenges can be overcome by updating a nonlinear reduced order model (ROM) rather than the full FE model. The updated ROM can be used to compute response statistics, stresses and ultimately life of the structure. However, one drawback of this approach is that one does not gain physical insight into which parameters in the FE model were in error, and so it is difficult to transfer the lessons learned to future FE models. This work explores the feasibility of updating a FE model to correlate with a ROM with a set of known parameters. The target ROM parameters consist of the linear stiffness and the nonlinear stiffness coefficients which have been identified previously by updating the ROM to correlate with experimental measurements. In this process an optimization routine is setup in which the free parameters are those of the FE model, such as boundary stiffness springs, imperfections, pre-stress, etc. The optimization procedure is wrapped around a ROM creation algorithm that iteratively tunes the FE parameters, creates a ROM, and evaluates the ROM parameters with respect to the target ROM in order to minimize the objective function. This work will demonstrate the procedure on a numerical example of a curved beam in order to verify its effectiveness on the nonlinear model updating.

Background Accurate biomechanical investigation of human intervertebral discs (IVDs) is difficult because of their complicated structural and material features. Aim To investigate probabilistic distributions of the biomechanical responses of the IVD by considering varying nonlinear structural and material properties using a stochastic finite element (FE) model. Methods A FE model of a L3-4 disc was reconstructed, including the nucleus pulposus (NP), annular matrix and fibers. A Monte Carlo method was used to randomly generate 500 sets of the nonlinear material properties and fiber orientations of the disc that were implemented into the FE model. The FE model was analyzed under seven loading conditions: a 500 N compressive force, a 7.5Nm moment simulating flexion, extension, left-right lateral bending, and left-right axial rotation, respectively. The distributions of the ranges of motion (ROMs), intradiscal pressures (IDP), fiber stresses and matrix strains of the disc were analyzed. Results Under the compressive load, the displacement varied between 0.29 mm and 0.76 mm. Under the 7.5Nm moment, the ROMs varied between 3.0° and 6.0° in primary rotations. The IDPs varied within 0.3 MPa under all the loading conditions. The maximal fiber stress (3.22 ± 0.64 MPa) and matrix strain (0.27 ± 0.12%) were observed under the flexion and extension moments, respectively. Conclusion The IVD biomechanics could be dramatically affected by the structural and material parameters used to construct the FE model. The stochastic FE model that includes the probabilistic distributions of the structural and material parameters provides a useful approach to analyze the statistical ranges of the biomechanical responses of the IVDs.

Carbon/carbon composites (C/Cs) are extensively utilized as structural materials and functional materials in the aerospace industry. Laser processing technology is an effective means of precision manufacturing C/Cs parts, with the advantages of no mechanical impact and high efficiency. Accurately predicting the material removal of C/Cs during laser processing is of great significance for the precision manufacturing of C/Cs parts. However, the numerical simulation models that can directly display the microstructure of C/Cs are still inadequate, and measuring the sublimation temperatures of two different phases of carbon is challenging. This paper establishes a three-dimensional microscopic heterogeneous finite element (FE) model of C/Cs, and the FE simulation of quasi-continuous wave (QCW) laser ablation of C/Cs is optimized using the restart method taking into account the residual temperature. Combining the optimized FE model, the material parameters of C/Cs are inverted using response surface methodology and genetic algorithm, resulting in the sublimation temperatures of the fiber phase being 4029.01 K and the matrix phase being 3481.86 K. After these parameters are substituted into the FE model, the resulting simulations are then compared with the experiments of QCW laser processing C/Cs, which reveals high correspondence between simulated morphology and experimental data, with the relative error of predicted ablation depth not exceeding 6.169%. The revised FE model can guide the laser processing of C/Cs, and the inverted material parameters can provide references for the theoretical study of the laser processing of C/Cs.

The aim of the present paper is the development of endodontic Cu-based single crystal Shape Memory Alloy (SMA) instruments in order to eliminate the antimicrobial and mechanical deficiencies observed with the conventional Nickel-Titane (NiTi) SMA files. A thermomechanical constitutive law, already developed and implemented in a finite element code by our research group, is adopted for the simulation of the single crystal SMA behavior. The corresponding material parameters were identified starting from experimental results for a tensile test at room temperature. A computer-aided design geometry has been achieved and considered for a finite element structural analysis of the endodontic Cu-based single crystal SMA files. They are meshed with tetrahedral continuum elements to improve the computation time and the accuracy of results. The geometric parameters tested in this study are the length of the active blade, the rod length, the pitch, the taper, the tip diameter, and the rod diameter. For each set of adopted parameters, a finite element model is built and tested in a combined bending-torsion loading in accordance with ISO 3630-1 norm. The numerical analysis based on finite element procedure allowed purposing an optimal geometry suitable for Cu-based single crystal SMA endodontic files. The same analysis was carried out for the classical NiTi SMA files and a comparison was made between the two kinds of files. It showed that Cu-based single crystal SMA files are less stiff than the NiTi files. The Cu-based endodontic files could be used to improve the root canal treatments. However, the finite element analysis brought out the need for further investigation based on experiments.

Finite element models (FEMs) including characteristic large deformations in highly nonlinear materials (hyperelasticity and coupled diffusive/convective transport of neutral mobile species) will allow quantitative study of in vivo tissues. Such FEMs will provide basic understanding of normal and pathological tissue responses and lead to optimization of local drug delivery strategies. We present a coupled porohyperelastic mass transport (PHEXPT) finite element approach developed using a commercially available ABAQUS finite element software. The PHEXPT transient simulations are based on sequential solution of the porohyperelastic (PHE) and mass transport (XPT) problems where an Eulerian PHE FEM is coupled to a Lagrangian XPT FEM using a custom-written FORTRAN program. The PHEXPT theoretical background is derived in the context of porous media transport theory and extended to ABAQUS finite element formulations. The essential assumptions needed in order to use ABAQUS are clearly identified in the derivation. Representative benchmark finite element simulations are provided along with analytical solutions (when appropriate). These simulations demonstrate the differences in transient and steady state responses including finite deformations, total stress, fluid pressure, relative fluid, and mobile species flux. A detailed description of important model considerations (e.g., material property functions and jump discontinuities at material interfaces) is also presented in the context of finite deformations. The ABAQUS-based PHEXPT approach enables the use of the available ABAQUS capabilities (interactive FEM mesh generation, finite element libraries, nonlinear material laws, pre- and postprocessing, etc.). PHEXPT FEMs can be used to simulate the transport of a relatively large neutral species (negligible osmotic fluid flux) in highly deformable hydrated soft tissues and tissue-engineered materials.

Central to the investigation of the biomechanics of traumatic brain injury (TBI) and the assessment of injury risk from head impact are finite element (FE) models of the human brain. However, many existing FE human brain models have been developed with simplified representations of the parenchyma, which may limit their applicability as an injury prediction tool. Recent advances in neuroimaging techniques and brain biomechanics provide new and necessary experimental data that can improve the biofidelity of FE brain models. In this study, the CAB-20MSym template model was developed, calibrated, and extensively verified. To implement material heterogeneity, a magnetic resonance elastography (MRE) template image was leveraged to define the relative stiffness gradient of the brain model. A multi-stage inverse FE (iFE) approach was used to calibrate the material parameters that defined the underlying non-linear deviatoric response by minimizing the error between model-predicted brain displacements and experimental displacement data. This process involved calibrating the infinitesimal shear modulus of the material using low-severity, low-deformation impact cases and the material non-linearity using high-severity, high-deformation cases from a dataset of in situ brain displacements obtained from cadaveric specimens. To minimize the geometric discrepancy between the FE models used in the iFE calibration and the cadaveric specimens from which the experimental data were obtained, subject-specific models of these cadaveric brain specimens were developed and used in the calibration process. Finally, the calibrated material parameters were extensively verified using independent brain displacement data from 33 rotational head impacts, spanning multiple loading directions (sagittal, coronal, axial), magnitudes (20–40 rad/s), durations (30–60 ms), and severity. Overall, the heterogeneous CAB-20MSym template model demonstrated good biofidelity with a mean overall CORA score of 0.63 ± 0.06 when compared to in situ brain displacement data. Strains predicted by the calibrated model under non-injurious rotational impacts in human volunteers (N = 6) also demonstrated similar biofidelity compared to in vivo measurements obtained from tagged magnetic resonance imaging studies. In addition to serving as an anatomically accurate model for further investigations of TBI biomechanics, the MRE-based framework for implementing material heterogeneity could serve as a foundation for incorporating subject-specific material properties in future models.

Dynamic delamination modelling using interface elements