[Geoqus] Sequentially Thermal Coupling

[Manooch Salehabadi]

Dear All,

 

Hi,

 

So far, I have made two different models (for studying oil well integrity in deep offshore environment) one Hydraulic-Mechanic (once again thanks for Wouter and Jose comments) and one Thermal. I have used CPE8RT element for Thermal analysis and CPE8RP for Hydraulic-Mechanic analysis.

Now I am trying sequentially coupling these models (similar to solved example of ABAQUS Manual) and look at pore pressure increasing due to thermal loading but I am facing a problem, which is different time scale for these two models.

I will be much appreciating for any comment, paper, reference or similar input file, which can guide me to sort out this problem.

 

Regards

Manooch 

 

The following is some parts of my H-M and Thermal models:

 

 

 

 

H-M MODEL

 

*Boundary

FAR-FIELD-BOUNDARY-NODES, 1, 1

** Name: BC-DIS-FAR-FIELD-2 Type: Displacement/Rotation

*Boundary

FAR-FIELD-BOUNDARY-NODES, 2, 2

*INITIAL CONDITIONS,TYPE=RATIO

TOTAL,0.33

*INITIAL CONDITIONS,TYPE=PORE PRESSURE

TOTAL,18

*INITIAL CONDITIONS,TYPE=STRESS

TOTAL,-6,-6,-6

** ----------------------------------------------------------------

** 

** STEP: Geostatic

** 

*Step, name=Geostatic

*Geostatic

** 

*MODEL CHANGE,REMOVE

CEMENT-CASING-ADD

** BOUNDARY CONDITIONS

** 

** Name: BC-DIS-CEMENT-OUTER-1 Type: Displacement/Rotation

*Boundary

CEMENT-OUTER-FACE-NODES, 1, 1

** Name: BC-DIS-CEMENT-OUTER-2 Type: Displacement/Rotation

*Boundary

CEMENT-OUTER-FACE-NODES, 2, 2

** Name: BC-POR-FAR-FIELD Type: Pore pressure

*Boundary

FAR-FIELD-BOUNDARY-NODES_PP_, 8, 8, 18.

** 

** OUTPUT REQUESTS

** 

*Restart, write, frequency=0

** 

** FIELD OUTPUT: F-Output-1

** 

*Output, field

*Element Output, directions=YES

ALPHA, CTSHR, DAMAGEC, DAMAGET, E, EE, ER, FLUVR, FLVEL, GELVR, IE, LE, MISESMAX, NE, PE, PEEQ

PEEQMAX, PEEQT, PEMAG, PEQC, PS, RD, S, SALPHA, SAT, SDEG, SE, SEE, SEP, SEPE, SPE, THE

TSHR, VE, VEEQ, VOIDR, VS, VVF, VVFG, VVFN

** 

** FIELD OUTPUT: F-Output-2

** 

*Node Output

CFF, POR, RVF, RVT, U, UR, UT, V

VR, VT

** 

** HISTORY OUTPUT: H-Output-1

** 

*Output, history, variable=PRESELECT

*End Step

** ----------------------------------------------------------------

** 

** STEP: DRILLING-RUNNING-CASING

** 

*Step, name=DRILLING-RUNNING-CASING,nlgeom, amplitude=RAMP, inc=1000

step 3

*Soils, consolidation, end=PERIOD, utol=5., creep=none

0.1, 10., 0.0001, 10., 

** 

*MODEL CHANGE,REMOVE

DRILLED-WELLBORE

*MODEL CHANGE,ADD

CEMENT-CASING-ADD

** LOADS

** 

** Name: Load-1   Type: Pressure

*Dsload

Surf-1, P, 20.

** 

** BOUNDARY CONDITIONS

** 

** Name: BC-DIS-CEMENT-OUTER-1 Type: Displacement/Rotation

*Boundary, op=NEW

** Name: BC-DIS-CEMENT-OUTER-2 Type: Displacement/Rotation

*Boundary, op=NEW

** Name: BC-DIS-DRILLED-SECTION Type: Displacement/Rotation

*Boundary, op=NEW

DRILLED-EXCEPT-WELLBORE-FACE, 1, 1

DRILLED-EXCEPT-WELLBORE-FACE, 2, 2

** Name: BC-DIS-DUMMY-NODES Type: Displacement/Rotation

*Boundary, op=NEW

DUMMY-NODES, 1, 1

DUMMY-NODES, 2, 2

** Name: BC-DIS-FAR-FIELD-1 Type: Displacement/Rotation

*Boundary, op=NEW

FAR-FIELD-BOUNDARY-NODES, 1, 1

** Name: BC-DIS-FAR-FIELD-2 Type: Displacement/Rotation

*Boundary, op=NEW

FAR-FIELD-BOUNDARY-NODES, 2, 2

** Name: BC-POR-DRILLED-SECTION Type: Pore pressure

*Boundary, op=NEW

DRILLED-EXCEPT-WELLBORE-FACE_PP_, 8, 8, 20.

** Name: BC-POR-DRILLING-MUD Type: Pore pressure

*Boundary, op=NEW

WELLBORE-FACE-NODES_PP_, 8, 8, 20.

** Name: BC-POR-FAR-FIELD Type: Pore pressure

*Boundary, op=NEW

FAR-FIELD-BOUNDARY-NODES_PP_, 8, 8, 18.

** 

** OUTPUT REQUESTS

**

*Restart, write, frequency=0

** 

** FIELD OUTPUT: F-Output-3

** 

*Output, field

*Element Output, directions=YES

ALPHA, CTSHR, DAMAGEC, DAMAGET, E, EE, ER, FLUVR, FLVEL, GELVR, IE, LE, MISESMAX, NE, PE, PEEQ

PEEQMAX, PEEQT, PEMAG, PEQC, PS, RD, S, SALPHA, SAT, SDEG, SE, SEE, SEP, SEPE, SPE, THE

TSHR, VE, VEEQ, VOIDR, VS, VVF, VVFG, VVFN

** 

** FIELD OUTPUT: F-Output-4

** 

*Node Output

CFF, POR, RVF, RVT, U, UR, UT, V

VR, VT

** 

** HISTORY OUTPUT: H-Output-2

** 

*Output, history, variable=PRESELECT

*End Step

** ----------------------------------------------------------------

** 

** STEP: AFTER-RUNNING-CASING

** 

*Step, name=AFTER-RUNNING-CASING,nlgeom, amplitude=RAMP, inc=1000

*Soils, end=ss ,utol=1.0, creep=none

10., 10000., 1.,10000.,0.01 

** 

** OUTPUT REQUESTS

** 

*Restart, write, frequency=0

** 

** FIELD OUTPUT: F-Output-3

** 

*Output, field

*Element Output, directions=YES

ALPHA, CTSHR, DAMAGEC, DAMAGET, E, EE, ER, FLUVR, FLVEL, GELVR, IE, LE, MISESMAX, NE, PE, PEEQ

PEEQMAX, PEEQT, PEMAG, PEQC, PS, RD, S, SALPHA, SAT, SDEG, SE, SEE, SEP, SEPE, SPE, THE

TSHR, VE, VEEQ, VOIDR, VS, VVF, VVFG, VVFN

** 

** FIELD OUTPUT: F-Output-4

** 

*Node Output

CFF, POR, RVF, RVT, U, UR, UT, V

VR, VT

** 

** HISTORY OUTPUT: H-Output-2

** 

*Output, history, variable=PRESELECT

*End Step

 

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

THERMAL MODEL

 

** PREDEFINED FIELDS

** 

** Name: Field-1   Type: Temperature

*Initial Conditions, type=TEMPERATURE

_PICKEDSET11, 288.

** ----------------------------------------------------------------

** 

** STEP: THERMAL

** 

*Step, name=THERMAL, amplitude=RAMP, inc=1000

*Coupled Temperature-Displacement, creep=none, deltmx=1.

1., 14400., 0.144, 14400.

** 

** BOUNDARY CONDITIONS

** 

** Name: Temp-BC-1 Type: Temperature

*Boundary

WELLBORE-FACE-WITH-CASING, 11, 11, 298.

** 

** OUTPUT REQUESTS

** 

*Restart, write, frequency=0

** 

** FIELD OUTPUT: F-Output-1

** 

*Output, field

*Node Output

NT, 

** 

** FIELD OUTPUT: F-Output-2

** 

*Element Output, directions=YES

TEMP, 

** 

** HISTORY OUTPUT: H-Output-1

** 

*Output, history, variable=PRESELECT

*End Step

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.ruhr-uni-bochum.de/pipermail/geoqus/attachments/20070505/8be8ed28/attachment.html>