• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.185-193
• /
• 1997

The order of the stress singularity that occurs at the termination of an interface between materials exhibiting bilinear stress-strain response under plane strain conditions has been calculated, The governing equation of elasticity together with traction-free boundary condition and interface continuity condition defines a two-point boundary value problem. The stress components near the free edge are assumed to be proportional to r/sup s-1/, with solutions existing only for certain values of s. Finding these values entails the solution of an eigenvalue problem. Because it has been impossible to integrate the differential equations analytically, the integration has been performed numerically with a shooting method coupled with a Newton improvement scheme.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.2 s.72
• /
• pp.213-219
• /
• 2006

Topology optimization problem requires numerous repeated evaluations of objective function and design sensitivity for elements within design domain with various density distributions. The recently proposed two-level condensation scheme(TLCS) is very promising for the construction of reduced system and for an accurate and efficient analysis concerned about eigenvalue and dynamic problems. We used the two-level dynamic condensation scheme for the analysis and sensitivity computation part in the structural topology optimization problem. The results of the topology optimization for the reduced system show the TLCS provides high accuracy and computation efficiency compared to the full scale system within engineering accuracy.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.310-323
• /
• 1990

A three-noded, with three degree-of-freedom at each node, in-plane curved beam element is formulated and employed in eigen-analysis of constant curvature beam. The conventional quadratic shape functions used in a three noded C .deg. type curved beam element produce such an undesirable large stiffness that a significant error is introduced in displacements and stresses. These phenomena are called 'Stiffness Locking Phenomena', which result from spurious strain energy due to inappropriate assumptions on independent isoparametric quadratic interpolation functions. Stiffness locking phenomena can be alleviated by using modified interpolation functions which get rid of spurious constraints of conventional interpolation functions. Eigenvalues and their modes as well as displacements and stresses may be locked because they are related to stiffness. Using modified curved beam element in eigenvalue problem of cantilever and arch, the property and performance of modified curved beam element are examined by numerical experimentations. In these eigen-analyses, mass matrices are calculated by using both modified and unmodified curved beam element, are compared with theoretical solutions. These comparisons show that the performance of the modified curved beam element is better than that of the unmodified curved beam element.

• Structural Engineering and Mechanics
• /
• v.64 no.3
• /
• pp.381-390
• /
• 2017

Present paper deals with the large amplitude flexural vibration of carbon nanotube reinforced composite (CNTRC) plates. Distribution of CNTs as reinforcements may be uniform or functionally graded (FG). The equivalent material properties of the composite media are obtained according to a refined rule of mixtures which contains efficiency parameters. To account for the large deformations, von $K{\acute{a}}rm{\acute{a}}n$ type of geometrical nonlinearity is included into the formulation. The matrix representation of the governing equations is obtained according to the Ritz method where the basic shape functions are written in terms of the Chebyshev polynomials. Time dependency of the problem is eliminated by means of the Galerkin method and the resulting nonlinear eigenvalue problem is solved employing a direct displacement control approach. Results are obtained for completely clamped and completely simply supported plates. Results are first validated for the especial cases of FG-CNTRC and cross-ply laminated plates. Afterwards, parametric studies are given for FG-CNTRC plates with different boundary conditions. It is shown that, nonlinear frequencies are highly dependent to the volume fraction and dispersion profiles of CNTs. Furthermore, mode redistribution is observed in both simply supported and clamped FG-CNTRC plates.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.651-659
• /
• 2002

Component Mode Synthesis (CMS) is a dynamic substructuring technique to get an approximate eigensolutions of large degree-of-freedom structures divisible into several components. But, In practice. most of large structures are modeled by different teams of engineers. and their respective finite element models often require different mesh resolutions. As a result, the finite element substructure models can be non-conforming and/or incompatible. In this work, A hybrid version of component mode synthesis using a localized lagrange multiplier to treat the non-conforming mesh problem was derived. Evolution Strategies (ESs) is a stochastic numerical optimization technique and has shown a robust performance for solving deterministic problems. An ESs conducts its search by processing a population of solutions for an optimization problem based on principles from natural evolution. An optimization example for raising the first natural frequency of a plate structure using beam stiffeners was presented using hybrid component mode synthesis and robust evolution strategies (RES) optimization technique. In the example. the design variables are the positions and lengths of beam stiffeners.

• Nuclear Engineering and Technology
• /
• v.27 no.3
• /
• pp.374-384
• /
• 1995

The mathematical adjoint solution of the Analytic Function Expansion (AFEN) method is found by solving the transposed matrix equation of AFEN nodal equation with only minor modification to the forward solution code AFEN. The perturbation calculations are then performed to estimate the change of reactivity by using the mathematical adjoint The adjoint calculational scheme in this study does not require the knowledge of the physical adjoint or the eigenvalue of the forward equation. Using the adjoint solutions, the exact and first-order perturbation calculations are peformed for the well-known benchmark problems (i.e., IAEA-2D benchmark problem and EPRI-9R benchmark problem). The results show that the mathematical adjoint flux calculated in the code is the correct adjoint solution of the AFEN method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.1
• /
• pp.66-76
• /
• 2003

The arbitrary V-notched crack problem is considered. The general expressions for the stress components on this problem are obtained as explicit series forms composed of independent unknown coefficients which are denoted by coefficients of eigenvector. For this results eigenvalue equation is performed first through introducing complex stress functions and applying the traction free boundary conditions. Next solving this equation, eigenvalues and corresponding eigenvectors are obtained respectively, and finally inserting these results into stress components, the general equations are obtained. These results are also shown to be applicable to the symmetric V-notched crack or straight crack. It can be shown that this solutions are composed of the linear combination of Mode I and Mode II solutions which are obtained from different characteristic equations, respectively. Through performing asymptotic analysis for stresses, the stress intensity factor is given as a closed form equipped with the unknown coefficients of eigenvector. In order to calculate the unknown coefficients. based on these general explicit equations, numerical programming using the overdetermined boundary collocation method which is algorithmed originally by Carpenter is also worked out. As this programming requires the input data, the commercial FE analysis for stresses is performed. From this study, for some V-notched problems, unknown coefficients can be calculated numerically and also fracture parameters are determined.

• Kyungpook Mathematical Journal
• /
• v.47 no.4
• /
• pp.455-471
• /
• 2007

The magnetohydrodynamic axisymmetric instability of a streaming gas jet surrounded by bounded fluid with non-uniform field has been developed. The problem is formulated, solved and the boundary conditions are applied across the interfaces. The eigenvalue relation is derived and discussed analytically and the results are confirmed numerically. Some reported works are recovered as limiting cases from the present general results. The streaming has a destabilizing effect for all short and long wavelengths. The capillary force is stabilizing for short wavelengths but it is destabilizing for long wavelengths. The axial magnetic fields interior the gas and fluid media are stabilizing. The transverse field is destabilizing for all wavelengths. The radii ratio of the gas and fluid cylinders plays an important role for stabilizing the model and made it more realistic one than the full liquid jet or/and the ordinary hollow jet. The numerical analysis clarify the stable and unstable domains based on different values of the various parameters of the problem.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.2 s.257
• /
• pp.221-230
• /
• 2007

An iterated improved reduced system (IIRS) procedure combined with sub-structuring scheme for nonclassically damped structural systems is presented. For dynamic analysis of such systems, complex eigenproperties are required to incorporate properly the nonclassical damping effect. In complex structural systems, the equations of motion are written in the state space from. Thus, the number of degrees of freedom of the new equations of motion and the size of the associated eigenvalue problem required to obtain the complex eigenvalues and eigenvectors are doubled. Iterated IRS method is an efficient reduction technique because the eigenproperties obtained in each iteration step improve the condensation matrix in the next iteration step. However, although this reduction technique reduces the size of problem drastically, it is not efficient to apply this technique to a single domain finite element model with degrees of freedom over several thousands. Therefore, for a practical application of the reduction method, accompanying sub-structuring scheme is necessary. In the present study, iterated IRS method combined with sub-structuring scheme for nonclssically damped structures is developed. Numerical examples demonstrate the convergence and the efficiency of a newly developed scheme.

• Wind and Structures
• /
• v.13 no.3
• /
• pp.235-256
• /
• 2010

The aeroelastic stability of bridge decks equipped with multiple tuned mass dampers is studied. The problem is attacked in the time domain, by representing self-excited loads with the aid of aerodynamic indicial functions approximated by truncated series of exponential filters. This approach allows to reduce the aeroelastic stability analysis in the form of a direct eigenvalue problem, by introducing an additional state variable for each exponential term adopted in the approximation of indicial functions. A general probabilistic framework for the optimal robust design of multiple tuned mass dampers is proposed, in which all possible sources of uncertainties can be accounted for. For the purposes of this study, the method is also simplified in a form which requires a lower computational effort and it is then applied to a general case study in order to analyze the control effectiveness of regular and irregular multiple tuned mass dampers. A special care is devoted to mistuning effects caused by random variations of the target frequency. Regular multiple tuned mass dampers are seen to improve both control effectiveness and robustness with respect to single tuned mass dampers. However, those devices exhibit an asymmetric behavior with respect to frequency mistuning, which may weaken their feasibility for technical applications. In order to overcome this drawback, an irregular multiple tuned mass damper is conceived which is based on unequal mass distribution. The optimal design of this device is finally pursued via a full domain search, which evidences a remarkable robustness against frequency mistuning, in the sense of the simplified design approach.