cyclic loading

  • 156 Views
  • Last Post 10 April 2020
  • Topic Is Solved
santrondela posted this 16 March 2020

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

 

With regard Sant

Order By: Standard | Newest | Votes
Wenlong posted this 18 March 2020

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

  • Liked by
  • santrondela
santrondela posted this 23 March 2020

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.

 

 

With regard Sant

peteroznewman posted this 23 March 2020

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.

  • Liked by
  • santrondela
santrondela posted this 23 March 2020

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.

 

With regard Sant

santrondela posted this 23 March 2020

With regard Sant

santrondela posted this 23 March 2020

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

With regard Sant

Wenlong posted this 23 March 2020

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 posted this 23 March 2020

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?


 

With regard Sant

Wenlong posted this 25 March 2020

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 posted this 27 March 2020

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.

With regard Sant

Wenlong posted this 27 March 2020

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 posted this 28 March 2020

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.

With regard Sant

santrondela posted this 28 March 2020

Hello Wenlong

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

is there anyway?

With regard Sant

santrondela posted this 31 March 2020

*------------------------------------------------------------------*
| |
| 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

With regard Sant

santrondela posted this 31 March 2020

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 solid\ansys\26.3.2020\3.57\_ProjectScratch\Scr20AA\
SET PARAMETER DIMENSIONS ON _WB_SOLVERFILES_DIR
TYPE=STRI DIMENSIONS= 248 1 1
PARAMETER _WB_SOLVERFILES_DIR(1) = G:\thesis solid\ansys\26.3.2020\3.57\without wire_files\dp0\SYS\MECH\
SET PARAMETER DIMENSIONS ON _WB_USERFILES_DIR
TYPE=STRI DIMENSIONS= 248 1 1
PARAMETER _WB_USERFILES_DIR(1) = G:\thesis solid\ansys\26.3.2020\3.57\without wire_files\user_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 *****

With regard Sant

santrondela posted this 31 March 2020

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

With regard Sant

santrondela posted this 31 March 2020

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
solid\ansys\26.3.2020\3.57\_ProjectScratch\Scr20AA\file.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.

With regard Sant

santrondela posted this 31 March 2020

*** 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
solid\ansys\26.3.2020\3.57\_ProjectScratch\Scr20AA\file.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.

With regard Sant

santrondela posted this 31 March 2020

****************************************
*** 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

With regard Sant

santrondela posted this 31 March 2020

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 solid\ansys\26.3.2020\3.57\_ProjectScratch\Scr20AA
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

With regard Sant

santrondela posted this 31 March 2020

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:///

With regard Sant

santrondela posted this 31 March 2020

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!

With regard Sant

Wenlong posted this 31 March 2020

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 posted this 01 April 2020

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 posted this 01 April 2020

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 posted this 01 April 2020

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 posted this 01 April 2020

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

Regards,

Wenlong

 

------------------------- Useful Links ----------------------------

 

santrondela posted this 01 April 2020

Thank you Wenlong

 

With regard Sant

vaibhavtaranekar posted this 02 April 2020

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 posted this 02 April 2020

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 ----------------------------

 

 

 

 

Show More Posts
Close