CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Structural Mechanics (https://www.cfd-online.com/Forums/structural-mechanics/)
-   -   Simple supported beam, incorrect solution (https://www.cfd-online.com/Forums/structural-mechanics/195729-simple-supported-beam-incorrect-solution.html)

optymista93 November 15, 2017 10:36

Simple supported beam, incorrect solution
 
Hello!

I'm new to ANSYS and don't quite get along with this software. I'm trying to do a simple analysis of beam shown below.

https://imgur.com/a/2MAku

I specified beam section with MM: Processor -> Sections -> Beam -> Common Sections and applied R=0.02257(which corresponds to A=0.0016), N=12, T=8.

I keep getting this incorrect solution where beam is bending the opposite way than it's supposed to.

https://imgur.com/a/wZaNT

Here's my code:

Code:

ENTER  /SHOW,DEVICE-NAME  TO ENABLE GRAPHIC DISPLAY
 ENTER  FINISH            TO LEAVE PREP7
 PRINTOUT KEY SET TO /GOPR (USE /NOPR TO SUPPRESS)

 ELEMENT TYPE      1 IS BEAM188      3-D 2-NODE BEAM
  KEYOPT( 1- 6)=        0      0      0        0      0      0
  KEYOPT( 7-12)=        0      0      0        0      0      0
  KEYOPT(13-18)=        0      0      0        0      0      0

 CURRENT NODAL DOF SET IS  UX    UY    UZ    ROTX  ROTY  ROTZ
  THREE-DIMENSIONAL MODEL

 *** PROPERTY TEMPERATURE TABLE    NUM. TEMPS=  1 ***
 SLOC=  1    0.000000

 PRINTOUT RESUMED BY /GOP

 PROPERTY TABLE EX    MAT=      1  NUM. POINTS=  1
 SLOC=  1    0.2100000E+12

 PROPERTY TABLE PRXY  MAT=      1  NUM. POINTS=  1
 SLOC=  1    0.3000000

 ALL CURRENT ANSYS DATA WRITTEN TO FILE NAME= Assignment5_try4.db
  FOR POSSIBLE RESUME FROM THIS POINT

 PRODUCE NODAL PLOT IN DSYS=  0

  INPUT SECTION ID NUMBER                1
  INPUT SECTION TYPE                BEAM
  INPUT BEAM SECTION SUBTYPE        Circular Solid
  INPUT BEAM SECTION NAME

  BEAM SECTION WITH SECTION ID NUMBER      1 IS OFFSET TO CENTROID

  SECTION ID NUMBER IS:        1
  BEAM SECTION TYPE IS:    Circular Solid
  BEAM SECTION NAME IS:
  COMPUTED BEAM SECTION DATA SUMMARY:
    Area                = 0.16001E-02
    Iyy                  = 0.20372E-06
    Iyz                  = 0.97194E-23
    Izz                  = 0.20372E-06
    Warping Constant    = 0.37199E-42
    Torsion Constant    = 0.40744E-06
    Centroid Y          = 0.41357E-18
    Centroid Z          =-0.15964E-17
    Shear Center Y      = 0.85597E-18
    Shear Center Z      = 0.33295E-18
    Shear Correction-yy  = 0.85698
    Shear Correction-yz  = 0.81254E-14
    Shear Correction-zz  = 0.85698

    Beam Section is offset to CENTROID of cross section

 KEYPOINT      1  X,Y,Z=  0.00000      0.00000      0.00000      IN CSYS=  0

 PLOT KEY POINTS FROM      1 TO      1 BY      1

 KEYPOINT      2  X,Y,Z=  10.0000      0.00000      0.00000      IN CSYS=  0

 PLOT KEY POINTS FROM      1 TO      2 BY      1

 KEYPOINT      3  X,Y,Z=  20.0000      0.00000      0.00000      IN CSYS=  0

 PLOT KEY POINTS FROM      1 TO      3 BY      1

 LINE CONNECTS KEYPOINTS      1    2
  LINE NO.=    1  KP1=      1  TAN1=  -1.0000  0.0000  0.0000
                  KP2=      2  TAN2=    1.0000  0.0000  0.0000

 LINE CONNECTS KEYPOINTS      2    3
  LINE NO.=    2  KP1=      2  TAN1=  -1.0000  0.0000  0.0000
                  KP2=      3  TAN2=    1.0000  0.0000  0.0000

 PLOT KEY POINTS FROM      1 TO      3 BY      1

 SET DIVISIONS ON ALL PICKED LINES
      TO  NDIV =  10,  SPACING RATIO =  1.000
 NDIV AND SPACE CAN BE OVERRIDDEN FOR CURVATURE OR PROXIMITY

 PLOT KEY POINTS FROM      1 TO      3 BY      1

 GENERATE NODES AND ELEMENTS  ON  ALL  PICKED  LINES

 NUMBER OF LINES MESHED    =        2
 MAXIMUM NODE NUMBER        =        21
 MAXIMUM ELEMENT NUMBER    =        20

 PRODUCE ELEMENT PLOT IN DSYS =  0

 PRODUCE ELEMENT PLOT IN DSYS =  0


 ***** ROUTINE COMPLETED *****  CP =        42.859



 *****  ANSYS SOLUTION ROUTINE  *****


 PERFORM A STATIC ANALYSIS
  THIS WILL BE A NEW ANALYSIS

 KP  NUMBERING KEY =  0

 LINE NUMBERING KEY =  0

 AREA NUMBERING KEY =  0

 VOLU NUMBERING KEY =  0

 NODE NUMBERING KEY =  1

 TABN NUMBERING KEY =  0

 SVAL NUMBERING KEY =  0

 DOMA NUMBERING KEY =  0

 NUMBER KEY SET TO  0  -1=NONE  0=BOTH  1=COLOR  2=NUMBER

 ELEM NUMBERING KEY =  0

 KP  NUMBERING KEY =  1

 LINE NUMBERING KEY =  0

 AREA NUMBERING KEY =  0

 VOLU NUMBERING KEY =  0

 NODE NUMBERING KEY =  1

 TABN NUMBERING KEY =  0

 SVAL NUMBERING KEY =  0

 DOMA NUMBERING KEY =  0

 NUMBER KEY SET TO  0  -1=NONE  0=BOTH  1=COLOR  2=NUMBER

 ELEM NUMBERING KEY =  0

 KP  NUMBERING KEY =  1

 LINE NUMBERING KEY =  0

 AREA NUMBERING KEY =  0

 VOLU NUMBERING KEY =  0

 NODE NUMBERING KEY =  1

 TABN NUMBERING KEY =  0

 SVAL NUMBERING KEY =  0

 DOMA NUMBERING KEY =  0

 NUMBER KEY SET TO  0  -1=NONE  0=BOTH  1=COLOR  2=NUMBER

 ELEM NUMBERING KEY =  1

 KP  NUMBERING KEY =  1

 LINE NUMBERING KEY =  0

 AREA NUMBERING KEY =  0

 VOLU NUMBERING KEY =  0

 NODE NUMBERING KEY =  1

 TABN NUMBERING KEY =  0

 SVAL NUMBERING KEY =  0

 DOMA NUMBERING KEY =  0

 NUMBER KEY SET TO  0  -1=NONE  0=BOTH  1=COLOR  2=NUMBER

 ELEM NUMBERING KEY =  0

 DISTANCE MULTIPLIER FOR WINDOW 1  0.7290

 DISTANCE FOR WINDOW 1  6.016

 PRODUCE ELEMENT PLOT IN DSYS =  0

 DISTANCE MULTIPLIER FOR WINDOW 1  1.372

 DISTANCE FOR WINDOW 1  8.252

 PRODUCE ELEMENT PLOT IN DSYS =  0


 SPECIFIED CONSTRAINT UX  FOR PICKED NODES
 REAL=  0.00000000      IMAG=  0.00000000
 ADDITIONAL DOFS=  UY

 ALL BOUNDARY CONDITION DISPLAY KEYS SET TO  1


 SPECIFIED CONSTRAINT UX  FOR PICKED NODES
 REAL=  0.00000000      IMAG=  0.00000000
 ADDITIONAL DOFS=  UY    UZ    ROTX  ROTY  ROTZ

 SPECIFIED SURFACE LOAD PRES FOR ALL PICKED BEAM ELEMENTS  LKEY = 2
    VALI,VALJ    =    -4000.0        -4000.0

 PRES LOAD SURFACE DISPLAY KEY =  1


 POINT LOAD ON ALL PICKED KEYPOINTS  LOAD LABEL = FY
      VALUES =    -15000.        0.0000

Here's my nodal solution:

Code:

PRINT U    NODAL SOLUTION PER NODE
 
  ***** POST1 NODAL DEGREE OF FREEDOM LISTING *****                           
 
  LOAD STEP=    0  SUBSTEP=    1                                           
  TIME=    1.0000      LOAD CASE=  0                                       
 
  THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE SYSTEM 
 
    NODE      UX          UY          UZ          USUM 
      1  0.0000      29.352      0.0000      29.352   
      2  0.0000      0.0000      0.0000      0.0000   
      3  0.0000      28.247      0.0000      28.247   
      4  0.0000      26.860      0.0000      26.860   
      5  0.0000      24.980      0.0000      24.980   
      6  0.0000      22.492      0.0000      22.492   
      7  0.0000      19.370      0.0000      19.370   
      8  0.0000      15.686      0.0000      15.686   
      9  0.0000      11.603      0.0000      11.603   
      10  0.0000      7.3786      0.0000      7.3786   
      11  0.0000      3.3628      0.0000      3.3628   
      12  0.0000      0.0000      0.0000      0.0000   
      13  0.0000      -2.3608      0.0000      2.3608   
      14  0.0000      -3.7286      0.0000      3.7286   
      15  0.0000      -4.2790      0.0000      4.2790   
      16  0.0000      -4.1876      0.0000      4.1876   
      17  0.0000      -3.6302      0.0000      3.6302   
      18  0.0000      -2.7823      0.0000      2.7823   
      19  0.0000      -1.8197      0.0000      1.8197   
      20  0.0000    -0.91800      0.0000      0.91800   
      21  0.0000    -0.25288      0.0000      0.25288   

 MAXIMUM ABSOLUTE VALUES
 NODE          0            1            0            1
 VALUE    0.0000      29.352      0.0000      29.352


Can someone please tell me what I'm doing wrong?
Thank You in advance!


All times are GMT -4. The time now is 20:16.