• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.32 no.4
• /
• pp.432-446
• /
• 1996

A stochastic Hamilton variational principle(SHVP) is formulated for dynamic problems of linear continuum. The SHVP allows incorporation of probabilistic distributions into the finite element analysis. The formulation is simplified by transformation of correlated random variables to a set of uncorrelated random variables through a standard eigenproblem. A procedure based on the Fourier analysis and synthesis is presented for eliminating secularities from the perturbation approach. In addition to, a method to analyse stochastic design sensitivity for structural dynamics is present. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element codes. An alternative form of the constraint functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The algorithms developed can readily be adapted to existing deterministic finite element codes. The numerical results for stochastic analysis by proceeding approach of cantilever, 2D-frame and 3D-frame illustrates in this paper.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.313-323
• /
• 2018

A numerical study is performed to investigate the impacts of thermal loading on the vibration and buckling of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) conical shells. Thermo-mechanical properties of constituents are considered to be temperature-dependent. Considering the shear deformation theory, the energy functional is derived, and applying the variational differential quadrature (VDQ) method, the mass and stiffness matrices are obtained. The shear correction factors are accurately calculated by matching the shear strain energy obtained from an exact three-dimensional distribution of the transverse shear stresses and shear strain energy related to the first-order shear deformation theory. Numerical results reveal that considering temperature-dependent material properties plays an important role in predicting the thermally induced vibration of FG-CNTRC conical shells, and neglecting this effect leads to considerable overestimation of the stiffness of the structure.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.5
• /
• pp.850-858
• /
• 1997

A new a posteriori error estimator is introduced and applied to variational inequalities occurring in problems of flow through porous media. In order to construct element-wise a posteriori error estimates the global error is localized by a special mixed formulation in which continuity conditions at interfaces are treated as constraints. This approach leads to error indicators which provide rigorous upper bounds of the element errors. A discussion of a compatibility condition for the well-posedness of the local error analysis problem is given. Two numerical examples are solved to check the compatibility of the local problems and convergence of the effectivity index both in a local and a global sense with respect to local refinements.

• East Asian mathematical journal
• /
• v.14 no.1
• /
• pp.63-76
• /
• 1998

we consider the hp-version to solve non-constant coefficients elliptic equations $-div(a{\nabla}u)=f$ with Dirichlet boundary conditions on a bounded polygonal domain $\Omega$ in $R^2$. In [6], M. Suri obtained an optimal error-estimate for the hp-version: ${\parallel}u-u^h_p{\parallel}_{1,\Omega}{\leq}Cp^{(\sigma-1)}h^{min(p,\sigma-1)}{\parallel}u{\parallel}_{\sigma,\Omega}$. This optimal result follows under the assumption that all integrations are performed exactly. In practice, the integrals are seldom computed exactly. The numerical quadrature rule scheme is needed to compute the integrals in the variational formulation of the discrete problem. In this paper we consider a family $G_p=\{I_m\}$ of numerical quadrature rules satisfying certain properties, which can be used for calculating the integrals. Under the numerical quadrature rules we will give the variational form of our non-constant coefficients elliptic problem and derive an error estimate of ${\parallel}u-\tilde{u}^h_p{\parallel}_{1,\Omega}$.

• Composites Research
• /
• v.30 no.6
• /
• pp.343-349
• /
• 2017

A multifield variational formulation is developed for the finite element (FE) cross-sectional analysis of composite beams. The cross-sectional warping displacements and sectional stresses are considered to be the primary variables through the application of Reissner's partially mixed principle. The warping displacements are modeled using generic FE shape functions with nonlinear distribution over the beam section. A generalized Timoshenko level stiffness matrix is derived which incorporates the effects of elastic couplings, transverse shear, and Poisson's deformations. The accuracy of the present analysis is validated for the stiffness constants and elastostatic responses of composite box beams which correlate well with the experimental data and other state-of-the-art approaches.

• Steel and Composite Structures
• /
• v.27 no.6
• /
• pp.789-801
• /
• 2018

A weak-form formulation of finite annular prism methods (FAPM) based on Reissner's mixed variational theorem (RMVT), is developed for the quasi three-dimensional (3D) static analysis of two-directional functionally graded (FG) circular plates with various boundary conditions and under mechanical loads. The material properties of the circular plate are assumed to obey either a two-directional power-law distribution of the volume fractions of the constituents through the radial-thickness surface or an exponential function distribution varying doubly exponentially through it. These FAPM solutions of the loaded FG circular plates with both simply-supported and clamped edges are in excellent agreement with the solutions obtained using the 3D analytical approach and two-dimensional advanced plate theories available in the literature.

• Structural Engineering and Mechanics
• /
• v.40 no.6
• /
• pp.813-827
• /
• 2011

This paper provides a numerical solution for a finite internally cracked plate using hybrid crack element method (HCE). In the formulation, an inclined crack is placed in any place of a rectangular element and the complex variable method is used. The complex potentials are expressed in a series form, and several undetermined coefficients are involved. The complex potentials for the cracked rectangle are first suggested in this paper. Based on a variational principle, the element stiffness matrix can be evaluated. The next steps are same as in the usual finite element method. Several numerical examples with computed stress intensity factor and T-stress are presented.

• Journal of Power System Engineering
• /
• v.10 no.4
• /
• pp.83-89
• /
• 2006

Curved beams are more efficient in transfer of loads than straight beams because the transfer is effected by bending, shear and membrane action. The finite element method is a versatile method for solving structural mechanics problems and curved beam problems have been solved using this method by many author. In this study, a new three-noded hybrid-mixed curved beam element is proposed to investigate the in-plane flexural vibration behavior of arches depending on the curvature, aspect ratio and boundary conditions, etc. The proposed element including the effect of shear deformation is based on the Hellinger-Reissner variational principle, and employs the quadratic displacement functions and consistent linear stress functions. The stress parameters are then eliminated from the stationary condition of the variational principle so that the standard stiffness equations are obtained. Several numerical examples confirm the accuracy of the proposed finite element and also show the dynamic behavior of arches with various shapes.

• Steel and Composite Structures
• /
• v.39 no.3
• /
• pp.291-306
• /
• 2021

Based on Reissner's mixed variational theorem (RMVT), the authors develop a semi-analytical finite element (FE) method for a three-dimensional (3D) bending analysis of nonhomogeneous orthotropic, complete and incomplete toroidal shells subjected to uniformly-distributed loads. In this formulation, the toroidal shell is divided into several finite annular prisms (FAPs) with quadrilateral cross-sections, where trigonometric functions and serendipity polynomials are used to interpolate the circumferential direction and meridian-radial surface variations in the primary field variables of each individual prism, respectively. The material properties of the toroidal shell are considered to be nonhomogeneous orthotropic over the meridianradial surface, such that homogeneous isotropic toroidal shells, laminated cross-ply toroidal shells, and single- and bi-directional functionally graded toroidal shells can be included as special cases in this work. Implementation of the current FAP methods shows that their solutions converge rapidly, and the convergent FAP solutions closely agree with the 3D elasticity solutions available in the literature.

• Earthquakes and Structures
• /
• v.16 no.1
• /
• pp.55-67
• /
• 2019

The finite solutions of deflection and the corresponding in-plane stress values of the graded sandwich shallow shell structure are computed in this research article via a higher-order polynomial shear deformation kinematics. The shell structural equilibrium equation is derived using the variational principle in association with a nine noded isoprametric element (nine degrees of freedom per node). The deflection values are computed via an own customized MATLAB code including the current formulation. The stability of the current finite element solutions including their accuracies have been demonstrated by solving different kind of numerical examples. Additionally, a few numerical experimentations have been conducted to show the influence of different design input parameters (geometrical and material) on the flexural strength of the graded sandwich shell panel including the geometrical configurations.