• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Advances in aircraft and spacecraft science
• /
• v.2 no.1
• /
• pp.1-16
• /
• 2015

In this paper, geometrical and mechanical approaches are proposed for the simulation of the draping of woven fabric onto complex parts. The geometrical discrete approach allows to define the ply shapes and fibres orientation in order to optimize the composite structural properties and the continuum meso-structural mechanical approach allows to take into account the mechanical properties of fibres and resin and the various dominating mode of deformation of woven fabrics during the forming process. Some numerical simulations of forming process are proposed and compared with the experimental results in order to demonstrate the efficiency of our approaches.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.4
• /
• pp.32-36
• /
• 2002

In this paper, a numerical methodology for health monitoring of weldment was proposed using finite element method coupled with continuum damage mechanics. The welding-induced residual stress distribution of T-joint weldment was calculated using a commercial finite element package SYSWELD+. The distribution of latent damage was evaluated from the stress and strain components taken as the output of a finite element calculation. Numerical examples were given to demonstrate the usefulness of this so-called "post-processing approach" in the case of welding-induced damage assessment.

• Transactions of Materials Processing
• /
• v.17 no.3
• /
• pp.161-169
• /
• 2008

In the present work, the two-back stress model is proposed and continuum damage mechanics (CDM) is incorporated into the plastic constitutive relation in order to describe the plastic deformation localization and the damage evolution in a deforming continuum body. Coupling between damage mechanics and isothermal rate independent plasticity is performed using the kinematic hardening rule, which in turn is formulated by combining the nonlinear Armstrong-Frederick rule and the Phillips rule. The numerical analyses are carried out within h deformation theory. It is noted that the damage evolution within a work piece accelerates the plastic deformation localization such that the material with lower hardening exponent results in a rapid shear band formation. Moreover, the results from the numerical analysis reflected closely with the micro-structures around the fractured regime. The effects of the various hardening parameters on deformation localization are also investigated. As the nonlinear strain rate description in the back stress evolution becomes dominant, the strain localization becomes intensified as well as the damage evolution.

• Structural Engineering and Mechanics
• /
• v.33 no.6
• /
• pp.657-674
• /
• 2009

In this paper a 3-D continuum damage mechanics formulation for composite laminates and its implementation into a finite element model that is based on the layer-wise laminate plate theory are described. In the damage formulation, each composite ply is treated as a homogeneous orthotropic material exhibiting orthotropic damage in the form of distributed microscopic cracks that are normal to the three principal material directions. The progressive damage of different angle ply composite laminates under quasi-static loading that exhibit the free edge effects are investigated. The effects of various numerical modeling parameters on the progressive damage response are investigated. It will be shown that the dominant damage mechanism in the lay-ups of [+30/-30]s and [+45/-45]s is matrix cracking. However, the lay-up of [+15/-15] may be delaminated in the vicinity of the edges and at $+{\theta}/-{\theta}$ layers interfaces.

• Geomechanics and Engineering
• /
• v.6 no.1
• /
• pp.17-31
• /
• 2014

This paper investigates the development of failure surfaces induced by an embankment on soft marine clay deposits and the characteristics of such surfaces through numerical simulations and its comparative study with monitoring results. It is well known that the factor of safety of embankment slopes is closely related to the vertical loading, including the height of the embankment. That is, an increase in the embankment height reduces the factor of safety. However, few studies have examined the relationship between the lateral movement of soft soil beneath the embankment and the factor of safety. In addition, no study has investigated the distribution of the pore pressure coefficient B value along the failure surface. This paper conducts a continuum analysis using finite difference methods to characterize the development of failure surfaces during embankment construction on soft marine clay deposits. The results of the continuum analysis for failure surfaces, stress, displacement, and the factor of safety can be used for the management of embankment construction. In failure mechanism, it has been validated that a large shear displacement causes change of stress and pore pressure along the failure surface. In addition, the pore pressure coefficient B value decreases along the failure surface as the embankment height increases. This means that the rate of change in stress is higher than that in pore pressure.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.511-521
• /
• 2002

A simple numerical modelling technique is proposed for estimating the shear stress distribution in beams of framed tube structures with multiple internal tubes. The structures arc analysed using a continuum approach in which each tube is individually modelled by a tube beam that accounts for the flexural and shear deformations, as well as the shear lag effects. The numerical analysis of shear stress is based on the mathematical analogy in conjunction with the elastic theory By simplifying assumptions regarding the form of strain distributions in external and internal tubes, the shear stress distributions are expressed in terms of a series of lineal functions of the second moments of area of the structures and the corresponding geometric and material properties, as well as the applied loads. Previous studies for axial stresses and shear lag phenomenon are further developed lot the numerical analysis of shear stresses in the tubes. The simplicity and accuracy of the proposed method are demonstrated through the solutions of throe numerical examples.

• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.519-540
• /
• 2007

In this paper, the energy-based plastic-damage model previously proposed by the authors [International Journal of Solids and Structures, 43(3-4): 583-612] is first simplified with an empirically defined evolution law for the irreversible strains, and then it is extended to its rate-dependent version to account for the strain rate effect. Regarding the energy dissipation by the motion of the structure under dynamic loadings, within the framework of continuum damage mechanics a new damping model is proposed and incorporated into the developed rate-dependent plastic-damage mode, leading to a unified constitutive model which is capable of directly considering the damping on the material scale. Pertinent computational aspects concerning the numerical implementation and the algorithmic consistent modulus for the unified model are also discussed in details, through which the dynamic nonlinear analysis of damping structures can be coped with by the same procedures as those without damping. The proposed unified plastic-damage model is verfied by the simulations of concrete specimens under different quasistatic and high rate straining loading conditions, and is then applied to the Koyna dam under earthquake motions. The numerical predictions agree fairly well with the results obtained from experimental tests and/or reported by other investigators, demonstrating its capability for reproducing most of the typical nonlinear performances of concrete under quasi-static and dynamic loading conditions.

• International Journal of Concrete Structures and Materials
• /
• v.2 no.2
• /
• pp.123-136
• /
• 2008

A new plastic-damage constitutive model applicable to lightweight concrete (LWC) and normal weight concrete (NWC) is proposed in this paper based on both continuum damage mechanics and plasticity theories. Two damage variables are used to represent tensile and compressive damage independently. The effective stress is computed in the Drucker-Prager multi-surface plasticity framework. The stress is then computed by multiplication of the damaged part and the effective part. The proposed model is coded as a user material subroutine and incorporated in a finite element analysis software. The constitutive integration algorithm is implemented by adopting the operator split involving elastic predictor, plastic corrector and damage corrector. The numerical study shows that the algorithm is efficient and robust in the finite element analysis. Experimental investigation is conducted to verify the proposed model involving both static and dynamic tests. The very good agreement between the numerical results and experimental results demonstrates the capability of the proposed model to capture the behaviors of LWC and NWC structures for static and impact loading.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.90-95
• /
• 2001

A continuum-based shape design sensitivity analysis (DSA) method is presented for the acoustic radiation from thin body. The normal derivative integral formulation is employed as an analysis formulation and differentiated directly by using material derivative to get the acoustic shape design sensitivity. In the acoustic sensitivity formulation, derivative coefficients of the structural normal velocities on the surface are required as the input. Thus, the shape design sensitivities of structural velocities on the surface with respect to the shape change are also calculated with continuum approach. A simple disk is considered as a numerical example to validate the accuracy and efficiency of the analytical shape design sensitivity equations derived in this research. This research should be very helpful to design an application involving thin body and to change its acoustic characteristics.