Users:Form Finding/Form Finding with CARAT++

From Carat++ Public Wiki
< Users:Form Finding(Difference between revisions)
Jump to: navigation, search
(Example of an Input Block)
(Benchmark examples)
 
(5 intermediate revisions by 2 users not shown)
Line 44: Line 44:
 
|elements that have to be formfound (i.e. that have no material contribution, but only geometrical stiffness), identified through their property ID               
 
|elements that have to be formfound (i.e. that have no material contribution, but only geometrical stiffness), identified through their property ID               
 
|-
 
|-
|colspan="3" style="background:#efefef;"| Optional Parameters
+
|colspan="3" style="background:#efefef;"| Optional parameters for full URS
 +
|-
 +
!HOMOTOPY_CURVE
 +
|LD-CURVE ''int'' 
 +
|Defines the homotopy factor (has to be between 0 and 1) for the different form-finding steps in a separate curve. If not used, the homotopy factor is set to 0, i.e. the reference problem is solved.
 +
|-
 +
|colspan="3" style="background:#efefef;"| Optional parameters for stress adaption / distortion control
 
|-
 
|-
 
!STRESS_ADAPTION
 
!STRESS_ADAPTION
Line 57: Line 63:
 
|''real''
 
|''real''
 
|Relaxation of the element distortion (Values between 0 and 1 are valid)
 
|Relaxation of the element distortion (Values between 0 and 1 are valid)
 +
|-
 +
|colspan="3" style="background:#efefef;"| Optional parameters for restart file
 +
|-
 +
!ELASTIC_RESTART_METHOD
 +
|''string'' 
 +
|Defines the update of elastic members for the restart file (NONE = everything is updated and results set to 0, FORCE = geometry is updated, results set to 0 and inital stresses added, DISPLACEMENT = geometry of elastic members is not updated, initial displacement is added and nodes with different neighbouring elements are doubled)
 +
|-
 +
|colspan="3" style="background:#efefef;"| Optional parameters for eXtended URS
 +
|-
 +
!X_URS
 +
|''bool'' 
 +
|Defines whether or not the xURS shall be applied to the computation resp. to the Membrane1 properties; Default is FALSE, leading to normal URS
 +
|-                     
 +
!X_URS_STEP
 +
|''int''
 +
|Defines the form-finding step from which on xURS is applied
 
|}       
 
|}       
  
Line 63: Line 85:
 
=== Example of an Input Block ===
 
=== Example of an Input Block ===
  
The input block shown below defines a formfinding analysis with the ID 1. The analysis performs 10 formfinding steps and stops the computation after 100 equilibrium iteration steps within each formfinding step if no convergence is achieved. The equilibrium accuracy of the analysis is defined as 10<sup>-6</sup>.
+
The input block shown below defines a formfinding analysis with the IDs 1,2 and 8. The analysis performs 10 formfinding steps and stops the computation after 100 equilibrium iteration steps within each formfinding step if no convergence is achieved. The equilibrium accuracy of the analysis is defined as 10<sup>-6</sup>.
  
 
The (optional) distortion control is enabled and used for elements releated to the elemnt properties with ID 1. The maximum allowed element distortion is 1.1 and the relaxation of the element distortion control is set to 0.  
 
The (optional) distortion control is enabled and used for elements releated to the elemnt properties with ID 1. The maximum allowed element distortion is 1.1 and the relaxation of the element distortion control is set to 0.  
Line 78: Line 100:
 
   MAX_ITER_EQUILIBRIUM = 100
 
   MAX_ITER_EQUILIBRIUM = 100
 
   EQUILIBRIUM_ACCURACY = 1e-06
 
   EQUILIBRIUM_ACCURACY = 1e-06
 +
  FORMFINDING_ELEMENTS = PROP_ID 1,2,8
 
   STRESS_ADAPTION = PROP_ID 1
 
   STRESS_ADAPTION = PROP_ID 1
 
   LAMBDA_MAX=1.1
 
   LAMBDA_MAX=1.1
 
   MUE_RELAX=0.0
 
   MUE_RELAX=0.0
 
</pre>
 
</pre>
 +
 +
=== Benchmark examples ===
 +
* form-finding of a four-point sail: ..\examples\benchmark_examples\analyses\formfinding_membrane1_I\cbm_4_point_Fofi.txt
 +
* form-finding of a catenoid with xURS: ..\examples\benchmark_examples\analyses\cbm_formfinding_xURS_Membrane1\cbm_formfinding_xURS_Membrane1.dat
  
 
== References ==
 
== References ==
  
 
<references/>
 
<references/>

Latest revision as of 11:19, 8 June 2017

The following table shows the important commands to define a form finding analysis in CARAT++ and gives an explanation for each parameter.

Compulsory Parameters
Parameter Values Description
PC-ANALYSIS int Analysis with the number int the number have to be followed by : FORMFINDING
DOMAIN EL-DOMAIN int Domain which should be calculated. The parameter EL-DOMAIN should be defined in the Element Block of the Inputfile
OUTPUT PC-OUT int Definition of the Output object. The Parameter PC-OUT should be defined in the PC-Problem Block of the Inputfile
SOLVER PC-SOLVER int Definition of the Solver object. The Parameter PC-SOLVER should be defined in the PC-Problem Block of the Inputfile
COMPCASE LD-COM int Definition of the Computation Case. The Parameter LD-COM should defined in the Load-Block of the Inputfile
FORMFINDING_STEP int Definition of the number of form finding steps
MAX_ITER_EQUILIBRIUM int Maximum of number of equilibrium iterations before CARAT++ stops the calculation
EQUILIBRIUM_ACCURACY real Equilibrium accuracy
FORMFINDING_ELEMENTS PROP_ID ints elements that have to be formfound (i.e. that have no material contribution, but only geometrical stiffness), identified through their property ID
Optional parameters for full URS
HOMOTOPY_CURVE LD-CURVE int Defines the homotopy factor (has to be between 0 and 1) for the different form-finding steps in a separate curve. If not used, the homotopy factor is set to 0, i.e. the reference problem is solved.
Optional parameters for stress adaption / distortion control
STRESS_ADAPTION PROP_ID ints Defines for which elements the distortion control is active through their property ID
LAMBDA_MAX real Maximum allowable element distortion in case of stress adaption
MUE_RELAX real Relaxation of the element distortion (Values between 0 and 1 are valid)
Optional parameters for restart file
ELASTIC_RESTART_METHOD string Defines the update of elastic members for the restart file (NONE = everything is updated and results set to 0, FORCE = geometry is updated, results set to 0 and inital stresses added, DISPLACEMENT = geometry of elastic members is not updated, initial displacement is added and nodes with different neighbouring elements are doubled)
Optional parameters for eXtended URS
X_URS bool Defines whether or not the xURS shall be applied to the computation resp. to the Membrane1 properties; Default is FALSE, leading to normal URS
X_URS_STEP int Defines the form-finding step from which on xURS is applied

For the better understanding of the syntax of this kind of analysis in the following an example of an input block for a form finding will be shown.

Example of an Input Block

The input block shown below defines a formfinding analysis with the IDs 1,2 and 8. The analysis performs 10 formfinding steps and stops the computation after 100 equilibrium iteration steps within each formfinding step if no convergence is achieved. The equilibrium accuracy of the analysis is defined as 10-6.

The (optional) distortion control is enabled and used for elements releated to the elemnt properties with ID 1. The maximum allowed element distortion is 1.1 and the relaxation of the element distortion control is set to 0. For a detailed discussion of the element distortion control see e.g. to [1] and [2].

PC-ANALYSIS 1: FORMFINDING
  DOMAIN = EL-DOMAIN 1
  OUTPUT = PC-OUT 1
  SOLVER = PC-SOLVER 1
  COMPCASE = LD-COM 1
  FORMFINDING_STEP = 10
  MAX_ITER_EQUILIBRIUM = 100
  EQUILIBRIUM_ACCURACY = 1e-06
  FORMFINDING_ELEMENTS = PROP_ID 1,2,8
  STRESS_ADAPTION = PROP_ID 1
  LAMBDA_MAX=1.1
  MUE_RELAX=0.0

Benchmark examples

  • form-finding of a four-point sail: ..\examples\benchmark_examples\analyses\formfinding_membrane1_I\cbm_4_point_Fofi.txt
  • form-finding of a catenoid with xURS: ..\examples\benchmark_examples\analyses\cbm_formfinding_xURS_Membrane1\cbm_formfinding_xURS_Membrane1.dat

References

  1. Linhard, J.: Numerisch-mechanische Betrachtung des Entwurfsprozesses von Membrantragwerken, Lehrstuhl für Statik, Technische Universität München, 2009
  2. Dieringer, F.: Numerical Methods for the Design and Analysis of Tensile Structures, Lehrstuhl für Statik, Technische Universität München, 2014




Whos here now:   Members 0   Guests 0   Bots & Crawlers 1
 
Personal tools
Content for Developers