General Mechanical

General Mechanical

cyclic loading

    • santrondela
      Subscriber

       can anybody help in the application of the cyclic loading? I'm trying to analyze dynamic effect of the structure WRT time.


       

    • Wenlong
      Ansys Employee

      Hi Santrondela,


      You can create in excel a table, with the first column your time and second column your load amplitude, then paste that table into the Load amplitude in Ansys Mechanical.


      Regards,


      Wenlong

    • santrondela
      Subscriber

      Thanks Wenlong.


      there is one more thing in which i need your help. which is when I'm applying gravity load and running for the solutions it take so long time.


      If you'll help help in it. it'll very helpful please reply fast. because I'm running out of time.


       


       

    • peteroznewman
      Subscriber

      You will have to provide a lot more detail about your model if you want help making it solve in less time.


      One way to make a Transient Structural model solve in less time is if you can accept a linear analysis. This means no nonlinear contact, no nonlinear material and no large deflections.  If you can accept those limitations, then you can use a Modal Superposition (MSUP) Transient Structural model. That will solve in a small fraction of the time it takes to solve a full Transient Structural model.

    • santrondela
      Subscriber

      Thank you  Mr. peteroznewman


      but thing is that I'm performing a non-linear study so I've to go through it anyhow.


      And sure whatever you need I'll provide you regard my geometry.


      first thing i want to mention that, it is imported from solidworks.


       

    • santrondela
      Subscriber

    • santrondela
      Subscriber

      here is the some info regard my geometry & I'm conducting a time-history analysis

    • Wenlong
      Ansys Employee

      Hi santrondela,


       


      So it is an RC structure under earthquake loading? What is the total number of your elements and how long does it take to solve? Do you have a damage model in your concrete material? How do the rebar and concrete interact with each other? Considering the geometric complexity and the nonlinearity, it could make sense it takes a long time to run. 


      Regards,


      Wenlong


       

    • santrondela
      Subscriber

      Thank you Wenlong for reply.


      yes it is a RC structure under eq loading & elements are 5087195, and it doesn't solve yet.


      I've not yet applied damage model but yes i do have one for now I'm applying standard gravity load and one side fixed and one end is friction-less and column upper face is under axial loading.


      after some changes I've run the meshing for the part it shows the warning "One or more surfaces or surface boundaries might not be oriented correctly."


      and after solving this error


      "At least one contact pair has no elements in it.  This may be due to mesh defeaturing.  Please modify defeaturing settings which are accessible on the Mesh object.  Alternatively, set the variable "contactAllowEmpty" to 1 in order to allow the solution to proceed with a warning which can be used to identify the offending pair(s)."


      i didn't understood why it is so?


    • Wenlong
      Ansys Employee

      Hi Sant,


      Wow, you have 5 million elements LOL! I won't be surprised at all that it takes a long time to run. Your concrete element mesh looks quite coarse, so I guess the mesh around the rebars is really small. Did you use solid elements or beam elements for the rebar? If using beam elements the contact between the rebar and concrete is a bit tricky (Please refer to this post: https://studentcommunity.ansys.com/thread/error-contact/), but it can probably make your element number drop down significantly. 


       


      Regards,


      Wenlong


       



      Useful Links



       

    • santrondela
      Subscriber

      Hello Wenlong


      I've reduced it up to 1076255 but still I'm facing the same problem.


      And about modal I've had done modeling  in solidworks.

    • Wenlong
      Ansys Employee

      Hi Sant,


      Oh, I am not asking what software you use to create the geometry, but what element type you use for the rebars, solid or beam? This is just out of curiosity as I'd like to know how do you bond the rebar and the concrete together because for beam elements, it would be a bit difficult. 


      For the mesh,  I would usually do a metrics check to make sure my mesh quality is good. Again may not be directly related to your problem but it is a practice I would do. You can find more information at: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v201/en/wb_msh/msh_metrics.html?q=mesh%20metrics . It can help you get an idea of whether you have some geometric defects (such as small faces or small edges) that caused an extremely fine mesh in some regions. 


      If your geometry do have some small faces or edges, you may use the "repair" tool in Spaceclaim to fix it before importing to Mechanical: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v201/en/spaceclaim/Discovery/user_manual/repair_overview.html?q=spaceclaim%20geometry%20repair


      Regards,


      Wenlong


       



      Useful Links



       


       


       


       

    • santrondela
      Subscriber

      An internal solution magnitude limit was exceeded. Please check your Environment for inappropriate load values or insufficient supports.  Also check that your mesh has more than 1 element in at least 2 directions if solid brick elements are present.  Please see the Troubleshooting section of the Help System for more information.

    • santrondela
      Subscriber

      Hello Wenlong


      actually i didn't use any command to connect both of them. what should i do about that?


      is there anyway?

    • santrondela
      Subscriber

      *


      *
      | |
      | W E L C O M E T O T H E A N S Y S (R) P R O G R A M |
      | |
      *
      *
      ***************************************************************
      * ANSYS Release 16.0 LEGAL NOTICES *
      ***************************************************************
      * *
      * Copyright 2014 SAS IP, Inc. All rights reserved. *
      * Unauthorized use, distribution or duplication is *
      * prohibited. *
      * *
      * Ansys is a registered trademark of ANSYS, Inc. or its *
      * subsidiaries in the United States or other countries. *
      * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
      * documentation CD or online help for the complete Legal *
      * Notice. *
      * *
      ***************************************************************
      * *
      * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION *
      * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY *
      * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. *
      * The software products and documentation are furnished by *
      * ANSYS, Inc. or its subsidiaries under a software license *
      * agreement that contains provisions concerning *
      * non-disclosure, copying, length and nature of use, *
      * compliance with exporting laws, warranties, disclaimers, *
      * limitations of liability, and remedies, and other *
      * provisions. The software products and documentation may be *
      * used, disclosed, transferred, or copied only in accordance *
      * with the terms and conditions of that software license *
      * agreement. *
      * *
      * ANSYS, Inc. is a UL registered *
      * ISO 9001:2008 company. *
      * *
      ***************************************************************
      * *
      * This product is subject to U.S. laws governing export and *
      * re-export. *
      * *
      * For U.S. Government users, except as specifically granted *
      * by the ANSYS, Inc. software license agreement, the use, *
      * duplication, or disclosure by the United States Government *
      * is subject to restrictions stated in the ANSYS, Inc. *
      * software license agreement and FAR 12.212 (for non-DOD *
      * licenses). *
      * *
      ***************************************************************
      Release 16.0
      Point Releases and Patches installed:
      ANSYS, Inc. Products Release 16.0
      ANSYS Mechanical Products Release 16.0
      ANSYS Customization Files for User Programmable Features Release 16.0
      ANSYS Autodyn Release 16.0
      ANSYS LS-DYNA Release 16.0
      ANSYS AIM Release 16.0
      ANSYS CFX (includes ANSYS CFD-Post) Release 16.0
      ANSYS Fluent (includes ANSYS CFD-Post) Release 16.0
      ANSYS TurboGrid Release 16.0
      ANSYS Polyflow (includes ANSYS CFD-Post) Release 16.0
      ANSYS CFD-Post only Release 16.0
      ANSYS ICEM CFD Release 16.0
      ANSYS Aqwa Release 16.0
      ANSYS Composite PrepPost Release 16.0
      ANSYS Icepak (includes ANSYS CFD-Post) Release 16.0
      ANSYS Remote Solve Manager Standalone Services Release 16.0
      Catia, Version 6 Release 16.0
      ANSYS, Inc. License Manager Release 16.0
      ***** ANSYS COMMAND LINE ARGUMENTS *****
      BATCH MODE REQUESTED (-b) = NOLIST
      INPUT FILE COPY MODE (-c) = COPY
      SHARED MEMORY PARALLEL REQUESTED
      SINGLE PROCESS WITH 2 THREADS REQUESTED
      TOTAL OF 2 CORES REQUESTED
      START-UP FILE MODE = NOREAD
      STOP FILE MODE = NOREAD
      RELEASE= Release 16.0 BUILD= 16.0 UP20141203 VERSION=WINDOWS x64
      CURRENT JOBNAME=file 11:51:08 MAR 31, 2020 CP= 1.469
      PARAMETER _DS_PROGRESS = 999.0000000
      /INPUT FILE= ds.dat LINE= 0
      DO NOT WRITE ELEMENT RESULTS INTO DATABASE
      *GET _WALLSTRT FROM ACTI ITEM=TIME WALL VALUE= 11.8525000
      TITLE=
      Solver

    • santrondela
      Subscriber

      without wire--Static Structural (A5)
      SET PARAMETER DIMENSIONS ON _WB_PROJECTSCRATCH_DIR
      TYPE=STRI DIMENSIONS= 248 1 1
      PARAMETER _WB_PROJECTSCRATCH_DIR(1) = G:thesis solidansys26.3.20203.57_ProjectScratchScr20AA
      SET PARAMETER DIMENSIONS ON _WB_SOLVERFILES_DIR
      TYPE=STRI DIMENSIONS= 248 1 1
      PARAMETER _WB_SOLVERFILES_DIR(1) = G:thesis solidansys26.3.20203.57without wire_filesdp0SYSMECH
      SET PARAMETER DIMENSIONS ON _WB_USERFILES_DIR
      TYPE=STRI DIMENSIONS= 248 1 1
      PARAMETER _WB_USERFILES_DIR(1) = G:thesis solidansys26.3.20203.57without wire_filesuser_files
      --- Data in consistent NMM units. See Solving Units in the help system for more
      MPA UNITS SPECIFIED FOR INTERNAL
      LENGTH = MILLIMETERS (mm)
      MASS = TONNE (Mg)
      TIME = SECONDS (sec)
      TEMPERATURE = CELSIUS (C)
      TOFFSET = 273.0
      FORCE = NEWTON (N)
      HEAT = MILLIJOULES (mJ)
      INPUT UNITS ARE ALSO SET TO MPA
      *** ANSYS - ENGINEERING ANALYSIS SYSTEM RELEASE Release 16.0 16.0 ***
      ANSYS Multiphysics
      00000000 VERSION=WINDOWS x64 11:51:09 MAR 31, 2020 CP= 1.609
      without wire--Static Structural (A5)
      ***** ANSYS ANALYSIS DEFINITION (PREP7) *****
      *********** Nodes for the whole assembly ***********
      *********** Elements for Body 1 "Solid" ***********
      *********** Elements for Body 2 "Solid" ***********
      *********** Elements for Body 3 "Solid" ***********
      *********** Elements for Body 4 "Solid" ***********
      *********** Elements for Body 5 "Solid" ***********
      *********** Elements for Body 6 "Solid" ***********
      *********** Elements for Body 7 "Solid" ***********
      *********** Elements for Body 8 "Solid" ***********
      *********** Elements for Body 9 "Solid" ***********
      *********** Elements for Body 10 "Solid" ***********
      *********** Elements for Body 11 "Solid" ***********
      *********** Elements for Body 12 "Solid" ***********
      *********** Elements for Body 13 "Solid" ***********
      *********** Elements for Body 14 "Solid" ***********
      *********** Elements for Body 15 "Solid" ***********
      *********** Elements for Body 16 "Solid" ***********
      *********** Elements for Body 17 "Solid" ***********
      *********** Elements for Body 18 "Solid" ***********
      *********** Elements for Body 19 "Solid" ***********
      *********** Elements for Body 20 "Solid" ***********
      *********** Elements for Body 21 "Solid" ***********
      *********** Send User Defined Coordinate System(s) ***********
      *********** Set Reference Temperature ***********
      *********** Send Materials ***********
      *********** Create Contact "Bonded - Multiple To Solid" ***********
      Real Constant Set For Above Contact Is 23 & 22
      *********** Fixed Supports ***********
      ******* Constant Zero Displacement X *******
      ******* Constant Zero Displacement Z *******
      *********** Define Pressure Using Surface Effect Elements ***********
      LIST DATA TABLE DMGE FOR MATERIAL 20
      *** ANSYS - ENGINEERING ANALYSIS SYSTEM RELEASE Release 16.0 16.0 ***
      ANSYS Multiphysics
      00000000 VERSION=WINDOWS x64 11:518 MAR 31, 2020 CP= 35.906
      without wire--Static Structural (A5)
      DMG EVLU (DMGE) Table For Material 20
      Material Prop Degradation
      1
      Temps 0.0000000e+000
      FIBR TEN 7.5100000e-001
      FIBR CMP 5.0000000e-001
      MATX TEN 8.0000000e-001
      MATX CMP 5.0000000e-001


      0.0000000e+000

      0.0000000e+000
      ***** ROUTINE COMPLETED ***** CP = 35.969
      --- Number of total nodes = 1485655
      --- Number of contact elements = 1030253
      --- Number of spring elements = 0
      --- Number of bearing elements = 0
      --- Number of solid elements = 682036
      --- Number of total elements = 1712289
      *GET _WALLBSOL FROM ACTI ITEM=TIME WALL VALUE= 11.8605556
      ****************************************************************************
      ************************* SOLUTION ********************************
      ****************************************************************************
      ***** ANSYS SOLUTION ROUTINE *****

    • santrondela
      Subscriber

      PERFORM A STATIC ANALYSIS
      THIS WILL BE A NEW ANALYSIS
      LARGE DEFORMATION ANALYSIS
      USE PRECONDITIONED CONJUGATE GRADIENT SOLVER
      CONVERGENCE TOLERANCE = 1.00000E-08
      MAXIMUM ITERATION = NumNode*DofPerNode* 1.0000
      MEMORY SAVING OPTION TURNED ON FOR PCG SOLVER
      CONTACT INFORMATION PRINTOUT LEVEL 1
      NLDIAG: Nonlinear diagnostics CONT option is set to ON.
      Writing frequency : each ITERATION.
      DEFINE RESTART CONTROL FOR LOADSTEP LAST
      AT FREQUENCY OF LAST AND NUMBER FOR OVERWRITE IS 0
      DELETE RESTART FILES OF ENDSTEP
      ****************************************************
      ******************* SOLVE FOR LS 1 ****************
      SELECT FOR ITEM=TYPE COMPONENT=
      IN RANGE 24 TO 24 STEP 1
      2117 ELEMENTS (OF 1712289 DEFINED) SELECTED BY ESEL COMMAND.
      SELECT ALL NODES HAVING ANY ELEMENT IN ELEMENT SET.
      4814 NODES (OF 1485655 DEFINED) SELECTED FROM
      2117 SELECTED ELEMENTS BY NSLE COMMAND.
      GENERATE SURFACE LOAD PRES ON SURFACE DEFINED BY ALL SELECTED NODES
      SET ACCORDING TO TABLE PARAMETER = _LOADVARI18224
      NUMBER OF PRES ELEMENT FACE LOADS STORED = 2117
      ALL SELECT FOR ITEM=NODE COMPONENT=
      IN RANGE 1 TO 1485655 STEP 1
      1485655 NODES (OF 1485655 DEFINED) SELECTED BY NSEL COMMAND.
      ALL SELECT FOR ITEM=ELEM COMPONENT=
      IN RANGE 1 TO 1712289 STEP 1
      1712289 ELEMENTS (OF 1712289 DEFINED) SELECTED BY ESEL COMMAND.
      *********** Create Acceleration ***********
      SET PARAMETER DIMENSIONS ON _ACELX TYPE=TABL DIMENSIONS= 2 1 1
      TIME
      PARAMETER _ACELX(1,0,1) = 0.000000000
      PARAMETER _ACELX(2,0,1) = 1.000000000
      PARAMETER _ACELX(1,1,1) = 0.000000000
      PARAMETER _ACELX(2,1,1) = 0.000000000
      SET PARAMETER DIMENSIONS ON _ACELY TYPE=TABL DIMENSIONS= 2 1 1
      TIME
      PARAMETER _ACELY(1,0,1) = 0.000000000
      PARAMETER _ACELY(2,0,1) = 1.000000000
      PARAMETER _ACELY(1,1,1) = 0.000000000
      PARAMETER _ACELY(2,1,1) = 9806.650000
      SET PARAMETER DIMENSIONS ON _ACELZ TYPE=TABL DIMENSIONS= 2 1 1
      TIME
      PARAMETER _ACELZ(1,0,1) = 0.000000000
      PARAMETER _ACELZ(2,0,1) = 1.000000000
      PARAMETER _ACELZ(1,1,1) = 0.000000000
      PARAMETER _ACELZ(2,1,1) = 0.000000000
      ACELX IS SET ACCORDING TO TABLE PARAMETER = _ACELX
      ACELY IS SET ACCORDING TO TABLE PARAMETER = _ACELY
      ACELZ IS SET ACCORDING TO TABLE PARAMETER = _ACELZ
      ALL SELECT FOR ITEM=ELEM COMPONENT=
      IN RANGE 1 TO 1712289 STEP 1
      1712289 ELEMENTS (OF 1712289 DEFINED) SELECTED BY ESEL COMMAND.
      PRINTOUT RESUMED BY /GOP
      USE AUTOMATIC TIME STEPPING THIS LOAD STEP
      USE 5 SUBSTEPS INITIALLY THIS LOAD STEP FOR ALL DEGREES OF FREEDOM
      FOR AUTOMATIC TIME STEPPING:
      USE 1000 SUBSTEPS AS A MAXIMUM
      USE 1 SUBSTEPS AS A MINIMUM
      TIME= 1.0000
      ERASE THE CURRENT DATABASE OUTPUT CONTROL TABLE.
      WRITE ALL ITEMS TO THE DATABASE WITH A FREQUENCY OF NONE
      FOR ALL APPLICABLE ENTITIES
      WRITE NSOL ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL

    • santrondela
      Subscriber

      FOR ALL APPLICABLE ENTITIES
      WRITE RSOL ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
      FOR ALL APPLICABLE ENTITIES
      WRITE STRS ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
      FOR ALL APPLICABLE ENTITIES
      WRITE EPEL ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
      FOR ALL APPLICABLE ENTITIES
      WRITE EPPL ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
      FOR ALL APPLICABLE ENTITIES
      WRITE AESO ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
      FOR ALL APPLICABLE ENTITIES
      NONLINEAR STABILIZATION CONTROL:
      KEY=OFF
      KEY TO TERMINATE RUN IF NO CONVERGENCE= 1
      DEGREE OF FREEDOM (DISPLACEMENT) LIMIT= 0.78886E-30
      CUMULATIVE ITERATION LIMIT= 0
      ELAPSED TIME LIMIT (SEC)= 0.0000 ( 0 HOURS 0 MIN)
      CPU TIME LIMIT (SEC)= 0.0000 ( 0 HOURS 0 MIN)
      *GET ANSINTER_ FROM ACTI ITEM=INT VALUE= 0.00000000
      *IF ANSINTER_ ( = 0.00000 ) NE
      0 ( = 0.00000 ) THEN
      *ENDIF
      ***** ANSYS SOLVE COMMAND *****
      *** WARNING *** CP = 42.953 TIME= 11:51:45
      Element shape checking is currently inactive. Issue SHPP,ON or
      SHPP,WARN to reactivate, if desired.
      *** NOTE *** CP = 57.188 TIME= 11:51:55
      The model data was checked and warning messages were found.
      Please review output or errors file ( G:thesis
      solidansys26.3.20203.57_ProjectScratchScr20AAfile.err ) for
      these warning messages.
      *** SELECTION OF ELEMENT TECHNOLOGIES FOR APPLICABLE ELEMENTS ***
      --- GIVE SUGGESTIONS AND RESET THE KEY OPTIONS ---
      ELEMENT TYPE 1 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 2 IS SOLID186. KEYOPT(2)=0 IS SUGGESTED AND HAS BEEN RESET.
      KEYOPT(1-12)= 0 0 0 0 0 0 0 0 0 0 0 0
      ELEMENT TYPE 3 IS SOLID186. KEYOPT(2)=0 IS SUGGESTED AND HAS BEEN RESET.
      KEYOPT(1-12)= 0 0 0 0 0 0 0 0 0 0 0 0
      ELEMENT TYPE 4 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 5 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 6 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 7 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 8 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 9 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 10 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 11 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 12 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 13 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 14 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 15 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 16 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 17 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 18 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 19 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 20 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.
      ELEMENT TYPE 21 IS SOLID187. IT IS NOT ASSOCIATED WITH FULLY INCOMPRESSIBLE
      HYPERELASTIC MATERIALS. NO SUGGESTION IS AVAILABLE AND NO RESETTING IS NEEDED.

    • santrondela
      Subscriber

      *** ANSYS - ENGINEERING ANALYSIS SYSTEM RELEASE Release 16.0 16.0 ***
      ANSYS Multiphysics
      00000000 VERSION=WINDOWS x64 11:51:56 MAR 31, 2020 CP= 57.922
      without wire--Static Structural (A5)
      S O L U T I O N O P T I O N S
      PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D
      DEGREES OF FREEDOM. . . . . . UX UY UZ
      ANALYSIS TYPE . . . . . . . . . . . . . . . . .STATIC (STEADY-STATE)
      OFFSET TEMPERATURE FROM ABSOLUTE ZERO . . . . . 273.15
      NONLINEAR GEOMETRIC EFFECTS . . . . . . . . . .ON
      EQUATION SOLVER OPTION. . . . . . . . . . . . .PCG
      MEMORY SAVING OPTION. . . . . . . . . . . .ON
      TOLERANCE. . . . . . . . . . . . . . . . . . 1.00000E-08
      PLASTIC MATERIAL PROPERTIES INCLUDED. . . . . .YES
      NEWTON-RAPHSON OPTION . . . . . . . . . . . . .PROGRAM CHOSEN
      GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC
      *** WARNING *** CP = 59.734 TIME= 11:51:58
      Material number 23 (used by element 682037 ) should normally have at
      least one MP or one TB type command associated with it. Output of
      energy by material may not be available.
      *** NOTE *** CP = 66.375 TIME= 11:52:06
      The step data was checked and warning messages were found.
      Please review output or errors file ( G:thesis
      solidansys26.3.20203.57_ProjectScratchScr20AAfile.err ) for
      these warning messages.
      *** NOTE *** CP = 66.375 TIME= 11:52:06
      This nonlinear analysis defaults to using the full Newton-Raphson
      solution procedure. This can be modified using the NROPT command.
      *** NOTE *** CP = 66.375 TIME= 11:52:06
      The conditions for direct assembly have been met. No .emat or .erot
      files will be produced.
      L O A D S T E P O P T I O N S
      LOAD STEP NUMBER. . . . . . . . . . . . . . . . 1
      TIME AT END OF THE LOAD STEP. . . . . . . . . . 1.0000
      AUTOMATIC TIME STEPPING . . . . . . . . . . . . ON
      INITIAL NUMBER OF SUBSTEPS . . . . . . . . . 5
      MAXIMUM NUMBER OF SUBSTEPS . . . . . . . . . 1000
      MINIMUM NUMBER OF SUBSTEPS . . . . . . . . . 1
      MAXIMUM NUMBER OF EQUILIBRIUM ITERATIONS. . . . 15
      STEP CHANGE BOUNDARY CONDITIONS . . . . . . . . NO
      STRESS-STIFFENING . . . . . . . . . . . . . . . ON
      TERMINATE ANALYSIS IF NOT CONVERGED . . . . . .YES (EXIT)
      NODAL DOF SOLUTION LIMIT. . . . . . . . . . . . 0.78886E-30
      CONVERGENCE CONTROLS. . . . . . . . . . . . . .USE DEFAULTS
      INERTIA LOADS X Y Z
      ACEL . . . . . . . . . . . ._ACELX _ACELY _ACELZ
      COPY INTEGRATION POINT VALUES TO NODE . . . . .YES, FOR ELEMENTS WITH
      ACTIVE MAT. NONLINEARITIES
      PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
      DATABASE OUTPUT CONTROLS
      ITEM FREQUENCY COMPONENT
      ALL NONE
      NSOL ALL
      RSOL ALL
      STRS ALL
      EPEL ALL
      EPPL ALL
      AESO ALL
      SOLUTION MONITORING INFO IS WRITTEN TO FILE=
      file.mntr
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      It is highly recommended to use the auto contact setting option by
      issuing CNCHECK,AUTO command for this problem in order to achieve
      better convergence.
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      Symmetric Deformable- deformable contact pair identified by real
      constant set 22 and contact element type 22 has been set up. The
      companion pair has real constant set ID 23. Both pairs should have
      the same behavior.
      ANSYS will deactivate the current pair and keep its companion pair,
      resulting in asymmetric contact.
      Contact algorithm: Augmented Lagrange method
      Contact detection at: Gauss integration point
      Contact stiffness factor FKN 10.000
      The resulting initial contact stiffness 0.14457E+06
      Default penetration tolerance factor FTOLN 0.10000
      The resulting penetration tolerance 0.49471
      Default opening contact stiffness OPSF will be used.
      Default tangent stiffness factor FKT 1.0000
      Default elastic slip factor SLTOL 0.50000E-02
      The resulting elastic slip 0.45970E-01
      Update contact stiffness at each iteration
      Default Max. friction stress TAUMAX 0.10000E+21
      Average contact surface length 9.1939
      Average contact pair depth 4.9471
      Default pinball region factor PINB 0.50000
      The resulting pinball region 2.4736
      Initial penetration/gap is excluded.
      Bonded contact (always) is defined.
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      No contact was detected for this contact pair.

    • santrondela
      Subscriber

      ****************************************
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      Symmetric Deformable- deformable contact pair identified by real
      constant set 23 and contact element type 22 has been set up. The
      companion pair has real constant set ID 22. Both pairs should have
      the same behavior.
      ANSYS will keep the current pair and deactivate its companion pair,
      resulting in asymmetric contact.
      Contact algorithm: Augmented Lagrange method
      Contact detection at: Gauss integration point
      Contact stiffness factor FKN 10.000
      The resulting initial contact stiffness 0.14457E+06
      Default penetration tolerance factor FTOLN 0.10000
      The resulting penetration tolerance 4.1502
      Default opening contact stiffness OPSF will be used.
      Default tangent stiffness factor FKT 1.0000
      Default elastic slip factor SLTOL 0.50000E-02
      The resulting elastic slip 0.27899
      Update contact stiffness at each iteration
      Default Max. friction stress TAUMAX 0.10000E+21
      Average contact surface length 55.797
      Average contact pair depth 41.502
      Default pinball region factor PINB 0.50000
      The resulting pinball region 20.751
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      One of the contact searching regions contains at least 7814 target
      elements. You may reduce the pinball radius.
      Initial penetration/gap is excluded.
      Bonded contact (always) is defined.
      *** NOTE *** CP = 742.922 TIME= 12:03:28
      Max. Initial penetration 160.848824 was detected between contact
      element 1703168 and target element 716308.
      You may move entire target surface by : x= 160.848824, y= 0, z= 0,to
      reduce initial penetration.
      *WARNING*: The detected penetration is larger than the pair based
      pinball. Please verify the pinball radius carefully.
      *WARNING*: The geometric gap/penetration may be too large. Increase
      pinball radius if it is a true geometric gap/penetration. Decrease
      pinball if it is a false one.
      *WARNING*: The initial penetration/gap is relatively large. Bonded/no
      separation option may cause an accuracy issue. Switch to MPC
      algorithm or you may use the CNCHECK,ADJUST command to move the
      contact nodes towards the target surface.
      ****************************************
      MAXIMUM NUMBER OF EQUILIBRIUM ITERATIONS HAS BEEN MODIFIED
      TO BE, NEQIT = 26, BY SOLUTION CONTROL LOGIC.
      *** NOTE *** CP = 760.297 TIME= 12:03:47
      The PCG solver has automatically set the level of difficulty for this
      model to 2.
      *************************************************
      SUMMARY FOR CONTACT PAIR IDENTIFIED BY REAL CONSTANT SET 22
      *** NOTE *** CP = 1020.969 TIME= 12:07:24
      Contact pair is inactive.
      *** NOTE *** CP = 1021.016 TIME= 12:07:24
      One of the contact searching regions contains at least 384 target
      elements. You may reduce the pinball radius (current value
      20.7507507) for contact pair identified by real constant set 23 to
      speed up contact searching.
      *********** PRECISE MASS SUMMARY ***********
      TOTAL RIGID BODY MASS MATRIX ABOUT ORIGIN
      Translational mass | Coupled translational/rotational mass
      6.8789 0.0000 0.0000 | 0.0000 2233.7 -15.061
      0.0000 6.8789 0.0000 | -2233.7 0.0000 5888.2
      0.0000 0.0000 6.8789 | 15.061 -5888.2 0.0000


      |

      | Rotational mass (inertia)
      | 0.36483E+07 -10322. -0.19125E+07
      | -10322. 0.12946E+08 -4844.8
      | -0.19125E+07 -4844.8 0.14656E+08
      TOTAL MASS = 6.8789
      The mass principal axes coincide with the global Cartesian axes
      CENTER OF MASS (X,Y,Z)= 855.99 2.1894 324.73
      TOTAL INERTIA ABOUT CENTER OF MASS
      0.29229E+07 2569.9 -429.81
      2569.9 0.71801E+07 45.811
      -429.81 45.811 0.96154E+07
      The inertia principal axes coincide with the global Cartesian axes
      *** MASS SUMMARY BY ELEMENT TYPE ***
      TYPE MASS
      1 0.315569
      2 0.115511E-01
      3 0.115511E-01
      4 0.150892
      5 0.248904E-01
      6 0.835647E-02
      7 0.835166E-02
      8 0.827129E-02
      9 0.827458E-02
      10 0.830644E-02

    • santrondela
      Subscriber

      11 0.835362E-02
      12 0.121986E-03
      13 0.121216E-03
      14 0.121862E-03
      15 0.121892E-03
      16 0.121814E-03
      17 0.121239E-03
      18 0.121878E-03
      19 0.121988E-03
      20 0.122401E-03
      21 6.31341
      Range of element maximum matrix coefficients in global coordinates
      Maximum = 526768871 at element 562821.
      Minimum = 117109.031 at element 1704882.
      *** ELEMENT MATRIX FORMULATION TIMES
      TYPE NUMBER ENAME TOTAL CP AVE CP
      1 342575 SOLID187 152.000 0.000444
      2 1540 SOLID186 2.297 0.001491
      3 1540 SOLID186 2.031 0.001319
      4 199080 SOLID187 93.984 0.000472
      5 42619 SOLID187 21.281 0.000499
      6 1737 SOLID187 0.672 0.000387
      7 1834 SOLID187 0.984 0.000537
      8 1700 SOLID187 1.000 0.000588
      9 1760 SOLID187 0.719 0.000408
      10 1724 SOLID187 0.641 0.000372
      11 1721 SOLID187 0.797 0.000463
      12 8539 SOLID187 3.266 0.000382
      13 6827 SOLID187 2.875 0.000421
      14 8175 SOLID187 3.406 0.000417
      15 8196 SOLID187 3.250 0.000397
      16 8256 SOLID187 3.344 0.000405
      17 6793 SOLID187 2.297 0.000338
      18 8203 SOLID187 3.234 0.000394
      19 8224 SOLID187 5.094 0.000619
      20 9386 SOLID187 4.125 0.000439
      21 11607 SOLID187 5.609 0.000483
      22 514068 CONTA174 91.656 0.000178
      23 514068 TARGE170 43.453 0.000085
      24 2117 SURF154 0.531 0.000251
      Time at end of element matrix formulation CP = 1021.03125.
      ALL CURRENT ANSYS DATA WRITTEN TO FILE NAME= file.rdb
      FOR POSSIBLE RESUME FROM THIS POINT
      FORCE CONVERGENCE VALUE = 0.8532E+06 CRITERION= 4353.
      curEqn= 420828 totEqn= 420828 Job CP sec= 1460.685
      Factor Done= 100% Factor Wall sec= 0.000 rate= 0.0 Mflops
      *** ERROR *** CP = 1070.719 TIME= 12:17:07
      Preconditioned conjugate gradient solver error level 1. Possibly, the
      model is unconstrained or additional iterations may be needed. Try
      running setting the multiplier MULT on the EQSLV command to greater
      than 1.0 (but less than 3.0).
      PRECONDITIONED SOLVER CP TIME = 251.812
      PRECONDITIONED SOLVER ELAPSED TIME = 251.812
      *** NOTE *** CP = 1071.172 TIME= 12:17:09
      During this loadstep the PCG iterative solver took more than 1000
      iterations to solve the system of equations. In the future it may be
      more efficient to choose a direct solver, such as the SPARSE solver,
      for this analysis.
      NUMBER OF WARNING MESSAGES ENCOUNTERED= 2
      NUMBER OF ERROR MESSAGES ENCOUNTERED= 1
      ***** PROBLEM TERMINATED BY INDICATED ERROR(S) OR BY END OF INPUT DATA *****
      *** WARNING *** CP = 1072.000 TIME= 12:17:21
      During this session the elapsed time exceeds the CPU time by 70%.
      Often this indicates either a lack of physical memory (RAM) required
      to efficiently handle this simulation or it indicates a particularly
      slow hard drive configuration. This simulation can be expected to run
      faster on identical hardware if additional RAM or a faster hard drive
      configuration is made available. For more details, please see the
      ANSYS Performance Guide which is part of the ANSYS Help system.
      +


      A N S Y S S T A T I S T I C S
      +
      Release: Release 16.0 Build: 16.0 Update: UP20141203 Platform: WINDOWS x64
      Date Run: 03/31/2020 Time: 12:17
      Windows Process ID: 11656
      Processor Model: Intel(R) Core(TM) i5-8250U CPU @ 1.60GHz
      Compiler: Intel(R) FORTRAN Compiler Version 14.0.0 (Build: 20140422)
      Intel(R) C/C++ Compiler Version 14.0.0 (Build: 20140422)
      Intel(R) Math Kernel Library Version 11.1.3 Product Build 20140917
      Total number of cores available : 8
      Number of physical cores available : 4
      Number of processes requested : 1
      Number of threads per process requested : 2
      Total number of cores requested : 2 (Shared Memory Parallel)
      GPU Acceleration: Not Requested
      Job Name: file
      Working Directory: G:thesis solidansys26.3.20203.57_ProjectScratchScr20AA
      Total CPU time for main thread : 926.9 seconds
      Total CPU time summed for all threads : 1072.0 seconds
      Elapsed time spent pre-processing model (/PREP7) : 28.8 seconds
      Elapsed time spent solution - preprocessing : 42.8 seconds

    • santrondela
      Subscriber

      Elapsed time spent computing solution : 1473.6 seconds
      Elapsed time spent solution - postprocessing : 0.0 seconds
      Elapsed time spent post-processing model (/POST1) : 0.0 seconds
      Equation solver computational rate : 122520.4 Mflops
      Maximum total memory used : 7786.0 MB
      Maximum total memory allocated : 8256.0 MB
      Maximum total memory available : 8 GB
      +


      E N D A N S Y S S T A T I S T I C S
      +
      *
      *
      | |
      | ANSYS RUN COMPLETED |
      | |
      |
      |
      | |
      | Ansys Release 16.0 Build 16.0 UP20141203 WINDOWS x64 |
      | |
      |
      |
      | |
      | Database Requested(-db) 1024 MB Scratch Memory Requested 1024 MB |
      | Maximum Database Used 2033 MB Maximum Scratch Memory Used 5753 MB |
      | |
      |
      |
      | |
      | CP Time (sec) = 1071.984 Time = 12:17:21 |
      | Elapsed Time (sec) = 1576.000 Date = 03/31/2020 |
      | |
      *
      *
      Page 8 of 8
      file:///

    • santrondela
      Subscriber

      Hello Wenlong


      These are the out put file can you please help me out in this because I'm running out of time. Please help me out!

    • Wenlong
      Ansys Employee

      Hi Sant,


      You may want to check this post: https://studentcommunity.ansys.com/thread/error-contact/. The model I created is very similar to yours. The difference is my rebars are modeled using beam elements, my mesh is more regular, and they are bonded through EREINF command. Using this method can greatly improve your mesh and thus reduce your simulation time. 


      Regards,


      Wenlong


       

    • vaibhavtaranekar
      Subscriber

      Hi Wenlong,


      Can you share your command to use EREINF? the ones in the link shared by you do not work and some error always comes. i want to join LINK180 with Solid185 but havn't been able to succesfully do so. Kindly help me out.

    • Wenlong
      Ansys Employee

      Hi 


      Please post the error message in another thread. Maybe titled as "EREINF error message", so that it can be better viewed by others. But here is the command I use: 


      finish
      /prep7

      !
      !
      ! Get the maximum element type id and section id
      !
      !
      *get, typid, elem, 0, typm
      *get, secid, secp, 0, num, max

      sectypeid = secid + 1 ! Create a new section id by shifting the max by 1
      etype_id = typid + 1 ! Create a new type id by shifting the max by 1

      !
      !
      ! Define the reinforcement section type and Mesh200 element type
      !
      !

      secarea = 50 ! Define a section area for the reinforcement

      sectype, sectypeid, reinf, discrete ! Define a discrete reinf section
      secdata, steel_mat_id, secarea, mesh ! Define the section properties

      et, etype_id, 200, 2 ! Define a new element type

      !
      !
      ! Generate Mesh200 elements
      !
      !

      cmsel, s, rebars, elem ! Select the name selection "rebars".
      ! The name selection has to be bodies instead of edges
      emodify, all, type, etype_id ! Change the element type to Mesh200
      emodify, all, secnum, sectypeid ! Change the section to discrete reinf

      !
      !
      ! Generate Reinf264 elements based on Mesh200 and solid elements
      !
      !

      cmsel, a, concrete, elem ! Add concrete to the selection
      EREINF ! Generate reinforcement elements

      esel,s,ename,,mesh200 ! Select the Mesh200 elements
      edele,all ! Delete the Mesh200 elements
      FINISH

      allsel,all
      /SOLU
      OUTRES,ALL,all

       


      Thanks,


      Regards,


      Wenlong


       

    • vaibhavtaranekar
      Subscriber

      Thanks for the update, i will try these and let you know soon. Do i have to define LINK180 elements normally as i do or these only work with beam elements?

    • Wenlong
      Ansys Employee

      It should work with both. I forgot to mention I have another command snippet added to the part.



      Regards,


      Wenlong


       



      Useful Links



       

    • santrondela
      Subscriber

      Thank you Wenlong


       

    • vaibhavtaranekar
      Subscriber

      It should work with both. I forgot to mention I have another command snippet added to the part.



      Regards,


      Wenlong



      Useful Links




       


      are you using default concrete and steel material from workbench or have you used command to set solid185 or any other material like that? can you share your file so that i can view the model?

    • Wenlong
      Ansys Employee

      ANSYS staff are not permitted to share or download models. I am using the default steel model for rebars and elastic material for concrete with solid 185. 


      Here is an example of a similar simulation where concrete is modelled using microplane damage model. This model is created completely using APDL commands: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v193/ans_tec/tecreinfconc.html and you can download the input.


      Regards,


      Wenlong


       



      Useful Links



       


       


       


       

    • vaibhavtaranekar
      Subscriber

      Thank you for your input, i was succesfully able to run the model using SOLID185 and REINF264. I have few queries regarding the REINF264 usage. 


      1. If there are multiple rebars having different sizes, then how to apply the command so that the secarea of REINF264 is as per the diameter of bar. 


      2. How to determine the stresses and strains in the REINF264 elements? Can we view the results in workbench instead of using APDL commands? 

    • santrondela
      Subscriber

      Hi Wenlong


      *** ERROR ***                           CP =      21.500   TIME= 14:21:21
       Element type 1 is not the same shape as LINK180.  Switching to a       
       different shape is not allowed while elements of type 1 exist.


      can You please help me out in this what should i do in this case


       

Viewing 32 reply threads
  • You must be logged in to reply to this topic.