site map | contact | login | Protection des données personnelles | Powered by FluxBB | réalisation artaban

You are not logged in.

- Topics: Active | Unanswered

Pages: **1**

Hello Jean Pierre,

Thank you for your quick answer.

I am trying to add bolts to another model that has been built using COQUE_3D.

To my anderstanding, COQUE_3D is adapted to geometries with high curvature.

Best regards,

jeanpierreaubry wrote:

hello

here are

the gmsh script to build the geometry

the .med filefrom it

the .comm filenote

2 differents ways of handling the bolt to plate joint

some lines have to be commented or uncommented according to the calculation typeif anything looks stange ask

jean pierre aubry

Dear Jean Pierre,

Thank you for the valuable files, they're really helped me.

I reproduce your case, it works just fine using DKT.

But when I change it to COQUE_3D. Of course i converted the mesh to quadratic and then QUAD9, also I converted the bolt lines to lineare mesh since POU_D_T works only with SEG2 elements.

You can find the figure of comparaison attached.

Thank you very much,

Regards

Hello,

I am facing exactly the same error on code_aster windows

!-----------------------------------------------------------------------------------------------!

! <EXCEPTION> <JEVEUX1_13> !

! !

! Le rÃ©pertoire de noms &CATA.GD.NOMGD $$XNUM contient 168 points d'entrÃ©e, la requÃªte !

! sur le numÃ©ro 0 est invalide. !

! Ce !

! message est un message d'erreur dÃ©veloppeur. !

! Contactez le support technique. !

!-----------------------------------------------------------------------------------------------!

it appears when I want to write results extracted from a base using the POURSUITE command.

Regards,

**vietazz**- Replies: 0

Hello everyone,

I am studying a very simple case of a 1m height cylinder of steel with a base using elastic analysis.

A pressure of 0.1 MPa has been applied on its lateral wall, nothing is applied on its base other than a fixation of the center.

I studied two radii 0.5m and 1m.

all cases have been studied using two methods: COQUE_3D and AXIS.

For r=0.5m, the displacement and the behavior is similar between COQUE_3D and AXIS (as you can see attached)

for r=1m, I couldn't get a convergence even if I refined more than once the mesh, I am really stuck.

```
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! Solver MUMPS: !
! The solution of the linear system is too vague: !
! Computed error: 0.00112041 !
! Acceptable error: 0.0002 !
! (RESI_RELA) !
! !
! Advices: !
! One can increase the value of the key word SOLVER/RESI_RELA. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
```

I tried to block the top too, but the results are non physical.

All the cases have been given in the attachment.

Any help would be very appreciated

Thank you in advance.

**vietazz**- Replies: 0

Hello everyone,

I am using THER_NON_LINE where the thermal conductivity depends in temperature.

I would like to know whether there is a way to export the conductivity to visualize it on ParaVis.

Thank you in advance,

Best regards

Hello,

Thank you for your reply.

I don't know how to get in touch with an expert to answer my question.

I'll be gratefull if you have any recommendations.

Best regards,

**vietazz**- Replies: 2

Hello everyone,

I succeeded to implement a thermal model of a moving plate (see the figure in the zip attached)

I also attached the case and the .comm to this thread.

My case is a non linear steady problem where the conductivity varies with temperature.

and since my plate is moving, I used THER_NON_LINE_MO. I created a field with CREA_CHAM that will hold the values of the velocity.

```
speed = CREA_CHAMP(
AFFE=_F(
NOM_CMP=('DX', 'DY'),
TOUT='OUI',
VALE=(0.1, 0.1)
),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_DEPL_R'
)
```

This field has being used to define the velocity in AFFE_CHAR_THER

```
load = AFFE_CHAR_THER(
CONVECTION=_F(
VITESSE=speed
),
FLUX_REP=_F(
FLUN=0.0,
GROUP_MA=('wall_b', 'wall_t')
),
MODELE=model,
TEMP_IMPO=(_F(
GROUP_MA=('wall_l', ),
TEMP=500.0
), _F(
GROUP_MA=('wall_r', ),
TEMP=100.0
))
)
```

The results were satisfying.

Now I would like to have a formula for the velocity (no longer constant).

So Here I defined a formula of my velocity (veloX,veloY)

```
veloX = FORMULE(
NOM_PARA=('X', 'Y'),
VALE='X/100'
)
```

So I changed the previously created field accordingly (I changed NOEU_DEPL_R to NOEU_DEPL_F, and entered my formulas) as shown below.

```
speed = CREA_CHAMP(
AFFE=_F(
NOM_CMP=('DX', 'DY'),
TOUT='OUI',
VALE_F=(veloX, veloY)
),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_DEPL_F'
)
```

The problem is when I run the simulation I get the following problem

```
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! le champ speed n'a pas la bonne grandeur : !
! grandeur autorise :DEPL_R !
! grandeur du champ :DEPL_F !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
```

After I have read in documentation of AFFE_CHAR_THER, I found the following :

SPEED =

Name of the field speed at the moment when calculation is carried out.

This field is a concept cham_no of type cham_no_depl_r. It must have been defined on all the

model for which one carries out calculation.

My question is this : is there a way in code_Aster to define a position dependent speed ?

Any help will be very appreciated.

regards

Hi,

I really struggled long time with this same problem.

I would like to know weather you succeeded to find a solution.

Any help will be appreciated.

Regards,

**vietazz**- Replies: 0

Bonjour tout le monde.

I am struggling to simulate steady heat transfer in a hollow rotating cylinder (a case similar to the Rotary Kiln Furnace). please find attached a representation of the case.

The angular velocity "omega" is constant, and the boundary conditions at the inner face are not uniform (two temperatures Tf and Tc are applied at different arcs of the inner face). The Fourier equation of temperature (to resolve) in a cylindrical coordinates (r,Theta) writes:

rho*cp*omega*dT/dTheta = (lambda) Laplacian(T)

(I didn't expand the Laplacian in (r,Theta) coordinates)

Is there a way to resolve that equation (case) using Code_Aster ?

I’ll be grateful to any suggestion

Thank you

Pages: **1**