Exporting only one animation to gltf from multiple objects each with an action I want to do an animated gltf file that shows how a set of molecules move.  Using python, I created an object for each atom, 1558 atoms in total and then I added keyframes  for updating the locations each frame.  Then I exported to gltf and I can see the molecules moving alright in the online viewer.  The problem is that there is an animation per each atom and I want only one animation in the file that includes all the atoms.

I tried putting all the actions in a single NLA strip but when I export the file I don't even see the strip in the viewer, not sure why (I do select the option for exporting NLA strips) and I still see the 1558 different animations.
Any help would be appreciated.
Update
If I create a few spheres and add keyframes manually, then everything is exported to a single animation as I want.  When I compare this to the file generated with python, I see that the F curves created for the objects actions are grouped into an ActionGroup called 'object tranfroms' but that group is not generated when I insert the keyframes by scripting.  I don't know if that grouping makes the GLTF have a single animation as I want.
How can I insert a keyframe programmatically  that exactly replicates what happens when you manually press I over the object to insert it?  I already tried adding an ActionGroup to the object but it is a read-only type and I think you can't add a F Curve to it.
Not sure why but I got it to work by grouping the channels using this code:
for obj in bpy.context.selected_objects:
    new_group=obj.animation_data.action.groups.new(&quot;ChGroup&quot;)
    for fc in obj.animation_data.action.fcurves:
        fc.group=new_group

And then exporting to gltf without selecting NLA Track:
bpy.ops.export_scene.gltf(export_format='GLB', filepath=filepath, export_nla_strips =False )

Now there is only one animation in the gltf (the color is different too but that's unrelated):