• Journal of the Korean Geotechnical Society
• /
• v.16 no.4
• /
• pp.191-199
• /
• 2000

In this research, the finite element analysis of piezocone penetration and dissipation tests has been conducted using the anisotropic elastoplastic-viscoplastic bounding surface model, virtual work equation, and theory of mixtures formulated in the Up[dated Lagrangian reference frame for the large deformation and finite strain nature of piezocone penetration. The formulated equations have been implemented into a finite element program. The cone resistance, excess pore water pressure, and dissipation of excess pore water pressure from the finite element analysis have been compared and investigated. An effective simulation could be performed with the use of the anisotropic and viscous soil model. The finite element formulations and the results are described in part 'I' and part 'II' respectively.

• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.125-143
• /
• 2014

This paper provides a new technique for solving the static analysis of arbitrarily shaped composite plates by using Strong Formulation Finite Element Method (SFEM). Several papers in literature by the authors have presented the proposed technique as an extension of the classic Generalized Differential Quadrature (GDQ) procedure. The present methodology joins the high accuracy of the strong formulation with the versatility of the well-known Finite Element Method (FEM). The continuity conditions among the elements is carried out by the compatibility or continuity conditions. The mapping technique is used to transform both the governing differential equations and the compatibility conditions between two adjacent sub-domains into the regular master element in the computational space. The numerical implementation of the global algebraic system obtained by the technique at issue is easy and straightforward. The main novelty of this paper is the application of the stress and strain recovery once the displacement parameters are evaluated. Computer investigations concerning a large number of composite plates have been carried out. SFEM results are compared with those presented in literature and a perfect agreement is observed.

• Computers and Concrete
• /
• v.6 no.4
• /
• pp.319-341
• /
• 2009

A new formulation based on lumped plasticity and inelastic hinges is presented in this paper for nonlinear analysis of Reinforced Concrete (RC) frame structures. Inelastic hinge behaviour is described using the principles of Continuum Damage Mechanics (CDM). Member formulation contains provisions to model stiffness degradation due to cracking of concrete and yielding of reinforcing steel. Depending on its nature, cracking is classified as concentrated or distributed. Concentrated cracking is accounted through a damage variable and its growth is defined based on strain energy principles. Presence of distributed flexural cracks in a member is taken care of by modelling it as non-prismatic. Plasticity theory supported by effective stress concept of CDM is applied to describe the post-yield response. Nonlinear quasi-static analysis is carried out on a RC column and a wide two-storey RC frame to verify the formulation. The column is subjected to constant axial load and monotonic lateral load while the frame is subjected to only lateral load. Computed results are compared with those due to experiments or other numerical methods to validate the performance of the formulation and also to highlight the contribution of distributed cracking on global response.

• International Journal of Quality Innovation
• /
• v.3 no.2
• /
• pp.174-182
• /
• 2002

System identification and controller formulation are essential in dynamic process control. In system identification, data for system identification are obtained, and then they are analyzed so that the system model of the process is built, identified, and diagnosed. In controller formulation, the control equation is derived based on the result of the system identification. There has been much theoretical research on system identification and controller formulation. These theories are very useful when they are appropriately applied. To our regret, however, these theories are not always effectively applied in practice because the engineers and the operators who manage the process often do not have the necessary understanding of required time series analysis methods. On the other hand, because of widespread use of statistical packages, system identification such as estimating ARMA models can be done with little understanding of time series analysis methods. Therefore, it might be said that the most theoretically difficult part in practice is the controller formulation. In this paper, lists of control equations are proposed as a useful tool for practitioners to use. The tool supports bridging the gap between theory and practice in dynamic process control. Also, for some models, the generalized control equations are obtained.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.1
• /
• pp.267-285
• /
• 2014

An upstream scheme based on the pseudostress-velocity mixed formulation is studied to solve convection-dominated Oseen equations. Lagrange multipliers are introduced to treat the trace-free constraint and the lowest order Raviart-Thomas finite element space on rectangular mesh is used. Error analysis for several quantities of interest is given. Particularly, first-order convergence in $L^2$ norm for the velocity is proved. Finally, numerical experiments for various cases are presented to show the efficiency of this method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.04a
• /
• pp.17-27
• /
• 1995

For the post-buckling and elasto-plastic analysis of shell structures, the total Lagrangian formulation is presented based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors in the iteration process and evaluating the total Green-Lagrange stain corresponding U total displacements. In the calculation of the stiffness matrix, the element formulation takes into account the effect of finite rotation increments by retaining second order rotation terms in the incremental displacement field. The selective or reduced integration scheme using the heterosis element is applied in order to overcome both shear locking phenomena and the zero energy mode. The load/displacement incremental scheme is adopted for geometric non-linear F .E. analysis. Based on such methodology, the computer program is developed and numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with references's results.

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.33-42
• /
• 2017

We present in this paper a finite element formulation for nonlinear torsional analysis of 3D beams with arbitrary composite cross-sections. Since the proposed formulation employs a continuum mechanics based beam element with kinematics enriched by the extended St. Venant solutions, it can precisely account higher order warping effect and its 3D couplings. We propose a numerical procedure to calculate the extended St. Venant equation and the twisting center of an arbitrary composite cross-section simultaneously. The accuracy and efficiency of the proposed formulation are thoroughly investigated through representative numerical examples.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.2
• /
• pp.199-209
• /
• 2012

An advanced aeroelastic formulation for flutter analyses is presented in this paper. Refined 1D structural models were coupled with the doublet lattice method, and the g-method was used for flutter analyses. Structural models were developed in the framework of the Carrera Unified Formulation (CUF). Higher-order 1D structural models were obtained by using Taylor-like expansions of the cross-section displacement field of the structure. The order (N) of the expansion was considered as a free parameter since it can be arbitrarily chosen as an input of the analysis. Convergence studies on the order of the structural model can be straightforwardly conducted in order to establish the proper 1D structural model for a given problem. Flutter analyses were conducted on several wing configurations and the results were compared to those from literature. Results show the enhanced capabilities of CUF 1D in dealing with the flutter analysis of typical wing structures with high accuracy and low computational costs.

• Journal of Mechanical Science and Technology
• /
• v.18 no.2
• /
• pp.193-202
• /
• 2004

This paper presents a computer method for the dynamic analysis of a system of rigid bodies in plane motion. The formulation rests upon the idea of replacing a rigid body by a dynamically equivalent constrained system of particles. Newton's second law is applied to study the motion of the resulting system of particles without introducing any rotational coordinates. A velocity transformation is used to transform the equations of motion to a reduced set. For an open-chain, this process automatically eliminates all of the non-working constraint forces and leads to an efficient integration of the equations of motion. For a closed-chain, suitable joints should be cut and few cut-joints constraint equations should be included. An example of a closed-chain is used to demonstrate the generality and efficiency of the proposed method.