@@ -170,7 +170,6 @@ defs = [
170
170
display(sort(unknowns(ssys), by=string))
171
171
172
172
prob = ODEProblem(ssys, defs, (0, 2))
173
- # sol.u[1] = [0.17054946565059462, -0.6015077878288565, 0.274732819217001]
174
173
sol = solve(prob, FBDF(autodiff=true))
175
174
@test SciMLBase.successful_retcode(sol)
176
175
rms(x) = sqrt(sum(abs2, x) / length(x))
@@ -259,7 +258,6 @@ display(sort(unknowns(ssys), by=string))
259
258
prob = ODEProblem(ssys, defs, (0, 4))
260
259
sol = solve(prob, FBDF(autodiff=true), initializealg = ShampineCollocationInit())
261
260
@test SciMLBase.successful_retcode(sol)
262
- # first(Multibody.render(model, sol, 0, show_axis=true, x=-1.5, y=0.0, z=0.0))
263
261
```
264
262
265
263
``` @example suspension
@@ -363,7 +361,6 @@ defs = [
363
361
# model.body_upright.prismatic_y.s => 0.17
364
362
# model.body_upright.prismatic_y.v => 0.14
365
363
366
-
367
364
# vec(ori(model.mass.frame_a).R .=> I(3))
368
365
# vec(ori(model.excited_suspension_r.suspension.r123.jointUSR.frame_a).R .=> I(3))
369
366
]
@@ -372,6 +369,8 @@ display(sort(unknowns(ssys), by=string))
372
369
373
370
prob = ODEProblem(ssys, defs, (0, 4))
374
371
sol = solve(prob, FBDF(autodiff=false), initializealg = ShampineCollocationInit())#, u0 = prob.u0 .+ 1e-6 .* randn.())
375
- @show SciMLBase.successful_retcode(sol)
376
- first(Multibody.render(model, sol, 0, show_axis=true, x=-1.5, y=0.0, z=0.0))
377
- ```
372
+ @test SciMLBase.successful_retcode(sol)
373
+ Multibody.render(model, sol, show_axis=false, x=-1.5, y=0.3, z=0.0, lookat=[0,0.1,0.0], timescale=3, filename="suspension_halfcar_wheels.gif") # Video
374
+ ```
375
+
376
+ ![ suspension with wheels] ( suspension_halfcar_wheels.gif )
0 commit comments