• Journal of the Optical Society of Korea
• /
• v.18 no.1
• /
• pp.82-88
• /
• 2014

A binary-coupled dipole approximation (BCDA) is described for designing metal nanoparticles with nonperiodic structures in one, two, and three dimensions. This method can be used to simulate the variation of near- and far-field properties through the interactions of metal nanoparticles. An advantage of this method is in its combination with the binary particle swarm optimization (BPSO) algorithm to find the best array of nanoparticles from all possible arrays. The BPSO algorithm has been used to design an array of plasmonic nanospheres to achieve maximum absorption, scattering, and extinction coefficient spectra. In BPSO, a swarm consists of a matrix with binary entries controlling the presence ('1') or the absence ('0') of nanospheres in the array. This approach is useful in optical applications such as solar cells, biosensors, and plasmonic nanoantennae, and optical cloaking.

• Journal of the Korean Institute of Telematics and Electronics S
• /
• v.35S no.5
• /
• pp.103-110
• /
• 1998

In this paper we propose various rotation block realization method of MIMO LBR system. For this, we first discuss how to realize a given LBR system in the form of LBR two-pair cascades. Then we examine several systematic extraction procedures which enable to realize each LBR two-pair in the form of rotation block. Then we consider which realization methods we select for each two-pair block in view of the dynamci range of the coefficients of the rotation blocks. Finally, we demonstrate the LBR realization method using a detailed numerical example that converts a 4th order IIR filter function into a 4*4 transform domain transfer matrix, and then realizes it in a rotation-block cascade.

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

The stochastic analysis of semi-infinite domain is presented using the weighted integral method, which is improved to include the higher order terms in expanding the displacement vector. To improve the weighted integral method, the Lagrangian remainder is taken into account in the expansion of the status variable with respect to the mean value of the random variables. In the resulting formulae only the 'proportionality coefficients' are introduced in the resulting equation, therefore no additional computation time and memory requirement is needed. The equations are applied in analyzing the semi-infinite domain. The results obtained by the improved weighted integral method are reasonable and are in good agreement with those of the Monte Carlo simulation. To model the semi-infinite domain, the Bettess's infinite element is adopted, where the theoretical decomposition of the strain-displacement matrix to calculate the deviatoric stiffness of the semi-infinite domains is introduced. The calculated value of mean and the covariance of the displacement are revealed to be larger than those given by the finite domain assumptions which is thought to be rational and should be considered in the design of structures on semi-infinite domains.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.683-691
• /
• 2010

In this paper, using an adjoint variable method, we develop a design sensitivity analysis(DSA) method applicable to heat conduction problems in steady state. Also, a topology design optimization method is developed using the developed DSA method. Design sensitivity expressions with respect to the thermal conductivity are derived. Since the already factorized system matrix is utilized to obtain the adjoint solution, the cost for the sensitivity computation is trivial. For the topology design optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of structures and allowable material volume respectively. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with finite difference ones, requiring less than 0.25% of CPU time for the finite differencing. Also, the topology optimization yields physical meaningful results.

• Wind and Structures
• /
• v.3 no.4
• /
• pp.291-302
• /
• 2000

A new finite element technique to solve the problem of wind and structure interactions is presented. Conventionally, wind analysis is performed on the Eulerian description in which the finite element mesh would not move in accordance with the wind flow. However, it is not the case in wind-structure interaction problems because nodes attached to the surface of structure should move with the displacement of structure. The arbitrary Lagrangian-Eulerian (ALE) method treats the mesh and flow independently, and allow the mesh to move. In this study, the analysis domain is divided into regions of the structure, air around the structure and the interface of two regions. To satisfy the compatibility and equilibrium conditions between separated regions and to carry out the efficient analysis, the rigid link is used. Also the equation of wind and that of structure are arranged in a single matrix equation.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.809-827
• /
• 2015

In this paper, seismic energy response of inelastic steel structures under earthquake excitations is investigated. For this purpose, a numerical procedure based on nonlinear dynamic analysis is developed by considering material, geometric and connection nonlinearities. Material nonlinearity is modeled by the inversion of Ramberg-Osgood equation. Nonlinearity caused by the interaction between the axial force and bending moment is also defined considering stability functions, while the geometric nonlinearity caused by axial forces is described using geometric stiffness matrix. Cyclic behaviour of steel connections is taken into account by employing independent hardening model. Dynamic equation of motion is solved by Newmark's constant acceleration method in the time history domain. Energy response analysis of space frames is performed by using this proposed numerical method. Finally, for the first time, the distribution of the different energy types versus time at the duration of the earthquake ground motion is obtained where in addition error analysis for the numerical solutions is carried out and plotted depending on the relative error calculated as a function of energy balance versus time.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.725-740
• /
• 2015

This study presents critical buckling load optimization of the axially graded layered uniform columns. In the first place, characteristic equations for the critical buckling loads for all boundary conditions are obtained using the transfer matrix method. Then, for each case, square of this equation is taken as a fitness function together with constraints. Due to explicitly unavailable objective function for the critical buckling loads as a function of segment length and volume fraction of the materials, especially for the column structures with higher segment numbers, initially, prescribed value is assumed for it and then the design variables satisfying constraints are searched using Differential Evolution (DE) optimization method coupled with eigen-value routine. For constraint handling, Exterior Penalty Function formulation is adapted to the optimization cycle. Different boundary conditions are considered. The results reveal that maximum increments in the critical buckling loads are attained about 20% for cantilevered and pinned-pinned end conditions and 18% for clamped-clamped case. Finally, the strongest column structure configurations will be determined. The scientific and statistical results confirmed efficiency, reliability and robustness of the Differential Evolution optimization method and it can be used in the similar problems which especially include transcendental functions.

• Wind and Structures
• /
• v.14 no.4
• /
• pp.359-381
• /
• 2011

This paper deals with the flutter instability problem of flexible bridge decks in the framework of bimodal-coupled aeroelastic bridge system analysis. Based on the analysis of coefficients of the polynomials deduced from the singularity conditions of an integral wind-structure impedance matrix, a set of simplified formulations for calculating the critical wind velocity and coupled frequency are presented. Several case studies are discussed and comparisons with available approximated approaches are made and presented, along with a conventional complex eigenvalue analysis and numerical results. From the results, it is found that the formulas that are presented in this study are applicable to a variety of bridge cross sections that are not only prone to coupled-mode but also to single-mode-dominated flutter.

• Structural Engineering and Mechanics
• /
• v.36 no.5
• /
• pp.529-544
• /
• 2010

Nonlinear equations of structures are generally solved numerically by the iterative solution of linear equations. However, this iterative procedure diverges when the tangent stiffness is ill-conditioned which occurs near limit points. In other words, a major challenge with simple iterative methods is failure caused by a singular or near singular Jacobian matrix. In this paper, using the Newton-Raphson algorithm based on Davidenko's equations, the iterations can traverse the limit point without difficulty. It is argued that the propose algorithm may be both more computationally efficient and more robust compared to the other algorithm when tracing path through severe nonlinearities such as those associated with structural collapse. Two frames are analyzed using the proposed algorithm and the results are compared with the previous methods. The ability of the proposed method, particularly for tracing the limit points, is demonstrated by those numerical examples.

• Steel and Composite Structures
• /
• v.12 no.4
• /
• pp.303-324
• /
• 2012

This paper presents the mode-shape analysis of the cross-ply laminated composite cylindrical shallow shells. First, the kinematic relations of strains and deformation are given. Then, using Hamilton's principle, governing differential equations are developed for a general curved shell. Finally, the stress-strain relation for the laminated, cross-ply composite shells are obtained. By using some simplifications and assuming Fourier series as a displacement field, the governed differential equations are solved by the matrix algebra for shallow shells. Employing the computer algebra system called MATHEMATICA; a computer program has been prepared for the solution. The results obtained by this solution are compared with the results obtained by (ANSYS and SAP2000) programs, in order to verify the accuracy and reliability of the solution presented.