Mercurial > repos > public > sbplib

changeset 617:4ced7d47bd1f feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Beam: Add ability to set one BC at a time
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 20 Oct 2017 15:56:10 +0200
parents 818d52d4928f
children c360bbecf260
files +scheme/Beam.m
diffstat 1 files changed, 38 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
diff -r 818d52d4928f -r 4ced7d47bd1f +scheme/Beam.m
--- a/+scheme/Beam.m	Wed Oct 11 16:01:22 2017 +0200
+++ b/+scheme/Beam.m	Fri Oct 20 15:56:10 2017 +0200
@@ -126,6 +126,44 @@
                     penalty{1} = -obj.Hi*tau;
                     penalty{1} = -obj.Hi*sig;
 
+                case 'e'
+                    alpha = obj.alpha;
+                    tuning = 1.1;
+
+                    tau1 = tuning * alpha/delt;
+                    tau4 = s*alpha;
+
+                    tau = tau1*e+tau4*d3;
+
+                    closure = obj.Hi*tau*e';
+                    penalty = -obj.Hi*tau;
+                case 'd1'
+                    alpha = obj.alpha;
+
+                    tuning = 1.1;
+
+                    sig2 = tuning * alpha/gamm;
+                    sig3 = -s*alpha;
+
+                    sig = sig2*d1+sig3*d2;
+
+                    closure = obj.Hi*sig*d1';
+                    penalty = -obj.Hi*sig;
+
+                case 'd2'
+                    a = obj.alpha;
+
+                    tau =  s*a*d1;
+
+                    closure = obj.Hi*tau*d2';
+                    penalty = -obj.Hi*tau;
+                case 'd3'
+                    a = obj.alpha;
+
+                    sig = -s*a*e;
+
+                    closure = obj.Hi*sig*d3';
+                    penalty = -obj.Hi*sig;
 
                 otherwise % Unknown, boundary condition
                     error('No such boundary condition: type = %s',type);