• Steel and Composite Structures
• /
• 제15권5호
• /
• pp.481-505
• /
• 2013

This paper focuses on thermal post-buckling analysis of functionally graded beams with temperature dependent physical properties by using the total Lagrangian Timoshenko beam element approximation. Material properties of the beam change in the thickness direction according to a power-law function. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces and therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In this study, the differences between temperature dependent and independent physical properties are investigated for functionally graded beams in detail in post-buckling case. With the effects of material gradient property and thermal load, the relationships between deflections, critical buckling temperature and maximum stresses of the beams are illustrated in detail in post-buckling case.

• Structural Engineering and Mechanics
• /
• 제44권1호
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• 한국입자에어로졸학회지
• /
• 제13권4호
• /
• pp.173-182
• /
• 2017

An Eulerian-Lagrangin approach is used to compute particle dispersion from a power plant chimney. For air flow, three-dimensional incompressible filtered Navier-Stokes equations are solved with a subgrid-scale model by integrating the Newton's equation, while the dispersed phase is solved in a Lagrangian framework. The velocity ratios between crossflow and a jet of 0.455 and 0.727 are considered. Flow fields and particle distribution of both cases are evaluated and compared. When the velocity ratio is 0.455, it demonstrates a Kelvin-Helmholtz vortex structure above the chimney caused by the interaction between crossflow and a jet, whereas the other case shows flow structures at the top of the chimney collapsed by fast crossflow. Also, complex wake structures cause different particle distributions behind the chimney. The case with the velocity ratio of 0.727 demonstrates strong particle concentration at the vortical region, whereas the case with the velocity ratio of 0.455 shows more dispersive particle distribution. The simulation result shows similar tendency to the experimental result.

• 대한기계학회논문집
• /
• 제12권5호
• /
• pp.992-1002
• /
• 1988

본 연구에서는 Monte Carlo 새뮬레이션 기법을 이용하여 지반가속도의 스펙트 럼 밀도함수(power spectral densities)로부터 여섯 성분의 지반가속도 시간이력곡선 을 얻고, 이들을 입력 데이터로 하여 운동방정식에 Newmark의 직접적분법을 이용하여 회전축-베어링 시스템의 응답상태벡터(response state vector)를 얻기로 한다. 충분 히 많은 수의 지반가속도 시간이력곡선을 시뮬레이션하고, 각 경우에 대응하는 응답상 태벡터들을 얻은 다음 일반적인 통계학 방법을 적용하여 평균함수, 표준편차 및 r.m.s (root mean square)등을 얻는다.

• Journal of Theoretical and Applied Mechanics
• /
• 제1권1호
• /
• pp.97-109
• /
• 1995

In this paper, a finite viscoelastic continuum model for rubber and its finite element analysis are presented. This finite viscoelatic model based on continuum mechanics is an extended model of Johnson and Wuigley's 1-D model. In this extended model, continuum based kinematic measures are rigorously defied and by using this kinematic measures, elastic stage law and flow rule are introduced. In kinematics, three configuration are introduced. In kinematics, three configuration are introduced. They are reference, current and virtual visco configurations. In elastic state law, it is assumed that at a certain time, there exists an elastic potential which describes the recoverable elastic energy. From this elastic potential, elastic state law is derived. The proposed flow rule is based on phenomenological observation. The flow rule gives precise relaxation response. In finite element approximation, mixed Lagrangian description is used, where total and similar method of updated Lagrangian descriptions are used together. This approach reduces the numerical job and gives simple nonlinear syatem of equations. To satisfy the incompressible condition, penalty-type modified Mooney-Rivlin energy function is adopted. By this method nearly incompressible condition is obtain the virtual visco configuration. For verification, uniaxial stretch tests are simulated for various stretch rates. The simulated results show good agreement with experiments. As a practical experiments. As a preactical example, pressurized rubber plate is simulated. The result shows finite viscoelastic effects clearly.

• 한국기계연구소 소보
• /
• 통권14호
• /
• pp.17-31
• /
• 1985

The Primary object of this study is to predict rigid carbody's vertical and pitch acceleration and/or displacement frequency response to vertical sinusoidal rail surface irregularities for any specified point of the carbody, and to verify the predictions by means of experiments. The developed computer program also calculates vertical and pith transmissibilities and acceleration spectra. This model can be used for first order analysis of ride behabior. it's main advantage is its simplicity and ease of use. This model can be used for first order analysis of ride quality behabior. It's main advantage is its simplicity and ease of use. The model was designed with 6 degreed of freedom. Equations of motion were derived by Lagrangian method. This calculation was applied to the vertical dynamic analysis in order to pursue a possible improvement of the dynamic performance of co-co locomotive, and results were very useful.

• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• Journal of applied mathematics & informatics
• /
• 제23권1_2호
• /
• pp.215-228
• /
• 2007

By using the variational calculus of fractional order, one derives a Hamilton-Jacobi equation and a Lagrangian variational approach to the optimal control of one-dimensional fractional dynamics with fractional cost function. It is shown that these two methods are equivalent, as a result of the Lagrange's characteristics method (a new approach) for solving non linear fractional partial differential equations. The key of this results is the fractional Taylor's series $f(x+h)=E_{\alpha}(h^{\alpha}D^{\alpha})f(x)$ where $E_{\alpha}(.)$ is the Mittag-Leffler function.

• Coupled systems mechanics
• /
• 제8권1호
• /
• pp.19-38
• /
• 2019

The influence of torsional rigidity of hinged flexible appendage on the linear dynamics of flexible spacecrafts with liquid on board was analyzed by considering the spacecraft's main body as a rigid tank, its flexible appendages as two elastically supported elastic beams, and the onboard liquid as an ideal liquid. The meniscus of the liquid free surface due to surface tension was considered. Using the Lagrangian of the spacecraft's main body (rigid tank), onboard liquid, and two beams (flexible appendages) in addition to assuming the system moved symmetrically, the coupled system frequency equations were obtained by applying the Rayleigh-Ritz method. The influence of the torsional rigidity of the flexible appendages on the spacecraft's coupled vibration characteristics was primary focus of investigation. It was found that coupled vibration modes especially that of appendage considerably changed with torsion spring parameter ${\kappa}_t$ of the flexible appendage. In addition, variation of the main body displacement with system parameters was investigated.

• 소성∙가공
• /
• 제4권3호
• /
• pp.222-232
• /
• 1995

The planar anisotripic FEM analysis for predicting earing profiles and draw-in amounts in the deep-drawing process is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-based unit vector and normal contact pressure. The consistent full set of governing relations, which is comprising euilbrium and geometric constraint equations, is appropriately linearized. Barlat's strain-rate potential is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and potential parameters. The linear triangular membrane elements are used for depicting the formed sheet. In the numerical simulations of deep drawing processes of a flat-top cylindrical cup for 2090-T3 aluminum alloy sheet show good agreement with experiments, although some discrepancies were observed in the directional trend of cup height and thickness strains.