Well, not-so-popular solution could be that you split one of the faces in a couple of surfaces, and create contact with different friction coeff. I know it is probably not the way you wanted, but could be your last resort..

Thanks for your answer Emir Nezirić . I already tried this solution with the script below, but the computation time has exploded...

This script creates, for each iteration, a contact between a random number of nodes of the upper part and the surface of the lower part (with a random friction coefficient).

In this case, I considered between 15 and 50 nodes of the upper part for each contact created.

With this range, around 10 contacts are created.

import random

model = ExtAPI.DataModel.Project.Model

#Nodes in contact of the upper body

IDs = [i for i in range(9780,10167)]

x=0

limite=0

limiteNM1=0

count = 0

for i in IDs:

mu = random.uniform(0.8, 1)

nbNodes = random.randint(15, 50)

limite += nbNodes

newContact=model.Connections.AddContactRegion()

newContact.ContactType=ContactType.Frictional

newContact.FrictionCoefficient = mu

newContact.DynamicCoefficient = mu

newContact.Behavior = ContactBehavior.Asymmetric

newContact.ContactFormulation=ContactFormulation.NormalLagrange

if limite > len(IDs):

limite = len(IDs)

#Creation d'une named selection contenant des noeud

ns = Model.AddNamedSelection()

ns.ScopingMethod=GeometryDefineByType.Worksheet

Criteria = ns.GenerationCriteria

Criteria.Add()

Criteria[0].Action = SelectionActionType.Add

Criteria[0].EntityType = SelectionType.MeshNode

Criteria[0].Criterion = SelectionCriterionType.NodeNumber

Criteria[0].Operator = SelectionOperatorType.RangeInclude

Criteria[0].LowerBound = IDs[limiteNM1]

Criteria[0].UpperBound = IDs[limite-1]

ns.Generate()

#Verif contact

contact = ExtAPI.DataModel.GetObjectsByName(model.NamedSelections.Children[x+1].Name)

if contact:

newContact.SourceLocation = contact[0]

target = ExtAPI.DataModel.GetObjectsByName('bas')

if target:

newContact.TargetLocation = target[0]

x = x+1

for i in range(limite-limiteNM1):

ligne = [i, mu]

muVector.append(ligne)

limiteNM1 = limite

break

elif i == IDs[0]:

#Creation d'une named selection contenant des noeud

ns = Model.AddNamedSelection()

ns.ScopingMethod=GeometryDefineByType.Worksheet

Criteria = ns.GenerationCriteria

Criteria.Add()

Criteria[0].Action = SelectionActionType.Add

Criteria[0].EntityType = SelectionType.MeshNode

Criteria[0].Criterion = SelectionCriterionType.NodeNumber

Criteria[0].Operator = SelectionOperatorType.RangeInclude

Criteria[0].LowerBound = IDs[limiteNM1]

Criteria[0].UpperBound = IDs[limite-1]

Criteria.Add()

Criteria[1].Action = SelectionActionType.Add

Criteria[1].EntityType = SelectionType.MeshNode

Criteria[1].Criterion = SelectionCriterionType.NodeNumber

Criteria[1].Operator = SelectionOperatorType.Equal

Criteria[1].Value = 9392

ns.Generate()

ExtAPI.DataModel.Tree.Refresh()

#Verif contact

contact = ExtAPI.DataModel.GetObjectsByName(model.NamedSelections.Children[x+1].Name)

if contact:

newContact.SourceLocation = contact[0]

target = ExtAPI.DataModel.GetObjectsByName('bas')

if target:

newContact.TargetLocation = target[0]

x = x+1

for i in range(count, count+nbNodes):

ligne = [i, mu]

muVector.append(ligne)

limiteNM1 = limite

count = count + nbNodes

else:

#Creation d'une named selection contenant des noeud

ns = Model.AddNamedSelection()

ns.ScopingMethod=GeometryDefineByType.Worksheet

Criteria = ns.GenerationCriteria

Criteria.Add()

Criteria[0].Action = SelectionActionType.Add

Criteria[0].EntityType = SelectionType.MeshNode

Criteria[0].Criterion = SelectionCriterionType.NodeNumber

Criteria[0].Operator = SelectionOperatorType.RangeInclude

Criteria[0].LowerBound = IDs[limiteNM1]

Criteria[0].UpperBound = IDs[limite-1]

ns.Generate()

ExtAPI.DataModel.Tree.Refresh()

#Verif contact

contact = ExtAPI.DataModel.GetObjectsByName(model.NamedSelections.Children[x+1].Name)

if contact:

newContact.SourceLocation = contact[0]

target = ExtAPI.DataModel.GetObjectsByName('bas')

if target:

newContact.TargetLocation = target[0]

x = x+1

limiteNM1 = limite

count = count + nbNodes