ME 1401 - FINITE ELEMENT ANALYSIS Two Marks Questions With Answers 2014

VII Semester Mechanical Engineering

ME 1401 - FINITE ELEMENT ANALYSIS

Short Q & A

1) What is meant by finite element analysis?

Finite element method is a numerical method for solving problems of engineering mathematical physics.

In the finite element method, instead of solving the problem for

the entire body in one operation, we formulate the equations for each finite element and combine them to obtain the solution of the whole body.

2) What is meant by finite element?

A small unit having definite shape of geometry and nodes is called finite element.

3) State the methods of engineering analysis.

There are three methods of engineering analysis. They are:

1. Experimental methods.

2. Analytical methods.

3. Numerical methods or approximate methods.

4) Give examples for the finite element.

1. One dimensional elements: (a) Truss element

(b) Bar, Beam element

2. Two dimensional elements (a) Triangular element

(b) Rectangular elements

3. Three dimensional elements (a) Tetrahedral elements.

(b) Hexahedral elements.

5) What is meant by node or Joint?

Each kind of finite element has a specific structural shape and is interconnected with the adjacent elements by nodal points or nodes. At the nodes, degrees of freedom are located. The forces

will act only at nodes and not at any other place in the element.

6) What do you mean by discretization?

Discretization is the basis of finite element method. The art of subdividing a structure into a convenient number of smaller components is known as discretization.

7) What are the types of boundary conditions?

There are two types of boundary conditions, they are

Primary boundary condition. Secondary boundary condition.

8) What are the three phases of finite element method.

The three phases are 1. Preprocessing

2. Analysis

3.Postprocessing

9) What is structural and non-structural problem?

Structural problem: In structural problems, displacement at each nodal point is obtained. By using these displacement solutions, stress and strain in each element can be calculated.

Non Structural problem: In non structural problem,temperatures or fluid pressure at each nodal point is obtained. By using these values, Properties such as heat flow, fluid flow, etcfor each element can be calculated.

10) What are the methods are generally associated with the finite element analysis?

The following two methods are generally associated with the finite element analysis. They are

1. Force method.

2. Displacement or stiffness method

11) Explain force method and stiffness method?

In force method, internal forces are considered as the unknowns of the problem. In displacement or stiffness metho, displacement of the node are considered as the of the problem. Among them two approaches,dislplacement method is desirable.

12) What is polynomial type of interpolation functions are mostly used in FEM?

The polynomial type of interpolation functions are mostly used due to the following reasons:

1. It is easy to formulate and computerize the finite element equations.

2. It is easy to perform differentiation or integration.

3. The accuracy of the results can be improved by increasing the order of the polynomial .

13) Name the variational methods.

1.Ritz method.

2. Rayleigh – Ritz method

14) Name the weighted residual methods.

1.Point collocation method.

2. subdomain collocation method.

3. Least square method

4. galerkin’s method

15) What is meant by post processing?

Analysis and evaluation of the solution results is referred to as post processing. Post processor computer programs help the user to interpret the results by displaying them in graphical form.

16) What is Rayleigh ritz method?

Rayleigh ritz method is a integral approach method which is useful for solving complex structural problems, encountered in finite element analysis. This method is possible only if a suitable functional is available.

17)What does assemblage mean?

The art of subdividing a structure into a convenient number of smaller components is known as discretization. These smaller components are then put together. The process of uniting the various elements together is called assemblage.

18)What is meant by DOF?

When the force or reaction acts at nodal point, node is subjected to deformation. The deformation includes displacement, rotations, and/or strains. These are collectively known as degrees of freedom (DOF).

19)What is aspect ratio?

Aspect ratio is defined as the ratio of the largest dimension of the element to the smallest dimension. In many cases, as the aspect ratio increases, the inaccuracy of the solution increases. The conclusion of many researches is that the aspect ratio should be close to unity as possible.

20)What is truss element?

The truss elements are the part of a truss structure linked together by point joints, which transmit only axial force to the element.

21) List the two advantages of post processing?

1.Required result can be obtained in graphical form.

2. Contour diagrams can be used to understand the solution easily and quickly.

22) If a displacement field in x direction ils given by u=2x2+4y2+6xy.Determine the strain in x direction.

U=2x2+4y2+6xy

Strain,e= δu/δx= 4x +6y

23) What are h and p versions of finite element method?

H version and p versions are used to improve the accuracy of the finite element method.

In h versions, the order of polynomial approximation for all elements is kept constant and the number of elements is increased.

In p version, the number of elements is maintained constant and the order of polynomial approximation of element is increased.

24) During discretization, mention the places where it is necessary to place a node

The following places are necessary to place a node during discretization process.

1. Concentrated load-acting point.

2. Cross section changing point

3. Different material interjunction point

4. Sudden change in load point.

25) What is the difference between static and dynamic analysis?

Static analysis: The solution of the problem does not vary with time is known as static analysis.

Example: Stress analysis on a beam.

Dynamic analysis: The solution of the problem varies with time is known as dynamic analysis.

Example: vibration analysis problems.

26) Name the four FEA softwares?

1. ANSYS

2. NASTRAN

3. COSMOS

4. NISA

27) Differentiate between global and local axes.

Local axes are established in an element. Since it is in the element level, they change with the change in orientation of the element. The direction differs from element to element.

Global axes are defined for the entire system. They are same in direction for all the elements even though the elements are differently oriented.

28) Distinguish between potential energy function and potential energy functional.

If a system has finite number of degrees of freedom (q1 q2 and q3)

then the potential energy is expressed as, π = f (q1 q2 and q3)

It is known as function.

If a system has infinite degrees of freedom, then the potential energy is expressed as,

π = ∫ f .( x, y, dy/dx, d2y/dx2,….)dx

It is known as functional.

29) What are the types of loading acting on the structure?

There are three types of loading acting on the body. They are:

1. Body force (f)

2. Traction force (T)

3. Point load (P)

30) Define body force (f).

A body force is distributed force acting on every elemental volume of the body.

Unit: Force per unit volume. Example: Self-weight due to gravity

31) Define traction force (T)

Traction force is defined as a distributed force acting on the surcace of the body

Unit: force per unit area.

Examples: Frictional resistance, viscous drag, surface shear etc.

32) What is point load (P)

Point load is force acting at a particular point, which causes displacement.

33) What are the basic steps involved in the finite element modeling.

Finite element modeling consists of the following:

1.Discretization of structure

2. Numbering of nodes.

34) What is discretization?

The art of subdividing a structure into a convenient number of smaller componentsis known as descretization.

35) What is the classification of co-ordinates?

The co ordinates are generally classified as follows:

1. Global co-ordinates

2. Local co-ordinates

3. Natural co-ordinates

36) What is Global co-ordinates?

The points in the entire structure are defined using co-ordinate system is known as global co-ordinate system.

clip_image002clip_image003clip_image004clip_image005clip_image006Example :

1 2 3 4 5

clip_image0071 2 3 4 5 6

37) What is natural co-ordinates?

A natural co-ordinate system is used too defined any point inside the element by a set of dimensionless numbers, whose magnitude never exceeds unity.This system is very useful in assembling of stiffness matrices.

38) Define shape function.

In finite element method,field variables within an element are generally expressed by the following approximate relation

N3 are the

Φ(x,y) = N1(x,y) φ1 + N2(x,y) φ2 +N3(x,y) φ3

Where φ1, φ2 and φ3 are the values of the field variables at the nodes and N1 ,N2 and

interpolation functions.

N1 ,N2 and N3 also called shape function because they are used to express the

geometry or shape of the element.

39) What are the characteristics of shape function?

The characterstics of shape function are as follows:

1. The shape function has unit value at one nodal point and zero value at other nodal points.

2. The sum of shape function is equal to one.

40) Why polynomial are generally used as shape function?

Polynomials are generally used as shape function due to the following reasons.

1. Differentiation and integration of polynomial are quite easy.

2. The accuracy of the results can be improved by increasing the order of the polynomial.

3. It is easy to formulate and computerize the finite element equations.

41) How do you calculate the size of the global stiffness matrix?

Global stiffness matrix size = Number of nodes X degree of freedom per node

1) Give th general expression for element stiffness matrix.

Stiffness matrix [k] = ∫ [B] T [D] [B] dv

[B] – Strain displacement matrix [row matrix]

[D] – Stress, Strain relationship matrix [Row matrix]

42) Write down the expression of stiffness matrix for one dimensional bar element.

clip_image008clip_image009clip_image010Stiffness matrix [k] = AE 1 -1

clip_image011L -1 1

A - Area of the bar element.

E –Youngs modulus of bar element

L – lenth of the bar element.

43) State the properties of a stiffness matrix.

The properties of a stiffness matix [k] are

1. It is symmetric matrix.

2. The sum of elements in any column must be equal to zero.

3. It is an unstable element.

So, the determinant is equal to zero.

44) Write down the general finite element equation.

General finite element equation is,

{F} = [K] {u}

{F} - Force vector [column matrix] [k] - Stiffness matrix [ row matix]

{u} - Degrees of freedom [coloumn matrix]

45) Write down the finite element equation for one dimensional two noded bar element.

clip_image012clip_image013clip_image014clip_image015clip_image016clip_image017The finite element equation for one-dimensional two noded bar element is,

F1

 

AE

 

1

-1

u1

F2

=

L

-1

 

1

u2

46) What is truss?

A truss is defined as a structure, made up of several bars, riveted or welded together.

47) State the assumptions are made while finding the forces in a truss.

The following assumptions are made while finding the forces in a truss.

1. All the members are pin jointed.

2. The truss is loaded only at the joints.

3. The self-wight of the members is neglected unless stated.

clip_image01848) Write down the expression of stiffness matrix for a truss element.

clip_image019Stiffness matrix [k] = Ae Ee l2 lm -l2 -lm

clip_image020Lm m2 -lm -m2

L -l2 -lm l2 lm

clip_image021clip_image022-lm -m2 lm m2

A – Area

E – Youngs modulus

le - Length of the element l,m – Direction cosines

49) Write down the expression of shape function N and displacement u for one- dimensional bar element.

For one dimensional bar element

Displacement function, u= N1 u1 + N2 U2

Where, shape function N1 = l-x / l

shape function N2= x/l

50) Define total potential energy.

The total potential energy π of an elastic body, is defined as the sum of total strain energy U and potential energy of the external forces,(W).

Total potential energy, π = Strain energy (U) + Potential energy of the external forces (W).

51) State the principle of minimum potential energy.

The principle of minimum potential energy states: Among all the displacement equations that satisfy internal compatibility and the boundary conditions, those that also satisfy the equations of equilibrium make the potential energy a minimum in a stable system.

52) What is the stationary property of total potential energy?

If a body is in equilibrium, its total potential energy π is stationary.

For stable equilibrium, δ2π >0 , other wise π is minimum for stable equilibrium. For neutral equilibrium, δ2π = 0. In this case π is unchanging.

For unstable equilibrium, δ2π < 0 , other wise π is maximum.

53) State the principle of virtual work?

A body is in equilibrium if the internal virtual work equals the external virtual work for every kinematical admissible displacement field.

54) Distinguish between essential boundary conditions and natural boundary conditions.

There are two types of boundary conditions. They are:

1. Primary boundary condition (or) Essential boundary condition

The boundary condition, which in terms of field variable, is known as primary boundary condition.

2. Secondary boundary condition or natural boundary conditions:

The boundary conditions, which are in the differential form of field variables, are known as secondary boundary condition.

Example: A bar is subjected to axial load as shown in fig.

U= 0

EA du/dx = p

In this problem, displacement u at node 1 = 0, that is primary boundary condition.

clip_image023EA du/dx = P, that is secondary boundary condition.

55) What are differences between boundary value problem and initial value problem?

The solution of differential equation is obtained for physical problems, which satisfies some specified conditions known as boundary conditions.

The differential equation together with these boundary conditions, subjected to a boundary value problem.

The differential equation together with initial conditions subjected to an initial value problem.

Examples: Boundary value problem. d2y/dx2 - a(x) dy/dx – b(x)y –c(x) = 0

with boundary conditions, y(m) = α and y(n) = β

initia l value problem, ax2 + bx + c =0

Boundary conditions: x(0) = 0

X(0) = 7

clip_image024

56) How do you define two-dimensional elements?

Two dimensional elements are defined by three or nodes in a two dimensional plane(ie x,y plane). The basic element useful for two dimensional analysis is the triangular element.

57) What is CST element?

Three-nodded triangular element is known Constant Strain Triangle (CST) which is shown in fig. it has six unknown displacement degrees of free (u1 v1 u2 v2, u3

v3 ). The element is called CST because it has a constant strain through it.

clip_image025V3

clip_image0263 u3

V2

2

U2

V1

1 U1

Merit:

Calculation of stiffness matrix is easier. Demerit:

The strain variation with in the element is considered as constant. So, the results will be poor.

58) What is LST element?

Six noded triangular elements are known as linear strain triangle (LST), which is shown in fig. It has twelve unknown displacement degrees of freedom. The displacement functions for the element are quadratic instead of linear as in the CST.

3

6 5

1 4 2

59) What is QST element?

Ten noded triangular elements is known as quadratic strain triangle (QST), which is shown if, fig. it is also called cubic displacement triangle.

9 10 6

clip_image0271 4 5 2

clip_image02860) What is meant by plane stress analysis?

Plane stress is defined to be a state in which the normal stress (σ) and shear stress directed perpendicular to the plane is assumed to be zero.

clip_image02961) Define plane strain analysis

Plane strain is defined to be a state of strain in which the strain normal to the xy plane and the shear strains are assumed to be zero.

clip_image030clip_image03162) Write a displacement function equation for CST element.

clip_image032Displacement function u = u (x,y) N1 0 N2 0 N3 0 U1

clip_image033clip_image013[1]clip_image034v(x,y) 0 N1 0 N2 0 N3 V1

U2

V2

U3

clip_image035V3

Where N1, N2, N3 are shape functions.

63) Write a strain-displacement matrix for CST element.

Strain displacement matrix for CST element is

Q1 0 q2 0 q3 0

clip_image036clip_image037[B] = 1/2A 0 r1 0 r2 0 r3

clip_image038clip_image039R1 q1 r2 q2 r3 q3

Where A = Area of the element

Q1 = y2 –y3 q2 = y3 –y1 q3 = y1-y2

R1= x3 –x2 r2= x1 –x3 r3= x2 –x1

64) Write down the stress strain relationship matrix for plane strain condition.

For plane strain problems, stress strain relationship matrix is,

(1-V) V 0

clip_image040clip_image041clip_image042clip_image043clip_image044[ D] = E 0 ( 1-V) 0 (1+V)(1-2V) 0 0 1-2V/2

65) Write down the stiffness matrix equation for two-dimensional CST element.

Stiffness matrix, [k] = [B]T [D] [B] A t

[B] = Strain displacement matrix

[D] = Stress strain matrix

A = Area of the element

T = Thickness of the element

66) Write down the stress strain relation ship matrix for plane stress condition.

clip_image045For plane stress problem, stress strain relation ship matrix is

clip_image046clip_image047clip_image0481-v v 0 [D] = E 0 1-v 0

(1+v)(1-2v) 0 0 1-2v/2

67) Write down the expression for the shape function for a constant strain triangular element.

For CST element,

Shape function, N1 = P1 + q1x + r1y

2A

N2 = p2+q2x+r2y

2A

Where P1= x2y3 – x3 y2

P2= x3y1 – x1 y3

P3= x1y2 – x2 y1

q1= y2 – y3

q2= y3 – y1

q3= y1 – y2

N3 = p3+q3x+r3y

2A

r1= x3 – x2 r2= x1 – x3 r3= x2 – x1

68) What is axisymmetric element?

Many three dimensional problems in engineering exhibit symmetry about an axis of rotation. Such types of problems are solved by a special two-dimensional element called the axisymmetric element.

69) What are the conditions for a problem to axisymmetric?

1. The problem domain must be symmetric about the axis of revolution.

2. All boundary conditions must be symmetric about the axis of revolution.

3. All loading condition must be symmetric about the axis of revolution.

70) Write down the displacement equation for an axisymmetric triangular element.

clip_image049clip_image050clip_image051clip_image052clip_image053clip_image054Displacement function, u(r,z) = u(r,z) N1 0 N2 0 N3 0 u1

clip_image042[1]w(r,z) 0 N1 0 N2 0 N3

w1

u2

w2 u3 w3

71) Write down the shape function for an axisymmetric triangular element.

clip_image055clip_image056clip_image057N1 = α1 x+ β1y +Ѵ1z

2A

N2 = α2 x+ β2y +Ѵ2z

2A

N3 = α3 x+ β3y +Ѵ3z

2A

Where α1 = r2z3 – r3 z2

α2 = r3z1 – r1 z3

α3 = r1z2 – r2 z1 β1= y2 – y3

β 2= z3 – z1

β 3= z1 – z2

Ѵ1= r3 – r2 Ѵ 2= r1 – r3 Ѵ 3= r2 – r1

72) Give the strain-displacement matrix equation for an axisymmetric triangular element.

Strain-displacement matrix,

β1 0 β 2 0 β 3

0

clip_image058[B]= 1/2A α1/r1+ β1 + Ѵ1z/r 0 α1/r1+ β1 + Ѵ1z/r 0 α1/r1+ β1 + Ѵ1z/r 0

0 Ѵ1 0 Ѵ2 0

Ѵ3

clip_image059Ѵ1 Ѵ2 Ѵ3

73) Write down the stress - strain relation ship matrix for an axisymmetric triangular element.

clip_image060clip_image0611-Ѵ Ѵ Ѵ 0

clip_image013[2]E Ѵ 1-Ѵ Ѵ 0

clip_image047[1]Stress- strain relationship matrix,[D] (1+ Ѵ)(1-2 Ѵ) Ѵ Ѵ 1- Ѵ

0

Ѵ/2

0 0 0 1-2

E – youngs modulus

V – poisson ratio

74) Give the stiffness matrix equation for an axisymmetric triangular element.

Stiffness matrix, [k] = 2π r A [B] T [D] [B] Where, co-ordinate r = r1 + r2 + r3 /3

A – area of the triangular element.

75) What are the ways in which a three dimensional problem can be reduced to a two dimensional approach.

1. Plane stress: One dimension is too small when compared to other two dimensions.

Example: gear – thickness is small

2. Plane strain: one dimension is too large when compard to other two dimensions example : Long pipe [length is long compared to diameter]

3. Axisymmetric : geometry is symmetric about the axis.

Example: cooling tower

76) Calculate the jacobian of the transformation J for the triangular element shown if fig.

clip_image062clip_image063z 1 (3,6)

clip_image064clip_image0651(2,3) 2(5,3)

R

r1 =2 ; r2 =5; r3 = 3 z1 = 3 ;z2 = 3; z3 = 6

clip_image066clip_image067clip_image068clip_image069j = r1 - r3 z1- z3 = -1 -3

clip_image070clip_image071clip_image072r2 - r3 z2- z3 2 -3

J = 3+6 = 9 units

77) What is the purpose of isoparametric elements?

It is difficult to represent the curved boundaries by straight edges finite elements. A large number of finite elements may be used to obtain reasonable resemblance between original body and assemblage. In order to overcome this drawback, isoparametric elements are used i.e. for problems involving curved boundaries; a family of elements known as “isoparametic elements” is used

78) Write down the shape function for 4 noded rectangular elements using natural co-ordinate system.

Shape functions: N1 = ¼ (1-Є) (1-η) N2 = ¼ (1+Є) (1-η)

N3 = ¼ (1+Є) (1+η)

N4 = ¼ (1-Є) (1+η)

Where, Є and η are natural co-ordinates.

79) Write down the jacobian matrix for four noded quadrilateral element.

clip_image073clip_image074Jacobian matrix,[ J ] = J11 J12

J21 J22

Where, J11 = ¼[-(1- η)x1 + (1- η)x2 + (1+η)x3 - (1+ η)x4 ] J12 = ¼[-(1- η)y1 + (1- η)y2 + (1+η)y3 - (1+ η)y4 ]

J21 = ¼[-(1- Є)x1 - (1+ Є)x2 + (1+ Є)x3 + (1- Є)x4 ] J22 = ¼[-(1- Є)y1 - (1+ Є)y2 + (1+ Є)y3 + (1- Є)y4 ]

Where, Є and η are natural co-ordinates.

x1 x2 x3 x4 y1 y2 y3 y4 are Cartesian coordinates.

80) Write down the stiffness matrix equation for four noded isoparametric quadrilateral elements.

1 1

Stiffness matrix, [k] = t ∫ ∫ [B]T [D] [B] X |J| X δ Є X δ η

-1 -1

Where t = thickness of the element

|J|=Determinant of the jacobian Є and η are natural co-ordinates [B] Strain displacement matrix

[D] Stress strain relationship matrix

81) Write down the element force vector equation for four noded quadrilateral elements.

clip_image075clip_image076Force vector, {F} = [N] T Fx

Fy

N is the shape function

Fx is the load or force in x direction

Fy is the force on y direction

82) Write down the Gaussian quadrature expression for numerical integration.

Gaussian quadrature expression

1 n

∫ f(x) dx = ∑ wi f(xi)

-1 i=1

Wi weight function

F(xi) values of the function at pre-determined sampling points.

83) Define super parametric element.

If the number nodes used for defining the geometry is more than number of nodes used for defining the displacements is known as super parametric element.

84) What is meant by sub parametric element?

If the number of nodes used for defining the geometry is less than number of nodes used for defining the desplacements known as isoparametric element.

85) What is meant by iso parametric element?

If the number of nodes used for defining the geometry is same as number of nodes used for defining the displacements is know as isoparametric element.

86) Is beam element an isoparametric element?

Beam element is not an isoparametric element since the geometry and displacements are defined by different order interpolation functions.

87) What is the difference between natural co-ordinates and simple natural co- ordinate?

A natural co-ordinate is one whose value lies between zero and one. Examples: L2 = x/l; l = (1-x/l)

Area co-ordinates : L1 = A1/A ; L2 = A2/A ; L3 = A3/A

A simple natural co ordinates is one whose value lies between -1 to +1

88) Give examples for essential (forced or geometric) and non-essential (natural)

boundary conditions.

The geometric boundary conditions are displacement, slopes, etc. the natural boundary conditions are bending moment, shear force,etc.

Descriptive type

89) Explain the general steps in FEA with the help of a flowchart?

90) A beam AB of span ‘L’ simply supported at ends and carrying a concentrated load

‘W’ at the center ‘C’. Determine the deflection at midspan by using the Rayleigh- Ritz method.

91) Solve the equations using Gauss- Elimination method

2x + 4y +2z = 15

2x + y + 2 z = -5

4x + y – 2z = 0

92) Describe the four types of weighted residual method.

93) Derive the stiffness matrix [K] for the truss element

94) Derive the shape function for one-dimensional bar element.

95) Using finite element, find the stress distribution in a uniformly tapering bar of circular cross sectional area 3cm2 and 2 cm2 at their ends, length 100mm, subjected to an axial tensile load of 50 N at smaller end and fixed at larger end. Take the value of Youngs modulus as 2 x

105N/mm2.

96) (i) Explain the Galerkin’s method.

(ii) Explain the Gaussian elimination.

97) Derive the constitutive matrix for 2D element.

98) Derive the expression for the stiffness matrix for an axisymmetric shell element.

99) Explain the terms plane stress and plane strain conditions. Give the constitutive laws for these cases.

100) Derive the element stiffness matrix for a linear isoparametric quadrilateral element.

101) Evaluate the integral by using Gaussian quadrature ∫ x2dx.

-------------------------------------------------------------------

Anna University, Chennai

Anna_University,_Chennai_logo

25 comments:

  1. really appreciate your effort, and i was surprised that there were no mistakes cool. Thank you

    ReplyDelete
  2. TESTIMONY ON HOW I GOT MY LOAN AMOUNT FROM A RELIABLE AND TRUSTED LOAN COMPANY LAST WEEK. Email for immediate response drbenjaminfinance@gmail.com

    Hello everyone, My name is Mrs. Carolin Glowski, I'm from Europe, am here to testify of how i got my loan from BENJAMIN LOAN FINANCE after i applied Two times from various loan lenders who claimed to be lenders right here this forum, i thought their lending where real and i applied but they never gave me loan until a friend of mine introduce me to {Dr. Benjamin Scarlet Owen} the C.E.O of BENJAMIN LOAN FINANCE who promised to help me with a loan of my desire and he really did as he promised without any form of delay, I never thought there are still reliable loan lenders until i met {Dr. Benjamin Scarlet Owen}, who really helped me with my loan and changed my life for the better. I don't know if you are in need of an urgent loan also, So feel free to contact Dr. Benjamin Scarlet Owen on his email address drbenjaminfinance@gmail.com


    THANKS

    ReplyDelete
  3. Hello. I want to advise you a very good and reliable site where you will find everything about .NET development. They are real professionals with experience all over the world.

    ReplyDelete
  4. What is the main advantage of using the Finite Element Method (FEM) in solving engineering problems? regards Telkom University

    ReplyDelete