RDforum.free.fr/Cast3M   Index du Forum
Site Web Cast3M
RDforum.free.fr/Cast3M
Forum d'entraide des utilisateurs de Cast3M
Répondre au sujet
Comment entrer des propriétés phy. variables par élément ?
Pascal Roux


Inscrit le: 02 Oct 2008
Messages: 3
Répondre en citant
Bonjour,

J'utilise une formulation poutre.
Je souhaiterais insérer des coefficients Young et des sections
VARIABLES par éléments.

j'ai vu dans la doc de MATER et CARAC que l'on peut entrer des
valeurs de type FLOTTANT, EVOL, ou MCHML pour les modules d'Young
et les sections.

supposons que je définisse un module d'Young et une section variable
avec une LISTREEL de n valeurs : n est le nombre d'éléments.

Soit x l'abscisse associée au centre de chaque élément :

e_de_x = ax + b ( avec la commande PROG LIN qui va bien )
s_de_x = cx + d


j'ai défini une ligne geo1 (type MAILLAGE) et un modèle mod1
je fais ceci :

you_var = MANU CHML geo1 YOU e_de_x;
sec_var = MANU CHML geo1 SECT s_de_x;

list you_var;
list sec_var;

Puis :
mat1 = MATER mod1 YOUN you_var NU 0.3;
car1 = CARA mod1 SECT sec_var INRY ... etc ;

Ok j'ai bien des coeff variables MAIS en sous type POINTEUR LISTREEL
Ils ne sont pas associés aux éléments. Ce qui plante ensuite dans le
calcul de la matrice de rigidité :

mat1 = mat1 ET car1 ; OK ça marche
rigi1 = RIGID mod1 mat1 ; ça plante :

***** ERREUR 769 ***** dans l'opérateur RIGI
2 sous zones ont le même maillage 47845 et le même constituant
(il met un large espace blanc) 47851


Donc une première question est :

comment utiliser MANU MCHML geo1 YOUN e_de_x;
pour obtenir un champ variable par éléments de type
réel et associé aux éléments de la géométrie geo1 ?
Est-ce que cela proviendrait de la commande PROG LIN
qui définit une liste de valeurs s'appuyant sur une liste
d'abscisses x n'ayant aucun rapport avec les éléments ?





Je fais une autre tentative en utilisant cette fois une variante
avec le mot clé REPA en s'appuyant sur le modèle mod1 :

mat1 = MANU CHML mod1 REPA YOUN e_de_x NU 0.3 RIGIDITE;
car1 = MANU CHML mod1 REPA SECT s_de_x INRY .... etc RIGIDITE;

list mat1;
list car1;

OK cette fois j'ai bien des coefficients variables par élément
de type réel et associé aux éléments.

mais ensuite, pour le calcul de la matrice de rigidité :

mat1 = mat1 ET car1; OK ça marche
rigi1 = RIGID mod1 mat1; ça plante avec le message :


**** ERREUR 130 ******** dans l'opérateur RIGI
L'objet de type CHAMELEM de sous type
(il met un large espace blanc) a déjà été donné


A noter que cela fonctionne SI l'on a QUE mat1 OU QUE car1
défini par MANU CHML mod1 REPA... RIGIDITE; mais PAS LES
DEUX A LA FOIS !

Soit on a :
mat1 = MANU CHML mod1 REPA YOUN e_de_x
car1 = CARAC mod1 SECT s_constante ....;

soit :

mat1 = MATER mod1 YOUN e_constant;
car1 = MANU CHML mod1 REPA SECT s_de_x...

mais PAS LES DEUX A LA FOIS !



Donc ma deuxième question est :

comment utiliser MANU MCHML mod1 REPA .... pour
obtenir des coefficients variables par éléments ?
Faut-il utiliser autre chose que le mot clé RIGIDITE à
la fin de la commande ?



Dernière question :

y a-t il une autre façon de procéder pour définir des
coefficients variable par éléments ?



Merci

----------------------------------------------
Pascal Roux
Post Doc en dynamique des structures à l'INRA
Voir le profil de l'utilisateurEnvoyer un message privé
lguenhael
Administrateur

Inscrit le: 25 Nov 2006
Messages: 146
Localisation: Laboratoire Roberval - UTC
Répondre en citant
Bonjour,

Un petit programme par cas considéré illustrant chacune de des erreurs rencontrées serait la bienvenue.

Cordialement.
Voir le profil de l'utilisateurEnvoyer un message privéEnvoyer un e-mail
Comment entrer des propriétés phy. variables par élément ?
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum
Vous ne pouvez pas joindre des fichiers
Vous ne pouvez pas télécharger des fichiers
Toutes les heures sont au format GMT + 2 Heures  
Page 1 sur 1  

  
  
 Répondre au sujet