Left

CivilFEM Online Help

Official Websites: CivilFEM | CivilFEM Licenses | Ingeciber | Current downloads
CivilFEM for Ansys: Install | Doc | Verif | Instala | Doc(ES) | PDF | PDF(ES) | Support | Soporte | Class3
CivilFEM Apps for Ansys Workbench: Install | Lic | FAQ | Concrete | Steel | Theory | PDF | Support | Soporte | Class3
CivilFEM powered by Marc: Install | Capabilities | Cap(ES) | Doc | Python | Theory | Verif | PDF | Support | Soporte | Class3
Right

CFVR4704 - CALTRANS: Seismic Analysis X-axis

Seismic Analysis on a 1 DOF beam according to CALTRANS. (soil profile type D)

This example compares spectrum accelerations given by CivilFEM to accelerations obtained when solving on a 1 DOF beam and to the correct CALTRANS spectum accelerations.
The model is 1 element beam with only 1 DOF, which direction the CivilFEM spectrum is set (the rest of the DOF are constrained). A different mass is applied at the free end of the beam for the beam period to be equal to the 20 CivilFEM spectrum points. After solving, the mass acceleration is obtained and compared to the CivilFEM spectrum acceleration and to the correct value. Therefore, this example checks the following results:

  • CivilFEM spectrum is defined correctly by comparing to the correct value
  • CivilFEM spectrum is applied correctly and response acceleration is equal to spectrum acceleration

Element types used in the model: BEAM4, MASS21

Needed CivilFEM Modules:
Geotechnical Module
Bridges and Civil Non Linearities Module
Advanced Prestressed Concrete Module

The example can be launched on an educational license.
KEYWORDS
Seismic
CALTRANS

Model Statistics
Number of elements 1
Number of nodes 2
Number of civil materials 1
Number of cross sections 1
Number of shell vertices 0

Log file: CFVR4704.DAT


   FINISH
  ~CFCLEAR,,1
  NomFile='CFVR4704'
  /TITLE, %NomFile%, CALTRANS: Seismic Analysis X-axis

! ----------------------------------------------------------------------
! Model definition
! ----------------------------------------------------------------------
! CivilFEM Setup: Code & Units
~UNITS,,LENG,CM
~UNITS,,TIME,S
~UNITS,,FORC,KP
~CODESEL,,,,,CALTRANS
/PREP7
! Preprocessor
! ----------------------------------------------------------------------
! Materials: A-42
 ~CFMP,1,LIB,STEEL,EA,A42
! Modify density (Rho = 0)
 ~CFMP,1,USER
 ~CFMP,1,DatGen,RHO,,0
! Element Types
  ET,1,BEAM4   ! Type 1: 3D Beam
  ET,2,MASS21  ! Type 2: Mass
! Sections
 ~SSECLIB,1,1,1,1   !IPE 80
 ~BMSHPRO,1,BEAM,1,1,,,4,,0,,

! Nodes
  L = 10	         ! L    : Bar length
  N, 1
  N,10,L

! Elements
  TYPE,1
  MAT ,1
  REAL,1
  EN,1,1,10
  EPLOT

! Parameters
  *GET,AR20,EX,1         		         ! Ex   : Elastic Modulus
  ~CFGET,AR22,SECTION,1,MECHPROP,A,,2    ! Area : Bar area
  ~CFGET,AR23,SECTION,1,MECHPROP,IZZ,,2  ! Izz  : Moment of inertia
  ~CFGET,AR24,SECTION,1,MECHPROP,IYY,,2  ! Iyy  : Moment of inertia
  pi=3.141592654

! Seismic Spectrum
  ~DEFSPEC,ALL,0.3,E1,0,5,1

! ----------------------------------------------------------------------
! DATA CHECK
! ----------------------------------------------------------------------
! Data comparision number
  NComp = 40
  NComp_ch = 0

! Marix dim.
  *DIM,LABEL,CHAR,Ncomp,1
  *DIM,LABEL_CH,CHAR,Ncomp_ch,1
  *DIM,VALUE,,Ncomp,3
  *DIM,VALUE_CH,CHAR,Ncomp_ch,3
  *DIM,TOLER,,Ncomp,2

*DO,ARG1,1,20,1
   /PREP7
   ~CFGET,AR26,SEISM,,SPEC,TH,,ARG1,1   ! Period
  ! Labels
   LABEL(ARG1,1) = 'Sa(%ARG1%)x'
   LABEL(20+ARG1,1) = 'Sa(%ARG1%)x'

   AR25=AR20*AR22/L
   AR27=AR25*(AR26/(2*pi))**2  	! Mass applied at end
   R,2,AR27
   TYPE,2
   MAT, 1
   REAL,2
   EN,2,10
   FINISH
   /SOLU
 ! Solution
 ! ---------------------------------------------------------------------
 ! Displacements
   D,1,,,,,,ALL
   D,10,,,,,,UY,UZ,ROTX,ROTY,ROTZ
   ~MODLSOL,1
   ~CMBMOD,HORIZONT,NONE,NONE,0,0.0001
 ! Query results
   *GET,ARG5,NODE,10,U,X
 ! Correct values
   VALUE(20-ARG1+1,1)=ARG5*AR25/AR27 			 ! Aceleration
 ! Obtained values
   ~CFGET,VALUE(20-ARG1+1,2),SEISM,,SPEC,SDH,,ARG1,1	 ! Y components of linear horizontal spectrum
   ~CFGET,VALUE(40-ARG1+1,2),SEISM,,SPEC,SDH,,ARG1,1	 ! Y components of linear horizontal spectrum
   ! CivilFEM gives the spectrum Normalized by gravity.
   ! To compare it must be changed to acceleration user units.
   VALUE(20-ARG1+1,2) = VALUE(20-ARG1+1,2)*981
   VALUE(40-ARG1+1,2) = VALUE(40-ARG1+1,2)*981

   /PREP7
   EDELE,2
*ENDDO

! Correct values
! ----------------------------------------------------------------------
   VALUE(40,1)=  92.2140
   VALUE(39,1)= 150.0930
   VALUE(38,1)= 276.6420
   VALUE(37,1)= 414.9630
   VALUE(36,1)= 680.8140
   VALUE(35,1)= 830.9070
   VALUE(34,1)= 923.1210
   VALUE(33,1)= 945.6840
   VALUE(32,1)= 924.5175
   VALUE(31,1)= 890.0080
   VALUE(30,1)= 800.6312
   VALUE(29,1)= 728.5564
   VALUE(28,1)= 677.5023
   VALUE(27,1)= 595.1941
   VALUE(26,1)= 535.4232
   VALUE(25,1)= 453.0806
   VALUE(24,1)= 358.0650
   VALUE(23,1)= 358.0650
   VALUE(22,1)= 358.0650
   VALUE(21,1)= 358.0650

! Warning and error tolerances
! ----------------------------------------------------------------------
   *DO,II,1,NComp
     TOLER(II, 1)= 0.1 $ TOLER(II, 2)= 0.1
   *ENDDO

! ----------------------------------------------------------------------
! Results Comparison
! ----------------------------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
Sa(1)x 358.07 358.071.000 0.1
Sa(2)x 358.06 358.071.000 0.1
Sa(3)x 358.06 358.071.000 0.1
Sa(4)x 358.07 358.071.000 0.1
Sa(5)x 453.12 453.121.000 0.1
Sa(6)x 535.43 535.431.000 0.1
Sa(7)x 595.17 595.171.000 0.1
Sa(8)x 677.48 677.481.000 0.1
Sa(9)x 728.59 728.591.000 0.1
Sa(10)x 800.59 800.591.000 0.1
Sa(11)x 889.96 889.961.000 0.1
Sa(12)x 924.5 924.491.000 0.1
Sa(13)x 945.68 945.681.000 0.1
Sa(14)x 923.12 923.121.000 0.1
Sa(15)x 830.91 830.911.000 0.1
Sa(16)x 680.81 680.811.000 0.1
Sa(17)x 414.96 414.961.000 0.1
Sa(18)x 276.64 276.641.000 0.1
Sa(19)x 150.09 150.091.000 0.1
Sa(20)x 92.214 92.2141.000 0.1
Sa(1)x 358.07 358.071.000 0.1
Sa(2)x 358.07 358.071.000 0.1
Sa(3)x 358.07 358.071.000 0.1
Sa(4)x 358.07 358.071.000 0.1
Sa(5)x 453.08 453.121.000 0.1
Sa(6)x 535.42 535.431.000 0.1
Sa(7)x 595.19 595.171.000 0.1
Sa(8)x 677.5 677.481.000 0.1
Sa(9)x 728.56 728.591.000 0.1
Sa(10)x 800.63 800.591.000 0.1
Sa(11)x 890.01 889.961.000 0.1
Sa(12)x 924.52 924.491.000 0.1
Sa(13)x 945.68 945.681.000 0.1
Sa(14)x 923.12 923.121.000 0.1
Sa(15)x 830.91 830.911.000 0.1
Sa(16)x 680.81 680.811.000 0.1
Sa(17)x 414.96 414.961.000 0.1
Sa(18)x 276.64 276.641.000 0.1
Sa(19)x 150.09 150.091.000 0.1
Sa(20)x 92.214 92.2141.000 0.1

Contains proprietary and confidential information of Ingeciber, S.A.