Open topic with navigation

Calculation stiffness of 2D element

Members which are not expected to be loaded above the level which would cause the tensile strength of the concrete to be exceeded anywhere within the member should be considered to be uncracked. Members which are expected to crack, but may not be fully cracked, will behave in a manner intermediate between the uncracked and fully cracked conditions. New stiffness (stiffness with taking into account cracking) is calculated in centre of each 2D element.

Two types of stiffness are calculated:

Short-term stiffness - is calculated using 28 days modulus of elasticity E_c = E_cm, it follows that value of stiffness is loaded directly from properties of the concrete material

Long-term stiffness - is calculated using effective E modulus based on creep coefficient for acting load, it follows E_c = E_c,eff = E_cm/(1+j).

Calculation effective modulus of elasticity is based on equation 5.27 in EN 1992-1-1, but instead of effective creep coefficient j_ef, only creep coefficient j is used

The following procedure is used for calculation stiffens of 2D element

• Calculation of principal stresses of 2D element for both surfaces

• Calculation angle of principal stresses at both surfaces

• Calculation of final value of principal stress

a=as1+ if s₁₊≥s_1-

a=as1- otherwise

• Recalculation internal forces to direction of principal stress a

where n_x,n_y,n_xy,m_x,m_y,m_xy are 2D forces in centre of 2D element

• Recalculation area of reinforcement to direction of of principal stress a

A_s(a) = A_s×cos²(a-a_s)

where A_s,a_s is area and angle of longitudinal reinforcement

• Calculation of non-linear stiffness in first principal direction according to procedure as for 1D element, see "Calculation stiffness of 1D element"
  • for rectangular cross-section (b =1m, h = thickness of 2D element in centre of gravity)
  • for internal forces N = n(a) , M_y= m(a) and M_z=0 according procedure as for 1D element
• Calculation of non-linear stiffness in second principal direction according to procedure as for 1D element, see "Calculation stiffness of 1D element"
  • for rectangular cross-section (b =1m, h = thickness of 2D element in centre of gravity)
  • for internal forces N = n(a+90) , M_y= m(a+90) and M_z=0 according procedure as for 1D element

The four type of stiffnesses is calculated for each 1D element and each dangerous combination:

The following stiffnesses are changes in stiffness matrix for 2D element:

D11 = EIy₁

D22 = EIy₂

D33 = 0.5×(1-m)×(D11×D22)^0.5

D44 = G×h/1.2

D55 = G×h/1.2

D12 = m• (D11×D22)^0.5

d11= EA₁

d22= EA₂

d33 = G×h

d12 =m× (d11×d22)^0.5

G is shear modulus of the concrete calculated according to formula G = 0.5×E_c/(1+m)

m is Poisson coefficient of the concrete loaded from material properties of the concrete

Eccentricity of stiifnness (distance between centre of gravity of concrete cross-section and centre of gravity of cracked transformed cross-section) is not taken into account in current version

Open topic with navigation