A reinforced concrete cylindrical roof (barrel vault) is subjected to a uniform load. It has a length of 50 feet, a thickness of 0.25 feet, a radius of 25 feet and an 80 degree arc length. A 90 lb/ft2 uniform downward load is to be applied.

Determine the vertical displacement at the center of the free edge i.e. at joint 49 of the model.



Due to symmetry of the structure and loading, a quarter of the cylindrical shell (as shown above) is used as a model. It consists of 6 x 6 mesh of shell elements.

The edge from joint 1 to 7 is simply supported and from 7 to 49 is free.

The boundary conditions along the edges from joint 1 to 43 and 43 to 49 are prescribed for symmetry. Here are the joint restraints specified for preserving symmetry:

Joints 1 to 7 XT YT ZR
Edge Joints 1 to 43 XT YR ZR
Joints 43 to 49 ZT XR YR
The material property of the shell is:
Modulus of Elasticity = 432,000,000 psf
Poisson's Ratio = 0.00
Shear Modulus = 216,000,000 psf
The uniform downward load of 90 lb/ft2 is equivalent to specifying a gravity loading in Global Y-direction with:
Weight Density = 360 lb/ft3
For a shell thickness of 0.25 ft and weight density of 360 lb/ft3, this combination gives the desired 90 lb/ft2 uniform loading in Global Y-direction.


The input for the model should be generated using a cylindrical coordinate system.

Reference

  • MacNeal, R.H. and Harder, R.C., "A Proposed Standard Set of Problems to Test Finite Element Accuracy," Finite Elements in Analysis and Design 1 (1985), pp.3-20, North-Holland.
  • Zienkiewicz, O. C., "The Finite Element Method," McGraw-Hill, 1977.

Result Comparison

Axial Z-Displacement in feet at the Support Joints 1-7

SAP90


Theory


LARSA


OPENBRIM

PERCENT DIFFERENCE

(LARSA vs. THEORY)

PERCENT DIFFERENCE

(OPENBRIM vs. THEORY)

Location=0.00(deg)

0.0000


0.0004


0.0000


0.0000

--

Location=6.67(deg)

0.0005


0.0009


0.0005


0.0005

80%

80%

Location=13.33(deg)


0.0018


0.0020

0.0018


0.0018

11%

11%

Location=20.00(deg)


0.0029


0.0030


0.0029


0.0029

3.45%

3.45%

Location=26.67(deg)

0.0024


0.0021


0.0024

0.0024

12.5%

12.5%

Location=33.33(deg)


-0.0017


-0.0016


-0.0017


-0.0017

5.88%

5.88%

Location=40.00(deg)

-0.0118


-0.0120


-0.0118


-0.0118

1.69%

1.69%


Vertical Y-Displacement in feet at the Center Joints 43-49

SAP90

Theory

LARSA

OPENBRIM

PERCENT DIFFERENCE

(LARSA vs. THEORY)

PERCENT DIFFERENCE

(OPENBRIM vs. THEORY)

Location=0.00(deg)

0.046

0.045


0.047


0.0463

4.26%

2.81%


Location=6.67(deg)

0.031


0.027


0.031


0.0309

12.9%

12.6%

Location=13.33(deg)

-0.013


-0.018


-0.013


-0.0131

38.5%

37.44%


Location=20.00(deg)

-0.078


-0.082


-0.079

-0.0786

3.8%

4.3%

Location=26.67(deg)

-0.155


-0.155


-0.157


-0.1557

1.3%

0.45%

Location=33.33(deg)

-0.234


-0.241


-0.236


-0.2341

2.1%

2.9%

Location=40.00(deg)

-0.307


-0.309


-0.309


-0.3072

0%

0.6%