Subsections

Initial Conditions

Initial conditions are entered as FEMGEN loads via the 'PROPERTY LOAD' or the 'PROPERTY INITIAL' command. The following initial conditions are supported:

 Initial Condition Mass Flow Rate Pore Pressure Void Ratio Rotating Velocity Saturation Stress Geostatic Stress Temperature Velocity Fluid Pressure

Note:

1.

The following sections describe the definition of these initial conditions.

Mass Flow Rate

There are 3 arguments, at least one of which one must be entered depending on the dimensionality of the elements to which the conditions apply.

 Argument(s) Meaning 1 Mass flow rate per unit area (x) 2 Mass flow rate per unit area (y) 3 Mass flow rate per unit area (z)

Example FEMGEN Command:

PROPERTY LOAD I-MFLOW S1 1 2 3

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=MASS FLOW RATE
S1  ,  1.000000    ,  2.000000    ,  3.000000

Pore Pressure

There is 1 argument, which must be entered.

 Argument Meaning 1 Fluid pore pressure

Example FEMGEN Command:

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=PORE PRESSURE
S1  ,  1.000000

Void Ratio

There are 4 arguments, all of which must be entered.

 Arguments Meaning 1 First value of void ratio 2 Vertical coordinate for the above value 3 Second value of void ratio 4 Vertical coordinate for the above value

Example FEMGEN Command:

PROPERTY LOAD I-RATIO S1 1 2 3 4

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=RATIO
S1  ,  1.000000    ,  2.000000    ,  3.000000    ,  4.000000

Rotating Velocity

There are 10 arguments, all of which must be entered.

 Arguments Meaning 1 Angular Velocity about axis a b. 2 Global x-component of translational velocity 3 Global y-component of translational velocity 4 Global z-component of translational velocity 5 Global x-component of point a 6 Global y-component of point a 7 Global z-component of point a 8 Global x-component of point b 9 Global y-component of point b 10 Global z-component of point b

Example FEMGEN Command:

PROPERTY LOAD I-RVEL S1 1 2 3 4 5 6 7 8 9 10

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=ROTATING VELOCITY
S1  ,  1.000000    ,  2.000000    ,  3.000000    ,  4.000000
5.000    , 6.000    , 7.000    , 8.000    , 9.000    , 10.00

Saturation

There is 1 argument, which must be entered.

 Arguments Meaning 1 Saturation value

Example FEMGEN Command:

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=SATURATION
S1  ,  1.000000

Stress

There are 6 arguments, at least one of which must be entered depending upon the number of stress components which are applicable to the elements to which the conditions apply.

 Arguments Meaning 1 First stress component 2 Second stress component 3 Third stress component 4 Fourth stress component 5 Fifth stress component 6 Sixth stress component

Example FEMGEN Command:

PROPERTY LOAD I-STRESS S1 1 2 3 4 5 6

Corresponding ABAQUS Input File:

*INITIAL CONDITIONS, TYPE=STRESS
S1  ,  1.00    ,  2.00    ,  3.00    ,  4.00    ,  5.00    ,  6.00

Geostatic Stress

There are 6 arguments, all of which must be entered.

 Arguments Meaning 1 First value of vertical component of effective stress 2 Vertical coordinate for the above value 3 Second value of vertical component of effective stress 4 Vertical coordinate for the above value 5 Coefficient of lateral stress (x-direction) 6 Coefficient of lateral stress (y-direction)

Example FEMGEN Command:

PROPERTY LOAD I-GSTRES S1 1 2 3 4 5 6

Corresponding ABAQUS Input File:

*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC
S1  ,  1.00    ,  2.00    ,  3.00    ,  4.00    ,  5.00    ,  6.00

Temperature

There are 7 arguments, of which at least one must be entered. Multiple temperature values are used in association with beam section points and through-thickness shell locations.

 Arguments Meaning 1 Temperature 2 Temperature 3 Temperature 4 Temperature 5 Temperature 6 Temperature 7 Temperature

Example FEMGEN Command:

PROPERTY LOAD I-TEMP S1 1 2 3 4 5

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=TEMPERATURE
S1      ,  1.000000    ,  2.000000    ,  3.000000    ,  4.000000    ,  5.000000

Velocity

There are 2 arguments, both of which must be entered.

 Arguments Meaning 1 Velocity 2 Degree of Freedom

Example FEMGEN Command:

PROPERTY LOAD I-VEL S1 50 1

Corresponding ABAQUS Input File:

*NSET, NSET=S1
1    2    3    4    5    6    7    8    9   10   11
*INITIAL CONDITIONS, TYPE=VELOCITY
S1  ,    1,  50.00000

Fluid Pressure

There is 1 argument, which must be entered.

 Argument Meaning 1 Fluid pressure

Example FEMGEN Command:

Corresponding ABAQUS Input File:

*NSET, NSET=P1
1
*INITIAL CONDITIONS, TYPE=FLUID PRESSURE
P1  ,  1.000000    Next: History Definition Up: Kinematic Constraints Previous: Kinematic Coupling Constraints

Femsys Limited
9/9/1999