• Coupled systems mechanics
• /
• 제7권6호
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Structural Engineering and Mechanics
• /
• 제70권5호
• /
• pp.551-562
• /
• 2019

In the present study, a suitable mathematical model considering parabolic transverse shear strains for dynamic analysis of laminated composite skew plates under different types of impulse and spatial loads was presented for the first time. The proposed mathematical model satisfies zero transverse shear strain at the top and bottom of the plate. On the basis of the cubic variation of thickness coordinate in in-plane displacement fields of the present mathematical model, a 2D finite element (FE) model was developed including skew transformations in the mathematical model. No shear correction factor is required in the present formulation and damping effect was also incorporated. This is the first FE implementation considering a cubic variation of thickness coordinate in in-plane displacement fields including skew transformations to solve the forced vibration problem of composite skew plates. The effect of transverse shear and rotary inertia was incorporated in the present model. The Newmark-${\beta}$ scheme was adapted to perform time integration from step to step. The $C^0$ FE formulation was implemented to overcome the problem of $C^1$ continuity associated with the cubic variation of thickness coordinate in in-plane displacement fields. The numerical studies showed that the present 2D FE model predicts the result close to the analytical results. Many new results varying different parameter such as skew angles, boundary conditions, etc. were presented.

• Earthquakes and Structures
• /
• 제22권3호
• /
• pp.245-261
• /
• 2022

Key questions to researchers interested in nonlinear analysis of skeletal structures are whether the distributed plasticity approach - albeit computationally demanding - is more reliable than the concentrated plasticity to adequately capture the extent and severity of the inelastic response, and whether force-based formulation is more efficient than displacement-based formulation without compromising accuracy. The present research focusing on performance-based seismic response of mid-span concrete bridges provides a pilot holistic investigation opting for some hands-on answers. OpenSees software is considered adopting different modeling techniques, viz. distributed plasticity (through either displacement-based or force-based elements) and concentrated plasticity via beam-with-hinges elements. The pros and cons of each are discussed based on nonlinear pushover analysis results, and fragility curves generated for various performance levels relying on incremental dynamic analyses under real earthquake records. Among prime conclusions, distributed plasticity modeling albeit inherently not relying on prior knowledge of plastic hinge length still somewhat depends on such information to ensure accurate results. For instance, displacement-based and force-based approaches secure optimal accuracy when dividing, for the former, the member into sub-elements, and satisfying, for the latter, a distance between any two consecutive integration points, close to the expected plastic hinge length. On the other hand, using beam-with-hinges elements is computationally more efficient relative to the distributed plasticity, yet with acceptable accuracy provided the user has prior reasonable estimate of the anticipated plastic hinge length. Furthermore, when intrusive performance levels (viz. life safety or collapse) are of concern, concentrated plasticity via beam-with-hinges ensures conservative predicted capacity of investigated bridge systems.

• Structural Engineering and Mechanics
• /
• 제52권1호
• /
• pp.115-136
• /
• 2014

In this study, an extended Kantorovich method, employing multi-term displacement functions, is applied to analyze the vibration problem of symmetrically laminated plates with arbitrary boundary conditions. The vibration behaviors of laminated plates are determined based on the variational principle of total energy minimization and the iterative Kantorovich method. The out-of-plane displacement is represented in the form of a series of a sum of products of functions in x and y directions. With a known function in the x or y directions, the formulation for the variation of total potential energy is transformed to a set of governing equations and a set of boundary conditions. The equations and boundary conditions are then numerically solved for the natural frequency and vibration mode shape. The solutions are verified with available solutions from the literature and solutions from the Ritz and finite element analysis. In most cases, the natural frequencies compare very well with the reference solutions. The vibration mode shapes are also very well modeled using the multi-term assumed displacement function in the terms of a power series. With the method used in this study, it is possible to solve the angle-ply plate problem, where the Kantorovich method with single-term displacement function is ineffective.

• Earthquakes and Structures
• /
• 제3권3_4호
• /
• pp.413-434
• /
• 2012

This paper presents a displacement-based seismic design method with damage control, in which the targets for the considered performance level are set as displacements and a damage distribution is proposed by the designer. The method is based on concepts of basic structural dynamics and of a reference single degree of freedom system associated to the fundamental mode with a bilinear behaviour. Based on the characteristics of this behaviour curve and on the requirements of modal spectral analysis, the stiffness and strength of the structural elements of the structure satisfying the target design displacement are calculated. The formulation of this method is presented together with the formulations of two other existing methods currently considered of practical interest. To illustrate the application of the proposed method, 5 reinforced concrete plane frames: 8, 17 and 25 storey regular, and 8 and 12 storey irregular in elevation. All frames are designed for a seismic demand defined by single earthquake record in order to compare the performances and damage distributions used as design targets with the corresponding results of the nonlinear step by step analyses of the designed structures subjected to the same seismic demand. The performances and damage distributions calculated with these analyses show a good agreement with those postulated as targets.

• 한국기계가공학회지
• /
• 제20권8호
• /
• pp.11-17
• /
• 2021

This paper describes the time rate of change of dynamic response of a cantilever beam inserted with a damping element, such as bonding, which is excited under a general force at various locations. A sensitivity analysis was performed in a finite element model to show that two types of second-order algebraic governing equations were used to predict the rate of change of dynamic displacement: one is related to the modal coordinate linked to a physical coordinate, and the other to the design parameter of the time rate of change of displacement. The sensitivity differential equation formulation includes more complicated terms compared with that of the undamped cantilever beam. The sensitivities of the dynamic response were observed by changing the location of the excitation force, displacement extraction, and cross-sectional area of the beam. The analytical results obtained by this suggested theory showed a relatively good agreement when compared with those obtained using the commercial finite element program. The suggested analysis procedure enables the prediction of the response sensitivity for any finite element model of the dynamic system.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• 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.

• International Journal of Aerospace System Engineering
• /
• 제3권1호
• /
• pp.10-16
• /
• 2016

A trigonometric Layerwise higher order shear deformation theory (TLHSDT) is developed and implemented for free vibration and buckling analysis of laminated composite and sandwich plates by analytical and finite element formulation. The present model assumes parabolic variation of out-plane stresses through the depth of the plate and also accomplish the zero transverse shear stresses over the surface of the plate. Thus a need of shear correction factor is obviated. The present zigzag model able to meet the transverse shear stress continuity and zigzag form of in-plane displacement continuity at the plate interfaces. Hence, botheration of shear correction coefficient is neglected. In the case of analytical method, the governing differential equation and boundary conditions are obtained from the principle of virtual work. For the finite element formulation, an efficient eight noded $C^0$ continuous isoparametric serendipity element is established and employed to examine the dynamic analysis. Like FSDT, the considered mathematical model possesses similar number of variables and which decides the present models computationally more effective. Several numerical predictions are carried out and results are compared with those of other existing numerical approaches.

• 한국전산구조공학회논문집
• /
• 제13권2호
• /
• pp.209-220
• /
• 2000

보강된 판 및 쉘구조의 기하학적 비선형해석을 수행하기 위하여, total lagrangian formulation에 근거한 증분 평형방정식을 적용하고, 강도행렬 산정시 회전각의 2차항을 포함시켜 기하학적 비선형 해석시 해의 수렴성을 향상시켰으며, 보강된 쉘 구조의 해석시 보강재를 쉘 요소로 모델링하고 주부재와 보강재의 연결점에서 일반적인 변환관계를 이용하였다. 등매개 쉘 유한요소의 단점인 locking 현상을 극복하기 위하여 가정 변형률장을 적용하여 감차적분 또는 선택적분시 나타날 수 있는 제로 에너지 모드를 제거하였다. 수치해석 예제를 통하여 가정 변형률장에 근거한 쉘유한요소에 대한 효율성 및 적용성을 확인하였다.

• Advances in nano research
• /
• 제4권4호
• /
• pp.251-264
• /
• 2016

In this paper, a simple and refined nonlocal hyperbolic higher-order beam theory is proposed for bending and vibration response of nanoscale beams. The present formulation incorporates the nonlocal scale parameter which can capture the small scale effect, and it considers both shear deformation and thickness stretching effects by a hyperbolic variation of all displacements across the thickness without employing shear correction factor. The highlight of this formulation is that, in addition to modeling the displacement field with only two unknowns, the thickness stretching effect (${\varepsilon}_z{\neq}0$) is also included in the present model. By utilizing the Hamilton's principle and the nonlocal differential constitutive relations of Eringen, the equations of motion of the nanoscale beam are reformulated. Verification studies demonstrate that the developed theory is not only more accurate than the refined nonlocal beam theory, but also comparable with the higher-order shear deformation theories which contain more number of unknowns. The theoretical formulation proposed herein may serve as a reference for nonlocal theories as applied to the static and dynamic responses of complex-nanobeam-system such as complex carbon nanotube system.