interaction.txt 2.16 KB
Newer Older
1
2
.. py:currentmodule:: pylmgc90.pre_lmgc

mozul's avatar
mozul committed
3
Interaction model and parameters definition
4
5
=========================================

mozul's avatar
mozul committed
6
One needs to define a tact_behav object with :py:func:`tact_behav` which corresponds to an
7
8
interaction model with its parameters.

mozul's avatar
mozul committed
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Meaning of the parameters

* friction: either constant *fric* or *stfr* static and *dyfr* dynamics
* restitution (newton): *rstn* normal, *rstt* tangential
* cohesion: *cohn* normal, *coht* tangential (unused with wet model)
* CZM parameters:

  * *cn* normal stiffness (pure mode I)
  * *ct* tangential stiffness (pure mode II)
  * *smax*  critical stress (with mixity it is pure mode I)
  * *w*  fracture energy (with mixity it is pure mode I)
  * *d* plastic displacement for TH (d=smax/cn means triangle)
  * *du*  rupture displacement for ABP
  * *phi* ratio between micro and macro fracture energy
  * *p* constant pressure in a broken CZM

25
26
27
28
29
30
31
32
33
34
35
36
Here is a list of existing models with the expected parameters:

RIGID/RIGID :
-------------

 * IQS_CLB: fric
 * IQS_CLB_g0: fric
 * IQS_DS_CLB: dyfr, stfr
 * IQS_WET_DS_CLB: cohn, coht, Wthk, dyfr, stfr
 * IQS_MOHR_DS_CLB: cohn, coht, dyfr, stfr
 * IQS_MAC_CZM: dyfr, stfr, cn, ct , b, w  
 * IQS_MAL_CZM: dyfr, stfr, cn, ct, smax, w
mozul's avatar
mozul committed
37
38
 * IQS_TH_CZM: dyfr, stfr, cn, ct, p, p0, smax, w , d
 * IQS_ABP_CZM: dyfr, stfr, cn, ct, p, smax, w , du, phi  
39
40
41
42
43
44
45
46
47
48
49
 * RST_CLB: rstn, rstt, fric

any/DEFORMABLE :
----------------

 * GAP_SGR_CLB: fric 
 * GAP_SGR_CLB_g0: fric
 * GAP_MOHR_DS_CLB: cohn, coht, dyfr, stfr
 * VEL_SGR_CLB: fric 
 * MAC_CZM: dyfr, stfr, cn, ct, b, w
 * MAL_CZM: dyfr, stfr, cn, ct, smax, w
mozul's avatar
mozul committed
50
51
52
53
54
55
 * MP_CZM: dyfr, stfr, cn, ct, p, p0, w
 * MSMP_CZM: dyfr, stfr, cn, ct, beta0, w
 * MP3_CZM: dyfr, stfr, cn, ct, p, p0, smax, w
 * TH_CZM: dyfr, stfr, cn, ct, p, p0, smax, w , d
 * ABP_CZM: dyfr, stfr, cn, ct, p, smax, w , du, phi
   
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
POINT/POINT:
------------

 * ELASTIC_WIRE: stiffness, prestrain
 * BRITTLE_ELASTIC_WIRE': stiffness, prestrain, Fmax
 * ELASTIC_ROD: stiffness, prestrain

any/any:
--------

 * COUPLED_DOF:
 * NORMAL_COUPLED_DOF:
 * ELASTIC_REPELL_CLB: stiffness, fric

**Example:** ::

  csas=tact_behav(name='gapc0', type='MAC_CZM', dyfr=0., stfr=0., cn=4.e+14, ct=3.e+14, b=0., w=60.)